Deciphering Abnormal Platelet Subpopulations in COVID-19, Sepsis and Systemic Lupus Erythematosus through Machine Learning and Single-Cell Transcriptomics.

This study focuses on understanding the transcriptional heterogeneity of activated platelets and its impact on diseases such as sepsis, COVID-19, and systemic lupus erythematosus (SLE). Recognizing the limited knowledge in this area, our research aims to dissect the complex transcriptional profiles of activated platelets to aid in developing targeted therapies for abnormal and pathogenic platelet subtypes. We analyzed single-cell transcriptional profiles from 47,977 platelets derived from 413 samples of patients with these diseases, utilizing Deep Neural Network (DNN) and eXtreme Gradient Boosting (XGB) to distinguish transcriptomic signatures predictive of fatal or survival outcomes. Our approach included source data annotations and platelet markers, along with SingleR and Seurat for comprehensive profiling. Additionally, we employed Uniform Manifold Approximation and Projection (UMAP) for effective dimensionality reduction and visualization, aiding in the identification of various platelet subtypes and their relation to disease severity and patient outcomes. Our results highlighted distinct platelet subpopulations that correlate with disease severity, revealing that changes in platelet transcription patterns can intensify endotheliopathy, increasing the risk of coagulation in fatal cases. Moreover, these changes may impact lymphocyte function, indicating a more extensive role for platelets in inflammatory and immune responses. This study identifies crucial biomarkers of platelet heterogeneity in serious health conditions, paving the way for innovative therapeutic approaches targeting platelet activation, which could improve patient outcomes in diseases characterized by altered platelet function.

B- and T-Lymphocyte Attenuator in Systemic Lupus Erythematosus Disease Pathogenesis.

B- and T-lymphocyte attenuator (BTLA; CD272) is an immunoglobulin superfamily member and part of a family of checkpoint inhibitory receptors that negatively regulate immune cell activation. The natural ligand for BTLA is herpes virus entry mediator (HVEM; TNFRSF14), and binding of HVEM to BTLA leads to attenuation of lymphocyte activation. In this study, we evaluated the role of BTLA and HVEM expression in the pathogenesis of systemic lupus erythematosus (SLE), a multisystem autoimmune disease. Peripheral blood mononuclear cells from healthy volunteers (N = 7) were evaluated by mass cytometry by time-of-flight to establish baseline expression of BTLA and HVEM on human lymphocytes compared with patients with SLE during a self-reported flare (N = 5). High levels of BTLA protein were observed on B cells, CD4+, and CD8+ T cells, and plasmacytoid dendritic cells in healthy participants. HVEM protein levels were lower in patients with SLE compared with healthy participants, while BTLA levels were similar between SLE and healthy groups. Correlations of BTLA-HVEM hub genes' expression with patient and disease characteristics were also analyzed using whole blood gene expression data from patients with SLE (N = 1,760) and compared with healthy participants (N = 60). HVEM, being one of the SLE-associated genes, showed an exceptionally strong negative association with disease activity. Several other genes in the BTLA-HVEM signaling network were strongly (negative or positive) correlated, while BTLA had a low association with disease activity. Collectively, these data provide a clinical rationale for targeting BTLA with an agonist in SLE patients with low HVEM expression.

Unusual structural and functional features of TpLRR/BspA-like LRR proteins.

Leucine Rich Repeat (LRR) domains, are present in hundreds of thousands of proteins across all kingdoms of life and are typically involved in protein-protein interactions and ligand recognition. LRR domains are classified into eight classes and when examined in three dimensions seven, of them form curved solenoid-like super-helices, also described as toruses, with a beta sheet on the concave (inside) and stacked alpha-helices on the convex (outside) of the torus. Here we present an overview of the least characterized 8th class of LRR proteins, the TpLRR-like LRRs, named after the Treponema pallidum protein Tp0225. Proteins from the TpLRR class differ from the proteins in all other known LRR classes by having a flipped curvature, with the beta sheet on the convex side of the torus and irregular secondary structure instead of helices on the opposite, now concave site. TpLRR proteins also present highly divergent sequence pattern of individual repeats and can associate with specific types of additional domains. Several of the characterized proteins from this class, specifically the BspA-like proteins, were found in human bacterial and protozoan pathogens, playing an important role in the interactions between the pathogens and the host immune system. In this paper we surveyed all existing experimental structures and selected AlphaFold models of the best-known proteins containing this class of LRR repeats, analyzing the relation between the pattern of conserved residues, specific structural features and functions of these proteins.

