An unusually high substitution rate in transplant-associated BK polyomavirus in vivo is further concentrated in HLA-C-bound viral peptides

Structural implications of BK polyomavirus sequence variations in the major viral capsid protein Vp1 and large T-antigen: a computational study

BK polyomavirus (BKPyV) is a double-stranded DNA virus causing nephropathy, hemorrhagic cystitis, and urothelial cancer in transplant patients. The BKPyV-encoded capsid protein Vp1 and large T-antigen (LTag) are key targets of neutralizing antibodies and cytotoxic T-cells, respectively. Our single-center data suggested that variability in Vp1 and LTag may contribute to failing BKPyV-specific immune control and impact vaccine design. We, therefore, analyzed all available entries in GenBank (1516 VP1; 742 LTAG) and explored potential structural effects using computational approaches. BKPyV-genotype (gt)1 was found in 71.18% of entries, followed by BKPyV-gt4 (19.26%), BKPyV-gt2 (8.11%), and BKPyV-gt3 (1.45%), but rates differed according to country and specimen type. Vp1-mutations matched a serotype different than the assigned one or were serotype-independent in 43%, 18% affected more than one amino acid. Notable Vp1-mutations altered antibody-binding domains, interactions with sialic acid receptors, or were predicted to change conformation. LTag-sequences were more conserved, with only 16 mutations detectable in more than one entry and without significant effects on LTag-structure or interaction domains. However, LTag changes were predicted to affect HLA-class I presentation of immunodominant 9mers to cytotoxic T-cells. These global data strengthen single center observations and specifically our earlier findings revealing mutant 9mer epitopes conferring immune escape from HLA-I cytotoxic T cells. We conclude that variability of BKPyV-Vp1 and LTag may have important implications for diagnostic assays assessing BKPyV-specific immune control and for vaccine design.

IMPORTANCE

Type and rate of amino acid variations in BKPyV may provide important insights into BKPyV diversity in human populations and an important step toward defining determinants of BKPyV-specific immunity needed to protect vulnerable patients from BKPyV diseases. Our analysis of BKPyV sequences obtained from human specimens reveals an unexpectedly high genetic variability for this double-stranded DNA virus that strongly relies on host cell DNA replication machinery with its proof reading and error correction mechanisms. BKPyV variability and immune escape should be taken into account when designing further approaches to antivirals, monoclonal antibodies, and vaccines for patients at risk of BKPyV diseases.

Evaluation of the inhibitory potential of bioactive compounds against SARS-CoV-2 by in silico approach

CONTEXT

The COVID-19 (coronavirus disease 19) pandemic brought on by the SARS-CoV-2 outbreak (severe acute respiratory syndrome coronavirus 2) has stimulated the exploration of various available chemical compounds that could be used to treat the infection. This has driven numerous researchers to investigate the antiviral potential of several bioactive compounds from medicinal plants due to their reduced adverse effects compared to chemicals. Some of the bioactive compounds used in folklore treatment strategies are reported as effective inhibitors against the proliferative and infective cycles of SARS-CoV-2. The secondary metabolites from plants are generally used to treat various diseases due to their intact medicinal properties. The present study analyzes the inhibitory potential of phytochemicals from medicinal plants like Sphaeranthus indicus, Lantana camara, and Nelumbo nucifera against SARS-CoV-2 by molecular docking.

METHODS

Ten druggable protein targets from SARS-CoV-2 are docked against the phytochemicals from the selected medicinal plants. The phytocompounds astragalin, isoquercetin, and 5-hydroxy-7-methoxy-6-c-glycosy flavone were found to have lower binding energy depicting their inhibitive potential compared with the reported inhibitors that are used in the treatment of SARS-CoV-2 infection. The phytocompounds found to have the least binding energy were selected for further analyses. To assess the compounds' potential as drugs, their ADMET characteristics were also examined. Sphaeranthus indicus, Lantana camara, and Nelumbo nucifera six possible compounds were separately screened for ADME and toxicity characteristics; then, the results were analyzed. To assess the impact of the phytocompound binding on the dynamics of SARS-CoV-2 ribonuclease protein NSP15, microsecond-level all atomistic molecular dynamics simulations were performed, and their dynamics were analyzed. Microsecond-level molecular dynamics simulations of both the ligands complexed with NSP15 revealed that the ligand induces allosteric effects on NSP15, which could lead to destabilization of NSP15 hexameric interface and loss of RNA binding. The low binding energy exhibited by the phytochemicals from Lantana camera, Sphaeranthus indicus, and Nelumbo nucifera against the protein targets of SARS-CoV-2 showed inhibitory potential by the selected molecules. Their predicted interference of the enzymes involved in the molecular mechanisms aiding the proliferation of SARS-CoV-2 indicated the inhibitive ability of the phytochemicals.

