The origins of SARS-CoV-2: A critical review.

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.

Generation and transmission of interlineage recombinants in the SARS-CoV-2 pandemic.

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.

Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa.

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, respectively1-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 function4. Here we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.

Natural selection in the evolution of SARS-CoV-2 in bats created a generalist virus and highly capable human pathogen.

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.

Repeatability and timing of tropical influenza epidemics.

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.

Selection Analysis Identifies Clusters of Unusual Mutational Changes in Omicron Lineage BA.1 That Likely Impact Spike Function.

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.

Benefits of near-universal vaccination and treatment access to manage COVID-19 burden in the United States.

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.

Using NS1 Flavivirus Protein Microarray to Infer Past Infecting Dengue Virus Serotype and Number of Past Dengue Virus Infections in Vietnamese Individuals.

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.

Viral CpG Deficiency Provides No Evidence That Dogs Were Intermediate Hosts for SARS-CoV-2.

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.

Optimal SARS-CoV-2 vaccine allocation using real-time attack-rate estimates in Rhode Island and Massachusetts.

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.

Improved Algorithmic Complexity for the 3SEQ Recombination Detection Algorithm.

Identifying recombinant sequences in an era of large genomic databases is challenging as it requires an efficient algorithm to identify candidate recombinants and parents, as well as appropriate statistical methods to correct for the large number of comparisons performed. In 2007, a computation was introduced for an exact nonparametric mosaicism statistic that gave high-precision P values for putative recombinants. This exact computation meant that multiple-comparisons corrected P values also had high precision, which is crucial when performing millions or billions of tests in large databases. Here, we introduce an improvement to the algorithmic complexity of this computation from O(mn3) to O(mn2), where m and n are the numbers of recombination-informative sites in the candidate recombinant. This new computation allows for recombination analysis to be performed in alignments with thousands of polymorphic sites. Benchmark runs are presented on viral genome sequence alignments, new features are introduced, and applications outside recombination analysis are discussed.

Multilocus sequence typing of Cryptococcus neoformans var. grubii from Laos in a regional and global context.

Cryptococcosis causes approximately 180 000 deaths each year in patients with human immunodeficiency virus (HIV). Patients with other forms of immunosuppression are also at risk, and disease is increasingly recognized in apparently immunocompetent individuals. Cryptococcus neoformans var. grubii, responsible for the majority of cases, is distributed globally. We used the consensus ISHAM Multilocus sequence typing (MLST) scheme to define the population structure of clinical C. neoformans var. grubii isolates from Laos (n = 81), which we placed into the global context using published MLST data from other countries (total N = 1047), including a reanalysis of 136 Vietnamese isolates previously reported. We observed a phylogeographical relationship in which the Laotian population was similar to its neighbor Thailand, being dominated (83%) by Sequence Types (ST) 4 and 6. This phylogeographical structure changed moving eastwards, with Vietnam's population consisting of an admixture of isolates dominated by the ST4/ST6 (35%) and ST5 (48%) lineages. The ST5 lineage is the predominant ST reported from China and East Asia, where it accounts for >90% of isolates. Analysis of genetic distance (Fst) between different populations of C. neoformans var. grubii supports this intermediate structure of the Vietnamese population. The pathogen and host diversity reported from Vietnam provide the strongest epidemiological evidence of the association between ST5 and HIV-uninfected patients. Regional anthropological genetic distances suggest diversity in the C. neoformans var. grubii population across Southeast Asia is driven by ecological rather than human host factors. Where the ST5 lineage is present, disease in HIV-uninfected patients is to be expected.

Poultry population dynamics and mortality risks in smallholder farms of the Mekong river delta region.

