Identification and characterization of novel alphacoronaviruses in Tadarida brasiliensis (Chiroptera, Molossidae) from Argentina: insights into recombination as a mechanism favoring bat coronavirus cross-species transmission

Bats are reservoirs of various coronaviruses that can jump between bat species or other mammalian hosts, including humans. This article explores coronavirus infection in three bat species (Tadarida brasiliensis, Eumops bonariensis, and Molossus molossus) of the family Molossidae from Argentina using whole viral metagenome analysis. Fecal samples of 47 bats from three semiurban or highly urbanized areas of the province of Santa Fe were investigated. After viral particle enrichment, total RNA was sequenced using the Illumina NextSeq 550 instrument; the reads were assembled into contigs and taxonomically and phylogenetically analyzed. Three novel complete Alphacoronavirus (AlphaCoV) genomes (Tb1-3) and two partial sequences were identified in T. brasiliensis (Tb4-5), and an additional four partial sequences were identified in M. molossus (Mm1-4). Phylogenomic analysis showed that the novel AlphaCoV clustered in two different lineages distinct from the 15 officially recognized AlphaCoV subgenera. Tb2 and Tb3 isolates appeared to be variants of the same virus, probably involved in a persistent infectious cycle within the T. brasiliensis colony. Using recombination analysis, we detected a statistically significant event in Spike gene, which was reinforced by phylogenetic tree incongruence analysis, involving novel Tb1 and AlphaCoVs identified in Eptesicus fuscus (family Vespertilionidae) from the U.S. The putative recombinant region is in the S1 subdomain of the Spike gene, encompassing the potential receptor-binding domain of AlphaCoVs. This study reports the first AlphaCoV genomes in molossids from the Americas and provides new insights into recombination as an important mode of evolution of coronaviruses involved in cross-species transmission. IMPORTANCE This study generated three novel complete AlphaCoV genomes (Tb1, Tb2, and Tb3 isolates) identified in individuals of Tadarida brasiliensis from Argentina, which showed two different evolutionary patterns and are the first to be reported in the family Molossidae in the Americas. The novel Tb1 isolate was found to be involved in a putative recombination event with alphacoronaviruses identified in bats of the genus Eptesicus from the U.S., whereas isolates Tb2 and Tb3 were found in different collection seasons and might be involved in persistent viral infections in the bat colony. These findings contribute to our knowledge of the global diversity of bat coronaviruses in poorly studied species and highlight the different evolutionary aspects of AlphaCoVs circulating in bat populations in Argentina.

Viral Metagenomic Data Analyses of Five New World Bat Species from Argentina: Identification of 35 Novel DNA Viruses

Bats are natural reservoirs of a variety of zoonotic viruses, many of which cause severe human diseases. Characterizing viruses of bats inhabiting different geographical regions is important for understanding their viral diversity and for detecting viral spillovers between animal species. Herein, the diversity of DNA viruses of five arthropodophagous bat species from Argentina was investigated using metagenomics. Fecal samples of 29 individuals from five species (Tadarida brasiliensis, Molossus molossus, Eumops bonariensis, Eumops patagonicus, and Eptesicus diminutus) living at two different geographical locations, were investigated. Enriched viral DNA was sequenced using Illumina MiSeq, and the reads were trimmed and filtered using several bioinformatic approaches. The resulting nucleotide sequences were subjected to viral taxonomic classification. In total, 4,520,370 read pairs were sequestered by sequencing, and 21.1% of them mapped to viral taxa. Circoviridae and Genomoviridae were the most prevalent among vertebrate viral families in all bat species included in this study. Samples from the T. brasiliensis colony exhibited lower viral diversity than samples from other species of New World bats. We characterized 35 complete genome sequences of novel viruses. These findings provide new insights into the global diversity of bat viruses in poorly studied species, contributing to prevention of emerging zoonotic diseases and to conservation policies for endangered species.

Cost-Effective Method to Perform SARS-CoV-2 Variant Surveillance: Detection of Alpha, Gamma, Lambda, Delta, Epsilon, and Zeta in Argentina

