Human leukocyte antigen (HLA) class I molecules play a crucial role in the development of a specific immune response to viral infections by presenting viral peptides at the cell surface where they will be further recognized by T cells. In the present manuscript, we explored whether HLA class I genotypes can be associated with the critical course of Coronavirus Disease-19 by searching possible connections between genotypes of deceased patients and their age at death. HLA-A, HLA-B, and HLA-C genotypes of n = 111 deceased patients with COVID-19 (Moscow, Russia) and n = 428 volunteers were identified with next-generation sequencing. Deceased patients were split into two groups according to age at the time of death: n = 26 adult patients aged below 60 and n = 85 elderly patients over 60. With the use of HLA class I genotypes, we developed a risk score (RS) which was associated with the ability to present severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) peptides by the HLA class I molecule set of an individual. The resulting RS was significantly higher in the group of deceased adults compared to elderly adults [p = 0.00348, area under the receiver operating characteristic curve (AUC ROC = 0.68)]. In particular, presence of HLA-A^*01:01 allele was associated with high risk, while HLA-A^*02:01 and HLA-A^*03:01 mainly contributed to low risk. The analysis of patients with homozygosity strongly highlighted these results: homozygosity by HLA-A^*01:01 accompanied early deaths, while only one HLA-A^*02:01 homozygote died before 60 years of age. Application of the constructed RS model to an independent Spanish patients cohort (n = 45) revealed that the score was also associated with the severity of the disease. The obtained results suggest the important role of HLA class I peptide presentation in the development of a specific immune response to COVID-19.

The interferon (IFN)-mediated innate immune response is the first line of defense against viruses. However, an IFN-stimulated gene, the adenosine deaminase acting on RNA 1 (ADAR1), favors the replication of several viruses. ADAR1 binds double-stranded RNA and converts adenosine to inosine by deamination. This form of editing makes duplex RNA unstable, thereby preventing IFN induction. To better understand how ADAR1 works at the cellular level, we generated cell lines that express exclusively either the IFN-inducible, cytoplasmic isoform ADAR1^p150, the constitutively expressed nuclear isoform ADAR1^p110, or no isoform. By comparing the transcriptome of these cell lines, we identified more than 150 polymerase II transcripts that are extensively edited, and we attributed most editing events to ADAR1^p150. Editing is focused on inverted transposable elements, located mainly within introns and untranslated regions, and predicted to form duplex RNA structures. Editing of these elements occurs also in primary human samples, and there is evidence for cross-species evolutionary conservation of editing patterns in primates and, to a lesser extent, in rodents. Whereas ADAR1^p150 rarely edits tightly encapsidated standard measles virus (MeV) genomes, it efficiently edits genomes with inverted repeats accidentally generated by a mutant MeV. We also show that immune activation occurs in fully ADAR1-deficient (ADAR1^KO) cells, restricting virus growth, and that complementation of these cells with ADAR1^p150 rescues virus growth and suppresses innate immunity activation. Finally, by knocking out either protein kinase R (PKR) or mitochondrial antiviral signaling protein (MAVS)—another protein controlling the response to duplex RNA—in ADAR1^KO cells, we show that PKR activation elicits a stronger antiviral response. Thus, ADAR1 prevents innate immunity activation by cellular transcripts that include extensive duplex RNA structures. The trade-off is that viruses take advantage of ADAR1 to elude innate immunity control. The innate immune response is a double-edged sword. It must protect the host from pathogens while avoiding accidental recognition of “self” molecular patterns, which can lead to autoimmune reactions. Double-stranded RNA is among the most potent inducers of cellular stress and interferon responses. We characterize here a mechanism that prevents autoimmune activation and show that an RNA virus, measles virus, can exploit it to elude innate immune responses. This mechanism relies on the enzyme adenosine deaminase acting on RNA 1 (ADAR1), which converts adenosine residues within duplex RNA structures to inosine. We identify duplex RNA structures in the 3′ untranslated regions of over 150 cellular transcripts and show that they are heavily edited in ADAR1-expressing cells. We detect the same type of editing in duplex RNA–forming defective genomes accidentally generated by measles virus. Loss of RNA editing causes strong innate immune responses and is detrimental to viral replication. Thus, by keeping the amount of duplex RNA in cells below an immune activation threshold, ADAR1 prevents autoimmunity while also favoring pathogens.