US National Institutes of Health Prioritization of SARS-CoV-2 Variants.

Since late 2020, SARS-CoV-2 variants have regularly emerged with competitive and phenotypic differences from previously circulating strains, sometimes with the potential to escape from immunity produced by prior exposure and infection. The Early Detection group is one of the constituent groups of the US National Institutes of Health National Institute of Allergy and Infectious Diseases SARS-CoV-2 Assessment of Viral Evolution program. The group uses bioinformatic methods to monitor the emergence, spread, and potential phenotypic properties of emerging and circulating strains to identify the most relevant variants for experimental groups within the program to phenotypically characterize. Since April 2021, the group has prioritized variants monthly. Prioritization successes include rapidly identifying most major variants of SARS-CoV-2 and providing experimental groups within the National Institutes of Health program easy access to regularly updated information on the recent evolution and epidemiology of SARS-CoV-2 that can be used to guide phenotypic investigations.

Gender differences in smoking-induced changes in the tumor immune microenvironment.

Both gender and smoking are correlated with prevalence and outcomes in many types of cancers. Tobacco smoke is a known carcinogen through its genotoxicity but can also affect cancer progression through its effect on the immune system. In this study, we aim to evaluate the hypothesis that the effects of smoking on the tumor immune microenvironment will be influenced differently by gender using large-scale analysis of publicly available cancer datasets. We used The Cancer Genomic Atlas (TCGA) datasets (n = 2724) to analyze effects of smoking on different cancer immune subtypes and the relative abundance of immune cell types between male and female cancer patients. We further validated our results by analyzing additional datasets, including Expression Project for Oncology (expO) bulk RNA-seq dataset (n = 1118) and single-cell RNA-seq dataset (n = 14). Results of our study indicate that in female patients, two immune subtypes, C1 and C2, are respectively over and under abundant in smokers vs. never smokers. In males, the only significant difference is underabundance of the C6 subtype in smokers. We identified gender-specific differences in the population of immune cell types between smokers and never smokers in all TCGA and expO cancer types. Increased plasma cell population was identified as the most consistent feature distinguishing smokers and never smokers, especially in current female smokers based on both TCGA and expO data. Our analysis of existing single-cell RNA-seq data further revealed that smoking differentially affects the gene expression profile of cancer patients based on the immune cell type and gender. In our analysis, female and male smokers show different smoking-induced patterns of immune cells in tumor microenvironment. Besides, our results suggest cancer tissues directly exposed to tobacco smoke undergo the most significant changes, but all other tissue types are affected as well. Findings of current study also indicate that changes in the populations of plasma cells and their correlations to survival outcomes are stronger in female current smokers, with implications for cancer immunotherapy of women smokers. In conclusion, results of this study can be used to develop personalized treatment plans for cancer patients who smoke, particularly women smokers, taking into account the unique immune cell profile of their tumors.

ViralVar: A Web Tool for Multilevel Visualization of SARS-CoV-2 Genomes.

The unprecedented growth of publicly available SARS-CoV-2 genome sequence data has increased the demand for effective and accessible SARS-CoV-2 data analysis and visualization tools. The majority of the currently available tools either require computational expertise to deploy them or limit user input to preselected subsets of SARS-CoV-2 genomes. To address these limitations, we developed ViralVar, a publicly available, point-and-click webtool that gives users the freedom to investigate and visualize user-selected subsets of SARS-CoV-2 genomes obtained from the GISAID public database. ViralVar has two primary features that enable: (1) the visualization of the spatiotemporal dynamics of SARS-CoV-2 lineages and (2) a structural/functional analysis of genomic mutations. As proof-of-principle, ViralVar was used to explore the evolution of the SARS-CoV-2 pandemic in the USA in pediatric, adult, and elderly populations (n > 1.7 million genomes). Whereas the spatiotemporal dynamics of the variants did not differ between these age groups, several USA-specific sublineages arose relative to the rest of the world. Our development and utilization of ViralVar to provide insights on the evolution of SARS-CoV-2 in the USA demonstrates the importance of developing accessible tools to facilitate and accelerate the large-scale surveillance of circulating pathogens.