SARS-CoV-2 variants with concerning characteristics have emerged since the end of 2020. Surveillance of SARS-CoV-2 variants was performed on a total of 4,851 samples from the capital city and 10 provinces of Argentina, during 51 epidemiological weeks (EWs) that covered the end of the first wave and the ongoing second wave of the COVID-19 pandemic in the country (EW 44/2020 to EW 41/2021). The surveillance strategy was mainly based on Sanger sequencing of a Spike coding region that allows the identification of signature mutations associated with variants. In addition, whole-genome sequences were obtained from 637 samples. The main variants found were Gamma and Lambda, and to a lesser extent, Alpha, Zeta, and Epsilon, and more recently, Delta. Whereas, Gamma dominated in different regions of the country, both Gamma and Lambda prevailed in the most populated area, the metropolitan region of Buenos Aires. The lineages that circulated on the first wave were replaced by emergent variants in a term of a few weeks. At the end of the ongoing second wave, Delta began to be detected, replacing Gamma and Lambda. This scenario is consistent with the Latin American variant landscape, so far characterized by a concurrent increase in Delta circulation and a stabilization in the number of cases. The cost-effective surveillance protocol presented here allowed for a rapid response in a resource-limited setting, added information on the expansion of Lambda in South America, and contributed to the implementation of public health measures to control the disease spread in Argentina.

A Preliminary Study of the Virome of the South American Free-Tailed Bats (Tadarida brasiliensis) and Identification of Two Novel Mammalian Viruses

Bats provide important ecosystem services as pollinators, seed dispersers, and/or insect controllers, but they have also been found harboring different viruses with zoonotic potential. Virome studies in bats distributed in Asia, Africa, Europe, and North America have increased dramatically over the past decade, whereas information on viruses infecting South American species is scarce. We explored the virome of Tadarida brasiliensis, an insectivorous New World bat species inhabiting a maternity colony in Rosario (Argentina), by a metagenomic approach. The analysis of five pooled oral/anal swab samples indicated the presence of 43 different taxonomic viral families infecting a wide range of hosts. By conventional nucleic acid detection techniques and/or bioinformatics approaches, the genomes of two novel viruses were completely covered clustering into the Papillomaviridae (Tadarida brasiliensis papillomavirus type 1, TbraPV1) and Genomoviridae (Tadarida brasiliensis gemykibivirus 1, TbGkyV1) families. TbraPV1 is the first papillomavirus type identified in this host and the prototype of a novel genus. TbGkyV1 is the first genomovirus reported in New World bats and constitutes a new species within the genus Gemykibivirus. Our findings extend the knowledge about oral/anal viromes of a South American bat species and contribute to understand the evolution and genetic diversity of the novel characterized viruses.

Assessing Gammapapillomavirus infections of mucosal epithelia with two broad-spectrum PCR protocols

Background

Human papillomaviruses (HPVs) have been divided into mucosal and cutaneous types according to their primary epithelial tissue tropism. However, recent studies showed the presence of several cutaneous types in mucosal lesions and healthy mucosa from different anatomical sites.

Methods

Here, the HPV prevalence and type-specific distribution were assessed in a variety of mucosal samples from 435 individuals using a combination of two established broad-spectrum primer systems: Gamma-PV PCR and CUT PCR.

Results

Overall HPV prevalence in anal canal swabs, cervical cancer biopsies, genital warts and oral swabs was 85, 47, 62 and 4%, respectively. In anal canal swabs, Alpha-PVs were most frequently found (59%), followed by Gamma- (37%) and Beta-PVs (4%). The prevalence and persistence of HPV infection in the anal canal of 226 individuals were further explored. Overall HPV, Gamma-PVs and multiple HPV infections were significantly higher in men vs. women (p = 0.034, p = 0.027 and p = 0.003, respectively); multiple HPV infections were more common in individuals ≤40 years (p = 0.05), and significantly higher prevalence of Gamma-PVs and multiple HPV infections was observed in HIV-1-positive vs. HIV-1-negative individuals (p = 0.003 and p = 0.04, respectively). Out of 21 patients with follow-up anal swabs, only one persistent infection with the same type (HPV58) was detected.

Conclusions

Our findings suggest that Gamma-PVs (except species Gamma-6) are ubiquitous viruses with dual muco-cutaneous tissue tropism. Anal canal Gamma-PV infections may be associated with sexual behavior and the host immune status. This study expands the knowledge on Gamma-PVs’ tissue tropism, providing valuable data on the characteristics of HPV infection in the anal canal.

High prevalence of Gammapapillomaviruses (Gamma-PVs) in pre-malignant cutaneous lesions of immunocompetent individuals using a new broad-spectrum primer system, and identification of HPV210, a novel Gamma-PV type

Genus Gammapapillomavirus (Gamma-PV) is the most diverse and largest clade within the Papillomaviridae family. A novel set of degenerate primers targeting the E1 gene was designed and further used in combination with the well-known CUT PCR assay to assess HPV prevalence and genus distribution in a variety of cutaneous samples from 448 immunocompetent individuals. General HPV, Gamma-PV and mixed infections prevalence were significantly higher in actinic keratosis with respect to benign and malignant neoplasms, respectively (p = 0.0047, p = 0.0172, p = 0.00001). Gamma-PVs were significantly more common in actinic keratosis biopsies than Beta- and Alpha-PVs (p = 0.002). The full-length genome sequence of a novel putative Gamma-PV type was amplified by 'hanging droplet' long-range PCR and cloned. The novel virus, designated HPV210, clustered within species Gamma-12. This study provides an additional tool enabling detection of HPV infections in skin and adds new insights about possible early roles of Gamma-PVs in the development of cutaneous malignant lesions.

