Tracking the genetic diversity of SARS-CoV-2 variants in Nicaragua throughout the COVID-19 Pandemic.

The global circulation of SARS-CoV-2 has been extensively documented, yet the dynamics within Central America, particularly Nicaragua, remain underexplored. This study characterizes the genomic diversity of SARS-CoV-2 in Nicaragua from March 2020 through December 2022, utilizing 1064 genomes obtained via next-generation sequencing. These sequences were selected nationwide and analyzed for variant classification, lineage predominance, and phylogenetic diversity. We employed both Illumina and Oxford Nanopore Technologies for all sequencing procedures. Results indicated a temporal and spatial shift in dominant lineages, initially from B.1 and A.2 in early 2020 to various Omicron subvariants towards the study's end. Significant lineage shifts correlated with changes in COVID-19 positivity rates, underscoring the epidemiological impact of variant dissemination. The comparative analysis with regional data underscored the low diversity of circulating lineages in Nicaragua and their delayed introduction compared to other countries in the Central American region. The study also linked specific viral mutations with hospitalization rates, emphasizing the clinical relevance of genomic surveillance. This research advances the understanding of SARS-CoV-2 evolution in Nicaragua and provide valuable information regarding its genetic diversity for public health officials in Central America. We highlight the critical role of ongoing genomic surveillance in identifying emergent lineages and informing public health strategies.

Introduction of New Dengue Virus Lineages of Multiple Serotypes after COVID-19 Pandemic, Nicaragua, 2022.

Major dengue epidemics throughout Nicaragua's history have been dominated by 1 of 4 dengue virus serotypes (DENV-1-4). To examine serotypes during the dengue epidemic in Nicaragua in 2022, we performed real-time genomic surveillance in-country and documented cocirculation of all 4 serotypes. We observed a shift toward co-dominance of DENV-1 and DENV-4 over previously dominant DENV-2. By analyzing 135 new full-length DENV sequences, we found that introductions underlay the resurgence: DENV-1 clustered with viruses from Ecuador in 2014 rather than those previously seen in Nicaragua; DENV-3, which last circulated locally in 2014, grouped instead with Southeast Asia strains expanding into Florida and Cuba in 2022; and new DENV-4 strains clustered within a South America lineage spreading to Florida in 2022. In contrast, DENV-2 persisted from the formerly dominant Nicaragua clade. We posit that the resurgence emerged from travel after the COVID-19 pandemic and that the resultant intensifying hyperendemicity could affect future dengue immunity and severity.

Phylodynamics of dengue virus 2 in Nicaragua leading up to the 2019 epidemic reveals a role for lineage turnover.

BACKGROUND: Dengue is a mosquito-borne viral disease posing a significant threat to public health. Dengue virus (DENV) evolution is often characterized by lineage turnover, which, along with ecological and immunological factors, has been linked to changes in dengue phenotype affecting epidemic dynamics. Utilizing epidemiologic and virologic data from long-term population-based studies (the Nicaraguan Pediatric Dengue Cohort Study and Nicaraguan Dengue Hospital-based Study), we describe a lineage turnover of DENV serotype 2 (DENV-2) prior to a large dengue epidemic in 2019. Prior to this epidemic, Nicaragua had experienced relatively low levels of DENV transmission from 2014 to 2019, a period dominated by chikungunya in 2014/15 and Zika in 2016. RESULTS: Our phylogenetic analyses confirmed that all Nicaraguan DENV-2 isolates from 2018 to 2019 formed their own clade within the Nicaraguan lineage of the Asian/American genotype. The emergence of the new DENV-2 lineage reflects a replacement of the formerly dominant clade presiding from 2005 to 2009, a lineage turnover marked by several shared derived amino acid substitutions throughout the genome. To elucidate evolutionary drivers of lineage turnover, we performed selection pressure analysis and reconstructed the demographic history of DENV-2. We found evidence of adaptive evolution by natural selection at the codon level as well as in branch formation. CONCLUSIONS: The timing of its emergence, along with a statistical signal of adaptive evolution and distinctive amino acid substitutions, the latest in the NS5 gene, suggest that this lineage may have increased fitness relative to the prior dominant DENV-2 strains. This may have contributed to the intensity of the 2019 DENV-2 epidemic, in addition to previously identified immunological factors associated with pre-existing Zika virus immunity.

