Epidemiological, virological and clinical characterization of a Dengue/Zika outbreak in the Caribbean region of Costa Rica 2017-2018.

The increase in incidence and geographical expansion of viruses transmitted by the Aedes mosquitoes, such as dengue (DENV) and zika (ZIKV) in the Americas, represents a burden for healthcare systems in tropical and subtropical regions. These and other under-detected arboviruses co-circulate in Costa Rica, adding additional complexity to their management due to their shared epidemiological behavior and similarity of symptoms in early stages. Since diagnostics of febrile illness is mostly based on clinical symptoms alone, we gathered acute-phase serum and urine from 399 samples of acute dengue-like cases from two healthcare facilities of Costa Rica, during an outbreak of arboviruses from July 2017 to May 2018, and tested them using molecular and serological methods. The analyses showed that of the clinically presumptive arbovirus cases that were reported, only 39.4% (n=153) of the samples were confirmed positive by RT-PCR to be DENV (DENV (10.3%), CHIKV (0.2%), ZIKV (27.3%), or mixed infections (1.5%). RT-PCR for other alphaviruses and flaviviruses, and PCR for Leptospira sp were negative. Furthermore, to assess flavivirus positivity in post-acute patients, the negative sera were tested against Dengue-IgM. 20% of sera were found positive, confounding even more the definitive number of cases, and emphasizing the need of several distinct diagnostic tools for accurate diagnostics. Molecular characterization of the prM and E genes from isolated viruses revealed that the American/Asian genotype of DENV-2 and the Asian lineage of ZIKV were circulating during this outbreak. Two different clades of DENV-2 American/Asian genotype were identified to co-circulate in the same region and a difference in the platelet and leukocyte count was noted between people infected with each clade, suggesting a putative distinct virulence. Our study sheds light on the necessity for healthcare strategies in managing arbovirus outbreaks, emphasizing the importance of comprehensive molecular and serological diagnostic approaches, as well as molecular characterization. This approach aids in enhancing our understanding of the clinical and epidemiological aspects of arboviral diseases during outbreaks. Our research highlights the need to strengthen training programs for health professionals and the need to increase research-based on laboratory evidence for diagnostic accuracy, guidance, development and implementation of public health interventions and epidemiological surveillance.

Occurrence and risk assessment of pharmaceuticals in hospital wastewater in Costa Rica.

This work aims to determine the occurrence, hazard and prioritization of pharmaceuticals from hospital wastewater in Costa Rica through the monitoring of 70 compounds and assessing their environmental risk through a hazard quotient approach (HQ). Moreover, the quantification of selected antibiotic resistance genes (ARGs) was conducted for the first time in this matrix in this geographical location. Thirty-four pharmaceuticals were detected, being caffeine, 1,7-dimethylxanthine, acetaminophen, ibuprofen, naproxen, ciprofloxacin and ketoprofen the most frequent (>50% of the samples). Eighteen pharmaceuticals exhibited high hazard (HQ ≥ 1), while five more showed medium hazard (1 > HQ ≥ 0.1). Prioritization, which also included frequency parameters, revealed caffeine, lovastatin, diphenhydramine, acetaminophen, ibuprofen, ciprofloxacin, and sildenafil as the compounds of major concern. Similarly, cumulative hazard per sample (ΣHQ) estimated high hazard towards aquatic organisms in every sample. All selected ARGs, except mcr-1 (polymyxin resistance), were detected. Among genes conferring resistance to beta-lactams, blaCTX-M and blaKPC were the most abundant, related to resistance to cephalosporins and carbapenems. Ecotoxicological evaluation showed mostly low toxicity towards Daphnia magna and Vibrio fischeri, contrary to the marked effect observed towards Lactuca sativa. These findings provide relevant and novel information on the risk posed by hospital wastewater and their pharmaceutical content in the Latin American environmental context.

Experimental infection of Artibeus lituratus bats and no detection of Zika virus in neotropical bats from French Guiana, Peru, and Costa Rica suggests a limited role of bats in Zika transmission.

Bats are important natural reservoir hosts of a diverse range of viruses that can be transmitted to humans and have been suggested to play an important role in the Zika virus (ZIKV) transmission cycle. However, the exact role of these animals as reservoirs for flaviviruses is still controversial. To further expand our understanding of the role of bats in the ZIKV transmission cycle in Latin America, we carried out an experimental infection in wild-caught Artibeus lituratus bats and sampled several free-living neotropical bats across three countries of the region. Experimental ZIKV infection was performed in wild-caught adult bats (4 females and 5 males). The most relevant findings were hemorrhages in the bladder, stomach and patagium. Significant histological findings included inflammatory infiltrate consisting of a predominance of neutrophils and lymphocytes, in addition to degeneration in the reproductive tract of males and females. This suggests that bat reproduction might be at some level affected by ZIKV. Leukopenia was also observed in some inoculated animals. Hemorrhages, genital alterations, and leukopenia are suggested to be caused by ZIKV; however, since these were wild-caught bats, we cannot exclude other agents. Detection of ZIKV by qPCR was observed at low concentrations in only two urine samples in two inoculated animals. All other animals and tissues tested were negative. Finally, no virus-neutralizing antibodies were found in any animal. To determine ZIKV infection in nature, the blood of a total of 2056 bats was sampled for ZIKV detection by qPCR. Most of the sampled individuals belonged to the genus Pteronotus sp. (23%), followed by the species Carollia sp. (17%), Anoura sp. (14%), and Molossus sp. (13.7%). No sample of any tested species was positive for ZIKV by qPCR. These results together suggest that bats are not efficient amplifiers or reservoirs of ZIKV and may not have an important role in ZIKV transmission dynamics.

Overview of the SARS-CoV-2 genotypes circulating in Latin America during 2021.

Latin America is one of the regions in which the COVID-19 pandemic has a stronger impact, with more than 72 million reported infections and 1.6 million deaths until June 2022. Since this region is ecologically diverse and is affected by enormous social inequalities, efforts to identify genomic patterns of the circulating SARS-CoV-2 genotypes are necessary for the suitable management of the pandemic. To contribute to the genomic surveillance of the SARS-CoV-2 in Latin America, we extended the number of SARS-CoV-2 genomes available from the region by sequencing and analyzing the viral genome from COVID-19 patients from seven countries (Argentina, Brazil, Costa Rica, Colombia, Mexico, Bolivia, and Peru). Subsequently, we analyzed the genomes circulating mainly during 2021 including records from GISAID database from Latin America. A total of 1,534 genome sequences were generated from seven countries, demonstrating the laboratory and bioinformatics capabilities for genomic surveillance of pathogens that have been developed locally. For Latin America, patterns regarding several variants associated with multiple re-introductions, a relatively low percentage of sequenced samples, as well as an increment in the mutation frequency since the beginning of the pandemic, are in line with worldwide data. Besides, some variants of concern (VOC) and variants of interest (VOI) such as Gamma, Mu and Lambda, and at least 83 other lineages have predominated locally with a country-specific enrichments. This work has contributed to the understanding of the dynamics of the pandemic in Latin America as part of the local and international efforts to achieve timely genomic surveillance of SARS-CoV-2.

Seasonal changes in the diversity, host preferences and infectivity of mosquitoes in two arbovirus-endemic regions of Costa Rica.

BACKGROUND: Mosquitoes are vectors of various arboviruses belonging to the genera Alphavirus and Flavivirus, and Costa Rica is endemic to several of them. The aim of this study was to describe and analyze the community structure of such vectors in Costa Rica. METHODS: Sampling was performed in two different coastal locations of Costa Rica with evidence of arboviral activity during rainy and dry seasons. Encephalitis vector surveillance traps, CDC female gravid traps and ovitraps were used. Detection of several arboviruses by Pan-Alpha and Pan-Flavi PCR was attempted. Blood meals were also identified. The Normalized Difference Vegetation Index (NDVI) was estimated for each area during the rainy and dry seasons. The Chao2 values for abundance and Shannon index for species diversity were also estimated. RESULTS: A total of 1802 adult mosquitoes belonging to 55 species were captured, among which Culex quinquefasciatus was the most caught species. The differences in NDVI were higher between seasons and between regions, yielding lower Chao-Sørensen similarity index values. Venezuelan equine encephalitis virus, West Nile virus and Madariaga virus were not detected at all, and dengue virus and Zika virus were detected in two separate Cx. quinquefasciatus specimens. The primary blood-meal sources were chickens (60%) and humans (27.5%). Both sampled areas were found to have different seasonal dynamics and population turnover, as reflected in the Chao2 species richness estimation values and Shannon diversity index. CONCLUSION: Seasonal patterns in mosquito community dynamics in coastal areas of Costa Rica have strong differences despite a geographical proximity. The NDVI influences mosquito diversity at the regional scale more than at the local scale. However, year-long continuous sampling is required to better understand local dynamics.

Clinical Profiles at the Time of Diagnosis of SARS-CoV-2 Infection in Costa Rica During the Pre-vaccination Period Using a Machine Learning Approach.

The clinical manifestations of COVID-19, caused by the SARS-CoV-2, define a large spectrum of symptoms that are mainly dependent on the human host conditions. In Costa Rica, more than 169,000 cases and 2185 deaths were reported during the year 2020, the pre-vaccination period. To describe the clinical presentations at the time of diagnosis of SARS-CoV-2 infection in Costa Rica during the pre-vaccination period, we implemented a symptom-based clustering using machine learning to identify clusters or clinical profiles at the population level among 18,974 records of positive cases. Profiles were compared based on symptoms, risk factors, viral load, and genomic features of the SARS-CoV-2 sequence. A total of 18 symptoms at time of diagnosis of SARS-CoV-2 infection were reported with a frequency > 1%, and those were used to identify seven clinical profiles with a specific composition of clinical manifestations. In the comparison between clusters, a lower viral load was found for the asymptomatic group, while the risk factors and the SARS-CoV-2 genomic features were distributed among all the clusters. No other distribution patterns were found for age, sex, vital status, and hospitalization. In conclusion, during the pre-vaccination time in Costa Rica, the symptoms at the time of diagnosis of SARS-CoV-2 infection were described in clinical profiles. The host co-morbidities and the SARS-CoV-2 genotypes are not specific of a particular profile, rather they are present in all the groups, including asymptomatic cases. In addition, this information can be used for decision-making by the local healthcare institutions (first point of contact with health professionals, case definition, or infrastructure). In further analyses, these results will be compared against the profiles of cases during the vaccination period. Supplementary Information: The online version contains supplementary material available at 10.1007/s43657-022-00058-x.

Serological Positivity against Selected Flaviviruses and Alphaviruses in Free-Ranging Bats and Birds from Costa Rica Evidence Exposure to Arboviruses Seldom Reported Locally in Humans.

Arboviruses have two ecological transmission cycles: sylvatic and urban. For some, the sylvatic cycle has not been thoroughly described in America. To study the role of wildlife in a putative sylvatic cycle, we sampled free-ranging bats and birds in two arbovirus endemic locations and analyzed them using molecular, serological, and histological methods. No current infection was detected, and no significant arbovirus-associated histological changes were observed. Neutralizing antibodies were detected against selected arboviruses. In bats, positivity in 34.95% for DENV-1, 16.26% for DENV-2, 5.69% for DENV-3, 4.87% for DENV-4, 2.43% for WNV, 4.87% for SLEV, 0.81% for YFV, 7.31% for EEEV, and 0.81% for VEEV was found. Antibodies against ZIKV were not detected. In birds, PRNT results were positive against WNV in 0.80%, SLEV in 5.64%, EEEV in 8.4%, and VEEV in 5.63%. An additional retrospective PRNT analysis was performed using bat samples from three additional DENV endemic sites resulting in a 3.27% prevalence for WNV and 1.63% for SLEV. Interestingly, one sample resulted unequivocally WNV positive confirmed by serum titration. These results suggest that free-ranging bats and birds are exposed to not currently reported hyperendemic-human infecting Flavivirus and Alphavirus; however, their role as reservoirs or hosts is still undetermined.

Serosurvey in Two Dengue Hyperendemic Areas of Costa Rica Evidence Active Circulation of WNV and SLEV in Peri-Domestic and Domestic Animals and in Humans.

Costa Rica harbors several flaviviruses, including Dengue (DENV), Zika (ZIKV), West Nile virus (WNV), and Saint Louis encephalitis virus (SLEV). While DENV and ZIKV are hyperendemic, previous research indicates restricted circulation of SLEV and WNV in animals. SLEV and WNV seroprevalence and high transmission areas have not yet been measured. To determine the extents of putative WNV and SLEV circulation, we sampled peri-domestic and domestic animals, humans, and mosquitoes in rural households located in two DENV and ZIKV hyperendemic regions during the rainy and dry seasons of 2017-2018 and conducted plaque reduction neutralization test assay for serology (PRNT) and RT-PCR for virus detection. In Cuajiniquil, serological evidence of WNV and SLEV was found in equines, humans, chickens, and wild birds. Additionally, five seroconversion events were recorded for WNV (2 equines), SLEV (1 human), and DENV-1 (2 humans). In Talamanca, WNV was not found, but serological evidence of SLEV circulation was recorded in equines, humans, and wild birds. Even though no active viral infection was detected, the seroconversion events recorded here indicate recent circulation of SLEV and WNV in these two regions. This study thus provides clear-cut evidence for WNV and SLEV presence in these areas, and therefore, they should be considered in arboviruses differential diagnostics and future infection prevention campaigns.

SARS-CoV-2 genomic surveillance in Costa Rica: Evidence of a divergent population and an increased detection of a spike T1117I mutation.

Genome sequencing is a key strategy in the surveillance of SARS-CoV-2, the virus responsible for the COVID-19 pandemic. Latin America is the hardest-hit region of the world, accumulating almost 20% of COVID-19 cases worldwide. In Costa Rica, from the first detected case on March 6th to December 31st almost 170,000 cases have been reported. We analyzed the genomic variability during the SARS-CoV-2 pandemic in Costa Rica using 185 sequences, 52 from the first months of the pandemic, and 133 from the current wave. Three GISAID clades (G, GH, and GR) and three PANGOLIN lineages (B.1, B.1.1, and B.1.291) were predominant, suggesting multiple re-introductions from other regions. The whole-genome variant calling analysis identified a total of 283 distinct nucleotide variants, following a power-law distribution with 190 single nucleotide mutations in a single sequence, and only 16 mutations were found in >5% sequences. These mutations were distributed through the whole genome. The prevalence of worldwide-found variant D614G in the Spike (98.9% in Costa Rica), ORF8 L84S (1.1%) is similar to what is found elsewhere. Interestingly, the frequency of mutation T1117I in the Spike has increased during the current pandemic wave beginning in May 2020 in Costa Rica, reaching 29.2% detection in the full genome analyses in November 2020. This variant has been observed in less than 1% of the GISAID reported sequences worldwide in 2020. Structural modeling of the Spike protein with the T1117I mutation suggests a potential effect on the viral oligomerization needed for cell infection, but no differences with other genomes on transmissibility, severity nor vaccine effectiveness are predicted. In conclusion, genome analyses of the SARS-CoV-2 sequences over the course of the COVID-19 pandemic in Costa Rica suggest the introduction of lineages from other countries and the detection of mutations in line with other studies, but pointing out the local increase in the detection of Spike-T1117I variant. The genomic features of this virus need to be monitored and studied in further analyses as part of the surveillance program during the pandemic.

Impaired pesticide removal and detoxification by biomixtures during the simulated pesticide application cycle of a tropical agricultural system.

Biopurification systems (BPS) or biobeds have been developed to attenuate point-source contamination due to inappropriate pesticide handling or disposal of agricultural wastewaters. The biomixture used for this strategy should be able to remove different active ingredients but its efficiency can vary due to the constant load of pesticides from crop application programs. For that reason, the performance of biomixtures in conditions that mimic the real pesticide treatment before their implementation in field settings should be assayed. This study aimed to evaluate the removal and detoxifying capacity of a previously formulated biomixture (coconut fiber, 50% v/v; compost, 25%; and soil pre-exposed to pesticides, 25%) during a simulated cycle of pesticide application (93 days) for potato production. The scheme included a first application of linuron followed by a weekly alternated treatment of the mixtures chlorpyrifos/metalaxyl and malathion/dimethomorph, and antibiotics at day 72. The biomixture showed efficient removal of linuron (half-life <15 days), and a fluctuating transformation rate for the other compounds. A constant and sustained removal was observed for malathion and methalaxyl. In contrast, lower efficiency and accumulation was described for chlorpyrifos and dimethomorph. Following antibiotic treatment, changes on pesticide removal were observed only in the case of chlorpyrifos, whose removal was slightly enhanced. Furthermore, acute toxicity assays showed limited detoxification of the matrix, especially when compounds began to accumulate. Summarizing, our experiments showed that the proposed biomixture does not support a proper removal of the pesticides during the simulated application cycle of potato production. Further optimization of a biopurification system is required to guarantee the successful elimination of pesticide combinations when applied in field conditions.

Neotropical bats that co-habit with humans function as dead-end hosts for dengue virus.

Several studies have shown Dengue Virus (DENV) nucleic acids and/or antibodies present in Neotropical wildlife including bats, suggesting that some bat species may be susceptible to DENV infection. Here we aim to elucidate the role of house-roosting bats in the DENV transmission cycle. Bats were sampled in households located in high and low dengue incidence regions during rainy and dry seasons in Costa Rica. We captured 318 bats from 12 different species in 29 households. Necropsies were performed in 205 bats to analyze virus presence in heart, lung, spleen, liver, intestine, kidney, and brain tissue. Histopathology studies from all organs showed no significant findings of disease or infection. Sera were analyzed by PRNT90 for a seroprevalence of 21.2% (51/241), and by PCR for 8.8% (28/318) positive bats for DENV RNA. From these 28 bats, 11 intestine samples were analyzed by RT-PCR. Two intestines were DENV RNA positive for the same dengue serotype detected in blood. Viral isolation from all positive organs or blood was unsuccessful. Additionally, viral load analyses in positive blood samples by qRT-PCR showed virus concentrations under the minimal dose required for mosquito infection. Simultaneously, 651 mosquitoes were collected using EVS-CO2 traps and analyzed for DENV and feeding preferences (bat cytochrome b). Only three mosquitoes were found DENV positive and none was positive for bat cytochrome b. Our results suggest an accidental presence of DENV in bats probably caused from oral ingestion of infected mosquitoes. Phylogenetic analyses suggest also a spillover event from humans to bats. Therefore, we conclude that bats in these urban environments do not sustain DENV amplification, they do not have a role as reservoirs, but function as epidemiological dead end hosts for this virus.

Neotropical primary bat cell lines show restricted dengue virus replication.

Dengue is the most widespread arboviral disease affecting humans. Bats are recognized carriers of emerging viral zoonoses and have been proposed as dengue reservoirs, since RNA/NS1 and/or antiviral antibodies have been detected. Yet, experimental inoculation of Artibeus bats failed to show virus replication. This conflicting results prevent drawing further conclusions of whether bats sustain dengue infection. To test bat cellular permissivity to dengue infection, we established primary bat embryonic cells from diverse organs and tissues of Artibeus jamaicensis, Molossus sinaloae, and Desmodus rotundus. We observed a limited serotype-, organ-, and bat species- specific dengue susceptibility. Only some Molossus-derived primary cells sustained poorly initial Dengue serotype-1 replication, though it was latter absent. To elucidate if Molossus bats may play a role in dengue replication, ecological or in vivo experiments must be performed. Taken together our results show that Dengue did not replicate efficiently in cell lines derived from Neotropical bat species.

Molecular Characterization of Two Major Dengue Outbreaks in Costa Rica.

Dengue virus (DENV) (Flavivirus, Flaviviridae) is a reemerging arthropod-borne virus with a worldwide circulation, transmitted mainly by Aedes aegypti and Aedes albopictus mosquitoes. Since the first detection of its main transmitting vector in 1992 and the invasion of DENV-1 in 1993, Costa Rica has faced dengue outbreaks yearly. In 2007 and 2013, Costa Rica experienced two of the largest outbreaks in terms of total and severe cases. To provide genetic information about the etiologic agents producing these outbreaks, we conducted phylogenetic analysis of viruses isolated from human samples. A total of 23 DENV-1 and DENV-2 sequences were characterized. These analyses signaled that DENV-1 genotype V and DENV-2 American/Asian genotype were circulating in those outbreaks. Our results suggest that the 2007 and 2013 outbreak viral strains of DENV-1 and DENV-2 originated from nearby countries and underwent in situ microevolution.