Next Article in Journal
Exploration of Preservation Methods for Utilizing Porcine Fetal-Organ-Derived Cells in Regenerative Medicine Research
Next Article in Special Issue
KRAS Mutation Subtypes and Their Association with Other Driver Mutations in Oncogenic Pathways
Previous Article in Journal
The Impact of Lens Epithelium-Derived Growth Factor p75 Dimerization on Its Tethering Function
Previous Article in Special Issue
TGFβ and Hippo Signaling Pathways Coordinate to Promote Acinar to Ductal Metaplasia in Human Pancreas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Significant Variations in Double-Stranded RNA Levels in Cultured Skin Cells

1
Biosciences Institute, Newcastle University, Newcastle upon Tyne NE2 4HH, UK
2
Fallujah College of Medicine, University of Fallujah, Al-Fallujah 31002, Iraq
3
Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen 40002, Thailand
4
Center of Excellence in Systems Microbiology, Faculty of Medicine, Chulalongkorn University, Bangkok 10330, Thailand
5
College of Medicine, University of Misan, Al-Sader Teaching Hospital, Amarah 62001, Iraq
6
Bioinformatics Support Unit, Medical School, Newcastle University, Newcastle upon Tyne NE2 4HH, UK
*
Author to whom correspondence should be addressed.
Cells 2024, 13(3), 226; https://doi.org/10.3390/cells13030226
Submission received: 27 November 2023 / Revised: 18 January 2024 / Accepted: 22 January 2024 / Published: 25 January 2024
(This article belongs to the Special Issue Molecular and Cellular Underpinnings of Cancer Vulnerability)

Abstract

:
Endogenous double-stranded RNA has emerged as a potent stimulator of innate immunity. Under physiological conditions, endogenous dsRNA is maintained in the cell nucleus or the mitochondria; however, if protective mechanisms are breached, it leaches into the cytoplasm and triggers immune signaling pathways. Ectopic activation of innate immune pathways is associated with various diseases and senescence and can trigger apoptosis. Hereby, the level of cytoplasmic dsRNA is crucial. We have enriched dsRNA from two melanoma cell lines and primary dermal fibroblasts, including a competing probe, and analyzed the dsRNA transcriptome using RNA sequencing. There was a striking difference in read counts between the cell lines and the primary cells, and the effect was confirmed by northern blotting and immunocytochemistry. Both mitochondria (10–20%) and nuclear transcription (80–90%) contributed significantly to the dsRNA transcriptome. The mitochondrial contribution was lower in the cancer cells compared to fibroblasts. The expression of different transposable element families was comparable, suggesting a general up-regulation of transposable element expression rather than stimulation of a specific sub-family. Sequencing of the input control revealed minor differences in dsRNA processing pathways with an upregulation of oligoadenylate synthase and RNP125 that negatively regulates the dsRNA sensors RIG1 and MDA5. Moreover, RT-qPCR, Western blotting, and immunocytochemistry confirmed the relatively minor adaptations to the hugely different dsRNA levels. As a consequence, these transformed cell lines are potentially less tolerant to interventions that increase the formation of endogenous dsRNA.

1. Introduction

The cycle of replication for most viruses involves double-stranded RNA (dsRNA) structures. Hence, the occurrence of cytoplasmic dsRNA is associated with viral infection and may trigger an innate immune response [1,2]. On the other hand, dsRNA is also produced from mitochondrial DNA and repetitive elements in nuclear DNA [2,3]. Mitochondrial DNA is transcribed in both directions, giving rise to complementary transcripts. Degradation of the noncoding sequences by SUV3 (mitochondrial ATP-dependent RNA helicase SUV3) and PNPase (Polynucleotide phosphorylase) reduces the levels of dsRNA and contains the remaining duplexes in the mitochondrial matrix. Functional impairment of SUV3 or PNPase due to mutations, drugs, or mitochondrial stress leads to an accumulation of dsRNA and leakage into the cytoplasm [3,4,5]. In the nucleus, dsRNA is predominantly formed from repetitive elements such as short and long interspersed elements (SINEs and LINEs) as well as endogenous retroviruses (ERVs) [2,6]. The transcripts from inverted repeats can fold back and form stem loops; alternatively, LINEs and ERVs can become bi-directionally transcribed and form hybrids in cis and -because of high copy numbers- in trans. Under physiological circumstances, transcription of the low complexity regions is minor with spurious dsRNA edited by ADAR1 p110 and contained in the nucleus. Increased transcription of repetitive elements in response to -for example- DNA de-methylation or a knockdown of ADAR1 leads to the export of dsRNA molecules into the cytoplasm [7].
DsRNA in the cytoplasm is recognized by various dsRNA-binding proteins that recognize structural differences between viral and endogenous transcripts to modulate a potential immune response [8]. The protein family of retinoic-acid-inducible gene-like receptors (RLRs) includes RIG1(retinoic-acid-inducible gene), MDA5 or IFIH1 (Interferon Induced with Helicase C Domain 1) and LGP-2/DHX58 (DExH-Box Helicase 58). Both RIG1 and MDA5 oligomerize upon dsRNA binding and signal via MAVS (Mitochondrial Antiviral Signaling Protein) to elicit an NF-κB (nuclear factor-kappa B)/IRF3/7 (interferon regulatory factors 3/7) mediated response. RIG1 recognizes short RNA hybrids with 5′ phosphate groups (up to ~300 bp), whereas MDA5 prefers long dsRNAs. Uncapped 5′ ends and extended perfectly complementary RNA duplexes are viral structural hallmarks; hence, RIG1/MDA5 predominantly recognize exogenous rather than endogenous dsRNA [2]. However, the discrimination can become blurred with increasing levels of cytosolic, endogenous dsRNA. LPG-2 displays regulatory effects rather than immune stimulation via interaction with dsRNA [9]. PKR (Protein kinase regulated by RNA, also known as eukaryotic translation initiation factor 2-alpha kinase 2, EIF2AK2), is an interferon-inducible protein that becomes activated by binding to dsRNA [10]. Binding of dsRNA induces dimerization of PKR and auto-phosphorylation. The activation leads to dissociation from the dsRNA and phosphorylation of the eukaryotic translation initiation factor 2 to inhibit global translation [11,12]. In addition, PKR activation leads to a type 1 interferon response via IRF3/7.
Endogenous dsRNA is implicated in various human diseases, potentially triggering an underlying inflammatory phenotype. For example, PKR is upregulated in Alzheimer’s disease [13], and its inhibition improves memory and cognition [14]. Conversely, PKR acts as a tumor suppressor and prevents oncogenesis via inhibition of translation [15]. Moreover, mutated RIG1 and MDA5 are associated with defects in innate immunity, confirming a tightly controlled balance between endogenous dsRNA formation and defense mechanisms [16].
Bidirectional transcription of endogenous retrovirus (ERV) stimulates endogenous dsRNA formation and interferon response, as shown in 63 cancer cell lines [17,18]. Since dsRNA-activated innate immunity may lead to apoptosis, hypomethylating drugs that promote dsRNA formation are widely used to treat hematological cancers and potentially other types of solid tissue tumors [19]. Similarly, induction of endogenous retrovirus transcription via HDAC inhibitors or demethylating drugs followed by dsRNA formation is harnessed to improve the action of checkpoint inhibitors [20].
A critical parameter that potentially shapes a dsRNA-mediated response to therapeutic or experimental interventions are endogenous levels of dsRNA and the related dsRNA sensor proteins [21]. Recent investigations into the dsRNA transcriptome of mouse liver and testis suggested significant organ-specific differences in both the amount and origin of dsRNA [22]. Given the importance of cell lines in establishing paradigms of cancer biology, we characterized the dsRNA transcriptome and the expression of sensor proteins in selected cellular models and found striking differences.

2. Materials and Methods

2.1. Cell Culture

