Molecular Evolutionary Analyses of the Spike Protein Gene and Spike Protein in the SARS-CoV-2 Omicron Subvariants.

To better understand the evolution of the SARS-CoV-2 Omicron subvariants, we performed molecular evolutionary analyses of the spike (S) protein gene/S protein using advanced bioinformatics technologies. First, time-scaled phylogenetic analysis estimated that a common ancestor of the Wuhan, Alpha, Beta, Delta variants, and Omicron variants/subvariants diverged in May 2020. After that, a common ancestor of the Omicron variant generated various Omicron subvariants over one year. Furthermore, a chimeric virus between the BM.1.1.1 and BJ.1 subvariants, known as XBB, diverged in July 2021, leading to the emergence of the prevalent subvariants XBB.1.5 and XBB.1.16. Next, similarity plot (SimPlot) data estimated that the recombination point (breakpoint) corresponded to nucleotide position 1373. As a result, XBB.1.5 subvariants had the 5' nucleotide side from the breakpoint as a strain with a BJ.1 sequence and the 3' nucleotide side as a strain with a BM.1.1.1 sequence. Genome network data showed that Omicron subvariants were genetically linked with the common ancestors of the Wuhan and Delta variants, resulting in many amino acid mutations. Selective pressure analysis estimated that the prevalent subvariants, XBB.1.5 and XBB.1.16, had specific amino acid mutations, such as V445P, G446S, N460K, and F486P, located in the RBD when compared with the BA.4 and BA.5 subvariants. Moreover, some representative immunogenicity-associated amino acid mutations, including L452R, F486V, R493Q, and V490S, were also found in these subvariants. These substitutions were involved in the conformational epitopes, implying that these mutations affect immunogenicity and vaccine evasion. Furthermore, these mutations were identified as positive selection sites. These results suggest that the S gene/S protein Omicron subvariants rapidly evolved, and mutations observed in the conformational epitopes may reduce the effectiveness of the current vaccine, including bivalent vaccines such as mRNA vaccines containing the BA.4/BA.5 subvariants.

Molecular evolutionary analyses of the fusion protein gene in human respirovirus 1.

Few evolutionary studies of the human respiratory virus (HRV) have been conducted, but most of them have focused on HRV3. In this study, the full-length fusion (F) genes in HRV1 strains collected from various countries were subjected to time-scaled phylogenetic, genome population size, and selective pressure analyses. Antigenicity analysis was performed on the F protein. The time-scaled phylogenetic tree using the Bayesian Markov Chain Monte Carlo method estimated that the common ancestor of the HRV1 F gene diverged in 1957 and eventually formed three lineages. Phylodynamic analyses showed that the genome population size of the F gene has doubled over approximately 80 years. Phylogenetic distances between the strains were short (< 0.02). No positive selection sites were detected for the F protein, whereas many negative selection sites were identified. Almost all conformational epitopes of the F protein, except one in each monomer, did not correspond to the neutralising antibody (NT-Ab) binding sites. These results suggest that the HRV1 F gene has constantly evolved over many years, infecting humans, while the gene may be relatively conserved. Mismatches between computationally predicted epitopes and NT-Ab binding sites may be partially responsible for HRV1 reinfection and other viruses such as HRV3 and respiratory syncytial virus.

Molecular Evolutionary Analyses of the Pseudomonas-Derived Cephalosporinase Gene.

Despite the increasing evidence of the clinical impact of Pseudomonas-derived cephalosporinase (PDC) sequence polymorphisms, the molecular evolution of its encoding gene, blaPDC, remains elusive. To elucidate this, we performed a comprehensive evolutionary analysis of blaPDC. A Bayesian Markov Chain Monte Carlo phylogenetic tree revealed that a common ancestor of blaPDC diverged approximately 4660 years ago, leading to the formation of eight clonal variants (clusters A-H). The phylogenetic distances within clusters A to G were short, whereas those within cluster H were relatively long. Two positive selection sites and many negative selection sites were estimated. Two PDC active sites overlapped with negative selection sites. In docking simulation models based on samples selected from clusters A and H, piperacillin was bound to the serine and the threonine residues of the PDC active sites, with the same binding mode for both models. These results suggest that, in P. aeruginosa, blaPDC is highly conserved, and PDC exhibits similar antibiotic resistance functionality regardless of its genotype.

Detailed Molecular Interactions between Respiratory Syncytial Virus Fusion Protein and the TLR4/MD-2 Complex In Silico.

Molecular interactions between respiratory syncytial virus (RSV) fusion protein (F protein) and the cellular receptor Toll-like receptor 4 (TLR4) and myeloid differentiation factor-2 (MD-2) protein complex are unknown. Thus, to reveal the detailed molecular interactions between them, in silico analyses were performed using various bioinformatics techniques. The present simulation data showed that the neutralizing antibody (NT-Ab) binding sites in both prefusion and postfusion proteins at sites II and IV were involved in the interactions between them and the TLR4 molecule. Moreover, the binding affinity between postfusion proteins and the TLR4/MD-2 complex was higher than that between prefusion proteins and the TLR4/MD-2 complex. This increased binding affinity due to conformational changes in the F protein may be able to form syncytium in RSV-infected cells. These results may contribute to better understand the infectivity and pathogenicity (syncytium formation) of RSV.

Molecular Evolution of the Pseudomonas aeruginosa DNA Gyrase gyrA Gene.

DNA gyrase plays important roles in genome replication in various bacteria, including Pseudomonasaeruginosa. The gyrA gene encodes the gyrase subunit A protein (GyrA). Mutations in GyrA are associated with resistance to quinolone-based antibiotics. We performed a detailed molecular evolutionary analyses of the gyrA gene and associated resistance to the quinolone drug, ciprofloxacin, using bioinformatics techniques. We produced an evolutionary phylogenetic tree using the Bayesian Markov Chain Monte Carlo (MCMC) method. This tree indicated that a common ancestor of the gene was present over 760 years ago, and the offspring formed multiple clusters. Quinolone drug-resistance-associated amino-acid substitutions in GyrA, including T83I and D87N, emerged after the drug was used clinically. These substitutions appeared to be positive selection sites. The molecular affinity between ciprofloxacin and the GyrA protein containing T83I and/or D87N decreased significantly compared to that between the drug and GyrA protein, with no substitutions. The rate of evolution of the gene before quinolone drugs were first used in the clinic, in 1962, was significantly lower than that after the drug was used. These results suggest that the gyrA gene evolved to permit the bacterium to overcome quinolone treatment.