Spatial and temporal distribution of SARS-CoV-2 diversity circulating in wastewater.

Wastewater-based epidemiology (WBE) has proven to be an effective tool for epidemiological surveillance of SARS-CoV-2 during the current COVID-19 pandemic. Furthermore, combining WBE together with high-throughput sequencing techniques can be useful for the analysis of SARS-CoV-2 viral diversity present in a given sample. The present study focuses on the genomic analysis of SARS-CoV-2 in 76 sewage samples collected during the three epidemiological waves that occurred in Spain from 14 wastewater treatment plants distributed throughout the country. The results obtained demonstrate that the metagenomic analysis of SARS-CoV-2 in wastewater allows the detection of mutations that define the B.1.1.7 lineage and the ability of the technique to anticipate the detection of certain mutations before they are detected in clinical samples. The study proves the usefulness of sewage sequencing to track Variants of Concern that can complement clinical testing to help in decision-making and in the analysis of the evolution of the pandemic.

Exploring the Diversity of the Human Blood Virome.

Metagenomics is greatly improving our ability to discover new viruses, as well as their possible associations with disease. However, metagenomics has also changed our understanding of viruses in general. The vast expansion of currently known viral diversity has revealed a large fraction of non-pathogenic viruses, and offers a new perspective in which viruses function as important components of many ecosystems. In this vein, studies of the human blood virome are often motivated by the search for new viral diseases, especially those associated with blood transfusions. However, these studies have revealed the common presence of apparently non-pathogenic viruses in blood, particularly human anelloviruses and, to a lower extent, human pegiviruses (HPgV). To shed light on the diversity of the human blood virome, we subjected pooled plasma samples from 587 healthy donors in Spain to a viral enrichment protocol, followed by massive parallel sequencing. This showed that anelloviruses were clearly the major component of the blood virome and showed remarkable diversity. In total, we assembled 332 complete or near-complete anellovirus genomes, 50 of which could be considered new species. HPgV was much less frequent, but we, nevertheless, recovered 17 different isolates that we subsequently used for characterizing the diversity of this virus. In-depth investigation of the human blood virome should help to elucidate the ecology of these viruses, and to unveil potentially associated diseases.

Hepatitis E virus in lettuce and water samples: A method-comparison study.

The hepatitis E virus (HEV), which is an increasing cause of acute viral hepatitis in Europe, is a zoonotic virus that is mainly transmitted through contaminated water, consumption of raw or undercooked meat from pigs or wild boar, blood transfusion, and organ transplantation. Although the role of HEV transmission through contaminated produce has not been confirmed, the presence of HEV has been reported in irrigation waters and in vegetables. The present study used a World Health Organization (WHO) international standard and clinical samples to evaluate the performance characteristics of three RT-qPCR assays for detection and quantification of HEV. Two of the evaluated assays provided good analytical sensitivity, as 250 international units (IU) per ml could be detected. Then, experiments focused on evaluating the elution conditions suitable for HEV release from vegetables, with the method proposed by the ISO 15216:2017 selected for evaluation in three types of fresh vegetables. The concentration method proposed by the ISO 15216:2017 combined with the RT-qPCR described by Schlosser et al. (2014) resulted in average HEV recoveries of 1.29%, 0.46%, and 3.95% in lettuce, spinach, and pepper, respectively, with an average detection limit of 1.47 × 105 IU/25 g. In naturally contaminated samples, HEV was detected in sewage only (10/14), while no detection was reported in lettuce (0/36) or in irrigation water samples (0/24).

EcPV2 DNA in equine genital squamous cell carcinomas and normal genital mucosa.

Squamous cell carcinoma (SCC) represents the most common genital malignant tumor in horses. Similar to humans, papillomaviruses (PVs) have been proposed as etiological agents and recently Equine papillomavirus type 2 (EcPV2) has been identified in a subset of genital SCCs. The goals of this study were (1) to determine the prevalence of EcPV2 DNA in tissue samples from equine genital SCCs, penile intraepithelial neoplasia (PIN) and penile papillomas, using EcPV2-specific PCR, (2) to examine the prevalence of latent EcPV2 infection in healthy genital mucosa and (3) to determine genetic variability within EcPV2 and to disentangle phylogenetic relationships of EcPV2 among PVs. EcPV2 DNA was detected in all but one penile SCC (15/16), in all PIN lesions (8/8) and penile papillomas (4/4). Additionally, EcPV2 DNA was demonstrated in one of two metastasized lymph nodes, one contact metastasis in the mouth, two vaginal and one anal lesion. In healthy horses, EcPV2 DNA was detected in 10% (4/39) of penile swabs but in none of vulvovaginal swabs (0/20). This study confirms the presence of EcPV2 DNA in equine genital SCCs and shows its involvement in anal lesions, a lymph node and contact metastases. Latent EcPV2 presence was also shown in normal male genital mucosa. We found that different EcPV2 variants cocirculate among horses and that EcPV2 is related to the Delta+Zeta PVs and is only a very distant relative of high-risk human PVs causing genital cancer. Thus, similar viral tropism and similar malignant outcome of the infection do not imply close evolutionary relationship.