Polyomavirus BK Genome Comparison Shows High Genetic Diversity in Kidney Transplant Recipients Three Months after Transplantation

BK polyomavirus (BKPyV) is a human DNA virus generally divided into twelve subgroups based on the genetic diversity of Viral Protein 1 (VP1). BKPyV can cause polyomavirus-associated nephropathy (PVAN) after kidney transplantation. Detection of BKPyV DNA in blood (viremia) is a source of concern and increase in plasma viral load is associated with a higher risk of developing PVAN. In this work, we looked for possible associations of specific BKPyV genetic features with higher plasma viral load in kidney transplant patients. We analyzed BKPyV complete genome in three-month samples from kidney recipients who developed viremia during their follow-up period. BKPyV sequences were obtained by next-generation sequencing and were de novo assembled using the new BKAnaLite pipeline. Based on the data from 72 patients, we identified 24 viral groups with unique amino acid sequences: three in the VP1 subgroup IVc2, six in Ib1, ten in Ib2, one in Ia, and four in II. In none of the groups did the mean plasma viral load reach a statistically significant difference from the overall mean observed at three months after transplantation. Further investigation is needed to better understand the link between the newly described BKPyV genetic variants and pathogenicity in kidney transplant recipients.

The past, current and future epidemiological dynamic of SARS-CoV-2

SARS-CoV-2, the agent of the COVID-19 pandemic, emerged in late 2019 in China, and rapidly spread throughout the world to reach all continents. As the virus expanded in its novel human host, viral lineages diversified through the accumulation of around two mutations a month on average. Different viral lineages have replaced each other since the start of the pandemic, with the most successful Alpha, Delta and Omicron variants of concern (VoCs) sequentially sweeping through the world to reach high global prevalence. Neither Alpha nor Delta was characterized by strong immune escape, with their success coming mainly from their higher transmissibility. Omicron is far more prone to immune evasion and spread primarily due to its increased ability to (re-)infect hosts with prior immunity. As host immunity reaches high levels globally through vaccination and prior infection, the epidemic is expected to transition from a pandemic regime to an endemic one where seasonality and waning host immunization are anticipated to become the primary forces shaping future SARS-CoV-2 lineage dynamics. In this review, we consider a body of evidence on the origins, host tropism, epidemiology, genomic and immunogenetic evolution of SARS-CoV-2 including an assessment of other coronaviruses infecting humans. Considering what is known so far, we conclude by delineating scenarios for the future dynamic of SARS-CoV-2, ranging from the good-circulation of a fifth endemic 'common cold' coronavirus of potentially low virulence, the bad-a situation roughly comparable with seasonal flu, and the ugly-extensive diversification into serotypes with long-term high-level endemicity.

Quasispecies of SARS-CoV-2 revealed by single nucleotide polymorphisms (SNPs) analysis

New SARS-CoV-2 mutants have been continuously indentified with enhanced transmission ever since its outbreak in early 2020. As an RNA virus, SARS-CoV-2 has a high mutation rate due to the low fidelity of RNA polymerase. To study the single nucleotide polymorphisms (SNPs) dynamics of SARS-CoV-2, 158 SNPs with high confidence were identified by deep meta-transcriptomic sequencing, and the most common SNP type was C > T. Analyses of intra-host population diversity revealed that intra-host quasispecies' composition varies with time during the early onset of symptoms, which implicates viral evolution during infection. Network analysis of co-occurring SNPs revealed the most abundant non-synonymous SNP 22,638 in the S glycoprotein RBD region and 28,144 in the ORF8 region. Furthermore, SARS-CoV-2 variations differ in an individual's respiratory tissue (nose, throat, BALF, or sputum), suggesting independent compartmentalization of SARS-CoV-2 populations in patients. The positive selection analysis of the SARS-CoV-2 genome uncovered the positive selected amino acid G251V on ORF3a. Alternative allele frequency spectrum (AAFS) of all variants revealed that ORF8 could bear alternate alleles with high frequency. Overall, the results show the quasispecies' profile of SARS-CoV-2 in the respiratory tract in the first two months after the outbreak.

