Amino Acid Substitution Rates Research Articles

QMix: An Efficient Program to Automatically Estimate Multi-Matrix Mixture Models for Amino Acid Substitution Process.

The single-matrix amino acid (AA) substitution models are widely used in phylogenetic analyses; however, they are unable to properly model the heterogeneity of AA substitution rates among sites. The multi-matrix mixture models can handle the site rate heterogeneity and outperform the single-matrix models. Estimating multi-matrix mixture models is a complex process and no computer program is available for this task. In this study, we implemented a computer program of the so-called QMix based on the algorithm of LG4X and LG4M with several enhancements to automatically estimate multi-matrix mixture models from large datasets. QMix employs QMaker algorithm instead of XRATE algorithm to accurately and rapidly estimate the parameters of models. It is able to estimate mixture models with different number of matrices and supports multi-threading computing to efficiently estimate models from thousands of genes. We re-estimate mixture models LG4X and LG4M from 1471 HSSP alignments. The re-estimated models (HP4X and HP4M) are slightly better than LG4X and LG4M in building maximum likelihood trees from HSSP and TreeBASE datasets. QMix program required about 10 hours on a computer with 18 cores to estimate a mixture model with four matrices from 200 HSSP alignments. It is easy to use and freely available for researchers.

Evolutionary analysis of seasonal influenza A viruses in Pakistan 2020-2023.

Influenza A viruses cause global health concerns due to their high amino acid substitution rates. They are linked to yearly seasonal epidemics and occasional pandemics. This study focused on sequencing influenza A virus strains in Pakistan. We analyzed the genetic characteristics of influenza A(H1N1)pdm09 and A(H3N2) viruses circulating in Pakistan from January 2020 to January 2023. Whole genome sequences from influenza A (n = 126) virus isolates were amplified and sequenced by the Oxford Nanopore (MinION) platform. The HA genes of influenza A(H1N1)pdm09 underwent amino acid substitutions at positions K54Q, A186T, Q189E, E224A, R259K, and K308R in sequenced samples. The HA genes of influenza A(H3N2) had amino acid substitutions at G53D, E83K, D104G, I140M, S205F, A212T, and K276R in the sequenced samples. Furthermore, the HA gene sequences of influenza A(H1N1)pdm09 in this study belonged to subclade 6B.1A.5a.2a. Similarly, the HA gene sequences of influenza A(H3N2) were classified under six subclades (3C.3a.1 and 3C.2a1b.2a [2, 2a.1, 2b, 2c, and 2a.3b]). Notably, amino acid substitutions in other gene segments of influenza A(H1N1)pdm09 and A(H3N2) were also found. These findings indicate influenza A(H1N1)pdm09 and A(H3N2) viruses co-circulated during the 2020-2023 influenza season in Pakistan. Continued surveillance is crucial for real-time monitoring of possible high-virulence variation and their relevance to existing vaccine strains.

Mitogenome-based phylogeny of the gastropod order Neomphalida points to multiple habitat shifts and a Pacific origin

Neomphalida is an order of gastropod molluscs with highly diverse morphology and a global distribution across various chemosynthesis-based ecosystems from organic falls to hot vents. The phylogenetic relationships of taxa within this order remain contentious, due to the rarity of material leading to a low taxonomic coverage and few genetic markers used. Neomphalida includes three families—Melanodrymiidae, Neomphalidae, and Peltospiridae—and molecular sequences are especially lacking in Melanodrymiidae. Here, we assembled a total of 11 mitogenomes covering these three families and 14 genus-level groups to reconstruct the most complete phylogeny of Neomphalida to date. Our current result recovered the monophyly of three families with maximum support and a likely interfamilial relationship of (Melanodrymiidae + Neomphalidae) + Peltospiridae. These indicate the possibility of habitat shifting from non-chemosynthetic deep sea to hot vent and then to sunken wood, accompanied by elevated mitogenome rearrangements and amino acid substitution rates in Melanodrymiidae. By mapping species distribution on the phylogeny, our findings suggest a Pacific origin of Neomphalida and multiple historical dispersal events of Peltospiridae to the Indian Ocean and at least once to the Atlantic.

An atlas of continuous adaptive evolution in endemic human viruses