A375 melanoma cells were kindly provided by M. Panagiotidis (Northumbria University, Newcastle, UK), and C8161 cells and fibroblasts were from PE Lovat (Newcastle University, Newcastle, UK). The cells were cultured in a DMEM high glucose medium (Sigma, Gillingham, UK) supplemented with 10% fetal calf serum (Sigma), 2 mM L-glutamine (Sigma), and 5 units/0.5 mg/mL penicillin/streptomycin (Sigma). Cells were grown at 37 °C in 5% CO2, and the medium was changed every 3 days. Passaging was performed when the cells reached 80–90% confluence. Cells were seeded in 6 well plates at a density of 250–500 × 103 cells per well, in 96 well plates at 5–10 × 103 cells, or in 50 mm dishes at 1 × 106 cells per dish. For immunocytochemistry, cells were seeded on coverslips in 12 well plates at a density of 80 × 103 cells per well.

2.2. RNA Extraction and RT-qPCR

RNA was extracted using the RNeasy kit (Qiagen, Manchester, UK) following the supplier’s instructions. The RNA was eluted from the columns with water and quantified using Nanodrop (Thermo Scientific, Paisley, UK). Samples were tested on 2% agarose gels in TAE buffer stained with Gel Green to ensure RNA integrity. For ᴄDNA synthesis, the Omniscript kit (Qiagen) with random hexamers (Qiagen) was used. Each 20 μL reaction contained 2 μg of total RNA in nuclease-free water up to 13 μL, 2 μL of dNTP (10 mM), 2 μL of 10× buffer, and 2 μL random hexamers (0.4 μg/mL) and 1 μL of Omniscript reverse transcriptase (200 U). Reactions were incubated at 37 °C for 1 h.
The expression profiling was performed by qPCR using a light cycler 480 (Roche, Burgess Hill, UK). Reactions were run in 96 well plates, each containing 0.5 μL of forward and reverse primers (10 μM), 1 μL of 1:5 diluted ᴄDNA, 3 μL of nuclease-free water, and 5 μL SYBR Green mix (Roche). The total reaction volume was 10 μL, and each condition was performed in three replicates. Actin was used as a reference gene. The cycling protocol was as follows: Activation of the polymerase at 95 °C for 10 min, followed by 45 cycles of 10 s at 95 °C, 20 s at 58 °C, and 5 s at 72 °C. Melting curve analysis was run to confirm the product. Primers were either designed using web-based tools or adapted from published studies and listed in Supplementary Table S1.

2.3. Western Blotting

Cells were grown in 50 mm Petri dishes followed by lysis in RIPA buffer (150 mM NaCl, 1.0% IGEPAL® CA-630, 0.5% sodium deoxycholate, 0.1% SDS, 50 mM Tris, pH 8.0) and protein quantification was performed using a BCA protein assay (Pierce, Paisley, UK). 30 μg of protein lysate were run on CriterionXT precast 3–8% gradient gels (Bio-Rad, Watford, UK) with a protein ladder (Bio-Rad) in 1× tricine acetate running buffer (Bio-Rad). Proteins were transferred to PVDF membranes (Bio-Rad) using the Trans-Blot Turbo system (Bio-Rad). Membranes were washed in PBS-T (1× PBS, 0.1% Tween) and then blocked in 5% (w/v) skimmed milk in TBS-T for 1 h at room temperature. Primary antibodies against PKR (Abcam, ab32052, Cambridge, UK), RIG1 (Abcam, ab180675), MDA5 (Abcam, ab126630), ADAR (Atlas, B115763, Stockholm, Sweden) were used at a dilution of 1:5000 and anti-actin (Rabbit, Sigma; A2066) at 1:10,000. All antibodies were diluted in a blocking solution and incubated with the membranes at 4 °C overnight. After two washes in TBS-T, the secondary antibody (goat anti-rabbit IgG-HRP, 4030-05, 1:5000, Southern Biotech, Upper Heyford, UK ) in blocking buffer was applied at room temperature for 1 h. After washes in PBS-T, ECL (Bio-Rad) was added, followed by imaging using a chemiluminescence documentation system.

2.4. Immunocytochemistry

Cells were grown on coverslips (13 mm), washed with 1 mL PBS, and fixed for 10 min in 4% paraformaldehyde in PBS (Affymetrix, High Wycombe, UK). Cells were washed three times with PBS-T for 5 min, followed by permeabilization with 0.25% Triton-X-100 (Sigma) in PBS for 15 min. Samples were washed three times in PBS-T and blocked with ready-made 10% normal goat serum (Life Technologies, Paisley, UK) for 1 h at room temperature. Primary antibodies were diluted in blocking solution: J2 (monoclonal mouse IgG2a kappa chain; Scicons, Susteren, The Netherlands) and TOM20 (Abcam, ab186734) at 1:200 dilution, PKR (monoclonal IgG rabbit anti-human, Abcam, ab32052), MDA5 (monoclonal IgG rabbit anti-human, Abcam, ab126630) at 1:1000. After incubation at RT for 1 h cells were washed (three times 5 min in PBS-T) followed by secondary antibody incubation: For J2, an Alexa fluor™ 488 coupled goat anti-mouse IgG2a (Thermo) was used, otherwise an Alexa fluor™ 594 goat anti-rabbit IgG H + L (Thermo). Both secondary antibodies were diluted in blocking buffer 1:1000 and incubated for 1 hr at RT in the dark. After three washes, the cells were stained with DAPI (Vector Laboratories, Newark, CA, USA) and mounted on slides with ProlongTM Glass Antifade hard mounting medium (Thermo). Cells were imaged using an LSM800 Airyscan confocal microscope (Zeiss, Jena, Germany).

2.5. Spike in Probe

Plasmids encoding the natural sense–antisense transcript pair (Slc34a2a and Slc34a2aas from zebrafish) [23] were linearized with XbaI and transcribed in vitro using the MEGAscript T7 Transcription Kit (Invitrogen, Paisley, UK). The transcripts are 2607 bases (sense, NM_131624) and 1371 bases (antisense, NR_002876.2) long and share 563 bp of complementarity over three exons. The resulting RNA was quantified and mixed in equimolar concentrations to a total concentration of 0.4 μg/μL. One microliter was diluted 500× with 0.1 M NaCl, heated to 70 °C, and gradually cooled to hybridize the two strands. One microliter of the spike-in probe (0.8 ng) was added to the testis homogenate before J2 binding (see below).

2.6. Double-Stranded RNA Immunopurification

In essence, a published protocol was followed [3,22]. Cells were scraped and homogenized in 220 µL of NP-40 lysis buffer (50 mM Tris pH 7.5, 150 mM NaCl, 5 mM EDTA, 1% NP-40, 0.5% Na-deoxycholine, 220 units RNasin) using a syringe with 25G needle. 10% of the homogenate was used as input control, and total RNA was isolated using Qiazol according to the supplier’s protocol. The volume of the remaining sample was adjusted to 1 mL per sample with NET2 + DOC buffer (50 mM Tris pH 7.5, 150 mM NaCl, 1 mM MgCl2, 0.5% Na-deoxycholine, 0.2 units/µL DNase I). J2 antibody (5 µg, Scicons) and the spike in the probe were added. Samples were rotated for 3 h at 4 °C, and then 100 µL of µMACS Protein G MicroBeads (Miltenyi Biotec, Tokyo, Japan) were added and incubated for 1 h. µMACS columns were equilibrated with NP-40 buffer (50 mM Tris pH 7.5, 150 mM NaCl, 5 mM EDTA, 1% NP-40) followed by loading and washing the samples with 300 µL of NP-40 and 3 × 250 µL of wash buffer (50 mM Tris pH 7.5, 1 M NaCl, 1 mM EDTA, 1% NP-40, 0.5% Na-deoxycholine, and 0.1% SDS). DsRNA was eluted with hot water and purified using Qiazol (Qiagen). RNA was reverse transcribed using the SMARTer Stranded Total RNA-Seq Kit v3—Pico Input (Takara, Tokyo, Japan) and sequenced on a NextSeq 550 System (Illumina, San Diego, CA, USA).

2.7. Data Analysis

