Does rapid sequence divergence preclude RNA structure conservation in vertebrates?

Accelerated evolution of any portion of the genome is of significant interest, potentially signaling positive selection of phenotypic traits and adaptation. Accelerated evolution remains understudied for structured RNAs, despite the fact that an RNA's structure is often key to its function. RNA structures are typically characterized by compensatory (structure-preserving) basepair changes that are unexpected given the underlying sequence variation, i.e., they have evolved through negative selection on structure. We address the question of how fast the primary sequence of an RNA can change through evolution while conserving its structure. Specifically, we consider predicted and known structures in vertebrate genomes. After careful control of false discovery rates, we obtain 13 de novo structures (and three known Rfam structures) that we predict to have rapidly evolving sequences-defined as structures where the primary sequences of human and mouse have diverged at least twice as fast (1.5 times for Rfam) as nearby neutrally evolving sequences. Two of the three known structures function in translation inhibition related to infection and immune response. We conclude that rapid sequence divergence does not preclude RNA structure conservation in vertebrates, although these events are relatively rare.

The identification and functional annotation of RNA structures conserved in vertebrates.

Structured elements of RNA molecules are essential in, e.g., RNA stabilization, localization, and protein interaction, and their conservation across species suggests a common functional role. We computationally screened vertebrate genomes for conserved RNA structures (CRSs), leveraging structure-based, rather than sequence-based, alignments. After careful correction for sequence identity and GC content, we predict ∼516,000 human genomic regions containing CRSs. We find that a substantial fraction of human-mouse CRS regions (1) colocalize consistently with binding sites of the same RNA binding proteins (RBPs) or (2) are transcribed in corresponding tissues. Additionally, a CaptureSeq experiment revealed expression of many of our CRS regions in human fetal brain, including 662 novel ones. For selected human and mouse candidate pairs, qRT-PCR and in vitro RNA structure probing supported both shared expression and shared structure despite low abundance and low sequence identity. About 30,000 CRS regions are located near coding or long noncoding RNA genes or within enhancers. Structured (CRS overlapping) enhancer RNAs and extended 3' ends have significantly increased expression levels over their nonstructured counterparts. Our findings of transcribed uncharacterized regulatory regions that contain CRSs support their RNA-mediated functionality.

Genes controlling the activation of natural killer lymphocytes are epigenetically remodeled in intestinal cells from germ-free mice.

Remodeling of the gut microbiota is implicated in various metabolic and inflammatory diseases of the gastrointestinal tract. We hypothesized that the gut microbiota affects the DNA methylation profile of intestinal epithelial cells (IECs) which could, in turn, alter intestinal function. In this study, we used mass spectrometry and methylated DNA capture to respectively investigate global and genome-wide DNA methylation of intestinal epithelial cells from germ-free (GF) and conventionally raised mice. In colonic IECs from GF mice, DNA was markedly hypermethylated. This was associated with a dramatic loss of ten-eleven-translocation activity, a lower DNA methyltransferase activity and lower circulating levels of the 1-carbon metabolite, folate. At the gene level, we found an enrichment for differentially methylated regions proximal to genes regulating the cytotoxicity of NK cells (false-discovery rate < 8.9E-6), notably genes regulating the cross-talk between NK cells and target cells, such as members of the NK group 2 member D ligand superfamily Raet. This distinct epigenetic signature was associated with a marked decrease in Raet1 expression and a loss of CD56+/CD45+ cells in the intestine of GF mice. Thus, our results indicate that altered activity of methylation-modifying enzymes in GF mice influences the IEC epigenome and modulates the crosstalk between IECs and NK cells. Epigenetic reprogramming of IECs may modulate intestinal function in diseases associated with altered gut microbiota.-Poupeau, A., Garde, C., Sulek, K., Citirikkaya, K., Treebak, J. T., Arumugam, M., Simar, D., Olofsson, L. E., Bäckhed, F., Barrès, R. Genes controlling the activation of natural killer lymphocytes are epigenetically remodeled in intestinal cells from germ-free mice.

Improved peptide-MHC class II interaction prediction through integration of eluted ligand and peptide affinity data.