BACKGROUND: Poultry farming is widely practiced by rural households in Vietnam and the vast majority of domestic birds are kept on small household farms. However, smallholder poultry production is constrained by several issues such as infectious diseases, including avian influenza viruses whose circulation remains a threat to public health. This observational study describes the demographic structure and dynamics of small-scale poultry farms of the Mekong river delta region. METHOD: Fifty three farms were monitored over a 20-month period, with farm sizes, species, age, arrival/departure of poultry, and farm management practices recorded monthly. RESULTS: Median flock population sizes were 16 for chickens (IQR: 10-40), 32 for ducks (IQR: 18-101) and 11 for Muscovy ducks (IQR: 7-18); farm size distributions for the three species were heavily right-skewed. Muscovy ducks were kept for long periods and outdoors, while chickens and ducks were farmed indoors or in pens. Ducks had a markedly higher removal rate (broilers: 0.14/week; layer/breeders: 0.05/week) than chickens and Muscovy ducks (broilers: 0.07/week; layer/breeders: 0.01-0.02/week) and a higher degree of specialization resulting in a substantially shorter life span. The rate of mortality due to disease did not differ much among species, with birds being less likely to die from disease at older ages, but frequency of disease symptoms differed by species. Time series of disease-associated mortality were correlated with population size for Muscovy ducks (Kendall's coefficient τ = 0.49, p-value < 0.01) and with frequency of outdoor grazing for ducks (τ = 0.33, p-value = 0.05). CONCLUSION: The study highlights some challenges to disease control in small-scale multispecies poultry farms. The rate of interspecific contact and overlap between flocks of different ages is high, making small-scale farms a suitable environment for pathogens circulation. Muscovy ducks are farmed outdoors with little investment in biosecurity and few inter-farm movements. Ducks and chickens are more at-risk of introduction of pathogens through movements of birds from one farm to another. Ducks are farmed in large flocks with high turnover and, as a result, are more vulnerable to disease spread and require a higher vaccination coverage to maintain herd immunity.

The decline of malaria in Vietnam, 1991-2014.

BACKGROUND: Despite the well-documented clinical efficacy of artemisinin-based combination therapy (ACT) against malaria, the population-level effects of ACT have not been studied thoroughly until recently. An ideal case study for these population-level effects can be found in Vietnam's gradual adoption of artemisinin in the 1990s. METHODS AND RESULTS: Analysis of Vietnam's national annual malaria reports (1991-2014) revealed that a 10% increase in artemisinin procurement corresponded to a 32.8% (95% CI 27.7-37.5%) decline in estimated malaria cases. There was no consistent national or regional effect of vector control on malaria. The association between urbanization and malaria was generally negative and sometimes statistically significant. CONCLUSIONS: The decline of malaria in Vietnam can largely be attributed to the adoption of artemisinin-based case management. Recent analyses from Africa showed that insecticide-treated nets had the greatest effect on lowering malaria prevalence, suggesting that the success of interventions is region-specific. Continuing malaria elimination efforts should focus on both vector control and increased access to ACT.

Modeling tuberculosis dynamics with the presence of hyper-susceptible individuals for Ho Chi Minh City from 1996 to 2015.

BACKGROUND: The depletion of CD4 cell is the underlying reason for TB hyper-susceptibility among people with HIV. Consequently, the trend of TB dynamics is usually hidden by the HIV outbreak. METHODS: Here, we aim to evaluate the trend of TB dynamics quantitatively by a simple mathematical model using the known prevalence of hyper-susceptible individuals in the population. In order to estimate the parameters governing transmission we fit this model in a maximum likelihood framework to both reported TB cases and data from samples tested with Interferon Gamma Assay from Ho Chi Minh City - a city with high TB transmission and strong synchronization between HIV/AIDS and TB dynamics. RESULTS: Our results show that TB transmission in HCMC has been declining among people without HIV; we estimate a 18% (95% CI: 9-25%) decline in the transmission parameter between 1996 and 2015. Furthermore, we show that co-infected patients have limited contribution to the transmission process. For hyper-susceptible individuals, our model suggests that the risk of a new active TB infection occurring is significantly higher than the risk of relapsed active TB, while this is not the case for people without hyper-susceptibility. CONCLUSIONS: The increase of TB notifications in Ho Chi Minh City from 1996 to 2008 is evitable when, as occurred, the number of hyper-susceptible individuals increased faster than the decrease of TB transmission rate. The sharp decrease in TB notifications observed in this city from 2008 to 2015 is the combined result of the decrease of TB transmission rate and the decrease of hyper-susceptible individuals in the population. For hyper-susceptible individuals, we propose that the reason for the reduced relapsed active TB risk is HIV treatment delay. According to HIV treatment guidelines issued by Vietnam's Ministry of Health, hyper-susceptible individuals usually have to wait until their CD4 cell count falls under 350 cells/µl to start ART. Once patients begin ART, they will remain on ART for the rest of their life and thus have greater protection against relapses of TB. We therefore hypothesize that the delay in using ART imposes considerable TB burden on HCMC despite the declining transmission process.

Impact of human mobility on the emergence of dengue epidemics in Pakistan.

The recent emergence of dengue viruses into new susceptible human populations throughout Asia and the Middle East, driven in part by human travel on both local and global scales, represents a significant global health risk, particularly in areas with changing climatic suitability for the mosquito vector. In Pakistan, dengue has been endemic for decades in the southern port city of Karachi, but large epidemics in the northeast have emerged only since 2011. Pakistan is therefore representative of many countries on the verge of countrywide endemic dengue transmission, where prevention, surveillance, and preparedness are key priorities in previously dengue-free regions. We analyze spatially explicit dengue case data from a large outbreak in Pakistan in 2013 and compare the dynamics of the epidemic to an epidemiological model of dengue virus transmission based on climate and mobility data from ∼40 million mobile phone subscribers. We find that mobile phone-based mobility estimates predict the geographic spread and timing of epidemics in both recently epidemic and emerging locations. We combine transmission suitability maps with estimates of seasonal dengue virus importation to generate fine-scale dynamic risk maps with direct application to dengue containment and epidemic preparedness.

K13 Propeller Mutations in Plasmodium falciparum Populations in Regions of Malaria Endemicity in Vietnam from 2009 to 2016.

The spread of artemisinin-resistant Plasmodium falciparum compromises the therapeutic efficacy of artemisinin combination therapies (ACTs) and is considered the greatest threat to current global initiatives to control and eliminate malaria. This is particularly relevant in Vietnam, where dihydroartemisinin-piperaquine (DP) is the recommended ACT for P. falciparum infection. The propeller domain gene of K13, a molecular marker of artemisinin resistance, was successfully sequenced in 1,060 P. falciparum isolates collected at 3 malaria hot spots in Vietnam between 2009 and 2016. Eight K13 propeller mutations (Thr474Ile, Tyr493His, Arg539Thr, Ile543Thr, Pro553Leu, Val568Gly, Pro574Leu, and Cys580Tyr), including several that have been validated to be artemisinin resistance markers, were found. The prevalences of K13 mutations were 29% (222/767), 6% (11/188), and 43% (45/105) in the Binh Phuoc, Ninh Thuan, and Gia Lai Provinces of Vietnam, respectively. Cys580Tyr became the dominant genotype in recent years, with 79.1% (34/43) of isolates in Binh Phuoc Province and 63% (17/27) of isolates in Gia Lai Province carrying this mutation. K13 mutations were associated with reduced ring-stage susceptibility to dihydroartemisinin (DHA) in vitro and prolonged parasite clearance in vivo An analysis of haplotypes flanking K13 suggested the presence of multiple strains with the Cys580Tyr mutation rather than a single strain expanding across the three sites.

Towards malaria elimination in Savannakhet, Lao PDR: mathematical modelling driven strategy design.

BACKGROUND: The number of Plasmodium falciparum malaria cases around the world has decreased substantially over the last 15 years, but with the spread of resistance against anti-malarial drugs and insecticides, this decline may not continue. There is an urgent need to consider alternative, accelerated strategies to eliminate malaria in countries like Lao PDR, where there are a few remaining endemic areas. A deterministic compartmental modelling tool was used to develop an integrated strategy for P. falciparum elimination in the Savannakhet province of Lao PDR. The model was designed to include key aspects of malaria transmission and integrated control measures, along with a user-friendly interface. RESULTS: Universal coverage was the foundation of the integrated strategy, which took the form of the deployment of community health workers who provided universal access to early diagnosis, treatment and long-lasting insecticidal nets. Acceleration was included as the deployment of three monthly rounds of mass drug administration targeted towards high prevalence villages, with the addition of three monthly doses of the RTS,S vaccine delivered en masse to the same high prevalence sub-population. A booster dose of vaccine was added 1 year later. The surveillance-as-intervention component of the package involved the screening and treatment of individuals entering the simulated population. CONCLUSIONS: In this modelling approach, the sequential introduction of a series of five available interventions in an integrated strategy was predicted to be sufficient to stop malaria transmission within a 3-year period. These interventions comprised universal access to early diagnosis and adequate treatment, improved access to long-lasting insecticidal nets, three monthly rounds of mass drug administration together with RTS,S vaccination followed by a booster dose of vaccine, and screening and treatment of imported cases.

The Community As the Patient in Malaria-Endemic Areas: Preempting Drug Resistance with Multiple First-Line Therapies.

Maciej F. Boni and colleagues propose deploying multiple first-line combination therapies against malaria within a community to delay drug-resistance evolution.