The quality of reads was assessed using FastQC, https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 11 March 2022), and adaptors were trimmed using Trimmomatic (version 0.3.6) [24]. The spike probe reads in the different samples were aligned to a fragment of zebrafish Chromosome 1 (Chr 1: 14,432,434–14,454,662) that contains the Slc34a2a gene and the related natural antisense transcript (Slc34a2aas) originating from the bidirectional Rbpja promoter using STAR version 2.5.2b [25]. All data sets were then quantified using Salmon [26], and expression differences between the input samples of the cell lines were established using DESeq2 [27]. To establish the dsRNA transcriptome reads from all samples were mapped to the reference genome (GRCm38.p5) using STAR. Peaks and coverage were established with BEDTools genomecov (-bg -ibam) and multicov (BEDTools suite version 2.28.0; [28]). Regions with coverage greater than or equal to five times more than the background were annotated using ChIPpeakAnno (version 3.19.4) [29], Figure 1.
SQuIRE (Software for Quantifying Interspersed Repeat Expression, version 0.9.9.9a-beta) is designed to quantify transposable element (TE) expression at the subfamily and locus level [30]. The required tools (STAR 2.5.3a, BEDTools 2.25.0, SAMTools 1.1, stringtie 1.3.3, DESeq2 1.16.1, R 3.4.1 and Python 2.7) were installed using Anaconda 4.12.0. SQuIRE Fetch was used to download repeat masker annotation files, chromosome information files, and the STAR index using the following command (SQuIRE Fetch --build hg38 --fasta --rmsk --chrom_info --gene -v). This step was followed by aligning the RNA-seq data to the reference genome using SQuIRE Map (SQuIRE Map -1 R1.fastq.gz -2 R2.fastq.gz -r 141 -n file_name -b hg38, where r equals the maximum read length). SQuIRE Count was then used with the following options (SQuIRE Count -r 141 -n file.bam -s 1 -b hg38, where –s is the library format). Differential gene expression analysis was performed using SQuIRE Call (SQuIRE Call -1 treatment.file -2 control.file -A treatment -B Control -o output -i SQuIRE_count_folder). The fastq files have been deposited to the Sequence Read Archive (SRA), bioproject PRJNA1053654.

3. Results

We used the established protocol, including the antibody J2, to enrich dsRNA under non-denaturing conditions from the melanoma cell lines A375 and C8161, as well as from skin fibroblasts [3,22]. To control for enrichment, in vitro synthesized dsRNA with stretches of perfect complementarity, bulges, and single-stranded ends were included. The recovered RNA, including an input control, was reverse transcribed and sequenced. After quality control, between 21–39 × 106 and 20–44 × 106, reads were obtained for dsRNA enriched samples and input RNA, respectively, of which 50–70% mapped to the human genome or the zebrafish probe. The amount of the added probe was in the low ng range per sample, controlling for dsRNA immune enrichment with J2 and serving as a carrier for samples with low intrinsic dsRNA. For unknown reasons, the immune enrichment for one fibroblast sample failed and was excluded from the analysis. A summary of the read statistics is given in Supplementary Table S2.
The alignment of the spike in reads confirmed the validity of the enrichment protocol and also underpinned the observation that complementary regions are processed with reduced efficiency during library preparation and sequencing (Figure 2). The significantly smaller number of reads mapping to the antisense strand—one would expect sense-antisense to be roughly equimolar—may indicate that one antisense strand binds to several sense-strand molecules or reduced stability of the antisense transcript.
There was a striking discrepancy in read numbers between A375 cells, C8161 cells, and fibroblasts, with the cancerous cell lines expressing significantly higher levels of dsRNA than the fibroblasts. In fact, dsRNA levels in fibroblasts were so low that the spike probe accounted for almost all isolated dsRNA (Supplementary Table S2). The low abundance of dsRNA in fibroblasts was confirmed using dot blot analysis and underpinned by immunocytochemistry, where no J2 staining was found (Figure 3). Hence, fibroblast dsRNA reads were excluded from the peak-calling analysis (see below), where a certain read density is required to establish the background and generate a meaningful number of peaks.
Reads were then mapped to the human genome, and expression was quantified using a published peak-calling pipeline to determine the presence of dsRNA-forming transcripts (Supplementary Figure S1) [22]. The dsRNA-associated reads mapped predominantly in promoter areas and in the 1000 bp downstream of genes. Fewer were mapped to exons, and less than 10% were each mapped to UTRs, introns, and intergenic regions. Hela cells [3] showed a different pattern, with most of the reads mapping to introns and intergenic regions (44.2% and 40.2%, respectively) followed by exons (7.1%) and less than 4% for the other biotypes (Supplementary Figure S1, upper panel). All input samples (A375, C8161, and fibroblasts) showed a comparable distribution of biotypes with read mapping predominantly to gene flanking regions (promoters, immediate downstream) and exons, whereas 5′ UTR, 3′ UTR, introns, and intergenic regions were associated with <7% of the reads.
Endogenous dsRNA is essentially produced from mitochondrial transcripts and repetitive regions of the genome. To characterize the dsRNA transcriptome further, reads from both dsRNA-enriched and input samples were aligned to the human mitochondrial genome; in addition, the contribution of repetitive elements to the transcriptome was analyzed using SQuIRE [30]. In all samples, approximately 10–15% of total mapped reads align to the mitochondrial genome (Figure 4). This proportion increased slightly (median 15 vs. 10.5%) in the dsRNA samples, though the change was not significant. The quantification and comparison of mitochondrial coverage between samples is problematic since the total number of mitochondrial DNA copies is unknown. Likewise, the dsRNA-enriched samples lack robust internal standards for normalization; consequently, the results need to be interpreted qualitatively.
The second common source of dsRNA in cells is repetitive DNA, i.e., SINEs, ERVs, and LINEs. Reads derived from repetitive DNA elements, both single-mapping and multi-mapping reads, were characterized and quantified using the software package SQuIRE [30]. The pipeline establishes an expression profile using uniquely mapping reads and assigns the multimappers accordingly to the relevant repeat elements. An arbitrary cut-off at 100 FPKM was set, and both TE expression and TE subfamilies of Alu elements, ERVs, and LINEs were determined. Normalization of these data is challenging because the dsRNA data sets lack parameters that could be used for normalization; hence, total read numbers were used. The input data of the different cell populations did not show striking differences in relative expression levels of the TEs, though trends become apparent. First, C8161 cells show higher relative levels of Alus, SINEs, and LINEs than A375 cells and fibroblasts (Figure 5A). This is also reflected by greater subfamily variety (Figure 5B). Moreover, both relative expression and variety of TEs are slightly enriched in the dsRNA-enriched samples, except for LINE subfamilies in fibroblasts, which show the opposite trend (Figure 5). This outlier is likely the result of low read numbers in the fibroblast samples, especially affecting values around the cut-off. The findings so far suggest that in all tested cells, a comparable repertoire of repetitive elements is transcribed, though at a strikingly different level. This raises the question whether the variations in dsRNA levels shape a specific cellular background regarding dsRNA sensors and innate immunity.
Cells may be exposed to dsRNA in the form of viruses or from endogenous sources, including repetitive elements and mitochondria. To mitigate dsRNA in the cytoplasm, cells express proteins that recognize RNA hybrids to either reduce their immunogenicity or trigger a signaling cascade activating innate immunity (Figure 6). The expression levels of these dsRNA sensors were tested using RT-qPCR, Western blotting (for practical reasons only in A375 and C8161 cells), immunocytochemistry, and targeted pathway analysis.
Levels of PKR, MDA5, and RIG1 mRNA were slightly but still significantly lower in A375 than in C8161 cells, whereas for ADAR1, no difference was found (Figure 6A). At the protein level, only MDA5 showed significant changes, though C8161 cells displayed higher expression. The same trend, higher transcript, and protein levels in C8161 were also apparent for PKR and RIG1, though the changes in protein were not significant (Figure 6B).
A375 cells are larger than C8161 cells and divide at a lower rate. Immunocytochemical staining with both cell lines revealed a punctate cytoplasmic pattern of dsRNA-related fluorescence (Figure 7, red staining). The dsRNA does not obviously co-localize with the sensor proteins ADAR1 and MDA5 (signal overlap between 0.2 and 0.4). RIG1 in both cell lines and PKR in A375 cells, however, showed co-localization with dsRNA, predominantly in a ring around the nucleus (R2 of around 0.6 or greater). Of note, A375 cells showed significantly higher co-localization of dsRNA with RIG1/PKR than C8161 cells, in line with the higher levels of dsRNA in A375 vs. C8161 cells. Staining of fibroblasts with the J2 antibody revealed no signal despite clear detection of all sensor proteins assessed (ADAR1, MDA5, PKR, pPKR, and also TOMM20 to visualize mitochondria), in line with the very low read counts in the RNAseq experiments (Figure 3C).
Next, we compared the expression values of dsRNA sensors and related proteins in the RNAseq input samples and performed a KEGG pathway analysis focussing on RIG1-like receptor signaling. Interestingly, 2′, 5′ oligoadenylate synthase (2′, 5′ OAS) is significantly upregulated in both melanoma cell lines as compared to fibroblasts (Figure 8A). 2′, 5′ OAS binds to dsRNA and activates RNaseL, leading to cleavage of viral and cellular single-stranded RNA. Sustained activation induces autophagy and apoptosis. In addition, RNaseL represses the transposition of LINE1 in cell lines and reduces LINE1 RNA [31]. The ring finger protein RNP125, an E3 ubiquitin-protein ligase, is also significantly upregulated in the two melanoma cell lines as compared to fibroblasts (Figure 8B). RNP125 ubiquinates both RIG1 and MDA5, which leads to their degradation and blunts innate immune signaling [32].
To conclude, the striking differences in dsRNA levels between fibroblasts and cancerous melanoma cell lines influence the expression pattern of dsRNA sensor proteins and potential regulators. However, the changes are relatively subtle.

