Molecular phylodynamics of fowl adenovirus serotype 11 and 8b from inclusion body hepatitis outbreaks.

Fowl adenovirus (FAdV) serotypes are involved in a variety of clinical manifestations in poultry and has resulted in substantial economic loss to the poultry farmers. Despite the endemicity of Inclusion body hepatitis (IBH) in South Asian countries, including India, its etiology is not well studied. In western India, the rural poultry flocks obtained from the vaccinated parents were experiencing disease outbreaks with substantial economic losses due to heavy outbreaks and mortality. Therefore, the study was conducted to decipher the molecular epidemiology of the FAdV from field outbreaks in western India. A total of 37 commercial broiler poultry flocks and 29 village poultry flocks of western India were visited during 2019 to 2021. Out of these, 19.14% flocks showed incidence of IBH during the age of 15 to 35 days. The mortality ranged from 3.3 percent to 55.28 percent. The samples were subjected for amplification of partial hexon gene covering loop 1 and loop 2. The results revealed 48.28% positivity by PCR. The sequence analysis identified 14 isolates as species D serotype 11 with 0.97 to 0.99% divergence and two as species E serotype 8b with 0.99% divergence. The FAdV-11 isolates showed amino acid substitutions D195N, T399A, N417S, and N496H. The amino acids I188 and N195 were conserved in FAdV-11. The molecular clock in Bayesian methods was used to determine most common ancestor. The isolates MH379249 and MH379248 were determined the most recent common ancestor for FAdV-11 and FAdV-8b isolates. The analysis suggested evolution of 10 FAdV-11 strains in 2012, and four FAdV-11 strains and two FAdV-8b strains in 2018.

Probability Distribution of Tree Age for the Simple Birth-Death Process, with Applications to Distributions of Number of Ancestral Lineages and Divergence Times for Pairs of Taxa in a Yule Tree.

In this contribution, a general expression is derived for the probability density of the time to the most recent common ancestor (TMRCA) of a simple birth-death tree, a widely used stochastic null-model of biological speciation and extinction, conditioned on the constant birth and death rates and number of extant lineages. This density is contrasted with a previous result which was obtained using a uniform prior for the time of origin. The new distribution is applied to two problems of phylogenetic interest. First, that of the probability of the number of taxa existing at any time in the past in a tree of a known number of extant species, and given birth and death rates, and second, that of determining the TMRCA of two randomly selected taxa in an unobserved tree that is produced by a simple birth-only, or Yule, process. In the latter case, it is assumed that only the rate of bifurcation (speciation) and the size, or number of tips, are known. This is shown to lead to a closed-form analytical expression for the probability distribution of this parameter, which is arrived at based on the known mathematical form of the age distribution of Yule trees of a given size and branching rate, which is derived here de novo, and a similar distribution which additionally is conditioned on tree age. The new distribution is the exact Yule prior for divergence times of pairs of taxa under the stated conditions and is potentially useful in statistical (Bayesian) inference studies of phylogenies.

Dates and Rates of Tick-Borne Encephalitis Virus-The Slowest Changing Tick-Borne Flavivirus.

We evaluated the temporal signal and substitution rate of tick-borne encephalitis virus (TBEV) using 276 complete open reading frame (ORF) sequences with known collection dates. According to a permutation test, the TBEV Siberian subtype (TBEV-S) data set has no temporal structure and cannot be applied for substitution rate estimation without other TBEV subtypes. The substitution rate obtained suggests that the common clade of TBEV (TBEV-common), including all TBEV subtypes and louping-ill virus (LIV), is characterized by the lowest rate (1.87 × 10-5 substitutions per site per year (s/s/y) or 1 nucleotide substitution per ORF per 4.9 years; 95% highest posterior density (HPD) interval, 1.3-2.4 × 10-5 s/s/y) among all tick-borne flaviviruses previously assessed. Within TBEV-common, the TBEV European subtype (TBEV-E) has the lowest substitution rate (1.3 × 10-5 s/s/y or 1 nucleotide substitution per ORF per 7.5 years; 95% HPD, 1.0-1.8 × 10-5 s/s/y) as compared with TBEV Far-Eastern subtype (3.0 × 10-5 s/s/y or 1 nucleotide substitution per ORF per 3.2 years; 95% HPD, 1.6-4.5 × 10-5 s/s/y). TBEV-common representing the species tick-borne encephalitis virus diverged 9623 years ago (95% HPD interval, 6373-13,208 years). The TBEV Baikalian subtype is the youngest one (489 years; 95% HPD, 291-697 years) which differs significantly by age from TBEV-E (848 years; 95% HPD, 596-1112 years), LIV (2424 years; 95% HPD, 1572-3400 years), TBEV-FE (1936 years, 95% HPD, 1344-2598 years), and the joint clade of TBEV-S (2505 years, 95% HPD, 1700-3421 years) comprising Vasilchenko, Zausaev, and Baltic lineages.