Through antigenic evolution, viruses such as seasonal influenza evade recognition by neutralizing antibodies. This means that a person with antibodies well tuned to an initial infection will not be protected against the same virus years later and that vaccine-mediated protection will decay. To expand our understanding of which endemic human viruses evolve in this fashion, we assess adaptive evolution across the genome of 28 endemic viruses spanning a wide range of viral families and transmission modes. Surface proteins consistently show the highest rates of adaptation, and ten viruses in this panel are estimated to undergo antigenic evolution to selectively fix mutations that enable the escape of prior immunity. Thus, antibody evasion is not an uncommon evolutionary strategy among human viruses, and monitoring this evolution will inform future vaccine efforts. Additionally, by comparing overall amino acid substitution rates, we show that SARS-CoV-2 is accumulating protein-coding changes at substantially faster rates than endemic viruses.

Protein Modelling Highlighted Key Catalytic Sites Involved in Position-Specific Glycosylation of Isoflavonoids.

Uridine diphosphate glycosyltransferases (UGTs) are known for promiscuity towards sugar acceptors, a valuable characteristic for host plants but not desirable for heterologous biosynthesis. UGTs characterized for the O-glycosylation of isoflavonoids have shown a variable efficiency, substrate preference, and OH site specificity. Thus, 22 UGTs with reported isoflavonoid O-glycosylation activity were analyzed and ranked for OH site specificity and catalysis efficiency. Multiple-sequence alignment (MSA) showed a 33.2% pairwise identity and 4.5% identical sites among selected UGTs. MSA and phylogenetic analysis highlighted a comparatively higher amino acid substitution rate in the N-terminal domain that likely led to a higher specificity for isoflavonoids. Based on the docking score, OH site specificity, and physical and chemical features of active sites, selected UGTs were divided into three groups. A significantly high pairwise identity (67.4%) and identical sites (31.7%) were seen for group 1 UGTs. The structural and chemical composition of active sites highlighted key amino acids that likely define substrate preference, OH site specificity, and glycosylation efficiency towards selected (iso)flavonoids. In conclusion, physical and chemical parameters of active sites likely control the position-specific glycosylation of isoflavonoids. The present study will help the heterologous biosynthesis of glycosylated isoflavonoids and protein engineering efforts to improve the substrate and site specificity of UGTs.

Evolutionary Study of DGAT1 Gene Indicates the Presence of Positive Selection in Mammals

This study was focused on the evolving role of DGAT1 using phylogenetic analysis. DGAT1 was retrieved from available databases of 37 mammalian genomes to detect signatures of positive selection in individual codons. Evolutionary analyses likelihood (REL), fast unconstrained Bayesian approximation (FUBAR), and Mixed-effects model evolution (MEME). REL detected 3 positive selection sites (at 109, 195, and 211 positions), and MEME detected 20 sites. Among them, 211 sites were confirmed as a strong positive selection of the DGAT1 gene. The codon model evolution (CME) model clustered amino acid substitution rates at three class levels (0.06, 0.20, and 0.40). Multiple hit tests showed that 20 of the DGAT1 triplet codes for amino acid formation were changed into new codes. In which, the GAC code for aspartate and GTC for valine had the maximum counts of 14 and 16. Strong evidence of positive selection of the DGAT1 gene in mammals was found.

Phylogenetic distribution, structural analysis and interaction of nucleotide excision repair proteins in cyanobacteria

Cyanobacteria are photosynthetic Gram-negative, oxygen evolving prokaryotes with cosmopolitan distribution. Ultraviolet radiation (UVR) and other abiotic stresses result in DNA lesions in cyanobacteria. Nucleotide excision repair (NER) pathway removes the DNA lesions produced by UVR to normal DNA sequence. In cyanobacteria, detailed knowledge about NER proteins is poorly studied. Therefore, we have studied the NER proteins in cyanobacteria. Analyses of 289 amino acids sequence from 77 cyanobacterial species have revealed the presence of a minimum of one copy of NER protein in their genome. Phylogenetic analysis of NER protein shows that UvrD has maximal rate of amino acid substitutions which resulted in increased branch length. The motif analysis shows that UvrABC proteins is more conserved than UvrD, Further, UvrA with UvrB protein interacts with each other and form stable complex which have DNA binding domain on the surface of the complex. UvrB also have DNA binding domain. Positive electrostatic potential was found in the DNA binding region, which is followed by negative and neutral electrostatic potential. Additionally, the surface accessibility values at the DNA strands of T5-T6 dimer binding site were maximal. Protein nucleotide interaction shows the strong binding of T5-T6 dimer with NER proteins of Synechocystis sp. PCC 6803. This process repairs the UV-induced DNA lesions in dark when photoreactivation is inactive. Regulation of NER proteins protect cyanobacterial genome and maintain the fitness of organism under different abiotic stresses.