Evolution and epidemiologic dynamics of dengue virus in Nicaragua during the emergence of chikungunya and Zika viruses.

Arthropod-borne viruses (arboviruses) comprise a significant and ongoing threat to human health, infecting hundreds of millions annually. Three such arboviruses include circumtropical dengue, Zika, and chikungunya viruses, exhibiting continuous emergence primarily via Aedes mosquito vectors. Nicaragua has experienced endemic dengue virus (DENV) transmission involving multiple serotypes since 1985, with chikungunya virus (CHIKV) reported in 2014-2015, followed by Zika virus (ZIKV) first reported in 2016. In order to identify patterns of genetic variation and selection pressures shaping the evolution of co-circulating DENV serotypes in light of the arrival of CHIKV and ZIKV, we employed whole-genome sequencing on an Illumina MiSeq platform of random-amplified total RNA libraries to characterize 42 DENV low-passage isolates, derived from viremic patients in Nicaragua between 2013 and 2016. Our approach also revealed clinically undetected co-infections with CHIKV. Of the three DENV serotypes (1, 2, and 3) co-circulating during our study, we uncovered distinct patterns of evolution using comparative phylogenetic inference. DENV-1 genetic variation was structured into two distinct co-circulating lineages with no evidence of positive selection in the origins of either lineage, suggesting they are equally fit. In contrast, the evolutionary history of DENV-2 was marked by positive selection, and a unique, divergent lineage correlated with high epidemic potential emerged in 2015 to drive an outbreak in 2016. DENV-3 genetic variation remained unstructured into lineages throughout the period of study. Thus, this study reveals insights into evolutionary and epidemiologic trends exhibited during the circulation of multiple arboviruses in Nicaragua.

Understanding Past Population Dynamics: Bayesian Coalescent-Based Modeling with Covariates.

Effective population size characterizes the genetic variability in a population and is a parameter of paramount importance in population genetics and evolutionary biology. Kingman's coalescent process enables inference of past population dynamics directly from molecular sequence data, and researchers have developed a number of flexible coalescent-based models for Bayesian nonparametric estimation of the effective population size as a function of time. Major goals of demographic reconstruction include identifying driving factors of effective population size, and understanding the association between the effective population size and such factors. Building upon Bayesian nonparametric coalescent-based approaches, we introduce a flexible framework that incorporates time-varying covariates that exploit Gaussian Markov random fields to achieve temporal smoothing of effective population size trajectories. To approximate the posterior distribution, we adapt efficient Markov chain Monte Carlo algorithms designed for highly structured Gaussian models. Incorporating covariates into the demographic inference framework enables the modeling of associations between the effective population size and covariates while accounting for uncertainty in population histories. Furthermore, it can lead to more precise estimates of population dynamics. We apply our model to four examples. We reconstruct the demographic history of raccoon rabies in North America and find a significant association with the spatiotemporal spread of the outbreak. Next, we examine the effective population size trajectory of the DENV-4 virus in Puerto Rico along with viral isolate count data and find similar cyclic patterns. We compare the population history of the HIV-1 CRF02_AG clade in Cameroon with HIV incidence and prevalence data and find that the effective population size is more reflective of incidence rate. Finally, we explore the hypothesis that the population dynamics of musk ox during the Late Quaternary period were related to climate change. [Coalescent; effective population size; Gaussian Markov random fields; phylodynamics; phylogenetics; population genetics.

Insights into the molecular evolution of Dengue virus type 4 in Puerto Rico over two decades of emergence.

Dengue has emerged globally as a major human health problem since the 1950s and is now the most important arboviral disease of humans, infecting nearly 400 million people annually. While some cases are asymptomatic, others can develop a febrile illness (dengue fever) or even progress to severe and fatal dengue. Dengue is caused by any of 4 closely related but distinct viruses, known as Dengue virus serotype 1 to 4 (DENV-1 to DENV-4) which are maintained in endemic transmission to humans in large urban centers of the tropics by Aedes mosquitoes. Since the early 1960s, Puerto Rico, a major metropolitan center in the Caribbean, has experienced increasingly larger and clinically more severe epidemics following the introduction of all four dengue serotypes. The first dengue hemorrhagic fever epidemic in 1986, and a particularly severe outbreak in 1998 were dominated by novel DENV-4 strains that evolved in Puerto Rico, replacing earlier strains and spreading throughout the region. Sequence characterization of 54 complete DENV-4 genomes and their comparative evolution against 74 previously published viral sequences from the region over several decades shows that DENV-4 strains from these periods were genetically distinct based on unique changes in the envelope and non-structural genes. Their replacement of earlier strains in Puerto Rico progressed rapidly, suggesting that strong natural selection played a role in their fixation. This study confirms that DENVs evolve through rapid lineage turnover driven in part by natural selection and genetic drift.

Phylogeography and population dynamics of dengue viruses in the Americas.

Changes in Dengue virus (DENV) disease patterns in the Americas over recent decades have been attributed, at least in part, to repeated introduction of DENV strains from other regions, resulting in a shift from hypoendemicity to hyperendemicity. Using newly sequenced DENV-1 and DENV-3 envelope (E) gene isolates from 11 Caribbean countries, along with sequences available on GenBank, we sought to document the population genetic and spatiotemporal transmission histories of the four main invading DENV genotypes within the Americas and investigate factors that influence the rate and intensity of DENV transmission. For all genotypes, there was an initial invasion phase characterized by rapid increases in genetic diversity, which coincided with the first confirmed cases of each genotype in the region. Rapid geographic dispersal occurred upon each genotype's introduction, after which individual lineages were locally maintained, and gene flow was primarily observed among neighboring and nearby countries. There were, however, centers of viral diversity (Barbados, Puerto Rico, Colombia, Suriname, Venezuela, and Brazil) that were repeatedly involved in gene flow with more distant locations. For DENV-1 and DENV-2, we found that a "distance-informed" model, which posits that the intensity of virus movement between locations is inversely proportional to the distance between them, provided a better fit than a model assuming equal rates of movement between all pairs of countries. However, for DENV-3 and DENV-4, the more stochastic "equal rates" model was preferred.

Isolation of a new strain of the flavivirus cell fusing agent virus in a natural mosquito population from Puerto Rico.

The genus Flavivirus contains approximately 70 single-stranded, positive-sense RNA viruses that are mosquito-borne, tick-borne or have no known vector. Two discoveries support previous suggestions of the existence of a large number of unsampled flaviviruses: (i) a new flavivirus, Kamiti River virus, was recently isolated from Kenyan mosquitoes, and (ii) sequences with high similarity to those of flaviviruses have been found integrated into the genome of Aedes mosquitoes, suggesting a past infection with a virus (or viruses) that has yet to be discovered. These sequences were related most closely to a flavivirus that infects insects alone, cell fusing agent virus (CFAV). CFAV was originally isolated in the laboratory from an Aedes aegypti cell line. To date, this virus had not been found in the wild. In the present study, over 40 isolates of a novel strain of CFAV were discovered from mature mosquitoes sampled from the wild in Puerto Rico. The viral strain was present in a range of mosquito species, including Aedes aegypti, Aedes albopictus and Culex sp., from numerous locations across the island and, importantly, in mosquitoes of both sexes, suggesting vertical transmission. Here, results from viral screening, and cell culture and molecular identification of the infected mosquitoes are presented. Experimental-infection tests were also conducted by using the original CFAV strain and a highly efficient reverse-transcription mechanism has been documented, in which initiation of copying occurs at the 3' terminus of either the genomic RNA or the intermediate of replication, potentially elucidating the mechanism by which flaviviral sequences may have integrated into mosquito genomes.