eSMC: a statistical model to infer admixture events from individual genomics data.

BACKGROUND: Inferring historical population admixture events yield essential insights in understanding a species demographic history. Methods are available to infer admixture events in demographic history with extant genetic data from multiple sources. Due to the deficiency in ancient population genetic data, there lacks a method for admixture inference from a single source. Pairwise Sequentially Markovian Coalescent (PSMC) estimates the historical effective population size from lineage genomes of a single individual, based on the distribution of the most recent common ancestor between the diploid's alleles. However, PSMC does not infer the admixture event. RESULTS: Here, we proposed eSMC, an extended PSMC model for admixture inference from a single source. We evaluated our model's performance on both in silico data and real data. We simulated population admixture events at an admixture time range from 5 kya to 100 kya (5 years/generation) with population admix ratio at 1:1, 2:1, 3:1, and 4:1, respectively. The root means the square error is [Formula: see text] kya for all experiments. Then we implemented our method to infer the historical admixture events in human, donkey and goat populations. The estimated admixture time for both Han and Tibetan individuals range from 60 kya to 80 kya (25 years/generation), while the estimated admixture time for the domesticated donkeys and the goats ranged from 40 kya to 60 kya (8 years/generation) and 40 kya to 100 kya (6 years/generation), respectively. The estimated admixture times were concordance to the time that domestication occurred in human history. CONCLUSION: Our eSMC effectively infers the time of the most recent admixture event in history from a single individual's genomics data. The source code of eSMC is hosted at https://github.com/zachary-zzc/eSMC .

Analysis of the genome of grapevine red blotch virus and related grabloviruses indicates diversification prior to the arrival of Vitis vinifera in North America.

In this study 163 complete whole-genome sequences of the emerging pathogen grapevine red blotch virus (GRBV; genus Grablovirus, family Geminiviridae) were used to reconstruct phylogenies using Bayesian analyses on time-tipped (heterochronous) data. Using different combinations of priors, Bayes factors identified heterochronous datasets (3×200 million chains) generated from strict clock and exponential tree priors as being the most robust. Substitution rates of 3.2×10-5 subsitutions per site per year (95% HPD 4.3-2.1×10-5) across the whole of the GRBV genome were estimated, suggesting ancestral GRBV diverged from ancestral wild Vitis latent virus 1 around 9 000 years ago, well before the first documented arrival of Vitis vinifera in North America. Whole-genome analysis of GRBV isolates in a single infected field-grown grapevine across 12 years identified 12 single nucleotide polymorphisms none of which were fixed substitutions: an observation not discordant with the in silico estimate. The substitution rate estimated here is lower than those estimated for other geminiviruses and is the first for a woody-host-infecting geminivirus.

Molecular phylodynamics of infectious bursal disease viruses.