The Structure of Evolutionary Model Space for Proteins across the Tree of Life

The factors that determine the relative rates of amino acid substitution during protein evolution are complex and known to vary among taxa. We estimated relative exchangeabilities for pairs of amino acids from clades spread across the tree of life and assessed the historical signal in the distances among these clade-specific models. We separately trained these models on collections of arbitrarily selected protein alignments and on ribosomal protein alignments. In both cases, we found a clear separation between the models trained using multiple sequence alignments from bacterial clades and the models trained on archaeal and eukaryotic data. We assessed the predictive power of our novel clade-specific models of sequence evolution by asking whether fit to the models could be used to identify the source of multiple sequence alignments. Model fit was generally able to correctly classify protein alignments at the level of domain (bacterial versus archaeal), but the accuracy of classification at finer scales was much lower. The only exceptions to this were the relatively high classification accuracy for two archaeal lineages: Halobacteriaceae and Thermoprotei. Genomic GC content had a modest impact on relative exchangeabilities despite having a large impact on amino acid frequencies. Relative exchangeabilities involving aromatic residues exhibited the largest differences among models. There were a small number of exchangeabilities that exhibited large differences in comparisons among major clades and between generalized models and ribosomal protein models. Taken as a whole, these results reveal that a small number of relative exchangeabilities are responsible for much of the structure of the "model space" for protein sequence evolution. The clade-specific models we generated may be useful tools for protein phylogenetics, and the structure of evolutionary model space that they revealed has implications for phylogenomic inference across the tree of life.

Molecular Evolution of RNA-Dependent RNA Polymerase Region in Norovirus Genogroup I.

Norovirus is the leading viral agent of gastroenteritis in humans. RNA-dependent RNA polymerase (RdRp) is essential in the replication of norovirus RNA. Here, we present a comprehensive evolutionary analysis of the norovirus GI RdRp gene. Our results show that the norovirus GI RdRp gene can be divided into three groups, and that the most recent common ancestor was 1484. The overall evolutionary rate of GI RdRp is 1.821 × 10-3 substitutions/site/year. Most of the amino acids of the GI RdRp gene were under negative selection, and only a few positively selected sites were recognized. Amino acid substitutions in the GI RdRp gene accumulated slowly over time. GI.P1, GI.P3 and GI.P6 owned the higher evolutionary rates. GI.P11 and GI.P13 had the faster accumulation rate of amino acid substitutions. GI.P2, GI.P3, GI.P4, GI.P6 and GI.P13 presented a strong linear evolution. These results reveal that the norovirus GI RdRp gene evolves conservatively, and that the molecular evolutionary characteristics of each P-genotype are diverse. Sequencing in RdRp and VP1 of norovirus should be advocated in the surveillance system to explore the effect of RdRp on norovirus activity.

Proteome-Wide Structural Computations Provide Insights into Empirical Amino Acid Substitution Matrices

The relative contribution of mutation and selection to the amino acid substitution rates observed in empirical matrices is unclear. Herein, we present a neutral continuous fitness-stability model, inspired by the Arrhenius law (qij=aije-ΔΔGij). The model postulates that the rate of amino acid substitution (i→j) is determined by the product of a pre-exponential factor, which is influenced by the genetic code structure, and an exponential term reflecting the relative fitness of the amino acid substitutions. To assess the validity of our model, we computed changes in stability of 14,094 proteins, for which 137,073,638 in silico mutants were analyzed. These site-specific data were summarized into a 20 square matrix, whose entries, ΔΔGij, were obtained after averaging through all the sites in all the proteins. We found a significant positive correlation between these energy values and the disease-causing potential of each substitution, suggesting that the exponential term accurately summarizes the fitness effect. A remarkable observation was that amino acids that were highly destabilizing when acting as the source, tended to have little effect when acting as the destination, and vice versa (source → destination). The Arrhenius model accurately reproduced the pattern of substitution rates collected in the empirical matrices, suggesting a relevant role for the genetic code structure and a tuning role for purifying selection exerted via protein stability.

Parallel evolution of opsin visual pigments in hawkmoths by tuning of spectral sensitivities during transition from a nocturnal to a diurnal ecology.

