Genetic Structuring of One of the Main Vectors of Sylvatic Yellow Fever: Haemagogus (Conopostegus) leucocelaenus (Diptera: Culicidae).

Genetic diversity and population structuring for the species Haemogogus leucocelaenus, a sylvatic vector of yellow fever virus, were found to vary with the degree of agricultural land use and isolation of fragments of Atlantic Forest in municipalities in the state of São Paulo where specimens were collected. Genotyping of 115 mitochondrial SNPs showed that the populations with the highest indices of genetic diversity (polymorphic loci and mean pairwise differences between the sequences) are found in areas with high levels of agricultural land use (northeast of the State). Most populations exhibited statistically significant negative values for the Tajima D and Fu FS neutrality tests, suggesting recent expansion. The results show an association between genetic diversity in this species and the degree of agricultural land use in the sampled sites, as well as signs of population expansion of this species in most areas, particularly those with the highest forest edge densities. A clear association between population structuring and the distance between the sampled fragments (isolation by distance) was observed: samples from a large fragment of Atlantic Forest extending along the coast of the state of São Paulo exhibited greater similarity with each other than with populations in the northwest of the state.

Gene coexpression network during ontogeny in the yellow fever mosquito, Aedes aegypti.

BACKGROUND: The behaviors and ontogeny of Aedes aegypti are closely related to the spread of diseases caused by dengue (DENV), chikungunya (CHIKV), Zika (ZIKV), and yellow fever (YFV) viruses. During the life cycle, Ae. aegypti undergoes drastic morphological, metabolic, and functional changes triggered by gene regulation and other molecular mechanisms. Some essential regulatory factors that regulate insect ontogeny have been revealed in other species, but their roles are still poorly investigated in the mosquito. RESULTS: Our study identified 6 gene modules and their intramodular hub genes that were highly associated with the ontogeny of Ae. aegypti in the constructed network. Those modules were found to be enriched in functional roles related to cuticle development, ATP generation, digestion, immunity, pupation control, lectins, and spermatogenesis. Additionally, digestion-related pathways were activated in the larvae and adult females but suppressed in the pupae. The integrated protein‒protein network also identified cilium-related genes. In addition, we verified that the 6 intramodular hub genes encoding proteins such as EcKinase regulating larval molt were only expressed in the larval stage. Quantitative RT‒PCR of the intramodular hub genes gave similar results as the RNA-Seq expression profile, and most hub genes were ontogeny-specifically expressed. CONCLUSIONS: The constructed gene coexpression network provides a useful resource for network-based data mining to identify candidate genes for functional studies. Ultimately, these findings will be key in identifying potential molecular targets for disease control.

Mosquito gene targeted RNAi studies for vector control.

Vector-borne diseases are serious public health concern. Mosquito is one of the major vectors responsible for the transmission of a number of diseases like malaria, Zika, chikungunya, dengue, West Nile fever, Japanese encephalitis, St. Louis encephalitis, and yellow fever. Various strategies have been used for mosquito control, but the breeding potential of mosquitoes is such tremendous that most of the strategies failed to control the mosquito population. In 2020, outbreaks of dengue, yellow fever, and Japanese encephalitis have occurred worldwide. Continuous insecticide use resulted in strong resistance and disturbed the ecosystem. RNA interference is one of the strategies opted for mosquito control. There are a number of mosquito genes whose inhibition affected mosquito survival and reproduction. Such kind of genes could be used as bioinsecticides for vector control without disturbing the natural ecosystem. Several studies have targeted mosquito genes at different developmental stages by the RNAi mechanism and result in vector control. In the present review, we included RNAi studies conducted for vector control by targeting mosquito genes at different developmental stages using different delivery methods. The review could help the researcher to find out novel genes of mosquitoes for vector control.

CRISPR-Cas9 Methods and Key Considerations in the Production of Aedes aegypti Mutant Strains.

Since the characterization of the CRISPR-Cas9 system in prokaryotes, it has become the prime choice in gene editing because of its exceptional flexibility, ease of use, high efficiency, and superior specificity. As a result, CRISPR-Cas9-mediated gene-editing technologies have enabled researchers not only to engineer transgenic animal strains with site-directed insertions more efficiently but also to generate desired mutants for previously intractable species. One such species is the invasive yellow fever mosquito, Aedes aegypti, which is notorious for its ability to transmit many blood-borne human pathogens. Methods for developing new transgenic strains of the yellow fever mosquito may aid in the effort to control its populations and provide significant benefits for the public. Here, we provide an overview of injection and noninjection methods for generating transgenic mosquitoes and also highlight important experimental design features.

RNA interference is essential to modulating the pathogenesis of mosquito-borne viruses in the yellow fever mosquito Aedes aegypti.

While it has long been known that the transmission of mosquito-borne viruses depends on the establishment of persistent and nonlethal infections in the invertebrate host, specific roles for the insects' antiviral immune pathways in modulating the pathogenesis of viral infections is the subject of speculation and debate. Here, we show that a loss-of-function mutation in the Aedes aegypti Dicer-2 (Dcr-2) gene renders the insect acutely susceptible to a disease phenotype upon infection with pathogens in multiple virus families associated with important human diseases. Additional interrogation of the disease phenotype demonstrated that the virus-induced pathology is controlled through a canonical RNA interference (RNAi) pathway, which functions as a resistance mechanism. These results suggest comparatively modest contributions of proposed tolerance mechanisms to the fitness of A. aegypti infected with these pathogens. Similarly, the production of virus-derived piwi-interacting RNAs (vpiRNAs) was not sufficient to prevent the pathology associated with viral infections in Dcr-2 null mutants, also suggesting a less critical, or potentially secondary, role for vpiRNAs in antiviral immunity. These findings have important implications for understanding the ecological and evolutionary interactions occurring between A. aegypti and the pathogens they transmit to human and animal hosts.

MEK/ERK activation plays a decisive role in Zika virus morphogenesis and release.

Brazil has experienced an increase in outbreaks caused by flaviviruses. The high incidence of dengue fever, the morbidity of Zika in children, and the high mortality of yellow fever have affected millions in recent years. Deciphering host-virus interactions is important for treating viral infections, and the mitogen-activated protein kinases (MAPK) are an interesting target because of their role in flavivirus replication. In particular, mitogen-activated protein kinase kinase (MEK), which targets extracellular-signal-regulated kinase (ERK), is necessary for dengue and yellow fever infections. In this study, we evaluated the role of the MEK/ERK pathway and the effect of the MEK inhibitor trametinib on the Asian ZIKV strain PE243 and the prototype African ZIKV strain MR766, addressing genome replication, morphogenesis, and viral release. ZIKV infection stimulated ERK phosphorylation in Vero cells at 12 and 18 hours postinfection (hpi). Trametinib showed sustained antiviral activity, inhibiting both ZIKV strains for at least four days, and electron microscopy showed probable inhibition of ZIKV morphogenesis. ZIKV PE243 can complete one cycle in Vero cells in 14 hours; genome replication was detected around 8 hpi, intracellular viral particles at 12 hpi, and extracellular progeny at 14 hpi. Treatments at 6-hour intervals showed that trametinib inhibited late stages of viral replication, and the titration of intra- or extracellular virions showed that the treatment especially affected viral morphogenesis and release. Thus, ZIKV stimulated ERK phosphorylation during viral morphogenesis and release, which correlated with trametinib inhibiting both the signaling pathway and viral replication.

CRISPR-Cas9 Genome Editing Uncovers the Mode of Action of Methoprene in the Yellow Fever Mosquito, Aedes aegypti.

Methoprene, a juvenile hormone (JH) analog, is widely used for insect control, but its mode of action is not known. To study methoprene action in the yellow fever mosquito, Aedes aegypti, the E93 (ecdysone-induced transcription factor) was knocked out using the CRISPR-Cas9 system. The E93 mutant pupae retained larval tissues similar to methoprene-treated insects. These insects completed pupal ecdysis and died as pupa. In addition, the expression of transcription factors, broad complex and Krüppel homolog 1 (Kr-h1), increased and that of programmed cell death (PCD) and autophagy genes decreased in E93 mutants. These data suggest that methoprene functions through JH receptor, methoprene-tolerant, and induces the expression of Kr-h1, which suppresses the expression of E93, resulting in a block in PCD and autophagy of larval tissues. Failure in the elimination of larval tissues and the formation of adult structures results in their death. These results answered long-standing questions on the mode of action of methoprene.

A chromosome-level assembly of the widely used Rockefeller strain of Aedes aegypti, the yellow fever mosquito.

Aedes aegypti is the vector of important human diseases, and genomic resources are crucial in facilitating the study of A. aegypti and its ecosystem interactions. Several laboratory-acclimated strains of this mosquito have been established, but the most used strain in toxicology studies is "Rockefeller," which was originally collected and established in Cuba 130 years ago. A full-length genome assembly of another reference strain, "Liverpool," was published in 2018 and is the reference genome for the species (AaegL5). However, genetic studies with the Rockefeller strain are complicated by the availability of only the Liverpool strain as the reference genome. Differences between Liverpool and Rockefeller have been known for decades, particularly in the expression of genes relevant to mosquito behavior and vector control (e.g. olfactory). These differences indicate that AaegL5 is likely not fully representative of the Rockefeller genome, presenting potential impediments to research. Here, we present a chromosomal-level assembly and annotation of the Rockefeller genome and a comparative characterization vs the Liverpool genome. Our results set the stage for a pan-genomic approach to understanding evolution and diversity within this important disease vector.

Real-Time Genomic Surveillance during the 2021 Re-Emergence of the Yellow Fever Virus in Rio Grande do Sul State, Brazil.

The 2021 re-emergence of yellow fever in non-human primates in the state of Rio Grande do Sul (RS), southernmost Brazil, resulted in the death of many howler monkeys (genus Alouatta) and led the state to declare a Public Health Emergency of State Importance, despite no human cases reported. In this study, near-complete genomes of yellow fever virus (YFV) recovered from the outbreak were sequenced and examined aiming at a better understanding of the phylogenetic relationships and the spatio-temporal dynamics of the virus distribution. Our results suggest that the most likely sequence of events involved the reintroduction of YFV from the state of São Paulo to RS through the states of Paraná and Santa Catarina, by the end of 2020. These findings reinforce the role of genomic surveillance in determining the pathways of distribution of the virus and in providing references for the implementation of preventive measures for populations in high risk areas.

Hormonal regulation of microRNA expression dynamics in the gut of the yellow fever mosquito Aedes aegypti.

The yellow fever mosquito Aedes aegypti is an obligatory blood feeder and a major arboviral disease vector, evoking severe public health concerns worldwide. In adult female mosquitoes, the gut is critical for blood digestion and pathogen entry. We aimed for a systematic exploration of microRNA expression dynamics in the gut during the gonadotrophic cycle. Small RNA libraries were constructed from female mosquito gut tissues at five time points. Unsupervised hierarchical clustering revealed three expression clusters (early, mid and late) peaking at sequential time points - 24, 48 and 72 h posteclosion. Differentially expressed miRNAs were identified at 24 h post-blood meal (PBM). Depletions of Methoprene-tolerant [Met; the juvenile hormone (JH) receptor] and Ecdysone receptor [EcR; the receptor to 20-hydroxyecdysone (20E)] were performed using dsRNA to these genes to investigate impacts on microRNA expressions. Our results suggest that Met-mediated signalling downregulates miRNA expression from the early cluster and upregulates that from the late cluster. EcR signalling either up- or downregulated miRNA levels at 24 h PBM, indicating a differential effect of this receptor in miRNA gene expression. Furthermore, miR-281, which is the most abundant miRNA in the gut tissue, is induced and repressed by Met- and EcR-mediated signalling, respectively. Systematic depletion using synthetic antagomir and phenotype examinations indicate that miR-281 is obligatory for the normal progression of blood digestion, ovarian development and reproduction. Collectively, this study unveils expression dynamics of microRNAs in the female gut tissue during the gonadotrophic cycle and demonstrates that they are affected by JH and 20E signalling.

Morphology and taxonomic status of Aedes aegypti populations across Senegal.

Aedes aegypti is the primary vector of dengue, Zika, yellow fever and chikungunya viruses to humans. In Africa, two subspecies, Ae. aegypti aegypti (Aaa) and Ae. aegypti formosus (Aaf) have been described. Until very recently, it was considered that the two forms were sympatric in East Africa and that only Aaf was present in Central and West Africa. However, recent data suggests that Aaa was also common in Senegal without any clear evidence of genetic differences with Aaf. This study was carried out in different Ae. aegypti populations from Senegal to better clarify their taxonomic status. The larvae, pupae and eggs were collected between July and September 2018 and reared individually to adult stage. For each population, F1 progeny from eggs laid by a single female F0 were reared as sibling samples. The number of pale scales on the first abdominal tergite (T1) and the basal part of the second tergite (T2) were counted. Individuals with no pale scale on T1 were classified as Aaf while those with at least one pale scale on this tergite were classified as Aaa. The morphological variations within families of Aaf were studied across 4 generations. In total, 2400 individuals constituting 240 families were identified, of which 42.5% were heterogeneous (families with both forms). Multivariate statistical analysis of variance including T1 and T2 data together showed that populations were significantly different from each other. Statistical analysis of T1 alone showed a similarity between populations from the southeast while variations were observed within northwest population. The analysis of family composition across generations showed the presence of Aaa and Aaf forms in each generation. The classification of Ae. aegypti into two subspecies is invalid in Senegal. Populations exhibit morphological polymorphism at the intra-family level that could have biological and epidemiological impacts.

Comparing the selective landscape of TLR7 and TLR8 across primates reveals unique sites under positive selection in Alouatta.

Among primates, susceptibility to yellow fever (YFV), a single-stranded (ss) RNA virus, ranges from complete resistance to high susceptibility. Howler monkeys (genus Alouatta) are the most susceptible to YFV. In order to identify Alouatta-specific genetic factors that may be responsible for their susceptibility, we collected skin samples from howler monkey museum specimens of the species A. caraya and A. guariba clamitans. We compared the rate of nonsynonymous to synonymous (dN/dS) changes of Toll-like receptor (TLR) 7 and TLR8, the two genes responsible for detecting all ssRNA viruses, across the Primate order. Overall, we found that the TLR7 gene is under stronger purifying selection in howler monkeys compared to other New World and Old World primates, but TLR8 is under the same selective pressure. When we evaluated dN/dS at each codon, we found six codons under positive selection in Alouatta TLR8 and two codons under positive selection in TLR7. The changes in TLR7 are unique to A. guariba clamitans and are found in functionally important regions likely to affect detection of ssRNA viruses by TLR7/TLR8, as well as downstream signaling. These amino acid differences in A. guariba clamitans may play a role in YFV susceptibility. These results have implications for identifying genetic factors affecting YFV susceptibility in primates.

Positively selected variants in functionally important regions of TLR7 in Alouatta guariba clamitans with yellow fever virus exposure in Northern Argentina.

BACKGROUND: In 2007-2009, a major yellow fever virus (YFV) outbreak in Northern Argentina decimated the local howler monkey (Alouatta) population. AIMS: To evaluate whether the surviving howler monkeys possess advantageous genetic variants inherited from monkeys alive prior to the YFV outbreak, we explored the relationship between Toll-like receptor (TLR) 7 and TLR8 gene variation and YFV susceptibility. METHODS: We used samples from Alouatta individuals in Misiones, Argentina alive before the YFV outbreak, individuals that died during the outbreak, and individuals that survived the outbreak and are alive today. We measured genetic divergence between Alouatta YFV exposure groups and evaluated Alouatta-specific substitutions for functional consequences. RESULTS: We did not find different allele frequencies in the post-YFV exposure Alouatta group compared to the pre-exposure group. We identified three nonsynonymous variants in TLR7 in Alouatta guariba clamitans. Two of these substitutions are under positive selection in functionally important regions of the gene. DISCUSSION AND CONCLUSIONS: Our results did not indicate that surviving howler monkey spossess advantageous genetic variants at greater frequency than those alive before the YFV outbreak. However, the positively selected unique coding differences in A. guariba clamitans are in the region important in pathogen detection which may affect YFV resistance. Morework is necessary to fully explore this hypothesis.

CREB-binding protein regulates metamorphosis and compound eye development in the yellow fever mosquito, Aedes aegypti.

Juvenile hormones (JH) and ecdysone coordinately regulate metamorphosis in Aedes aegypti. We studied the function of an epigenetic regulator and multifunctional transactivator, CREB binding protein (CBP) in A. aegypti. RNAi-mediated knockdown of CBP in Ae. aegypti larvae resulted in suppression of JH primary response gene, Krüppel-homolog 1 (Kr-h1), and induction of primary ecdysone response gene, E93, resulting in multiple effects including early metamorphosis, larval-pupal intermediate formation, mortality and inhibition of compound eye development. RNA sequencing identified hundreds of genes, including JH and ecdysone response genes regulated by CBP. In the presence of JH, CBP upregulates Kr-h1 by acetylating core histones at the Kr-h1 promoter and facilitating the recruitment of JH receptor and other proteins. CBP suppresses metamorphosis regulators, EcR-A, USP-A, BR-C, and E93 through the upregulation of Kr-h1 and E75A. CBP regulates the expression of core eye specification genes including those involved in TGF-ß and EGFR signaling. These studies demonstrate that CBP is an essential player in JH and 20E action and regulates metamorphosis and compound eye development in Ae. aegypti.

