Next Article in Journal
Cloning and Functional Analysis of BcMYB101 Gene Involved in Leaf Development in Pak Choi (Brassica rapa ssp. Chinensis)
Next Article in Special Issue
The Feasibility of Host Transcriptome Profiling as a Diagnostic Tool for Microbial Etiology in Childhood Cancer Patients with Febrile Neutropenia
Previous Article in Journal
Neuroprotective Properties of Resveratrol and Its Derivatives—Influence on Potential Mechanisms Leading to the Development of Alzheimer’s Disease
Previous Article in Special Issue
A Meta-Analysis of Multiple Whole Blood Gene Expression Data Unveils a Diagnostic Host-Response Transcript Signature for Respiratory Syncytial Virus
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 8
10.3390/ijms21082748
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

RNA-Seq Data-Mining Allows the Discovery of Two Long Non-Coding RNA Biomarkers of Viral Infection in Humans

by
Ruth Barral-Arca
1,2,3,
Alberto Gómez-Carballa
1,2,3,
Miriam Cebey-López
1,2,3,
María José Currás-Tuala
1,2,3,
Sara Pischedda
1,2,3,
Sandra Viz-Lasheras
1,2,3,
Xabier Bello
1,2,3,
Federico Martinón-Torres
3,4 and
Antonio Salas
1,2,3,*
1
Unidade de Xenética, Instituto de Ciencias Forenses (INCIFOR), Facultade de Medicina, Universidade de Santiago de Compostela, 15782 Galicia, Spain
2
GenPoB Research Group, Instituto de Investigaciones Sanitarias (IDIS), Hospital Clínico Universitario de Santiago (SERGAS), 15706 Galicia, Spain
3
Genetics, Vaccines and Pediatric Infectious Diseases Research Group (GENVIP), Instituto de Investigación Sanitaria de Santiago (IDIS) and Universidad de Santiago de Compostela (USC), 15706 Galicia, Spain
4
Translational Pediatrics and Infectious Diseases, Department of Pediatrics, Hospital Clínico Universitario de Santiago de Compostela (SERGAS), 15706 Galicia, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(8), 2748; https://doi.org/10.3390/ijms21082748
Submission received: 10 March 2020 / Revised: 3 April 2020 / Accepted: 11 April 2020 / Published: 15 April 2020
(This article belongs to the Special Issue Host Infectomics in the Childhood)
Browse Figures
Versions Notes

Abstract

:
There is a growing interest in unraveling gene expression mechanisms leading to viral host invasion and infection progression. Current findings reveal that long non-coding RNAs (lncRNAs) are implicated in the regulation of the immune system by influencing gene expression through a wide range of mechanisms. By mining whole-transcriptome shotgun sequencing (RNA-seq) data using machine learning approaches, we detected two lncRNAs (ENSG00000254680 and ENSG00000273149) that are downregulated in a wide range of viral infections and different cell types, including blood monocluclear cells, umbilical vein endothelial cells, and dermal fibroblasts. The efficiency of these two lncRNAs was positively validated in different viral phenotypic scenarios. These two lncRNAs showed a strong downregulation in virus-infected patients when compared to healthy control transcriptomes, indicating that these biomarkers are promising targets for infection diagnosis. To the best of our knowledge, this is the very first study using host lncRNAs biomarkers for the diagnosis of human viral infections.

1. Introduction

The majority of erroneous antibiotic prescriptions occur in virus-infected patients, for which antibiotics offer no benefit except for mixed viral/bacterial infections [1,2]. Inadequate use of antibiotics not only is an economic burden for the sanitary systems—especially for developing countries—but also it increases the risk of adverse events and, more importantly, the emergence of resistant bacteria [2]. According to the World Health Organization (WHO), antibiotics overuse is one of the biggest threats to human health nowadays [3]. The development of a fast and accurate diagnostic testing to early distinguish viral from bacterial infections in clinical settings and hospitals would facilitate a reduction in the overuse of wide-spectrum antibiotics, helping physicians make the right decisions and fight the appearance of antibiotic-resistant bacteria.
The gold standard procedure to detect the presence of bacterial infections is bacterial culture, but this technique has important limitations. Firstly, obtaining results from cultures usually takes 48–72 h, a timeframe that might be inadequate for decision-making in terms of antibiotic prescriptions to children with suspected infection. Secondly, bacterial cultures have limited sensitivity, to the extent that, e.g., failure to detect causal microorganisms occurs in 50% of pneumonia patients in critical care units [4,5]. According to Rozon et al. [6], the sensitivity and specificity of a Gram stain from a good-quality sample in the diagnosis of pneumococcal pneumonia and H. influenzae ranged 57–82%. Besides, samples obtained from accessible sites (e.g., blood) can be useless for pathogen detection in cases where the infection is located in inaccessible anatomical sites. Third, viruses are usually detected in healthy children, and also in children with bacterial infections. For instance, in a large international collaborative study, Martinón-Torres et al. [7] recently reported that pathogen detection has been of limited help to distinguish viral from bacterial infection.
In 2007, Ramilo et al. [8] showed that RNA from blood leukocytes of children with life-threatening viral or bacterial infections bears host-pathogen specific transcriptional signatures. This finding represented a shift of paradigm in the study of infectious diseases since, for the first time, the host, and not the pathogen, became the focus of attention. Since then, RNA analysis has arisen as a powerful screening tool to find diagnostic biomarkers that may be used to develop new tests that overcome the limitations of bacterial culture. Recently, several studies have been exploring host-specific transcriptomic biomarkers that may allow distinguishing between viral and bacterial infections or pathogen-specific signatures [9,10,11,12,13,14,15,16]. Related to transcriptional signatures, there are also several studies relating host genetic susceptibility factors to infectious diseases [17,18,19,20,21].
Thanks to the recent advances of transcriptomics and bioinformatics, new long non-coding RNAs (lncRNAs) are being discovered and characterized each year. LncRNAs are transcripts with lengths exceeding 200 nucleotides that do not translate into proteins. Despite their abundance, the role of these molecules in infectious diseases has not been studied yet.
The present study sought viral lncRNAs host biomarkers that allow discriminating viral from healthy controls and assessed whether these genes could be the basis for a new diagnostic tool that could be implemented in the clinical routine.

2. Results

2.1. Evaluation of lncRNAs in the Context of Viral Infection