The present study was conducted to study the molecular phylodynamics of the Indian field IBDVs. A total of 13 organized commercial poultry farms and 29 village poultry flocks were recruited in the study. The broiler flocks showed 15.25-60.18% mortality, followed by 12.4% in improved native poultry varieties and 5% in indigenous birds. The 664 bp hypervariable VP2 gene fragment of Western and Central Indian vvIBDVs showed 97.14-98.79 and 94.49-96.69% identity to Pakistani and South Indian vvIBDVs, respectively. An isolate was 99.54% identical to the Ventri-Plus vaccine strain, while three IBDVs showed maximum identity with the Georgia strain. Out of 22, 19 strains showed typical vvIBDV amino acid signature, while three strains showed substitutions specific to classical IBDVs. Central Indian vvIBDVs showed conserved substitutions at N212D and E300A, which can be used as a regional marker. Phylogenetic genogrouping placed global IBDVs into seven genogroups based upon virulence and geographical distribution. Nineteen field vvIBDVs were placed in the G3 genogroup, and the other three were grouped with classical IBDVs in G1 genogroup. A nucleotide span from 584 to 1248 covering VP2 hypervariable fragment was found suitable for correct genogrouping of field IBDVs. The Bayesian evolutionary analysis showed tMRCA of the year 2009 for 8 Western Indian vvIBDVs with vvIBDV from Pakistan. Central Indian vvIBDVs were evolved in the year 1991 from BD-3 and PY12 strains of vvIBDVs from Bangladesh and Pondicherry, respectively. An isolate showed evolution in year 2010 from the Nigerian ABIC strain, while three classical strains showed tMRCA of the year 2009 with the Georgia strain as a recent common ancestor.

YMrCA: Improving Y-chromosomal ancestor time estimation for DNA kinship research.

The Y-chromosome is a valuable kinship indicator in family history and forensic research. To reconstruct genealogies, the time to the most recent common ancestor (tMRCA) between paternal relatives can be estimated through Y-STR analysis. Existing models are the stepwise mutation model (SMM, only one-step Y-STR changes) and the infinite allele model (IAM, new allele per Y-STR change). In this study, these mutation models and all existing tMRCA calculators were validated through a genetic-genealogy database containing 1,120 biologically related genealogical pairs confirmed by 46 Y-STRs with known tMRCA (18,109 generations). Consistent under- and overestimation and broad confidence intervals were observed, leading to dubious tMRCA estimates. This is because they do not include individual mutation rates or multi-step changes and ignore hidden multiple, back, or parallel modifications. To improve tMRCA estimation, we developed a user-friendly calculator, the "YMrCA", including all previously mentioned mutation characteristics. After extensive validation, we observed that the YMrCA calculator demonstrated a promising performance. The YMrCA yields a significantly higher tMRCA success rate (96%; +20%) and a lower tMRCA error (7; -3) compared to the mutation models and all online tMRCA calculators. Therefore, YMrCA offers the next step towards more objective tMRCA estimation for DNA kinship research.

Exploitation of a novel slowly mutating Y-STRs set and evaluation of slowly mutating Y-STRs plus Y-SNPs typing strategy in forensic genetics and evolutionary research.

The Y-chromosome short tandem repeats (Y-STRs) loci with different mutation rates existing in the Y chromosome non-recombination region (NRY) allow to be applied in human forensics, genealogical researches, historical investigations and evolutionary studies. Currently, there is a high demand for pedigree search to narrow the scope of crime investigations. However, the commonly used Y-STRs kits generally contain Y-STRs with high mutation rates that could cause individuals from the same pedigree to display different haplotypes. Herein, we put forward a new strategy of Slowly Mutating (SM) Y-STRs plus Y-SNPs typing, which could not only improve the resolution and accuracy of pedigree search, but also be applicable to evolutionary research. First, we developed a nine SM Y-STRs assay by evaluating their mutation rates in 210 pedigrees. Then the gene diversity and efficiency of the SM Y-STRs and 172 Y-SNPs sets were investigated by 2304 unrelated males from 24 populations. Furthermore, network and time estimation analyses were performed to evaluate the new strategy's capability to reconstruct phylogenetic tree and reliability to infer the time to the most recent common ancestor (TMRCA). The nine SM Y-STRs assay even had a higher resolution and a comparable capacity of revealing population genetic differentiation compared to 172 Y-SNPs system. This new strategy could optimize the phylogenetic tree generated by commonly used Y-STR panels and obtain a quite consistent time estimations with the published dating.

Potential of large "first generation" human-to-human transmission of 2019-nCoV.

