IFITM3 rs12252 polymorphism association with COVID-19 severity and mortality in a Brazilian sample: an update and a meta-analysis.

This study investigated the association between the IFITM3 rs12252 polymorphism and the severity and mortality of COVID-19 in hospitalized Brazilian patients. A total of 102 COVID-19 patients were included, and the outcomes of interest were defined as death and the need for mechanical ventilation. Genotypes were assessed using Taqman probes. No significant associations were found between the rs12252 polymorphism and COVID-19 outcomes in the original sample, both for death and the need for mechanical ventilation. A meta-analysis, incorporating previous studies that used death as a severity indicator, revealed no association in the allelic and C-recessive models. However, due to the rarity of the T allele and its absence in the sample, further replication studies in larger and more diverse populations are needed to clarify the role of rs12252 in COVID-19 prognosis.

Bioinformatics and system biology approaches to determine the connection of SARS-CoV-2 infection and intrahepatic cholangiocarcinoma.

INTRODUCTION: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causal agent of coronavirus disease 2019 (COVID-19), has infected millions of individuals worldwide, which poses a severe threat to human health. COVID-19 is a systemic ailment affecting various tissues and organs, including the lungs and liver. Intrahepatic cholangiocarcinoma (ICC) is one of the most common liver cancer, and cancer patients are particularly at high risk of SARS-CoV-2 infection. Nonetheless, few studies have investigated the impact of COVID-19 on ICC patients. METHODS: With the methods of systems biology and bioinformatics, this study explored the link between COVID-19 and ICC, and searched for potential therapeutic drugs. RESULTS: This study identified a total of 70 common differentially expressed genes (DEGs) shared by both diseases, shedding light on their shared functionalities. Enrichment analysis pinpointed metabolism and immunity as the primary areas influenced by these common genes. Subsequently, through protein-protein interaction (PPI) network analysis, we identified SCD, ACSL5, ACAT2, HSD17B4, ALDOA, ACSS1, ACADSB, CYP51A1, PSAT1, and HKDC1 as hub genes. Additionally, 44 transcription factors (TFs) and 112 microRNAs (miRNAs) were forecasted to regulate the hub genes. Most importantly, several drug candidates (Periodate-oxidized adenosine, Desipramine, Quercetin, Perfluoroheptanoic acid, Tetrandrine, Pentadecafluorooctanoic acid, Benzo[a]pyrene, SARIN, Dorzolamide, 8-Bromo-cAMP) may prove effective in treating ICC and COVID-19. CONCLUSION: This study is expected to provide valuable references and potential drugs for future research and treatment of COVID-19 and ICC.

Identification of a membrane-associated element (MAE) in the C-terminal region of SARS-CoV-2 nsp6 that is essential for viral replication.