The Rise of SARS-CoV-2 Variants and the Role of Convalescent Plasma Therapy for Management of Infections

Novel therapies for the treatment of COVID-19 are continuing to emerge as the SARS-Cov-2 pandemic progresses. PCR remains the standard benchmark for initial diagnosis of COVID-19 infection, while advances in immunological profiling are guiding clinical treatment. The SARS-Cov-2 virus has undergone multiple mutations since its emergence in 2019, resulting in changes in virulence that have impacted on disease severity globally. The emergence of more virulent variants of SARS-Cov-2 remains challenging for effective disease control during this pandemic. Major variants identified to date include B.1.1.7, B.1.351; P.1; B.1.617.2; B.1.427; P.2; P.3; B.1.525; and C.37. Globally, large unvaccinated populations increase the risk of more and more variants arising. With successive waves of COVID-19 emerging, strategies that mitigate against community transmission need to be implemented, including increased vaccination coverage. For treatment, convalescent plasma therapy, successfully deployed during recent Ebola outbreaks and for H1N1 influenza, can increase survival rates and improve host responses to viral challenge. Convalescent plasma is rich with cytokines (IL-1β, IL-2, IL-6, IL-17, and IL-8), CCL2, and TNFα, neutralizing antibodies, and clotting factors essential for the management of SARS-CoV-2 infection. Clinical trials can inform and guide treatment policy, leading to mainstream adoption of convalescent therapy. This review examines the limited number of clinical trials published, to date that have deployed this therapy and explores clinical trials in progress for the treatment of COVID-19.

Al-Ameer H, E. Al Tamam S, Alashqar R, Al bawareed O, Gushchina Y, Samy Abousenna M, Ayyash AM, Sharab A, M. Alnaimat S, Imraish A, W. Akash M, R. Hadi N. Correlates of SARS-CoV-2 Variants on Deaths, Case Incidence and Case Fatality Ratio among the Continents for the Period of 1 December 2020 to 15 March 2021

The outbreak of coronavirus disease 2019 (COVID-19), by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has quickly developed into a worldwide pandemic. Mutations in the SARS-CoV-2 genome may affect various aspects of the disease including fatality ratio. In this study, 553,518 SARS-CoV-2 genome sequences isolated from patients from continents for the period 1 December 2020 to 15 March 2021 were comprehensively analyzed and a total of 82 mutations were identified concerning the reference sequence. In addition, associations between the mutations and the case fatality ratio (CFR), cases per million and deaths per million, were examined. The mutations having the highest frequencies among different continents were Spike_D614G and NSP12_P323L. Among the identified mutations, NSP2_T153M, NSP14_I42V and Spike_L18F mutations showed a positive correlation to CFR. While the NSP13_Y541C, NSP3_T73I and NSP3_Q180H mutations demonstrated a negative correlation to CFR. The Spike_D614G and NSP12_P323L mutations showed a positive correlation to deaths per million. The NSP3_T1198K, NS8_L84S and NSP12_A97V mutations showed a significant negative correlation to deaths per million. The NSP12_P323L and Spike_D614G mutations showed a positive correlation to the number of cases per million. In contrast, NS8_L84S and NSP12_A97V mutations showed a negative correlation to the number of cases per million. In addition, among the identified clades, none showed a significant correlation to CFR. The G, GR, GV, S clades showed a significant positive correlation to deaths per million. The GR and S clades showed a positive correlation to number of cases per million. The clades having the highest frequencies among continents were G, followed by GH and GR. These findings should be taken into consideration during epidemiological surveys of the virus and vaccine development.

In Vivo Generation of BK and JC Polyomavirus Defective Viral Genomes in Human Urine Samples Associated with Higher Viral Loads