Our machine learning model detected two lncRNAs of viral infection, namely, ENSG00000273149 and ENSG00000254680. According to Ensembl (www.ensembl.org), ENSG00000273149 is a recently described transcript located in chromosome 13 and anti-sense transcribed with respect to the translationally-controlled tumor protein gene (TPT1). As ENSG00000273149 is a recently described transcript not much is known about it, except for the fact that is has been described to be associated with vascular stenosis [22]. ENSG00000254680 is located in chromosome 11; since it is also a recently described transcript, there is no published information about it.
We evaluated if patients clustered according to their disease status (viral infection, bacterial infection, and healthy controls) when applying the Viral Score (VS; see methods) to different groups represented by several pathogens and tissues. For this comparison, we generated boxplots with overlapping one-dimensional scatter plots, where high VS values indicate healthy status, whereas a lower VS indicates viral infection. This analysis showed that in most of the scenarios (Figure 1), there are significant differences in the VS of affected children by viral infection versus healthy controls (1A,1C,1D). The only exception occurs in the Mexican group (Figure 1B) where the boxplot shows some overlap in the first and third quartiles (75th percentile) of the boxplots for virus-infected patients versus controls. This could be explained by the fact that the Mexican samples represent a mixed group of children with mild and severe infections [23].
The diagnostic accuracy of the 2-lncRNA test was evaluated by ROC analysis (Figure 2, Table 1). For all the tested scenarios, the ROC curves indicate that the accuracy of the test based on the two lncRNAs is very high (AUC > 89%) when comparing viral infection against healthy controls in all the datasets. The accuracy drops moderately when comparing bacterial versus viral infection (AUC = 73%), and it drops slightly more when comparing bacteria versus controls (AUC = 70%):
  • Influenza versus controls (PRJNA230906), AUC = 100% (CI95%: 100–100%)
  • Dengue versus controls (GSE98859), AUC = 100% (CI95%: 100–100%)
  • Enterovirus/coxsakie virus (GSE94551) versus controls, AUC = 100% (CI95%: 100–100%)
  • RV versus controls (PRJNA325575), AUC = 100% (CI95%: 100–100%)
  • RV/norovirus/adenovirus versus controls (GSE69529), AUC = 89% (CI95%: 83–95%)
  • E. coli/Salmonella/Shigella versus controls (GSE69529), AUC = 70% (CI95%: 62–79%)
  • E. coli/Salmonella/Shigella versus RV/norovirus/adenovirus, AUC = 73% (CI95%: 66–80%)
  • Varicella zoster versus controls (GSE121385), AUC = 100% (CI95%: 100–100%).

2.2. The efficiency of lncRNAs versus Other Minimal Transcriptomic Signatures

We compared the performance of the two lncRNAs signature with the 2-transcript signature (IFI44L and FAM89) described by Herberg et al. [13] using microarray data, which was also identified as a minimum transcriptomic signature for distinguishing between viral and bacterial infections.
We found that the power to differentiate between viral and bacterial infection in our groups was very similar for the two signatures: AUClncRNA = 0.7311 [CI95%: 0.6615–0.8006] versus AUCIFI44L + FAM89 = 0.745 [CI95%: 0.6743–0.8157]. However, when comparing against healthy controls the performance of the two lncRNA is slightly better (Table 1).
We also evaluated the possibility of combining the information provided by the two lncRNAs with the expression level of IFI44L alone [24], and IFI44L + FAM89A. We found that the prediction accuracy in terms of AUC did not improve significantly (Table 1):
  • Influenza versus controls (PRJNA230906), AUClncRNAs + IFI44L + FAM89A = 100% (CI95%: 100–100%),
  • Dengue versus controls (GSE98859), AUClncRNAs + IFI44L + FAM89A = 100% (CI95%: 100–100%),
  • Enterovirus/coxsakie virus (GSE94551) versus controls, AUClncRNAs + IFI44L + FAM89A = NE; meaning that at least one of the genes of the signature was not expressed (NE) and the score could not be calculated,
  • RV versus healthy controls (PRJNA325575), AUClncRNAs + IFI44L + FAM89A = 100% (CI95%: 100–100%),
  • RV/norovirus/adenovirus versus controls (GSE69529), AUClncRNAs + IFI44L + FAM89A = 85% (CI95%: 77–92%),
  • E. coli/Salmonella/Shigella versus controls (GSE69529), AUClncRNAs + IFI44L + FAM89A = 75% (CI95%: 68–82%),
  • E. coli/Salmonella/Shigella versus RV/norovirus/adenovirus, AUClncRNAs + IFI44L + FAM89A = 68% (CI95%: 58–78%),
  • Varicella zoster versus controls (GSE121385), AUClncRNAs + IFI44L + FAM89A = 100% (CI95%: 100–100%).
In agreement with recent findings [24], the addition of FAM89A to the signature did not significantly improve (and even worsened) the signature performance.

3. Discussion