The coronavirus disease 2019 (COVID-19) pandemic, caused by the novel coronavirus severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2), has rapidly spread worldwide since its emergence in late 2019. Its ongoing evolution poses challenges for antiviral drug development. Coronavirus nsp6, a multiple-spanning transmembrane protein, participates in the biogenesis of the viral replication complex, which accommodates the viral replication-transcription complex. The roles of its structural domains in viral replication are not well studied. Herein, we predicted the structure of the SARS-CoV-2 nsp6 protein using AlphaFold2 and identified a highly folded C-terminal region (nsp6C) downstream of the transmembrane helices. The enhanced green fluorescent protein (EGFP)-fused nsp6C was found to cluster in the cytoplasm and associate with membranes. Functional mapping identified a minimal membrane-associated element (MAE) as the region from amino acids 237 to 276 (LGV-KLL), which is mainly composed of the α-helix H1 and the α-helix H2; the latter exhibits characteristics of an amphipathic helix (AH). Mutagenesis studies and membrane flotation experiments demonstrate that AH-like H2 is required for MAE-mediated membrane association. This MAE was functionally conserved across MERS-CoV, HCoV-OC43, HCoV-229E, HCoV-HKU1, and HCoV-NL63, all capable of mediating membrane association. In a SARS-CoV-2 replicon system, mutagenesis studies of H2 and replacements of H1 and H2 with their homologous counterparts demonstrated requirements of residues on both sides of the H2 and properly paired H1-H2 for MAE-mediated membrane association and viral replication. Notably, mutations I266A and K274A significantly attenuated viral replication without dramatically affecting membrane association, suggesting a dual role of the MAE in viral replication: mediating membrane association as well as participating in protein-protein interactions.IMPORTANCESevere acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) assembles a double-membrane vesicle (DMV) by the viral non-structural proteins for viral replication. Understanding the mechanisms of the DMV assembly is of paramount importance for antiviral development. Nsp6, a multiple-spanning transmembrane protein, plays an important role in the DMV biogenesis. Herein, we predicted the nsp6 structure of SARS-CoV-2 and other human coronaviruses using AlphaFold2 and identified a putative membrane-associated element (MAE) in the highly conserved C-terminal regions of nsp6. Experimentally, we verified a functionally conserved minimal MAE composed of two α-helices, the H1, and the amphipathic helix-like H2. Mutagenesis studies confirmed the requirement of H2 for MAE-mediated membrane association and viral replication and demonstrated a dual role of the MAE in viral replication, by mediating membrane association and participating in residue-specific interactions. This functionally conserved MAE may serve as a novel anti-viral target.

