Integration of individualized and population-level molecular epidemiology data to model COVID-19 outcomes.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with enhanced transmissibility and immune escape have emerged periodically throughout the coronavirus disease 2019 (COVID-19) pandemic, but the impact of these variants on disease severity has remained unclear. In this single-center, retrospective cohort study, we examined the association between SARS-CoV-2 clade and patient outcome over a two-year period in Chicago, Illinois. Between March 2020 and March 2022, 14,252 residual diagnostic specimens were collected from SARS-CoV-2-positive inpatients and outpatients alongside linked clinical and demographic metadata, of which 2,114 were processed for viral whole-genome sequencing. When controlling for patient demographics and vaccination status, several viral clades were associated with risk for hospitalization, but this association was negated by the inclusion of population-level confounders, including case count, sampling bias, and shifting standards of care. These data highlight the importance of integrating non-virological factors into disease severity and outcome models for the accurate assessment of patient risk.

Resurgence of SARS-CoV-2 Delta after Omicron variant superinfection in an immunocompromised pediatric patient.

BACKGROUND: Persistent SARS-CoV-2 infection in immunocompromised hosts is thought to contribute to viral evolution by facilitating long-term natural selection and viral recombination in cases of viral co-infection or superinfection. However, there are limited data on the longitudinal intra-host population dynamics of SARS-CoV-2 co-infection/superinfection, especially in pediatric populations. Here, we report a case of Delta-Omicron superinfection in a hospitalized, immunocompromised pediatric patient. METHODS: We conducted Illumina whole genome sequencing (WGS) for longitudinal specimens to investigate intra-host dynamics of SARS-CoV-2 strains. Topoisomerase PCR cloning of Spike open-reading frame and Sanger sequencing of samples was performed for four specimens to validate the findings. Analysis of publicly available SARS-CoV-2 sequence data was performed to investigate the co-circulation and persistence of SARS-CoV-2 variants. RESULTS: Results of WGS indicate the patient was initially infected with the SARS-CoV-2 Delta variant before developing a SARS-CoV-2 Omicron variant superinfection, which became predominant. Shortly thereafter, viral loads decreased below the level of detection before resurgence of the original Delta variant with no residual trace of Omicron. After 54 days of persistent infection, the patient tested negative for SARS-CoV-2 but ultimately succumbed to a COVID-19-related death. Despite protracted treatment with remdesivir, no antiviral resistance mutations emerged. These results indicate a unique case of persistent SARS-CoV-2 infection with the Delta variant interposed by a transient superinfection with the Omicron variant. Analysis of publicly available sequence data suggests the persistence and ongoing evolution of Delta subvariants despite the global predominance of Omicron, potentially indicative of continued transmission in an unknown population or niche. CONCLUSION: A better understanding of SARS-CoV-2 intra-host population dynamics, persistence, and evolution during co-infections and/or superinfections will be required to continue optimizing patient care and to better predict the emergence of new variants of concern.

Genome-wide screens reveal shared and strain-specific genes that facilitate enteric colonization by Klebsiella pneumoniae.

Gastrointestinal (GI) colonization by Klebsiella pneumoniae is a risk factor for subsequent infection as well as transmission to other patients. Additionally, colonization is achieved by many strain types that exhibit high diversity in genetic content. Thus, we aimed to study strain-specific requirements for K. pneumoniae GI colonization by applying transposon insertion sequencing to three classical clinical strains: a carbapenem-resistant strain, an extended-spectrum beta-lactamase-producing strain, and a non-epidemic antibiotic-susceptible strain. The transposon insertion libraries were screened in a murine model of GI colonization. At 3 days post-inoculation, 27 genes were required by all three strains for colonization. Isogenic deletion mutants for three genes/operons (acrA, carAB, and tatABCD) confirmed colonization defects in each of the three strains. Additionally, deletion of acrA reduced bile tolerance in vitro, while complementation restored both bile tolerance in vitro and colonization ability in vivo. Transposon insertion sequencing suggested that some genes were more important for the colonization of one strain than the others. For example, deletion of the sucrose porin-encoding gene scrY resulted in a colonization defect in the carbapenemase-producing strain but not in the extended-spectrum beta-lactamase producer or the antibiotic-susceptible strain. These findings demonstrate that classical K. pneumoniae strains use both shared and strain-specific strategies to colonize the mouse GI tract. IMPORTANCE Klebsiella pneumoniae is a common cause of difficult-to-treat infections due to its propensity to express resistance to many antibiotics. For example, carbapenem-resistant K. pneumoniae has been named an urgent threat by the United States Centers for Disease Control and Prevention. Gastrointestinal colonization in patients with K. pneumoniae has been linked to subsequent infection, making it a key process to control in the prevention of multidrug-resistant infections. However, the bacterial factors which contribute to K. pneumoniae colonization are not well understood. Additionally, individual strains exhibit large amounts of genetic diversity, begging the question of whether some colonization factors are strain dependent. This study identifies the enteric colonization factors of three classical strains using transposon mutant screens to define a core colonization program for K. pneumoniae as well as detecting strain-to-strain differences in colonization strategies.