Defective viral genomes (DVGs) are parasitic viral sequences containing point mutations, deletions, or duplications that might interfere with replication. DVGs are often associated with viral passage at high multiplicities of infection in culture systems but have been increasingly reported in clinical specimens. To date however, only RNA viruses have been shown to contain DVGs in clinical specimens. Here, using direct deep sequencing with multiple library preparation strategies and confirmatory digital droplet PCR (ddPCR) of urine samples taken from immunosuppressed individuals, we show that clinical BK polyomavirus (BKPyV) and JC polyomavirus (JCPyV) strains contain widespread genomic rearrangements across multiple loci that likely interfere with viral replication. BKPyV DVGs were derived from BKPyV genotypes Ia, Ib-1, and Ic. The presence of DVGs was associated with specimens containing higher viral loads but never reached clonality, consistent with a model of parasitized replication. These DVGs persisted during clinical infection as evidenced in two separate pairs of samples containing BK virus collected from the same individual up to 302 days apart. In a separate individual, we observed the generation of DVGs after a 57.5-fold increase in viral load. In summary, by extending the presence of DVGs in clinical specimens to DNA viruses, we demonstrate the ubiquity of DVGs in clinical virology.

IMPORTANCE Defective viral genomes (DVGs) can have a significant impact on the production of infectious virus particles. DVGs have only been identified in cultured viruses passaged at high multiplicities of infection and RNA viruses collected from clinical specimens; no DNA virus in the wild has been shown to contain DVGs. Here, we identified BK and JC polyomavirus DVGs in clinical urine specimens and demonstrated that these DVGs are more frequently identified in samples with higher viral loads. The strains containing DVGs had rearrangements throughout their genomes, with the majority affecting genes required for viral replication. Longitudinal analysis showed that these DVGs can persist during an infection but do not reach clonality within the chronically infected host. Our identification of polyomavirus DVGs suggests that these parasitic sequences exist across the many classes of viruses capable of causing human disease.

Genetics and genomics of SARS-CoV-2: A review of the literature with the special focus on genetic diversity and SARS-CoV-2 genome detection

The outbreak of 2019-novel coronavirus disease (COVID-19), caused by SARS-CoV-2, started in late 2019; in a short time, it has spread rapidly all over the world. Although some possible antiviral and anti-inflammatory medications are available, thousands of people are dying daily. Well-understanding of the SARS-CoV-2 genome is not only essential for the development of new treatments/vaccines, but it also can be used for improving the sensitivity and specificity of current approaches for virus detection. Accordingly, we reviewed the most critical findings related to the genetics of the SARS-CoV-2, with a specific focus on genetic diversity and reported mutations, molecular-based diagnosis assays, using interfering RNA technology for the treatment of patients, and genetic-related vaccination strategies. Additionally, considering the unanswered questions or uncertainties in these regards, different topics were discussed.

The Role of HLA and KIR Immunogenetics in BK Virus Infection after Kidney Transplantation

BK virus (BKV) is a polyomavirus with high seroprevalence in the general population with an unremarkable clinical presentation in healthy people, but a potential for causing serious complications in immunosuppressed transplanted patients. Reactivation or primary infection in kidney allograft recipients may lead to allograft dysfunction and subsequent loss. Currently, there is no widely accepted specific treatment for BKV infection and reduction of immunosuppressive therapy is the mainstay therapy. Given this and the sequential appearance of viruria-viremia-nephropathy, screening and early detection are of utmost importance. There are numerous risk factors associated with BKV infection including genetic factors, among them human leukocyte antigens (HLA) and killer cell immunoglobulin-like receptors (KIR) alleles have been shown to be the strongest so far. Identification of patients at risk for BKV infection would be useful in prevention or early action to reduce morbidity and progression to frank nephropathy. Assessment of risk involving HLA ligands and KIR genotyping of recipients in the pre-transplant or early post-transplant period might be useful in clinical practice. This review summarizes current knowledge of the association between HLA, KIR and BKV infection and potential future directions of research, which might lead to optimal utilization of these genetic markers.

No evidence for increased transmissibility from recurrent mutations in SARS-CoV-2

