Identification of Leishmania infantum Epidemiology, Drug Resistance and Pathogenicity Biomarkers with Nanopore Sequencing

The emergence of drug-resistant strains of the parasite Leishmania infantum infecting dogs and humans represents an increasing threat. L. infantum genomes are complex and unstable with extensive structural variations, ranging from aneuploidies to multiple copy number variations (CNVs). These CNVs have recently been validated as biomarkers of Leishmania concerning virulence, tissue tropism, and drug resistance. As a proof-of-concept to develop a novel diagnosis platform (LeishGenApp), four L. infantum samples from humans and dogs were nanopore sequenced. Samples were epidemiologically typed within the Mediterranean L. infantum group, identifying members of the JCP5 and non-JCP5 subgroups, using the conserved region (CR) of the maxicircle kinetoplast. Aneuploidies were frequent and heterogenous between samples, yet only chromosome 31 tetrasomy was common between all the samples. A high frequency of aneuploidies was observed for samples with long passage history (MHOM/TN/80/IPT-1), whereas fewer were detected for samples maintained in vivo (MCRI/ES/2006/CATB033). Twenty-two genes were studied to generate a genetic pharmacoresistance profile against miltefosine, allopurinol, trivalent antimonials, amphotericin, and paromomycin. MHOM/TN/80/IPT-1 and MCRI/ES/2006/CATB033 displayed a genetic profile with potential resistance against miltefosine and allopurinol. Meanwhile, MHOM/ES/2016/CATB101 and LCAN/ES/2020/CATB102 were identified as potentially resistant against paromomycin. All four samples displayed a genetic profile for resistance against trivalent antimonials. Overall, this proof-of-concept revealed the potential of nanopore sequencing and LeishGenApp for the determination of epidemiological, drug resistance, and pathogenicity biomarkers in L. infantum.

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.

Considerable escape of SARS-CoV-2 Omicron to antibody neutralization

The SARS-CoV-2 Omicron variant was first identified in November 2021 in Botswana and South Africa^1-3. It has since spread to many countries and is expected to rapidly become dominant worldwide. The lineage is characterized by the presence of around 32 mutations in spike-located mostly in the N-terminal domain and the receptor-binding domain-that may enhance viral fitness and enable antibody evasion. Here we isolated an infectious Omicron virus in Belgium from a traveller returning from Egypt. We examined its sensitivity to nine monoclonal antibodies that have been clinically approved or are in development⁴, and to antibodies present in 115 serum samples from COVID-19 vaccine recipients or individuals who have recovered from COVID-19. Omicron was completely or partially resistant to neutralization by all monoclonal antibodies tested. Sera from recipients of the Pfizer or AstraZeneca vaccine, sampled five months after complete vaccination, barely inhibited Omicron. Sera from COVID-19-convalescent patients collected 6 or 12 months after symptoms displayed low or no neutralizing activity against Omicron. Administration of a booster Pfizer dose as well as vaccination of previously infected individuals generated an anti-Omicron neutralizing response, with titres 6-fold to 23-fold lower against Omicron compared with those against Delta. Thus, Omicron escapes most therapeutic monoclonal antibodies and, to a large extent, vaccine-elicited antibodies. However, Omicron is neutralized by antibodies generated by a booster vaccine dose.

Considerable escape of SARS-CoV-2 Omicron to antibody neutralization