4. Discussion

We report the dsRNA transcriptome of two melanoma cell lines and compare it to primary dermal fibroblasts; published data from HeLa cells were analyzed in parallel as a benchmark. Our results suggest that the cancerous cell lines express significantly higher levels of endogenous dsRNA than the primary cells. This discrepancy led to minor changes in dsRNA sensor proteins that eventually shape or suppress a cellular innate immune response.
The impact of intrinsic dsRNA on innate immune signaling has been widely established in the context of various diseases, most prominently in the context of autoimmune disease and cancer and recently also in cellular senescence [33,34,35,36]. The focus of these reports is on dsRNA-associated proteins and downstream signaling rather than the dsRNA itself. However, the level and nature of endogenous dsRNA will shape the various protective mechanisms that feed into the activation of innate immunity and, as a consequence, how a cell reacts to changes in dsRNA due to disease or drug interventions [16].
Cancerous cells show increased expression of both LINEs and SINEs due to transcriptional activation and copy number amplification [37]. On the other hand, mitochondrial dysfunction -another hallmark of cancer cells- promotes the leakage of dsRNA into the cytoplasm and activation of the dsRNA sensor proteins PKR and MDA5 [3,38]. Increased levels of dsRNA in the two investigated melanoma cell lines (A375 and C8161) compared to dermal fibroblasts could, therefore, be expected, though the scale of the difference is still striking. At this point, we cannot predict whether different dsRNA levels are attributable to variations between normal and cancer cells or between cell lines and primary cultured cells; comparable data are lacking. Should our observation apply to other cancer versus un-transformed cells, alternative options to enhance the efficacy of immunotherapy may open. Rather than using adjuvants (Azacitidine, for example) to increase the already significant levels of dsRNA in cancer cells, pathways that clear dsRNA (ADAR, PNPase, or autophagy [39]) could be targeted [40]. The fact that non-cancerous cells express much less dsRNA may make them less susceptible to the drug with potentially milder side effects. Both Azacitidine and inhibition of ADAR1 (or knockdown under experimental conditions) lead to an increase in cytosolic dsRNA and activation of the OAS/RNaseL pathway [41].
Comparative expression and pathway analysis have established increased expression of OAS in the two melanoma cell lines, an unexpected finding in light of increased levels of dsRNA in these cells. However, the OAS/RNaseL pathway also plays a crucial role in cellular RNA turnover, and enhanced expression could be due to roles unrelated to dsRNA [42]. Accordingly, OAS is upregulated in a large number of primary tumors and enhances resistance against therapy-induced DNA damage [43]. In contrast, the ubiquitin ligase RNP125, which is also upregulated in the melanoma cell lines, triggers the degradation of MDA5 and RIG1, making the cells more tolerant towards a high dsRNA load (Figure 8B).
The dsRNA transcriptome has been investigated in various cell lines and whole organs, though most of these efforts focused on relative changes in physiological vs. pathophysiological conditions, for example [44,45,46,47]. The approach we applied here, including a competing dsRNA probe, produced an estimate of net levels of dsRNA in the three different cell samples. We did not count the cells before immune enrichment, though the striking difference in read numbers and the lack of dsRNA signals on dot blots and immunocytochemistry validate the reported discrepancy between cell lines and primary dermal fibroblasts. We have investigated the source of the dsRNA and confirmed the important contribution of mitochondrial transcripts and repetitive elements. The relative scale of mitochondrial dsRNA expression remains to be determined since the number of mitochondria and mitochondrial genome copies may vary in the different cell models. We also attempted to determine the proportion of light to heavy strands in response to oxidative stress by RT-qPCR as a proxy for dsRNA formation. Despite evidence that oxidative stress leads to the release of mitochondrial dsRNA and activation of PKR [39], we were unable to establish consistent changes in the proportion of light to heavy strand RNA in response to H2O2 (Supplementary Figure S3).
Transcripts from repetitive elements were assessed using SQuIRE, and a trend of slight enrichment of TEs in the dsRNA samples was found. Transposon expression has been studied in detail in various cancers [48]. Their expression is greatly enhanced by histone deacetylases and DNA methyltransferase inhibitors, and many derive from unannotated, cryptic start sites [49,50]. Moreover, for acute myeloid leukemia, the ratio of TE [1] expression versus gene expression may have prognostic value for cancer progression. Of note, the low repeat cohort (low dsRNA) showed wide-ranging repression of immune response genes and poor prognosis. In contrast, the high-repeat patients (high dsRNA) showed reduced expression of tumor-related pathways and increased apoptosis [51]. Importantly, de-repression of transposons has been found in a wide range of diseases and is potentially responsible for the underlying chronic inflammatory phenotype [36,52,53,54]. TE expression and dsRNA also occur in pluripotent stem cells [55] and during hematopoietic development, though in the immune-privileged environment, other protective mechanisms appear to play a role, such as RNA methylation (m6A) [56].
The varied expression of dsRNA in cell lines did not cause strikingly different RNA and protein levels; the observed differences were mostly non-significant. These minor changes agree with the pathway analysis of differentially expressed genes between the three cell models (A375, C8161, and skin fibroblasts). None of the investigated dsRNA sensors was significantly changed; however, the ubiquitin ligase RNP125 that negatively regulates RIG1 and MDA5 was elevated in both transformed cell lines that have higher levels of dsRNA. Moreover, OAS is also upregulated to provide protection against high levels of dsRNA, which makes them highly responsive to increases in dsRNA, for example, induced by Azacitidine [41]. The finely balanced complexity of the dsRNA response is also documented by the significant expression changes of dsRNA sensors observed when cells are challenged with polyI:C or demethylating agents. The dsRNA stress caused significant expression changes depending on the nature of the cell lines and the particular stressors applied [57], thesis Dr Shaymaa Sadeq, Newcastle University). Evidence so far suggests that significant differences in dsRNA levels between individual tissues or specific cells within a tissue do exist [58] and may have physiological relevance. Though the nature of cells (healthy or pathological) or cell culture conditions could also have a determining effect on the dsRNA transcriptome, the underlying dsRNA sensors and the related immune signaling pathways, in essence, call for further research in the field.
Comprehensive characterization of the dsRNA transcriptome is hampered by significant technical difficulties. First, the immune enrichment depends essentially on lysis conditions. The breakdown of cellular compartments will mix and potentially hybridize RNAs that were not bona fide dsRNA molecules. Moreover, enzymatic processes, particularly reverse transcription, are negatively affected by dsRNA structures and may lead to a strand bias (Figure 2). Finally, bioinformatic analysis is challenging because of missing internal standards. For example, the quantification of mitochondrial reads is ambiguous because neither the number of mitochondria per cell nor the DNA molecules per mitochondria are usually known for a specific cell. Additionally, the quantification of nuclear reads is challenging if datasets from different biological systems were included. This may be achievable with a set of thoroughly tested spike-in probes added at the very start of the J2 enrichment procedure.