COVID-19 is caused by the coronavirus SARS-CoV-2, which jumped into the human population in late 2019 from a currently uncharacterised animal reservoir. Due to this recent association with humans, SARS-CoV-2 may not yet be fully adapted to its human host. This has led to speculations that SARS-CoV-2 may be evolving towards higher transmissibility. The most plausible mutations under putative natural selection are those which have emerged repeatedly and independently (homoplasies). Here, we formally test whether any homoplasies observed in SARS-CoV-2 to date are significantly associated with increased viral transmission. To do so, we develop a phylogenetic index to quantify the relative number of descendants in sister clades with and without a specific allele. We apply this index to a curated set of recurrent mutations identified within a dataset of 46,723 SARS-CoV-2 genomes isolated from patients worldwide. We do not identify a single recurrent mutation in this set convincingly associated with increased viral transmission. Instead, recurrent mutations currently in circulation appear to be evolutionary neutral and primarily induced by the human immune system via RNA editing, rather than being signatures of adaptation. At this stage we find no evidence for significantly more transmissible lineages of SARS-CoV-2 due to recurrent mutations.

Emergence of genomic diversity and recurrent mutations in SARS-CoV-2

SARS-CoV-2 is a SARS-like coronavirus of likely zoonotic origin first identified in December 2019 in Wuhan, the capital of China's Hubei province. The virus has since spread globally, resulting in the currently ongoing COVID-19 pandemic. The first whole genome sequence was published on January 5 2020, and thousands of genomes have been sequenced since this date. This resource allows unprecedented insights into the past demography of SARS-CoV-2 but also monitoring of how the virus is adapting to its novel human host, providing information to direct drug and vaccine design. We curated a dataset of 7666 public genome assemblies and analysed the emergence of genomic diversity over time. Our results are in line with previous estimates and point to all sequences sharing a common ancestor towards the end of 2019, supporting this as the period when SARS-CoV-2 jumped into its human host. Due to extensive transmission, the genetic diversity of the virus in several countries recapitulates a large fraction of its worldwide genetic diversity. We identify regions of the SARS-CoV-2 genome that have remained largely invariant to date, and others that have already accumulated diversity. By focusing on mutations which have emerged independently multiple times (homoplasies), we identify 198 filtered recurrent mutations in the SARS-CoV-2 genome. Nearly 80% of the recurrent mutations produced non-synonymous changes at the protein level, suggesting possible ongoing adaptation of SARS-CoV-2. Three sites in Orf1ab in the regions encoding Nsp6, Nsp11, Nsp13, and one in the Spike protein are characterised by a particularly large number of recurrent mutations (>15 events) which may signpost convergent evolution and are of particular interest in the context of adaptation of SARS-CoV-2 to the human host. We additionally provide an interactive user-friendly web-application to query the alignment of the 7666 SARS-CoV-2 genomes.

Amino acid sequence associated with bacteriophage recombination site helps to reveal genes potentially acquired through horizontal gene transfer

Background

Horizontal gene transfer, i.e. the acquisition of genetic material from nonparent organism, is considered an important force driving species evolution. Many cases of horizontal gene transfer from prokaryotes to eukaryotes have been registered, but no transfer mechanism has been deciphered so far, although viruses were proposed as possible vectors in several studies. In agreement with this idea, in our previous study we discovered that in two eukaryotic proteins bacteriophage recombination site (AttP) was adjacent to the regions originating via horizontal gene transfer. In one of those cases AttP site was present inside the introns of cysteine-rich repeats. In the present study we aimed to apply computational tools for finding multiple horizontal gene transfer events in large genome databases. For that purpose we used a sequence of cysteine-rich repeats to identify genes potentially acquired through horizontal transfer.

Results

HMMER remote similarity search significantly detected 382 proteins containing cysteine-rich repeats. All of them, except 8 sequences, belong to eukaryotes. In 124 proteins the presence of conserved structural domains was predicted. In spite of the fact that cysteine-rich repeats are found almost exclusively in eukaryotic proteins, many predicted domains are most common for prokaryotes or bacteriophages. Ninety-eight proteins out of 124 contain typical prokaryotic domains. In those cases proteins were considered as potentially originating via horizontal transfer. In addition, HHblits search revealed that two domains of the same fungal protein, Glycoside hydrolase and Peptidase M15, have high similarity with proteins of two different prokaryotic species, hinting at independent horizontal gene transfer events.

Conclusions