The SARS-CoV-2 Omicron variant was first identified in November 2021 in Botswana and South Africa^1-3. It has since then spread to many countries and is expected to rapidly become dominant worldwide. The lineage is characterized by the presence of about 32 mutations in the spike, located mostly in the N-terminal domain (NTD) and the receptor binding domain (RBD), which may enhance viral fitness and allow antibody evasion. Here, we isolated an infectious Omicron virus in Belgium, from a traveller returning from Egypt. We examined its sensitivity to 9 monoclonal antibodies (mAbs) clinically approved or in development⁴, and to antibodies present in 115 sera from COVID-19 vaccine recipients or convalescent individuals. Omicron was totally or partially resistant to neutralization by all mAbs tested. Sera from Pfizer or AstraZeneca vaccine recipients, sampled 5 months after complete vaccination, barely inhibited Omicron. Sera from COVID-19 convalescent patients collected 6 or 12 months post symptoms displayed low or no neutralizing activity against Omicron. Administration of a booster Pfizer dose as well as vaccination of previously infected individuals generated an anti-Omicron neutralizing response, with titers 6 to 23 fold lower against Omicron than against Delta. Thus, Omicron escapes most therapeutic monoclonal antibodies and to a large extent vaccine-elicited antibodies. Omicron remains however neutralized by antibodies generated by a booster vaccine dose.

Exploration of the Ixodes ricinus virosphere unveils an extensive virus diversity including novel coltiviruses and other reoviruses

Recent metagenomics studies have revealed several tick species to host a variety of previously undiscovered RNA viruses. Ixodes ricinus, which is known to be a vector for many viral, bacterial, and protozoan pathogens, is the most prevalent tick species in Europe. For this study, we decided to investigate the virosphere of Belgian I. ricinus ticks. High-throughput sequencing of tick pools collected from six different sampling sites revealed the presence of viruses belonging to many different viral orders and families, including Mononegavirales, Bunyavirales, Partitiviridae, and Reoviridae. Of particular interest was the detection of several new reoviruses, two of which cluster together with members of the genus Coltivirus. This includes a new strain of Eyach virus, a known causative agent of tick-borne encephalitis. All genome segments of this new strain are highly similar to those of previously published Eyach virus genomes, except for the fourth segment, encoding VP4, which is markedly more dissimilar, potentially indicating the occurrence of a genetic reassortment. Further polymerase chain reaction-based screening of over 230 tick pools for 14 selected viruses showed that most viruses could be found in all six sampling sites, indicating the wide spread of these viruses throughout the Belgian tick population. Taken together, these results illustrate the role of ticks as important virus reservoirs, highlighting the need for adequate tick control measures.

Human Diversity of Killer Cell Immunoglobulin-Like Receptors and Human Leukocyte Antigen Class I Alleles and Ebola Virus Disease Outcomes

We investigated the genetic profiles of killer cell immunoglobulin-like receptors (KIRs) in Ebola virus–infected patients. We studied the relationship between KIR–human leukocyte antigen (HLA) combinations and the clinical outcomes of patients with Ebola virus disease (EVD). We genotyped KIRs and HLA class I alleles using DNA from uninfected controls, EVD survivors, and persons who died of EVD. The activating 2DS4–003 and inhibitory 2DL5 genes were significantly more common among persons who died of EVD; 2DL2 was more common among survivors. We used logistic regression analysis and Bayesian modeling to identify 2DL2, 2DL5, 2DS4–003, HLA-B-Bw4-Thr, and HLA-B-Bw4-Ile as probably having a significant relationship with disease outcome. Our findings highlight the importance of innate immune response against Ebola virus and show the association between KIRs and the clinical outcome of EVD.

A Phylodynamic Workflow to Rapidly Gain Insights into the Dispersal History and Dynamics of SARS-CoV-2 Lineages

Since the start of the COVID-19 pandemic, an unprecedented number of genomic sequences of SARS-CoV-2 have been generated and shared with the scientific community. The unparalleled volume of available genetic data presents a unique opportunity to gain real-time insights into the virus transmission during the pandemic, but also a daunting computational hurdle if analyzed with gold-standard phylogeographic approaches. To tackle this practical limitation, we here describe and apply a rapid analytical pipeline to analyze the spatiotemporal dispersal history and dynamics of SARS-CoV-2 lineages. As a proof of concept, we focus on the Belgian epidemic, which has had one of the highest spatial densities of available SARS-CoV-2 genomes. Our pipeline has the potential to be quickly applied to other countries or regions, with key benefits in complementing epidemiological analyses in assessing the impact of intervention measures or their progressive easement.