To investigate the genetic diversity, time origin, and evolutionary history of the 2019-nCoV outbreak in China and Thailand, a total of 12 genome sequences of the virus with known sampling date (24 December 2019 and 13 January 2020) and geographic location (primarily Wuhan city, Hubei Province, China, but also Bangkok, Thailand) were analyzed. Phylogenetic and likelihood-mapping analyses of these genome sequences were performed. On the basis of our results, the star-like signal and topology of 2019-nCoV may be indicative of potentially large "first generation" human-to-human virus transmission. We estimated that 2019-nCoV likely originated in Wuhan on 9 November 2019 (95% credible interval: 25 September 2019 and 19 December 2019), and that Wuhan is the major hub for the spread of the 2019-nCoV outbreak in China and elsewhere. Our results could be useful for designing effective prevention strategies for 2019-nCoV in China and beyond.

Evolutionary history, potential intermediate animal host, and cross-species analyses of SARS-CoV-2.

To investigate the evolutionary history of the recent outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in China, a total of 70 genomes of virus strains from China and elsewhere with sampling dates between 24 December 2019 and 3 February 2020 were analyzed. To explore the potential intermediate animal host of the SARS-CoV-2 virus, we reanalyzed virome data sets from pangolins and representative SARS-related coronaviruses isolates from bats, with particular attention paid to the spike glycoprotein gene. We performed phylogenetic, split network, transmission network, likelihood-mapping, and comparative analyses of the genomes. Based on Bayesian time-scaled phylogenetic analysis using the tip-dating method, we estimated the time to the most recent common ancestor and evolutionary rate of SARS-CoV-2, which ranged from 22 to 24 November 2019 and 1.19 to 1.31 × 10-3 substitutions per site per year, respectively. Our results also revealed that the BetaCoV/bat/Yunnan/RaTG13/2013 virus was more similar to the SARS-CoV-2 virus than the coronavirus obtained from the two pangolin samples (SRR10168377 and SRR10168378). We also identified a unique peptide (PRRA) insertion in the human SARS-CoV-2 virus, which may be involved in the proteolytic cleavage of the spike protein by cellular proteases, and thus could impact host range and transmissibility. Interestingly, the coronavirus carried by pangolins did not have the RRAR motif. Therefore, we concluded that the human SARS-CoV-2 virus, which is responsible for the recent outbreak of COVID-19, did not come directly from pangolins.

Transmission dynamics and evolutionary history of 2019-nCoV.

To investigate the time origin, genetic diversity, and transmission dynamics of the recent 2019-nCoV outbreak in China and beyond, a total of 32 genomes of virus strains sampled from China, Thailand, and the USA with sampling dates between 24 December 2019 and 23 January 2020 were analyzed. Phylogenetic, transmission network, and likelihood-mapping analyses of the genome sequences were performed. On the basis of the likelihood-mapping analysis, the increasing tree-like signals (from 0% to 8.2%, 18.2%, and 25.4%) over time may be indicative of increasing genetic diversity of 2019-nCoV in human hosts. We identified three phylogenetic clusters using the Bayesian inference framework and three transmission clusters using transmission network analysis, with only one cluster identified by both methods using the above genome sequences of 2019-nCoV strains. The estimated mean evolutionary rate for 2019-nCoV ranged from 1.7926 × 10-3 to 1.8266 × 10-3 substitutions per site per year. On the basis of our study, undertaking epidemiological investigations and genomic data surveillance could positively impact public health in terms of guiding prevention efforts to reduce 2019-nCOV transmission in real-time.

Increased predominance of HIV-1 CRF01_AE and its recombinants in the Philippines.