Genome-wide screens reveal shared and strain-specific genes that facilitate enteric colonization by Klebsiella pneumoniae.

Gastrointestinal (GI) colonization by Klebsiella pneumoniae is a risk factor for subsequent infection as well as transmission to other patients. Additionally, colonization is achieved by many strain types that exhibit high diversity in genetic content. Thus, we aimed to study strain-specific requirements for K. pneumoniae GI colonization by applying transposon insertion sequencing to three classical clinical strains: a carbapenem-resistant strain, an extended-spectrum beta-lactamase producing strain, and a non-epidemic antibiotic-susceptible strain. The transposon insertion libraries were screened in a murine model of GI colonization. At three days post-inoculation, 27 genes were required by all three strains for colonization. Isogenic deletion mutants for three genes/operons (acrA, carAB, tatABCD) confirmed colonization defects in each of the three strains. Additionally, deletion of acrA reduced bile tolerance in vitro, while complementation restored both bile tolerance in vitro and colonization ability in vivo. Transposon insertion sequencing suggested that some genes were more important for colonization of one strain than the others. For example, deletion of the sucrose porin-encoding gene scrY resulted in a colonization defect in the carbapenemase-producing strain but not in the extended-spectrum beta-lactamase producer or the antibiotic-susceptible strain. These findings demonstrate that classical K. pneumoniae strains use both shared and strain-specific strategies to colonize the mouse GI tract. IMPORTANCE: Klebsiella pneumoniae is a common cause of difficult-to-treat infections due to its propensity to express resistance to many antibiotics. For example, carbapenem-resistant K. pneumoniae (CR-Kp) has been named an urgent threat by the United States Centers for Disease Control and Prevention. Gastrointestinal colonization of patients with K. pneumoniae has been linked to subsequent infection, making it a key process to control in prevention of multidrug-resistant infections. However, the bacterial factors which contribute to K. pneumoniae colonization are not well understood. Additionally, individual strains exhibit large amounts of genetic diversity, begging the question of whether some colonization factors are strain-dependent. This study identifies the enteric colonization factors of 3 classical strains using transposon mutant screens to define a core colonization program for K. pneumoniae as well as detecting strain-to-strain differences in colonization strategies.

Quantity of SARS-CoV-2 RNA copies exhaled per minute during natural breathing over the course of COVID-19 infection.

SARS-CoV-2 is spread through exhaled breath of infected individuals. A fundamental question in understanding transmission of SARS-CoV-2 is how much virus an individual is exhaling into the environment while they breathe, over the course of their infection. Research on viral load dynamics during COVID-19 infection has focused on internal swab specimens, which provide a measure of viral loads inside the respiratory tract, but not on breath. Therefore, the dynamics of viral shedding on exhaled breath over the course of infection are poorly understood. Here, we collected exhaled breath specimens from COVID-19 patients and used RTq-PCR to show that numbers of exhaled SARS-CoV-2 RNA copies during COVID-19 infection do not decrease significantly until day 8 from symptom-onset. COVID-19-positive participants exhaled an average of 80 SARS-CoV-2 viral RNA copies per minute during the first 8 days of infection, with significant variability both between and within individuals, including spikes over 800 copies a minute in some patients. After day 8, there was a steep drop to levels nearing the limit of detection, persisting for up to 20 days. We further found that levels of exhaled viral RNA increased with self-rated symptom-severity, though individual variation was high. Levels of exhaled viral RNA did not differ across age, sex, time of day, vaccination status or viral variant. Our data provide a fine-grained, direct measure of the number of SARS-CoV-2 viral copies exhaled per minute during natural breathing-including 312 breath specimens collected multiple times daily over the course of infection-in order to fill an important gap in our understanding of the time course of exhaled viral loads in COVID-19.

Systems analysis reveals differential expression of endocervical genes in African women randomized to DMPA-IM, LNG implant or cu-IUD.