Resistin Concentration in Early Sepsis and All-Cause Mortality at a Safety-Net Hospital in Riverside County.

Background: Sepsis mortality has remained unchanged for greater than a decade, and early recognition continues to be the most important factor in mortality outcome. Plasma resistin concentration is increased in sepsis, but its mechanism and clinical relevance is unclear. As one function, resistin interacts with toll-like receptor 4 in competition with lipopolysaccharide, a main component of the gram-negative bacterial cell wall. It is not known if the type of infection leading to sepsis influences resistin production. The objective of this study was to investigate whether 1) early plasma resistin concentration can predict mortality, 2) elevated plasma resistin concentration is associated with clinical disease severity scores, such as SOFA, mSOFA and APACHE II, and 3) plasma resistin concentrations differ between gram negative versus other etiologies of sepsis. Methods: This was an exploratory study in the framework of a prospective observational design. Peripheral venous blood samples were obtained from subjects admitted to the intensive care unit at clinical recognition of sepsis (0 hour) and at 6 and 24 hours. Vasopressor utilization was not a requirement for inclusion. Plasma was analyzed for resistin concentration by ELISA. Cytokine concentrations including IL-6, IL-8, and IL-10 were determined by cytokine bead array. Cytokine data were evaluated against publicly available sepsis RNA expression datasets to compare protein versus RNA expression levels in predicting clinical disease state. Clinical data were collected from electronic health records for clinical severity index calculations and context for interpretation of resistin and cytokine concentrations. Subjects were followed up to 60 days, or until death, whichever came first. Statistical analysis was completed with R package and SPSS software. Results: Resistin levels were elevated in subjects admitted to the intensive care unit with sepsis. Four-hundred subjects were screened with 45 subjects included in the final analysis. Thirteen of 45 patients were non-survivors. Mortality within 60 days correlated with significantly higher resistin concentrations than in survivors. A resistin concentration of >126 ng/mL at clinical recognition of sepsis and >197 ng/mL within the first 24 hours were associated with mortality within 60 days with an area under the curve of 0.82 and 0.88, respectively. Most subjects with resistin concentration greater than these threshold values were deceased prior to 30 days. Resistin concentrations correlated with SOFA, mSOFA, and APACHE II scores in addition to having association with increases in inflammatory and sepsis biomarkers. These associations were validated with analysis of RNA expression datasets. Conclusion: Plasma resistin concentrations of >126 ng/mL at clinical recognition of sepsis and >197 ng/mL within the first 24 hours of clinical sepsis recognition are associated with all-cause mortality. Resistin concentration within this timeframe also has comparable mortality association to well-validated clinical severity indices of SOFA, mSOFA, and APACHE II scores.

Increased Frequency of Indels in Hypervariable Regions of SARS-CoV-2 Proteins-A Possible Signature of Adaptive Selection.

Most attention in the surveillance of evolving SARS-CoV-2 genome has been centered on nucleotide substitutions in the spike glycoprotein. We show that, as the pandemic extends into its second year, the numbers and ratio of genomes with in-frame insertions and deletions (indels) increases significantly, especially among the variants of concern (VOCs). Monitoring of the SARS-CoV-2 genome evolution shows that co-occurrence (i.e., highly correlated presence) of indels, especially deletions on spike N-terminal domain and non-structural protein 6 (NSP6) is a shared feature in several VOCs such as Alpha, Beta, Delta, and Omicron. Indels distribution is correlated with spike mutations associated with immune escape and growth in the number of genomes with indels coincides with the increasing population resistance due to vaccination and previous infections. Indels occur most frequently in the spike, but also in other proteins, especially those involved in interactions with the host immune system. We also showed that indels concentrate in regions of individual SARS-CoV-2 proteins known as hypervariable regions (HVRs) that are mostly located in specific loop regions. Structural analysis suggests that indels remodel viral proteins' surfaces at common epitopes and interaction interfaces, affecting the virus' interactions with host proteins. We hypothesize that the increased frequency of indels, the non-random distribution of them and their independent co-occurrence in several VOCs is another mechanism of response to elevated global population immunity.

What the protein data bank tells us about the evolutionary conservation of protein conformational diversity.