Sestrins induce natural killer function in senescent-like CD8+ T cells.

Aging is associated with remodeling of the immune system to enable the maintenance of life-long immunity. In the CD8+ T cell compartment, aging results in the expansion of highly differentiated cells that exhibit characteristics of cellular senescence. Here we found that CD27-CD28-CD8+ T cells lost the signaling activity of the T cell antigen receptor (TCR) and expressed a protein complex containing the agonistic natural killer (NK) receptor NKG2D and the NK adaptor molecule DAP12, which promoted cytotoxicity against cells that expressed NKG2D ligands. Immunoprecipitation and imaging cytometry indicated that the NKG2D-DAP12 complex was associated with sestrin 2. The genetic inhibition of sestrin 2 resulted in decreased expression of NKG2D and DAP12 and restored TCR signaling in senescent-like CD27-CD28-CD8+ T cells. Therefore, during aging, sestrins induce the reprogramming of non-proliferative senescent-like CD27-CD28-CD8+ T cells to acquire a broad-spectrum, innate-like killing activity.

Real-time Metagenomic Analysis of Undiagnosed Fever Cases Unveils a Yellow Fever Outbreak in Edo State, Nigeria.

Fifty patients with unexplained fever and poor outcomes presented at Irrua Specialist Teaching Hospital (ISTH) in Edo State, Nigeria, an area endemic for Lassa fever, between September 2018 - January 2019. After ruling out Lassa fever, plasma samples from these epidemiologically-linked cases were sent to the African Centre of Excellence for Genomics of Infectious Diseases (ACEGID), Redeemer's University, Ede, Osun State, Nigeria, where we carried out metagenomic sequencing which implicated yellow fever virus (YFV) as the etiology of this outbreak. Twenty-nine of the 50 samples were confirmed positive for YFV by reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), 14 of which resulted in genome assembly. Maximum likelihood phylogenetic analysis revealed that these YFV sequences formed a tightly clustered clade more closely related to sequences from Senegal than sequences from earlier Nigerian isolates, suggesting that the YFV clade responsible for this outbreak in Edo State does not descend directly from the Nigerian YFV outbreaks of the last century, but instead reflects a broader diversity and dynamics of YFV in West Africa. Here we demonstrate the power of metagenomic sequencing for identifying ongoing outbreaks and their etiologies and informing real-time public health responses, resulting in accurate and prompt disease management and control.

[Development of inactivated cultural yellow fever vaccine].

INTRODUCTION: The only currently available live vaccine against yellow fever (YF) based on chicken embryos infected with an attenuated 17D strain of the YF virus is one of the most effective vaccine preparations. However, the live vaccine is associated with "viscerotropic syndrome" (approximately 0.4 cases per 100 000 vaccinated). Therefore, the development and introduction of highly purified inactivated vaccine against YF is intended to ensure the maximum safety of vaccination against one of the most common human viral diseases.Goals and objectives. Development and evaluation of immunogenicity of the cultural inactivated vaccine against YF at the laboratory model level. MATERIAL AND METHODS: Adaptation of 17D strain of YF virus to Vero cell culture, cultivation, removal of cellular DNA, inactivation with ß-propiolactone, concentration, chromatographic purification, determination of protein and antigen of YF virus, assessment of immunogenicity in mice in parallel with commercial live vaccine. RESULTS AND DISCUSSION: Immunogenicity: the determination of specific antibodies of class G (IgG) and virus neutralizing antibodies in the sera of immunized mice showed high level of antibodies exceeding that of immunized with commercial live vaccine. The optimal dose of antigen in the vaccine (total protein) was 50 µg/ml (5 µg/0.1 ml -dose and volume per 1 vaccination of mice). Thus, the laboratory version of cultural inactivated vaccine against YF is as effective (and even superior) as the commercial live vaccine. CONCLUSION: Laboratory version of cultural inactivated vaccine against YF, which is not inferior in immunogenicity (in animal model) to commercial live vaccine, has been developed.

Genetic evidence for the origin of Aedes aegypti, the yellow fever mosquito, in the southwestern Indian Ocean.