A measles-vectored vaccine candidate expressing prefusion-stabilized SARS-CoV-2 spike protein brought to phase I/II clinical trials: protection of African green monkeys from COVID-19 disease.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged at the end of 2019 and is responsible for the largest human pandemic in 100 years. Thirty-four vaccines are currently approved for use worldwide, and approximately 67% of the world population has received a complete primary series of one, yet countries are dealing with new waves of infections, variant viruses continue to emerge, and breakthrough infections are frequent secondary to waning immunity. Here, we evaluate a measles virus (MV)-vectored vaccine expressing a stabilized prefusion SARS-CoV-2 spike (S) protein (MV-ATU3-S2PΔF2A; V591) with demonstrated immunogenicity in mouse models (see companion article [J. Brunet, Z. Choucha, M. Gransagne, H. Tabbal, M.-W. Ku et al., J Virol 98:e01693-23, 2024, https://doi.org/10.1128/jvi.01693-23]) in an established African green monkey model of disease. Animals were vaccinated with V591 or the control vaccine (an equivalent MV-vectored vaccine with an irrelevant antigen) intramuscularly using a prime/boost schedule, followed by challenge with an early pandemic isolate of SARS-CoV-2 at 56 days post-vaccination. Pre-challenge, only V591-vaccinated animals developed S-specific antibodies that had virus-neutralizing activity as well as S-specific T cells. Following the challenge, V591-vaccinated animals had lower infectious virus and viral (v) RNA loads in mucosal secretions and stopped shedding virus in these secretions earlier. vRNA loads were lower in these animals in respiratory and gastrointestinal tract tissues at necropsy. This correlated with a lower disease burden in the lungs as quantified by PET/CT at early and late time points post-challenge and by pathological analysis at necropsy.IMPORTANCESevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the largest human pandemic in 100 years. Even though vaccines are currently available, countries are dealing with new waves of infections, variant viruses continue to emerge, breakthrough infections are frequent, and vaccine hesitancy persists. This study uses a safe and effective measles vaccine as a platform for vaccination against SARS-CoV-2. The candidate vaccine was used to vaccinate African green monkeys (AGMs). All vaccinated AGMs developed robust antigen-specific immune responses. After challenge, these AGMs produced less virus in mucosal secretions, for a shorter period, and had a reduced disease burden in the lungs compared to control animals. At necropsy, lower levels of viral RNA were detected in tissue samples from vaccinated animals, and the lungs of these animals lacked the histologic hallmarks of SARS-CoV-2 disease observed exclusively in the control AGMs.

Molecular basis of hippopotamus ACE2 binding to SARS-CoV-2.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a wide range of hosts, including hippopotami, which are semi-aquatic mammals and phylogenetically closely related to Cetacea. In this study, we characterized the binding properties of hippopotamus angiotensin-converting enzyme 2 (hiACE2) to the spike (S) protein receptor binding domains (RBDs) of the SARS-CoV-2 prototype (PT) and variants of concern (VOCs). Furthermore, the cryo-electron microscopy (cryo-EM) structure of the SARS-CoV-2 PT S protein complexed with hiACE2 was resolved. Structural and mutational analyses revealed that L30 and F83, which are specific to hiACE2, played a crucial role in the hiACE2/SARS-CoV-2 RBD interaction. In addition, comparative and structural analysis of ACE2 orthologs suggested that the cetaceans may have the potential to be infected by SARS-CoV-2. These results provide crucial molecular insights into the susceptibility of hippopotami to SARS-CoV-2 and suggest the potential risk of SARS-CoV-2 VOCs spillover and the necessity for surveillance. IMPORTANCE: The hippopotami are the first semi-aquatic artiodactyl mammals wherein SARS-CoV-2 infection has been reported. Exploration of the invasion mechanism of SARS-CoV-2 will provide important information for the surveillance of SARS-CoV-2 in hippopotami, as well as other semi-aquatic mammals and cetaceans. Here, we found that hippopotamus ACE2 (hiACE2) could efficiently bind to the RBDs of the SARS-CoV-2 prototype (PT) and variants of concern (VOCs) and facilitate the transduction of SARS-CoV-2 PT and VOCs pseudoviruses into hiACE2-expressing cells. The cryo-EM structure of the SARS-CoV-2 PT S protein complexed with hiACE2 elucidated a few critical residues in the RBD/hiACE2 interface, especially L30 and F83 of hiACE2 which are unique to hiACE2 and contributed to the decreased binding affinity to PT RBD compared to human ACE2. Our work provides insight into cross-species transmission and highlights the necessity for monitoring host jumps and spillover events on SARS-CoV-2 in semi-aquatic/aquatic mammals.

ISMI-VAE: A deep learning model for classifying disease cells using gene expression and SNV data.

Various studies have linked several diseases, including cancer and COVID-19, to single nucleotide variations (SNV). Although single-cell RNA sequencing (scRNA-seq) technology can provide SNV and gene expression data, few studies have integrated and analyzed these multimodal data. To address this issue, we introduce Interpretable Single-cell Multimodal Data Integration Based on Variational Autoencoder (ISMI-VAE). ISMI-VAE leverages latent variable models that utilize the characteristics of SNV and gene expression data to overcome high noise levels and uses deep learning techniques to integrate multimodal information, map them to a low-dimensional space, and classify disease cells. Moreover, ISMI-VAE introduces an attention mechanism to reflect feature importance and analyze genetic features that could potentially cause disease. Experimental results on three cancer data sets and one COVID-19 data set demonstrate that ISMI-VAE surpasses the baseline method in terms of both effectiveness and interpretability and can effectively identify disease-causing gene features.

A positive selection at binding site 501 in the B.1 lineage might have triggered the highly infectious sub-lineages of SARS-CoV-2.

The descendants of the B lineage are the most predominant variants among the SARS-CoV-2 virus due to the incorporation of new mutations augmenting the infectivity of the virus. There is a substantial increase in the transition transversion bias, nucleotide diversity and purifying selection on the spike protein in the descendants of the B lineage of the SARS-CoV-2 virus on a temporal scale. A strong bias for C-to-U substitutions is found in the genes encoding spike protein in this lineage. The positive selection has operated on the spike gene of B lineages and its sub-lineages. The B.1 lineage has undergone positive selection on site 501 of the receptor binding domain ultimately reflected in a key substitution N501Y in its three descendant lineages namely B.1.1.7, B.1.351 and P.1. The intensity of purifying selection on the multiple sites of the spike gene has increased substantially in the sub-lineages of B.1 in a timescale. The binding site 501 on the spike protein in B lineage is found to coevolve with other amino acid sites. This study sheds light on the evolutionary trajectory of the B lineage into highly infectious descendants in the recent past under the influence of positive and purifying selection exerted by natural immunity and vaccination of the host.

Conserved structures and dynamics in 5'-proximal regions of Betacoronavirus RNA genomes.

Betacoronaviruses are a genus within the Coronaviridae family of RNA viruses. They are capable of infecting vertebrates and causing epidemics as well as global pandemics in humans. Mitigating the threat posed by Betacoronaviruses requires an understanding of their molecular diversity. The development of novel antivirals hinges on understanding the key regulatory elements within the viral RNA genomes, in particular the 5'-proximal region, which is pivotal for viral protein synthesis. Using a combination of cryo-electron microscopy, atomic force microscopy, chemical probing, and computational modeling, we determined the structures of 5'-proximal regions in RNA genomes of Betacoronaviruses from four subgenera: OC43-CoV, SARS-CoV-2, MERS-CoV, and Rousettus bat-CoV. We obtained cryo-electron microscopy maps and determined atomic-resolution models for the stem-loop-5 (SL5) region at the translation start site and found that despite low sequence similarity and variable length of the helical elements it exhibits a remarkable structural conservation. Atomic force microscopy imaging revealed a common domain organization and a dynamic arrangement of structural elements connected with flexible linkers across all four Betacoronavirus subgenera. Together, these results reveal common features of a critical regulatory region shared between different Betacoronavirus RNA genomes, which may allow targeting of these RNAs by broad-spectrum antiviral therapeutics.

Global patterns of phylogenetic diversity and transmission of bat coronavirus.

Bats are reservoirs for multiple coronaviruses (CoVs). However, the phylogenetic diversity and transmission of global bat-borne CoVs remain poorly understood. Here, we performed a Bayesian phylogeographic analysis based on 3,594 bat CoV RdRp gene sequences to study the phylogenetic diversity and transmission of bat-borne CoVs and the underlying driving factors. We found that host-switching events occurred more frequently for α-CoVs than for ß-CoVs, and the latter was highly constrained by bat phylogeny. Bat species in the families Molossidae, Rhinolophidae, Miniopteridae, and Vespertilionidae had larger contributions to the cross-species transmission of bat CoVs. Regions of eastern and southern Africa, southern South America, Western Europe, and Southeast Asia were more frequently involved in cross-region transmission events of bat CoVs than other regions. Phylogenetic and geographic distances were the most important factors limiting CoV transmission. Bat taxa and global geographic hotspots associated with bat CoV phylogenetic diversity were identified, and bat species richness, mean annual temperature, global agricultural cropland, and human population density were strongly correlated with the phylogenetic diversity of bat CoVs. These findings provide insight into bat CoV evolution and ecological transmission among bat taxa. The identified hotspots of bat CoV evolution and transmission will guide early warnings of bat-borne CoV zoonotic diseases.

Human coronaviruses: The emergence of SARS-CoV-2 and management of COVID-19.

To date, a total of seven human coronaviruses (HCoVs) have been identified, all of which are important respiratory pathogens. Recently, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19), has led to a global pandemic causing millions of infections and deaths. Here, we summarize the discovery and fundamental virology of HCoVs, discuss their zoonotic transmission and highlight the weak species barrier of SARS-CoV-2. We also discuss the possible origins of SARS-CoV-2 variants of concern identified to date and discuss the experimental challenges in characterizing mutations of interest and propose methods to circumvent them. As the COVID-19 treatment and prevention landscape rapidly evolves, we summarize current therapeutics and vaccines, and their implications on SARS-CoV-2 variants. Finally, we explore how interspecies transmission of SARS-CoV-2 may drive the emergence of novel strains, how disease severity may evolve and how COVID-19 will likely continue to burden healthcare systems globally.

Genomic Comparisons of Alphacoronaviruses and Betacoronaviruses from Korean Bats.

Coronaviruses are well known as a diverse family of viruses that affect a wide range of hosts. Since the outbreak of severe acute respiratory syndrome, a variety of bat-associated coronaviruses have been identified in many countries. However, they do not represent all the specific geographic locations of their hosts. In this study, full-length genomes representing newly identified bat coronaviruses in South Korea were obtained using an RNA sequencing approach. The analysis, based on genome structure, conserved replicase domains, spike gene, and nucleocapsid genes revealed that bat Alphacoronaviruses are from three different viral species. Among them, the newly identified B20-97 strain may represent a new putative species, closely related to PEDV. In addition, the newly-identified MERS-related coronavirus exhibited shared genomic nucleotide identities of less than 76.4% with other Merbecoviruses. Recombination analysis and multiple alignments of spike and RBD amino acid sequences suggested that this strain underwent recombination events and could possibly use hDPP4 molecules as its receptor. The bat SARS-related CoV B20-50 is unlikely to be able to use hACE2 as its receptor and lack of an open reading frame in ORF8 gene region. Our results illustrate the diversity of coronaviruses in Korean bats and their evolutionary relationships. The evolution of the bat coronaviruses related ORF8 accessory gene is also discussed.

First detection of Hedgehog coronavirus 1 in Poland.

Hedgehogs are common in the majority of European countries and are known to host various pathogens, including viruses. The recent discovery of MERS-related coronaviruses (CoVs) in hedgehogs from Germany, France, the UK, China, and Italy suggests that hedgehogs may represent a wild reservoir of betacoronaviruses. This study reports the first detection and characterization of novel betacoronovirus, subgenus Merbecovirus in wild hedgehogs in Poland. The CoV RNA was detected in 10 out of 40 hedgehogs' rectal swabs and in 1 out of 18 samples of the lung. No viral RNA was identified in the duodenum and kidney. There was no significant relationship between clinical status, gender, hedgehogs' age, and coronaviral RNA detection. Phylogenetic analysis showed that CoVs detected in our study grouped together with other representatives of Hedgehog coronavirus 1 species identified in Western Europe. Our findings provide further evidence that hedgehogs are a natural reservoir of Merbecovirus. Considering the high mutation rate of CoVs and their potential for crossing interspecies barriers, the proper management of hedgehogs admitted to wildlife rehabilitation centres is needed. It cannot be excluded that merbecovirus strains detected in hedgehogs may recombine with other CoVs leading to new viruses with potential for interspecies transmission.

Analysis of SARS-CoV-2 synonymous codon usage evolution throughout the COVID-19 pandemic.

SARS-CoV-2, the seventh coronavirus known to infect humans, can cause severe life-threatening respiratory pathologies. To better understand SARS-CoV-2 evolution, genome-wide analyses have been made, including the general characterization of its codons usage profile. Here we present a bioinformatic analysis of the evolution of SARS-CoV-2 codon usage over time using complete genomes collected since December 2019. Our results show that SARS-CoV-2 codon usage pattern is antagonistic to, and it is getting farther away from that of the human host. Further, a selection of deoptimized codons over time, which was accompanied by a decrease in both the codon adaptation index and the effective number of codons, was observed. All together, these findings suggest that SARS-CoV-2 could be evolving, at least from the perspective of the synonymous codon usage, to become less pathogenic.

The unresolved question on COVID-19 virus origin: The three cards game?

The ongoing discussion about the real origin of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) feeds acrimonious debates. Where did SARS-CoV-2 come from? Was SARS-CoV-2 transmitted in the wild from an animal to a person before exploding in Wuhan or was it an engineered virus that escaped from research or a laboratory in Wuhan? Right now, we still don't know enough whether SARS-CoV-2 is human-made or not, and lab-leak theories remain essentially speculative. Many recent studies have pointed out several plausible scenarios. Anyhow, currently, even if suspicions by some about the possibility of lab-leak hypothesis still remain, the consensus view is that the pandemic probably started from a natural source and, to determine the real origin of the SARS-CoV-2 virus, further research is needed.