Proteins sample a multitude of different conformations by undergoing small- and large-scale conformational changes that are often intrinsic to their functions. Information about these changes is often captured in the Protein Data Bank by the apparently redundant deposition of independent structural solutions of identical proteins. Here, we mine these data to examine the conservation of large-scale conformational changes between homologous proteins. This is important for both practical reasons, such as predicting alternative conformations of a protein by comparative modeling, and conceptual reasons, such as understanding the extent of conservation of different features in evolution. To study this question, we introduce a novel approach to compare conformational changes between proteins by the comparison of their difference distance maps (DDMs). We found that proteins undergoing similar conformational changes have similar DDMs and that this similarity could be quantified by the correlation between the DDMs. By comparing the DDMs of homologous protein pairs, we found that large-scale conformational changes show a high level of conservation across a broad range of sequence identities. This shows that conformational space is usually conserved between homologs, even relatively distant ones.

A Genomic Island of Vibrio cholerae Encodes a Three-Component Cytotoxin with Monomer and Protomer Forms Structurally Similar to Alpha-Pore-Forming Toxins.

Alpha-pore-forming toxins (α-PFTs) are secreted by many species of bacteria, including Escherichia coli, Aeromonas hydrophila, and Bacillus thuringiensis, as part of their arsenal of virulence factors, and are often cytotoxic. In particular, for α-PFTs, the membrane-spanning channel they form is composed of hydrophobic α-helices. These toxins oligomerize at the surface of target cells and transition from a soluble to a protomer state in which they expose their hydrophobic regions and insert into the membrane to form a pore. The pores may be composed of homooligomers of one component or heterooligomers with two or three components, resulting in bi- or tripartite toxins. The multicomponent α-PFTs are often expressed from a single operon. Recently, motility-associated killing factor A (MakA), an α-PFT, was discovered in Vibrio cholerae. We report that makA is found on the V. cholerae GI-10 genomic island within an operon containing genes for two other potential α-PFTs, MakB and MakE. We determined the X-ray crystal structures for MakA, MakB, and MakE and demonstrated that all three are structurally related to the α-PFT family in the soluble state, and we modeled their protomer state based on the α-PFT AhlB from A. hydrophila. We found that MakA alone is cytotoxic at micromolar concentrations. However, combining MakA with MakB and MakE is cytotoxic at nanomolar concentrations, with specificity for J774 macrophage cells. Our data suggest that MakA, -B, and -E are α-PFTs that potentially act as a tripartite pore-forming toxin with specificity for phagocytic cells. IMPORTANCE The bacterium Vibrio cholerae causes gastrointestinal, wound, and skin infections. The motility-associated killing factor A (MakA) was recently shown to be cytotoxic against colon, prostate, and other cancer cells. However, at the outset of this study, the capacity of MakA to damage cells in combination with other Mak proteins encoded in the same operon had not been elucidated. We determined the structures of three Mak proteins and established that they are structurally related to the α-PFTs. Compared to MakA alone, the combination of all three toxins was more potent specifically in mouse macrophages. This study highlights the idea that the Mak toxins are selectively cytotoxic and thus may function as a tripartite toxin with cell type specificity.

Defining the risk of SARS-CoV-2 variants on immune protection.

The global emergence of many severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants jeopardizes the protective antiviral immunity induced after infection or vaccination. To address the public health threat caused by the increasing SARS-CoV-2 genomic diversity, the National Institute of Allergy and Infectious Diseases within the National Institutes of Health established the SARS-CoV-2 Assessment of Viral Evolution (SAVE) programme. This effort was designed to provide a real-time risk assessment of SARS-CoV-2 variants that could potentially affect the transmission, virulence, and resistance to infection- and vaccine-induced immunity. The SAVE programme is a critical data-generating component of the US Government SARS-CoV-2 Interagency Group to assess implications of SARS-CoV-2 variants on diagnostics, vaccines and therapeutics, and for communicating public health risk. Here we describe the coordinated approach used to identify and curate data about emerging variants, their impact on immunity and effects on vaccine protection using animal models. We report the development of reagents, methodologies, models and notable findings facilitated by this collaborative approach and identify future challenges. This programme is a template for the response to rapidly evolving pathogens with pandemic potential by monitoring viral evolution in the human population to identify variants that could reduce the effectiveness of countermeasures.

Novel putative polyethylene terephthalate (PET) plastic degrading enzymes from the environmental metagenome.