The growth rate of new HIV infections in the Philippines was the fastest of any countries in the Asia-Pacific region between 2010 and 2016. To date, HIV-1 subtyping results in the Philippines have been determined by characterizing only partial viral genome sequences. It is not known whether recombination occurs in the majority of unsequenced genome regions. Near-full-length genome (NFLG) sequences were obtained by amplifying two overlapping half genomes from plasma samples collected between 2015 and 2017 from 23 newly diagnosed infected individuals in the Philippines. Phylogenetic analysis showed that the newly characterized sequences were CRF01_AE (14), subtype B (3), CRF01/B recombinants (5) and a CRF01/CRF07/B recombinant (1). All 14 CRF01_AE formed a tight cluster, suggesting that they were derived from a single introduction. The time to the most recent common ancestor (tMRCA) for CRF01_AE in the Philippines was 1995 (1992-1998), about 10-15 years later than that of CRF01_AE in China and Thailand. All five CRF01/B recombinants showed distinct recombination patterns, suggesting ongoing recombination between the two predominant circulating viruses. The identification of partial CRF07_BC sequences in one CRF01/CRF07/B recombinant, not reported previously in the Philippines, indicated that CRF07_BC may have been recently introduced into that country from China, where CRF07_BC is prevalent. Our results show that the major epidemic strains may have shifted to an increased predominance of CRF01_AE and its recombinants, and that other genotypes such as CRF07_BC may have been introduced into the Philippines.

Genetic characterization of bovine coronavirus in Vietnam.

A maximum clade credibility tree constructed using the full-length spike (S) and hemagglutinin-esterase genes revealed that Vietnamese Bovine coronavirus (BCoV) strains belong to a single cluster (C1); therefore, they might share a common origin with Cuban and Chinese BCoV strains. The omega values of cluster 1 (C1) and cluster 2 (C2) were 0.15734 and 0.11613, respectively, and naive empirical bayes analysis identified two amino acid positions (179 and 501) in the S protein in C1 and three amino acid positions (113, 501, and 525) in that of C2 that underwent positive selection (p > 99%). The evolutionary rate of C1 was estimated to be 7.6206 × 10-4 substitutions/site/year, and the most recent common ancestor (tMRCA) of Vietnamese BCoVs was estimated to date back to 1962 (95% HPD 1950-1973). The effective population sizes of C1 and C2 underwent a rapid reduction after 2000 and 2004, respectively.

Is Hepatitis Delta infections important in Brazil?

BACKGROUND: The Hepatitis Delta Virus (HDV) can increase the incidence of fulminant hepatitis. For this infection occurs, the host must also be infected with Hepatitis B Virus. Previous studies demonstrated the endemicity and near exclusivity of this infection in the Amazon region, and as a consequence of the difficulty in accessing this area we used dried blood spots (DBS) in sample collection. The aims of this study were to investigate the presence of recombination, to analyze the epidemiology, ancestry and evolutionary pressures on HDV in Brazil. METHODS: Blood samples from 50 individuals were collected using dried-blood spots (DBS 903, Whatman), and sent via regular mail to Retrovirology Laboratory from Federal University of São Paulo, where the samples were processed. In the analysis the following software were used: PhyML, RDP, BEAST, jModelTest and CODEML. RESULTS: Our results confirm the prevalence of HDV-3 in the Amazon region of Brazil, with the absence of inter-genotypic recombination. It was identified a positive selection in probable epitopes of HDV on B lymphocytes that might indicate that the virus is changing to escape the humoral response of the host. The analysis of the time of the most common ancestor demonstrated the exponential growth of this virus in late 1970s that lasted until 1995, after which it remained constant. It was also observed a probable founder effect in two cities, which demonstrate the need to focus on prevention methods against HBV/HDV infection. CONCLUSION: We confirmed the prevalence of HDV-3 in the Amazon region of Brazil, without inter-genotypic recombination. The analysis of the time of the most common ancestor showed that this infection remain constant in the studied area. Taking into account the probable founder effect established in the cities of Rio Branco and Porto Velho, a focus on preventive methods is recommended against these infections.

Population Growth Combined with Wide Offspring Distributions can Increase Fixation Rate and Reduce Genetic Diversity.

