A new era of virus bioinformatics

Revealing the hierarchical structure of microbial communities

Measuring the dynamics of microbial communities results in high-dimensional measurements of taxa abundances over time and space, which is difficult to analyze due to complex changes in taxonomic compositions. This paper presents a new method to investigate and visualize the intrinsic hierarchical community structure implied by the measurements. The basic idea is to identify significant intersection sets, which can be seen as sub-communities making up the measured communities. Using the subset relationship, the intersection sets together with the measurements form a hierarchical structure visualized as a Hasse diagram. Chemical organization theory (COT) is used to relate the hierarchy of the sets of taxa to potential taxa interactions and to their potential dynamical persistence. The approach is demonstrated on a data set of community data obtained from bacterial 16S rRNA gene sequencing for samples collected monthly from four groundwater wells over a nearly 3-year period (n = 114) along a hillslope area. The significance of the hierarchies derived from the data is evaluated by showing that they significantly deviate from a random model. Furthermore, it is demonstrated how the hierarchy is related to temporal and spatial factors; and how the idea of a core microbiome can be extended to a set of interrelated core microbiomes. Together the results suggest that the approach can support developing models of taxa interactions in the future.

Identification of allograft inflammatory factor-1 suppressing the progression and indicating good prognosis of osteosarcoma

BACKGROUND

Osteosarcoma is one of the most common cancers worldwide. Intense efforts have been made to elucidate the pathogeny, but the mechanisms of osteosarcoma are still not well understood. We aimed to investigate the potential biomarker, allograft inflammatory factor-1 (AIF1), affecting the progression and prognosis of osteosarcoma.

METHODS

Three microarray datasets were downloaded from GEO datasets and one was obtained from the TCGA dataset. The differentially expressed genes (DEGs) were identified. GO and KEGG functional enrichment analyses of overlapped DEGs were performed. The PPI network of overlapped DEGs was constructed by STRING and visualized with Cytoscape. Overall survival (OS) and Metastasis free survival (MFS) were analyzed from GSE21257. Finally, the effect of the most relevant core gene affecting the progression of osteosarcoma was examined in vitro.

RESULTS

One hundred twenty six DEGs were identified, consisting of 65 upregulated and 61 downregulated genes. Only AIF1 was significantly associated with OS and MFS. It was found that AIF1 could be enriched into the NF-κB signaling pathway. GSEA and ssGSEA analyses showed that AIF1 was associated with the immune invasion of tumors. Cell experiments showed that AIF1 was underexpressed in osteosarcoma cell lines, while the malignant propriety was attenuated after overexpressing the expression of AIF1. Moreover, AIF1 also affects the expression of the NF-κB pathway.

CONCLUSION

In conclusion, DEGs and hub genes identified in the present study help us understand the molecular mechanisms underlying the carcinogenesis and progression of osteosarcoma, and provide candidate targets for diagnosis and treatment of osteosarcoma.

Bioinformatics analysis identifies a key gene HLA_DPA1 in severe influenza-associated immune infiltration

BACKGROUND

Severe influenza is a serious global health issue that leads to prolonged hospitalization and mortality on a significant scale. The pathogenesis of this infectious disease is poorly understood. Therefore, this study aimed to identify the key genes associated with severe influenza patients necessitating invasive mechanical ventilation.

METHODS

The current study utilized two publicly accessible gene expression profiles (GSE111368 and GSE21802) from the Gene Expression Omnibus database. The research focused on identifying the genes exhibiting differential expression between severe and non-severe influenza patients. We employed three machine learning algorithms, namely the Least Absolute Shrinkage and Selection Operator regression model, Random Forest, and Support Vector Machine-Recursive Feature Elimination, to detect potential key genes. The key gene was further selected based on the diagnostic performance of the target genes substantiated in the dataset GSE101702. A single-sample gene set enrichment analysis algorithm was applied to evaluate the participation of immune cell infiltration and their associations with key genes.

RESULTS

A total of 44 differentially expressed genes were recognized; among them, we focused on 10 common genes, namely PCOLCE2, HLA_DPA1, LOC653061, TDRD9, MPO, HLA_DQA1, MAOA, S100P, RAP1GAP, and CA1. To ensure the robustness of our findings, we employed overlapping LASSO regression, Random Forest, and SVM-RFE algorithms. By utilizing these algorithms, we were able to pinpoint the aforementioned 10 genes as potential biomarkers for distinguishing between both cases of influenza (severe and non-severe). However, the gene HLA_DPA1 has been recognized as a crucial factor in the pathological condition of severe influenza. Notably, the validation dataset revealed that this gene exhibited the highest area under the receiver operating characteristic curve, with a value of 0.891. The use of single-sample gene set enrichment analysis has provided valuable insights into the immune responses of patients afflicted with severe influenza that have further revealed a categorical correlation between the expression of HLA_DPA1 and lymphocytes.

CONCLUSION

The findings indicated that the HLA_DPA1 gene may play a crucial role in the immune-pathological condition of severe influenza and could serve as a promising therapeutic target for patients infected with severe influenza.

CDK1-SRC Interaction-Dependent Transcriptional Activation of HSP90AB1 Promotes Antitumor Immunity in Hepatocellular Carcinoma

This study aimed to analyze multiomics data and construct a regulatory network involving kinases, transcription factors, and immune genes in hepatocellular carcinoma (HCC) prognosis. The researchers used transcriptomic, proteomic, and clinical data from TCGA and GEO databases to identify immune genes associated with HCC. Statistical analysis, meta-analysis, and protein-protein interaction analyses were performed to identify key immune genes and their relationships. In vitro and in vivo experiments validated the CDK1-SRC-HSP90AB1 network's effects on HCC progression and antitumor immunity. A prognostic risk model was developed using clinicopathological features and immune infiltration. The immune genes LPA, BIRC5, HSP90AB1, ROBO1, and CCL20 were identified as the key prognostic factors. The CDK1-SRC-HSP90AB1 network promoted HCC cell proliferation and migration, with HSP90AB1 being transcriptionally activated by the CDK1-SRC interaction. Manipulating SRC or HSP90AB1 reversed the effects of CDK1 and SRC on HCC. The CDK1-SRC-HSP90AB1 network also influenced HCC tumor formation and antitumor immunity. Overall, this study highlights the importance of the CDK1-SRC-HSP90AB1 network as a crucial immune-regulatory network in the HCC prognosis.