Enterovirus co-infections and onychomadesis after hand, foot, and mouth disease, Spain, 2008.

Hand, foot, and mouth disease (HFMD), a common disease caused by enteroviruses (EVs), usually affects children. Clustered and sporadic HFMD cases, followed by onychomadesis (nail shedding), occurred during summer and fall 2008 in Valencia, Spain. Fecal samples from onychomadesis patients, who did or did not have previous HFMD, and from healthy children exposed to onychomadesis patients tested positive for EV. The complete viral protein 1 capsid gene sequence was obtained for typing and phylogenetic analysis. Two EV serotypes, coxsackievirus A10 and coxsackievirus B1 (CVB1), were mainly detected as a monoinfection or co-infection in a childcare center where an onychomadesis outbreak occurred. On the basis of our results, and detection of CVB1 in 2 other contemporary onychomadesis outbreaks in childcare centers in Spain, we propose that mixed infection of an EV serotype that causes HFMD, plus the serotype CVB1, could explain the emergence after HFMD of onychomadesis, a rare and late complication.

Modular organizations of novel cetacean papillomaviruses.

The phylogenetic position of cetacean papillomaviruses (PVs: Omikron-PVs and Upsilon-PVs) varies depending on the region of the genome analysed. They cluster together with Alpha-PVs when analysing early genes and with Xi-PVs and Phi-PVs when analysing late genes. We cloned and sequenced the complete genomes of five novel PVs, sampled from genital and oesophageal lesions of free-ranging cetaceans: Delphinus delphis (DdPV1), Lagenorhynchus acutus (TtPV3 variant), and Phocoena phocoena (PphPV1, PphPV2, and PphPV3). Using Maximum Likelihood and Bayesian approaches, all cetacean PVs constituted a monophyletic group with Alpha-, Omega-, and Dyodelta-PVs as inferred from E1-E2 early genes analyses, thus matching the shared phenotype of mucosal tropism. However, cetacean PVs, with the exception of PphPV3, were the closest relatives of Xi-PVs and Phi-PVs in L2-L1 late genes analyses, isolated from cow and goat, thus reflecting the close relationship between Cetacea and Artiodactyla. Our results are compatible with a recombination between ancestral PVs infecting the Cetartiodactyla lineage. Our study supports a complex evolutionary scenario with multiple driving forces for PV diversification, possibly including recombination and also interspecies transmission.

The role of positive selection in hepatitis C virus.

Hepatitis C virus (HCV) is a major health problem worldwide, infecting an estimated 170 million people. In this study, we have employed a large data set of sequences (14,654 sequences from between 25 and 100 clone sequences per analyzed region and per patient) from 67 patients infected with HCV genotype 1 (23 subtype 1a and 44 subtype 1b). For all patients, a sample prior to combined therapy with alpha interferon plus ribavirin was available, whereas for some patients additional samples after 6 or 12 months of treatment were also available. Twenty-seven patients responded to treatment (12 subtype 1a and 15 subtype 1b) and forty patients did not respond to treatment (11 subtype 1a vs. 29 subtype 1b). Two regions of the HCV genome were analyzed, one compressing the hypervariable regions (HVR1, HVR2 and HVR3) of the envelope 2 glycoprotein and another one including the interferon sensitive determining region (ISDR) and the V3 domain of the NS5A protein. Previously (Cuevas, J.M., Torres-Puente, M., Jiménez-Hernández, N., Bracho, M.A., García-Robles, I., Wrobel, B., Carnicer, F., del Olmo, J., Ortega, E., Moya, A., González-Candelas, F., 2008b. Genetic variability of hepatitis C virus before and after combined therapy of interferon plus ribavirin. Plos One 3 (8), e3058), several amino acid positions in both regions analyzed were detected to be under positive selection. Here, we have compared the amino acid composition of each positively selected position between responder and non-responder patients for both subtypes. If we exclude some non-conclusive cases, no clear differences were detected in any case. In conclusion, identifying specific positions as completely discriminatory of treatment response seems to be a difficult task. Our results, in concordance with previous studies, suggest that HCV evasion strategies are more likely based on a global increased variability, which would yield combinations of mutations with an increased resistance, than on the fixation of specific amino acids conferring resistance to antiviral treatment or immune response. In this sense, the particular systemic response from each patient could play an essential role in determining the outcome of the antiviral treatment.

