Network meta-analysis of transcriptome expression changes in different manifestations of dengue virus infection.

BACKGROUND: Several studies have been performed to study transcriptome profiles after dengue virus infections with partly different results. Due to slightly different settings of the individual studies, different genes and enriched gene sets are reported in these studies. The main aim of this network meta-analysis was to aggregate a selection of these studies to identify genes and gene sets that are more generally associated with dengue virus infection, i.e. with less dependence on the individual study settings. METHODS: We performed network meta-analysis by different approaches using publicly available gene expression data of five selected studies from the Gene Expression Omnibus database. The study network includes dengue fever (DF), hemorrhagic fever (DHF), shock syndrome (DSS) patients as well as convalescent and healthy control individuals. After data merging and missing value imputation, study-specific batch effects were removed. Pairwise differential expression analysis and subsequent gene-set enrichment analysis were performed between the five study groups. Furthermore, mutual information networks were derived from the top genes of each group comparison, and the separability between the three patient groups was studied by machine learning models. RESULTS: From the 10 possible pairwise group comparisons in the study network, six genes (IFI27, TPX2, CDT1, DTL, KCTD14 and CDCA3) occur with a noticeable frequency among the top listed genes of each comparison. Thus, there is an increased evidence that these genes play a general role in dengue virus infections. IFI27 and TPX2 have also been highlighted in the context of dengue virus infection by other studies. A few of the identified gene sets from the network meta-analysis overlap with findings from the original studies. Mutual information networks yield additional genes for which the observed pairwise correlation is different between the patient groups. Machine learning analysis shows a moderate separability of samples from the DF, DHF and DSS groups (accuracy about 80%). CONCLUSIONS: Due to an increased sample size, the network meta-analysis could reveal additional genes which are called differentially expressed between the studied groups and that may help to better understand the molecular basis of this disease.

A resampling strategy for studying robustness in virus detection pipelines.

Next-generation sequencing is regularly used to identify viral sequences in DNA or RNA samples of infected hosts. A major step of most pipelines for virus detection is to map sequence reads against known virus genomes. Due to small differences between the sequences of related viruses, and due to several biological or technical errors, mapping underlies uncertainties. As a consequence, the resulting list of detected viruses can lack robustness. A new approach for generating artificial sequencing reads together with a strategy of resampling from the original findings is proposed that can help to assess the robustness of the originally identified list of viruses. From the original mapping result in form of a SAM file, a set of statistical distributions are derived. These are used in the resampling pipeline to generate new artificial reads which are again mapped versus the reference genomes. By summarizing the resampling procedure, the analyst receives information about whether the presence of a particular virus in the sample gains or losses evidence, and thus about the robustness of the original mapping list but also that of individual viruses in this list. To judge robustness, several indicators are derived from the resampling procedure such as the correlation between original and resampling read counts, or the statistical detection of outliers in the differences of read counts. Additionally, graphical illustrations of read count shifts via Sankey diagrams are provided. To demonstrate the use of the new approach, the resampling approach is applied to three real-world data samples, one of them with laboratory-confirmed Influenza sequences, and to artificially generated data where virus sequences have been spiked into the sequencing data of a host. By applying the resampling pipeline, several viruses drop from the original list while new viruses emerge, showing robustness of those viruses that remain in the list. The evaluation of the new approach shows that the resampling approach is helpful to analyze the viral content of a biological sample, to rate the robustness of original findings and to better show the overall distribution of findings. The method is also applicable to other virus detection pipelines based on read mapping.