Natural history of human papillomavirus infection of sun-exposed healthy skin of immunocompetent individuals over three climatic seasons and identification of HPV209, a novel betapapillomavirus

We present the first longitudinal study reporting the natural history of human papillomavirus (HPV) infection in sun-exposed skin of healthy individuals living in a geographical area in which solar UV radiation is influenced by the ozone content of the atmosphere. During three climatic seasons, skin swab samples were obtained from 78 healthy individuals and the prevalence of cutaneous HPVs was assessed with broad-spectrum FAP and CUT primers and determined at 54, 45 and 47 % in spring, summer and winter, respectively. Frequencies of mixed HPV infections were significantly higher in spring with respect to summer and winter (P=0.02). Seventy-one different HPV types/putative types were identified. While 62 volunteers were HPV-infected in at least one season, 23 had persistent infections. β-PVs (β-1) were the most prevalent and persistent. Age was associated with both the infection status (P=0.01) and the type of HPV infection (no infection, indeterminate/transient, persistent P=0.02). The molecular/phylogenetic analysis of the newly identified β-PV, officially designated as HPV209, showed that the virus has a typical genomic organization of cutaneous HPVs with five early (E6, E7, E1, E2 and E4) and two late genes (L2 and L1), which clusters to the species β-2. This provides useful data on cutaneous HPV infections in high UV-exposed regions.

Characterization of novel human papillomavirus types 157, 158 and 205 from healthy skin and recombination analysis in genus γ-Papillomavirus

Gammapapillomavirus (γ-PV) is a diverse and rapidly expanding genus, currently consisting of 79 fully characterized human PV (HPV) types. In this study, three novel types, HPV157, HPV158 and HPV205, obtained from healthy sun-exposed skin of two immunocompetent individuals, were amplified by the "Hanging droplet" long PCR technique, cloned, sequenced and characterized. HPV157, HPV158 and HPV205 genomes comprise 7154-bp, 7192-bp and 7298-bp, respectively, and contain four early (E1, E2, E6 and E7) and two late genes (L1 and L2). Phylogenetic analysis of the L1 ORF placed all novel types within the γ-PV genus: HPV157 was classified as a new member of species γ-12 while HPV158 and HPV205 belong to species γ-1. We then explored potential recombination events in genus γ-PV with the RDP4 program in a dataset of 74 viruses (71 HPV types with available full-length genomes and the 3 novel types). Two events, both located in the E1 ORF, met the inclusion criterion (p-values <0.05 with at least four methods) and persisted in different ORF combinations: an inter-species recombination in species γ-8 (major and minor parents: species γ-24 and γ-11, respectively), and an intra-species recombination in species γ-7 (recombinant strain: HPV170; major and minor parents: HPV-109 and HPV-149, respectively). These findings were confirmed by phylogenetic tree incongruence analysis. An additional incongruence was found in members of species γ-9 but it was not detected by the RDP4. This report expands our knowledge of the family Papillomaviridae and provides for the first time in silico evidence of recombination in genus γ-PV.

Evidence for homologous recombination in Chikungunya Virus

Chikungunya Virus (CHIKV), a mosquito-transmitted alphavirus, causes acute fever and joint pain in humans. Recently, endemic CHIKV infection outbreaks have jeopardized public health in wider geographical regions. Here, we analyze the phylogenetic associations of CHIKV and explore the potential recombination events on 152 genomic isolates deposited in GenBank database. The CHIKV genotypes [West African, Asian, East/Central/South African (ECSA)], and a clear division of ECSA clade into three sub-groups (I-II-III), were defined by Bayesian analysis; similar results were obtained using E1 gene sequences. A nucleotide identity-based approach is provided to facilitate CHIKV classification within ECSA clade. Using seven methods to detect recombination, we found a statistically significant event (p-values range: 1.14×10(-7)-4.45×10(-24)) located within the nsP3 coding region. This finding was further confirmed by phylogenetic networks (PHI Test, p=0.004) and phylogenetic tree incongruence analysis. The recombinant strain, KJ679578/India/2011 (ECSA III), derives from viruses of ECSA III and ECSA I. Our study demonstrates that recombination is an additional mechanism of genetic diversity in CHIKV that might assist in the cross-species transmission process.