Combined therapy of interferon plus ribavirin promotes multiple adaptive solutions in hepatitis C virus.

Hepatitis C virus (HCV) presents several regions involved potentially in evading antiviral treatment and host immune system. Two regions, known as PKR-BD and V3 domains, have been proposed to be involved in resistance to interferon. Additionally, hypervariable regions in the envelope E2 glycoprotein are also good candidates to participate in evasion from the immune system. In this study, we have used a cohort of 22 non-responder patients to combined therapy (interferon alpha-2a plus ribavirin) for which samples obtained just before initiation of therapy and after 6 or/and 12 months of treatment were available. A range of 25-100 clones per patient, genome region and time sample were obtained. The predominant amino acid sequences for each time sample and patient were determined. Next, the sequences of the PKR-BD and V3 domains and the hypervariable regions from different time samples were compared for each patient. The highest levels of variability were detected at the three hypervariable regions of the E2 protein and, to a lower extent, at the V3 domain of the NS5A protein. However, no clear patterns of adaptation to the host immune system or to antiviral treatment were detected. In summary, although high levels of variability are correlated to viral adaptive response, antiviral treatment does not seem to promote convergent adaptive changes. Consequently, other regions must be involved in evasion strategies likely based on a combination of multiple mechanisms, in which pools of changes along the HCV genome could confer viruses the ability to overcome strong selective pressures.

Complete genome of a European hepatitis C virus subtype 1g isolate: phylogenetic and genetic analyses.

BACKGROUND: Hepatitis C virus isolates have been classified into six main genotypes and a variable number of subtypes within each genotype, mainly based on phylogenetic analysis. Analyses of the genetic relationship among genotypes and subtypes are more reliable when complete genome sequences (or at least the full coding region) are used; however, so far 31 of 80 confirmed or proposed subtypes have at least one complete genome available. Of these, 20 correspond to confirmed subtypes of epidemic interest. RESULTS: We present and analyse the first complete genome sequence of a HCV subtype 1g isolate. Phylogenetic and genetic distance analyses reveal that HCV-1g is the most divergent subtype among the HCV-1 confirmed subtypes. Potential genomic recombination events between genotypes or subtype 1 genomes were ruled out. We demonstrate phylogenetic congruence of previously deposited partial sequences of HCV-1g with respect to our sequence. CONCLUSION: In light of this, we propose changing the current status of its subtype-specific designation from provisional to confirmed.

Hepatitis C virus and the controversial role of the interferon sensitivity determining region in the response to interferon treatment.

The degree of variability of the interferon sensitivity determining region (ISDR) in the hepatitis C virus (HCV) genome has been postulated to predict the response to interferon therapy, mainly in patients infected with subtype 1b, although this prediction has been the subject of a long controversy. This prediction has been tested by analyzing a cohort of 67 Spanish patients infected with HCV genotype 1, 23 of which were infected with subtype 1a and 44 with subtype 1b. A sample previous to therapy with alpha-interferon plus ribavirin was obtained and several clones (between 25 and 96) including the ISDR were sequenced from each patient. A significant correlation between mutations at the ISDR and response to treatment for subtype 1b patients, but not for those infected with subtype 1a, has been detected. Although the results suggest that the same relationship holds true for subtype 1a, lack of statistical power because of the small sample size of this subtype prevented firmer conclusions. However, identical ISDR sequences were found in responder and non-responder patients, suggesting that the stability of the ISDR sequence can occasionally help HCV to evade interferon therapy, although this is not a sufficient condition. More complex interactions, including the ISDR or not, are likely to exist and govern the HCV response to interferon treatment.

Contribution of insertions and deletions to the variability of hepatitis C virus populations.

Little is known about the potential effects of insertions and deletions (indels) on the evolutionary dynamics of hepatitis C virus (HCV). In fact, the consequences of indels on antiviral treatment response are a field of investigation completely unexplored. Here, an extensive sequencing project was undertaken by cloning and sequencing serum samples from 25 patients infected with HCV subtype 1a and 48 patients with subtype 1b. For 23 patients, samples obtained after treatment with alpha interferon plus ribavirin were also available. Two genome fragments containing the hypervariable regions in the envelope 2 glycoprotein and the PKR-BD domain in NS5A were sequenced, yielding almost 16 000 sequences. Our results show that insertions are quite rare, but they are often present in biologically relevant domains of the HCV genome. Moreover, their frequency distributions between different time samples reflect the quasispecies dynamics of HCV populations. Deletions seem to be subject to negative selection.