Aa Substitutions Research Articles

Diversity, functional classification and genotyping of SHV β-lactamases in Klebsiella pneumoniae.

Interpreting the phenotypes of bla SHV alleles in Klebsiella pneumoniae genomes is complex. Whilst all strains are expected to carry a chromosomal copy conferring resistance to ampicillin, they may also carry mutations in chromosomal bla SHV alleles or additional plasmid-borne bla SHV alleles that have extended-spectrum β-lactamase (ESBL) activity and/or β-lactamase inhibitor (BLI) resistance activity. In addition, the role of individual mutations/a changes is not completely documented or understood. This has led to confusion in the literature and in antimicrobial resistance (AMR) gene databases [e.g. the National Center for Biotechnology Information (NCBI) Reference Gene Catalog and the β-lactamase database (BLDB)] over the specific functionality of individual sulfhydryl variable (SHV) protein variants. Therefore, the identification of ESBL-producing strains from K. pneumoniae genome data is complicated. Here, we reviewed the experimental evidence for the expansion of SHV enzyme function associated with specific aa substitutions. We then systematically assigned SHV alleles to functional classes (WT, ESBL and BLI resistant) based on the presence of these mutations. This resulted in the re-classification of 37 SHV alleles compared with the current assignments in the NCBI's Reference Gene Catalog and/or BLDB (21 to WT, 12 to ESBL and 4 to BLI resistant). Phylogenetic and comparative genomic analyses support that (i) SHV-1 (encoded by bla SHV-1) is the ancestral chromosomal variant, (ii) ESBL- and BLI-resistant variants have evolved multiple times through parallel substitution mutations, (iii) ESBL variants are mostly mobilized to plasmids and (iv) BLI-resistant variants mostly result from mutations in chromosomal bla SHV. We used matched genome-phenotype data from the KlebNET-GSP AMR Genotype-Phenotype Group to identify 3999 K. pneumoniae isolates carrying one or more bla SHV alleles but no other acquired β-lactamases to assess genotype-phenotype relationships for bla SHV. This collection includes human, animal and environmental isolates collected between 2001 and 2021 from 24 countries. Our analysis supports that mutations at Ambler sites 238 and 179 confer ESBL activity, whilst most omega-loop substitutions do not. Our data also provide support for the WT assignment of 67 protein variants, including 8 that were noted in public databases as ESBL. These eight variants were reclassified as WT because they lack ESBL-associated mutations, and our phenotype data support susceptibility to third-generation cephalosporins (SHV-27, SHV-38, SHV-40, SHV-41, SHV-42, SHV-65, SHV-164 and SHV-187). The approach and results outlined here have been implemented in Kleborate v2.4.1 (a software tool for genotyping K. pneumoniae), whereby known and novel bla SHV alleles are classified based on causative mutations. Kleborate v2.4.1 was updated to include ten novel protein variants from the KlebNET-GSP dataset and all alleles in public databases as of November 2023. This study demonstrates the power of sharing AMR phenotypes alongside genome data to improve the understanding of resistance mechanisms.

Epidemiology and molecular evolution of GI.1 sapovirus in the recent era.

Sapovirus (SaV) infection is increasing worldwide. Herein, we provided evidence of a significant increase in SaV infection in Japan during 2010-2022, primarily due to the considerable (p = 0.0003) rise of the GI.1 genotype. Furthermore, we found that all major and minor SaV outbreaks in Japan, including the largest SaV outbreak in 2021-2022, were caused by the GI.1 genotype. Therefore, to get insight into the underlying molecular mechanism behind this rising trend of the SaV GI.1 type, we selected 15 SaV GI.1 outbreak strains for complete genome analysis through next-generation sequencing. Phylogenetically, our strains remained clustered in different branches in lineages I and II among the GI.1 genotype. We showed all amino acid (aa) substitutions in different open reading frames (ORFs) in these strains. Importantly, we have demonstrated that the strains involved in the largest SaV outbreak in Japan in 2021-2022 belonged to lineage II and possessed the third ORF. We have identified some unique aa mutations in these major outbreak strains in the NS1 and NS6-NS7 regions that are thought to be associated with viral pathogenicity, cell tropism, and epidemiological competence. Thus, in addition to enriching the database of SaV's complete sequences, this study provides insights into its important mutations.

Polymorphism Identification in the Coding Sequences (ORFs) of the Porcine Pregnancy-Associated Glycoprotein 2-like Gene Subfamily in Pigs.