5. Conclusions

To conclude, we established a significant difference in dsRNA levels between two melanoma cell lines and primary dermal fibroblasts. However, the differences only lead to subtle changes in dsRNA sensing and signaling proteins. These findings may explain the discrepancies in susceptibility to dsRNA stress observed in different cell types. Our results highlight the importance of intrinsic dsRNA in shaping innate immunity and call for a comprehensive characterization of the dsRNA landscape in model systems such as cell lines, primary cells, or tissues. Moreover, experimental strategies that investigate the dsRNA transcriptome need to be standardized to allow for the comparison of different studies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells13030226/s1, Figure S1: Expression analysis using a peak annotation pipeline; Figure S2: Protein blots for Figure 6B; Figure S3: Mitochondrial dsRNA; Table S1: Primers for dsRNA sensors; Table S2: Statistics of RNA sequencing data. References [3,59,60,61,62] are cited in Supplementary Materials.

Author Contributions

Conceptualization, S.S., S.A.-H. and A.W.; methodology, S.S., S.A.-H., S.C. and S.R.; software, S.C., J.C. and S.A.-H.; formal analysis, S.S., S.A.-H., S.C., S.R. and A.W.; data curation, S.S., S.A.-H., S.C. and A.W.; writing—original draft preparation, A.W., S.S., S.A.-H. and S.C.; writing—review and editing, A.W., S.S., S.A.-H. and S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly funded by The Northern Counties Kidney Research Fund, Grant 18.011 to A.W., and the Iraqi Ministry of Higher Education (S.S. and S.A.-H.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

RNAseq data generated and discussed in this project have been submitted to SRA (accession number PRJNA1053654).

Acknowledgments

We would like to thank Penny Lovat (Newcastle University) for the C8161 cell line and dermal fibroblasts, as well as Mihalis Panagiotidis (Northumbria University) for A375 cells.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Chiappinelli, K.B.; Strissel, P.L.; Desrichard, A.; Li, H.; Henke, C.; Akman, B.; Hein, A.; Rote, N.S.; Cope, L.M.; Snyder, A.; et al. Inhibiting DNA Methylation Causes an Interferon Response in Cancer via dsRNA Including Endogenous Retroviruses. Cell 2015, 162, 974–986. [Google Scholar] [CrossRef]
  2. Hur, S. Double-Stranded RNA Sensors and Modulators in Innate Immunity. Annu. Rev. Immunol. 2019, 37, 349–375. [Google Scholar] [CrossRef] [PubMed]
  3. Dhir, A.; Dhir, S.; Borowski, L.S.; Jimenez, L.; Teitell, M.; Rotig, A.; Crow, Y.J.; Rice, G.I.; Duffy, D.; Tamby, C.; et al. Mitochondrial double-stranded RNA triggers antiviral signalling in humans. Nature 2018, 560, 238–242. [Google Scholar] [CrossRef] [PubMed]
  4. Kim, S.; Ku, Y.; Ku, J.; Kim, Y. Evidence of Aberrant Immune Response by Endogenous Double-Stranded RNAs: Attack from Within. Bioessays 2019, 41, e1900023. [Google Scholar] [CrossRef] [PubMed]
  5. Sadeq, S.; Al-Hashimi, S.; Cusack, C.M.; Werner, A. Endogenous Double-Stranded RNA. Non-Coding RNA 2021, 7, 15. [Google Scholar] [CrossRef]
  6. Chen, Y.G.; Hur, S. Cellular origins of dsRNA, their recognition and consequences. Nat. Rev. Mol. Cell Biol. 2022, 23, 286–301. [Google Scholar] [CrossRef]
  7. Datta, R.; Adamska, J.Z.; Bhate, A.; Li, J.B. A-to-I RNA editing by ADAR and its therapeutic applications: From viral infections to cancer immunotherapy. Wiley Interdiscip. Rev. RNA 2023, 15, e1817. [Google Scholar] [CrossRef]
  8. Wang, Q.; Carmichael, G.G. Effects of length and location on the cellular response to double-stranded RNA. Microbiol. Mol. Biol. Rev. 2004, 68, 432–452. [Google Scholar] [CrossRef]
  9. Bruns, A.M.; Pollpeter, D.; Hadizadeh, N.; Myong, S.; Marko, J.F.; Horvath, C.M. ATP hydrolysis enhances RNA recognition and antiviral signal transduction by the innate immune sensor, laboratory of genetics and physiology 2 (LGP2). J. Biol. Chem. 2013, 288, 938–946. [Google Scholar] [CrossRef]
  10. Pfaller, C.K.; Li, Z.; George, C.X.; Samuel, C.E. Protein kinase PKR and RNA adenosine deaminase ADAR1: New roles for old players as modulators of the interferon response. Curr. Opin. Immunol. 2011, 23, 573–582. [Google Scholar] [CrossRef]
  11. Lemaire, P.A.; Anderson, E.; Lary, J.; Cole, J.L. Mechanism of PKR Activation by dsRNA. J. Mol. Biol. 2008, 381, 351–360. [Google Scholar] [CrossRef]
  12. García, M.A.; Gil, J.; Ventoso, I.; Guerra, S.; Domingo, E.; Rivas, C.; Esteban, M. Impact of protein kinase PKR in cell biology: From antiviral to antiproliferative action. Microbiol. Mol. Biol. Rev. 2006, 70, 1032–1060. [Google Scholar] [CrossRef]
  13. Ohno, M. Roles of eIF2α kinases in the pathogenesis of Alzheimer’s disease. Front. Mol. Neurosci. 2014, 7, 22. [Google Scholar] [CrossRef]
  14. Zhu, P.J.; Huang, W.; Kalikulov, D.; Yoo, J.W.; Placzek, A.N.; Stoica, L.; Zhou, H.; Bell, J.C.; Friedlander, M.J.; Krnjević, K.; et al. Suppression of PKR promotes network excitability and enhanced cognition by interferon-γ-mediated disinhibition. Cell 2011, 147, 1384–1396. [Google Scholar] [CrossRef] [PubMed]
  15. Benkirane, M.; Neuveut, C.; Chun, R.F.; Smith, S.M.; Samuel, C.E.; Gatignol, A.; Jeang, K.T. Oncogenic potential of TAR RNA binding protein TRBP and its regulatory interaction with RNA-dependent protein kinase PKR. EMBO J. 1997, 16, 611–624. [Google Scholar] [CrossRef] [PubMed]
  16. de Reuver, R.; Maelfait, J. Novel insights into double-stranded RNA-mediated immunopathology. Nat. Rev. Immunol. 2023. ahead of print. [Google Scholar] [CrossRef] [PubMed]
  17. Atkins, M.B.; Lotze, M.T.; Dutcher, J.P.; Fisher, R.I.; Weiss, G.; Margolin, K.; Abrams, J.; Sznol, M.; Parkinson, D.; Hawkins, M.; et al. High-dose recombinant interleukin 2 therapy for patients with metastatic melanoma: Analysis of 270 patients treated between 1985 and 1993. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 1999, 17, 2105–2116. [Google Scholar] [CrossRef] [PubMed]
  18. Chen, L.; Li, Y.; Lin, C.H.; Chan, T.H.; Chow, R.K.; Song, Y.; Liu, M.; Yuan, Y.F.; Fu, L.; Kong, K.L.; et al. Recoding RNA editing of AZIN1 predisposes to hepatocellular carcinoma. Nat. Med. 2013, 19, 209–216. [Google Scholar] [CrossRef] [PubMed]
  19. Tsai, H.C.; Li, H.; Van Neste, L.; Cai, Y.; Robert, C.; Rassool, F.V.; Shin, J.J.; Harbom, K.M.; Beaty, R.; Pappou, E.; et al. Transient low doses of DNA-demethylating agents exert durable antitumor effects on hematological and epithelial tumor cells. Cancer Cell 2012, 21, 430–446. [Google Scholar] [CrossRef] [PubMed]
  20. Covre, A.; Coral, S.; Nicolay, H.; Parisi, G.; Fazio, C.; Colizzi, F.; Fratta, E.; Di Giacomo, A.M.; Sigalotti, L.; Natali, P.G.; et al. Antitumor activity of epigenetic immunomodulation combined with CTLA-4 blockade in syngeneic mouse models. Oncoimmunology 2015, 4, e1019978. [Google Scholar] [CrossRef]
  21. Wolpaw, A.J.; Grossmann, L.D.; Dessau, J.L.; Dong, M.M.; Aaron, B.J.; Brafford, P.A.; Volgina, D.; Pascual-Pasto, G.; Rodriguez-Garcia, A.; Uzun, Y.; et al. Epigenetic state determines inflammatory sensing in neuroblastoma. Proc. Natl. Acad. Sci. USA 2022, 119, e2102358119. [Google Scholar] [CrossRef]
  22. Werner, A.; Clark, J.E.; Samaranayake, C.; Casement, J.; Zinad, H.S.; Sadeq, S.; Al-Hashimi, S.; Smith, M.; Kotaja, N.; Mattick, J.S. Widespread formation of double-stranded RNAs in testis. Genome Res. 2021, 7, 1174–1186. [Google Scholar] [CrossRef] [PubMed]
  23. Nalbant, P.; Boehmer, C.; Dehmelt, L.; Wehner, F.; Werner, A. Functional characterization of a Na+-phosphate cotransporter (NaPi-II) from zebrafish and identification of related transcripts. J. Physiol. 1999, 520 Pt 1, 79–89. [Google Scholar] [CrossRef] [PubMed]
  24. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  25. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2012, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  26. Patro, R.; Duggal, G.; Love, M.I.; Irizarry, R.A.; Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 2017, 14, 417–419. [Google Scholar] [CrossRef] [PubMed]
  27. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed]
  28. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef]
  29. Zhu, L.J.; Gazin, C.; Lawson, N.D.; Pages, H.; Lin, S.M.; Lapointe, D.S.; Green, M.R. ChIPpeakAnno: A Bioconductor package to annotate ChIP-seq and ChIP-chip data. BMC Bioinform. 2010, 11, 237. [Google Scholar] [CrossRef]
  30. Yang, W.R.; Ardeljan, D.; Pacyna, C.N.; Payer, L.M.; Burns, K.H. SQuIRE reveals locus-specific regulation of interspersed repeat expression. Nucleic Acids Res. 2019, 47, e27. [Google Scholar] [CrossRef]
  31. Zhang, A.; Dong, B.; Doucet, A.J.; Moldovan, J.B.; Moran, J.V.; Silverman, R.H. RNase L restricts the mobility of engineered retrotransposons in cultured human cells. Nucleic Acids Res. 2013, 42, 3803–3820. [Google Scholar] [CrossRef] [PubMed]
  32. Arimoto, K.; Takahashi, H.; Hishiki, T.; Konishi, H.; Fujita, T.; Shimotohno, K. Negative regulation of the RIG-I signaling by the ubiquitin ligase RNF125. Proc. Natl. Acad. Sci. USA 2007, 104, 7500–7505. [Google Scholar] [CrossRef] [PubMed]
  33. Dunker, W.; Ye, X.; Zhao, Y.; Liu, L.; Richardson, A.; Karijolich, J. TDP-43 prevents endogenous RNAs from triggering a lethal RIG-I-dependent interferon response. Cell Rep. 2021, 35, 108976. [Google Scholar] [CrossRef] [PubMed]
  34. de Reuver, R.; Dierick, E.; Wiernicki, B.; Staes, K.; Seys, L.; De Meester, E.; Muyldermans, T.; Botzki, A.; Lambrecht, B.N.; Van Nieuwerburgh, F.; et al. ADAR1 interaction with Z-RNA promotes editing of endogenous double-stranded RNA and prevents MDA5-dependent immune activation. Cell Rep. 2021, 36, 109500. [Google Scholar] [CrossRef] [PubMed]
  35. Mullani, N.; Porozhan, Y.; Mangelinck, A.; Rachez, C.; Costallat, M.; Batsché, E.; Goodhardt, M.; Cenci, G.; Mann, C.; Muchardt, C. Reduced RNA turnover as a driver of cellular senescence. Life Sci. Alliance 2021, 4, e202000809. [Google Scholar] [CrossRef] [PubMed]
  36. Tossberg, J.T.; Heinrich, R.M.; Farley, V.M.; Crooke, P.S., 3rd; Aune, T.M. Adenosine-to-Inosine RNA Editing of Alu Double-Stranded (ds)RNAs Is Markedly Decreased in Multiple Sclerosis and Unedited Alu dsRNAs Are Potent Activators of Proinflammatory Transcriptional Responses. J. Immunol. 2020, 205, 2606–2617. [Google Scholar] [CrossRef]
  37. Gualtieri, A.; Andreola, F.; Sciamanna, I.; Sinibaldi-Vallebona, P.; Serafino, A.; Spadafora, C. Increased expression and copy number amplification of LINE-1 and SINE B1 retrotransposable elements in murine mammary carcinoma progression. Oncotarget 2013, 4, 1882–1893. [Google Scholar] [CrossRef]
  38. Kim, Y.; Park, J.; Kim, S.; Kim, M.; Kang, M.G.; Kwak, C.; Kang, M.; Kim, B.; Rhee, H.W.; Kim, V.N. PKR Senses Nuclear and Mitochondrial Signals by Interacting with Endogenous Double-Stranded RNAs. Mol. Cell 2018, 71, 1051–1063.e1056. [Google Scholar] [CrossRef]
  39. Kim, S.; Lee, K.; Choi, Y.S.; Ku, J.; Kim, H.; Kharbash, R.; Yoon, J.; Lee, Y.S.; Kim, J.H.; Lee, Y.J.; et al. Mitochondrial double-stranded RNAs govern the stress response in chondrocytes to promote osteoarthritis development. Cell Rep. 2022, 40, 111178. [Google Scholar] [CrossRef]
  40. Kung, C.P.; Cottrell, K.A.; Ryu, S.; Bramel, E.R.; Kladney, R.D.; Bao, E.A.; Freeman, E.C.; Sabloak, T.; Maggi, L., Jr.; Weber, J.D. Evaluating the therapeutic potential of ADAR1 inhibition for triple-negative breast cancer. Oncogene 2021, 40, 189–202. [Google Scholar] [CrossRef]
  41. Banerjee, S.; Gusho, E.; Gaughan, C.; Dong, B.; Gu, X.; Holvey-Bates, E.; Talukdar, M.; Li, Y.; Weiss, S.R.; Sicheri, F.; et al. OAS-RNase L innate immune pathway mediates the cytotoxicity of a DNA-demethylating drug. Proc. Natl. Acad. Sci. USA 2019, 116, 5071–5076. [Google Scholar] [CrossRef]
  42. Prangley, E.; Korennykh, A. 2-5A-Mediated decay (2-5AMD): From antiviral defense to control of host RNA. Crit. Rev. Biochem. Mol. Biol. 2022, 57, 477–491. [Google Scholar] [CrossRef] [PubMed]
  43. Weichselbaum, R.R.; Ishwaran, H.; Yoon, T.; Nuyten, D.S.; Baker, S.W.; Khodarev, N.; Su, A.W.; Shaikh, A.Y.; Roach, P.; Kreike, B.; et al. An interferon-related gene signature for DNA damage resistance is a predictive marker for chemotherapy and radiation for breast cancer. Proc. Natl. Acad. Sci. USA 2008, 105, 18490–18495. [Google Scholar] [CrossRef] [PubMed]
  44. Wiatrek, D.M.; Candela, M.E.; Sedmík, J.; Oppelt, J.; Keegan, L.P.; O’Connell, M.A. Activation of innate immunity by mitochondrial dsRNA in mouse cells lacking p53 protein. RNA 2019, 25, 713–726. [Google Scholar] [CrossRef] [PubMed]
  45. Shadle, S.C.; Bennett, S.R.; Wong, C.J.; Karreman, N.A.; Campbell, A.E.; van der Maarel, S.M.; Bass, B.L.; Tapscott, S.J. DUX4-induced bidirectional HSATII satellite repeat transcripts form intranuclear double-stranded RNA foci in human cell models of FSHD. Hum. Mol. Genet. 2019, 28, 3997–4011. [Google Scholar] [CrossRef] [PubMed]
  46. Arnaiz, E.; Miar, A.; Dias, A.G.; Prasad, N.; Schulze, U.; Waithe, D.; Rehwinkel, J.; Harris, A. Hypoxia regulates endogenous double-stranded RNA production via reduced mitochondrial DNA transcription. Front. Oncol. 2021, 11, 779739. [Google Scholar] [CrossRef]
  47. Gao, Y.; Chen, S.; Halene, S.; Tebaldi, T. Transcriptome-wide quantification of double-stranded RNAs in live mouse tissues by dsRIP-Seq. STAR Protoc. 2021, 2, 100366. [Google Scholar] [CrossRef]
  48. Burns, K.H. Transposable elements in cancer. Nat. Rev. Cancer 2017, 17, 415–424. [Google Scholar] [CrossRef]
  49. Brocks, D.; Schmidt, C.R.; Daskalakis, M.; Jang, H.S.; Shah, N.M.; Li, D.; Li, J.; Zhang, B.; Hou, Y.; Laudato, S.; et al. DNMT and HDAC inhibitors induce cryptic transcription start sites encoded in long terminal repeats. Nat. Genet. 2017, 49, 1052–1060. [Google Scholar] [CrossRef]
  50. Colombo, A.R.; Triche, T., Jr.; Ramsingh, G. Transposable Element Expression in Acute Myeloid Leukemia Transcriptome and Prognosis. Sci. Rep. 2018, 8, 16449. [Google Scholar] [CrossRef]
  51. Onishi-Seebacher, M.; Erikson, G.; Sawitzki, Z.; Ryan, D.; Greve, G.; Lübbert, M.; Jenuwein, T. Repeat to gene expression ratios in leukemic blast cells can stratify risk prediction in acute myeloid leukemia. BMC Med. Genom. 2021, 14, 166. [Google Scholar] [CrossRef]
  52. Rostami, M.R.; Bradic, M. The derepression of transposable elements in lung cells is associated with the inflammatory response and gene activation in idiopathic pulmonary fibrosis. Mob. DNA 2021, 12, 14. [Google Scholar] [CrossRef] [PubMed]
  53. Ochoa, E.; Ramirez, P.; Gonzalez, E.; De Mange, J.; Ray, W.J.; Bieniek, K.F.; Frost, B. Pathogenic tau-induced transposable element-derived dsRNA drives neuroinflammation. Sci. Adv. 2023, 9, eabq5423. [Google Scholar] [CrossRef] [PubMed]
  54. Dhillon, P.; Mulholland, K.A.; Hu, H.; Park, J.; Sheng, X.; Abedini, A.; Liu, H.; Vassalotti, A.; Wu, J.; Susztak, K. Increased levels of endogenous retroviruses trigger fibroinflammation and play a role in kidney disease development. Nat. Commun. 2023, 14, 559. [Google Scholar] [CrossRef] [PubMed]
  55. Fu, H.; Zhang, W.; Li, N.; Yang, J.; Ye, X.; Tian, C.; Lu, X.; Liu, L. Elevated retrotransposon activity and genomic instability in primed pluripotent stem cells. Genome Biol. 2021, 22, 201. [Google Scholar] [CrossRef]
  56. Gao, Y.; Vasic, R.; Song, Y.; Teng, R.; Liu, C.; Gbyli, R.; Biancon, G.; Nelakanti, R.; Lobben, K.; Kudo, E.; et al. m(6)A Modification Prevents Formation of Endogenous Double-Stranded RNAs and Deleterious Innate Immune Responses during Hematopoietic Development. Immunity 2020, 52, 1007–1021.e1008. [Google Scholar] [CrossRef]
  57. Guo, R.; Lu, S.Y.; Ma, J.X.; Wang, Q.L.; Zhang, L.; Tang, L.Y.; Shen, Y.; Shen, C.L.; Wang, J.J.; Lu, L.M.; et al. RIG-I acts as a tumor suppressor in melanoma via regulating the activation of the MKK/p38MAPK signaling pathway. Hum. Cell 2022, 35, 1071–1083. [Google Scholar] [CrossRef]
  58. Thomsen, C.; Røge, R.; Fred, Å.; Wanders, A. Immunohistochemical detection of double-stranded RNA in formalin-fixed paraffin-embedded tissue. APMIS 2023, 131, 197–205. [Google Scholar] [CrossRef]
  59. Zhang, N.N.; Shen, S.H.; Jiang, L.J.; Zhang, W.; Zhang, H.X.; Sun, Y.P.; Li, X.Y.; Huang, Q.H.; Ge, B.X.; Chen, S.J.; et al. RIG-I plays a critical role in negatively regulating granulocytic proliferation. Proc. Natl. Acad. Sci. USA 2008, 105, 10553–10558. [Google Scholar] [CrossRef]
  60. Nasirudeen, A.M.; Wong, H.H.; Thien, P.; Xu, S.; Lam, K.P.; Liu, D.X. RIG-I, MDA5 and TLR3 synergistically play an important role in restriction of dengue virus infection. PLoS Negl. Trop. Dis. 2011, 5, e926. [Google Scholar] [CrossRef]
  61. Wang, Q.; Hui, H.; Guo, Z.; Zhang, W.; Hu, Y.; He, T.; Tai, Y.; Peng, P.; Wang, L. ADAR1 regulates ARHGAP26 gene expression through RNA editing by disrupting miR-30b-3p and miR-573 binding. RNA 2013, 19, 1525–1536. [Google Scholar] [CrossRef] [PubMed]
  62. Strydom, E.; Pietersen, G. Development of a strand-specific RT-PCR to detect the positive sense replicative strand of Soybean blotchy mosaic virus. J. Virol. Methods 2018, 259, 39–44. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Bioinformatics analysis pipeline. The different samples are indicated at the top, and the software packages used are in blue. In addition to our own data sets from A375 cells, C8161 cells, and skin fibroblasts, previously published data from HeLa cells [3] was analyzed in parallel as an additional control. The zebrafish reference (*) encompasses bases 14,431,900–14,443,060 of chromosome 1 of the zebrafish genome assembly GRCz11.