Several plastic degrading enzymes have been described in the literature, most notably PETases that are capable of hydrolyzing polyethylene terephthalate (PET) plastic. One of them, the PETase from Ideonella sakaiensis, a bacterium isolated from environmental samples within a PET bottle recycling site, was a subject of extensive studies. To test how widespread PETase functionality is in other bacterial communities, we used a cascade of BLAST searches in the JGI metagenomic datasets and showed that close homologs of I. sakaiensis PETase can also be found in other metagenomic environmental samples from both human-affected and relatively pristine sites. To confirm their classification as putative PETases, we verified that the newly identified proteins have the PETase sequence signatures common to known PETases and that phylogenetic analyses group them with the experimentally characterized PETases. Additionally, docking analysis was performed in order to further confirm the functional assignment of the putative environmental PETases.

Dynamic changes in human single-cell transcriptional signatures during fatal sepsis.

Systemic infections, especially in patients with chronic diseases, may result in sepsis: an explosive, uncoordinated immune response that can lead to multisystem organ failure with a high mortality rate. Patients with similar clinical phenotypes or sepsis biomarker expression upon diagnosis may have different outcomes, suggesting that the dynamics of sepsis is critical in disease progression. A within-subject study of patients with Gram-negative bacterial sepsis with surviving and fatal outcomes was designed and single-cell transcriptomic analyses of peripheral blood mononuclear cells (PBMC) collected during the critical period between sepsis diagnosis and 6 h were performed. The single-cell observations in the study are consistent with trends from public datasets but also identify dynamic effects in individual cell subsets that change within hours. It is shown that platelet and erythroid precursor responses are drivers of fatal sepsis, with transcriptional signatures that are shared with severe COVID-19 disease. It is also shown that hypoxic stress is a driving factor in immune and metabolic dysfunction of monocytes and erythroid precursors. Last, the data support CD52 as a prognostic biomarker and therapeutic target for sepsis as its expression dynamically increases in lymphocytes and correlates with improved sepsis outcomes. In conclusion, this study describes the first single-cell study that analyzed short-term temporal changes in the immune cell populations and their characteristics in surviving or fatal sepsis. Tracking temporal expression changes in specific cell types could lead to more accurate predictions of sepsis outcomes and identify molecular biomarkers and pathways that could be therapeutically controlled to improve the sepsis trajectory toward better outcomes.

The interplay of SARS-CoV-2 evolution and constraints imposed by the structure and functionality of its proteins.

The unprecedented pace of the sequencing of the SARS-CoV-2 virus genomes provides us with unique information about the genetic changes in a single pathogen during ongoing pandemic. By the analysis of close to 200,000 genomes we show that the patterns of the SARS-CoV-2 virus mutations along its genome are closely correlated with the structural and functional features of the encoded proteins. Requirements of foldability of proteins' 3D structures and the conservation of their key functional regions, such as protein-protein interaction interfaces, are the dominant factors driving evolutionary selection in protein-coding genes. At the same time, avoidance of the host immunity leads to the abundance of mutations in other regions, resulting in high variability of the missense mutation rate along the genome. "Unexplained" peaks and valleys in the mutation rate provide hints on function for yet uncharacterized genomic regions and specific protein structural and functional features they code for. Some of these observations have immediate practical implications for the selection of target regions for PCR-based COVID-19 tests and for evaluating the risk of mutations in epitopes targeted by specific antibodies and vaccine design strategies.

ModFlex: Towards Function Focused Protein Modeling.