Pregnancy-associated glycoproteins (PAGs) are a polygenic family with many scattered genes and pseudogenes resulting from the duplication or fusion of a pseudogene with expression beginning in the trophoblast during the peri-implantation period and continuing in the trophectoderm. In this study, single-nucleotide polymorphism (SNP) and insertion/deletion (InDels) in the open reading frame (nine exons) of crossbreed pigs are reported for the first time. Novel SNPs/InDels were researched using genomic DNA templates isolated from the leukocytes of crossbreed pigs (N = 25), which were amplified, gel-out-purified, and sequenced. Sixteen SNPs and one InDel (g.6961_6966 Ins TGCCAA) were identified in the crossbreed pigs. In silico analysis revealed that among 16 SNPs, only 10 SNPs cause amino acid (aa) substitutions, and InDel codes asparagine (N298) and alanine (A299). The results provide a novel broad-based database (main pattern) that will be critical for future research into the possible correlations between the SNP genotypes of the pPAG2-L subfamily in pigs of various breeds whose reproductive traits are known.

Molecular and phylogenetic analysis of transmissible gastroenteritis virus strain VET-16, isolated from piglets in Vietnam.

Porcine transmissible gastroenteritis virus (TGEV) is a major pathogen that causes viral enteritis and severe diarrhea in newborn piglets. TGEV strains have been isolated in the USA, Europe, and China, and their molecular characteristics are well known. However, there have been few reports of molecular analysis of TGEV strains isolated in Southeast Asia. In 2016, we isolated TGEV strain VET-16 from fecal samples collected from piglets in Vietnam and determined its complete genome sequence by Sanger sequencing. We found that, while the full genome of the VET-16 strain was 92.4-99.9% identical to those of other TGEV strains, the ORF3 gene showed very little sequence similarity. Phylogenetic analysis suggested that the VET-16 strain belongs to the Purdue subgroup. Comparison of the predicted amino acid (aa) sequence of the spike protein of strain VET-16 with those of other TGEV strains revealed three aa substitutions (V378L, S379T, and D380N) and a 3-aa insertion (F383_F387insWEK) in antigenic site D of the VET-16 strain. Also, a single aa deletion (∆F1413) was found in the transmembrane domain of the spike gene of VET-16. Like the ORF3 gene from the TGEV Miller M60 vaccine strain, the VET-16 strain has a large deletion (∆725 nt) in the ORF3 gene. Previous studies have suggested that these mutations in the spike and ORF3 genes might be associated with a reduction in pathogenicity. The data from this study will facilitate further genetic analysis and research into the evolution of TGEV in pigs in Vietnam.

Genetic Sequence and Pathogenicity of Infectious Bursal Disease Virus in Chickens in Egypt During 2017-2021.

The continued circulation of infectious bursal disease virus (IBDV) in Egypt, despite the use of various vaccines, is a serious problem that requires continuous detection of IBDV. In the current study, real-time reverse transcriptase polymerase chain reaction testing of 100 diseased chicken flocks during 2017-2021 revealed the presence of very virulent IBDV (vvIBDV) in 67% of the flocks, non-vvIBDV in 11%, and a mixture of both vvIBDV and non-vvIBDV in 4%. Twenty-nine IBDV isolates were submitted for partial sequencing of the viral protein 2 hypervariable region (VP2-HVR), and 27 isolates were confirmed to be genogroup A3 (vvIBDV) with 96.3%-98.5% similarity to the global A3 (vvIBDV) and 88.9%-97% similarity to genogroup A1 vaccine strains. The remaining two isolates were non-vvIBDV and showed 91.1% and 100% identity with classical genogroup A1 strains, respectively. Furthermore, the sequence and phylogenetic analysis of VP1 (amino acids 33-254) of two selected isolates of A3, 5/2017 and 98/2021, clustered them as B2, vvIBDV-like, strains with high similarity (99.5%) to four Egyptian, 99% to Chinese and European, and 97.7% to Chinese and Polish vvIBDV isolates. Experimental infection of commercial broiler chickens with two vvIBDV-A3B2 isolates (5/2017 and 98/2021) showed no mortality despite typical tissue lesions, clear histopathological changes, and strong ELISA antibody response. Isolate 98/2021 was more pathogenic, as confirmed by histopathology, whereas isolate 5/2017 induced a stronger serological response. In conclusion, vvIBDV (A3B2) strains with two amino acid (aa) substitutions in VP1 as V141I and V234I as well as VP2 as Y220F and G254S are still circulating in Egypt.

Functional Characterization of RNA Silencing Suppressor Encoded by Cotton Leafroll Dwarf Virus