Light environments differ dramatically between day and night. The transition between diurnal and nocturnal visual ecology has happened repeatedly throughout evolution in many species. However, the molecular mechanism underlying the evolution of vision in recent diurnal-nocturnal transition is poorly understood. Here, we focus on hawkmoths (Lepidoptera: Sphingidae) to address this question by investigating five nocturnal and five diurnal species. We performed RNA-sequencing analysis and identified opsin genes corresponding to the ultraviolet (UV), short-wavelength (SW) and long-wavelength (LW)-absorbing visual pigments. We found no significant differences in the expression patterns of opsin genes between the nocturnal and diurnal species. We then constructed the phylogenetic trees of hawkmoth species and opsins. The diurnal lineages had emerged at least three times from the nocturnal ancestors. The evolutionary rates of amino acid substitutions in the three opsins differed between the nocturnal and diurnal species. We found an excess number of parallel amino acid substitutions in the opsins in three independent diurnal lineages. The numbers were significantly more than those inferred from neutral evolution, suggesting that positive selection acted on these parallel substitutions. Moreover, we predicted the visual pigment absorption spectra based on electrophysiologically determined spectral sensitivity in two nocturnal and two diurnal species belonging to different clades. In the diurnal species, the LW pigments shift 10 nm towards shorter wavelengths, and the SW pigments shift 10 nm in the opposite direction. Taken together, our results suggest that parallel evolution of opsins may have enhanced the colour discrimination properties of diurnal hawkmoths in ambient light.

Molecular characterisation of Lymphocystis Disease Virus detected from Sea Bream in Turkey

Lymphocystis disease is a common, chronic and benign infection of fish. The disease has a worldwide distribution and has been reported from more than 125 marine and freshwater sites. In this study, molecular diagnosis and characterisation of Lymphocystis Disease Virus were performed from samples collected from sea bream fry showing lymphocystis disease symptoms in a farm located in the Aegean Sea. As a result of the phylogenetic analysis, the virus was classified as genotype VII. According to the similarity of nucleotide sequence, SA1.ETun.2011 virus reported from Tunisia was the closest related virus. The rate of amino acid substitution between genotypes was quite high, while amongst the genotype VII isolates substitution rates were low. When the similarity between the genotypes was investigated using the nucleotide sequences, the closest and most distant genotypes to the genotype VII were V and I, respectively.

Contrasting modes of mitochondrial genome evolution in sister taxa of wood-eating marine bivalves (Teredinidae and Xylophagaidae).

The bivalve families Teredinidae and Xylophagaidae include voracious consumers of wood in shallow-water and deep-water marine environments, respectively. The taxa are sister clades whose members consume wood as food with the aid of intracellular cellulolytic endosymbionts housed in their gills. This combination of adaptations is found in no other group of animals and was likely present in the common ancestor of both families. Despite these commonalities, the two families have followed dramatically different evolutionary paths with respect to anatomy, life history, and distribution. Here, we present 42 new mitochondrial genome sequences from Teredinidae and Xylophagaidae and show that distinct trajectories have also occurred in the evolution and organization of their mitochondrial genomes. Teredinidae display significantly greater rates of amino acid substitution but absolute conservation of protein-coding gene order, whereas Xylophagaidae display significantly less amino acid change but have undergone numerous and diverse changes in genome organization since their divergence from a common ancestor. As with many bivalves, these mitochondrial genomes encode 2 ribosomal RNAs, 12 protein-coding genes, and 22 tRNAs; atp8 was not detected. We further show that their phylogeny, as inferred from amino acid sequences of 12 concatenated mitochondrial protein-coding genes, is largely congruent with those inferred from their nuclear genomes based on 18S and 28S ribosomal RNA sequences. Our results provide a robust phylogenetic framework to explore the tempo and mode of mitochondrial genome evolution and offer directions for future phylogenetic and taxonomic studies of wood-boring bivalves.

Parthenogenesis doubles the rate of amino acid substitution in whiptail mitochondria.

Sexual reproduction is ubiquitous in the natural world, suggesting that sex must have extensive benefits to overcome the cost of males compared to asexual reproduction. One hypothesized advantage of sex with strong theoretical support is that sex plays a role in removing deleterious mutations from the genome. Theory predicts that transitions to asexuality should lead to the suppression of recombination and segregation and, in turn, weakened natural selection, allowing for the accumulation of slightly deleterious mutations. We tested this prediction by estimating the dN/dS ratios in asexual vertebrate lineages in the genus Aspidoscelis using whole mitochondrial genomes from seven asexual and five sexual species. We found higher dN/dS ratios in asexual Aspidoscelis species, indicating that asexual whiptails accumulate nonsynonymous substitutions due to weaker purifying selection. Additionally, we estimated nucleotide diversity and found that asexuals harbor significantly less diversity. Thus, despite their recent origins, slightly deleterious mutations accumulated rapidly enough in asexual lineages to be detected. We provide empirical evidence to corroborate the connection between asexuality and increased amino acid substitutions in asexual vertebrate lineages.