LinearTurboFold: Linear-time global prediction of conserved structures for RNA homologs with applications to SARS-CoV-2.

The constant emergence of COVID-19 variants reduces the effectiveness of existing vaccines and test kits. Therefore, it is critical to identify conserved structures in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes as potential targets for variant-proof diagnostics and therapeutics. However, the algorithms to predict these conserved structures, which simultaneously fold and align multiple RNA homologs, scale at best cubically with sequence length and are thus infeasible for coronaviruses, which possess the longest genomes (∼30,000 nt) among RNA viruses. As a result, existing efforts on modeling SARS-CoV-2 structures resort to single-sequence folding as well as local folding methods with short window sizes, which inevitably neglect long-range interactions that are crucial in RNA functions. Here we present LinearTurboFold, an efficient algorithm for folding RNA homologs that scales linearly with sequence length, enabling unprecedented global structural analysis on SARS-CoV-2. Surprisingly, on a group of SARS-CoV-2 and SARS-related genomes, LinearTurboFold's purely in silico prediction not only is close to experimentally guided models for local structures, but also goes far beyond them by capturing the end-to-end pairs between 5' and 3' untranslated regions (UTRs) (∼29,800 nt apart) that match perfectly with a purely experimental work. Furthermore, LinearTurboFold identifies undiscovered conserved structures and conserved accessible regions as potential targets for designing efficient and mutation-insensitive small-molecule drugs, antisense oligonucleotides, small interfering RNAs (siRNAs), CRISPR-Cas13 guide RNAs, and RT-PCR primers. LinearTurboFold is a general technique that can also be applied to other RNA viruses and full-length genome studies and will be a useful tool in fighting the current and future pandemics.