Computing all persistent subspaces of a reaction-diffusion system

An algorithm is presented for computing a reaction-diffusion partial differential equation (PDE) system for all possible subspaces that can hold a persistent solution of the equation. For this, all possible sub-networks of the underlying reaction network that are distributed organizations (DOs) are identified. Recently it has been shown that a persistent subspace must be a DO. The algorithm computes the hierarchy of DOs starting from the largest by a linear programming approach using integer cuts. The underlying constraints use elementary reaction closures as minimal building blocks to guarantee local closedness and global self-maintenance, required for a DO. Additionally, the algorithm delivers for each subspace an affiliated set of organizational reactions and minimal compartmentalization that is necessary for this subspace to persist. It is proved that all sets of organizational reactions of a reaction network, as already DOs, form a lattice. This lattice contains all potentially persistent sets of reactions of all constrained solutions of reaction-diffusion PDEs. This provides a hierarchical structure of all persistent subspaces with regard to the species and also to the reactions of the reaction-diffusion PDE system. Here, the algorithm is described and the corresponding Python source code is provided. Furthermore, an analysis of its run time is performed and all models from the BioModels database as well as further examples are examined. Apart from the practical implications of the algorithm the results also give insights into the complexity of solving reaction-diffusion PDEs.

Revolutionizing Medical Microbiology: How Molecular and Genomic Approaches Are Changing Diagnostic Techniques

Molecular and genomic approaches have revolutionized medical microbiology by offering faster and more accurate diagnostic techniques for infectious diseases. Traditional methods, which include culturing microbes and biochemical testing, are time-consuming and may not detect antibiotic-resistant strains. In contrast, molecular and genomic methods, including polymerase chain reaction (PCR)-based techniques and whole-genome sequencing, provide rapid and precise detection of pathogens, early-stage diseases, and antibiotic-resistant strains. These approaches have advantages such as high sensitivity and specificity, the potential for targeted therapies, and personalized medicine. However, implementing molecular and genomic techniques faces challenges related to cost, equipment, expertise, and data analysis. Ethical and legal considerations regarding patient privacy and genetic data usage also arise. Nonetheless, the future of medical microbiology lies in the widespread adoption of molecular and genomic approaches, which can lead to improved patient outcomes and the identification of antibiotic-resistant strains. Continued advancements, education, and exploration of ethical implications are necessary to fully harness the potential of molecular and genomic techniques in medical microbiology.

The International Virus Bioinformatics Meeting 2023

The 2023 International Virus Bioinformatics Meeting was held in Valencia, Spain, from 24-26 May 2023, attracting approximately 180 participants worldwide. The primary objective of the conference was to establish a dynamic scientific environment conducive to discussion, collaboration, and the generation of novel research ideas. As the first in-person event following the SARS-CoV-2 pandemic, the meeting facilitated highly interactive exchanges among attendees. It served as a pivotal gathering for gaining insights into the current status of virus bioinformatics research and engaging with leading researchers and emerging scientists. The event comprised eight invited talks, 19 contributed talks, and 74 poster presentations across eleven sessions spanning three days. Topics covered included machine learning, bacteriophages, virus discovery, virus classification, virus visualization, viral infection, viromics, molecular epidemiology, phylodynamic analysis, RNA viruses, viral sequence analysis, viral surveillance, and metagenomics. This report provides rewritten abstracts of the presentations, a summary of the key research findings, and highlights shared during the meeting.

Identification of key genes in spontaneous cerebral hemorrhage and prevention of disease damage: LASSO and SVM regression

Prevention is more important than treatment, and the incidence of intracerebral hemorrhage can be effectively reduced by intervening on the risk factors of intracerebral hemorrhage. By studying the risk factors of spontaneous intracerebral hemorrhage, we can identify the risk factors to achieve the target of treatment and prevention. Through the use of the Least Absolute Shrinkage and Selection Operator (LASSO) and the Support Vector Machine (SVM), the two essential SICH-related genes, NUAK1 and ERO1L, were eliminated from consideration. A Venn analysis was performed, and based on the two important modules, it found that SICH was related with four critical genes: VCM1, CRNDE, COL6A2, and HSPB6. One gene (NUAK1) was dramatically downregulated in the illness group compared to the control group, whereas three essential genes (ERO1L, VCAM1, and COL6A2) were significantly upregulated in the disease group. In the end, the genes ERO1L, VCAM1, COL6A2, and NUAK1 were shown to be the most important ones for SICH. It is anticipated that these genes will become novel biomarkers as well as targets for the development of new pharmacotherapies for SICH.

Navigating the Landscape: A Comprehensive Review of Current Virus Databases

Viruses are abundant and diverse entities that have important roles in public health, ecology, and agriculture. The identification and surveillance of viruses rely on an understanding of their genome organization, sequences, and replication strategy. Despite technological advancements in sequencing methods, our current understanding of virus diversity remains incomplete, highlighting the need to explore undiscovered viruses. Virus databases play a crucial role in providing access to sequences, annotations and other metadata, and analysis tools for studying viruses. However, there has not been a comprehensive review of virus databases in the last five years. This study aimed to fill this gap by identifying 24 active virus databases and included an extensive evaluation of their content, functionality and compliance with the FAIR principles. In this study, we thoroughly assessed the search capabilities of five database catalogs, which serve as comprehensive repositories housing a diverse array of databases and offering essential metadata. Moreover, we conducted a comprehensive review of different types of errors, encompassing taxonomy, names, missing information, sequences, sequence orientation, and chimeric sequences, with the intention of empowering users to effectively tackle these challenges. We expect this review to aid users in selecting suitable virus databases and other resources, and to help databases in error management and improve their adherence to the FAIR principles. The databases listed here represent the current knowledge of viruses and will help aid users find databases of interest based on content, functionality, and scope. The use of virus databases is integral to gaining new insights into the biology, evolution, and transmission of viruses, and developing new strategies to manage virus outbreaks and preserve global health.