Factors That Affect the Rates of Adaptive and Nonadaptive Evolution at the Gene Level in Humans and Chimpanzees.

The rate of amino acid substitution has been shown to be correlated to a number of factors including the rate of recombination, the age of the gene, the length of the protein, mean expression level, and gene function. However, the extent to which these correlations are due to adaptive and nonadaptive evolution has not been studied in detail, at least not in hominids. We find that the rate of adaptive evolution is significantly positively correlated to the rate of recombination, protein length and gene expression level, and negatively correlated to gene age. These correlations remain significant when each factor is controlled for in turn, except when controlling for expression in an analysis of protein length; and they also generally remain significant when biased gene conversion is taken into account. However, the positive correlations could be an artifact of population size contraction. We also find that the rate of nonadaptive evolution is negatively correlated to each factor, and all these correlations survive controlling for each other and biased gene conversion. Finally, we examine the effect of gene function on rates of adaptive and nonadaptive evolution; we confirm that virus-interacting proteins (VIPs) have higher rates of adaptive and lower rates of nonadaptive evolution, but we also demonstrate that there is significant variation in the rate of adaptive and nonadaptive evolution between GO categories when removing VIPs. We estimate that the VIP/non-VIP axis explains about 5–8 fold more of the variance in evolutionary rate than GO categories.

Dynamic Patterns of Sex Chromosome Evolution in Neognath Birds: Many Independent Barriers to Recombination at the ATP5F1A Locus

Avian sex chromosomes evolved after the divergence of birds and crocodilians from their common ancestor, so they are younger than the better-studied chromosomes of mammals. It has long been recognized that there may have been several stages to the evolution of avian sex chromosomes. For example, the CHD1 undergoes recombination in paleognaths but not neognaths. Genome assemblies have suggested that there may be variation in the timing of barriers to recombination among Neognathae, but there remains little understanding of the extent of this variability. Here, we look at partial sequences of ATP5F1A, which is on the avian Z and W chromosomes. It is known that recombination of this gene has independently ceased in Galliformes, Anseriformes, and at least five neoavian orders, but whether there are other independent cessations of recombination among Neoaves is not understood. We analyzed a combination of data extracted from published chromosomal-level genomes with data collected using PCR and cloning to identify Z and W copies in 22 orders. Our results suggest that there may be at least 19 independent cessations of recombination within Neognathae, and 3 clades that may still be undergoing recombination (or have only recently ceased recombination). Analyses of ATP5F1A protein sequences revealed an increased amino acid substitution rate for W chromosome gametologs, suggesting relaxed purifying selection on the W chromosome. Supporting this hypothesis, we found that the increased substitution rate was particularly pronounced for buried residues, which are expected to be more strongly constrained by purifying selection. This highlights the dynamic nature of avian sex chromosomes, and that this level of variation among clades means they should be a good system to understand sex chromosome evolution.

Rapid sequence evolution is associated with genetic incompatibilities in the plastid Clp complex.

Replacing the native clpP1 gene in the Nicotiana plastid genome with homologs from different donor species showed that the extent of genetic incompatibilities depended on the rate of sequence evolution. The plastid caseinolytic protease (Clp) complex plays essential roles in maintaining protein homeostasis and comprises both plastid-encoded and nuclear-encoded subunits. Despite the Clp complex being retained across green plants with highly conserved protein sequences in most species, examples of extremely accelerated amino acid substitution rates have been identified in numerous angiosperms. The causes of these accelerations have been the subject of extensive speculation but still remain unclear. To distinguish among prevailing hypotheses and begin to understand the functional consequences of rapid sequence divergence in Clp subunits, we used plastome transformation to replace the native clpP1 gene in tobacco (Nicotiana tabacum) with counterparts from another angiosperm genus (Silene) that exhibits a wide range in rates of Clp protein sequence evolution. We found that antibiotic-mediated selection could drive a transgenic clpP1 replacement from a slowly evolving donor species (S. latifolia) to homoplasmy but that clpP1 copies from Silene species with accelerated evolutionary rates remained heteroplasmic, meaning that they could not functionally replace the essential tobacco clpP1 gene. These results suggest that observed cases of rapid Clp sequence evolution are a source of epistatic incompatibilities that must be ameliorated by coevolutionary responses between plastid and nuclear subunits.