Aedes aegypti is among the best-studied mosquitoes due to its critical role as a vector of human pathogens and ease of laboratory rearing. Until now, this species was thought to have originated in continental Africa, and subsequently colonized much of the world following the establishment of global trade routes. However, populations of this mosquito on the islands in the southwestern Indian Ocean (SWIO), where the species occurs with its nearest relatives referred to as the Aegypti Group, have received little study. We re-evaluated the evolutionary history of Ae. aegypti and these relatives, using three data sets: nucleotide sequence data, 18,489 SNPs and 12 microsatellites. We found that: (a) the Aegypti Group diverged 16 MYA (95% HPD: 7-28 MYA) from its nearest African/Asian ancestor; (b) SWIO populations of Ae. aegypti are basal to continental African populations; (c) after diverging 7 MYA (95% HPD: 4-15 MYA) from its nearest formally described relative (Ae. mascarensis), Ae. aegypti moved to continental Africa less than 85,000 years ago, where it recently (<1,000 years ago) split into two recognized subspecies Ae. aegypti formosus and a human commensal, Ae. aegypti aegypti; (d) the Madagascar samples form a clade more distant from all other Ae. aegypti than the named species Ae. mascarensis, implying that Madagascar may harbour a new cryptic species; and (e) there is evidence of introgression between Ae. mascarensis and Ae. aegypti on Réunion, and between the two subspecies elsewhere in the SWIO, a likely consequence of recent introductions of domestic Ae. aegypti aegypti from Asia.

Zika and Flavivirus Shell Disorder: Virulence and Fetal Morbidity.

Zika virus (ZIKV) was first discovered in 1947 in Africa. Since then, sporadic ZIKV infections of humans have been reported in Africa and Asia. For a long time, this virus was mostly unnoticed due to its mild symptoms and low fatality rates. However, during the 2015-2016 epidemic in Central and South America, when millions of people were infected, it was discovered that ZIKV causes microcephaly in the babies of mothers infected during pregnancy. An examination of the M and C proteins of the ZIKV shell using the disorder predictor PONDR VLXT revealed that the M protein contains relatively high disorder levels comparable only to those of the yellow fever virus (YFV). On the other hand, the disorder levels in the C protein are relatively low, which can account for the low case fatality rate (CFR) of this virus in contrast to the more virulent YFV, which is characterized by high disorder in its C protein. A larger variation was found in the percentage of intrinsic disorder (PID) in the C protein of various ZIKV strains. Strains of African lineage are characterized by higher PIDs. Using both in vivo and in vitro experiments, laboratories have also previously shown that strains of African origin have a greater potential to inflict higher fetal morbidity than do strains of Asian lineage, with dengue-2 virus (DENV-2) having the least potential. Strong correlations were found between the potential to inflict fetal morbidity and shell disorder in ZIKV (r2 = 0.9) and DENV-2 (DENV-2 + ZIKV, r2 = 0.8). A strong correlation between CFR and PID was also observed when ZIKV was included in an analysis of sets of shell proteins from a variety of flaviviruses (r2 = 0.8). These observations have potential implications for antiviral vaccine development and for the design of cancer therapeutics in terms of developing therapeutic viruses that penetrate hard-to-reach organs.

Analysis By Deep Sequencing of Discontinued Neurotropic Yellow Fever Vaccine Strains.

Deep sequencing of live-attenuated viral vaccines has focused on vaccines in current use. Here we report characterization of a discontinued live yellow fever (YF) vaccine associated with severe adverse events. The French neurotropic vaccine (FNV) strain of YF virus was derived empirically in 1930 by 260 passages of wild-type French viscerotropic virus (FVV) in mouse brain. The vaccine was administered extensively in French-speaking Africa until discontinuation in 1982, due to high rates of post-vaccination encephalitis in children. Using rare archive strains of FNV, viral RNAs were sequenced and analyzed by massively parallel, in silico methods. Diversity and specific population structures were compared in reference to the wild-type parental strain FVV, and between the vaccine strains themselves. Lower abundance of polymorphism content was observed for FNV strains relative to FVV. Although the vaccines were of lower diversity than FVV, heterogeneity between the vaccines was observed. Reversion to wild-type identity was variably observed in the FNV strains. Specific population structures were recovered from vaccines with neurotropic properties; loss of neurotropism in mice was associated with abundance of wild-type RNA populations. The analysis provides novel sequence evidence that FNV is genetically unstable, and that adaptation of FNV contributed to the neurotropic adverse phenotype.