Differences between Omicron SARS-CoV-2 RBD and other variants in their ability to interact with cell receptors and monoclonal antibodies

SARS-CoV-2 remains a health threat with the continuous emergence of new variants. This work aims to expand the knowledge about the SARS-CoV-2 receptor-binding domain (RBD) interactions with cell receptors and monoclonal antibodies (mAbs). By using constant-pH Monte Carlo simulations, the free energy of interactions between the RBD from different variants and several partners (Angiotensin-Converting Enzyme-2 (ACE2) polymorphisms and various mAbs) were predicted. Computed RBD-ACE2-binding affinities were higher for two ACE2 polymorphisms (rs142984500 and rs4646116) typically found in Europeans which indicates a genetic susceptibility. This is amplified for Omicron (BA.1) and its sublineages BA.2 and BA.3. The antibody landscape was computationally investigated with the largest set of mAbs so far in the literature. From the 32 studied binders, groups of mAbs were identified from weak to strong binding affinities (e.g. S2K146). These mAbs with strong binding capacity and especially their combination are amenable to experimentation and clinical trials because of their high predicted binding affinities and possible neutralization potential for current known virus mutations and a universal coronavirus.Communicated by Ramaswamy H. Sarma.

Inhibition of apoptosis through AKT-mTOR pathway in ovarian cancer and renal cancer

OBJECTIVE

Ovarian cancer and renal cancer are malignant tumors; however, the relationship between TTK Protein Kinase (TTK), AKT-mTOR pathway and ovarian cancer, renal cancer remains unclear.

METHODS

Download GSE36668 and GSE69428 from Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was performed. Created protein-protein interaction (PPI) network. Used Gene Ontology analysis (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) for functional enrichment analysis. Gene Set Enrichment Analysis (GSEA) analysis and survival analysis were performed. Created animal model for western blot analysis. Gene Expression Profiling Interactive Analysis (GEPIA) was performed to explore the role of TTK on the overall survival of renal cancer.

RESULTS

GO showed that DEGs were enriched in anion and small molecule binding, and DNA methylation. KEGG analysis presented that they mostly enriched in cholesterol metabolism, type 1 diabetes, sphingolipid metabolism, ABC transporters, etc., TTK, mTOR, p-mTOR, AKT, p-AKT, 4EBP1, p-4EBP1 and Bcl-2 are highly expressed in ovarian cancer, Bax, Caspase3 are lowly expressed in ovarian cancer, cell apoptosis is inhibited, leading to deterioration of ovarian cancer. Furthermore, the TTK was not only the hub biomarker of ovarian cancer, but also one significant hub gene of renal cancer, and its expression was up-regulated in the renal cancer. Compared with the renal cancer patients with low expression of TTK, the patients with high expression of TTK have the poor overall survival (P = 0.0021).

CONCLUSION

TTK inhibits apoptosis through AKT-mTOR pathway, worsening ovarian cancer. And TTK was also one significant hub biomarker of renal cancer.

CCNB2 as a potential biomarker of bladder cancer via the high throughput technology

Bladder cancer and oral squamous cell carcinoma (OSCC) seriously affect people's health. However, the relationship between bladder cancer and OSCC remains unclear. Got GSE138206, GSE146483, GSE184616, and bladder cancer datasets GSE65635, GSE100926 from Gene Expression Omnibus database. Weighted gene co-expression network analysis was used to identify the significant module. Functional enrichment analysis was performed via the Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes. Furthermore, the Gene Set Enrichment Analysis was also used to complete the enrichment analysis. Comparative Toxicogenomics Database found most relevant diseases to core genes. TargetScan is used to forecast analysis of microRNA and target genes. In Gene Ontology analysis, differentially expressed genes were mostly concentrated in cell differentiation, extrallular region, structural molecule activity, and actin binding. In Kyoto Encyclopedia of Genes and Genomes analysis, the differentially expressed genes were mainly enriched in PI3K-Akt signaling pathway, pathway in cancer, and extracellular matrix-receptor interaction. Seven hub genes (cyclin B2 [CCNB2], TK1, CDC20, PCNA, CKS1B, CDCA5, MCM4) were obtained. Hub genes (CCNB2, CDC20) are highly expressed in OSCC and bladder cancer samples. CCNB2 was one common oncogene of bladder cancer and OSCC.

Towards an optimal monoclonal antibody with higher binding affinity to the receptor-binding domain of SARS-CoV-2 spike proteins from different variants

A highly efficient and robust multiple scales in silico protocol, consisting of atomistic Molecular Dynamics (MD), coarse-grain (CG) MD, and constant-pH CG Monte Carlo (MC), has been developed and used to study the binding affinities of selected antigen-binding fragments of the monoclonal antibody (mAbs) CR3022 and several of its here optimized versions against 11 SARS-CoV-2 variants including the wild type. Totally 235,000 mAbs structures were initially generated using the RosettaAntibodyDesign software, resulting in top 10 scored CR3022-like-RBD complexes with critical mutations and compared to the native one, all having the potential to block virus-host cell interaction. Of these 10 finalists, two candidates were further identified in the CG simulations to be the best against all SARS-CoV-2 variants. Surprisingly, all 10 candidates and the native CR3022 exhibited a higher affinity for the Omicron variant despite its highest number of mutations. The multiscale protocol gives us a powerful rational tool to design efficient mAbs. The electrostatic interactions play a crucial role and appear to be controlling the affinity and complex building. Studied mAbs carrying a more negative total net charge show a higher affinity. Structural determinants could be identified in atomistic simulations and their roles are discussed in detail to further hint at a strategy for designing the best RBD binder. Although the SARS-CoV-2 was specifically targeted in this work, our approach is generally suitable for many diseases and viral and bacterial pathogens, leukemia, cancer, multiple sclerosis, rheumatoid, arthritis, lupus, and more.