Longitudinal monitoring in Cambodia suggests higher circulation of alpha and betacoronaviruses in juvenile and immature bats of three species.

Recent studies suggest that coronaviruses circulate widely in Southeast Asian bat species and that the progenitors of the SARS-Cov-2 virus could have originated in rhinolophid bats in the region. Our objective was to assess the diversity and circulation patterns of coronavirus in several bat species in Southeast Asia. We undertook monthly live-capture sessions and sampling in Cambodia over 17 months to cover all phases of the annual reproduction cycle of bats and test specifically the association between their age and CoV infection status. We additionally examined current information on the reproductive phenology of Rhinolophus and other bat species presently known to occur in mainland southeast China, Vietnam, Laos and Cambodia. Results from our longitudinal monitoring (573 bats belonging to 8 species) showed an overall proportion of positive PCR tests for CoV of 4.2% (24/573) in cave-dwelling bats from Kampot and 4.75% (22/463) in flying-foxes from Kandal. Phylogenetic analysis showed that the PCR amplicon sequences of CoVs (n = 46) obtained clustered in Alphacoronavirus and Betacoronavirus. Interestingly, Hipposideros larvatus sensu lato harbored viruses from both genera. Our results suggest an association between positive detections of coronaviruses and juvenile and immature bats in Cambodia (OR = 3.24 [1.46-7.76], p = 0.005). Since the limited data presently available from literature review indicates that reproduction is largely synchronized among rhinolophid and hipposiderid bats in our study region, particularly in its more seasonal portions (above 16° N), this may lead to seasonal patterns in CoV circulation. Overall, our study suggests that surveillance of CoV in insectivorous bat species in Southeast Asia, including SARS-CoV-related coronaviruses in rhinolophid bats, could be targeted from June to October for species exhibiting high proportions of juveniles and immatures during these months. It also highlights the need to develop long-term longitudinal surveys of bats and improve our understanding of their ecology in the region, for both biodiversity conservation and public health reasons.