NCBI's Virus Discovery Codeathon: Building "FIVE" -The Federated Index of Viral Experiments API Index

Viruses represent important test cases for data federation due to their genome size and the rapid increase in sequence data in publicly available databases. However, some consequences of previously decentralized (unfederated) data are lack of consensus or comparisons between feature annotations. Unifying or displaying alternative annotations should be a priority both for communities with robust entry representation and for nascent communities with burgeoning data sources. To this end, during this three-day continuation of the Virus Hunting Toolkit codeathon series (VHT-2), a new integrated and federated viral index was elaborated. This Federated Index of Viral Experiments (FIVE) integrates pre-existing and novel functional and taxonomy annotations and virus–host pairings. Variability in the context of viral genomic diversity is often overlooked in virus databases. As a proof-of-concept, FIVE was the first attempt to include viral genome variation for HIV, the most well-studied human pathogen, through viral genome diversity graphs. As per the publication of this manuscript, FIVE is the first implementation of a virus-specific federated index of such scope. FIVE is coded in BigQuery for optimal access of large quantities of data and is publicly accessible. Many projects of database or index federation fail to provide easier alternatives to access or query information. To this end, a Python API query system was developed to enhance the accessibility of FIVE.

Preparedness needs research: How fundamental science and international collaboration accelerated the response to COVID-19

The first cluster of patients suffering from coronavirus disease 2019 (COVID-19) was identified on December 21, 2019, and as of July 29, 2020, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections have been linked with 664,333 deaths and number at least 16,932,996 worldwide. Unprecedented in global societal impact, the COVID-19 pandemic has tested local, national, and international preparedness for viral outbreaks to the limits. Just as it will be vital to identify missed opportunities and improve contingency planning for future outbreaks, we must also highlight key successes and build on them. Concomitant to the emergence of a novel viral disease, there is a 'research and development gap' that poses a threat to the overall pace and quality of outbreak response during its most crucial early phase. Here, we outline key components of an adequate research response to novel viral outbreaks using the example of SARS-CoV-2. We highlight the exceptional recent progress made in fundamental science, resulting in the fastest scientific response to a major infectious disease outbreak or pandemic. We underline the vital role of the international research community, from the implementation of diagnostics and contact tracing procedures to the collective search for vaccines and antiviral therapies, sustained by unique information sharing efforts.

BKTyper: Free Online Tool for Polyoma BK Virus VP1 and NCCR Typing

Human BK polyomavirus (BKPyV) prevalence has been increasing due to the introduction of more potent immunosuppressive agents in transplant recipients, and its clinical interest. BKPyV has been linked mostly to polyomavirus-associated hemorrhagic cystitis, in allogenic hematopoietic stem cell transplant, and polyomavirus-associated nephropathy in kidney transplant patients. BKPyV is a circular double-stranded DNA virus that encodes for seven proteins, of which Viral Protein 1 (VP1), the major structural protein, has been extensively used for genotyping. BKPyV also contains the noncoding control region (NCCR), configured by five repeat blocks (OPQRS) known to be highly repetitive and diverse, and linked to viral infectivity and replication. BKPyV genetic diversity has been mainly studied based on the NCCR and VP1, due to the high occurrence of BKPyV-associated diseases in transplant patients and their clinical implications. Here BKTyper is presented, a free online genotyper for BKPyV, based on a VP1 genotyping and a novel algorithm for NCCR block identification. VP1 genotyping is based on a modified implementation of the BK typing and grouping regions (BKTGR) algorithm, providing a maximum-likelihood phylogenetic tree using a custom internal BKPyV database. Novel NCCR block identification relies on a minimum of 12-bp motif recognition and a novel sorting algorithm. A graphical representation of the OPQRS block organization is provided.

Accounting for population structure reveals ambiguity in the Zaire Ebolavirus reservoir dynamics

