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.

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.

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.

The COVID-19 pandemic has lately been driven by Omicron. This work aimed to study the dynamics of SARS-CoV-2 Omicron lineages during the third and fourth waves of COVID-19 in Argentina. Molecular surveillance was performed on 3431 samples from Argentina, between EW44/2021 and EW31/2022. Sequencing, phylogenetic and phylodynamic analyses were performed. A differential dynamic between the Omicron waves was found. The third wave was associated with lineage BA.1, characterized by a high number of cases, very fast displacement of Delta, doubling times of 3.3 days and a low level of lineage diversity and clustering. In contrast, the fourth wave was longer but associated with a lower number of cases, initially caused by BA.2, and later by BA.4/BA.5, with doubling times of about 10 days. Several BA.2 and BA.4/BA.5 sublineages and introductions were detected, although very few clusters with a constrained geographical distribution were observed, suggesting limited transmission chains. The differential dynamic could be due to waning immunity and an increase in population gatherings in the BA.1 wave, and a boosted population (for vaccination or recent prior immunity for BA.1 infection) in the wave caused by BA2/BA.4/BA.5, which may have limited the establishment of the new lineages.

The second wave of COVID-19 occurred in South America in early 2021 and was mainly driven by Gamma and Lambda variants. In this study, we aimed to describe the emergence and local genomic diversity of the SARS-CoV-2 Lambda variant in Argentina, from its initial entry into the country until its detection ceased. Molecular surveillance was conducted on 9356 samples from Argentina between October 2020 and April 2022, and sequencing, phylogenetic, and phylogeographic analyses were performed. Our findings revealed that the Lambda variant was first detected in Argentina in January 2021 and steadily increased in frequency until it peaked in April 2021, with continued detection throughout the year. Phylodynamic analyses showed that at least 18 introductions of the Lambda variant into the country occurred, with nine of them having evidence of onward local transmission. The spatial–-temporal reconstruction showed that Argentine clades were associated with Lambda sequences from Latin America and suggested an initial diversification in the Metropolitan Area of Buenos Aires before spreading to other regions in Argentina. Genetic analyses of genome sequences allowed us to describe the mutational patterns of the Argentine Lambda sequences and detect the emergence of rare mutations in an immunocompromised patient. Our study highlights the importance of genomic surveillance in identifying the introduction and geographical distribution of the SARS-CoV-2 Lambda variant, as well as in monitoring the emergence of mutations that could be involved in the evolutionary leaps that characterize variants of concern.