Although effective contraceptives are crucial for preventing unintended pregnancies, evidence suggests that their use may perturb the female genital tract (FGT). A comparative analysis of the effects of the most common contraceptives on the FGT have not been evaluated in a randomized clinical trial setting. Here, we evaluated the effect of three long-acting contraceptive methods: depot medroxyprogesterone acetate(DMPA-IM), levonorgestrel(LNG) implant, and a copper intrauterine device (Cu-IUD), on the endocervical host transcriptome in 188 women from the Evidence for Contraceptive Options and HIV Outcomes Trial (ECHO) trial. Cu-IUD usage showed the most extensive transcriptomic changes, and was associated with inflammatory and anti-viral host responses. DMPA-IM usage was enriched for pathways associated with T cell responses. LNG implant had the mildest effect on endocervical gene expression, and was associated with growth factor signaling. These data provide a mechanistic basis for the diverse influence that varying contraceptives have on the FGT.

Large library docking for novel SARS-CoV-2 main protease non-covalent and covalent inhibitors.

Antiviral therapeutics to treat SARS-CoV-2 are needed to diminish the morbidity of the ongoing COVID-19 pandemic. A well-precedented drug target is the main viral protease (MPro ), which is targeted by an approved drug and by several investigational drugs. Emerging viral resistance has made new inhibitor chemotypes more pressing. Adopting a structure-based approach, we docked 1.2 billion non-covalent lead-like molecules and a new library of 6.5 million electrophiles against the enzyme structure. From these, 29 non-covalent and 11 covalent inhibitors were identified in 37 series, the most potent having an IC50 of 29 and 20 µM, respectively. Several series were optimized, resulting in low micromolar inhibitors. Subsequent crystallography confirmed the docking predicted binding modes and may template further optimization. While the new chemotypes may aid further optimization of MPro inhibitors for SARS-CoV-2, the modest success rate also reveals weaknesses in our approach for challenging targets like MPro versus other targets where it has been more successful, and versus other structure-based techniques against MPro itself.

Genomics of Aminoglycoside Resistance in Pseudomonas aeruginosa Bloodstream Infections at a United States Academic Hospital.

Pseudomonas aeruginosa frequently becomes resistant to aminoglycosides by the acquisition of aminoglycoside modifying enzyme (AME) genes and the occurrence of mutations in the mexZ, fusA1, parRS, and armZ genes. We examined resistance to aminoglycosides in a collection of 227 P. aeruginosa bloodstream isolates collected over 2 decades from a single United States academic medical institution. Resistance rates of tobramycin and amikacin were relatively stable over this time, while the resistance rates of gentamicin were somewhat more variable. For comparison, we examined resistance rates to piperacillin-tazobactam, cefepime, meropenem, ciprofloxacin, and colistin. Resistance rates to the first four antibiotics were also stable, although uniformly higher for ciprofloxacin. Colistin resistance rates were initially quite low, rose substantially, and then began to decrease at the end of the study. Clinically relevant AME genes were identified in 14% of isolates, and mutations predicted to cause resistance were relatively common in the mexZ and armZ genes. In a regression analysis, resistance to gentamicin was associated with the presence of at least one gentamicin-active AME gene and significant mutations in mexZ, parS, and fusA1. Resistance to tobramycin was associated with the presence of at least one tobramycin-active AME gene. An extensively drug-resistant strain, PS1871, was examined further and found to contain five AME genes, most of which were within clusters of antibiotic resistance genes embedded in transposable elements. These findings demonstrate the relative contributions of aminoglycoside resistance determinants to P. aeruginosa susceptibilities at a United States medical center. IMPORTANCE Pseudomonas aeruginosa is frequently resistant to multiple antibiotics, including aminoglycosides. The rates of resistance to aminoglycosides in bloodstream isolates collected over 2 decades at a United States hospital remained constant, suggesting that antibiotic stewardship programs may be effective in countering an increase in resistance. Mutations in the mexZ, fusA1, parR, pasS, and armZ genes were more common than acquisition of genes encoding aminoglycoside modifying enzymes. The whole-genome sequence of an extensively drug resistant isolate indicates that resistance mechanisms can accumulate in a single strain. Together, these results suggest that aminoglycoside resistance in P. aeruginosa remains problematic and confirm known resistance mechanisms that can be targeted for the development of novel therapeutics.

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.

Lactic acid from vaginal microbiota enhances cervicovaginal epithelial barrier integrity by promoting tight junction protein expression.