Host miRNAs are known as important regulators of virus replication and pathogenesis. They can interact with various viruses through several possible mechanisms including direct binding of viral RNA. Identification of human miRNAs involved in coronavirus-host interplay becomes important due to the ongoing COVID-19 pandemic. In this article we performed computational prediction of high-confidence direct interactions between miRNAs and seven human coronavirus RNAs. As a result, we identified six miRNAs (miR-21-3p, miR-195-5p, miR-16-5p, miR-3065-5p, miR-424-5p and miR-421) with high binding probability across all analyzed viruses. Further bioinformatic analysis of binding sites revealed high conservativity of miRNA binding regions within RNAs of human coronaviruses and their strains. In order to discover the entire miRNA-virus interplay we further analyzed lungs miRNome of SARS-CoV infected mice using publicly available miRNA sequencing data. We found that miRNA miR-21-3p has the largest probability of binding the human coronavirus RNAs and being dramatically up-regulated in mouse lungs during infection induced by SARS-CoV.

Interactions of the extracellular matrix (ECM) and cellular receptors constitute one of the crucial pathways involved in colorectal cancer progression and metastasis. With the use of bioinformatics analysis, we comprehensively evaluated the prognostic information concentrated in the genes from this pathway. First, we constructed a ECM–receptor regulatory network by integrating the transcription factor (TF) and 5’-isomiR interaction databases with mRNA/miRNA-seq data from The Cancer Genome Atlas Colon Adenocarcinoma (TCGA-COAD). Notably, one-third of interactions mediated by 5’-isomiRs was represented by noncanonical isomiRs (isomiRs, whose 5’-end sequence did not match with the canonical miRBase version). Then, exhaustive search-based feature selection was used to fit prognostic signatures composed of nodes from the network for overall survival prediction. Two reliable prognostic signatures were identified and validated on the independent The Cancer Genome Atlas Rectum Adenocarcinoma (TCGA-READ) cohort. The first signature was made up by six genes, directly involved in ECM–receptor interaction: AGRN, DAG1, FN1, ITGA5, THBS3, and TNC (concordance index 0.61, logrank test p = 0.0164, 3-years ROC AUC = 0.68). The second hybrid signature was composed of three regulators: hsa-miR-32-5p, NR1H2, and SNAI1 (concordance index 0.64, logrank test p = 0.0229, 3-years ROC AUC = 0.71). While hsa-miR-32-5p exclusively regulated ECM-related genes (COL1A2 and ITGA5), NR1H2 and SNAI1 also targeted other pathways (adhesion, cell cycle, and cell division). Concordant distributions of the respective risk scores across four stages of colorectal cancer and adjacent normal mucosa additionally confirmed reliability of the models.

Development of novel approaches for regulating the expression of angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2) is becoming increasingly important within the context of the ongoing COVID-19 pandemic since these enzymes play a crucial role in cell infection. In this work we searched for putative ACE2 and TMPRSS2 expression regulation networks mediated by various miRNA isoforms (isomiR) across different human organs using publicly available paired miRNA/mRNA-sequencing data from The Cancer Genome Atlas (TCGA) project. As a result, we identified several miRNA families targeting ACE2 and TMPRSS2 genes in multiple tissues. In particular, we found that lysine-specific demethylase 5B (JARID1B), encoded by the KDM5B gene, can indirectly affect ACE2 / TMPRSS2 expression by repressing transcription of hsa-let-7e / hsa-mir-125a and hsa-mir-141 / hsa-miR-200 miRNA families which are targeting these genes.

Rapidly appearing SARS-CoV-2 mutations can affect T cell epitopes, which can help the virus to evade either CD8 or CD4 T-cell responses. We developed T-cell COVID-19 Atlas (T-CoV, https://t-cov.hse.ru) – the comprehensive web portal, which allows one to analyze how SARS-CoV-2 mutations alter the presentation of viral peptides by HLA molecules. The data are presented for common virus variants and the most frequent HLA class I and class II alleles. Binding affinities of HLA molecules and viral peptides were assessed with accurate in silico methods. The obtained results highlight the importance of taking HLA alleles diversity into account: mutation-mediated alterations in HLA-peptide interactions were highly dependent on HLA alleles. For example, we found that the essential number of peptides tightly bound to HLA-B*07:02 in the reference Wuhan variant ceased to be tight binders for the Indian (Delta) and the UK (Alpha) variants. In summary, we believe that T-CoV will help researchers and clinicians to predict the susceptibility of individuals with different HLA genotypes to infection with variants of SARS-CoV-2 and/or forecast its severity.