Best Molecular Tools to Investigate Coronavirus Diversity in Mammals: A Comparison.

Coronaviruses (CoVs) are widespread and highly diversified in wildlife and domestic mammals and can emerge as zoonotic or epizootic pathogens and consequently host shift from these reservoirs, highlighting the importance of veterinary surveillance. All genera can be found in mammals, with α and ß showing the highest frequency and diversification. The aims of this study were to review the literature for features of CoV surveillance in animals, to test widely used molecular protocols, and to identify the most effective one in terms of spectrum and sensitivity. We combined a literature review with analyses in silico and in vitro using viral strains and archive field samples. We found that most protocols defined as pan-coronavirus are strongly biased towards α- and ß-CoVs and show medium-low sensitivity. The best results were observed using our new protocol, showing LoD 100 PFU/mL for SARS-CoV-2, 50 TCID50/mL for CaCoV, 0.39 TCID50/mL for BoCoV, and 9 ± 1 log2 ×10-5 HA for IBV. The protocol successfully confirmed the positivity for a broad range of CoVs in 30/30 field samples. Our study points out that pan-CoV surveillance in mammals could be strongly improved in sensitivity and spectrum and propose the application of a new RT-PCR assay, which is able to detect CoVs from all four genera, with an optimal sensitivity for α-, ß-, and γ-.

