Search

Background: Malaria transmission in Tanzania has declined significantly over the last 2 decades due to scaled-up control interventions. However, recent confirmation of artemisinin partial resistance (ART-R) in Kagera region in northwest Tanzania threatens the ongoing efforts to eliminate malaria in the country. This study was conducted according to the World Health Organization recommendation to generate evidence of malaria burden in areas with confirmed ART-R as the first step before developing a response strategy to the resistance., Methods: We assessed the local burden of malaria in Kagera region by geospatial analysis, using data collected retrospectively from health facilities and community surveys from 2015 to 2023 to identify malaria hot spots., Results: From 2017 to 2023, a total of 8 124 363 suspected malaria cases were reported by health facilities, and 2 983 717 (36.7% [95% range across wards, 22.7%-50.7%]) tested positive by rapid diagnostic tests. Test positivity rates were similar among patients aged <5 years (33.1% [95% range, 19.7%-46.5%]) and those aged ≥5 years (33.7% [21.0%-46.5%]). The malaria prevalence was 10.0% (95% range across wards, 5.1%-14.9% [n = 84 999 of 853 761]) in pregnant women and 26.1% (11.7%-40.6% [n = 3409 of 13 065]) in schoolchildren. Despite high temporal variations, we identified hot spots and cold spots, including persistently high burden in 69 of 192 wards (35.9%)., Conclusions: The malaria burden in Kagera exhibited high temporal and spatial heterogeneity, with schoolchildren showing the highest prevalence. This demographic pattern underlines the need for targeted interventions and provides evidence for developing an ART-R response for the region., Competing Interests: Potential conflicts of interest. All authors: No reported conflicts., (© The Author(s) 2024. Published by Oxford University Press on behalf of Infectious Diseases Society of America.)

Dengue is endemic in Vietnam with circulation of all four serotypes (DENV1-4) all year-round. It is hard to estimate the disease's true serotype-specific transmission patterns from cases due to its high asymptomatic rate, low reporting rate and complex immunity and transmission dynamics. Seroprevalence studies have been used to great effect for understanding patterns of dengue transmission. We tested 991 population serum samples (ages 1-30 years, collected 2013 to 2017), 531 from Ho Chi Minh City and 460 from Khanh Hoa in Vietnam, using a flavivirus protein microarray assay. By applying our previously developed inference framework to the antibody profiles from this assay, we can (1) determine proportions of a population that have not been infected or infected, once, or more than once, and (2) infer the infecting serotype in those infected once. With these data, we then use mathematical models to estimate the force of infection (FOI) for all four DENV serotypes in HCMC and KH over 35 years up to 2017. Models with time-varying or serotype-specific DENV FOI assumptions fit the data better than constant FOI. Annual dengue FOI ranged from 0.005 (95%CI: 0.003-0.008) to 0.201 (95%CI: 0.174-0.228). FOI varied across serotypes, higher for DENV1 (95%CI: 0.033-0.048) and DENV2 (95%CI: 0.018-0.039) than DENV3 (95%CI: 0.007-0.010) and DENV4 (95%CI: 0.010-0.016). The use of the PMA on serial age-stratified cross-sectional samples increases the amount of information on transmission and population immunity, and should be considered for future dengue serological surveys, particularly to understand population immunity given vaccines with differential efficacy against serotypes, however, there remains limits to what can be inferred even using this assay., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Phuong et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Influenza virus seasonality, synchronicity, and vaccine supply differ substantially between temperate and tropical settings, and optimal vaccination strategy may differ on this basis. Many national vaccine recommendations focus on high-risk groups, elderly populations, and healthcare workers despite previous analyses demonstrating broad benefits to vaccinating younger high-contact age groups. In this study, we parameterized an age-structured nonseasonal asynchronous epidemiological model of influenza virus transmission for a tropical low-income setting. We evaluated timing and age allocation of vaccines across vaccine supplies ranging from 10 to 90% using decade-based age groups. Year-round vaccination was beneficial when compared with more concentrated annual vaccine distribution. When targeting a single age group for vaccine prioritization, maximum vaccine allocation to the 10-19 high-contact age group minimized annual influenza mortality for all but one vaccine supply. When evaluating across all possible age allocations, optimal strategies always allocated a plurality of vaccines to school-age children (10-19). The converse, however, was not true as not all strategies allocating a plurality to children aged 10-19 minimized mortality. Allocating a high proportion of vaccine supply to the 10-19 age group is necessary but not sufficient to minimize annual mortality as distribution of remaining vaccine doses to other age groups also needs to be optimized. Strategies focusing on indirect benefits (vaccinating children) showed higher variance in mortality outcomes than strategies focusing on direct benefits (vaccinating the elderly). However, the indirect benefit approaches showed a lower mean mortality and a lower minimum mortality than vaccination focused on the elderly., (© The Author(s) 2024. Published by Oxford University Press on behalf of National Academy of Sciences.)

Human influenza virus evolves to escape neutralization by polyclonal antibodies. However, we have a limited understanding of how the antigenic effects of viral mutations vary across the human population and how this heterogeneity affects virus evolution. Here, we use deep mutational scanning to map how mutations to the hemagglutinin (HA) proteins of two H3N2 strains, A/Hong Kong/45/2019 and A/Perth/16/2009, affect neutralization by serum from individuals of a variety of ages. The effects of HA mutations on serum neutralization differ across age groups in ways that can be partially rationalized in terms of exposure histories. Mutations that were fixed in influenza variants after 2020 cause greater escape from sera from younger individuals compared with adults. Overall, these results demonstrate that influenza faces distinct antigenic selection regimes from different age groups and suggest approaches to understand how this heterogeneous selection shapes viral evolution., Competing Interests: Declaration of interests J.D.B. is on the scientific advisory boards of Apriori Bio, Invivyd, Aerium Therapeutics, and the Vaccine Company. J.D.B., A.N.L., and F.C.W. receive royalty payments as inventors on Fred Hutch licensed patents related to viral deep mutational scanning. H.Y.C. reports consulting with Ellume, Pfizer, and the Bill and Melinda Gates Foundation. She has served on advisory boards for Vir, Merck, and Abbvie. She has conducted CME teaching with Medscape, Vindico, and Clinical Care Options. She has received research funding from Gates Ventures and support and reagents from Ellume and Cepheid outside of the submitted work. J.A.E. receives institutional funding from AstraZeneca, Merck, GlaxoSmithKline, and Pfizer. She is a consultant for Abbvie, AstraZeneca, GlaxoSmithKline, Meissa Vaccines, Moderna, Pfizer, and SanofiPasteur. A.L.G. reports contract testing from Abbott, Cepheid, Novavax, Pfizer, Janssen, and Hologic and research support from Gilead outside of the submitted work., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Historically Plasmodium falciparum has followed a pattern of drug resistance first appearing in low-transmission settings before spreading to high-transmission settings. Several features of low-transmission regions are hypothesized as explanations: higher chance of symptoms and treatment seeking, better treatment access, less within-host competition among clones and lower rates of recombination. Here, we test whether importation of drug-resistant parasites is more likely to lead to successful emergence and establishment in low-transmission or high-transmission periods of the same epidemiological setting, using a spatial, individual-based stochastic model of malaria and drug-resistance evolution calibrated for Burkina Faso. Upon controlling for the timing of importation of drug-resistant genotypes and examination of key model variables, we found that drug-resistant genotypes imported during the low-transmission season were (i) more susceptible to stochastic extinction due to the action of genetic drift, and (ii) more likely to lead to establishment of drug resistance when parasites are able to survive early stochastic loss due to drift. This implies that rare importation events are more likely to lead to establishment if they occur during a high-transmission season, but that constant importation (e.g. neighbouring countries with high levels of resistance) may produce a greater risk during low-transmission periods.

Delaying and slowing antimalarial drug resistance evolution is a priority for malaria-endemic countries. Until novel therapies become available, the mainstay of antimalarial treatment will continue to be artemisinin-based combination therapy (ACT). Deployment of different ACTs can be optimized to minimize evolutionary pressure for drug resistance by deploying them as a set of co-equal multiple first-line therapies (MFT) rather than rotating therapies in and out of use. Here, we consider one potential detriment of MFT policies, namely, that the simultaneous deployment of multiple ACTs could drive the evolution of different resistance alleles concurrently and that these resistance alleles could then be brought together by recombination into double-resistant or triple-resistant parasites. Using an individual-based model, we compare MFT and cycling policies in malaria transmission settings ranging from 0.1% to 50% prevalence. We define a total risk measure for multi-drug resistance (MDR) by summing the area under the genotype-frequency curves (AUC) of double- and triple-resistant genotypes. When prevalence ≥ 1%, total MDR risk ranges from statistically similar to 80% lower under MFT policies than under cycling policies, irrespective of whether resistance is imported or emerges de novo. At 0.1% prevalence, there is little statistical difference in MDR risk between MFT and cycling., (© 2024. The Author(s).)

Artemisinin combination therapies (ACTs) are highly effective at treating uncomplicated Plasmodium falciparum malaria, but the emergence of the new pfkelch13 R561H mutation in Rwanda, associated with delayed parasite clearance, suggests that interventions are needed to slow its spread. Using a Rwanda-specific spatial calibration of an individual-based malaria model, we evaluate 26 strategies aimed at minimizing treatment failures and delaying the spread of R561H after 3, 5 and 10 years. Lengthening ACT courses and deploying multiple first-line therapies (MFTs) reduced treatment failures after 5 years when compared to the current approach of a 3-d course of artemether-lumefantrine. The best among these options (an MFT policy) resulted in median treatment failure counts that were 49% lower and a median R561H allele frequency that was 0.15 lower than under baseline. New approaches to resistance management, such as triple ACTs or sequential courses of two different ACTs, were projected to have a larger impact than longer ACT courses or MFT; these were associated with median treatment failure counts in 5 years that were 81-92% lower than the current approach. A policy response to currently circulating artemisinin-resistant genotypes in Africa is urgently needed to prevent a population-wide rise in treatment failures., (© 2023. The Author(s).)

Introduction: It is well known that influenza and other respiratory viruses are wintertime-seasonal in temperate regions. However, respiratory disease seasonality in the tropics is less well understood. In this study, we aimed to characterise the seasonality of influenza-like illness (ILI) and influenza virus in Ho Chi Minh City, Vietnam., Methods: We monitored the daily number of ILI patients in 89 outpatient clinics from January 2010 to December 2019. We collected nasal swabs and tested for influenza from a subset of clinics from May 2012 to December 2019. We used spectral analysis to describe the periodic signals in the system. We evaluated the contribution of these periodic signals to predicting ILI and influenza patterns through lognormal and gamma hurdle models., Results: During 10 years of community surveillance, 66 799 ILI reports were collected covering 2.9 million patient visits; 2604 nasal swabs were collected, 559 of which were PCR-positive for influenza virus. Both annual and nonannual cycles were detected in the ILI time series, with the annual cycle showing 8.9% lower ILI activity (95% CI 8.8% to 9.0%) from February 24 to May 15. Nonannual cycles had substantial explanatory power for ILI trends (ΔAIC=183) compared with all annual covariates (ΔAIC=263) in lognormal regression. Near-annual signals were observed for PCR-confirmed influenza but were not consistent over time or across influenza (sub)types. The explanatory power of climate factors for ILI and influenza virus trends was weak., Conclusion: Our study reveals a unique pattern of respiratory disease dynamics in a tropical setting influenced by both annual and nonannual drivers, with influenza dynamics showing near-annual periodicities. Timing of vaccination campaigns and hospital capacity planning may require a complex forecasting approach., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2023. Re-use permitted under CC BY. Published by BMJ.)

Background: As we continue the fourth year of the COVID-19 epidemic, SARS-CoV-2 infections still cause high morbidity and mortality in the United States. During 2020-2022, COVID-19 was one of the leading causes of death in the United States and by far the leading cause among infectious diseases. Vaccination uptake remains low despite this being an effective burden reducing intervention. The development of COVID-19 therapeutics provides hope for mitigating severe clinical outcomes. This modeling study examines combined strategies of vaccination and treatment to reduce the burden of COVID-19 epidemics over the next decade., Methods: We use a validated mathematical model to evaluate the reduction of incident cases, hospitalized cases, and deaths in the United States through 2033 under various levels of vaccination and treatment coverage. We assume that future seasonal transmission patterns for COVID-19 will be similar to those of influenza virus and account for the waning of infection-induced immunity and vaccine-induced immunity in a future with stable COVID-19 dynamics. Due to uncertainty in the duration of immunity following vaccination or infection, we consider three exponentially distributed waning rates, with means of 365 days (1 year), 548 days (1.5 years), and 730 days (2 years). We also consider treatment failure, including rebound frequency, as a possible treatment outcome., Results: As expected, universal vaccination is projected to eliminate transmission and mortality. Under current treatment coverage (13.7%) and vaccination coverage (49%), averages of 81,000-164,600 annual reported deaths, depending on duration of immunity, are expected by the end of this decade. Annual mortality in the United States can be reduced below 50,000 per year with 52-80% annual vaccination coverage and below 10,000 annual deaths with 59-83% annual vaccination coverage, depending on duration of immunity. Universal treatment reduces hospitalizations by 88.6% and deaths by 93.1% under current vaccination coverage. A reduction in vaccination coverage requires a comparatively larger increase in treatment coverage in order for hospitalization and mortality levels to remain unchanged., Conclusions: Adopting universal vaccination and universal treatment goals in the United States will likely lead to a COVID-19 mortality burden below 50,000 deaths per year, a burden comparable to that of influenza virus., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Swine pathogens have a long history of zoonotic transmission to humans, occasionally leading to sustained outbreaks or pandemics. Through a retrospective epidemiological study of swine populations in China, we describe novel lineages of porcine hemagglutinating encephalomyelitis virus (PHEV) complex coronaviruses (CoVs) that cause exclusively respiratory symptoms with no signs of the neurological symptoms typically associated with classical PHEV infection. Through large-scale epidemiological surveillance, we show that these novel lineages have circulated in at least eight provinces in southeastern China. Phylogenetic and recombination analyses of twenty-four genomes identified two major viral lineages causing respiratory symptoms with extensive recombination within them, between them, and between classical PHEV and the novel respiratory variant PHEV (rvPHEV) lineages. Divergence times among the sampled lineages in the PHEV virus complex date back to 1886-1958 (mean estimate 1928), with the two major rvPHEV lineages separating approximately 20 years later. Many rvPHEV viruses show amino acid substitutions at the carbohydrate-binding site of hemagglutinin esterase (HE) and/or have lost the cysteine required for HE dimerization. This resembles the early adaptation of human CoVs, where HE lost its hemagglutination ability to adapt to growth in the human respiratory tract. Our study represents the first report of the evolutionary history of rvPHEV circulating in swine and highlights the importance of characterizing CoV diversity and recombination in swine to identify pathogens with outbreak potential that could threaten swine farming., Competing Interests: None declared., (© The Author(s) 2023. Published by Oxford University Press.)

Increasing levels of artemisinin and partner drug resistance threaten malaria control and elimination globally. Triple artemisinin-based combination therapies (TACTs) which combine artemisinin derivatives with two partner drugs are efficacious and well tolerated in clinical trials, including in areas of multidrug-resistant malaria. Whether early TACT adoption could delay the emergence and spread of antimalarial drug resistance is a question of vital importance. Using two independent individual-based models of Plasmodium falciparum epidemiology and evolution, we evaluated whether introduction of either artesunate-mefloquine-piperaquine or artemether-lumefantrine-amodiaquine resulted in lower long-term artemisinin-resistance levels and treatment failure rates compared with continued ACT use. We show that introduction of TACTs could significantly delay the emergence and spread of artemisinin resistance and treatment failure, extending the useful therapeutic life of current antimalarial drugs, and improving the chances of malaria elimination. We conclude that immediate introduction of TACTs should be considered by policy makers in areas of emerging artemisinin resistance., (© 2023. The Author(s).)

Mass drug administration (MDA) with antimalarials has been shown to reduce prevalence and interrupt transmission in small populations, in populations with reliable access to antimalarial drugs, and in populations where sustained improvements in diagnosis and treatment are possible. In addition, when MDA is effective it eliminates both drug-resistant parasites and drug-sensitive parasites, which has the long-term benefit of extending the useful therapeutic life of first-line therapies for all populations, not just the focal population where MDA was carried out. However, in order to plan elimination measures effectively, it is necessary to characterize the conditions under which failed MDA could exacerbate resistance. We use an individual-based stochastic model of Plasmodium falciparum transmission to evaluate this risk for MDA using dihydroartemisinin-piperaquine (DHA-PPQ), in populations where access to antimalarial treatments may not be uniformly high and where re-importation of drug-resistant parasites may be common. We find that artemisinin-resistance evolution at the kelch13 locus can be accelerated by MDA when all three of the following conditions are met: (1) strong genetic bottlenecking that falls short of elimination, (2) re-importation of artemisinin-resistant genotypes, and (3) continued selection pressure during routine case management post-MDA. Accelerated resistance levels are not immediate but follow the rebound of malaria cases post-MDA, if this is allowed to occur. Crucially, resistance is driven by the selection pressure during routine case management post-MDA and not the selection pressure exerted during the MDA itself. Second, we find that increasing treatment coverage post-MDA increases the probability of local elimination in low-transmission regions (prevalence < 2%) in scenarios with both low and high levels of drug-resistance importation. This emphasizes the importance of planning for and supporting high coverage of diagnosis and treatment post-MDA., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Nguyen et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Much of the world experiences influenza in yearly recurring seasons, particularly in temperate areas. These patterns can be considered repeatable if they occur predictably and consistently at the same time of year. In tropical areas, including southeast Asia, timing of influenza epidemics is less consistent, leading to a lack of consensus regarding whether influenza is repeatable. This study aimed to assess repeatability of influenza in Vietnam, with repeatability defined as seasonality that occurs at a consistent time of year with low variation. We developed a mathematical model incorporating parameters to represent periods of increased transmission and then fitted the model to data collected from sentinel hospitals throughout Vietnam as well as four temperate locations. We fitted the model for individual (sub)types of influenza as well as all combined influenza throughout northern, central, and southern Vietnam. Repeatability was evaluated through the variance of the timings of peak transmission. Model fits from Vietnam show high variance (sd = 64-179 days) in peak transmission timing, with peaks occurring at irregular intervals and throughout different times of year. Fits from temperate locations showed regular, annual epidemics in winter months, with low variance in peak timings (sd = 32-57 days). This suggests that influenza patterns are not repeatable or seasonal in Vietnam. Influenza prevention in Vietnam therefore cannot rely on anticipation of regularly occurring outbreaks., Competing Interests: The authors declare that they have no competing interests., (Copyright: © 2023 Servadio et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Nowadays, the ability to make data-driven decisions in public health is of utmost importance. To achieve this, it is necessary for modelers to comprehend the impact of models on the future state of healthcare systems. Compartmental models are a valuable tool for making informed epidemiological decisions, and the proper parameterization of these models is crucial for analyzing epidemiological events. This work evaluated the use of compartmental models in conjunction with Particle Swarm Optimization (PSO) to determine optimal solutions and understand the dynamics of Dengue epidemics. The focus was on calculating and evaluating the rate of case reproduction, R 0 , for the Republic of Panama. Three compartmental models were compared: Susceptible-Infected-Recovered (SIR), Susceptible-Exposed-Infected-Recovered (SEIR), and Susceptible-Infected-Recovered Human-Susceptible-Infected Vector (SIR Human-SI Vector, SIR-SI). The models were informed by demographic data and Dengue incidence in the Republic of Panama between 1999 and 2022, and the susceptible population was analyzed. The SIR, SEIR, and SIR-SI models successfully provided R 0 estimates ranging from 1.09 to 1.74. This study provides, to the best of our understanding, the first calculation of R 0 for Dengue outbreaks in the Republic of Panama., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2023 The Author(s).)

Background: Ongoing tetanus cases and sporadic outbreaks of vaccine-preventable diseases associated with routine vaccination programmes remain problems in many low and middle-income countries, including Vietnam. With no human-to-human transmission or natural immunity, tetanus antibody levels indicate both individual risk of tetanus and gaps in vaccination programmes., Methods: To investigate gaps in immunity to tetanus in Vietnam, a country with a historically high level of tetanus vaccination coverage, tetanus antibodies were measure by ELISA from samples selected from a long-term serum bank, established for the purposes of general-population seroepidemiological investigations in southern Vietnam. Samples were selected from 10 provinces, focussing on age-groups targeted by national vaccination programmes for infants and pregnant women (Expanded Programme on Immunization, EPI, and Maternal and Neonatal Tetanus, MNT)., Results: Antibodies were measured from a total of 3864 samples. Highest tetanus antibody concentrations occurred in children under 4 years old, over 90 % of whom had protective levels. Approximately 70 % of children aged 7-12 years had protective antibody concentrations although there was variation among provinces. For infants and children, there were no significant differences in tetanus protection between males and females, but for adults aged 20-35 years, in five of the ten provinces surveyed, protection against tetanus was higher in females (p < 0.05) who are eligible for booster doses under the MNT programme. In seven of ten provinces, antibody concentrations were inversely related to age (p < 0.01) and protection of older individuals was generally low., Conclusion: Widespread immunity to tetanus toxoid is seen in infants and young children consistent with the high coverage rates reported for diptheria tetanus toxoid and pertussis (DTP) in Vietnam. However, the lower antibody concentrations seen in older children and men suggest reduced immunity to tetanus in populations not targeted by EPI and MNT programmes., Competing Interests: Declaration of Competing Interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Louise Thwaites reports financial support was provided by Wellcome Trust. OUCRU affiliated authors reports financial support was provided by Wellcome Trust., (Copyright © 2023 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Disease emergence in livestock is a product of environment, epidemiology and economic forces. The environmental factors contributing to novel pathogen emergence in humans have been studied extensively, but the two-way relationship between farm microeconomics and outbreak risk has received comparably little attention. We introduce a game-theoretic model where farmers produce and sell two goods, one of which (e.g. pigs, poultry) is susceptible to infection by a pathogen. We model market and epidemiological effects at both the individual farm level and the community level. We find that in the case of low demand elasticity for livestock meat, the presence of an animal pathogen causing production losses can lead to a bistable system where two outcomes are possible: (i) successful disease control or (ii) maintained disease circulation, where farmers slaughter their animals at a low rate, face substantial production losses, but maintain large herds because of the appeal of high meat prices. Our observations point to the potentially critical effect of price elasticity of demand for livestock products on the success or failure of livestock disease control policies. We show the potential epidemiological benefits of (i) policies aimed at stabilizing livestock product prices, (ii) subsidies for alternative agricultural activities during epidemics, and (iii) diversifying agricultural production and sources of proteins available to consumers., (© 2023 The Authors.)

We describe the MalariaGEN Pf7 data resource, the seventh release of Plasmodium falciparum genome variation data from the MalariaGEN network.  It comprises over 20,000 samples from 82 partner studies in 33 countries, including several malaria endemic regions that were previously underrepresented.  For the first time we include dried blood spot samples that were sequenced after selective whole genome amplification, necessitating new methods to genotype copy number variations.  We identify a large number of newly emerging crt mutations in parts of Southeast Asia, and show examples of heterogeneities in patterns of drug resistance within Africa and within the Indian subcontinent.  We describe the profile of variations in the C-terminal of the csp gene and relate this to the sequence used in the RTS,S and R21 malaria vaccines.  Pf7 provides high-quality data on genotype calls for 6 million SNPs and short indels, analysis of large deletions that cause failure of rapid diagnostic tests, and systematic characterisation of six major drug resistance loci, all of which can be freely downloaded from the MalariaGEN website., Competing Interests: No competing interests were disclosed., (Copyright: © 2023 MalariaGEN et al.)

Malaria due to the Plasmodium falciparum parasite remains a threat to human health despite eradication efforts and the development of anti-malarial treatments, such as artemisinin combination therapies. Human movement and migration have been linked to the propagation of malaria on national scales, highlighting the need for the incorporation of human movement in modeling efforts. Spatially couped individual-based models have been used to study how anti-malarial resistance evolves and spreads in response to drug policy changes; however, as the spatial scale of the model increases, the challenges associated with modeling of movement also increase. In this paper we discuss the development, calibration, and validation of a movement model in the context of a national-scale, spatial, individual-based model used to study the evolution of drug resistance in the malaria parasite., (© 2023. The Author(s).)

Treatment of symptomatic malaria became a routine component of the clinical and public health response to malaria after the second world war. However, all antimalarial drugs deployed against malaria eventually generated enough drug resistance that they had to be removed from use. Chloroquine, sulfadoxine-pyrimethamine, and mefloquine are well known examples of antimalarial drugs to which resistance did and still does ready evolve. Artemisinin-based combination therapies (ACTs) are currently facing the same challenge as artemisinin resistance is widespread in Southeast Asia and emerging in Africa. Here, I review some aspects of drug-resistance management in malaria that influence the strength of selective pressure on drug-resistant malaria parasites, as well as an approach we can take in the future to avoid repeating the common mistake of deploying a new drug and waiting for drug resistance and treatment failure to arrive. A desirable goal of drug-resistance management is to reduce selection pressure without reducing the overall percentage of patients that are treated. This can be achieved by distributing multiple first-line therapies (MFT) simultaneously in the population for the treatment of uncomplicated falciparum malaria, thereby keeping treatment levels high but the overall selection pressure exerted by each individual therapy low. I review the primary reasons that make MFT a preferred resistance management option in many malaria-endemic settings, and I describe two exceptions where caution and additional analyses may be warranted before deploying MFT. MFT has shown to be feasible in practice in many endemic settings. The continual improvement and increased coverage of genomic surveillance in malaria may allow countries to implement custom MFT strategies based on their current drug-resistance profiles., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (© 2022 Boni.)

Background: Artemisinin-resistant genotypes of Plasmodium falciparum have now emerged a minimum of six times on three continents despite recommendations that all artemisinins be deployed as artemisinin combination therapies (ACTs). Widespread resistance to the non-artemisinin partner drugs in ACTs has the potential to limit the clinical and resistance benefits provided by combination therapy. We aimed to model and evaluate the long-term effects of high levels of partner-drug resistance on the early emergence of artemisinin-resistant genotypes., Methods: Using a consensus modelling approach, we used three individual-based mathematical models of Plasmodium falciparum transmission to evaluate the effects of pre-existing partner-drug resistance and ACT deployment on the evolution of artemisinin resistance. Each model simulates 100 000 individuals in a particular transmission setting (malaria prevalence of 1%, 5%, 10%, or 20%) with a daily time step that updates individuals' infection status, treatment status, immunity, genotype-specific parasite densities, and clinical state. We modelled varying access to antimalarial drugs if febrile (coverage of 20%, 40%, or 60%) with one primary ACT used as first-line therapy: dihydroartemisinin-piperaquine (DHA-PPQ), artesunate-amodiaquine (ASAQ), or artemether-lumefantrine (AL). The primary outcome was time until 0·25 580Y allele frequency for artemisinin resistance (the establishment time)., Findings: Higher frequencies of pre-existing partner-drug resistant genotypes lead to earlier establishment of artemisinin resistance. Across all models, a 10-fold increase in the frequency of partner-drug resistance genotypes on average corresponded to loss of artemisinin efficacy 2-12 years earlier. Most reductions in time to artemisinin resistance establishment were observed after an increase in frequency of the partner-drug resistance genotype from 0·0 to 0·10., Interpretation: Partner-drug resistance in ACTs facilitates the early emergence of artemisinin resistance and is a major public health concern. Higher-grade partner-drug resistance has the largest effect, with piperaquine resistance accelerating the early emergence of artemisinin-resistant alleles the most. Continued investment in molecular surveillance of partner-drug resistant genotypes to guide choice of first-line ACT is paramount., Funding: Schmidt Science Fellowship in partnership with the Rhodes Trust; Bill & Melinda Gates Foundation; Wellcome Trust., Competing Interests: Declaration of interests We declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license. Published by Elsevier Ltd.. All rights reserved.)

Importance: In emergency epidemic and pandemic settings, public health agencies need to be able to measure the population-level attack rate, defined as the total percentage of the population infected thus far. During vaccination campaigns in such settings, public health agencies need to be able to assess how much the vaccination campaign is contributing to population immunity; specifically, the proportion of vaccines being administered to individuals who are already seropositive must be estimated., Objective: To estimate population-level immunity to SARS-CoV-2 through May 31, 2021, in Rhode Island, Massachusetts, and Connecticut., Design, Setting, and Participants: This observational case series assessed cases, hospitalizations, intensive care unit occupancy, ventilator occupancy, and deaths from March 1, 2020, to May 31, 2021, in Rhode Island, Massachusetts, and Connecticut. Data were analyzed from July 2021 to November 2021., Exposures: COVID-19-positive test result reported to state department of health., Main Outcomes and Measures: The main outcomes were statistical estimates, from a bayesian inference framework, of the percentage of individuals as of May 31, 2021, who were (1) previously infected and vaccinated, (2) previously uninfected and vaccinated, and (3) previously infected but not vaccinated., Results: At the state level, there were a total of 1 160 435 confirmed COVID-19 cases in Rhode Island, Massachusetts, and Connecticut. The median age among individuals with confirmed COVID-19 was 38 years. In autumn 2020, SARS-CoV-2 population immunity (equal to the attack rate at that point) in these states was less than 15%, setting the stage for a large epidemic wave during winter 2020 to 2021. Population immunity estimates for May 31, 2021, were 73.4% (95% credible interval [CrI], 72.9%-74.1%) for Rhode Island, 64.1% (95% CrI, 64.0%-64.4%) for Connecticut, and 66.3% (95% CrI, 65.9%-66.9%) for Massachusetts, indicating that more than 33% of residents in these states were fully susceptible to infection when the Delta variant began spreading in July 2021. Despite high vaccine coverage in these states, population immunity in summer 2021 was lower than planned owing to an estimated 34.1% (95% CrI, 32.9%-35.2%) of vaccines in Rhode Island, 24.6% (95% CrI, 24.3%-25.1%) of vaccines in Connecticut, and 27.6% (95% CrI, 26.8%-28.6%) of vaccines in Massachusetts being distributed to individuals who were already seropositive., Conclusions and Relevance: These findings suggest that future emergency-setting vaccination planning may have to prioritize high vaccine coverage over optimized vaccine distribution to ensure that sufficient levels of population immunity are reached during the course of an ongoing epidemic or pandemic.

Among the 30 nonsynonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (1) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (2) interactions of Spike with ACE2 receptors, and (3) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any virus within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and, in combination with other mutations, adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron overall previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The SARS-CoV-2 epidemic in southern Africa has been characterized by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, while the second and third waves were driven by the Beta (B.1.351) and Delta (B.1.617.2) variants, respectively 1-3 . In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron, B.1.1.529) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, which are predicted to influence antibody neutralization and spike function 4 . Here we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity., (© 2022. The Author(s).)

Artemisinin combination therapies (ACTs) are the WHO-recommended first-line therapies for uncomplicated Plasmodium falciparum malaria. The emergence and spread of artemisinin-resistant genotypes is a major global public health concern due to the increased rate of treatment failures that result. This is particularly germane for WHO designated 'high burden to high impact' (HBHI) countries, such as Burkina Faso, where there is increased emphasis on improving guidance, strategy, and coordination of local malaria response in an effort to reduce the prevalence of P. falciparum malaria. To explore how the increased adoption of ACTs may affect the HBHI malaria setting of Burkina Faso, we added spatial structure to a validated individual-based stochastic model of P. falciparum transmission and evaluated the long-term effects of increased ACT use. We explored how de novo emergence of artemisinin-resistant genotypes, such as pfkelch13 580Y, may occur under scenarios in which private-market drugs are eliminated or multiple first-line therapies (MFT) are deployed. We found that elimination of private market drugs would result in lower treatment failures rates (between 11.98% and 12.90%) when compared to the status quo (13.11%). However, scenarios incorporating MFT with equal deployment of artemether-lumefantrine (AL) and dihydroartemisinin-piperaquine (DHA-PPQ) may accelerate near-term drug resistance (580Y frequency ranging between 0.62 to 0.84 in model year 2038) and treatment failure rates (26.69% to 34.00% in 2038), due to early failure and substantially reduced treatment efficacy resulting from piperaquine-resistant genotypes. A rebalanced MFT approach (90% AL, 10% DHA-PPQ) results in approximately equal long-term outcomes to using AL alone but may be difficult to implement in practice., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2022 Zupko et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

State-level reopenings in late spring 2020 facilitated the resurgence of severe acute respiratory syndrome coronavirus 2 transmission. Here, we analyze age-structured case, hospitalization, and death time series from three states-Rhode Island, Massachusetts, and Pennsylvania-that had successful reopenings in May 2020 without summer waves of infection. Using 11 daily data streams, we show that from spring to summer, the epidemic shifted from an older to a younger age profile and that elderly individuals were less able to reduce contacts during the lockdown period when compared to younger individuals. Clinical case management improved from spring to summer, resulting in fewer critical care admissions and lower infection fatality rate. Attack rate estimates through 31 August 2020 are 6.2% [95% credible interval (CI), 5.7 to 6.8%] of the total population infected for Rhode Island, 6.7% (95% CI, 5.4 to 7.6%) in Massachusetts, and 2.7% (95% CI, 2.5 to 3.1%) in Pennsylvania.

The present analysis uses the data of confirmed incidence of dengue cases in the metropolitan region of Panama from 1999 to 2017 and climatic variables (air temperature, precipitation, and relative humidity) during the same period to determine if there exists a correlation between these variables. In addition, we compare the predictive performance of two regression models (SARIMA, SARIMAX) and a recurrent neural network model (RNN-LSTM) on the dengue incidence series. For this data from 1999-2014 was used for training and the three subsequent years of incidence 2015-2017 were used for prediction. The results show a correlation coefficient between the climatic variables and the incidence of dengue were low but statistical significant. The RMSE and MAPE obtained for the SARIMAX and RNN-LSTM models were 25.76, 108.44 and 26.16, 59.68, which suggest that any of these models can be used to predict new outbreaks. Although, it can be said that there is a limited role of climatic variables in the outputs the models. The value of this work is that it helps understand the behaviour of cases in a tropical setting as is the Metropolitan Region of Panama City, and provides the basis needed for a much needed early alert system for the region.

The relationship between age and seroprevalence can be used to estimate the annual attack rate of an infectious disease. For pathogens with multiple serologically distinct strains, there is a need to describe composite exposure to an antigenically variable group of pathogens. In this study, we assay 24,402 general-population serum samples, collected in Vietnam between 2009 to 2015, for antibodies to eleven human influenza A strains. We report that a principal components decomposition of antibody titer data gives the first principal component as an appropriate surrogate for seroprevalence; this results in annual attack rate estimates of 25.6% (95% CI: 24.1% - 27.1%) for subtype H3 and 16.0% (95% CI: 14.7% - 17.3%) for subtype H1. The remaining principal components separate the strains by serological similarity and associate birth cohorts with their particular influenza histories. Our work shows that dimensionality reduction can be used on human antibody profiles to construct an age-seroprevalence relationship for antigenically variable pathogens., (© 2021. The Author(s).)

We present evidence for multiple independent origins of recombinant SARS-CoV-2 viruses sampled from late 2020 and early 2021 in the United Kingdom. Their genomes carry single-nucleotide polymorphisms and deletions that are characteristic of the B.1.1.7 variant of concern but lack the full complement of lineage-defining mutations. Instead, the remainder of their genomes share contiguous genetic variation with non-B.1.1.7 viruses circulating in the same geographic area at the same time as the recombinants. In four instances, there was evidence for onward transmission of a recombinant-origin virus, including one transmission cluster of 45 sequenced cases over the course of 2 months. The inferred genomic locations of recombination breakpoints suggest that every community-transmitted recombinant virus inherited its spike region from a B.1.1.7 parental virus, consistent with a transmission advantage for B.1.1.7's set of mutations., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Since the first reports of a novel severe acute respiratory syndrome (SARS)-like coronavirus in December 2019 in Wuhan, China, there has been intense interest in understanding how severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in the human population. Recent debate has coalesced around two competing ideas: a "laboratory escape" scenario and zoonotic emergence. Here, we critically review the current scientific evidence that may help clarify the origin of SARS-CoV-2., Competing Interests: Declaration of interests E.C.H. is an honorary visiting professor at Fudan University (Shanghai Public Health Clinical Center), Shanghai, China and, from 2014–2020, was a guest professor at the Chinese Center for Disease Control and Prevention, Beijing, China. These affiliations are only used in papers co-authored with Prof. Yong-Zhen Zhang (Shanghai Public Health Clinical Center) and involve no formal appointment, no duties, and no remuneration nor research funding. J.O.W. receives funding from the U.S. Centers for Disease Control and Prevention (ongoing) via grants and contracts to his institution unrelated to this research. S.R.W. consults for Immunome and Ocugen. A.R., A.L.R., M.F.B., S.A.G., and K.G.A. have received consulting fees and compensation for expert testimony on SARS-CoV-2 and the COVID-19 pandemic. R.F.G. is co-founder of Zalgen Labs., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: When three SARS-CoV-2 vaccines came to market in Europe and North America in the winter of 2020-2021, distribution networks were in a race against a major epidemiological wave of SARS-CoV-2 that began in autumn 2020. Rapid and optimized vaccine allocation was critical during this time. With 95% efficacy reported for two of the vaccines, near-term public health needs likely require that distribution is prioritized to the elderly, health care workers, teachers, essential workers, and individuals with comorbidities putting them at risk of severe clinical progression., Methods: We evaluate various age-based vaccine distributions using a validated mathematical model based on current epidemic trends in Rhode Island and Massachusetts. We allow for varying waning efficacy of vaccine-induced immunity, as this has not yet been measured. We account for the fact that known COVID-positive cases may not have been included in the first round of vaccination. And, we account for age-specific immune patterns in both states at the time of the start of the vaccination program. Our analysis assumes that health systems during winter 2020-2021 had equal staffing and capacity to previous phases of the SARS-CoV-2 epidemic; we do not consider the effects of understaffed hospitals or unvaccinated medical staff., Results: We find that allocating a substantial proportion (>75%) of vaccine supply to individuals over the age of 70 is optimal in terms of reducing total cumulative deaths through mid-2021. This result is robust to different profiles of waning vaccine efficacy and several different assumptions on age mixing during and after lockdown periods. As we do not explicitly model other high-mortality groups, our results on vaccine allocation apply to all groups at high risk of mortality if infected. A median of 327 to 340 deaths can be avoided in Rhode Island (3444 to 3647 in Massachusetts) by optimizing vaccine allocation and vaccinating the elderly first. The vaccination campaigns are expected to save a median of 639 to 664 lives in Rhode Island and 6278 to 6618 lives in Massachusetts in the first half of 2021 when compared to a scenario with no vaccine. A policy of vaccinating only seronegative individuals avoids redundancy in vaccine use on individuals that may already be immune, and would result in 0.5% to 1% reductions in cumulative hospitalizations and deaths by mid-2021., Conclusions: Assuming high vaccination coverage (>28%) and no major changes in distancing, masking, gathering size, hygiene guidelines, and virus transmissibility between 1 January 2021 and 1 July 2021 a combination of vaccination and population immunity may lead to low or near-zero transmission levels by the second quarter of 2021., (© 2021. The Author(s).)

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination., Competing Interests: No competing interests were disclosed., (Copyright: © 2021 MalariaGEN et al.)

Tetanus arises from wound contamination with Clostridium tetani, but approximately one fifth of patients have no discernable entry wound. Clostridium tetani is culturable from animal feces, suggesting the gastrointestinal tract could be an endogenous reservoir or direct-entry portal, but human data are lacking. In this study of 101 Vietnamese adults with tetanus and 29 hospitalized control subjects, admission stool samples were cultured for C. tetani. Anti-tetanus toxin antibodies were measured by ELISA. Clostridium tetani toxigenicity was evaluated using polymerase chain reaction and sequencing. Toxigenic C. tetani was cultured from stool samples in 50 of 100 (50%) tetanus cases and 12 of 28 (42.9%) control subjects (P = 0.50), and stool samples of 44 of 85 (52.4%) tetanus cases with clinically identified wounds compared with 6 of 15 (47.6%) patients without clinically identified wounds (P = 0.28). Nine of 12 (75%) control subjects with toxigenic C. tetani in their stool samples lacked protective antibody concentrations. These findings fail to show evidence of an association between gastrointestinal C. tetani and tetanus infection, but emphasize the importance of increasing vaccination coverage.

Background: In recent years, researchers have had an increased focus on multiplex microarray assays, in which antibodies are measured against multiple related antigens, for use in seroepidemiological studies to infer past transmission., Methods: We assess the performance of a flavivirus microarray assay for determining past dengue virus (DENV) infection history in a dengue-endemic setting, Vietnam. We tested the microarray on samples from 1 and 6 months postinfection from DENV-infected patients (infecting serotype was determined using reverse-transcription polymerase chain reaction during acute, past primary, and secondary infection assessed using plaque reduction neutralization tests 6 months postinfection)., Results: Binomial models developed to discriminate past primary from secondary infection using the protein microarray (PMA) titers had high area under the curve (0.90-0.97) and accuracy (0.84-0.86). Multinomial models developed to identify most recent past infecting serotype using PMA titers performed well in those with past primary infection (average test set: κ = 0.85, accuracy of 0.92) but not those with past secondary infection (κ = 0.24, accuracy of 0.45)., Conclusions: Our results suggest that the microarray will be useful in seroepidemiological studies aimed at classifying the past infection history of individuals (past primary vs secondary and serotype of past primary infections) and thus inferring past transmission intensity of DENV in dengue-endemic settings. Future work to validate these models should be undertaken in different transmission settings and with samples later after infection., (© The Author(s) 2020. Published by Oxford University Press for the Infectious Diseases Society of America.)

Virus host shifts are generally associated with novel adaptations to exploit the cells of the new host species optimally. Surprisingly, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has apparently required little to no significant adaptation to humans since the start of the Coronavirus Disease 2019 (COVID-19) pandemic and to October 2020. Here we assess the types of natural selection taking place in Sarbecoviruses in horseshoe bats versus the early SARS-CoV-2 evolution in humans. While there is moderate evidence of diversifying positive selection in SARS-CoV-2 in humans, it is limited to the early phase of the pandemic, and purifying selection is much weaker in SARS-CoV-2 than in related bat Sarbecoviruses. In contrast, our analysis detects evidence for significant positive episodic diversifying selection acting at the base of the bat virus lineage SARS-CoV-2 emerged from, accompanied by an adaptive depletion in CpG composition presumed to be linked to the action of antiviral mechanisms in these ancestral bat hosts. The closest bat virus to SARS-CoV-2, RmYN02 (sharing an ancestor about 1976), is a recombinant with a structure that includes differential CpG content in Spike; clear evidence of coinfection and evolution in bats without involvement of other species. While an undiscovered "facilitating" intermediate species cannot be discounted, collectively, our results support the progenitor of SARS-CoV-2 being capable of efficient human-human transmission as a consequence of its adaptive evolutionary history in bats, not humans, which created a relatively generalist virus., Competing Interests: The authors have declared that no competing interests exist.

There are outstanding evolutionary questions on the recent emergence of human coronavirus SARS-CoV-2 including the role of reservoir species, the role of recombination and its time of divergence from animal viruses. We find that the sarbecoviruses-the viral subgenus containing SARS-CoV and SARS-CoV-2-undergo frequent recombination and exhibit spatially structured genetic diversity on a regional scale in China. SARS-CoV-2 itself is not a recombinant of any sarbecoviruses detected to date, and its receptor-binding motif, important for specificity to human ACE2 receptors, appears to be an ancestral trait shared with bat viruses and not one acquired recently via recombination. To employ phylogenetic dating methods, recombinant regions of a 68-genome sarbecovirus alignment were removed with three independent methods. Bayesian evolutionary rate and divergence date estimates were shown to be consistent for these three approaches and for two different prior specifications of evolutionary rates based on HCoV-OC43 and MERS-CoV. Divergence dates between SARS-CoV-2 and the bat sarbecovirus reservoir were estimated as 1948 (95% highest posterior density (HPD): 1879-1999), 1969 (95% HPD: 1930-2000) and 1982 (95% HPD: 1948-2009), indicating that the lineage giving rise to SARS-CoV-2 has been circulating unnoticed in bats for decades.

Due to the scope and impact of the COVID-19 pandemic there exists a strong desire to understand where the SARS-CoV-2 virus came from and how it jumped species boundaries to humans. Molecular evolutionary analyses can trace viral origins by establishing relatedness and divergence times of viruses and identifying past selective pressures. However, we must uphold rigorous standards of inference and interpretation on this topic because of the ramifications of being wrong. Here, we dispute the conclusions of Xia (2020. Extreme genomic CpG deficiency in SARS-CoV-2 and evasion of host antiviral defense. Mol Biol Evol. doi:10.1093/molbev/masa095) that dogs are a likely intermediate host of a SARS-CoV-2 ancestor. We highlight major flaws in Xia's inference process and his analysis of CpG deficiencies, and conclude that there is no direct evidence for the role of dogs as intermediate hosts. Bats and pangolins currently have the greatest support as ancestral hosts of SARS-CoV-2, with the strong caveat that sampling of wildlife species for coronaviruses has been limited., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)