Cysteine-rich repeats in eukaryotic proteins are usually accompanied by conserved domains typical for prokaryotes or bacteriophages. These proteins, containing both cysteine-rich repeats, and characteristic prokaryotic domains, might represent multiple independent horizontal gene transfer events from prokaryotes to eukaryotes. We believe that the presence of bacteriophage recombination site inside cysteine-rich repeat coding sequence may facilitate horizontal genes transfer. Thus computational approach, described in the present study, can help finding multiple sequences originated from horizontal transfer in eukaryotic genomes.

Molecular epidemiology and phylodynamics of goose haemorrhagic polyomavirus

Goose haemorrhagic polyomavirus (GHPV, or Anser anser polyomavirus 1) is a small dsDNA virus of the Polyomaviridae family. The virus infects the internal organs causing haemorrhagic nephritis and enteritis of geese that may be fatal for goslings. In this study, GHPV positivity was examined in goose and duck samples collected in Hungary between 2005 and 2019. In this period, 384 of the investigated 1,111 specimens were diagnosed as GHPV-positive by PCR assay. Twenty-two GHPV genomes were sequenced and subjected to phylogenetic and evolutionary analysis. Based on the sequence data, the mean evolutionary rates were estimated 6.57 × 10^-6 -5.82 × 10^-5  s/s/y for both GHPV complete genomes and individual genes, with negative selection acting on each gene. When GHPV VP1 sequences originating from wild birds were also included in the analyses, the nt and aa mutations inflated the substitution rate to 1.54 × 10^-4  s/s/y that may imply adaptation of the virus to novel host species. Our data suggested the co-circulation of various GHPV strains in Hungarian goose farms; the source of these may be persistently infected domesticated or migratory wild birds. Detection and characterization of GHPV in wild birds and domestic waterfowls may help to elaborate new strategies for more effective disease control and prevention.

BK polyomavirus diversity-Why viral variation matters

BK polyomavirus (BKPyV or BKV) is a non-enveloped, circular double-stranded DNA virus that may exceed 80% seroprevalence in adults. BKV infection typically occurs during childhood, and the majority of adults are latently infected. While BKV infection is rarely associated with clinical disease in most individuals, in immunosuppressed individuals, reactivation may cause kidney (BK-associated nephropathy) or bladder (hemorrhagic cystitis and ureteral stenosis) injury. No antiviral therapies have been approved for the treatment of BKV infection. Reducing immunosuppression is the most effective therapy, although this is not feasible in many patients. Thus, a robust understanding of viral pathogenesis and viral diversity remains important for the development of future therapeutic strategies. Studies of BKV diversity are quite sparse compared to other common viral infections; thus, much of our understanding of BVK variability and evolution relies heavily analogous studies of other viruses such as HIV or viral hepatitis. We provide a comprehensive review of BKV diversity at the population and individual level with careful consideration of how viral variability may impact viral replication, pathogenesis, tropism, and protein function. We also discuss a number of outstanding questions related to BK virus diversity that should be explored rigorously in future studies.

Viral quasispecies

Viral quasispecies refers to a population structure that consists of extremely large numbers of variant genomes, termed mutant spectra, mutant swarms or mutant clouds. Fueled by high mutation rates, mutants arise continually, and they change in relative frequency as viral replication proceeds. The term quasispecies was adopted from a theory of the origin of life in which primitive replicons) consisted of mutant distributions, as found experimentally with present day RNA viruses. The theory provided a new definition of wild type, and a conceptual framework for the interpretation of the adaptive potential of RNA viruses that contrasted with classical studies based on consensus sequences. Standard clonal analyses and deep sequencing methodologies have confirmed the presence of myriads of mutant genomes in viral populations, and their participation in adaptive processes. The quasispecies concept applies to any biological entity, but its impact is more evident when the genome size is limited and the mutation rate is high. This is the case of the RNA viruses, ubiquitous in our biosphere, and that comprise many important pathogens. In virology, quasispecies are defined as complex distributions of closely related variant genomes subjected to genetic variation, competition and selection, and that may act as a unit of selection. Despite being an integral part of their replication, high mutation rates have an upper limit compatible with inheritable information. Crossing such a limit leads to RNA virus extinction, a transition that is the basis of an antiviral design termed lethal mutagenesis.