According to the Encyclopedia of DNA Elements (ENCODE; [25]) pilot project, most of the human genome (74–93%) transcribes in multiple RNAs, and the vast majority of them does not display protein coding capacity [26,27]. New sequence technologies (e.g., RNA-seq) have recently led to the discovery of many previously unknown no protein-coding lncRNAs. Even though these RNAs are encoded by an important proportion of the genome, most of their functions are still unknown [28]. This is because, unlike mRNAs or proteins, presently the function of lncRNAs cannot be predicted by looking at the sequence or structure. Their function is most likely related to chromatin modification, transcriptional regulation, and post-transcriptional regulation [29]. It is currently widely accepted that lncRNAs are involved in biological and pathological processes such as X dosage compensation in mammals, regulation of the immune response, gene imprinting, regulation of the cell cycle, telomere length, etc. [28,30]. The strongest evidence to support the role of host lncRNAs in the regulation of innate and adaptative immune system has come from animal models of infection and/or disease [30]. Thus, Peng et al. [31] observed that coronavirus infection in lung tissue altered globally the lncRNAs levels in several mouse strains. A few years later, experiments by Gomez et al. [32] carried out in transgenic mice infected with Salmonella and Theiler’s virus found that the upregulation of the lncRNA NeST increased the clearance of Salmonella infection, but also reduced resistance to the mouse Theiler’s picornavirus. According to their results, NeST induces interferon γ (IFN-γ), a cytokine that is critical for innate and adaptive immunity against viral, some bacterial and protozoal infections expression in T cells. Almost simultaneously, Zhang et al. [33] reported that the knockdown of the lncRNA NEAT1 enhances HIV-1 virus production. A year later, Imamura et al. [34] found that the NEAT1 promoted IL8 expression in response to viral infection in human cells. Nevertheless, it has been reported that a major problem in the study of lncRNA in animal models is the lack of evolutionary conservation of the lncRNAs between species, which constitutes a major barrier in extrapolating results from animal models to humans [30].
There is a growing body of evidence suggesting that lncRNAs play a role in the host susceptibility and defense against viral infections [30]. The present study used a data-driven approach based on comparing the whole transcriptome of healthy children and children suffering a viral infection. Following this method, we found two lncRNAs (ENSG00000254680 and ENSG00000273149) downregulated upon viral infection whose expression level can be used to detect the presence of viral infection. These two lncRNAs perfectly distinguish healthy controls from viral infected patients in a broad sense, which suggests that they may play a role in the host defense against viral infections or in the host susceptibility to infection. These two lncRNA molecules are abundant enough to be detected in a variety of different biological scenarios as seen in boxplots of Figure 2.
Little information is available for these two lncRNAs. The antisense transcript ENSG00000273149 is also known as AL138963.3, and it is located in chromosome 13 (http://www.ensembl.org/Homo_sapiens/Gene/Summary?g=ENSG00000273149;r=13:45340039-45341183;t=ENST00000610057). To the best of our knowledge, the only reference about it in the literature describes its upregulation in the tissue of arteriovenous fistula [22].
The lncRNA ENSG00000254680 (also known as AC079329.1) was found to be downregulated in viral infection. It is located in chromosome 11, and it also has been recently described (http://www.ensembl.org/Homo_sapiens/Gene/Summary?g=ENSG00000254680; r = 11:12261426-12263173; t = ENST00000527288). The only reference existing in the literature for this lncRNA corresponds to a patent describing methods to reduce T cells exhaustion, indicating that AC079329.1 is downregulated by c-Jun gene expression [35]. This description is in line with our observation because the c-Jun protein has proven to be very similar to an avian sarcoma virus 17 protein, which directly recognizes specific DNA sequences to regulate its expression [36].
The possibility of using only a few transcripts to diagnose viral infections renders these two lncRNAs attractive biomarkers to design a rapid point of care test, or even a qPCR-based assay, applicable in hospital settings [24,37]. It is worth noting that there already exists a commercial test that uses lncRNAs for diagnosis and prognosis purposes: the ExoDx™ Prostate (IntelliScore) urine test detects three biomarkers associated with aggressive prostate cancer, one of them is a lncRNA transcript [38,39].
Even so, the use of expression biomarkers, such as lncRNAs in diagnosis settings, is still limited because the RT-qPCR, the gold standard method for gene expression measurement, requires access to a laboratory, a thermocycler, and relatively sophisticated data analysis [24,38]. Nevertheless, recent advances in transcriptomic and portable technologies (e.g., Nanopore [4]), suggest that in the next few years we will probably see an explosion in the use of mRNAs and lncRNAs expression signatures as a point of care diagnostic tools for many pathologies [40,41,42,43].
In all of the phenotypic scenarios examined, our results suggest that the two described lncRNAs provide a signal strong enough to identify viral infections in a broad sense. The major limitation of the present study is the low sample size used to build the classification model; unfortunately, there are only a few RNA-seq studies on infectious diseases available [23,26]. Despite the limited amount of information available, our sample sizes are in line with other RNA-seq studies [44]. It is important to note that our meta-analysis could not benefit from the vast datasets publicly available on commercial microarrays because the lncRNAs found in our study were not included in these datasets.
Further research is thus needed to evaluate the accuracy of these non-coding biomarkers in different clinical scenarios, including different severities, the evolution of the signal according to the time from the onset of the disease, more microorganisms, scenarios of co-infection [45], etc. However, taken together, the present results suggest that the 2-lncRNAs signature has a good prediction capacity, comparable to the coding RNAs described in the literature to date.
Finally, the fact that the lncRNAs signature shows very good performance in different human population groups, tissues, and families of virus (such as Herpesviridae, Flaviviridae, Orthomyxoviridae, and Reoviridae), suggests that these biomarkers are most likely related to a molecular mechanism related to host response or susceptibility against viral infection that is well-preserved from an evolutionary point of view. Unravelling their role may allow the discovery of unknown pathogenic pathways and drug targets, which might eventually lead to the discovery of wide-spectrum antiviral drugs.

4. Conclusions

The present study represents a stepping stone to the ultimate goal of understanding virus–host interaction mechanisms in viral parthenogenesis. As far as we are aware, this is the first study to have found host lncRNAs with potential as viral diagnostic biomarkers or therapeutic targets. We have shown that two lncRNAs are downregulated during viral infections in blood and different types of cells compared to healthy control samples. Moving the present results to translational medicine and the future use of these markers as diagnostic tools in clinical settings would require preliminary testing on a wider range of well-selected samples from cultured cells and clinical samples. Moreover, further functional studies are needed to unravel the mechanisms by which these lncRNAs act during the infection.

5. Material and Methods

5.1. Sample Groups

Six groups of patients or infected cells and uninfected controls were analyzed in the present study (Figure 3).
The Spanish group consists of 18 western European children [46]. Blood samples were collected between 2013–2014 at the Hospital Clínico Universitario of Santiago de Compostela (Galicia; northwest Spain); they include: (i) 6 rotavirus (RV)-infected children that required medical attention at two different time-points, namely, acute (during medical attendance) and convalescent phases; and (ii) 6 healthy controls (with all the vaccines of the Spanish immunization schedule up to date but no RV vaccine); Table S1. A subset of these controls and infected children was previously explored in a separate study [26].
The Mexican group consists of 255 blood samples of healthy controls (n = 35), and patients with acute diarrhea caused by different bacterial/viral pathogens: RV (n = 60), norovirus (n = 7), adenovirus (n = 13), Salmonella (n = 42), Shigella (n = 38) and different strains of E. coli (n = 60); Gene Expression Omnibus (GEO) accession number: GSE69529 (more detailed information about the samples is provided in [47]).
The Chinese group comprises four blood samples collected from patients with H7N9 infection (n = 2) and healthy controls (n = 2); GEO accession number: PRJNA230906 (see [48] for additional information on the samples). This study analyzed the transcriptomic response of individuals to avian influenza virus (H7N9) infection using whole blood from infected and healthy adults from China using RNA-seq.
The Varicella Zoster infected fibroblast group is composed of six samples of Varicella Zoster Virus (VZV)-infected human dermal fibroblasts cell line (HDF) infected with different strains or vaccines (Suduvax® and Varivix®); uninfected cells (n = 1), cells infected with wildtype strains (n = 2), and cells exposed to the vaccine (n = 3); GEO accession number: GSE121385 (see [49] for additional information on the samples). This study focused on analyzing the effects of the wild type VZV and different vaccines using cell cultures; more specifically, primary human dermal fibroblasts infected with wild-type VZV and attenuated varicella strains.
The mononuclear cells group comprises five blood mononuclear cell samples infected with dengue virus (n = 4) and one uninfected cell sample (n = 1); GEO accession number: GSE98859. This study analyzed changes in the transcriptome of peripheral blood mononuclear cells (PBMCs) when undergoing dengue infection.
The human umbilical vein endothelial cells (HUVECs) group contains samples from enterovirus patients (n = 2), coxsackievirus patients (n = 2), and healthy controls (n = 2); GEO accession number: GSE94551 (see [50] for additional information on the samples).

5.2. Data Processing and Statistical Analysis

Quality control of total RNA, libraries preparation and whole transcriptome sequencing (RNA-seq) of the Spanish group (Table S1) were carried out following the protocol described before in [26]. All the RNA-seq raw data files of the present study were preprocessed and normalized as described in [23].
To explore the predictive value of lncRNAs in the context of viral infections, we first detected differentially expressed lncRNAs (meaning transcripts with biotype ‘lncRNA’ according to ENSEMBL) by comparing RV-children (acute and convalescent) versus controls. We focused on the transcripts belonging to the transcript biotype lncRNA according to BiomaRt [51] (Table S2). For this purpose, we employed the R package Deseq2 using a negative binomial generalized linear model, and including age and sex as covariates of the model [52].
We filtered the lncRNAs to obtain a list of the best candidates according to the following thresholds: P-adjusted < 0.05, and log2 FoldChange > |2|.
We then used a machine-learning algorithm to select the lncRNAs that are more differentially expressed. We applied the variable selection algorithm elastic net implemented in the R package glmnet to the list of filtered lncRNA, using the Spanish group (Table S1) as the training group, and the Mexican dataset GSE69529 as the test group. As the elastic net algorithm usually yields a model that includes many genes, (thus rendering its translation to clinical settings difficult) we looked for the most informative genes (minimal signature) among the ones filtered. The parameters needed for the calculation of elastic net were estimated using 10-fold cross-validation. Following a machine learning approach, a single-hidden-layer neural network model was fitted with the R package nnet [36], obtaining a 2-transcripts (ENSG00000273149 and ENSG00000254680) viral signature. Viral score (VS) was calculated analogously as in Herberg et al. [13], namely:
VS = log2(exprs(ENSG00000273149)) + log2(exprs(ENSG00000254680))
We next evaluated the performance of these markers to detect viral infections in a broader sense. For this purpose, we used other GEO datasets (PRJNA230906, GSE121385, GSE98859, and GSE94551) including controls and viral patient’s RNA-seq data from blood and different cell types samples. While the first study (PRJNA230906,) analyzed the transcriptome of individuals suffering a viral infection, the other three studies focused on analyzing host response to viral infections using cell cultures or single cell types.
To evaluate the diagnostic accuracy of the lncRNAs signature in all the included studies, we used receiver operating characteristic (ROC) curves and the area under the curves (AUC) using the pROC [53] package in R. The threshold value, defined as the point on the ROC curve that maximized sensitivity and specificity, was calculated using the R package OptimalCutPoints [54]. The calculation of the confidence intervals for sensitivity and specificity was based on a stratified bootstrap resampling.

Supplementary Materials

The Supplementary Materials are available online at https://www.mdpi.com/1422-0067/21/8/2748/s1.

Author Contributions

Conceptualization, A.S., R.B.-A. and F.M.-T.; methodology, R.B.-A., M.C.-L., M.J.C.-T., S.P., and S.V.-L.; software, formal analysis and data curation, A.S., R.B.-A., A.G.-C., and X.B.; writing—original draft preparation, A.S. and R.B.-A.; writing—review, editing and supervision, all authors; funding acquisition, A.S. and F.M.-T. All authors have read and agreed to the published version of the manuscript.

Funding

This study received support from the Instituto de Salud Carlos III: project GePEM (Instituto de Salud Carlos III(ISCIII)/PI16/01478/Cofinanciado FEDER), DIAVIR (Instituto de Salud Carlos III(ISCIII)/DTS19/00049/Cofinanciado FEDER; Proyecto de Desarrollo Tecnológico en Salud) and Resvi-Omics (Instituto de Salud Carlos III(ISCIII)/PI19/01039/Cofinanciado FEDER) and project BI-BACVIR (PRIS-3; Agencia de Conocimiento en Salud (ACIS)—Servicio Gallego de Salud (SERGAS)—Xunta de Galicia; Spain) given to A.S.; and project ReSVinext (Instituto de Salud Carlos III(ISCIII)/PI16/01569/Cofinanciado FEDER), and Enterogen (Instituto de Salud Carlos III(ISCIII)/ PI19/01090/Cofinanciado FEDER) given to F.M.-T.

Acknowledgments

We would like to acknowledge CESGA (Supercomputing Centre of Galicia, Santiago de Compostela, Spain) for its supercomputing availability, web hosting, and support. We would also like to thank Peter Plinsley for kindly granting us access to transcriptome information on the norovirus and adenovirus samples.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Gill, J.M.; Fleischut, P.; Haas, S.; Pellini, B.; Crawford, A.; Nash, D.B. Use of antibiotics for adult upper respiratory infections in outpatient settings: A national ambulatory network study. Fam. Med. 2006, 38, 349–354. [Google Scholar]
  2. Fleming-Dutra, K.E.; Hersh, A.L.; Shapiro, D.J.; Bartoces, M.; Enns, E.A.; File, T.M., Jr.; Finkelstein, J.A.; Gerber, J.S.; Hyun, D.Y.; Linder, J.A.; et al. Prevalence of inappropriate antibiotic prescriptions among us ambulatory care visits, 2010–2011. JAMA 2016, 315, 1864–1873. [Google Scholar] [CrossRef] [Green Version]
  3. Zaman, S.B.; Hussain, M.A.; Nye, R.; Mehta, V.; Mamun, K.T.; Hossain, N. A review on antibiotic resistance: Alarm bells are ringing. Cureus 2017, 9, e1403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Charalampous, T.; Kay, G.L.; Richardson, H.; Aydin, A.; Baldan, R.; Jeanes, C.; Rae, D.; Grundy, S.; Turner, D.J.; Wain, J.; et al. Nanopore metagenomics enables rapid clinical diagnosis of bacterial lower respiratory infection. Nat. Biotechnol. 2019, 37, 783–792. [Google Scholar] [CrossRef] [PubMed]
  5. Cookson, W.; Cox, M.J.; Moffatt, M.F. New opportunities for managing acute and chronic lung infections. Nat. Rev. Microbiol 2018, 16, 111–120. [Google Scholar] [CrossRef] [PubMed]
  6. Roson, B.; Carratala, J.; Verdaguer, R.; Dorca, J.; Manresa, F.; Gudiol, F. Prospective study of the usefulness of sputum gram stain in the initial approach to community-acquired pneumonia requiring hospitalization. Clin. Infect. Dis. 2000, 31, 869–874. [Google Scholar] [CrossRef]
  7. Martinón-Torres, F.; Salas, A.; Rivero-Calle, I.; Cebey-López, M.; Pardo-Seco, J.; Herberg, J.A.; Boeddha, N.P.; Klobassa, D.S.; Secka, F.; Paulus, S.; et al. Life-threatening infections in children in europe (the euclids project): A prospective cohort study. Lancet Child Adolesc. Health 2018, 2, 404–414. [Google Scholar] [CrossRef]
  8. Ramilo, O.; Allman, W.; Chung, W.; Mejias, A.; Ardura, M.; Glaser, C.; Wittkowski, K.M.; Piqueras, B.; Banchereau, J.; Palucka, A.K.; et al. Gene expression patterns in blood leukocytes discriminate patients with acute infections. Blood 2007, 109, 2066–2077. [Google Scholar] [CrossRef] [Green Version]
  9. van Houten, C.B.; de Groot, J.A.H.; Klein, A.; Srugo, I.; Chistyakov, I.; de Waal, W.; Meijssen, C.B.; Avis, W.; Wolfs, T.F.W.; Shachor-Meyouhas, Y.; et al. A host-protein based assay to differentiate between bacterial and viral infections in preschool children (opportunity): A double-blind, multicentre, validation study. Lancet Infect. Dis. 2017, 17, 431–440. [Google Scholar] [CrossRef]
  10. Srugo, I.; Klein, A.; Stein, M.; Golan-Shany, O.; Kerem, N.; Chistyakov, I.; Genizi, J.; Glazer, O.; Yaniv, L.; German, A.; et al. Validation of a novel assay to distinguish bacterial and viral infections. Pediatrics 2017, 140. [Google Scholar] [CrossRef] [Green Version]
  11. Mahajan, P.; Kuppermann, N.; Mejias, A.; Suarez, N.; Chaussabel, D.; Casper, T.C.; Smith, B.; Alpern, E.R.; Anders, J.; Atabaki, S.M.; et al. Association of rna biosignatures with bacterial infections in febrile infants aged 60 days or younger. JAMA 2016, 316, 846–857. [Google Scholar] [CrossRef] [PubMed]
  12. Kaforou, M.; Herberg, J.A.; Wright, V.J.; Coin, L.J.M.; Levin, M. Diagnosis of bacterial infection using a 2-transcript host rna signature in febrile infants 60 days or younger. JAMA 2017, 317, 1577–1578. [Google Scholar] [CrossRef] [PubMed]
  13. Herberg, J.A.; Kaforou, M.; Wright, V.J.; Shailes, H.; Eleftherohorinou, H.; Hoggart, C.J.; Cebey-López, M.; Carter, M.J.; Janes, V.A.; Gormley, S.; et al. Diagnostic test accuracy of a 2-transcript host rna signature for discriminating bacterial vs viral infection in febrile children. JAMA 2016, 316, 835–845. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Bhattacharya, S.; Rosenberg, A.F.; Peterson, D.R.; Grzesik, K.; Baran, A.M.; Ashton, J.M.; Gill, S.R.; Corbett, A.M.; Holden-Wiltse, J.; Topham, D.J.; et al. Transcriptomic biomarkers to discriminate bacterial from nonbacterial infection in adults hospitalized with respiratory illness. Sci. Rep. 2017, 7, 6548. [Google Scholar] [CrossRef] [PubMed]
  15. Sampson, D.L.; Fox, B.A.; Yager, T.D.; Bhide, S.; Cermelli, S.; McHugh, L.C.; Seldon, T.A.; Brandon, R.A.; Sullivan, E.; Zimmerman, J.J.; et al. A four-biomarker blood signature discriminates systemic inflammation due to viral infection versus other etiologies. Sci. Rep. 2017, 7, 2914. [Google Scholar] [CrossRef] [PubMed]
  16. Barral-Arca, R.; Gómez-Carballa, A.; Cebey-López, M.; Bello, X.; Martinón-Torres, F.; Salas, A. A meta-analysis of multiple whole blood gene expression data unveils a diagnostic host-response transcript signature for respiratory syncytial virus. Int. J. Mol. Sci. 2020, 21, 1831. [Google Scholar] [CrossRef] [Green Version]
  17. Salas, A.; Pardo-Seco, J.; Barral-Arca, R.; Cebey-López, M.; Gómez-Carballa, A.; Rivero-Calle, I.; Pischedda, S.; Curras-Tuala, M.J.; Amigo, J.; Gómez-Rial, J.; et al. Whole exome sequencing identifies new host genomic susceptibility factors in empyema caused by streptococcus pneumoniae in children: A pilot study. Genes (Basel) 2018, 9, 240. [Google Scholar] [CrossRef] [Green Version]
  18. Borghini, L.; Png, E.; Binder, A.; Wright, V.J.; Pinnock, E.; de Groot, R.; Hazelzet, J.; Emonts, M.; Van der Flier, M.; Schlapbach, L.J.; et al. Identification of regulatory variants associated with genetic susceptibility to meningococcal disease. Sci. Rep. 2019, 9, 6966. [Google Scholar] [CrossRef]
  19. Davila, S.; Wright, V.J.; Khor, C.C.; Sim, K.S.; Binder, A.; Breunis, W.B.; Inwald, D.; Nadel, S.; Betts, H.; Carrol, E.D.; et al. Genome-wide association study identifies variants in the cfh region associated with host susceptibility to meningococcal disease. Nat. Genet. 2010, 42, 772–776. [Google Scholar]
  20. Martinon-Torres, F.; Png, E.; Khor, C.C.; Davila, S.; Wright, V.J.; Sim, K.S.; Vega, A.; Fachal, L.; Inwald, D.; Nadel, S.; et al. Natural resistance to meningococcal disease related to cfh loci: Meta-analysis of genome-wide association studies. Sci. Rep. 2016, 6, 35842. [Google Scholar] [CrossRef] [Green Version]
  21. Salas, A.; Pardo-Seco, J.; Cebey-López, M.; Gómez-Carballa, A.; Obando-Pacheco, P.; Rivero-Calle, I.; Curras-Tuala, M.J.; Amigo, J.; Gómez-Rial, J.; Martinón-Torres, F.; et al. Whole exome sequencing reveals new candidate genes in host genomic susceptibility to respiratory syncytial virus disease. Sci. Rep. 2017, 7, 15888. [Google Scholar] [CrossRef] [PubMed]
  22. Yan, Y.; Ye, W.; Chen, Q.; Yang, L.; Zhang, L.; Liu, Y.; Zhou, X.; Wang, G. Differential expression profile of long non-coding rna in the stenosis tissue of arteriovenous fistula. Gene 2018, 664, 127–138. [Google Scholar] [CrossRef] [PubMed]
  23. Barral-Arca, R.; Pardo-Seco, J.; Martinón-Torres, F.; Salas, A. A 2-transcript host cell signature distinguishes viral from bacterial diarrhea and it is influenced by the severity of symptoms. Sci. Rep. 2018, 8, 8043. [Google Scholar] [CrossRef] [PubMed]
  24. Gómez-Carballa, A.; Cebey-López, M.; Pardo-Seco, J.; Barral-Arca, R.; Rivero-Calle, I.; Pischedda, S.; Curras-Tuala, M.J.; Gómez-Rial, J.; Barros, F.; Martinón-Torres, F.; et al. A qpcr expression assay of ifi44l gene differentiates viral from bacterial infections in febrile children. Sci. Rep. 2019, 9, 11780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Consortium, E.P.; Birney, E.; Stamatoyannopoulos, J.A.; Dutta, A.; Guigo, R.; Gingeras, T.R.; Margulies, E.H.; Weng, Z.; Snyder, M.; Dermitzakis, E.T.; et al. Identification and analysis of functional elements in 1% of the human genome by the encode pilot project. Nature 2007, 447, 799–816. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Salas, A.; Marco-Puche, G.; Trivino, J.C.; Gomez-Carballa, A.; Cebey-Lopez, M.; Rivero-Calle, I.; Vilanova-Trillo, L.; Rodriguez-Tenreiro, C.; Gomez-Rial, J.; Martinon-Torres, F. Strong down-regulation of glycophorin genes: A host defense mechanism against rotavirus infection. Infect. Genet. Evol. 2016, 44, 403–411. [Google Scholar] [CrossRef]
  27. Carninci, P.; Yasuda, J.; Hayashizaki, Y. Multifaceted mammalian transcriptome. Curr. Opin. Cell Biol. 2008, 20, 274–280. [Google Scholar] [CrossRef]
  28. Geisler, S.; Coller, J. Rna in unexpected places: Long non-coding rna functions in diverse cellular contexts. Nat. Rev. Mol. Cell Biol. 2013, 14, 699–712. [Google Scholar] [CrossRef] [Green Version]
  29. Mercer, T.R.; Dinger, M.E.; Mattick, J.S. Long non-coding rnas: Insights into functions. Nat. Rev. Genet. 2009, 10, 155–159. [Google Scholar] [CrossRef]
  30. Heward, J.A.; Lindsay, M.A. Long non-coding rnas in the regulation of the immune response. Trends Immunol. 2014, 35, 408–419. [Google Scholar] [CrossRef] [Green Version]
  31. Peng, X.; Gralinski, L.; Armour, C.D.; Ferris, M.T.; Thomas, M.J.; Proll, S.; Bradel-Tretheway, B.G.; Korth, M.J.; Castle, J.C.; Biery, M.C.; et al. Unique signatures of long noncoding rna expression in response to virus infection and altered innate immune signaling. mBio 2010, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Gomez, J.A.; Wapinski, O.L.; Yang, Y.W.; Bureau, J.F.; Gopinath, S.; Monack, D.M.; Chang, H.Y.; Brahic, M.; Kirkegaard, K. The nest long ncrna controls microbial susceptibility and epigenetic activation of the interferon-gamma locus. Cell 2013, 152, 743–754. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Zhang, Q.; Chen, C.Y.; Yedavalli, V.S.; Jeang, K.T. Neat1 long noncoding rna and paraspeckle bodies modulate hiv-1 posttranscriptional expression. mBio 2013, 4, e00596–e00612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Imamura, K.; Imamachi, N.; Akizuki, G.; Kumakura, M.; Kawaguchi, A.; Nagata, K.; Kato, A.; Kawaguchi, Y.; Sato, H.; Yoneda, M.; et al. Long noncoding rna neat1-dependent sfpq relocation from promoter region to paraspeckle mediates il8 expression upon immune stimuli. Mol. Cell 2014, 53, 393–406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Mackall, C.; Lynn, R.; Weber, E.; Sotillo, E. Compositions and methods for inhibiting t cell exhaustion. Patent application no. 16/221,155, 14 December 2018. [Google Scholar]
  36. Maki, Y.; Bos, T.J.; Davis, C.; Starbuck, M.; Vogt, P.K. Avian sarcoma virus 17 carries the jun oncogene. Proc. Natl. Acad. Sci. USA 1987, 84, 2848–2852. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Gliddon, H.D.; Herberg, J.A.; Levin, M.; Kaforou, M. Genome-wide host rna signatures of infectious diseases: Discovery and clinical translation. Immunology 2018, 153, 171–178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Byron, S.A.; Van Keuren-Jensen, K.R.; Engelthaler, D.M.; Carpten, J.D.; Craig, D.W. Translating rna sequencing into clinical diagnostics: Opportunities and challenges. Nat. Rev. Genet. 2016, 17, 257–271. [Google Scholar] [CrossRef]
  39. Du, Z.; Fei, T.; Verhaak, R.G.; Su, Z.; Zhang, Y.; Brown, M.; Chen, Y.; Liu, X.S. Integrative genomic analyses reveal clinically relevant long noncoding rnas in human cancer. Nat. Struct. Mol. Biol. 2013, 20, 908–913. [Google Scholar] [CrossRef] [Green Version]
  40. Tang, H.; Wu, Z.; Zhang, J.; Su, B. Salivary lncrna as a potential marker for oral squamous cell carcinoma diagnosis. Mol. Med. Rep. 2013, 7, 761–766. [Google Scholar] [CrossRef] [Green Version]
  41. Duan, W.; Du, L.; Jiang, X.; Wang, R.; Yan, S.; Xie, Y.; Yan, K.; Wang, Q.; Wang, L.; Zhang, X.; et al. Identification of a serum circulating lncrna panel for the diagnosis and recurrence prediction of bladder cancer. Oncotarget 2016, 7, 78850–78858. [Google Scholar] [CrossRef] [Green Version]
  42. Zhou, M.; Diao, Z.; Yue, X.; Chen, Y.; Zhao, H.; Cheng, L.; Sun, J. Construction and analysis of dysregulated lncrna-associated cerna network identified novel lncrna biomarkers for early diagnosis of human pancreatic cancer. Oncotarget 2016, 7, 56383–56394. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Yang, G.; Lu, X.; Yuan, L. Lncrna: A link between rna and cancer. Biochim. Biophys. Acta 2014, 1839, 1097–1109. [Google Scholar] [CrossRef] [PubMed]
  44. Conesa, A.; Madrigal, P.; Tarazona, S.; Gómez-Cabrero, D.; Cervera, A.; McPherson, A.; Szczesniak, M.W.; Gaffney, D.J.; Elo, L.L.; Zhang, X.; et al. A survey of best practices for rna-seq data analysis. Genome Biol. 2016, 17, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Cebey-López, M.; Herberg, J.; Pardo-Seco, J.; Gómez-Carballa, A.; Martinon-Torres, N.; Salas, A.; Martinón-Sánchez, J.M.; Gormley, S.; Sumner, E.; Fink, C.; et al. Viral co-infections in pediatric patients hospitalized with lower tract acute respiratory infections. PLoS ONE 2015, 10, e0136526. [Google Scholar] [CrossRef] [Green Version]
  46. Barral-Arca, R.; Pardo-Seco, J.; Bello, X.; Martinón-Torres, F.; Salas, A. Ancestry patterns inferred from massive rnaseq data. RNA 2019, 25, 857–868. [Google Scholar] [CrossRef]
  47. DeBerg, H.A.; Zaidi, M.B.; Altman, M.C.; Khaenam, P.; Gersuk, V.H.; Campos, F.D.; Perez-Martinez, I.; Meza-Segura, M.; Chaussabel, D.; Banchereau, J.; et al. Shared and organism-specific host responses to childhood diarrheal diseases revealed by whole blood transcript profiling. PLoS ONE 2018, 13, e0192082. [Google Scholar] [CrossRef] [Green Version]
  48. Mei, B.; Ding, X.; Xu, H.Z.; Wang, M.T. Global gene expression changes in human peripheral blood after h7n9 infection. Gene 2014, 551, 255–260. [Google Scholar] [CrossRef]
  49. Oh, S.J.; Lim, S.; Song, M.J.; Ahn, J.H.; Lee, C.H.; Shin, O.S. Whole transcriptome analyses reveal differential mrna and microrna expression profiles in primary human dermal fibroblasts infected with clinical or vaccine strains of varicella zoster virus. Pathogens 2019, 8, 183. [Google Scholar] [CrossRef] [Green Version]
  50. Song, J.; Hu, Y.; Li, J.; Zheng, H.; Wang, J.; Guo, L.; Ning, R.; Li, H.; Yang, Z.; Fan, H.; et al. Different microrna profiles reveal the diverse outcomes induced by ev71 and ca16 infection in human umbilical vein endothelial cells using high-throughput sequencing. PLoS ONE 2017, 12, e0177657. [Google Scholar] [CrossRef]
  51. Durinck, S.; Moreau, Y.; Kasprzyk, A.; Davis, S.; De Moor, B.; Brazma, A.; Huber, W. Biomart and bioconductor: A powerful link between biological databases and microarray data analysis. Bioinform 2005, 21, 3439–3440. [Google Scholar] [CrossRef] [Green Version]
  52. 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] [Green Version]
  53. Robin, X.; Turck, N.; Hainard, A.; Tiberti, N.; Lisacek, F.; Sanchez, J.C.; Muller, M. Proc: An open-source package for r and s+ to analyze and compare roc curves. BMC Bioinform. 2011, 12, 77. [Google Scholar] [CrossRef] [PubMed]
  54. López-Ratón, M.; Rodríguez-Álvarez, M.X.; Cadarso-Suárez, C.; Gude-Sampedro, F. Optimalcutpoints: An r package for selecting optimal cutpoints in diagnostic tests. J. Stat. Softw. 2014, 61, 1–36. [Google Scholar] [CrossRef] [Green Version]
Ijms 21 02748 g001 550
Figure 1. Classification performance based on the two lncRNA VS, considering different viral pathogens and studies. Box and whisker plots of VS values: (A) Discovery group (Spanish group); (B) Test group GSE69529; (C) External validation group GSE98859; (D) External validation group GSE121385; (E) External validation group GSE94551; and (F) External validation group PRJNA230906. For all box plots, the horizontal lines in boxes indicate median of the groups; the lower and upper sides of boxes interquartile ranges and the whiskers < 1 times the interquartile range. On the x-axis is the sample status and on the y-axis the VS calculated as log2 of the sum of counts of our 2-transcript diagnosis model. The acronyms stand for: RV (Rotavirus), VZV (Varicella Zoster Virus), diffuse-adhering Escherichia coli (DAEC), enteropathogenic Escherichia coli (EPEC), and enteroaggregative Escherichia coli (EAEC).
Figure 1. Classification performance based on the two lncRNA VS, considering different viral pathogens and studies. Box and whisker plots of VS values: (A) Discovery group (Spanish group); (B) Test group GSE69529; (C) External validation group GSE98859; (D) External validation group GSE121385; (E) External validation group GSE94551; and (F) External validation group PRJNA230906. For all box plots, the horizontal lines in boxes indicate median of the groups; the lower and upper sides of boxes interquartile ranges and the whiskers < 1 times the interquartile range. On the x-axis is the sample status and on the y-axis the VS calculated as log2 of the sum of counts of our 2-transcript diagnosis model. The acronyms stand for: RV (Rotavirus), VZV (Varicella Zoster Virus), diffuse-adhering Escherichia coli (DAEC), enteropathogenic Escherichia coli (EPEC), and enteroaggregative Escherichia coli (EAEC).
Ijms 21 02748 g001
Ijms 21 02748 g002 550
Figure 2. Evaluation of VS test performance. AUC (Area Under the Curve) values are provided in Table 1. ROC curves for VS: (A) Discovery group (Spanish group); (B) Test group GSE69529; (C) External validation group GSE98859; (D) External validation group GSE121385; (E) External validation group GSE94551; and (F) External validation group PRJNA230906. For all figures, the horizontal lines in boxes indicate median of the groups; the lower and upper sides of boxes interquartile ranges and the whiskers < 1 times the interquartile range. On the x-axis is the sample status and on the y-axis the VS calculated as log2 of the sum of counts of our 2-transcript diagnosis model.
Figure 2. Evaluation of VS test performance. AUC (Area Under the Curve) values are provided in Table 1. ROC curves for VS: (A) Discovery group (Spanish group); (B) Test group GSE69529; (C) External validation group GSE98859; (D) External validation group GSE121385; (E) External validation group GSE94551; and (F) External validation group PRJNA230906. For all figures, the horizontal lines in boxes indicate median of the groups; the lower and upper sides of boxes interquartile ranges and the whiskers < 1 times the interquartile range. On the x-axis is the sample status and on the y-axis the VS calculated as log2 of the sum of counts of our 2-transcript diagnosis model.
Ijms 21 02748 g002
Ijms 21 02748 g003 550
Figure 3. Diagram of the study design. The acronyms stand for: Rotateq® (RV5); rotavirus (RV); diffuse-adhering Escherichia coli (DAEC), enteropathogenic Escherichia coli (EPEC) and enteroaggregative Escherichia coli (EAEC), and Varicella Zoster Virus (VZV).
Figure 3. Diagram of the study design. The acronyms stand for: Rotateq® (RV5); rotavirus (RV); diffuse-adhering Escherichia coli (DAEC), enteropathogenic Escherichia coli (EPEC) and enteroaggregative Escherichia coli (EAEC), and Varicella Zoster Virus (VZV).
Ijms 21 02748 g003
Table
Table 1. Area under the curve (AUC) values and viral score (VS) for different pathogens and studies. Figures in round brackets show the CI95% (calculated using 2000 bootstrap replicates). No expression (NE) means that at least one of the transcripts of the signature was not expressed. CC: Chinese group; SC: Spanish group; VZF: Varicella Zoster infected fibroblasts; MNC: mononuclear cells with dengue, EC: human umbilical vein endothelial cells enterovirus group, MC: Mexican group.
Table 1. Area under the curve (AUC) values and viral score (VS) for different pathogens and studies. Figures in round brackets show the CI95% (calculated using 2000 bootstrap replicates). No expression (NE) means that at least one of the transcripts of the signature was not expressed. CC: Chinese group; SC: Spanish group; VZF: Varicella Zoster infected fibroblasts; MNC: mononuclear cells with dengue, EC: human umbilical vein endothelial cells enterovirus group, MC: Mexican group.
Database ComparisonnAUC Sensitivity Specificity VS
2-LncRNA
CCPRJNA230906virus vs. control41 (1–1)112.569
MNCGSE98859virus vs. control51 (1–1)114.939
ECGSE94551virus vs. control61 (1–1)114.043
SCPRJNA325575virus vs. control181 (1–1)117.535
MCGSE69529virus vs. control1150.891 (0.833–0.950)0.8570.8388.350
MCGSE69529virus vs. bacteria2200.731 (0.662–0.800)0.8140.6008.113
MCGSE69529bacteria vs. control1750.707 (0.621–0.793)0.650 0.7438.392
VZFGSE121385virus vs. control61 (1–1)1115.015
2-LncRNA + IFI44L
CCPRJNA230906virus vs. control41 (1–1)111.980
MNCGSE98859virus vs. control51 (1–1)112.596
ECGSE94551virus vs. controlNENENENENE
SCPRJNA325575virus vs. control181 (1–1)113.227
MCGSE69529virus vs. control1150.857 (0.784–0.929)0.800 0.8506.133
MCGSE69529virus vs. bacteria2200.724 (0.651–0.797)0.8550.5254.101
MCGSE69529bacteria vs. control1750.723 (0.623–0.824)0.7900.7146.883
VZFGSE121385virus vs. control61 (1–1)119.670
2-LncRNA + IFI44L + FAM89A
CCPRJNA230906virus vs. control41 (1–1)112.328
MNCGSE98859virus vs. control51 (1–1)114.899
ECGSE94551virus vs. controlNENENENENE
SCPRJNA325575virus vs. control181 (1–1)1111.640
MCGSE69529virus vs. control1150.850 (0.774–0.927)0.8000.8638.247
MCGSE69529virus vs. bacteria2200.753 (0.684–0.822)0.6300.7887.469
MCGSE69529bacteria vs. control1750.681 (0.579–0.783)0.7100.7148.773
VZFGSE121385virus vs. control61 (1–1)119.679
FAM89A + IFI44L
CCPRJNA230906virus vs. control41 (1–1)11−2.650
MNCGSE98859virus vs. control50.625 (0.500–0.875)10.250−0.091
ECGSE94551virus vs. controlNENENENENE
SCPRJNA325575virus vs. control181 (1–1)114.715
MCGSE69529virus vs. control1150.829 (0.749–0.910)0.7710.825−3.764
MCGSE69529virus vs. bacteria2200.745 (0.674–0.816)0.8130.638−4.834
MCGSE69529bacteria vs. control1750.647 (0.544–0.749)0.7480.629−2.924
VZFGSE121385virus vs. control61 (1–1)11−2.270

Share and Cite

MDPI and ACS Style

Barral-Arca, R.; Gómez-Carballa, A.; Cebey-López, M.; Currás-Tuala, M.J.; Pischedda, S.; Viz-Lasheras, S.; Bello, X.; Martinón-Torres, F.; Salas, A. RNA-Seq Data-Mining Allows the Discovery of Two Long Non-Coding RNA Biomarkers of Viral Infection in Humans. Int. J. Mol. Sci. 2020, 21, 2748. https://doi.org/10.3390/ijms21082748

AMA Style

Barral-Arca R, Gómez-Carballa A, Cebey-López M, Currás-Tuala MJ, Pischedda S, Viz-Lasheras S, Bello X, Martinón-Torres F, Salas A. RNA-Seq Data-Mining Allows the Discovery of Two Long Non-Coding RNA Biomarkers of Viral Infection in Humans. International Journal of Molecular Sciences. 2020; 21(8):2748. https://doi.org/10.3390/ijms21082748

Chicago/Turabian Style

Barral-Arca, Ruth, Alberto Gómez-Carballa, Miriam Cebey-López, María José Currás-Tuala, Sara Pischedda, Sandra Viz-Lasheras, Xabier Bello, Federico Martinón-Torres, and Antonio Salas. 2020. "RNA-Seq Data-Mining Allows the Discovery of Two Long Non-Coding RNA Biomarkers of Viral Infection in Humans" International Journal of Molecular Sciences 21, no. 8: 2748. https://doi.org/10.3390/ijms21082748

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