Figure 1. Bioinformatics analysis pipeline. The different samples are indicated at the top, and the software packages used are in blue. In addition to our own data sets from A375 cells, C8161 cells, and skin fibroblasts, previously published data from HeLa cells [3] was analyzed in parallel as an additional control. The zebrafish reference (*) encompasses bases 14,431,900–14,443,060 of chromosome 1 of the zebrafish genome assembly GRCz11.
Cells 13 00226 g001
Figure 2. Reads mapping to the zebrafish Slc34a2a locus from the dsRNA spike in the probe. The scale for sense (Slc34a2a) is 290,000, and for the antisense (Slc34a2a(as)), it is 15,000.
Figure 2. Reads mapping to the zebrafish Slc34a2a locus from the dsRNA spike in the probe. The scale for sense (Slc34a2a) is 290,000, and for the antisense (Slc34a2a(as)), it is 15,000.
Cells 13 00226 g002
Figure 3. Fibroblasts express minimal amounts of dsRNA. DsRNA was assessed using the J2 antibody with dot blots and immunocytochemistry. (A) Three technical replicates of RNA were spotted on nitrocellulose probed with J2 and visualized using an HRP-coupled secondary antibody and CSPD star (Roche) detection. (B) Quantification of the luminescence signal in (A), three technical and two biological replicates (ANOVA with Tukey’s post hoc test). (C) Pictures 1–7, fibroblast cells; picture 8 (bottom right), A375 cells. The cells were fixed and stained with combinations of antibodies against dsRNA (J2, red) and dsRNA binding proteins (ADAR, RIG1, MDA5, PKR, and phosphoPKR) as well as TOMM20 (mitochondria) and DAPI (nucleus). In A375 control cells, the J2 signal is clearly visible throughout the cytoplasm, whereas it is undetectable in fibroblast.
Figure 3. Fibroblasts express minimal amounts of dsRNA. DsRNA was assessed using the J2 antibody with dot blots and immunocytochemistry. (A) Three technical replicates of RNA were spotted on nitrocellulose probed with J2 and visualized using an HRP-coupled secondary antibody and CSPD star (Roche) detection. (B) Quantification of the luminescence signal in (A), three technical and two biological replicates (ANOVA with Tukey’s post hoc test). (C) Pictures 1–7, fibroblast cells; picture 8 (bottom right), A375 cells. The cells were fixed and stained with combinations of antibodies against dsRNA (J2, red) and dsRNA binding proteins (ADAR, RIG1, MDA5, PKR, and phosphoPKR) as well as TOMM20 (mitochondria) and DAPI (nucleus). In A375 control cells, the J2 signal is clearly visible throughout the cytoplasm, whereas it is undetectable in fibroblast.
Cells 13 00226 g003
Figure 4. Summary of reads aligned to the mitochondrial genome. (A) Input and dsRNA reads of the three cell lines assessed, A375, C8161, and HeLa [3], as well as input from fibroblasts. (B) Total read numbers mapping to mitochondria, heavy (gene rich) versus light strand (gene poor). (C) Percentage of heavy (H) vs. light (L) strand mapping reads (t-test).
Figure 4. Summary of reads aligned to the mitochondrial genome. (A) Input and dsRNA reads of the three cell lines assessed, A375, C8161, and HeLa [3], as well as input from fibroblasts. (B) Total read numbers mapping to mitochondria, heavy (gene rich) versus light strand (gene poor). (C) Percentage of heavy (H) vs. light (L) strand mapping reads (t-test).
Cells 13 00226 g004
Figure 5. Expression of TE in A375, C8161, and fibroblast cells. (A) Total counts of normalized reads in FPKM mapping to Alu, ERV, and LINE elements. (B) Different TE families of Alu, ERV, and LINE elements in the three cell lines, A375, C8161, and fibroblasts. * p < 0.05; ** p < 0.01; ANOVA with Tukey’s post hoc test.
Figure 5. Expression of TE in A375, C8161, and fibroblast cells. (A) Total counts of normalized reads in FPKM mapping to Alu, ERV, and LINE elements. (B) Different TE families of Alu, ERV, and LINE elements in the three cell lines, A375, C8161, and fibroblasts. * p < 0.05; ** p < 0.01; ANOVA with Tukey’s post hoc test.
Cells 13 00226 g005
Figure 6. Double-stranded RNA sensors. (A) mRNA levels of dsRNA sensors in A375 and C8161 cells assessed using RT-qPCR. Actin was used as a reference. **** p < 0.001, ns, not significant, t-test. Of note, increased ΔCt values reflect lower mRNA expression levels. (B) Protein levels of dsRNA sensors assessed using Western blotting. ** p < 0.01, ns, not significant, t-test. Blots are provided in Supplementary Figure S2. (C) Schematic diagram of dsRNA triggered pathways. Left, PKR recognizes short perfect RNA hybrids, dimerizes, and undergoes autophosphorylation. Phospho PKR induces global translational arrest via phosphorylation of eIF2α and leads to NF-κB mediated interferon signaling. Middle, the two Rig-like receptors (RLR) RIG1and MDA5 bind dsRNA with different specificity; RIG1 binds 5′ phosphate groups and polymerizes along dsRNA, whereas MDA5 binds long (500–1000bp) RNA hybrids. Both activate MAVS and lead to NF-κB mediated immune signaling and apoptosis. Right, dsRNA from nuclear DNA becomes edited by ADAR1 and is prevented from exiting the nucleus. Should the barrier be breached, MDA5 and PKR may become activated.
Figure 6. Double-stranded RNA sensors. (A) mRNA levels of dsRNA sensors in A375 and C8161 cells assessed using RT-qPCR. Actin was used as a reference. **** p < 0.001, ns, not significant, t-test. Of note, increased ΔCt values reflect lower mRNA expression levels. (B) Protein levels of dsRNA sensors assessed using Western blotting. ** p < 0.01, ns, not significant, t-test. Blots are provided in Supplementary Figure S2. (C) Schematic diagram of dsRNA triggered pathways. Left, PKR recognizes short perfect RNA hybrids, dimerizes, and undergoes autophosphorylation. Phospho PKR induces global translational arrest via phosphorylation of eIF2α and leads to NF-κB mediated interferon signaling. Middle, the two Rig-like receptors (RLR) RIG1and MDA5 bind dsRNA with different specificity; RIG1 binds 5′ phosphate groups and polymerizes along dsRNA, whereas MDA5 binds long (500–1000bp) RNA hybrids. Both activate MAVS and lead to NF-κB mediated immune signaling and apoptosis. Right, dsRNA from nuclear DNA becomes edited by ADAR1 and is prevented from exiting the nucleus. Should the barrier be breached, MDA5 and PKR may become activated.
Cells 13 00226 g006
Figure 7. Immunohistochemical staining of dsRNA and dsRNA sensor proteins. The two cell lines, A375 and C8161, were assessed with the indicated antibodies against MDA5, RIG1, ADAR, PKR, and dsRNA (J2). The left panel shows nuclear DAPI staining (blue), followed by the specific dsRNA sensor proteins (green) and dsRNA (J2) in red. The panel on the right displays the composite. Scale bars are 10 µm (left panel). In the right panel, the spatial overlap of signals is compared between A375 and C8161 cells (R2). * p < 0.05, ** p < 0.01, **** p < 0.001, ns, not significant, t-test.
Figure 7. Immunohistochemical staining of dsRNA and dsRNA sensor proteins. The two cell lines, A375 and C8161, were assessed with the indicated antibodies against MDA5, RIG1, ADAR, PKR, and dsRNA (J2). The left panel shows nuclear DAPI staining (blue), followed by the specific dsRNA sensor proteins (green) and dsRNA (J2) in red. The panel on the right displays the composite. Scale bars are 10 µm (left panel). In the right panel, the spatial overlap of signals is compared between A375 and C8161 cells (R2). * p < 0.05, ** p < 0.01, **** p < 0.001, ns, not significant, t-test.
Cells 13 00226 g007
Figure 8. Expression differences in components of dsRNA signaling pathways between cell lines. (A) Expression differences as established using DESeq of selected dsRNA sensors, including MDA5, PKR, ADAR, RIG1, OAS, RNaseL, and MAVS. * adjusted p-value < 0.05; *** adjusted p-value < 0.005. (B) Components enriched in the KEGG pathway ‘Rig-like receptor signaling’ in C8161 cells versus fibroblasts and A375 versus fibroblasts. +P, phosphorylation; +U, ubiquitination; −U, deubiquitination. Adjusted p-value < 0.05.
Figure 8. Expression differences in components of dsRNA signaling pathways between cell lines. (A) Expression differences as established using DESeq of selected dsRNA sensors, including MDA5, PKR, ADAR, RIG1, OAS, RNaseL, and MAVS. * adjusted p-value < 0.05; *** adjusted p-value < 0.005. (B) Components enriched in the KEGG pathway ‘Rig-like receptor signaling’ in C8161 cells versus fibroblasts and A375 versus fibroblasts. +P, phosphorylation; +U, ubiquitination; −U, deubiquitination. Adjusted p-value < 0.05.
Cells 13 00226 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sadeq, S.; Chitcharoen, S.; Al-Hashimi, S.; Rattanaburi, S.; Casement, J.; Werner, A. Significant Variations in Double-Stranded RNA Levels in Cultured Skin Cells. Cells 2024, 13, 226. https://doi.org/10.3390/cells13030226

AMA Style

Sadeq S, Chitcharoen S, Al-Hashimi S, Rattanaburi S, Casement J, Werner A. Significant Variations in Double-Stranded RNA Levels in Cultured Skin Cells. Cells. 2024; 13(3):226. https://doi.org/10.3390/cells13030226

Chicago/Turabian Style

Sadeq, Shaymaa, Suwalak Chitcharoen, Surar Al-Hashimi, Somruthai Rattanaburi, John Casement, and Andreas Werner. 2024. "Significant Variations in Double-Stranded RNA Levels in Cultured Skin Cells" Cells 13, no. 3: 226. https://doi.org/10.3390/cells13030226

APA Style

Sadeq, S., Chitcharoen, S., Al-Hashimi, S., Rattanaburi, S., Casement, J., & Werner, A. (2024). Significant Variations in Double-Stranded RNA Levels in Cultured Skin Cells. Cells, 13(3), 226. https://doi.org/10.3390/cells13030226

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop