H19 influenza A virus exhibits species-specific MHC class II receptor usage.

Avian influenza A virus (IAV) surveillance in Northern California, USA, revealed unique IAV hemagglutinin (HA) genome sequences in cloacal swabs from lesser scaups. We found two closely related HA sequences in the same duck species in 2010 and 2013. Phylogenetic analyses suggest that both sequences belong to the recently discovered H19 subtype, which thus far has remained uncharacterized. We demonstrate that H19 does not bind the canonical IAV receptor sialic acid (Sia). Instead, H19 binds to the major histocompatibility complex class II (MHC class II), which facilitates viral entry. Unlike the broad MHC class II specificity of H17 and H18 from bat IAV, H19 exhibits a species-specific MHC class II usage that suggests a limited host range and zoonotic potential. Using cell lines overexpressing MHC class II, we rescued recombinant H19 IAV. We solved the H19 crystal structure and identified residues within the putative Sia receptor binding site (RBS) that impede Sia-dependent entry.

Influence of Sex on Respiratory Syncytial Virus Genotype Infection Frequency and Nasopharyngeal Microbiome.

Respiratory syncytial virus (RSV) has a significant health burden in children, older adults, and the immunocompromised. However, limited effort has been made to identify emergence of new RSV genotypes' frequency of infection and how the combination of nasopharyngeal microbiome and viral genotypes impact RSV disease outcomes. In an observational cohort designed to capture the first infant RSV infection, we employed multi-omics approaches to sequence 349 RSV complete genomes and matched nasopharyngeal microbiomes, during which the 2012/2013 season was dominated by RSV-A, whereas 2013 and 2014 was dominated by RSV-B. We found non-G-72nt-duplicated RSV-A strains were more frequent in male infants (P = 0.02), whereas G-72nt-duplicated genotypes (which is ON1 lineage) were seen equally in both males and females. DESeq2 testing of the nasal microbiome showed Haemophilus was significantly more abundant in infants with RSV-A infection compared to infants with RSV-B infection (adjusted P = 0.002). In addition, the broad microbial clustering of the abundant genera was significantly associated with infant sex (P = 0.03). Overall, we show sex differences in infection by RSV genotype and host nasopharyngeal microbiome, suggesting an interaction between host genetics, virus genotype, and associated nasopharyngeal microbiome. IMPORTANCE Respiratory syncytial virus (RSV) is one of the leading causes of lower respiratory tract infections in young children and is responsible for high hospitalization rates and morbidity in infants and the elderly. To understand how the emergence of RSV viral genotypes and viral-respiratory microbiome interactions contribute to infection frequency and severity, we utilized an observational cohort designed to capture the first infant RSV infection we employed multi-omics approaches to sequence 349 RSV complete genomes and matched nasopharyngeal microbiomes. We found non-G-72nt-duplicated RSV-A genotypes were more frequent in male infants, whereas G-72nt-duplicated RSV-A strains (ON1 lineage) were seen equally in both males and females. Microbiome analysis show Haemophilus was significantly more abundant in infants with RSV-A compared to infants with RSV-B infection and the microbial clustering of the abundant genera was associated with infant sex. Overall, we show sex differences in RSV genotype-nasopharyngeal microbiome, suggesting an interaction host genetics-virus-microbiome interaction.

Ferroptosis of tumour neutrophils causes immune suppression in cancer.

Ferroptosis is a non-apoptotic form of regulated cell death that is triggered by the discoordination of regulatory redox mechanisms culminating in massive peroxidation of polyunsaturated phospholipids. Ferroptosis inducers have shown considerable effectiveness in killing tumour cells in vitro, yet there has been no obvious success in experimental animal models, with the notable exception of immunodeficient mice1,2. This suggests that the effect of ferroptosis on immune cells remains poorly understood. Pathologically activated neutrophils (PMNs), termed myeloid-derived suppressor cells (PMN-MDSCs), are major negative regulators of anti-tumour immunity3-5. Here we found that PMN-MDSCs in the tumour microenvironment spontaneously die by ferroptosis. Although decreasing the presence of PMN-MDSCs, ferroptosis induces the release of oxygenated lipids and limits the activity of human and mouse T cells. In immunocompetent mice, genetic and pharmacological inhibition of ferroptosis abrogates suppressive activity of PMN-MDSCs, reduces tumour progression and synergizes with immune checkpoint blockade to suppress the tumour growth. By contrast, induction of ferroptosis in immunocompetent mice promotes tumour growth. Thus, ferroptosis is a unique and targetable immunosuppressive mechanism of PMN-MDSCs in the tumour microenvironment that can be pharmacologically modulated to limit tumour progression.

Evolutionary features of a prolific subtype of avian influenza A virus in European waterfowl.

Avian influenza A virus (AIV) is ubiquitous in waterfowl and is detected annually at high prevalence in waterfowl during the Northern Hemisphere autumn. Some AIV subtypes are globally common in waterfowl, such as H3N8, H4N6, and H6N2, and are detected in the same populations at a high frequency, annually. In order to investigate genetic features associated to the long-term maintenance of common subtypes in migratory ducks, we sequenced 248 H4 viruses isolated across 8 years (2002-9) from mallards (Anas platyrhynchos) sampled in southeast Sweden. Phylogenetic analyses showed that both H4 and N6 sequences fell into three distinct lineages, structured by year of isolation. Specifically, across the 8 years of the study, we observed lineage replacement, whereby a different HA lineage circulated in the population each year. Analysis of deduced amino acid sequences of the HA lineages illustrated key differences in regions of the globular head of hemagglutinin that overlap with established antigenic sites in homologous hemagglutinin H3, suggesting the possibility of antigenic differences among these HA lineages. Beyond HA, lineage replacement was common to all segments, such that novel genome constellations were detected across years. A dominant genome constellation would rapidly amplify in the duck population, followed by unlinking of gene segments as a result of reassortment within 2-3 weeks following introduction. These data help reveal the evolutionary dynamics exhibited by AIV on both annual and decadal scales in an important reservoir host.

Resistance to Durvalumab and Durvalumab plus Tremelimumab Is Associated with Functional STK11 Mutations in Patients with Non-Small Cell Lung Cancer and Is Reversed by STAT3 Knockdown.

Mutations in the STK11 (LKB1) gene regulate resistance to PD-1/PD-L1 blockade. This study evaluated this association in patients with nonsquamous non-small cell lung cancer (NSCLC) enrolled in three phase I/II trials. STK11 mutations were associated with resistance to the anti-PD-L1 antibody durvalumab (alone/with the anti-CTLA4 antibody tremelimumab) independently of KRAS mutational status, highlighting STK11 as a potential driver of resistance to checkpoint blockade. Retrospective assessments of tumor tissue, whole blood, and serum revealed a unique immune phenotype in patients with STK11 mutations, with increased expression of markers associated with neutrophils (i.e., CXCL2, IL6), Th17 contexture (i.e., IL17A), and immune checkpoints. Associated changes were observed in the periphery. Reduction of STAT3 in the tumor microenvironment using an antisense oligonucleotide reversed immunotherapy resistance in preclinical STK11 knockout models. These results suggest that STK11 mutations may hinder response to checkpoint blockade through mechanisms including suppressive myeloid cell biology, which could be reversed by STAT3-targeted therapy. SIGNIFICANCE: Patients with nonsquamous STK11-mutant (STK11mut) NSCLC are less likely than STK11 wild-type (STK11wt) patients to respond to anti-PD-L1 ± anti-CTLA4 immunotherapies, and their tumors show increased expression of genes and cytokines that activate STAT3 signaling. Preclinically, STAT3 modulation reverses this resistance, suggesting STAT3-targeted agents as potential combination partners for immunotherapies in STK11mut NSCLC.This article is highlighted in the In This Issue feature, p. 2659.

A novel automated SARS-CoV-2 saliva PCR test protects a global asymptomatic workforce.

Regular PCR testing of nasopharyngeal swabs from symptomatic individuals for SARS-CoV-2 virus has become the established method by which health services are managing the COVID-19 pandemic. Businesses such as AstraZeneca have also prioritised voluntary asymptomatic testing to keep workplaces safe and maintain supply of essential medicines to patients. We describe the development of an internal automated SARS-CoV-2 testing programme including the transformative introduction of saliva as an alternative sample type.

Molecular characterization of respiratory syncytial viruses circulating in a paediatric cohort in Amman, Jordan.

Respiratory syncytial viruses (RSVs) are an important cause of mortality worldwide and a major cause of respiratory tract infections in children, driving development of vaccine candidates. However, there are large gaps in our knowledge of the local evolutionary and transmission dynamics of RSVs, particularly in understudied regions such as the Middle East. To address this gap, we sequenced the complete genomes of 58 RSVA and 27 RSVB samples collected in a paediatric cohort in Amman, Jordan, between 2010 and 2013. RSVA and RSVB co-circulated during each winter epidemic of RSV in Amman, and each epidemic comprised multiple independent viral introductions of RSVA and RSVB. However, RSVA and RSVB alternated in dominance across years, potential evidence of immunological interactions. Children infected with RSVA tended to be older than RSVB-infected children [30 months versus 22.4 months, respectively (P value = 0.02)], and tended to developed bronchopneumonia less frequently than those with RSVB, although the difference was not statistically significant (P value = 0.06). Differences in spatial patterns were investigated, and RSVA lineages were often identified in multiple regions in Amman, whereas RSVB introductions did not spread beyond a single region of the city, although these findings were based on small sample sizes. Multiple RSVA genotypes were identified in Amman, including GA2 viruses as well as three viruses from the ON1 sub-genotype that emerged in 2009 and are now the dominant genotype circulating worldwide. As vaccine development advances, further sequencing of RSV is needed to understand viral ecology and transmission, particularly in under-studied locations.

Chemically Defined, High-Density Insect Cell-Based Expression System for Scalable AAV Vector Production.

The recombinant adeno-associated virus (AAV) vector is one of the most utilized viral vectors in gene therapy due to its robust, long-term in vivo transgene expression and low toxicity. One major hurdle for clinical AAV applications is large-scale manufacturing. In this regard, the baculovirus-based AAV production system is highly attractive due to its scalability and predictable biosafety. Here, we describe a simple method to improve the baculovirus-based AAV production using the ExpiSf Baculovirus Expression System with a chemically defined medium for suspension culture of high-density ExpiSf9 cells. Baculovirus-infected ExpiSf9 cells produced up to 5 × 1011 genome copies of highly purified AAV vectors per 1 mL of suspension culture, which is up to a 19-fold higher yield than the titers we obtained from the conventional Sf9 cell-based system. When mice were administered the same dose of AAV vectors, we saw comparable transduction efficiency and biodistributions between the vectors made in ExpiSf9 and Sf9 cells. Thus, the ExpiSf Baculovirus Expression System would support facile and scalable AAV manufacturing amenable for preclinical and clinical applications.

Author Correction: Microbiome disturbance and resilience dynamics of the upper respiratory tract during influenza A virus infection.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Microbiome disturbance and resilience dynamics of the upper respiratory tract during influenza A virus infection.

Infection with influenza can be aggravated by bacterial co-infections, which often results in disease exacerbation. The effects of influenza infection on the upper respiratory tract (URT) microbiome are largely unknown. Here, we report a longitudinal study to assess the temporal dynamics of the URT microbiomes of uninfected and influenza virus-infected humans and ferrets. Uninfected human patients and ferret URT microbiomes have stable healthy ecostate communities both within and between individuals. In contrast, infected patients and ferrets exhibit large changes in bacterial community composition over time and between individuals. The unhealthy ecostates of infected individuals progress towards the healthy ecostate, coinciding with viral clearance and recovery. Pseudomonadales associate statistically with the disturbed microbiomes of infected individuals. The dynamic and resilient microbiome during influenza virus infection in multiple hosts provides a compelling rationale for the maintenance of the microbiome homeostasis as a potential therapeutic target to prevent IAV associated bacterial co-infections.

Amino Acid Substitutions in Positions 385 and 393 of the Hydrophobic Region of VP4 May Be Associated with Rotavirus Attenuation and Cell Culture Adaptation.

Rotaviruses (RVs) are the leading cause of the acute viral gastroenteritis in young children and livestock animals worldwide. Although live attenuated vaccines have been applied to control RV infection for many years, the underlying mechanisms of RV attenuation following cell culture adaption are unknown. To study these mechanisms at the genomic level, we have sequenced and conducted a comparative analysis of two virulent human (Wa, G1P[8] and M, G3P[8]) and two virulent porcine (Gottfried, G4P[6] and OSU, G5P[7]) RV strains maintained in gnotobiotic piglets for 22, 11, 12 and 9 serial passages, respectively, with their attenuated counterparts serially passaged in MA-104 cell cultures for 25, 43, 54 and 43 passages, respectively. We showed that most of the mutations were clustered in the VP4 gene, with a relatively high nonsynonymous substitution rate (81.2%). Moreover, two amino acid substitutions observed in the VP4 gene were conserved between two or more strain pairs. D385N substitution was found in M, Wa and Gottfried strains, and another one, S471H/L was present in Wa and Gottfried strains. Importantly, D385 was reported previously in another study and may be involved in regulation of virus entry. Of interest, although no 385 substitution was found in OSU strains, the attenuated OSU strain contained a unique D393H substitution within the same VP4 hydrophobic domain. Collectively, our data suggest that the VP4 hydrophobic region may play an important role in RV attenuation and aa385 and aa393 may represent potential targets for RV vaccine development using reverse genetics and site-specific mutagenesis.

The tumor inflammation signature (TIS) is associated with anti-PD-1 treatment benefit in the CERTIM pan-cancer cohort.

BACKGROUND: The 18-gene tumor inflammation signature (TIS) is a clinical research assay that enriches for clinical benefit to immune checkpoint blockade. We evaluated its ability to predict clinical benefit of immunotherapy in cancer patients treated with PD-1 checkpoint inhibitors in routine clinical care. METHODS: The CERTIM cohort is a prospective cohort which includes patients receiving immune checkpoint inhibitors in Cochin University hospital. RNA extracted from 58 archival formalin fixed paraffin embedded tumor blocks (including 38 lung cancers, 5 melanomas, 10 renal carcinomas, 4 urothelial carcinomas and 1 colon carcinoma) was hybridized to a beta version of the NanoString® PanCancer IO360™ CodeSet using nCounter® technology. Gene expression signatures were correlated with tumor responses (by RECIST criteria) and overall survival. PD-L1 immunostaining on tumor cells was assessed in 37 non-small cell lung cancer (NSCLC) samples and tumor mutational burden (TMB) measured by whole exome sequencing in 19 of these. RESULTS: TIS scores were significantly associated with complete or partial response to anti-PD-1 treatment in the whole cohort (odds ratio = 2.64, 95% CI [1.4; 6.0], p = 0.008), as well as in the NSCLC population (odds ratio = 3.27, 95% CI [1.2; 11.6], p = 0.03). Patients whose tumor had a high TIS score (upper tertile) showed prolonged overall survival compared to patients whose tumor had lower TIS scores, both in the whole cohort (hazard ratio = 0.37, 95% CI [0.18, 0.76], p = 0.005) and in the NSCLC population (hazard ratio = 0.36, 95% CI [0.14, 0.90], p = 0.02). In the latter, the TIS score was independent from either PD-L1 staining on tumor cells (spearman coefficient 0.2) and TMB (spearman coefficient - 0.2). CONCLUSIONS: These results indicate that validated gene expression assay measuring the level of tumor microenvironment inflammation such as TIS, are accurate and independent predictive biomarkers and can be easily implemented in the clinical practice.

Erratum: Evolutionary dynamics and molecular epidemiology of West Nile virus in New York State: 1999-2015.

[This corrects the article DOI: 10.1093/ve/vez020.].

Sequencing and analysis of globally obtained human parainfluenza viruses 1 and 3 genomes.

Human Parainfluenza viruses (HPIV) type 1 and 3 are important causes of respiratory tract infections in young children globally. HPIV infections do not confer complete protective immunity so reinfections occur throughout life. Since no effective vaccine is available for the two virus subtypes, comprehensive understanding of HPIV-1 and HPIV-3 genetic and epidemic features is important for diagnosis, prevention, and treatment of HPIV-1 and HPIV-3 infections. Relatively few whole genome sequences are available for both HPIV-1 and HPIV-3 viruses, so our study sought to provide whole genome sequences from multiple countries to further the understanding of the global diversity of HPIV at a whole-genome level. We collected HPIV-1 and HPIV-3 samples and isolates from Argentina, Australia, France, Mexico, South Africa, Switzerland, and USA from the years 2003-2011 and sequenced the genomes of 40 HPIV-1 and 75 HPIV-3 viruses with Sanger and next-generation sequencing with the Ion Torrent, Illumina, and 454 platforms. Phylogenetic analysis showed that the HPIV-1 genome is evolving at an estimated rate of 4.97 × 10-4 mutations/site/year (95% highest posterior density 4.55 × 10-4 to 5.38 × 10-4) and the HPIV-3 genome is evolving at a similar rate (3.59 × 10-4 mutations/site/year, 95% highest posterior density 3.26 × 10-4 to 3.94 × 10-4). There were multiple genetically distinct lineages of both HPIV-1 and 3 circulating on a global scale. Further surveillance and whole-genome sequencing are greatly needed to better understand the spatial dynamics of these important respiratory viruses in humans.

Evolutionary dynamics and molecular epidemiology of West Nile virus in New York State: 1999-2015.

Following its introduction into New York State (NYS) in 1999, West Nile virus (WNV; Flavivirus, Flaviviridae) underwent a rapid expansion throughout the USA and into Canada and Latin America. WNV has been characterized as being evolutionarily stable, with weak geographic structure, a dominance of purifying selection and limited adaptive change. We analyzed all available full-genome WNV sequences, focusing on the 543 available sequences from NYS, which included 495 newly sequenced 2000-15 isolates. In addition, we analyzed deep-sequencing data from 317 of these isolates. While our data are generally in agreement with the limited pace of evolutionary change and broad geographic and temporal mixing identified in other studies, we have identified some important exceptions. Most notably, there are 14 codons which demonstrated evidence of positive selection as determined by multiple models, including some positions with evidence of selection in NYS exclusively. Coincident with increased WNV activity, genotypes possessing one or more of these mutations, designated NY01, NY07, and NY10, have increased in prevalence in recent years and displaced historic strains. In addition, we have found a geographical bias with many of these mutations, which suggests selective pressures and adaptations could be regional. Lastly, our deep-sequencing data suggest both increased overall diversity in avian tissue isolates relative to mosquito isolates and multiple non-synonymous minority variants that are both host-specific and retained over time and space. Together, these data provide novel insight into the evolutionary pressures on WNV and the need for continued genetic surveillance and characterization of emergent strains.

Evaluation of ultra-low input RNA sequencing for the study of human T cell transcriptome.

Deeper understanding of T cell biology is crucial for the development of new therapeutics. Human naïve T cells have low RNA content and their numbers can be limiting; therefore we set out to determine the parameters for robust ultra-low input RNA sequencing. We performed transcriptome profiling at different cell inputs and compared three protocols: Switching Mechanism at 5' End of RNA Template technology (SMART) with two different library preparation methods (Nextera and Clontech), and AmpliSeq technology. As the cell input decreased the number of detected coding genes decreased with SMART, while stayed constant with AmpliSeq. However, SMART enables detection of non-coding genes, which is not feasible for AmpliSeq. The detection is dependent on gene abundance, but not transcript length. The consistency between technical replicates and cell inputs was comparable across methods above 1 K but highly variable at 100 cell input. Sensitivity of detection for differentially expressed genes decreased dramatically with decreased cell inputs in all protocols, support that additional approaches, such as pathway enrichment, are important for data interpretation at ultra-low input. Finally, T cell activation signature was detected at 1 K cell input and above in all protocols, with AmpliSeq showing better detection at 100 cells.

Multiple Introductions and Antigenic Mismatch with Vaccines May Contribute to Increased Predominance of G12P[8] Rotaviruses in the United States.

Rotavirus is the leading global cause of diarrheal mortality for unvaccinated children under 5 years of age. The outer capsid of rotavirus virions consists of VP7 and VP4 proteins, which determine viral G and P types, respectively, and are primary targets of neutralizing antibodies. Successful vaccination depends upon generating broadly protective immune responses following exposure to rotaviruses presenting a limited number of G- and P-type antigens. Vaccine introduction resulted in decreased rotavirus disease burden but also coincided with the emergence of uncommon G and P genotypes, including G12. To gain insight into the recent predominance of G12P[8] rotaviruses in the United States, we evaluated 142 complete rotavirus genome sequences and metadata from 151 clinical specimens collected in Nashville, TN, from 2011 to 2013 through the New Vaccine Surveillance Network. Circulating G12P[8] strains were found to share many segments with other locally circulating strains but to have distinct constellations. Phylogenetic analyses of G12 sequences and their geographic sources provided evidence for multiple separate introductions of G12 segments into Nashville, TN. Antigenic epitopes of VP7 proteins of G12P[8] strains circulating in Nashville, TN, differ markedly from those of vaccine strains. Fully vaccinated children were found to be infected with G12P[8] strains more frequently than with other rotavirus genotypes. Multiple introductions and significant antigenic mismatch may in part explain the recent predominance of G12P[8] strains in the United States and emphasize the need for continued monitoring of rotavirus vaccine efficacy against emerging rotavirus genotypes.IMPORTANCE Rotavirus is an important cause of childhood diarrheal disease worldwide. Two immunodominant proteins of rotavirus, VP7 and VP4, determine G and P genotypes, respectively. Recently, G12P[8] rotaviruses have become increasingly predominant. By analyzing rotavirus genome sequences from stool specimens obtained in Nashville, TN, from 2011 to 2013 and globally circulating rotaviruses, we found evidence of multiple introductions of G12 genes into the area. Based on sequence polymorphisms, VP7 proteins of these viruses are predicted to present themselves to the immune system very differently than those of vaccine strains. Many of the sick children with G12P[8] rotavirus in their diarrheal stools also were fully vaccinated. Our findings emphasize the need for continued monitoring of circulating rotaviruses and the effectiveness of the vaccines against strains with emerging G and P genotypes.

Author Correction: Genomic analysis of a large set of currently-and historically-important human adenovirus pathogens.

The original publication of this article [1] was missing the below author that made contributions to the research and the published article.

Where do all the subtypes go? Temporal dynamics of H8-H12 influenza A viruses in waterfowl.

Influenza A virus (IAV) is ubiquitous in waterfowl. In the northern hemisphere IAV prevalence is highest during the autumn and coincides with a peak in viral subtype diversity. Although haemagglutinin subtypes H1-H12 are associated with waterfowl hosts, subtypes H8-H12 are detected very infrequently. To better understand the role of waterfowl in the maintenance of these rare subtypes, we sequenced H8-H12 viruses isolated from Mallards (Anas platyrhynchos) from 2002 to 2009. These rare viruses exhibited varying ecological and phylodynamic features. The Eurasian clades of H8 and H12 phylogenies were dominated by waterfowl sequences; mostly viruses sequenced in this study. H11, once believed to be a subtype that infected charadriiformes (shorebirds), exhibited patterns more typical of common virus subtypes. Finally, subtypes H9 and H10, which have maintained lineages in poultry, showed markedly different patterns: H10 was associated with all possible NA subtypes and this drove HA lineage diversity within years. Rare viruses belonging to subtypes H8-H12 were highly reassorted, indicating that these rare subtypes are part of the broader IAV pool. Our results suggest that waterfowl play a role in the maintenance of these rare subtypes, but we recommend additional sampling of non-traditional hosts to better understand the reservoirs of these rare viruses.

Differing epidemiological dynamics of Chikungunya virus in the Americas during the 2014-2015 epidemic.

Chikungunya virus (CHIKV) has been detected sporadically since the 1950s and includes three distinct co-circulating genotypes. In late 2013, the Asian genotype of CHIKV was responsible for the Caribbean outbreak (CO) that rapidly became an epidemic throughout the Americas. There is a limited understanding of the molecular evolution of CHIKV in the Americas during this epidemic. We sequenced 185 complete CHIKV genomes collected mainly from Nicaragua in Central America and Florida in the United States during the 2014-2015 Caribbean/Americas epidemic. Our comprehensive phylogenetic analyses estimated the epidemic history of the Asian genotype and the recent Caribbean outbreak (CO) clade, revealed considerable genetic diversity within the CO clade, and described different epidemiological dynamics of CHIKV in the Americas. Specifically, we identified multiple introductions in both Nicaragua and Florida, with rapid local spread of viruses in Nicaragua but limited autochthonous transmission in Florida in the US. Our phylogenetic analysis also showed phylogeographic clustering of the CO clade. In addition, we identified the significant amino acid substitutions that were observed across the entire Asian genotype during its evolution and examined amino acid changes that were specific to the CO clade. Deep sequencing analysis identified specific minor variants present in clinical specimens below-consensus levels. Finally, we investigated the association between viral phylogeny and geographic/clinical metadata in Nicaragua. To date, this study represents the largest single collection of CHIKV complete genomes during the Caribbean/Americas epidemic and significantly expands our understanding of the emergence and evolution of CHIKV CO clade in the Americas.