Evolving ribonucleocapsid assembly/packaging signals in the genomes of the human and animal coronaviruses: targeting, transmission and evolution

A world-wide COVID-19 pandemic intensified strongly the studies of molecular mechanisms related to the coronaviruses. The origin of coronaviruses and the risks of human-to-human, animal-to-human and human-to-animal transmission of coronaviral infections can be understood only on a broader evolutionary level by detailed comparative studies. In this paper, we studied ribonucleocapsid assembly-packaging signals (RNAPS) in the genomes of all seven known pathogenic human coronaviruses, SARS-CoV, SARS-CoV-2, MERS-CoV, HCoV-OC43, HCoV-HKU1, HCoV-229E and HCoV-NL63 and compared them with RNAPS in the genomes of the related animal coronaviruses including SARS-Bat-CoV, MERS-Camel-CoV, MHV, Bat-CoV MOP1, TGEV and one of camel alphacoronaviruses. RNAPS in the genomes of coronaviruses were evolved due to weakly specific interactions between genomic RNA and N proteins in helical nucleocapsids. Combining transitional genome mapping and Jaccard correlation coefficients allows us to perform the analysis directly in terms of underlying motifs distributed over the genome. In all coronaviruses, RNAPS were distributed quasi-periodically over the genome with the period about 54 nt biased to 57 nt and to 51 nt for the genomes longer and shorter than that of SARS-CoV, respectively. The comparison with the experimentally verified packaging signals for MERS-CoV, MHV and TGEV proved that the distribution of particular motifs is strongly correlated with the packaging signals. We also found that many motifs were highly conserved in both characters and positioning on the genomes throughout the lineages that make them promising therapeutic targets. The mechanisms of encapsidation can affect the recombination and co-infection as well.Communicated by Ramaswamy H. Sarma.

Coexpression Network Analysis-Based Identification of Critical Genes Differentiating between Latent and Active Tuberculosis

Methods

Three Gene Expression Omnibus (GEO) microarray datasets (GSE19491, GSE98461, and GSE152532) were downloaded, with GSE19491 and GSE98461 then being merged to form a training dataset. Hub genes capable of differentiating between ATB and LTBI were then identified through differential expression analyses and a WGCNA analysis of this training dataset. Receiver operating characteristic (ROC) curves were then used to gauge to the diagnostic accuracy of these hub genes in the test dataset (GSE152532). Gene expression-based immune cell infiltration and the relationship between such infiltration and hub gene expression were further assessed via a single-sample gene set enrichment analysis (ssGSEA).

Results

In total, 485 differentially expressed genes were analyzed, with the WGCNA approach yielding 8 coexpression models. Of these, the black module was the most closely correlated with ATB. In total, five hub genes (FBXO6, ATF3, GBP1, GBP4, and GBP5) were identified as potential biomarkers associated with LTBI progression to ATB based on a combination of differential expression and LASSO analyses. The area under the ROC curve values for these five genes ranged from 0.8 to 0.9 in the test dataset, and ssGSEA revealed the expression of these genes to be negatively correlated with lymphocyte activity but positively correlated with myeloid and inflammatory cell activity.

Conclusion

The five hub genes identified in this study may play a novel role in tuberculosis-related immunopathology and offer value as novel biomarkers differentiating LTBI from ATB.

Co-expression network analysis identifies potential candidate hub genes in severe influenza patients needing invasive mechanical ventilation

BACKGROUND

Influenza is a contagious disease that affects people of all ages and is linked to considerable mortality during epidemics and occasional outbreaks. Moreover, effective immunological biomarkers are needed for elucidating aetiology and preventing and treating severe influenza. Herein, we aimed to evaluate the key genes linked with the disease severity in influenza patients needing invasive mechanical ventilation (IMV). Three gene microarray data sets (GSE101702, GSE21802, and GSE111368) from blood samples of influenza patients were made available by the Gene Expression Omnibus (GEO) database. The GSE101702 and GSE21802 data sets were combined to create the training set. Hub indicators for IMV patients with severe influenza were determined using differential expression analysis and Weighted correlation network analysis (WGCNA) from the training set. The receiver operating characteristic curve (ROC) was also used to evaluate the hub genes from the test set's diagnostic accuracy. Different immune cells' infiltration levels in the expression profile and their correlation with hub gene markers were examined using single-sample gene set enrichment analysis (ssGSEA).

RESULTS

In the present study, we evaluated a total of 447 differential genes. WGCNA identified eight co-expression modules, with the red module having the strongest correlation with IMV patients. Differential genes were combined to obtain 3 hub genes (HLA-DPA1, HLA-DRB3, and CECR1). The identified genes were investigated as potential indicators for patients with severe influenza who required IMV using the least absolute shrinkage and selection operator (LASSO) approach. The ROC showed the diagnostic value of the three hub genes in determining the severity of influenza. Using ssGSEA, it has been revealed that the expression of key genes was negatively correlated with neutrophil activation and positively associated with adaptive cellular immune response.

CONCLUSION

We evaluated three novel hub genes that could be linked to the immunopathological mechanism of severe influenza patients who require IMV treatment and could be used as potential biomarkers for severe influenza prevention and treatment.

Identification of a 5-lncRNA-Based Signature for Immune Characteristics and Prognosis of Lung Squamous Cell Carcinoma and Verification of the Function of lncRNA SPATA41

Lung squamous cell carcinoma (LUSC) is one of the most lethal cancers worldwide. Traditional tumor-node-metastasis (TNM) staging system has many insufficiencies in predicting immune characteristics, overall survival (OS), and prognosis of LUSC. LncRNA is currently found involved in tumor development and effectively predicts tumor prognosis. We screened potential tumor-related lncRNAs for immune characteristics and constructed a nomogram combining lncRNA and traditional clinical indicators for prognosis prediction. We obtained the large-scale gene expression profiles of samples from 492 LUSC patients in The Cancer Genome Atlas database. SPATA41, AL034550.2, AP003721.2, AC106786.1, and AC078889.1 were finally screened to construct a 5-lncRNA-based signature. The risk score of the signature divided patients into subgroups of high-risk and low-risk with significant differences in OS. Their area under the curve (AUC) reached more than 0.70 in 1, 3, and 5 years. In addition, compared with the high-risk subgroup, the low-risk subgroup exhibited a remarkably favorable prognosis and TME score, along with a higher immune infiltration score and lower TIDE score. The signature also significantly related to chemotherapy response, especially in cisplatin, vinorelbine, and paclitaxel. Importantly, the nomogram we constructed had good reliability with the assessment of the calibration chart and consistency index (c-index). GO and KEGG enrichment analysis indicated that co-expression mRNAs of the 5 lncRNAs were mainly focused on RNA splicing, DNA replication, and protein serine/threonine kinase activity. Functional assays demonstrated that SPATA41, one of the five OS-related lncRNAs, regulated invasion, migration, proliferation, and programmed death in vitro. In summary, our 5-lncRNA-based signature has a good performance in predicting immune characteristics and prognosis of LUSC patients.

Women in the European Virus Bioinformatics Center

Viruses are the cause of a considerable burden to human, animal and plant health, while on the other hand playing an important role in regulating entire ecosystems. The power of new sequencing technologies combined with new tools for processing "Big Data" offers unprecedented opportunities to answer fundamental questions in virology. Virologists have an urgent need for virus-specific bioinformatics tools. These developments have led to the formation of the European Virus Bioinformatics Center, a network of experts in virology and bioinformatics who are joining forces to enable extensive exchange and collaboration between these research areas. The EVBC strives to provide talented researchers with a supportive environment free of gender bias, but the gender gap in science, especially in math-intensive fields such as computer science, persists. To bring more talented women into research and keep them there, we need to highlight role models to spark their interest, and we need to ensure that female scientists are not kept at lower levels but are given the opportunity to lead the field. Here we showcase the work of the EVBC and highlight the achievements of some outstanding women experts in virology and viral bioinformatics.

The International Virus Bioinformatics Meeting 2022

The International Virus Bioinformatics Meeting 2022 took place online, on 23-25 March 2022, and has attracted about 380 participants from all over the world. The goal of the meeting was to provide a meaningful and interactive scientific environment to promote discussion and collaboration and to inspire and suggest new research directions and questions. The participants created a highly interactive scientific environment even without physical face-to-face interactions. This meeting is a focal point to gain an insight into the state-of-the-art of the virus bioinformatics research landscape and to interact with researchers in the forefront as well as aspiring young scientists. The meeting featured eight invited and 18 contributed talks in eight sessions on three days, as well as 52 posters, which were presented during three virtual poster sessions. The main topics were: SARS-CoV-2, viral emergence and surveillance, virus-host interactions, viral sequence analysis, virus identification and annotation, phages, and viral diversity. This report summarizes the main research findings and highlights presented at the meeting.

Up State of the SARS-COV-2 Spike Homotrimer Favors an Increased Virulence for New Variants

The COVID-19 pandemic has spread worldwide. However, as soon as the first vaccines-the only scientifically verified and efficient therapeutic option thus far-were released, mutations combined into variants of SARS-CoV-2 that are more transmissible and virulent emerged, raising doubts about their efficiency. This study aims to explain possible molecular mechanisms responsible for the increased transmissibility and the increased rate of hospitalizations related to the new variants. A combination of theoretical methods was employed. Constant-pH Monte Carlo simulations were carried out to quantify the stability of several spike trimeric structures at different conformational states and the free energy of interactions between the receptor-binding domain (RBD) and angiotensin-converting enzyme II (ACE2) for the most worrying variants. Electrostatic epitopes were mapped using the PROCEEDpKa method. These analyses showed that the increased virulence is more likely to be due to the improved stability to the S trimer in the opened state, in which the virus can interact with the cellular receptor, ACE2, rather than due to alterations in the complexation RBD-ACE2, since the difference observed in the free energy values was small (although more attractive in general). Conversely, the South African/Beta variant (B.1.351), compared with the SARS-CoV-2 wild type (wt), is much more stable in the opened state with one or two RBDs in the up position than in the closed state with three RBDs in the down position favoring the infection. Such results contribute to understanding the natural history of disease and indicate possible strategies for developing new therapeutic molecules and adjusting the vaccine doses for higher B-cell antibody production.

The SARS-CoV-2 SSHHPS Recognized by the Papain-like Protease

Viral proteases are highly specific and recognize conserved cleavage site sequences of ∼6-8 amino acids. Short stretches of homologous host-pathogen sequences (SSHHPS) can be found spanning the viral protease cleavage sites. We hypothesized that these sequences corresponded to specific host protein targets since >40 host proteins have been shown to be cleaved by Group IV viral proteases and one Group VI viral protease. Using PHI-BLAST and the viral protease cleavage site sequences, we searched the human proteome for host targets and analyzed the hit results. Although the polyprotein and host proteins related to the suppression of the innate immune responses may be the primary targets of these viral proteases, we identified other cleavable host proteins. These proteins appear to be related to the virus-induced phenotype associated with Group IV viruses, suggesting that information about viral pathogenesis may be extractable directly from the viral genome sequence. Here we identify sequences cleaved by the SARS-CoV-2 papain-like protease (PLpro) in vitro within human MYH7 and MYH6 (two cardiac myosins linked to several cardiomyopathies), FOXP3 (an X-linked Treg cell transcription factor), ErbB4 (HER4), and vitamin-K-dependent plasma protein S (PROS1), an anticoagulation protein that prevents blood clots. Zinc inhibited the cleavage of these host sequences in vitro. Other patterns emerged from multispecies sequence alignments of the cleavage sites, which may have implications for the selection of animal models and zoonosis. SSHHPS/nsP is an example of a sequence-specific post-translational silencing mechanism.

Comparative Transcriptome Profiling of Human and Pig Intestinal Epithelial Cells after Porcine Deltacoronavirus Infection

Porcine deltacoronavirus (PDCoV) is an emerging infectious disease of swine with zoonotic potential. Phylogenetic analysis suggests that PDCoV originated recently from a host-switching event between birds and mammals. Little is known about how PDCoV interacts with its differing hosts. Human-derived cell lines are susceptible to PDCoV infection. Herein, we compare the gene expression profiles of an established host swine cells to potential emerging host human cells after infection with PDCoV. Cell lines derived from intestinal lineages were used to reproduce the primary sites of viral infection in the host. Porcine intestinal epithelial cells (IPEC-J2) and human intestinal epithelial cells (HIEC) were infected with PDCoV. RNA-sequencing was performed on total RNA extracted from infected cells. Human cells exhibited a more pronounced response to PDCoV infection in comparison to porcine cells with more differentially expressed genes (DEGs) in human, 7486, in comparison to pig cells, 1134. On the transcriptional level, the adoptive host human cells exhibited more DEGs in response to PDCoV infection in comparison to the primary pig host cells, where different types of cytokines can control PDCoV replication and virus production. Key immune-associated DEGs and signaling pathways are shared between human and pig cells during PDCoV infection. These included genes related to the NF-kappa-B transcription factor family, the interferon (IFN) family, the protein-kinase family, and signaling pathways such as the apoptosis signaling pathway, JAK-STAT signaling pathway, inflammation/cytokine–cytokine receptor signaling pathway. MAP4K4 was unique in up-regulated DEGs in humans in the apoptosis signaling pathway. While similarities exist between human and pig cells in many pathways, our research suggests that the adaptation of PDCoV to the porcine host required the ability to down-regulate many response pathways including the interferon pathway. Our findings provide an important foundation that contributes to an understanding of the mechanisms of PDCoV infection across different hosts. To our knowledge, this is the first report of transcriptome analysis of human cells infected by PDCoV.

Structure and Hierarchy of SARS-CoV-2 Infection Dynamics Models Revealed by Reaction Network Analysis

This work provides a mathematical technique for analyzing and comparing infection dynamics models with respect to their potential long-term behavior, resulting in a hierarchy integrating all models. We apply our technique to coupled ordinary and partial differential equation models of SARS-CoV-2 infection dynamics operating on different scales, that is, within a single organism and between several hosts. The structure of a model is assessed by the theory of chemical organizations, not requiring quantitative kinetic information. We present the Hasse diagrams of organizations for the twelve virus models analyzed within this study. For comparing models, each organization is characterized by the types of species it contains. For this, each species is mapped to one out of four types, representing uninfected, infected, immune system, and bacterial species, respectively. Subsequently, we can integrate these results with those of our former work on Influenza-A virus resulting in a single joint hierarchy of 24 models. It appears that the SARS-CoV-2 models are simpler with respect to their long term behavior and thus display a simpler hierarchy with little dependencies compared to the Influenza-A models. Our results can support further development towards more complex SARS-CoV-2 models targeting the higher levels of the hierarchy.

The International Virus Bioinformatics Meeting 2020

The International Virus Bioinformatics Meeting 2020 was originally planned to take place in Bern, Switzerland, in March 2020. However, the COVID-19 pandemic put a spoke in the wheel of almost all conferences to be held in 2020. After moving the conference to 8-9 October 2020, we got hit by the second wave and finally decided at short notice to go fully online. On the other hand, the pandemic has made us even more aware of the importance of accelerating research in viral bioinformatics. Advances in bioinformatics have led to improved approaches to investigate viral infections and outbreaks. The International Virus Bioinformatics Meeting 2020 has attracted approximately 120 experts in virology and bioinformatics from all over the world to join the two-day virtual meeting. Despite concerns being raised that virtual meetings lack possibilities for face-to-face discussion, the participants from this small community created a highly interactive scientific environment, engaging in lively and inspiring discussions and suggesting new research directions and questions. The meeting featured five invited and twelve contributed talks, on the four main topics: (1) proteome and RNAome of RNA viruses, (2) viral metagenomics and ecology, (3) virus evolution and classification and (4) viral infections and immunology. Further, the meeting featured 20 oral poster presentations, all of which focused on specific areas of virus bioinformatics. This report summarizes the main research findings and highlights presented at the meeting.

On the interactions of the receptor-binding domain of SARS-CoV-1 and SARS-CoV-2 spike proteins with monoclonal antibodies and the receptor ACE2

A new betacoronavirus named SARS-CoV-2 has emerged as a new threat to global health and economy. A promising target for both diagnosis and therapeutics treatments of the new disease named COVID-19 is the coronavirus (CoV) spike (S) glycoprotein. By constant-pH Monte Carlo simulations and the PROCEEDpKa method, we have mapped the electrostatic epitopes for four monoclonal antibodies and the angiotensin-converting enzyme 2 (ACE2) on both SARS-CoV-1 and the new SARS-CoV-2 S receptor binding domain (RBD) proteins. We also calculated free energy of interactions and shown that the S RBD proteins from both SARS viruses binds to ACE2 with similar affinities. However, the affinity between the S RBD protein from the new SARS-CoV-2 and ACE2 is higher than for any studied antibody previously found complexed with SARS-CoV-1. Based on physical chemical analysis and free energies estimates, we can shed some light on the involved molecular recognition processes, their clinical aspects, the implications for drug developments, and suggest structural modifications on the CR3022 antibody that would improve its binding affinities for SARS-CoV-2 and contribute to address the ongoing international health crisis.

Origin and evolution of emerging Liao ning Virus (genus Seadornavirus, family Reoviridae)

Background

Liao ning virus (LNV) is a member of the genus Seadornavirus, family Reoviridae and has been isolated from kinds of vectors in Asia and Australia. However, there are no systematic studies describe the molecular genetic evolution and migration of LNVs. With the development of bioinformatics, viral genetic data combining the information of virus isolation time and locations could be integrated to infer the virus evolution and spread in nature.

Methods

Here, a phylogenetic and phylogeographic analysis using Bayesian Markov chain Monte Carlo simulations was conducted on the LNVs isolated from a variety of vectors during 1990–2014 to identify the evolution and migration patterns of LNVs.

Results

The results demonstrated that the LNV could be divided into 3 genotypes, of which genotype 1 mainly composed of LNVs isolated from Australia during 1990 to 2014 and the original LNV strain (LNV-NE97–31) isolated from Liaoning province in northern China in 1997, genotype 2 comprised of the isolates all from Xinjiang province in western China and genotype 3 consisted the isolates from Qinghai and Shanxi province of central China. LNVs emerged about 272 years ago and gradually evolved into three lineages in the order genotype 1, genotype 2 and genotype 3. Following phylogeographic analysis, it shows genotype 1 LNVs transmitted from Australia (113°E-153°E,10°S-42°S) to Liaoning province (118°E-125°E,38°N-43°N) in Northeast Asian continent then further spread across the central part of China to western China (75°E-95°E,35°N-50°N).

Conclusion

LNVs were initially isolated from Liaoning province of China in the Northeast Asia, however, the present study revealed that LNVs were first appeared in Australia in the South Pacific region and transmitted to mainland China then rapidly spread across China and evolved three different genotypes. The above results suggested that LNV had the characteristics of long-distance transmission and there were great genetic diversity existed in the LNV population. Notably, current information of 80 strains of LNVs are limited. It is of great importance to strengthen the surveillance of LNVs to explore its real origin in nature and monitoring of the LNVs’ population variation and maintain vigilance to avoid LNV breaking through the species barrier and further clarify its relationship to human and animal infection.

Improving risk assessment of the emergence of novel influenza A viruses by incorporating environmental surveillance

Reassortment is an evolutionary mechanism by which influenza A viruses (IAV) generate genetic novelty. Reassortment is an important driver of host jumps and is widespread according to retrospective surveillance studies. However, predicting the epidemiological risk of reassortant emergence in novel hosts from surveillance data remains challenging. IAV strains persist and co-occur in the environment, promoting co-infection during environmental transmission. These conditions offer opportunity to understand reassortant emergence in reservoir and spillover hosts. Specifically, environmental RNA could provide rich information for understanding the evolutionary ecology of segmented viruses, and transform our ability to quantify epidemiological risk to spillover hosts. However, significant challenges with recovering and interpreting genomic RNA from the environment have impeded progress towards predicting reassortant emergence from environmental surveillance data. We discuss how the fields of genomics, experimental ecology and epidemiological modelling are well positioned to address these challenges. Coupling quantitative disease models and natural transmission studies with new molecular technologies, such as deep-mutational scanning and single-virus sequencing of environmental samples, should dramatically improve our understanding of viral co-occurrence and reassortment. We define observable risk metrics for emerging molecular technologies and propose a conceptual research framework for improving accuracy and efficiency of risk prediction. This article is part of the theme issue 'Dynamic and integrative approaches to understanding pathogen spillover'.

The Third Annual Meeting of the European Virus Bioinformatics Center

The Third Annual Meeting of the European Virus Bioinformatics Center (EVBC) took place in Glasgow, United Kingdom, 28-29 March 2019. Virus bioinformatics has become central to virology research, and advances in bioinformatics have led to improved approaches to investigate viral infections and outbreaks, being successfully used to detect, control, and treat infections of humans and animals. This active field of research has attracted approximately 110 experts in virology and bioinformatics/computational biology from Europe and other parts of the world to attend the two-day meeting in Glasgow to increase scientific exchange between laboratory- and computer-based researchers. The meeting was held at the McIntyre Building of the University of Glasgow; a perfect location, as it was originally built to be a place for "rubbing your brains with those of other people", as Rector Stanley Baldwin described it. The goal of the meeting was to provide a meaningful and interactive scientific environment to promote discussion and collaboration and to inspire and suggest new research directions and questions. The meeting featured eight invited and twelve contributed talks, on the four main topics: (1) systems virology, (2) virus-host interactions and the virome, (3) virus classification and evolution and (4) epidemiology, surveillance and evolution. Further, the meeting featured 34 oral poster presentations, all of which focused on specific areas of virus bioinformatics. This report summarizes the main research findings and highlights presented at the meeting.

Interpreting Viral Deep Sequencing Data with GLUE

Using deep sequencing technologies such as Illumina's platform, it is possible to obtain reads from the viral RNA population revealing the viral genome diversity within a single host. A range of software tools and pipelines can transform raw deep sequencing reads into Sequence Alignment Mapping (SAM) files. We propose that interpretation tools should process these SAM files, directly translating individual reads to amino acids in order to extract statistics of interest such as the proportion of different amino acid residues at specific sites. This preserves per-read linkage between nucleotide variants at different positions within a codon location. The samReporter is a subsystem of the GLUE software toolkit which follows this direct read translation approach in its processing of SAM files. We test samReporter on a deep sequencing dataset obtained from a cohort of 241 UK HCV patients for whom prior treatment with direct-acting antivirals has failed; deep sequencing and resistance testing have been suggested to be of clinical use in this context. We compared the polymorphism interpretation results of the samReporter against an approach that does not preserve per-read linkage. We found that the samReporter was able to properly interpret the sequence data at resistance-associated locations in nine patients where the alternative approach was equivocal. In three cases, the samReporter confirmed that resistance or an atypical substitution was present at NS5A position 30. In three further cases, it confirmed that the sofosbuvir-resistant NS5B substitution S282T was absent. This suggests the direct read translation approach implemented is of value for interpreting viral deep sequencing data.