Major histocompatibility complex (MHC) class II antigen presentation is a key component in eliciting a CD4+ T cell response. Precise prediction of peptide-MHC (pMHC) interactions has thus become a cornerstone in defining epitope candidates for rational vaccine design. Current pMHC prediction tools have, so far, primarily focused on inference from in vitro binding affinity. In the current study, we collate a large set of MHC class II eluted ligands generated by mass spectrometry to guide the prediction of MHC class II antigen presentation. We demonstrate that models developed on eluted ligands outperform those developed on pMHC binding affinity data. The predictive performance can be further enhanced by combining the eluted ligand and pMHC affinity data in a single prediction model. Furthermore, by including ligand data, the peptide length preference of MHC class II can be accurately learned by the prediction model. Finally, we demonstrate that our model significantly outperforms the current state-of-the-art prediction method, NetMHCIIpan, on an external dataset of eluted ligands and appears superior in identifying CD4+ T cell epitopes.

Identification of two microRNA nodes as potential cooperative modulators of liver metabolism.

AIM: Hepatic insulin resistance is a hallmark of type 2 diabetes and non-alcoholic fatty liver disease. Dysregulation of microRNA (miRNA) expression in insulin-resistant livers might coordinate impaired hepatic metabolic function. Here, we aimed to discover miRNAs and their downstream targets involved in hepatic insulin resistance. METHODS: We determined miRNA expression profiles by small RNA sequencing of two mouse models of impaired hepatic insulin action: high-fat diet-induced obesity and liver-specific insulin receptor knockout. Conversely, we assessed the hepatic miRNA expression profile after treatment with the antidiabetic hormone, fibroblast growth factor 21 (FGF21). Ontology analysis of predicted miRNA gene targets was performed to identify regulated gene pathways. Target enrichment analysis and miRNA mimic overexpression in vitro were used to identify unified protein targets of nodes of regulated miRNAs. RESULTS: We identified an array of miRNA species regulated by impaired liver insulin action or after fibroblast growth factor 21 treatment. Ontology analysis of predicted miRNA gene targets identified pathways controlling hepatic energy metabolism and insulin sensitivity. We identified a node of two miRNAs downregulated in the livers of liver-specific insulin receptor knockout mice, miR-883b and miR-205, which positively regulate the expression of transcription factor zinc finger E-box-binding homeobox 1 (ZBED1). We found another node of two miRNAs upregulated in the livers of fibroblast growth factor 21-treated mice, miR-155-3p and miR-1968-5p, which canonically downregulates the caveola component, polymerase I and transcript release factor (PTRF), a gene previously implicated in hepatic energy metabolism. CONCLUSIONS: This study identifies two nodes of coregulated miRNAs that might coordinately control hepatic energy metabolism in states of insulin resistance.

miR-18b overexpression identifies mantle cell lymphoma patients with poor outcome and improves the MIPI-B prognosticator.

Recent studies show that mantle cell lymphoma (MCL) express aberrant microRNA (miRNA) profiles; however, the clinical effect of miRNA expression has not previously been examined and validated in large prospective homogenously treated cohorts. We performed genome-wide miRNA microarray profiling of 74 diagnostic MCL samples from the Nordic MCL2 trial (screening cohort). Prognostic miRNAs were validated in diagnostic MCL samples from 94 patients of the independent Nordic MCL3 trial (validation cohort). Three miRNAs (miR-18b, miR-92a, and miR-378d) were significantly differentially expressed in patients who died of MCL in both cohorts. MiR-18b was superior to miR-92a and miR-378d in predicting high risk. Thus, we generated a new biological MCL International Prognostic Index (MIPI-B)-miR prognosticator, combining expression levels of miR-18b with MIPI-B data. Compared to the MIPI-B, this prognosticator improved identification of high-risk patients with regard to cause-specific, overall, and progression-free survival. Transfection of 2 MCL cell lines with miR-18b decreased their proliferation rate without inducing apoptosis, suggesting that miR-18b may render MCL cells resistant to chemotherapy by decelerating cell proliferation. We conclude that overexpression of miR-18b identifies patients with poor prognosis in 2 large prospective MCL cohorts and adds prognostic information to the MIPI-B. MiR-18b may reduce the proliferation rate of MCL cells as a mechanism of chemoresistance.

Global hypomethylation is an independent prognostic factor in diffuse large B cell lymphoma.

Global hypomethylation has been linked to disease progression in several cancers, but has not been reported for Diffuse Large B Cell Lymphoma (DLBCL). This study aimed to assess global methylation in DLBCL and describe its prognostic value. Mean LINE1 methylation, a validated surrogate measure for global methylation, was measured in DNA from 67 tumor biopsies. Additionally, cell-free circulating DNA (cfDNA) in plasma samples from 74 patients was tested to assess the feasibility of global hypomethylation as a biomarker in liquid biopsies. LINE1 methylation was assessed using a commercially available kit, based on pyrosequencing of PCR amplified bisulfite-treated DNA. Global hypomethylation was detected in a subset of cases and was associated with poor overall survival in both tumor biopsies (P = .001) and cfDNA (P = .009). It was the strongest risk factor in multivariate analysis in both biopsies (HR: 10.65, CI: 2.03-55.81, P = .005) and cfDNA (HR: 11.87, CI: 2.80-50.20, P = .001), outperforming conventional clinical risk factors. Finally, hierarchical cluster analyses were performed for the cfDNA samples using previously published gene-specific methylation data. This analysis shows that global hypomethylation co-occurs with other epigenetic abnormalities, including DAPK1 promoter hypermethylation. In conclusion, we have shown that global hypomethylation is strongly associated with poor survival in DLBCL both when present in tumor biopsy DNA and when detected in plasma cfDNA, and has potential for clinical application as a prognostic biomarker.

15-year follow-up of the Second Nordic Mantle Cell Lymphoma trial (MCL2): prolonged remissions without survival plateau.

In recent decades, the prognosis of Mantle Cell Lymphoma (MCL) has been significantly improved by intensified first-line regimens containing cytarabine, rituximab and consolidation with high-dose-therapy and autologous stem cell transplantation. One such strategy is the Nordic MCL2 regimen, developed by the Nordic Lymphoma Group. We here present the 15-year updated results of the Nordic MCL2 study after a median follow-up of 11·4 years: For all patients on an intent-to-treat basis, the median overall and progression-free survival was 12·7 and 8·5 years, respectively. The MCL International Prognostic Index (MIPI), biological MIPI, including Ki67 expression (MIPI-B) and the MIPI-B including mIR-18b expression (MIPI-B-miR), in particular, significantly divided patients into distinct risk groups. Despite very long response durations of the low and intermediate risk groups, we observed a continuous pattern of relapse and the survival curves never reached a plateau. In conclusion, despite half of the patients being still alive and 40% in first remission after more than 12 years, we still see an excess disease-related mortality, even among patients experiencing long remissions. Even though we consider the Nordic regimen as a very good choice of regimen, we recommend inclusion in prospective studies to explore the benefit of novel agents in the frontline treatment of MCL.

Microbial biofilm as a smart material.

Microbial biofilm colonies will in many cases form a smart material capable of responding to external threats dependent on their size and internal state. The microbial community accordingly switches between passive, protective, or attack modes of action. In order to decide which strategy to employ, it is essential for the biofilm community to be able to sense its own size. The sensor designed to perform this task is termed a quorum sensor, since it only permits collective behaviour once a sufficiently large assembly of microbes have been established. The generic quorum sensor construct involves two genes, one coding for the production of a diffusible signal molecule and one coding for a regulator protein dedicated to sensing the signal molecules. A positive feedback in the signal molecule production sets a well-defined condition for switching into the collective mode. The activation of the regulator involves a slow dimerization, which allows low-pass filtering of the activation of the collective mode. Here, we review and combine the model components that form the basic quorum sensor in a number of Gram-negative bacteria, e.g., Pseudomonas aeruginosa.

Personalized therapy with peptide-based neoantigen vaccine (EVX-01) including a novel adjuvant, CAF®09b, in patients with metastatic melanoma.

The majority of neoantigens arise from unique mutations that are not shared between individual patients, making neoantigen-directed immunotherapy a fully personalized treatment approach. Novel technical advances in next-generation sequencing of tumor samples and artificial intelligence (AI) allow fast and systematic prediction of tumor neoantigens. This study investigates feasibility, safety, immunity, and anti-tumor potential of the personalized peptide-based neoantigen vaccine, EVX-01, including the novel CD8+ T-cell inducing adjuvant, CAF®09b, in patients with metastatic melanoma (NTC03715985). The AI platform PIONEERTM was used for identification of tumor-derived neoantigens to be included in a peptide-based personalized therapeutic cancer vaccine. EVX-01 immunotherapy consisted of 6 administrations with 5-10 PIONEERTM-predicted neoantigens as synthetic peptides combined with the novel liposome-based Cationic Adjuvant Formulation 09b (CAF®09b) to strengthen T-cell responses. EVX-01 was combined with immune checkpoint inhibitors to augment the activity of EVX-01-induced immune responses. The primary endpoint was safety, exploratory endpoints included feasibility, immunologic and objective responses. This interim analysis reports the results from the first dose-level cohort of five patients. We documented a short vaccine manufacturing time of 48-55 days which enabled the initiation of EVX-01 treatment within 60 days from baseline biopsy. No severe adverse events were observed. EVX-01 elicited long-lasting EVX-01-specific T-cell responses in all patients. Competitive manufacturing time was demonstrated. EVX-01 was shown to be safe and able to elicit immune responses targeting tumor neoantigens with encouraging early indications of a clinical and meaningful antitumor efficacy, warranting further study.

Bacteria-host transcriptional response during endothelial invasion by Staphylococcus aureus.

Staphylococcus aureus is the cause of serious vascular infections such as sepsis and endocarditis. These infections are notoriously difficult to treat, and it is believed that the ability of S. aureus to invade endothelial cells and persist intracellularly is a key mechanism for persistence despite ongoing antibiotic treatment. Here, we used dual RNA sequencing to study the simultaneous transcriptional response of S. aureus and human endothelial cells during in vitro infections. We revealed discrete and shared differentially expressed genes for both host and pathogen at the different stages of infection. While the endothelial cells upregulated genes involved in interferon signalling and antigen presentation during late infection, S. aureus downregulated toxin expression while upregulating genes related to iron scavenging. In conclusion, the presented data provide an important resource to facilitate functional investigations into host-pathogen interaction during S. aureus invasive infection and a basis for identifying novel drug target sites.

Endurance exercise training-responsive miR-19b-3p improves skeletal muscle glucose metabolism.

Skeletal muscle is a highly adaptable tissue and remodels in response to exercise training. Using short RNA sequencing, we determine the miRNA profile of skeletal muscle from healthy male volunteers before and after a 14-day aerobic exercise training regime. Among the exercise training-responsive miRNAs identified, miR-19b-3p was selected for further validation. Overexpression of miR-19b-3p in human skeletal muscle cells increases insulin signaling, glucose uptake, and maximal oxygen consumption, recapitulating the adaptive response to aerobic exercise training. Overexpression of miR-19b-3p in mouse flexor digitorum brevis muscle enhances contraction-induced glucose uptake, indicating that miR-19b-3p exerts control on exercise training-induced adaptations in skeletal muscle. Potential targets of miR-19b-3p that are reduced after aerobic exercise training include KIF13A, MAPK6, RNF11, and VPS37A. Amongst these, RNF11 silencing potentiates glucose uptake in human skeletal muscle cells. Collectively, we identify miR-19b-3p as an aerobic exercise training-induced miRNA that regulates skeletal muscle glucose metabolism.

Transcriptomic and epigenomics atlas of myotubes reveals insight into the circadian control of metabolism and development.

Aim: Innate circadian rhythms are critical for optimal tissue-specific functions, including skeletal muscle, a major insulin-sensitive tissue responsible for glucose homeostasis. We determined whether transcriptional oscillations are associated with CpG methylation changes in skeletal muscle. Materials & methods: We performed rhythmicity analysis on the transcriptome and CpG methylome of circadian synchronized myotubes. Results: We identified several transcripts and CpG-sites displaying oscillatory behavior, which were enriched with Gene Ontology terms related to metabolism and development. Oscillating CpG methylation was associated with rhythmic expression of 31 transcripts. Conclusion: Although circadian oscillations may be regulated by rhythmic DNA methylation, strong rhythmic associations between transcriptome and CpG methylation were not identified. This resource constitutes a transcriptomic/epigenomic atlas of skeletal muscle and regulation of circadian rhythms.

Thermostability profiling of MHC-bound peptides: a new dimension in immunopeptidomics and aid for immunotherapy design.

The features of peptide antigens that contribute to their immunogenicity are not well understood. Although the stability of peptide-MHC (pMHC) is known to be important, current assays assess this interaction only for peptides in isolation and not in the context of natural antigen processing and presentation. Here, we present a method that provides a comprehensive and unbiased measure of pMHC stability for thousands of individual ligands detected simultaneously by mass spectrometry (MS). The method allows rapid assessment of intra-allelic and inter-allelic differences in pMHC stability and reveals profiles of stability that are broader than previously appreciated. The additional dimensionality of the data facilitated the training of a model which improves the prediction of peptide immunogenicity, specifically of cancer neoepitopes. This assay can be applied to any cells bearing MHC or MHC-like molecules, offering insight into not only the endogenous immunopeptidome, but also that of neoepitopes and pathogen-derived sequences.

Exercise training alters the genomic response to acute exercise in human adipose tissue.

AIM: To determine the genomic mechanisms by which adipose tissue responds to acute and chronic exercise. METHODS: We profiled the transcriptomic and epigenetic response to acute exercise in human adipose tissue collected before and after endurance training. RESULTS: Although acute exercises were performed at same relative intensities, the magnitude of transcriptomic changes after acute exercise was reduced by endurance training. DNA methylation remodeling induced by acute exercise was more prominent in trained versus untrained state. We found an overlap between gene expression and DNA methylation changes after acute exercise for 32 genes pre-training and six post-training, notably at adipocyte-specific genes. CONCLUSION: Training status differentially affects the epigenetic and transcriptomic response to acute exercise in human adipose tissue.

RAIN: RNA-protein Association and Interaction Networks.

Protein association networks can be inferred from a range of resources including experimental data, literature mining and computational predictions. These types of evidence are emerging for non-coding RNAs (ncRNAs) as well. However, integration of ncRNAs into protein association networks is challenging due to data heterogeneity. Here, we present a database of ncRNA-RNA and ncRNA-protein interactions and its integration with the STRING database of protein-protein interactions. These ncRNA associations cover four organisms and have been established from curated examples, experimental data, interaction predictions and automatic literature mining. RAIN uses an integrative scoring scheme to assign a confidence score to each interaction. We demonstrate that RAIN outperforms the underlying microRNA-target predictions in inferring ncRNA interactions. RAIN can be operated through an easily accessible web interface and all interaction data can be downloaded.Database URL: http://rth.dk/resources/rain.

Muscle Contraction Induces Acute Hydroxymethylation of the Exercise-Responsive Gene Nr4a3.

Exercise training triggers numerous positive adaptations through the regulation of genes controlling muscle structure and function. Epigenetic modifications, including DNA methylation, participate in transcriptional activation by allowing the recruitment of the transcription machinery to gene promoters. Exercise induces dynamic DNA demethylation at gene promoters; however, the contribution of the demethylation precursor hydroxymethylcytosine is unknown. Given the evanescent nature of hydroxymethylcytosine, a muscle contraction model that allows for the collection of samples that are repeatedly stimulated over time is required to determine whether contraction-induced demethylation is preceded by changes in the hydroxymethylcytosine level. Here, we established an acute skeletal muscle contraction model to mimic the effects of acute exercise on gene expression. We used this model to investigate the effect of muscle contraction on DNA demethylation and hydroxymethylation. First, we performed an acute exercise study in healthy humans to identify an exercise-responsive gene that we could study in culture. We identified the nuclear receptor subfamily 4 group A member 3 (Nr4a3) gene with the highest fold-expression increase after acute exercise. We then refined an electrical pulse stimulation (EPS) protocol that could induce expression of the Nr4a3 gene in C2C12 myotubes. Using targeted bisulfite sequencing, we found that in response to EPS, a region of the Nr4a3 promoter is rapidly demethylated at 60 min and re-methylated at 120 min. Of interest, hydroxymethylation of the differentially methylated region of Nr4a3 promoter after EPS was elevated immediately after EPS, with lowest levels reached at 60 min after EPS. In conclusion, we have established a cell culture-based protocol to mimic the acute transcriptional responses to exercise. Furthermore, we provide insight into the mechanism by which the exercise-responsive gene Nr4a3 is demethylated after muscle contraction.