BACKGROUND: Women with a cervicovaginal microbiota dominated by Lactobacillus spp. are at reduced risk of acquiring sexually transmitted infections including HIV, but the biological mechanisms involved remain poorly defined. Here, we performed metaproteomics on vaginal swab samples from young South African women (n = 113) and transcriptomics analysis of cervicovaginal epithelial cell cultures to examine the ability of lactic acid, a metabolite produced by cervicovaginal lactobacilli, to modulate genital epithelial barrier function. RESULTS: Compared to women with Lactobacillus-depleted microbiota, women dominated by vaginal lactobacilli exhibit higher abundance of bacterial lactate dehydrogenase, a key enzyme responsible for lactic acid production, which is independently associated with an increased abundance of epithelial barrier proteins. Physiological concentrations of lactic acid enhance epithelial cell culture barrier integrity and increase intercellular junctional molecule expression. CONCLUSIONS: These findings reveal a novel ability of vaginal lactic acid to enhance genital epithelial barrier integrity that may help prevent invasion by sexually transmitted pathogens. Video abstract.

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.

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.

Investigation into the bacterial diversity of sediment samples obtained from Berg River, Western Cape, South Africa.

This study used conventional culturing and 16S rRNA metagenomics analyses to assess the diversity of bacterial communities in sediment samples obtained from the Berg River, Western Cape, South Africa. Samples were collected from six points: a residential and recreational area, an industrial area, an informal residential settlement, a point next to a wastewater treatment plant (WWTP), a pumping station, and a residential and agricultural farming area along the river. High bacterial counts recorded on general selective and differential culture media signify substantial microbial contamination along the sampling sites. The most prevalent bacterial phyla detected (through metagenomics analyses) along the sampling sites were Proteobacteria (61%), Planctomycetes (9.5%), Firmicutes (7.8%), Bacteroidetes (5%), Acidobacteria (4.6%), and Actinobacteria (4.6%). Some members of the identified predominant bacterial phyla, genera, and classes are important public health bacteria that have been implicated in human diseases and outbreaks, while some others are metal or hydrocarbon tolerant, indicating possible significant environmental pollution. Notable human pathogenic genera such as Bacillus, Clostridium, Shigella, Legionella, Mycobacterium, and Pseudomonas were identified in varying percentages at five of the six sampling areas. Fecal contamination was particularly rife at all residential areas, with the informal housing area being the most notably polluted. Diverse functional pathways were predicted for identified bacteria, such as those associated with different chronic and infectious human diseases as well as those related to hydrocarbon and metal remediation. The point next to a WWTP contained vastly diverse groups of bacterial contaminants as well as the most abundant pathway identities and titles.

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.

Microbial function and genital inflammation in young South African women at high risk of HIV infection.

BACKGROUND: Female genital tract (FGT) inflammation is an important risk factor for HIV acquisition. The FGT microbiome is closely associated with inflammatory profile; however, the relative importance of microbial activities has not been established. Since proteins are key elements representing actual microbial functions, this study utilized metaproteomics to evaluate the relationship between FGT microbial function and inflammation in 113 young and adolescent South African women at high risk of HIV infection. Women were grouped as having low, medium, or high FGT inflammation by K-means clustering according to pro-inflammatory cytokine concentrations. RESULTS: A total of 3186 microbial and human proteins were identified in lateral vaginal wall swabs using liquid chromatography-tandem mass spectrometry, while 94 microbial taxa were included in the taxonomic analysis. Both metaproteomics and 16S rRNA gene sequencing analyses showed increased non-optimal bacteria and decreased lactobacilli in women with FGT inflammatory profiles. However, differences in the predicted relative abundance of most bacteria were observed between 16S rRNA gene sequencing and metaproteomics analyses. Bacterial protein functional annotations (gene ontology) predicted inflammatory cytokine profiles more accurately than bacterial relative abundance determined by 16S rRNA gene sequence analysis, as well as functional predictions based on 16S rRNA gene sequence data (p < 0.0001). The majority of microbial biological processes were underrepresented in women with high inflammation compared to those with low inflammation, including a Lactobacillus-associated signature of reduced cell wall organization and peptidoglycan biosynthesis. This signature remained associated with high FGT inflammation in a subset of 74 women 9 weeks later, was upheld after adjusting for Lactobacillus relative abundance, and was associated with in vitro inflammatory cytokine responses to Lactobacillus isolates from the same women. Reduced cell wall organization and peptidoglycan biosynthesis were also associated with high FGT inflammation in an independent sample of ten women. CONCLUSIONS: Both the presence of specific microbial taxa in the FGT and their properties and activities are critical determinants of FGT inflammation. Our findings support those of previous studies suggesting that peptidoglycan is directly immunosuppressive, and identify a possible avenue for biotherapeutic development to reduce inflammation in the FGT. To facilitate further investigations of microbial activities, we have developed the FGT-DB application that is available at http://fgtdb.org/ . Video Abstract.

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.

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.