The P0 proteins encoded by Cotton leafroll dwarf virus (CLRDV) serve as viral suppressors of RNA silencing (VSR). CLRDV P0 proteins share a VSR-associated F-box-like motif. This study investigated the role of specific amino acid (aa) residues within the F-box-like motif, elucidating their implications for VSR potency, pathogenicity, and intracellular localization. Analysis of various single aa substitution mutants within the F-box-like motif to understand their effect on P0 proteins’ VSR potency revealed that specific mutations significantly affect P0’s ability to suppress RNA silencing, emphasizing the crucial role of the F-box-like motif. Subcellular localization examinations indicated that the P0 proteins associate with the nucleus and endoplasmic reticulum, which may be linked to their VSR function. Also, the induction of hypersensitive response (HR)-like symptoms and the accumulation of reactive oxygen species by P0 proteins was observed, suggesting their role as pathogenicity factors. These results indicated that specific mutations significantly affect the functions of multifaceted P0 proteins, highlighting the F-box-like motif’s crucial role. Our study highlights the importance of further molecular investigations to elucidate how different CLRDV strains manifest various disease symptoms and severity. This is crucial considering the global economic importance of cotton and the potential emergence of more threatening CLRDV isolates.

Abstract IA011: Polymerase-based bypass of atypical UV photoproducts