Ancient rapid functional differentiation and fixation of the duplicated members in rice Dof genes after whole genome duplication.

Whole genome duplication (WGD) in plants is typically followed by genomic downsizing, where large portions of the new genome are lost. Whether this downsizing is accompanied by increased or decreased evolutionary rates of the remaining genes is poorly known, not least because homeolog pairings are often obscured by chromosomal rearrangement. Here, we use the newly published genome from a sedge, namely Kobresia littledalei, and CRISPR/Cas-9 editing to investigate how the Rho WGD event 70 million years ago (MYA) affected transcription factor evolutionary rates, fates, and function in rice (Oryza sativa) and sorghum (Sorghum bicolor). We focus on the 30-member DNA-binding with one zinc finger (Dof) transcription factor family in both crops due to their agronomic importance. Using the known speciation dates of rice from Kobresia (97 MYA) and sorghum (50 MYA), we find that rates of amino acid substitution in the critical Dof domain region were over twofold higher during the 20-million-year period following the WGD than before or afterward. Through comparison of synteny blocks, we report that at least 11% of Dof genes were purged from 70 to 50 MYA, while only 6% have been lost in the most recent 50-million-year interval. CRISPR/Cas9 editing revealed widespread fitness-related defects in flowering and lack of redundancy of paired members, as well as significant differences in expression between gene pairs. Together these findings demonstrate the strength of Dof genes as a model for deep evolutionary study and offer one of the most detailed portraits yet of the Rho WGD impact on a gene lineage.

Temporal spread and evolution of SARS-CoV-2 in the second pandemic wave in Brazil.

Severe acute respiratory syndrome coronavirus‐2 (SARS‐CoV‐2) pandemic spread rapidly and this scenario is concerning in South America, mainly in Brazil that presented more than 21 million coronavirus disease 2019 cases and 590 000 deaths. The recent emergence of novel lineages carrying several mutations in the spike protein has raised additional public health concerns worldwide. The present study describes the temporal spreading and evolution of SARS‐CoV2 in the beginning of the second pandemic wave in Brazil, highlighting the fast dissemination of the two major concerning variants (P.1 and P.2). A total of 2507 SARS‐CoV‐2 whole‐genome sequences (WGSs) with available information from the country (Brazil) and sampling date (July 2020–February 2021), were obtained and the frequencies of the lineages were evaluated in the period of the growing second pandemic wave. The results demonstrated the increasing prevalence of P.1 and P.2 lineages in the period evaluated. P.2 lineage was first detected in the middle of 2020, but a high increase occurred only in the last trimester of this same year and the spreading to all Brazilian regions. P.1 lineage emerged even later, first in the North region in December 2020 and really fast dissemination to all other Brazilian regions in January and February 2021. All SARS‐CoV‐2 WGSs of P.1 and P.2 were further separately evaluated with a Bayesian approach. The rates of nucleotide and amino acid substitutions were statistically higher in P.1 than P.2 (p < 0.01). The phylodynamic analysis demonstrated that P.2 gradually spread in all the country from September 2020 to January 2021, while P.1 disseminated even faster from December 2020 to February 2021. Skyline plots of both lineages demonstrated a slight rise in the spreading for P.2 and exponential growth for P.1. In conclusion, these data demonstrated that the P.1 (recently renamed as Gamma) and P.2 lineages have predominated in the second pandemic wave due to the very high spreading across all geographic regions in Brazil at the end of 2020 and beginning of 2021.

Genetic Characterizations and Molecular Evolution of the Measles Virus Genotype B3's Hemagglutinin (H) Gene in the Elimination Era.

Measles virus (MeV) genotype B3 is one globally significant circulating genotype. Here, we present a systematic description of long-term evolutionary characterizations of the MeV genotype B3’s hemagglutinin (H) gene in the elimination era. Our results show that the B3 H gene can be divided into two main sub-genotypes, and the highest intra-genotypic diversity was observed in 2004. MeV genotype B3’s H gene diverged in 1976; its overall nucleotide substitution rate is estimated to be 5.697 × 10−4 substitutions/site/year, and is slowing down. The amino acid substitution rate of genotype B3’s H gene is also decreasing, and the mean effective population size has been in a downward trend since 2000. Selection pressure analysis only recognized a few sites under positive selection, and the number of positive selection sites is getting smaller. All of these observations may reveal that genotype B3’s H gene is not under strong selection pressure, and is becoming increasingly conservative. MeV H-gene or whole-genome sequencing should be routine, so as to better elucidate the molecular epidemiology of MeV in the future.