Ebolaviruses pose a substantial threat to wildlife populations and to public health in Africa. Evolutionary analyses of virus genome sequences can contribute significantly to elucidate the origin of new outbreaks, which can help guide surveillance efforts. The reconstructed between-outbreak evolutionary history of Zaire ebolavirus so far has been highly consistent. By removing the confounding impact of population growth bursts during local outbreaks on the free mixing assumption that underlies coalescent-based demographic reconstructions, we find-contrary to what previous results indicated-that the circulation dynamics of Ebola virus in its animal reservoir are highly uncertain. Our findings also accentuate the need for a more fine-grained picture of the Ebola virus diversity in its reservoir to reliably infer the reservoir origin of outbreak lineages. In addition, the recent appearance of slower-evolving variants is in line with latency as a survival mechanism and with bats as the natural reservoir host.

NCBI's Virus Discovery Hackathon: Engaging Research Communities to Identify Cloud Infrastructure Requirements

A wealth of viral data sits untapped in publicly available metagenomic data sets when it might be extracted to create a usable index for the virological research community. We hypothesized that work of this complexity and scale could be done in a hackathon setting. Ten teams comprised of over 40 participants from six countries, assembled to create a crowd-sourced set of analysis and processing pipelines for a complex biological data set in a three-day event on the San Diego State University campus starting 9 January 2019. Prior to the hackathon, 141,676 metagenomic data sets from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) were pre-assembled into contiguous assemblies (contigs) by NCBI staff. During the hackathon, a subset consisting of 2953 SRA data sets (approximately 55 million contigs) was selected, which were further filtered for a minimal length of 1 kb. This resulted in 4.2 million (Mio) contigs, which were aligned using BLAST against all known virus genomes, phylogenetically clustered and assigned metadata. Out of the 4.2 Mio contigs, 360,000 contigs were labeled with domains and an additional subset containing 4400 contigs was screened for virus or virus-like genes. The work yielded valuable insights into both SRA data and the cloud infrastructure required to support such efforts, revealing analysis bottlenecks and possible workarounds thereof. Mainly: (i) Conservative assemblies of SRA data improves initial analysis steps; (ii) existing bioinformatic software with weak multithreading/multicore support can be elevated by wrapper scripts to use all cores within a computing node; (iii) redesigning existing bioinformatic algorithms for a cloud infrastructure to facilitate its use for a wider audience; and (iv) a cloud infrastructure allows a diverse group of researchers to collaborate effectively. The scientific findings will be extended during a follow-up event. Here, we present the applied workflows, initial results, and lessons learned from the hackathon.

Human cytomegalovirus genomics and transcriptomics through the lens of next-generation sequencing: revision and future challenges

The human cytomegalovirus (HCMV) genome was sequenced by hierarchical shotgun almost 30 years ago. Over these years, low and high passaged strains have been sequenced, improving, albeit still far from complete, the understanding of the coding potential, expression dynamics and diversity of wild-type HCMV strains. Next-generation sequencing (NGS) platforms have enabled a huge advancement, facilitating the comparison of differentially passaged strains, challenging diagnostics and research based on a single or reduced gene set genotyping. In addition, it allowed to link genetic features to different viral phenotypes as for example, correlating large genomic re-arrangements to viral attenuation or different mutations to antiviral resistance and cell tropism. NGS platforms provided the first high-resolution experiments to HCMV dynamics, allowing the study of intra-host viral population structures and the description of rare transcriptional events. Long-read sequencing has recently become available, helping to identify new genomic re-arrangements, partially accounting for the genetic variability displayed in clinical isolates, as well as, in changing the understanding of the HCMV transcriptome. Better knowledge of the transcriptome resulted in a vast number of new splicing events and alternative transcripts, although most of them still need additional validation. This review summarizes the sequencing efforts reached so far, discussing its approaches and providing a revision and new nuances on HCMV sequence variability in the sequencing field.