Utilizing the VirIdAl Pipeline to Search for Viruses in the Metagenomic Data of Bat Samples.

According to various estimates, only a small percentage of existing viruses have been discovered, naturally much less being represented in the genomic databases. High-throughput sequencing technologies develop rapidly, empowering large-scale screening of various biological samples for the presence of pathogen-associated nucleotide sequences, but many organisms are yet to be attributed specific loci for identification. This problem particularly impedes viral screening, due to vast heterogeneity in viral genomes. In this paper, we present a new bioinformatic pipeline, VirIdAl, for detecting and identifying viral pathogens in sequencing data. We also demonstrate the utility of the new software by applying it to viral screening of the feces of bats collected in the Moscow region, which revealed a significant variety of viruses associated with bats, insects, plants, and protozoa. The presence of alpha and beta coronavirus reads, including the MERS-like bat virus, deserves a special mention, as it once again indicates that bats are indeed reservoirs for many viral pathogens. In addition, it was shown that alignment-based methods were unable to identify the taxon for a large proportion of reads, and we additionally applied other approaches, showing that they can further reveal the presence of viral agents in sequencing data. However, the incompleteness of viral databases remains a significant problem in the studies of viral diversity, and therefore necessitates the use of combined approaches, including those based on machine learning methods.

Neutral evolution test of the spike protein of SARS-CoV-2 and its implications in the binding to ACE2.

As the SARS-CoV-2 has spread and the pandemic has dragged on, the virus continued to evolve rapidly resulting in the emergence of new highly transmissible variants that can be of public health concern. The evolutionary mechanisms that drove this rapid diversity are not well understood but neutral evolution should open the first insight. The neutral theory of evolution states that most mutations in the nucleic acid sequences are random and they can be fixed or disappear by purifying selection. Herein, we performed a neutrality test to better understand the selective pressures exerted over SARS-CoV-2 spike protein from homologue proteins of Betacoronavirus, as well as to the spikes from human clinical isolates of the virus. Specifically, Tyr and Asn have higher occurrence rates on the Receptor Binding Domain (RBD) and in the overall sequence of spike proteins of Betacoronavirus, whereas His and Arg have lower occurrence rates. The in vivo evolutionary phenomenon of SARS-CoV-2 shows that Glu, Lys, Phe, and Val have the highest probability of occurrence in the emergent viral particles. Amino acids that have higher occurrence than the expected by the neutral control, are favorable and are fixed in the sequence while the ones that have lower occurrence than expected, influence the stability and/or functionality of the protein. Our results show that most unique mutations either for SARS-CoV-2 or its variants of health concern are under selective pressures, which could be related either to the evasion of the immune system, increasing the virus' fitness or altering protein - protein interactions with host proteins. We explored the consequences of those selected mutations in the structure and protein - protein interaction with the receptor. Altogether all these forces have shaped the spike protein and the continually evolving variants.