Proteolytic cleavage of host proteins by the Group IV viral proteases of Venezuelan equine encephalitis virus and Zika virus

The alphaviral nonstructural protein 2 (nsP2) cysteine proteases (EC 3.4.22.-) are essential for the proteolytic processing of the nonstructural (ns) polyprotein and are validated drug targets. A common secondary role of these proteases is to antagonize the effects of interferon (IFN). After delineating the cleavage site motif of the Venezuelan equine encephalitis virus (VEEV) nsP2 cysteine protease, we searched the human genome to identify host protein substrates. Here we identify a new host substrate of the VEEV nsP2 protease, human TRIM14, a component of the mitochondrial antiviral-signaling protein (MAVS) signalosome. Short stretches of homologous host-pathogen protein sequences (SSHHPS) are present in the nonstructural polyprotein and TRIM14. A 25-residue cyan-yellow fluorescent protein TRIM14 substrate was cleaved in vitro by the VEEV nsP2 protease and the cleavage site was confirmed by tandem mass spectrometry. A TRIM14 cleavage product also was found in VEEV-infected cell lysates. At least ten other Group IV (+)ssRNA viral proteases have been shown to cleave host proteins involved in generating the innate immune responses against viruses, suggesting that the integration of these short host protein sequences into the viral protease cleavage sites may represent an embedded mechanism of IFN antagonism. This interference mechanism shows several parallels with those of CRISPR/Cas9 and RNAi/RISC, but with a protease recognizing a protein sequence common to both the host and pathogen. The short host sequences embedded within the viral genome appear to be analogous to the short phage sequences found in a host's CRISPR spacer sequences. To test this algorithm, we applied it to another Group IV virus, Zika virus (ZIKV), and identified cleavage sites within human SFRP1 (secreted frizzled related protein 1), a retinal G_s alpha subunit, NT5M, and Forkhead box protein G1 (FOXG1) in vitro. Proteolytic cleavage of these proteins suggests a possible link between the protease and the virus-induced phenotype of ZIKV. The algorithm may have value for selecting cell lines and animal models that recapitulate virus-induced phenotypes, predicting host-range and susceptibility, selecting oncolytic viruses, identifying biomarkers, and de-risking live virus vaccines. Inhibitors of the proteases that utilize this mechanism may both inhibit viral replication and alleviate suppression of the innate immune responses.

SilentMutations (SIM): A tool for analyzing long-range RNA-RNA interactions in viral genomes and structured RNAs

•

We developed a tool to analyze the effect of multiple point mutations on the secondary structures of two interacting viral RNAs.

•

Our tool simulates destructive and compensatory mutants of two key regions from a single-stranded RNA.

•

The simulated mutants can be utilized for the combinatorial in vitro analysis of RNA–RNA interactions.

•

We predicted potential mutants for in vitro validation experiments of influenza A virus and hepatitis C virus interactions.

A single nucleotide change in the coding region can alter the amino acid sequence of a protein. In consequence, natural or artificial sequence changes in viral RNAs may have various effects not only on protein stability, function and structure but also on viral replication.

In recent decades, several tools have been developed to predict the effect of mutations in structured RNAs such as viral genomes or non-coding RNAs. Some tools use multiple point mutations and also take coding regions into account. However, none of these tools was designed to specifically simulate the effect of mutations on viral long-range interactions.

Here, we developed SilentMutations (SIM), an easy-to-use tool to analyze the effect of multiple point mutations on the secondary structures of two interacting viral RNAs. The tool can simulate disruptive and compensatory mutants of two interacting single-stranded RNAs. This allows a fast and accurate assessment of key regions potentially involved in functional long-range RNA–RNA interactions and will eventually help virologists and RNA-experts to design appropriate experiments.

SIM only requires two interacting single-stranded RNA regions as input. The output is a plain text file containing the most promising mutants and a graphical representation of all interactions.

We applied our tool on two experimentally validated influenza A virus and hepatitis C virus interactions and we were able to predict potential double mutants for in vitro validation experiments.

The source code and documentation of SIM are freely available at github.com/desiro/silentMutations.

Advances in Analyzing Virus-Induced Alterations of Host Cell Splicing

Alteration of host cell splicing is a common feature of many viral infections which is underappreciated because of the complexity and technical difficulty of studying alternative splicing (AS) regulation. Recent advances in RNA sequencing technologies revealed that up to several hundreds of host genes can show altered mRNA splicing upon viral infection. The observed changes in AS events can be either a direct consequence of viral manipulation of the host splicing machinery or result indirectly from the virus-induced innate immune response or cellular damage. Analysis at a higher resolution with single-cell RNAseq, and at a higher scale with the integration of multiple omics data sets in a systems biology perspective, will be needed to further comprehend this complex facet of virus-host interactions.

Bioinformatics Meets Virology: The European Virus Bioinformatics Center's Second Annual Meeting

The Second Annual Meeting of the European Virus Bioinformatics Center (EVBC), held in Utrecht, Netherlands, focused on computational approaches in virology, with topics including (but not limited to) virus discovery, diagnostics, (meta-)genomics, modeling, epidemiology, molecular structure, evolution, and viral ecology. The goals of the Second Annual Meeting were threefold: (i) to bring together virologists and bioinformaticians from across the academic, industrial, professional, and training sectors to share best practice; (ii) to provide a meaningful and interactive scientific environment to promote discussion and collaboration between students, postdoctoral fellows, and both new and established investigators; (iii) to inspire and suggest new research directions and questions. Approximately 120 researchers from around the world attended the Second Annual Meeting of the EVBC this year, including 15 renowned international speakers. This report presents an overview of new developments and novel research findings that emerged during the meeting.