Genetic diversity grows with the population size in most neutral evolution models. Empirical evidence of large populations with limited diversity has been proposed to be due either to genetic bottlenecks or to selection. An alternative explanation is that the limited diversity is a result of rare reproduction events. Indeed, recent estimates of the offspring number distribution highlight the role of large reproduction events. We here show that in a large class of models containing such rare events, genetic diversity decreases as the population size increases, in neutral evolution models. For many realistic offspring number distributions, the contribution of rare events to the dynamics grows with the population size. In the context of genetic diversity, these rare events induce a decrease in the time to the most recent common ancestor and in the genetic heterogeneity as the population grows. This phenomenon may explain the observed rapid fixation of genes in large populations, in the absence of observable selection or bottlenecks.

Evolution of the viral hemorrhagic septicemia virus: divergence, selection and origin.

Viral hemorrhagic septicemia virus (VHSV) is an economically significant rhabdovirus that affects an increasing number of freshwater and marine fish species. Extensive studies have been conducted on the molecular epizootiology, genetic diversity, and phylogeny of VHSV. However, there are discrepancies between the reported estimates of the nucleotide substitution rate for the G gene and the divergence times for the genotypes. Herein, Bayesian coalescent analyses were conducted to the time-stamped entire coding sequences of the six VHSV genes. Rate estimates based on the G gene indicated that the marine genotypes/subtypes might not all evolve slower than their major European freshwater counterpart. Age calculations on the six genes revealed that the first bifurcation event of the analyzed isolates might have taken place within the last 300 years, which was much younger than previously thought. Selection analyses suggested that two codons of the G gene might be positively selected. Surveys of codon usage bias showed that the P, M and NV genes exhibited genotype-specific variations. Furthermore, we proposed that VHSV originated from the Pacific Northwest of North America.

Dating ancient splits in phylogenetic trees, with application to the human-Neanderthal split.

BACKGROUND: We tackle the problem of estimating species TMRCAs (Time to Most Recent Common Ancestor), given a genome sequence from each species and a large known phylogenetic tree with a known structure (typically from one of the species). The number of transitions at each site from the first sequence to the other is assumed to be Poisson distributed, and only the parity of the number of transitions is observed. The detailed phylogenetic tree contains information about the transition rates in each site. We use this formulation to develop and analyze multiple estimators of the species' TMRCA. To test our methods, we use mtDNA substitution statistics from the well-established Phylotree as a baseline for data simulation such that the substitution rate per site mimics the real-world observed rates. RESULTS: We evaluate our methods using simulated data and compare them to the Bayesian optimizing software BEAST2, showing that our proposed estimators are accurate for a wide range of TMRCAs and significantly outperform BEAST2. We then apply the proposed estimators on Neanderthal, Denisovan, and Chimpanzee mtDNA genomes to better estimate their TMRCA with modern humans and find that their TMRCA is substantially later, compared to values cited recently in the literature. CONCLUSIONS: Our methods utilize the transition statistics from the entire known human mtDNA phylogenetic tree (Phylotree), eliminating the requirement to reconstruct a tree encompassing the specific sequences of interest. Moreover, they demonstrate notable improvement in both running speed and accuracy compared to BEAST2, particularly for earlier TMRCAs like the human-Chimpanzee split. Our results date the human - Neanderthal TMRCA to be [Formula: see text] years ago, considerably later than values cited in other recent studies.

Paleopteran molecular clock: Time drift and recent acceleration.

Applying BEAST v1.10.4, we constructed a Bayesian Inference tree comprising 322 taxa, primarily representing Paleoptera (Odonata and Ephemeroptera; Pterygota), Zygentoma and Archaeognatha (Apterygota; paraphyly), and Neoptera (Plecoptera; Pterygota), based on a 2685 bp sequence dataset. Our analyses revealed that robust dating required the incorporation of both Quaternary and pre-Quaternary dates. To achieve this, our dating incorporated a 1.55 Ma (Quaternary) geological event (the formation of the Ryukyu Islands) and a set of chronologically well-founded fossil dates, spanning from up to 400 Ma (Devonian) for the stem Archaeognatha, 320 Ma (Carboniferous) for the crown of Paleoptera, 300 Ma (Carboniferous) for the crown Ephemeroptera, and 280 Ma (Permian) for the crown Odonata, down to 1.76 Ma (Quaternary) for Calopteryx japonica, encompassing a total of 22 calibration points (events: 6, fossils: 16; Quaternary: 7, pre-Quaternary: 15). The resulting dated tree aligns with previous research, albeit with some dates being overestimated. This overestimation was mainly due to the lack of Quaternary calibration and the exclusive dependence on pre-Quaternary calibration, though the application of maximum age constraints also played a role. Our minimum age dating demonstrates that the molecular clock did not uniformly progress, rendering rate dating an inapplicable approach. We observed that the base substitution rate is time-dependent, with an exponential increase evident from around 20 Ma (Miocene) to the present time, exceeding an order of magnitude. The extensive radiation and speciation of Insecta and Paleoptera, potentially resulting from the severe climatic changes associated with the Quaternary, including the commencement of glacial and interglacial cycles, may have significantly contributed to this increase in base substitution rates. Additionally, we identified a potential peak in base substitution rates during the Carboniferous period, around 320 million years ago, possibly corresponding to the Late Paleozoic Ice Age.