High-resolution kinetics and modeling of hydrogen peroxide degradation in live cells.

Although the role of oxidative stress factors and their regulation is well studied, the temporal dynamics of stress recovery is still poorly understood. In particular, measuring the kinetics of stress recovery in the first minutes after acute exposure provides a powerful technique for assessing the role of regulatory proteins or enzymes through the use of mutant backgrounds. This project endeavors to screen the temporal dynamics of intracellular oxidant levels in live cells as a function of gene deletion in the budding yeast, Saccharomyces cerevisiae. Using the detailed time dynamics of extra- and intra-cellular peroxide we have developed a mathematical model that describes two distinct kinetic processes, an initial rapid degradation in the first 10-20min followed by a slower process. Using this model, a qualitative comparison allowed us to assign the dependence of temporal events to genetic factors. Surprisingly, we found that the deletion of transcription factors Yap1p or Skn7p was sufficient to disrupt the establishment of the second degradation phase but not the initial phase. A better fundamental understanding of the role protective factors play in the recovery from oxidative stress may lead to strategies for protecting or sensitizing cell to this stress.

Ionizing Radiation Potentiates High-Fat Diet-Induced Insulin Resistance and Reprograms Skeletal Muscle and Adipose Progenitor Cells.

Exposure to ionizing radiation increases the risk of chronic metabolic disorders such as insulin resistance and type 2 diabetes later in life. We hypothesized that irradiation reprograms the epigenome of metabolic progenitor cells, which could account for impaired metabolism after cancer treatment. C57Bl/6 mice were treated with a single dose of irradiation and subjected to high-fat diet (HFD). RNA sequencing and reduced representation bisulfite sequencing were used to create transcriptomic and epigenomic profiles of preadipocytes and skeletal muscle satellite cells collected from irradiated mice. Mice subjected to total body irradiation showed alterations in glucose metabolism and, when challenged with HFD, marked hyperinsulinemia. Insulin signaling was chronically disrupted in skeletal muscle and adipose progenitor cells collected from irradiated mice and differentiated in culture. Epigenomic profiling of skeletal muscle and adipose progenitor cells from irradiated animals revealed substantial DNA methylation changes, notably for genes regulating the cell cycle, glucose/lipid metabolism, and expression of epigenetic modifiers. Our results show that total body irradiation alters intracellular signaling and epigenetic pathways regulating cell proliferation and differentiation of skeletal muscle and adipose progenitor cells and provide a possible mechanism by which irradiation used in cancer treatment increases the risk for metabolic disease later in life.

Tumor suppressor microRNAs are downregulated in myelodysplastic syndrome with spliceosome mutations.

Spliceosome mutations are frequently observed in patients with myelodysplastic syndromes (MDS). However, it is largely unknown how these mutations contribute to the disease. MicroRNAs (miRNAs) are small noncoding RNAs, which have been implicated in most human cancers due to their role in post transcriptional gene regulation. The aim of this study was to analyze the impact of spliceosome mutations on the expression of miRNAs in a cohort of 34 MDS patients. In total, the expression of 76 miRNAs, including mirtrons and splice site overlapping miRNAs, was accurately quantified using reverse transcriptase quantitative PCR. The majority of the studied miRNAs have previously been implicated in MDS. Stably expressed miRNA genes for normalization of the data were identified using GeNorm and NormFinder algorithms. High-resolution melting assays covering all mutational hotspots within SF3B1, SRSF2, and U2AF1 (U2AF35) were developed, and all detected mutations were confirmed by Sanger sequencing. Overall, canonical miRNAs were downregulated in spliceosome mutated samples compared to wild-type (P = 0.002), and samples from spliceosome mutated patients clustered together in hierarchical cluster analyses. Among the most downregulated miRNAs were several tumor-suppressor miRNAs, including several let-7 family members, miR-423, and miR-103a. Finally, we observed that the predicted targets of the most downregulated miRNAs were involved in apoptosis, hematopoiesis, and acute myeloid leukemia among other cancer- and metabolic pathways. Our data indicate that spliceosome mutations may play an important role in MDS pathophysiology by affecting the expression of tumor suppressor miRNA genes involved in the development and progression of MDS.