Abstract UV light is a ubiquitous environmental mutagen that is an etiological cause of skin cancers like melanoma. While the genomes of human cutaneous melanomas have high mutation burdens primarily consisting of C to T substitutions in TC and CC dinucleotides (consistent with cyclobutane pyrimidine dimer formation), many of the driver mutations contributing to disease progression involve other types of substitutions at non-dipyrimidine sequences. Our whole genome sequencing of UVB or UVC irradiated yeast has revealed non-canonical UV induced substitution types. These mutations included TA to TT and AC to TT substitutions that commonly occur as activating BRAF mutations in human melanomas. UV treatment and sequencing of deletion mutants of DNA polymerase eta, the protein mutated in human XPV patients, indicated that this polymerase does not impact TA to TT mutagenesis, but contributes to the formation of AC to TT substitutions. Surprisingly, loss of pol eta revealed additional CA to AA substitutions, highlighting another UV-induced mutation class possibly created by a non-canonical photoproduct. All novel mutation classes in pol eta deficient yeast displayed transcriptional asymmetry, consistent with the repair of the cognate lesion by transcription-coupled nucleotide excision repair. However, UV-induced mutations lacked replication asymmetry in both wild type and rad30Δ yeast, suggesting TLS occurs nearly equally on leading and lagging strands. Both TA to TT and AC to TT substitutions were dependent on DNA polymerase zeta, and altered by a mutant of DNA polymerase epsilon, suggesting the involvement of both TLS and replicative polymerases in UV-induced mutagenesis. These results emphasize the importance of non-canonical UV lesions in driving melanomagenesis, provide mechanisms of how these lesions result in mutations, and suggest the variable accuracy of TLS polymerases across from atypical UV photoproducts may alter the trajectory of skin cancer evolution in XPV patients. Citation Format: Steven A. Roberts. Polymerase-based bypass of atypical UV photoproducts [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: DNA Damage Repair: From Basic Science to Future Clinical Application; 2024 Jan 9-11; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2024;84(1 Suppl):Abstract nr IA011.

Emergence of SARS-CoV-2 Variants Are Induced by Coinfections With Dengue.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that emerged in late 2019 has accumulated a series of point mutations and evolved into several variants of concern (VOCs), some of which are more transmissible and potentially more severe than the original strain. The most notable VOCs are Alpha, Beta, Gamma, Delta, and Omicron, which have spread to various parts of the world. This study conducted surveillance in Jashore, Bangladesh to identify the prevalence of SARS-CoV-2 coinfected with dengue virus and their genomic effect on the emergence of VOCs. A hospital-based COVID-19 surveillance from June to August, 2021 identified 9 453 positive patients in the surveillance area. The study enrolled 572 randomly selected COVID-19-positive patients, of which 11 (2%) had dengue viral coinfection. Whole genome sequences of SARS-CoV-2 were analyzed and compared between coinfection positive and negative group. In addition, we extracted 185 genome sequences from GISAID to investigate the cross-correlation function between SARS-CoV-2 mutations and VOC; multiple ARIMAX(p,d,q) models were developed to estimate the average number of amino acid (aa) substitution among different SARS-CoV-2 VOCs. The results of the study showed that the coinfection group had an average of 30.6 (±1.7) aa substitutions in SARS-CoV-2, whereas the dengue-negative COVID-19 group had that average of 25.6 (±1.8; P < .01). The coinfection group showed a significant difference of aa substitutions in open reading frame (ORF) and N-protein when compared to dengue-negative group (P = .03). Our ARIMAX models estimated that the emergence of SARS-CoV-2 variants Delta required additional 9 to 12 aa substitutions than Alpha, Beta, or Gamma variant. The emergence of Omicron accumulated additional 19 (95% confidence interval [CI]: 15.74, 21.95) aa substitution than Delta. Increased number of point mutations in SARS-CoV-2 genome identified from coinfected cases could be due to the compromised immune function of host and induced adaptability of pathogens during coinfections. As a result, new variants might be emerged when series of coinfection events occur during concurrent two epidemics.

Whole genome sequencing of outbreak strains from 2017 to 2018 reveals an endemic clade of dengue 1 virus in Cameroon

ABSTRACT Dengue fever is expanding as a global public health threat including countries within Africa. For the past few decades, Cameroon has experienced sporadic cases of arboviral infections including dengue fever. Here, we conducted genomic analyses to investigate the origin and phylogenetic profile of Cameroon DENV-1 outbreak strains and predict the impact of emerging therapeutics on these strains. Bayesian and maximum-likelihood phylogenetic inference approaches were employed in virus evolutionary analyses. An in silico analysis was performed to assess the divergence in immunotherapeutic and vaccine targets in the new genomes. Six complete DENV-1 genomes were generated from 50 samples that met a clinical definition for DENV infection. Phylogenetic analyses revealed that the strains from the current study belong to a sub-lineage of DENV-1 genotype V and form a monophyletic taxon with a 2012 strain from Gabon. The most recent common ancestor (TMRCA) of the Cameroon and Gabon strains was estimated to have existed around 2008. Comparing our sequences to the vaccine strains, 19 and 15 amino acid (aa) substitutions were observed in the immuno-protective prM-E protein segments of the Dengvaxia® and TetraVax-DV-TV003 vaccines, respectively. Epitope mapping revealed mismatches in aa residues at positions E155 and E161 located in the epitope of the human anti-DENV-1 monoclonal antibody HMAb 1F4. The new DENV strains constitute a conserved genomic pool of viruses endemic to the Central African region that needs prospective monitoring to track local viral evolution. Further work is needed to ascertain the performance of emerging therapeutics in DENV strains from the African region.

The stem region of group A transferase is crucial for its specificity, and its alteration promotes heterologous Forssman synthase activity

Some stem region mutants of human blood group A transferase (hAT) possess Forssman synthase (FS) activity, but very little is known about the mechanisms responsible for this enzymatic crosstalk. We performed confocal microscopy and image analysis to determine whether different intra-Golgi localization was accountable for this acquired activity. We also performed structural modeling and mutational and normal mode analyses. We introduced new mutations in the stem region and tested its FS and AT activities. No differences in subcellular localization were found between hAT and FS-positive mutants. AlphaFold models of hAT and mFS (mouse Forssman synthase) showed that the hAT stem region has a tether-like stem region, while in mFS, it encircles its catalytic domain. In silico analysis of FS-positive mutants indicated that stem region mutations induced structural changes, decreasing interatomic interactions and mobility of hAT that correlated with FS activity. Several additional mutations introduced in that region also bestowed FS activity without altering the AT activity: hAT 37–55 aa substitution by mFS 34–52, 37–55 aa deletion, and missense mutations: S46P, Q278Y, and Q286M. Stem region structure, mobility, and interactions are crucial for hAT specificity. Moreover, stem region mutations can lead to heterologous Forssman activity without changes in the catalytic machinery.

First Report of Resistance-Breaking Variants of Tomato Spotted Wilt Virus (TSWV) Infecting Tomatoes with the Sw-5 Tospovirus-Resistance Gene in North Carolina.

Widespread use of tomato cultivars with the Sw-5 resistance gene has led to the emergence of resistance-breaking (RB) strains of tomato spotted wilt virus across the globe. In June of 2022, tomato spotted wilt (TSW) symptoms were observed at two farms (A and B, within 15 miles of each other) in Rowan County, NC on several commercial TSW resistant tomato cultivars (all heterozygous for the Sw-5 gene). At farm A, ~10% of plants had symptomatic foliage with ~30% of fruit with symptoms, while at farm B, up to 50% of plants had symptomatic foliage with ~80% of fruit with symptoms. Visual symptoms included stunting, severe leaf curling and bronzing, necrotic lesions on leaves, petioles and stems, and concentric ring spots on fruit (Supplementary Fig. 1). TSWV ImmunoStrips (AgDia, Elkhart, IN) and reverse-transcription (RT)-PCR with NSm primers (di Rienzo et al 2018) confirmed the presence of TSWV in 12 symptomatic plants sampled across the two farms. Primers designed to detect Impatiens necrotic spot virus, groundnut ringspot virus, tomato chlorotic spot virus, tomato chlorosis virus, alfalfa mosaic virus, and tomato necrotic streak virus (ilarvirus, Badillo et al., 2016) failed to generate amplicons of the expected size from cDNA generated from these field samples. The amplicons from full-length NSm cDNA were sequenced from independent, single-leaflet isolates from the TSWV-positive plants (three from farm A, nine from farm B) with the expectation of finding an amino acid (aa) substitution associated with the Sw-5 RB phenotype identified previously in CA (C118Y, Batuman et al. 2017) or Spain (C118Y and T120N, Lopez et al. 2011). All three nucleotide sequences from farm A contained the NSm C118Y substitution reported in CA. All three sequences were 99% identical (including the C118Y mutation) to NCBI GenBank accession KU179600.1, a TSWV isolate collected from GA in 2014 with no cultivar information reported. The nine nucleotide sequences from farm B contained neither of the two previously reported aa substitutions associated with the RB phenotype. Instead, all contained a D122G substitution within a conserved region of the TSWV NSm protein reported to be involved in direct interaction with the Sw-5 protein (Zhu et al 2017). Likewise, Huang et al (2021) generated a D122A mutation in TSWV-NSm, resulting in failure to elicit a Sw-5 mediated hypersensitive response. Three NSm sequences retrieved from GenBank contained the D122G substitution (AY848921.1, HM015516.1, KU179582.1), however, this mutation was not implicated directly with RB phenotypes (Ciuffo et al., 2005; Lopez et al., 2011; Marshall, 2016). The RB phenotype was confirmed with the NC variants on 'Mountain Merit' (Sw-5) by two means of virus inoculation: mechanical, rub-inoculation with extracted sap from infected plants, and thrips transmission assays with lab colony-maintained, Frankliniella occidentalis, the western flower thrips. Symptomatic leaf tissue obtained from these inoculation assays tested positive for TSWV by DAS-ELISA (AgDia, Elkhart, IN) and RT-PCR with NSm primers, providing definitive evidence of the occurrence of RB-TSWV at both farms, and subsequent sequencing confirmed the C118Y and D122G substitutions. This report warrants further investigation of the putative origins, prevalence and epidemiological implications of RB-TSWV variants in NC tomato production, and the development of new sources of resistance to TSWV.

Search of a genomic sequence database for potential novel blood group antigens: Investigation into why some amino acid substitutions are not immunogenic.

Polypeptide blood group antigens are typically identified through investigation of the antibodies they induce. Human genome sequence databases are a new tool to identify AA substitutions that potentially create blood group antigens. The Erythrogene genomic sequence database was searched for missense mutations not known to be blood group antigens in the extracellular domains of selected RBC proteins in European populations. Any mutations found with prevalence of 1%-90% and not known to have induced antibodies in transfusion practice were analyzed using protein structural analysis and epitope prediction programs to determine why they apparently lack immunogenicity. Thirteen missense mutations not known to create blood group antigens were identified in the extracellular domains of Kell, BCAM, and RhD proteins, but not in RhCE, Urea Transporter 1 (Kidd), Atypical Chemokine Receptor 1 (Duffy), glycophorin A or glycophorin B. While 11 of the 13 mutations had low prevalence (<1%), a Kell Ser726Pro substitution and a BCAM Val196Ile substitution had predicted phenotype prevalences of 43.2% and 5.7%, respectively. Ser726Pro had multiple properties of a linear B-cell epitope, but possible suboptimal protein location for B-cell receptor binding and limited T-cell epitope possibilities. Val196Ile was not predicted to be in a linear B-cell epitope. Multiple potential new blood group antigens of low prevalence were identified. Whether they are antigenic remains to be determined. Two higher prevalence variants in Kell and BCAM are unlikely antigens, otherwise their antibodies presumably would already have been identified. Possible reasons for their poor immunogenicity were identified.

A novel allelic donkey β-lactoglobulin I protein isoform generated by a non-AUG translation initiation codon is associated with a nonsynonymous SNP

β-Lactoglobulin I (β-LG I) is one of the most important whey proteins in donkey milk. However, to our knowledge, there has been no study focusing on the full nucleotide sequences of this gene (BLG I). Current investigation of donkey BLG I gene is very limited with only 2 variants (A and B) characterized so far at the protein level. Recently, a new β-LG I variant, with a significantly higher mass (+1,915 Da) than known variants has been detected. In this study, we report the whole nucleotide sequence of the BLG I gene from 2 donkeys, whose milk samples are characterized by the β-LG I SDS-PAGE band with a normal electrophoretic mobility (18,514.25 Da, β-LG I B1 form) the first, and by the presence of a unique β-LG I band with a higher electrophoretic mobility (20,428.5 Da, β-LG I D form) the latter. A high genetic variability was found all over the 2 sequenced BLG I alleles. In particular, 16 polymorphic sites were found in introns, one in the 5' flanking region, 3 SNPs in the 5' untranslated region and one SNP in the coding region (g.458G > A) located at the 40th nucleotide of exon 2 and responsible for the AA substitutions p.Asp28 > Asn in the mature protein. Two SNPs (g.920-922CAC > TGT and g.1871G/A) were genotyped in 93 donkeys of 2 Italian breeds (60 Ragusana and 33 Amiatina, respectively) and the overall frequencies of g.920-922CAC and g.1871A were 0.3065 and 0.043, respectively. Only the rare allele g.1871A was observed to be associated with the slower migrating β-LG I. Considering this genetic diversity and those found in the database, it was possible to deduce at least 5 different alleles (BLG I A, B, B1, C, D) responsible for 4 potential β-LG I translations. Among these alleles, B1 and D are those characterized in the present research, with the D allele of real novel identification. Haplotype data analysis suggests an evolutionary pathway of donkey BLG I gene and a possible phylogenetic map is proposed. Analyses of mRNA secondary structure showed relevant changes in the structures, as consequence of the g.1871G > A polymorphism, that might be responsible for the recognition of an alternative initiation site providing an additional signal peptide. The extension of 19 AA sequence to the mature protein, corresponding to the canonical signal peptide with an additional alanine residue, is sufficient to provide the observed molecular weight of the slower migrating β-LG I encoded by the BLG I D allele.

Sewage-based surveillance shows presence of Klebsiella pneumoniae resistant against last resort antibiotics in the population in Bergen, Norway.

The aim of this study was to understand the prevalence of antibiotic resistance in Klebsiella pneumoniae present in the population in Bergen city, Norway using city-scale sewage-based surveillance, as well as the potential spread of K. pneumoniae into the marine environment through treated sewage. From a total of 30 sewage samples collected from five different sewage treatment plants (STPs), 563 presumptive K. pneumoniae isolates were obtained on Simmons Citrate Agar with myo-Inositol (SCAI) plates, and 44 presumptive K. pneumoniae isolates on SCAI plates with cefotaxime. Colistin resistance was observed in 35 isolates, while cefotaxime resistance and tigecycline resistance was observed in only five isolates each, out of 563 presumptive K. pneumoniae isolates. All 44 isolates obtained on cefotaxime-containing plates were multidrug-resistant, with 25% (n=11) showing resistance against tigecycline. Clinically important acquired antibiotic resistance genes (ARGs), like blaCTX-M-14, blaCTX-M-15, qnrS1, aac(3)-IIe, tet(A), and sul1, were detected in several sequenced Klebsiella spp. isolates (n=53). All sequenced colistin-resistant isolates (n=13) had a mutation in the mgrB gene with nucleotide substitution at position C88T creating a premature stop codon. All sequenced tigecycline-resistant isolates (n=4) harbored a Tet(A) variant with 22 amino acid (aa) substitutions compared to the reference protein. The sequenced K. pneumoniae isolates (n=44) belonged to 22 different sequence types (STs) with ST730 (29.5%) as most prevalent, followed by pathogenic ST307 (11.4%). Virulence factors, including aerobactin (iutA), enterobactin (entABCDEFS and fepABCDG), salmochelin (iro), and yersiniabactin (ybt) were detected in several sequenced K. pneumoniae isolates, suggesting pathogenicity potential. Heavy metal resistance genes were common in sequenced K. pneumoniae isolates (n=44) with silver (silABCEFPRS) and copper (pcoABDRS) resistance genes present in 79.5% of the isolates. Sewage-based surveillance can be a useful tool for understanding antibiotic resistance in pathogens present within a population and to provide up-to date information on the current resistance situation. Our study presents a framework for population-based surveillance of resistance in K. pneumoniae.

Unveiling mutation effects on the structural dynamics of the main protease from SARS-CoV-2 with hybrid simulation methods

The main protease of SARS-CoV-2 (called Mpro or 3CLpro) is essential for processing polyproteins encoded by viral RNA. Several Mpro mutations were found in SARS-CoV-2 variants, which are related to higher transmissibility, pathogenicity, and resistance to neutralization antibodies. Macromolecules adopt several favored conformations in solution depending on their structure and shape, determining their dynamics and function. In this study, we used a hybrid simulation method to generate intermediate structures along the six lowest frequency normal modes and sample the conformational space and characterize the structural dynamics and global motions of WT SARS-CoV-2 Mpro and 48 mutations, including mutations found in P.1, B.1.1.7, B.1.351, B.1.525 and B.1.429+B.1.427 variants. We tried to contribute to the elucidation of the effects of mutation in the structural dynamics of SARS-CoV-2 Mpro. A machine learning analysis was performed following the investigation regarding the influence of the K90R, P99L, P108S, and N151D mutations on the dimeric interface assembling of the SARS-CoV-2 Mpro. The parameters allowed the selection of potential structurally stable dimers, which demonstrated that some single surface aa substitutions not located at the dimeric interface (K90R, P99L, P108S, and N151D) are able to induce significant quaternary changes. Furthermore, our results demonstrated, by a Quantum Mechanics method, the influence of SARS-CoV-2 Mpro mutations on the catalytic mechanism, confirming that only one of the chains of the WT and mutant SARS-CoV-2 Mpros are prone to cleave substrates. Finally, it was also possible to identify the aa residue F140 as an important factor related to the increasing enzymatic reactivity of a significant number of SARS-CoV-2 Mpro conformations generated by the normal modes-based simulations.

Prevalence of HIV-1 Natural Polymorphisms and Integrase-Resistance-Associated Mutations in African Children.

Integrase inhibitors (INIs) are a potent option for HIV treatment. Limited data exist on INI resistance in West Africa, particularly in children living with HIV/AIDS. We determined the prevalence of integrase gene polymorphisms and the frequency of naturally occurring amino acid (aa) substitutions at positions associated with INI resistance. Dried blood spot (DBS) samples were obtained from one hundred and seven (107) HIV-1-infected children aged less than 15 years old in two West African countries, Benin and Mali. All children were naïve to INI treatment, 56 were naïve to anti-retroviral therapy (ART), and 51 had received ART. Genetic sequencing of HIV integrase was successful in 75 samples. The aa changes at integrase positions associated with INI resistance were examined according to the Stanford HIV Genotypic Resistance database. The median ages were 2.6 and 10 years for ART-naïve and -treated children, respectively. The most common subtypes observed were CRF02_AG (74.7%) followed by CRF06_cpx (20%). No major INI-resistance mutations at positions 66, 92, 121, 143, 147, 148, 155, and 263 were detected. The most prevalent INI accessory resistance mutations were: L74I/M (14/75, 18.6%) followed by E157Q (8/75, 10.6%), G163E/N/T/Q (5/75, 6.6%), Q95A/H/P (2/75, 2.6%), and T97A (4/75, 5.3%). Other substitutions observed were M50I/L/P, H51E/P/S/Q, I72V, T112V, V201I, and T206S. Polymorphisms at positions which may influence the genetic barrier and/or drive the selection of specific INI-resistance pathways were detected. However, no transmitted drug resistance (TDR) to INI was detected among samples of INI-naïve patients. These findings support the use of this treatment class for children with HIV-1, particularly in West Africa.

Genomic and comparative clinico-pathological assessment of two Pakistani pigeon-derived newcastle disease virus sub-genotypes XXI.1.1 and XXI.1.2 isolated in 2017

Pigeon paramyxovirus type-1 (PPMV-1) is an antigenic-variant of Newcastle disease virus (NDV) which is associated with infection in Columbidae family. In this study, we isolated two pigeon-derived strains pi/Pak/Lhr/SA_1/17 (designed as SA_1) and pi/Pak/Lhr/SA_2/17 (designed as SA_2) from diseased pigeons collected in Punjab province in 2017. We performed the whole genome, phylogenetic analysis and comparative clinico-pathological evaluation of two viruses in pigeons. Phylogenetic analysis based on fusion (F) gene and complete genome sequences showed that SA_1 belonged to sub-genotype XXI.1.1 and SA_2 clustered in sub-genotype XXI.1.2. SA_1 and SA_2 viruses contributed to morbidity and mortality in pigeons. Remarkably, although the two viruses resulted in comparatively similar pattern of pathogenesis and replication ability in various tissues of infected pigeons, SA_2 could cause more severe histopathological lesions and had comparatively high replication ability in pigeons than SA_1. Moreover, pigeons infected with SA_2 had higher shedding efficiency than that of pigeons infected with SA_1. Moreover, several aa substitutions in the major functional domains of the F and HN proteins might be contributed to the pathogenic differences between the two isolates in pigeons. Overall, these findings provide us with important insight into the epidemiology and evolution of PPMV-1 in Pakistan and laid the foundation for the further elucidation of the mechanism underlying the pathogenic difference of PPMV-1 in pigeons.

A novel C-terminal modification method enhanced the yield of human papillomavirus L1 or chimeric L1-L2 virus-like particles in the baculovirus system.

Human papillomavirus (HPV) major capsid protein L1 virus-like particles (VLPs) produced in the baculovirus system showed excellent safety and immunogenicity, but the relatively high production cost stands as a substantial barrier to extensive commercialization, especially in producing multivalent vaccines. Here, a novel method, C-terminal basic amino acid (aa) substitution, was developed for increasing VLP and chimeric VLP (cVLP) production in this system. A series of mutants of five HPV types, including three L1 VLPs (6L1, 11L1, and 52L1) and two L1-L2 cVLPs (16L1-33L2, 58L1-16L2), were constructed. We found that most mutants exhibited higher protein expression in Sf9 cells, among which the yields of the superior mutants, 6L1CS4, 11L1CS3, 52L1m4∆N13CS1, 16L1-33L2 CS1, and 58L1-16L2 CS3, were up to 40, 35, 20, 35, and 60mg/L, which respectively increased by 4.2-, 7.3-, 5-, 2.5-, and 3.4-fold, and they also showed robust immunogenicity and great stabilities. Additionally, we found that the increased level of steady-state mRNA may play a crucial role in promoting L1 protein expression. Our results demonstrated that this novel method was cost-effective and can be used to reduce the production costs of L1 VLPs and L1-L2 cVLPs to develop broadly protective and affordable multivalent HPV vaccines.

Variability of the HCV core region and host genetic and epigenetic factors can predict the response to pegylated interferon/ribavirin therapy in genotype 1b hepatitis C patients from Serbia

Variations in the hepatitis C virus (HCV) core sequence have been related to disease progression and response to antiviral therapy. Previously we showed that the methylation status of RASSF1A and p16 genes, and IL28B genotypes affects the response to pegylated interferon/ribavirin (PEG-IFN/RBV) therapy. Herein we investigated whether amino acid (aa) substitutions in the HCV core region alone or in combination with IL28B genotypes and RASSF1A/p16 methylation affect the response to PEG-IFN/RBV therapy and liver disease progression. Among 29 examined patients, we found no association between single aa substitutions and response to therapy. However, we observed that patients with the HCV core aa substitution at position 75 and CT/TT IL28B genotypes were non-responders (NR), (P=0.023). Moreover, these patients had unmethylated RASSF1A. In contrast, most patients (75%) with aa substitutions at position 91 and CC IL28B genotype achieved sustained virologic response (SVR), (P=0.030), and 70% of them had methylated RASSF1A gene. Our results suggest that combined analysis of aa substitutions in the core protein, the IL28B rs12979860 polymorphism, and the methylation status of the RASSF1A gene may help in predicting treatment response to PEG-IFN/RBV in genotype 1b chronic hepatitis C patients.

Evolution of avian infectious bronchitis vaccine viruses reisolated from vaccinated broilers

Genetic and phenotypic diversity among infectious bronchitis virus (IBVs) is a key element of the epidemiology of avian infectious bronchitis. Virtually, all studies have focused on the evolution of wild viruses. Live-attenuated vaccine strains deserve a special attention. Two strains of IBV, 2T/17 and 16T/17, detected by quantitative reverse transcription polymerase chain reaction (qRT-PCR) in samples from H120 vaccinated chickens showing respiratory signs, were genetically analysed targeting the S1 gene. The total RNA was extracted by the classical Trizol method. The whole S1 gene was amplified in an end-point RT-PCR using a specific primer pair, before being sequenced. The obtained sequences were compared to those of vaccine strains and reference viruses and finally subjected to a phylogenetic analysis. The two isolates turned out to be vaccine strains as evidenced by the monophyletic cluster they formed with the H120 vaccine strain (100% bootstrap). Also, they shared with this vaccine strain 99.60% and 99.47% of nucleotides (nt) for strain 2T/17 and 16T/17 respectively. The sequence alignment revealed, for each of the two strains, one non-synonymous nt substitution outside the hypervariable regions: C56T and T39G leading to two amino acid (aa) substitutions: Alv19Val and Cys13Trp, respectively for the strain 2T/17 and 16T/17. The second hypervariable region of strain 2T/17 displayed 2 nt substitutions at positions 345 and 347, resulting in an aa substitution Lys116Thr. As for 16T/17 strain, it harboured two non-synonymous nt substitutions (T353G and T380G), leading to two aa substitutions: Val118Gly and Met127Arg. The chromatogram analysis revealed minor and major peaks at all positions where nt changes were observed. Nucleotides associated with minor peaks matched with those observed in the sequence of the H120 vaccine strain at the same positions, indicating subpopulation selection. Only the substitution T353G results in a mutation as evidenced by the absence of minor peak at this position. Viral subpopulations contained in IBV live-attenuated vaccines are susceptible to be in vivo selected and emerge as persistent pathotype with reversion to virulence. This draws attention to the need of monitoring the evolution of IBV vaccine strains.