Bayesian Molecular Dating Analyses Combined with Mutational Profiling Suggest an Independent Origin and Evolution of SARS-CoV-2 Omicron BA.1 and BA.2 Sub-Lineages.

The ongoing evolution of severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) has resulted in the recent emergence of a highly divergent variant of concern (VOC) defined as Omicron or B.1.1.529. This VOC is of particular concern because it has the potential to evade most therapeutic antibodies and has undergone a sustained genetic evolution, resulting in the emergence of five distinct sub-lineages. However, the evolutionary dynamics of the initially identified Omicron BA.1 and BA.2 sub-lineages remain poorly understood. Herein, we combined Bayesian phylogenetic analysis, mutational profiling, and selection pressure analysis to track the virus's genetic changes that drive the early evolutionary dynamics of the Omicron. Based on the Omicron dataset chosen for the improved temporal signals and sampled globally between November 2021 and January 2022, the most recent common ancestor (tMRCA) and substitution rates for BA.1 were estimated to be that of 18 September 2021 (95% highest posterior density (HPD), 4 August-22 October 2021) and 1.435 × 10-3 (95% HPD = 1.021 × 10-3 - 1.869 × 10-3) substitution/site/year, respectively, whereas 3 November 2021 (95% highest posterior density (HPD) 26 September-28 November 2021) and 1.074 × 10-3 (95% HPD = 6.444 × 10-4 - 1.586 × 10-3) substitution/site/year were estimated for the BA.2 sub-lineage. The findings of this study suggest that the Omicron BA.1 and BA.2 sub-lineages originated independently and evolved over time. Furthermore, we identified multiple sites in the spike protein undergoing continued diversifying selection that may alter the neutralization profile of BA.1. This study sheds light on the ongoing global genomic surveillance and Bayesian molecular dating analyses to better understand the evolutionary dynamics of the virus and, as a result, mitigate the impact of emerging variants on public health.

SARS-CoV-2 Lineage A.27: New Data from African Countries and Dynamics in the Context of the COVID-19 Pandemic.

SARS-CoV-2 is constantly evolving with lineages emerging and others eclipsing. Some lineages have an important epidemiological impact and are known as variants of interest (VOIs), variants under monitoring (VUMs) or variants of concern (VOCs). Lineage A.27 was first defined as a VUM since it holds mutations of concern. Here, we report additional lineage A.27 data and sequences from five African countries and describe the molecular characteristics, and the genetic history of this lineage worldwide. Based on the new sequences investigated, the most recent ancestor (tMRCA) of lineage A.27 was estimated to be from April 2020 from Niger. It then spread to Europe and other parts of the world with a peak observed between February and April 2021. The detection rate of A.27 then decreased with only a few cases reported during summer 2021. The phylogenetic analysis revealed many sub-lineages. Among them, one was defined by the substitution Q677H in the spike (S) gene, one was defined by the substitution D358N in the nucleoprotein (N) gene and one was defined by the substitution A2143V in the ORF1b gene. This work highlights the importance of molecular characterization and the timely submission of sequences to correctly describe the circulation of particular strains in order to be proactive in monitoring the pandemic.