There is a wide, and continuously widening, gap between the number of proteins known only by their amino acid sequence versus those structurally characterized by direct experiment. To close this gap, we mostly rely on homology-based inference and modeling to reason about the structures of the uncharacterized proteins by using structures of homologous proteins as templates. With the rapidly growing size of the Protein Data Bank, there are often multiple choices of templates, including multiple sets of coordinates from the same protein. The substantial conformational differences observed between different experimental structures of the same protein often reflect function related structural flexibility. Thus, depending on the questions being asked, using distant homologs, or coordinate sets with lower resolution but solved in the appropriate functional form, as templates may be more informative. The ModFlex server (https://modflex.org/) addresses this seldom mentioned gap in the standard homology modeling approach by providing the user with an interface with multiple options and tools to select the most relevant template and explore the range of structural diversity in the available templates. ModFlex is closely integrated with a range of other programs and servers developed in our group for the analysis and visualization of protein structural flexibility and divergence.

Identification and phylogenetic analysis of RNA binding domain abundant in apicomplexans or RAP proteins.

The RNA binding domain abundant in apicomplexans (RAP) is a protein domain identified in a diverse group of proteins, called RAP proteins, many of which have been shown to be involved in RNA binding. To understand the expansion and potential function of the RAP proteins, we conducted a hidden Markov model based screen among the proteomes of 54 eukaryotes, 17 bacteria and 12 archaea. We demonstrated that the domain is present in closely and distantly related organisms with particular expansions in Alveolata and Chlorophyta, and are not unique to Apicomplexa as previously believed. All RAP proteins identified can be decomposed into two parts. In the N-terminal region, the presence of variable helical repeats seems to participate in the specific targeting of diverse RNAs, while the RAP domain is mostly identified in the C-terminal region and is highly conserved across the different phylogenetic groups studied. Several conserved residues defining the signature motif could be crucial to ensure the function(s) of the RAP proteins. Modelling of RAP domains in apicomplexan parasites confirmed an ⍺/ß structure of a restriction endonuclease-like fold. The phylogenetic trees generated from multiple alignment of RAP domains and full-length proteins from various distantly related eukaryotes indicated a complex evolutionary history of this family. We further discuss these results to assess the potential function of this protein family in apicomplexan parasites.

The interplay of SARS-CoV-2 evolution and constraints imposed by the structure and functionality of its proteins.

Fast evolution of the SARS-CoV-2 virus provides us with unique information about the patterns of genetic changes in a single pathogen in the timescale of months. This data is used extensively to track the phylodynamic of the pandemic's spread and its split into distinct clades. Here we show that the patterns of SARS-CoV-2 virus mutations along its genome are closely correlated with the structural features of the coded proteins. We show that the foldability of proteins' 3D structures and conservation of their functions are the universal factors driving evolutionary selection in protein-coding genes. Insights from the analysis of mutation distribution in the context of the SARS-CoV-2 proteins' structures and functions have practical implications including evaluating potential antigen epitopes or selection of primers for PCR-based COVID-19 tests.

The crystal structure of nsp10-nsp16 heterodimer from SARS-CoV-2 in complex with S-adenosylmethionine.

SARS-CoV-2 is a member of the coronaviridae family and is the etiological agent of the respiratory Coronavirus Disease 2019. The virus has spread rapidly around the world resulting in over two million cases and nearly 150,000 deaths as of April 17, 2020. Since no treatments or vaccines are available to treat COVID-19 and SARS-CoV-2, respiratory complications derived from the infections have overwhelmed healthcare systems around the world. This virus is related to SARS-CoV-1, the virus that caused the 2002-2004 outbreak of Severe Acute Respiratory Syndrome. In January 2020, the Center for Structural Genomics of Infectious Diseases implemented a structural genomics pipeline to solve the structures of proteins essential for coronavirus replication-transcription. Here we show the first structure of the SARS-CoV-2 nsp10-nsp16 2'-O-methyltransferase complex with S-adenosylmethionine at a resolution of 1.80 Å. This heterodimer complex is essential for capping viral mRNA transcripts for efficient translation and to evade immune surveillance.

FATCAT 2.0: towards a better understanding of the structural diversity of proteins.

FATCAT 2.0 server (http://fatcat.godziklab.org/), provides access to a flexible protein structure alignment algorithm developed in our group. In such an alignment, rotations and translations between elements in the structure are allowed to minimize the overall root mean square deviation (RMSD) between the compared structures. This allows to effectively compare protein structures even if they underwent structural rearrangements in different functional forms, different crystallization conditions or as a result of mutations. The major update for the server introduces a new graphical interface, much faster database searches and several new options for visualization of the structural differences between proteins.

Coronavirus3D: 3D structural visualization of COVID-19 genomic divergence.

MOTIVATION: As the COVID-19 pandemic is spreading around the world, the SARS-CoV-2 virus is evolving with mutations that potentially change and fine-tune functions of the proteins coded in its genome. RESULTS: Coronavirus3D website integrates data on the SARS-CoV-2 virus mutations with information about 3D structures of its proteins, allowing users to visually analyze the mutations in their 3D context. AVAILABILITY AND IMPLEMENTATION: Coronavirus3D server is freely available at https://coronavirus3d.org.