Molecular evolution of dengue 2 virus in Puerto Rico: positive selection in the viral envelope accompanies clade reintroduction.

Dengue virus is a circumtropical, mosquito-borne flavivirus that infects 50-100 million people each year and is expanding in both range and prevalence. Of the four co-circulating viral serotypes (DENV-1 to DENV-4) that cause mild to severe febrile disease, DENV-2 has been implicated in the onset of dengue haemorrhagic fever (DHF) in the Americas in the early 1980s. To identify patterns of genetic change since DENV-2's reintroduction into the region, molecular evolution in DENV-2 from Puerto Rico (PR) and surrounding countries was examined over a 20 year period of fluctuating disease incidence. Structural genes (over 20 % of the viral genome), which affect viral packaging, host-cell entry and immune response, were sequenced for 91 DENV-2 isolates derived from both low- and high-prevalence years. Phylogenetic analyses indicated that DENV-2 outbreaks in PR have been caused by viruses assigned to subtype IIIb, originally from Asia. Variation amongst DENV-2 viruses in PR has since largely arisen in situ, except for a lineage-replacement event in 1994 that appears to have non-PR New World origins. Although most structural genes have remained relatively conserved since the 1980s, strong evidence was found for positive selection acting on a number of amino acid sites in the envelope gene, which have also been important in defining phylogenetic structure. Some of these changes are exhibited by the multiple lineages present in 1994, during the largest Puerto Rican outbreak of dengue, suggesting that they may have altered disease dynamics, although their functional significance will require further investigation.

Invasion and maintenance of dengue virus type 2 and type 4 in the Americas

Dengue virus type 4 (DENV-4) was first reported in the Americas in 1981, where it caused epidemics of dengue fever throughout the region. In the same year, the region’s first epidemic of dengue hemorrhagic fever was reported, caused by an Asian strain of dengue virus type 2 (DENV-2) that was distinct from the American subtype circulating previously. Despite the importance of these epidemics, little is known about the rates or determinants of viral spread among island and mainland populations or their directions of movement. We employed a Bayesian coalescent approach to investigate the transmission histories of DENV-2 and DENV-4 since their introduction in 1981 and a parsimony method to assess patterns of strain migration. For both viruses there was an initial invasion phase characterized by an exponential increase in the number of DENV lineages, after which levels of genetic diversity remained constant despite reported fluctuations in DENV-2 and DENV-4 activity. Strikingly, viral lineage numbers increased far more rapidly for DENV-4 than DENV-2, indicative of a more rapid rate of exponential population growth in DENV-4 or a higher rate of geographic dispersal, allowing this virus to move more effectively among localities. We propose that these contrasting dynamics may reflect underlying differences in patterns of host immunity. Despite continued gene flow along particular transmission routes, the overall extent of viral traffic was less than expected under panmixis. Hence, DENV in the Americas has a clear geographic structure that maintains viral diversity between outbreaks.

Phylogeography and molecular evolution of dengue 2 in the Caribbean basin, 1981-2000.

We sequenced the envelope (E) genes of 59 DEN-2 isolates collected from ten Caribbean islands, six South American countries, and two Central American countries between 1981 and 2000, a period characterized by hyperendemicity and increased incidence of severe dengue. Fifty-two isolates belonged to "American/Asian" subtype IIIb, possessing a characteristic polar residue at envelope aa position 390 (N [n = 48] or S [n = 4]) common to that group. Six isolates from Trinidad (1981), Honduras (1991 [4]), and El Salvador (1987) fell into the "Native American" subtype V (D at aa 390), and one from Honduras (1986) belonged to "Asian" subtype I. The data suggest that after its first isolation in the Caribbean in 1981, genotype IIIb spread throughout the Americas and effectively replaced subtype V throughout the Caribbean basin. The strain also evolved into several distinct lineages, based on substitutions in the E glycoprotein (amino acids 91 and 131), two of which were still in circulation in 2000. Interestingly, a molecular clock did not fit the data well, suggesting that other sources of rate variation, such as differential selection or differences in effective population sizes, may exist among lineages. Our results indicate the importance of large temporal- and geographical-scale phylogenetic studies in understanding disease dynamics, particularly where replacements between regions can occur.

Phylogeography and molecular evolution of dengue 2 in the Caribbean basin, 1981–2000

We sequenced the envelope (E) genes of 59 DEN-2 isolates collected from ten Caribbean islands, six South American countries, and two Central American countries between 1981 and 2000, a period characterized by hyperendemicity and increased incidence of severe dengue. Fifty-two isolates belonged to “American/Asian” subtype IIIb, possessing a characteristic polar residue at envelope aa position 390 (N [n = 48] or S [n = 4]) common to that group. Six isolates from Trinidad (1981), Honduras (1991 [4]), and El Salvador (1987) fell into the “Native American” subtype V (D at aa 390), and one from Honduras (1986) belonged to “Asian” subtype I. The data suggest that after its first isolation in the Caribbean in 1981, genotype IIIb spread throughout the Americas and effectively replaced subtype V throughout the Caribbean basin. The strain also evolved into several distinct lineages, based on substitutions in the E glycoprotein (amino acids 91 and 131), two of which were still in circulation in 2000. Interestingly, a molecular clock did not fit the data well, suggesting that other sources of rate variation, such as differential selection or differences in effective population sizes, may exist among lineages. Our results indicate the importance of large temporal- and geographical-scale phylogenetic studies in understanding disease dynamics, particularly where replacements between regions can occur.

Selection-driven evolution of emergent dengue virus.

In the last four decades the incidence of dengue fever has increased 30-fold worldwide, and over half the world's population is now threatened with infection from one or more of four co-circulating viral serotypes (DEN-1 through DEN-4). To determine the role of viral molecular evolution in emergent disease dynamics, we sequenced 40% of the genome of 82 DEN-4 isolates collected from Puerto Rico over the 20 years since the onset of endemic dengue on the island. Isolates were derived from years with varying levels of DEN-4 prevalence. Over our sampling period there were marked evolutionary shifts in DEN-4 viral populations circulating in Puerto Rico; viral lineages were temporally clustered and the most common genotype at a particular sampling time often arose from a previously rare lineage. Expressed changes in structural genes did not appear to drive this lineage turnover, even though these regions include primary determinants of viral antigenic properties. Instead, recent dengue evolution can be attributed in part to positive selection on the nonstructural gene 2A (NS2A), whose functions may include replication efficiency and antigenicity. During the latest and most severe DEN-4 epidemic in Puerto Rico, in 1998, viruses were distinguished by three amino acid changes in NS2A that were fixed far faster than expected by drift alone. Our study therefore demonstrates viral genetic turnover within a focal population and the potential importance of adaptive evolution in viral epidemic expansion.

Molecular evolution and phylogeny of dengue type 4 virus in the Caribbean

We sequenced the E gene and adjacent prM/M and NS1 junctions (1940 bp) of 48 Dengue-4 (DEN-4) isolates collected between 1981 and 1999 from 8 Caribbean islands and from 7 South and Central American countries. Phylogenetic analysis confirms a single introduction in the early 1980s and a high degree of gene flow resulting in a pattern of evolution defined more by time period than geographic origin, especially within the Caribbean basin. A modern Caribbean clade consisting of four distinct lineages has arisen, comprised of isolates from Caribbean islands and nearby regions of South America. This clade is defined by three amino acid substitutions in the E (aa 163 and 351) and NS1 (aa 52) proteins. These findings highlight the importance of migration and gene flow in dengue viral change and suggest that efforts to understand disease dynamics in the Caribbean basin need to focus at regional, rather than local scales.