FUBAR: a fast, unconstrained bayesian approximation for inferring selection

More than mcr: canonical plasmid- and transposon-encoded mobilized colistin resistance genes represent a subset of phosphoethanolamine transferases

Mobilized colistin resistance genes (mcr) may confer resistance to the last-resort antimicrobial colistin and can often be transmitted horizontally. mcr encode phosphoethanolamine transferases (PET), which are closely related to chromosomally encoded, intrinsic lipid modification PET (i-PET; e.g., EptA, EptB, CptA). To gain insight into the evolution of mcr within the context of i-PET, we identified 69,814 MCR-like proteins present across 256 bacterial genera (obtained by querying known MCR family representatives against the National Center for Biotechnology Information [NCBI] non-redundant protein database via protein BLAST). We subsequently identified 125 putative novel mcr-like genes, which were located on the same contig as (i) ≥1 plasmid replicon and (ii) ≥1 additional antimicrobial resistance gene (obtained by querying the PlasmidFinder database and NCBI's National Database of Antibiotic Resistant Organisms, respectively, via nucleotide BLAST). At 80% amino acid identity, these putative novel MCR-like proteins formed 13 clusters, five of which represented putative novel MCR families. Sequence similarity and a maximum likelihood phylogeny of mcr, putative novel mcr-like, and ipet genes indicated that sequence similarity was insufficient to discriminate mcr from ipet genes. A mixed-effect model of evolution (MEME) indicated that site- and branch-specific positive selection played a role in the evolution of alleles within the mcr-2 and mcr-9 families. MEME suggested that positive selection played a role in the diversification of several residues in structurally important regions, including (i) a bridging region that connects the membrane-bound and catalytic periplasmic domains, and (ii) a periplasmic loop juxtaposing the substrate entry tunnel. Moreover, eptA and mcr were localized within different genomic contexts. Canonical eptA genes were typically chromosomally encoded in an operon with a two-component regulatory system or adjacent to a TetR-type regulator. Conversely, mcr were represented by single-gene operons or adjacent to pap2 and dgkA, which encode a PAP2 family lipid A phosphatase and diacylglycerol kinase, respectively. Our data suggest that eptA can give rise to "colistin resistance genes" through various mechanisms, including mobilization, selection, and diversification of genomic context and regulatory pathways. These mechanisms likely altered gene expression levels and enzyme activity, allowing bona fide eptA to evolve to function in colistin resistance.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

Evaluating the role of selection in the evolution of mitochondrial genomes of aboriginal peoples of Siberia

Studies of the nature of mitochondrial DNA (mtDNA) variability in human populations have shown that protein-coding genes are under negative (purifying) selection, since their mutation spectra are characterized by a pronounced predominance of synonymous substitutions over non-synonymous ones (Ka/Ks < 1). Meanwhile, a number of studies have shown that the adaptation of populations to various environmental conditions may be accompanied by a relaxation of negative selection in some mtDNA genes. For example, it was previously found that in Arctic populations, negative selection is relaxed in the mitochondrial ATP6 gene, which encodes one of the subunits of ATP synthase. In this work, we performed a Ka/Ks analysis of mitochondrial genes in large samples of three regional population groups in Eurasia: Siberia (N = 803), Western Asia/Transcaucasia (N = 753), and Eastern Europe (N = 707). The main goal of this work is to search for traces of adaptive evolution in the mtDNA genes of aboriginal peoples of Siberia represented by populations of the north (Koryaks, Evens) and the south of Siberia and the adjacent territory of Northeast China (Buryats, Barghuts, Khamnigans). Using standard Ka/Ks analysis, it was found that all mtDNA genes in all studied regional population groups are subject to negative selection. The highest Ka/Ks values in different regional samples were found in almost the same set of genes encoding subunits of ATP synthase (ATP6, ATP8), NADH dehydrogenase complex (ND1, ND2, ND3), and cytochrome bc1 complex (CYB). The highest Ka/Ks value, indicating a relaxation of negative selection, was found in the ATP6 gene in the Siberian group. The results of the analysis performed using the FUBAR method (HyPhy software package) and aimed at searching for mtDNA codons under the influence of selection also showed the predominance of negative selection over positive selection in all population groups. In Siberian populations, nucleotide sites that are under positive selection and associated with mtDNA haplogroups were registered not in the north (which is expected under the assumption of adaptive evolution of mtDNA), but in the south of Siberia.

Infectivity of three Mayaro Virus geographic isolates in human cell lines

Mayaro virus (MAYV) is an emergent arthropod-borne virus that causes an acute febrile illness accompanied by arthralgia, similar to chikungunya virus. Increasing urbanization of MAYV outbreaks in the Americas has led to concerns for geographic expansion and spillover. Given the potential importance of this pathogen, we sought to fill critical gaps in knowledge regarding MAYV infectivity and geographic variation. This study describes the cytopathogenicity of MAYV in human dermal fibroblasts, human skeletal muscle satellite cells, human embryonic kidney cells (HEK), peripherally derived human macrophages, and Vero cells. We found that regional differences between these viruses do not affect replication kinetics, with high titers peaking at 37 h post infection. MAYV-U, did however, cause the most cytopathic effect in a time-dependent manner. Compared to the other two prototypic isolates, MAYV-U harbors unique mutations in the E2 protein, D60G and S205F, that are likely to interact with the host cell receptor and could affect infectivity. We further demonstrate that pre-treatment of cells with interferon-β inhibited viral replication in a dose-dependent manner. Together, these findings advance our understanding of MAYV infection of human target cells and provide initial data regarding variation according to geography.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Host-specific sensing of coronaviruses and picornaviruses by the CARD8 inflammasome

Hosts have evolved diverse strategies to respond to microbial infections, including the detection of pathogen-encoded proteases by inflammasome-forming sensors such as NLRP1 and CARD8. Here, we find that the 3CL protease (3CLpro) encoded by diverse coronaviruses, including Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), cleaves a rapidly evolving region of human CARD8 and activates a robust inflammasome response. CARD8 is required for cell death and the release of pro-inflammatory cytokines during SARS-CoV-2 infection. We further find that natural variation alters CARD8 sensing of 3CLpro, including 3CLpro-mediated antagonism rather than activation of megabat CARD8. Likewise, we find that a single nucleotide polymorphism (SNP) in humans reduces CARD8's ability to sense coronavirus 3CLpros and, instead, enables sensing of 3C proteases (3Cpro) from select picornaviruses. Our findings demonstrate that CARD8 is a broad sensor of viral protease activities and suggests that CARD8 diversity contributes to inter- and intraspecies variation in inflammasome-mediated viral sensing and immunopathology.

The emergence, impact, and evolution of human metapneumovirus variants from 2014 to 2021 in Spain

BACKGROUND

Human metapneumovirus (HMPV) is an important aetiologic agent of respiratory tract infection (RTI). This study aimed to describe the prevalence, genetic diversity, and evolutionary dynamics of HMPV.

METHODS

Laboratory-confirmed HMPV were characterised based on partial-coding G gene sequences with MEGA.v6.0. WGS was performed with Illumina, and evolutionary analyses with Datamonkey and Nextstrain.

RESULTS

HMPV prevalence was 2.5%, peaking in February-April and with an alternation in the predominance of HMPV-A and -B until the emergence of SARS-CoV-2, not circulating until summer and autumn-winter 2021, with a higher prevalence and with the almost only circulation of A2c_111dup. G and SH proteins were the most variable, and 70% of F protein was under negative selection. Mutation rate of HMPV genome was 6.95 ×10^-4 substitutions/site/year.

CONCLUSION

HMPV showed a significant morbidity until the emergence of SARS-CoV-2 pandemic in 2020, not circulating again until summer and autumn 2021, with a higher prevalence and with almost the only circulation of A2c_111dup, probably due to a more efficient immune evasion mechanism. The F protein showed a very conserved nature, supporting the need for steric shielding. The tMRCA showed a recent emergence of the A2c variants carrying duplications, supporting the importance of virological surveillance.

The plastome reveals new insights into the evolutionary and domestication history of peonies in East Asia

BACKGROUD

Paeonia holds considerable value in medicinal, ornamental horticultural, and edible oil industries, but the incomplete state of phylogenetic research in this genus poses a challenge to the effective conservation and development of wild germplasm, and also impedes the practical utilization of existing cultivars. Due to its uniparental inheritance and lack of recombination, the plastome (i.e., plastid genome), which is a valuable molecular marker for phylogenetic analyses, is characterized by an appropriate rate of nucleotide evolution.

METHODS

In this study, 10 newly assembled data and available reported data were combined to perform a comparative genomics and phylogenetics analysis of 63 plastomes of 16 Paeonia species, primarily from East Asia, which is the origin and diversity center of Paeonia.

RESULTS

Ranging between 152,153 and 154,405 bp, most plastomes displayed a conserved structure and relatively low nucleotide diversity, except for six plastomes, which showed obvious IR construction or expansion. A total of 111 genes were annotated in the Paeonia plastomes. Four genes (rpl22, rps3, rps19 and ycf1) showed different copy numbers among accessions while five genes (rpl36, petN, psbI, rpl33 and psbJ) showed strong codon usage biases (ENC < 35). Additional selection analysis revealed that no genes were under positive selection during the domestication of tree peony cultivars whereas four core photosynthesis-related genes (petA, psaA, psaB and rbcL) were under positive selection in herbaceous peony cultivars. This discovery might contribute to the wide adaption of these cultivars. Two types of molecular markers (SSR and SNP) were generated from the 63 plastomes. Even though SSR was more diverse than SNP, it had a weaker ability to delimit Paeonia species than SNP. The reconstruction of a phylogenetic backbone of Paeonia in East Asia revealed significant genetic divergence within the P. ostii groups. Evidence also indicated that the majority of P. suffruticosa cultivars had a maternal origin, from P. ostii. The results of this research also suggest that P. delavayi var. lutea, which likely resulted from hybridization with P. ludlowii, should be classified as a lineage within the broader P. delavayi group.

CONCLUSIONS

Overall, this study's research findings suggest that the Paeonia plastome is highly informative for phylogenetic and comparative genomic analyses, and could be useful in future research related to taxonomy, evolution, and domestication.

Adaptive Selection in the Evolution of Aquaglyceroporins in Mammals

Aquaporins (AQPs) are integral membrane proteins responsible for water transport across cellular membranes in both prokaryotes and eukaryotes. A subfamily of AQPs, known as aquaglyceroporins (AQGPs), facilitate the transport of small solutes such as glycerol, water, and other solutes across cellular membranes. These proteins are involved in a variety of physiological processes, such as organogenesis, wound healing, and hydration. Although AQPs have been studied extensively in different species, their conservation patterns, phylogenetic relationships, and evolution in mammals remain unexplored. In the present study, 119 AQGP coding sequences from 31 mammalian species were analysed to identify conserved residues, gene organisation, and most importantly, the nature of AQGP gene selection. Repertoire analysis revealed the absence of AQP7, 9, and 10 genes in certain species of Primates, Rodentia, and Diprotodontia, although not all three genes were absent in a single species. Two Asparagine-Proline-Alanine (NPA) motifs located at the N- and C-terminal ends, aspartic acid (D) residues, and the ar/R region were conserved in AQP3, 9, and 10. Six exons encoding the functional MIP domain of AQGP genes were found to be conserved across mammalian species. Evolutionary analysis indicated signatures of positive selection in AQP7, 9, and 10 amongst different mammalian lineages. Furthermore, substitutions of certain amino acids located close to critical residues may alter AQGP functionality, which is crucial for substrate selectivity, pore formation, and transport efficiency required for the maintenance of homeostasis in different mammalian species.

Evolution of insect innate immunity through domestication of bacterial toxins

Toxin cargo genes are often horizontally transferred by phages between bacterial species and are known to play an important role in the evolution of bacterial pathogenesis. Here, we show how these same genes have been horizontally transferred from phage or bacteria to animals and have resulted in novel adaptations. We discovered that two widespread bacterial genes encoding toxins of animal cells, cytolethal distending toxin subunit B (cdtB) and apoptosis-inducing protein of 56 kDa (aip56), were captured by insect genomes through horizontal gene transfer from bacteria or phages. To study the function of these genes in insects, we focused on Drosophila ananassae as a model. In the D. ananassae subgroup species, cdtB and aip56 are present as singular (cdtB) or fused copies (cdtB::aip56) on the second chromosome. We found that cdtB and aip56 genes and encoded proteins were expressed by immune cells, some proteins were localized to the wasp embryo's serosa, and their expression increased following parasitoid wasp infection. Species of the ananassae subgroup are highly resistant to parasitoid wasps, and we observed that D. ananassae lines carrying null mutations in cdtB and aip56 toxin genes were more susceptible to parasitoids than the wild type. We conclude that toxin cargo genes were captured by these insects millions of years ago and integrated as novel modules into their innate immune system. These modules now represent components of a heretofore undescribed defense response and are important for resistance to parasitoid wasps. Phage or bacterially derived eukaryotic toxin genes serve as macromutations that can spur the instantaneous evolution of novelty in animals.

Genomic epidemiology of the rotavirus G2P[4] strains in coastal Kenya pre- and post-rotavirus vaccine introduction, 2012-8

The introduction of rotavirus vaccines into the national immunization programme in many countries has led to a decline in childhood diarrhoea disease burden. Coincidentally, the incidence of some rotavirus group A (RVA) genotypes has increased, which may result from non-vaccine-type replacement. Here, we investigate the evolutionary genomics of rotavirus G2P[4] which has shown an increase in countries that introduced the monovalent Rotarix® vaccine. We examined sixty-three RVA G2P[4] strains sampled from children (aged below 13 years) admitted to Kilifi County Hospital, coastal Kenya, pre- (2012 to June 2014) and post-(July 2014 to 2018) rotavirus vaccine introduction. All the sixty-three genome sequences showed a typical DS-1-like genome constellation (G2-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2). Pre-vaccine G2 sequences predominantly classified as sub-lineage IVa-3 and co-circulated with low numbers of sub-lineage IVa-1 strains, whereas post-vaccine G2 sequences mainly classified into sub-lineage IVa-3. In addition, in the pre-vaccine period, P[4] sub-lineage IVa strains co-circulated with low numbers of P[4] lineage II strains, but P[4] sub-lineage IVa strains predominated in the post-vaccine period. On the global phylogeny, the Kenyan pre- and post-vaccine G2P[4] strains clustered separately, suggesting that different virus populations circulated in the two periods. However, the strains from both periods exhibited conserved amino acid changes in the known antigenic epitopes, suggesting that replacement of the predominant G2P[4] cluster was unlikely a result of immune escape. Our findings demonstrate that the pre- and post-vaccine G2P[4] strains circulating in Kilifi, coastal Kenya, differed genetically but likely were antigenically similar. This information informs the discussion on the consequences of rotavirus vaccination on rotavirus diversity.

Phylodynamics and Molecular Mutations of the Hemagglutinin Affecting Global Transmission and Host Adaptation of H5Nx Viruses

Highly pathogenic avian influenza (HPAI) H5 viruses have circulated globally causing incidental human infection with a substantial pandemic threat. The present study investigated the molecular evolution and phylodynamics of hemagglutinin (HA) in avian and human-isolated H5Nx viruses globally circulating since 2000. We investigated the dynamics of amino acid substitution in the HA sequences of avian and human H5Nx viruses and performed a phylogenetic analysis. Our study found that the H5Nx lineages dominantly expanded since 2000 and diverged into multiple sublineages with unique genetic mutations. P185S mutation in HA became a molecular characteristic of dominant H5Nx viruses throughout clades 2.3.4.1 to 2.3.4.4 (2.3.4.1-4). The key mutations, ΔE130 and I155T, and potential N-linked glycosylation at residues 128, 144, and 159 in the HA gene of human-isolated viruses possibly contributed to both the individual and population levels of the H5 evolution and the host adaptation. Our analysis detected heterogeneity in amino acid sites under positive selection in the HA gene of clades 2.3.4.1-4. Accumulated mutations in the HA protein may potentially affect not only the genetic and antigenic diversity of HPAI H5Nx viruses but also increase the functional compatibility with NA subtypes. Given the global spread and incessantly occurring HA mutations of H5Nx viruses, our results emphasize the importance of early identification of HA mutations as well as the need for a comprehensive assessment of H5Nx variants in terms of pandemic preparedness.

Molecular characterization of human respiratory syncytial virus in Seoul, South Korea, during 10 consecutive years, 2010-2019

Respiratory syncytial virus (RSV) is the leading cause of lower respiratory tract infections and hospitalization in infants and young children. Here, we analyzed the genetic diversity of RSV using partial G gene sequences in 84 RSV-A and 78 RSV- B positive samples collected in Seoul, South Korea, for 10 consecutive years, from 2010 to 2019. Our phylogenetic analysis revealed that RSV-A strains were classified into either the ON1 (80.9%) or NA1 (19.0%) genotypes. On the other hand, RSV-B strains demonstrated diversified clusters within the BA genotype. Notably, some sequences designated as BA-SE, BA-SE1, and BA-DIS did not cluster with previously identified BA genotypes in the phylogenetic trees. Despite this, they did not meet the criteria for the assignment of a new genotype based on recent classification methods. Selection pressure analysis identified three positive selection sites (amino acid positions 273, 274, and 298) in RSV-A, and one possible positive selection site (amino acid position 296) in RSV-B, respectively. The mean evolutionary rates of Korean RSV-A from 1999 to 2019 and RSV-B strains from 1991 and 2019 were estimated at 3.51 × 10-3 nucleotides (nt) substitutions/site/year and 3.32 × 10-3 nt substitutions/site/year, respectively. The population dynamics in the Bayesian skyline plot revealed fluctuations corresponding to the emergence of dominant strains, including a switch of the dominant genotype from NA1 to ON1. Our study on time-scaled cumulative evolutionary analysis contributes to a better understanding of RSV epidemiology at the local level in South Korea.

Evaluation of codon usage patterns and molecular evolution dynamics in Japanese encephalitis virus: An integrated bioinformatics approach

In the recent survey, Japanese encephalitis (JE) is one of the most common mosquito-borne diseases, accounting for ∼30% of fatalities. The outbreaks of the JE virus (JEV) suggests that exhaustive study is essential for the prevention and management of the disease. The disease mainly spreads in humans and pigs by the vector: mosquito; as this is a major concern, this study had employed various bioinformatics tools to investigate the codon usage bias, evolutionary inference and selection pressure analysis of the Japanese encephalitis virus disease. The results indicated that the JE virus was biased and natural selection was the main factor shaping the codon usage that was determined and confirmed with the Nc, neutrality, PR2 plots and correlation analysis. The evolutionary analysis revealed that the virus had a substitution rate of 1.54 × 10-4 substitution/site/year and the tMRCA was found to be in 1723. The transmission of the virus in the map found transmissions mostly from China and transmitted across Asia and Africa. The selection pressure analysis employed three methods which had 969th codon site as diversifying site and had many purifying sites that shows the virus had evolved rapidly. The inferences from this study would aid people to employ this methodology on various diseases and also perform insilico studies in the field of vaccinology and immunoinformatics.

Molecular Characteristics and Pathogenicity of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) 1 in Taiwan during 2019-2020

Two variants of porcine reproductive and respiratory syndrome virus (PRRSV), PRRSV 1 and PRRSV 2, have caused abortion in pregnant sows and respiratory distress in nursery pigs worldwide. PRRSV 2 has been thoroughly researched in Taiwan since 1993; however, the first case of PRRSV 1 was not reported until late 2018. To decipher the genetic characteristics of PRRSV 1 in Taiwan, open reading frame 5 (ORF5) genes of PRRSV 1 strains collected from 11 individual pig farms in 2019-2020 were successfully sequenced. All Taiwanese ORF5 sequences were closely related to Spanish-like PRRSV strains, which are considered to share a common evolutionary origin with the strain used for the PRRSV 1 vaccine. Analyses of amino acid (aa) and non-synonymous substitutions showed that genetic variations resulted in numerously specific codon mutations scattered across the neutralizing epitopes within the ORF5 gene. The PRRSV 1 challenge experiment disclosed the pathogenetic capability of the NPUST2789 isolate in nursery pigs. These findings provide comprehensive knowledge of the molecular diversity of the PRRSV 1 variant in local Taiwanese fields and facilitate the development of suitable immunization programs against this disease.

Genomic characterization and molecular evolution of SARS-CoV-2 in Rio Grande do Sul State, Brazil

SARS-CoV-2 is the virus responsible for the COVID-19 and has afflicted the world since the end of 2019. Different lineages have been discovered and the Gamma lineage, which started the second wave of infections, was first described in Brazil, one of the most affected countries by pandemic. Therefore, this study analyzed SARS-CoV-2 sequenced genomes from Esteio city in Rio Grande do Sul, Southern Brazil. We also comparatively analyzed genomes of the two first years of the pandemic from Rio Grande do Sul state for understanding their genomic and evolutionary patterns. The phylogenomic analysis showed monophyletic groups for Alpha, Gamma, Delta and Omicron, as well as for other circulating lineages in the state. Molecular evolutionary analysis identified several sites under adaptive selection in membrane and nucleocapsid proteins which could be related to a prevalent stabilizing effect on membrane protein structure, as well as majoritarily destabilizing effects on C-terminal nucleocapsid domain.

Natural selection and convergent evolution of the HOX gene family in Carnivora

HOX genes play a central role in the development and regulation of limb patterns. For mammals in the order Carnivora, limbs have evolved in different forms, and there are interesting cases of phenotypic convergence, such as the pseudothumb of the giant and red pandas, and the flippers or specialized limbs of the pinnipeds and sea otter. However, the molecular bases of limb development remain largely unclear. Here, we studied the molecular evolution of the HOX9 ~ 13 genes of 14 representative species in Carnivora and explored the molecular evolution of other HOX genes. We found that only one limb development gene, HOXC10, underwent convergent evolution between giant and red pandas and was thus an important candidate gene related to the development of pseudothumbs. No signals of amino acid convergence and natural selection were found in HOX9 ~ 13 genes between pinnipeds and sea otter, but there was evidence of positive selection and rapid evolution in four pinniped species. Overall, few HOX genes evolve via natural selection or convergent evolution, and these could be important candidate genes for further functional validation. Our findings provide insights into potential molecular mechanisms of the development of specialized pseudothumbs and flippers (or specialized limbs).

Evolution of tetraspanin antigens in the zoonotic Asian blood fluke Schistosoma japonicum

BACKGROUND

Despite successful control efforts in China over the past 60 years, zoonotic schistosomiasis caused by Schistosoma japonicum remains a threat with transmission ongoing and the risk of localised resurgences prompting calls for a novel integrated control strategy, with an anti-schistosome vaccine as a core element. Anti-schistosome vaccine development and immunisation attempts in non-human mammalian host species, intended to interrupt transmission, and utilising various antigen targets, have yielded mixed success, with some studies highlighting variation in schistosome antigen coding genes (ACGs) as possible confounders of vaccine efficacy. Thus, robust selection of target ACGs, including assessment of their genetic diversity and antigenic variability, is paramount. Tetraspanins (TSPs), a family of tegument-surface antigens in schistosomes, interact directly with the host's immune system and are promising vaccine candidates. Here, for the first time to our knowledge, diversity in S. japonicum TSPs (SjTSPs) and the impact of diversifying selection and sequence variation on immunogenicity in these protiens were evaluated.

METHODS

SjTSP sequences, representing parasite populations from seven provinces across China, were gathered by baiting published short-read NGS data and were analysed using in silico methods to measure sequence variation and selection pressures and predict the impact of selection on variation in antigen protein structure, function and antigenic propensity.

RESULTS

Here, 27 SjTSPs were identified across three subfamilies, highlighting the diversity of TSPs in S. japonicum. Considerable variation was demonstrated for several SjTSPs between geographical regions/provinces, revealing that episodic, diversifying positive selection pressures promote amino acid variation/variability in the large extracellular loop (LEL) domain of certain SjTSPs. Accumulating polymorphisms in the LEL domain of SjTSP-2, -8 and -23 led to altered structural, functional and antibody binding characteristics, which are predicted to impact antibody recognition and possibly blunt the host's ability to respond to infection. Such changes, therefore, appear to represent a mechanism utilised by S. japonicum to evade the host's immune system.

CONCLUSION

Whilst the genetic and antigenic geographic variability observed amongst certain SjTSPs could present challenges to vaccine development, here we demonstrate conservation amongst SjTSP-1, -13 and -14, revealing their likely improved utility as efficacious vaccine candidates. Importantly, our data highlight that robust evaluation of vaccine target variability in natural parasite populations should be a prerequisite for anti-schistosome vaccine development.

Analysis of sequence diversity in Plasmodium falciparum glutamic acid-rich protein (PfGARP), an asexual blood stage vaccine candidate

Glutamic acid-rich protein of Plasmodium falciparum (PfGARP) binds to erythrocyte band 3 and may enhance cytoadherence of infected erythrocytes. Naturally acquired anti-PfGARP antibodies could confer protection against high parasitemia and severe symptoms. While whole genome sequencing analysis has suggested high conservation in this locus, little is known about repeat polymorphism in this vaccine candidate antigen. Direct sequencing was performed from the PCR-amplified complete PfGARP gene of 80 clinical isolates from four malaria endemic provinces in Thailand and an isolate from a Guinean patient. Publicly available complete coding sequences of this locus were included for comparative analysis. Six complex repeat (RI-RVI) and two homopolymeric glutamic acid repeat (E1 and E2) domains were identified in PfGARP. The erythrocyte band 3-binding ligand in domain RIV and the epitope for mAB7899 antibody eliciting in vitro parasite killing property were perfectly conserved across isolates. Repeat lengths in domains RIII and E1-RVI-E2 seemed to be correlated with parasite density of the patients. Sequence variation in PfGARP exhibited genetic differentiation across most endemic areas of Thailand. Phylogenetic tree inferred from this locus has shown that most Thai isolates formed closely related lineages, suggesting local expansion/contractions of repeat-encoding regions. Positive selection was observed in non-repeat region preceding domain RII which corresponded to a helper T cell epitope predicted to be recognized by a common HLA class II among Thai population. Predicted linear B cell epitopes were identified in both repeat and non-repeat domains. Besides length variation in some repeat domains, sequence conservation in non-repeat regions and almost all predicted immunogenic epitopes have suggested that PfGARP-derived vaccine may largely elicit strain-transcending immunity.

Emergence, phylogeography, and adaptive evolution of mpox virus

Mpox (Monkeypox) is a zoonotic disease caused by mpox virus (MPXV). A multi-country MPXV outbreak in non-endemic demographics was identified in May 2022. A systematic evaluation of MPXV evolutionary trajectory and genetic diversity could be a timely addition to the MPXV diagnostics and prophylaxis. Herein, we integrated a systematic evolution analysis including phylogenomic and phylogeographic, followed by an in-depth analysis of the adaptive evolution and amino acid variations in type I interferon binding protein (IFNα/βBP). Mutations in IFNα/βBP protein may impair its binding capacity, affecting the MPXV immune evasion strategy. Based on the equilibrated data, we found an evolutionary rate of 7.75 × 10 - 5 substitutions/site/year, and an earlier original time (2021.25) of the clade IIb. We further discovered significant genetic variations in MPXV genomes from different regions and obtained six plausible spread trajectories from its intricate viral flow network, implying that North America might have acted as a bridge for the spread of MPXV from Africa to other continents. We identified two amino acids under positive selection in the Rifampicin resistance protein and extracellular enveloped virus (EEV) type-I membrane glycoprotein, indicating a role in adaptive evolution. Our research sheds light on the emergence, dispersal, and adaptive evolution of MPXV, providing theoretical support for mitigating and containing its expansion.

Proposed taxonomic revision of Giardia duodenalis

Giardia duodenalis, Giardia enterica, Giardia intestinalis and Giardia lamblia are the synonyms for a species complex of 8-11 phylogenetically distinct species of Giardia infecting a broad range of animals including humans. Retrospective alignment of 8409 gene sequences from 3 loci confirmed host associations of Assemblages and sub-Assemblages within this species complex and molecular species delimitation testing confirmed that the Assemblages and sub-Assemblages AI and AII should be recognised as distinct species. It is recommended to synonymise the Assemblages with historic species descriptions based on host associations and consider descriptions for new species where no corresponding description exists. Synonyms, Giardia duodenalis, Giardia intestinalis and Giardia enterica, to be removed from synonymy: synonymise "Giardia duodenalis-Assemblage AI" syn. n. to Giardia duodenalis (Davaine, 1875), Kofoid and Christansen, 1915, synonymise "Giardia duodenalis-Assemblage AII" syn. n. to Giardia intestinalis (Lambl, 1859; Blanchard, 1885), Alexeieff, 1914 and synonymise "Giardia duodenalis-Assemblage B" syn. n. to Giardia enterica (Grassi, 1881a), Kofoid, 1920. Host specific Assemblages synonymised: synonymise canid-associated "Giardia duodenalis-Assemblage C" syn. n. to Giardia canisHegner, 1922; synonymise artiodactyl-associated "Giardia duodenalis-Assemblage E" syn. n. to Giardia bovisFantham, 1921; synonymise feline-associated "Giardia duodenalis-Assemblage F" syn. n. to Giardia catiDeschiens, 1925; and synonymise rodent-associated "Giardia duodenalis-Assemblage G" syn. n. to Giardia simoniLavier, 1924. New description for parasite type infecting specific host: canid-associated "Giardia duodenalis-Assemblage D" named Giardia lupus, sp. n. (LSID: urn:lsid:zoobank.org:act:1651A8CB-CBA8-40D9-AB59-D4AB11AC18A3). New proposed names and descriptions for consideration for parasite types infecting specific hosts: cervid-associated "Giardia duodenalis-sub-Assemblage AIII" for consideration "cervus" and Pinnipedia-associated "Giardia duodenalis-Assemblage H" for consideration "pinnipedis".

Molecular evolution analysis of three species gyroviruses in China from 2018 to 2019

Gyrovirus (GyV) is a widespread ssDNA virus with a high population diversity, and several of its species, including the chicken anemia virus (CAV), gyrovirus galga 1 (GyG1), and gyrovirus homsa 1 (GyH1), have been shown to be pathogenic to poultry. The evolution of these viruses, however, is still unclear. Our study analyzed epidemiology and molecular evolution of three species of GyVs (CAV, GyG1, and GyH1) from 2018 to 2019 in China. The survey results indicated that GyV was widespread in China. It is vital to consider the coinfections among the three species of GyV. The phylogenetic analysis showed that CAV was divided into three clades and GyG1 and GyH1 were divided into two clades. Based on the recombination analysis, CAV and GyG1 had similar recombination regions associated with viral replication and transcription. Furthermore, the substitution rates for CAV and GyG1 were approximately 6.09 × 10^-4 and 2.784 × 10^-4 nucleotides per site per year, respectively. The high substitution rate and recombination were the main factors for the high diversity of GyVs. Unfortunately, GyH1 strains have not been discovered in enough numbers to allow evolutionary analysis. The GyVs had several positively selected sites, possibly related to their potential to escape the host immune response. In summary, our study provides insights into the time of origin, evolution rate, and recombination of GyV for assessing their evolutionary process and genetic diversity.

Pan-Chromosome and Comparative Analysis of Agrobacterium fabrum Reveal Important Traits Concerning the Genetic Diversity, Evolutionary Dynamics, and Niche Adaptation of the Species

Agrobacterium fabrum has been critical for the development of plant genetic engineering and agricultural biotechnology due to its ability to transform eukaryotic cells. However, the gene composition, evolutionary dynamics, and niche adaptation of this species is still unknown. Therefore, we established a comparative genomic analysis based on a pan-chromosome data set to evaluate the genetic diversity of A. fabrum. Here, 25 A. fabrum genomes were selected for analysis by core genome phylogeny combined with the average nucleotide identity (ANI), amino acid identity (AAI), and in silico DNA-DNA hybridization (DDH) values. An open pan-genome of A. fabrum exhibits genetic diversity with variable accessorial genes as evidenced by a consensus pan-genome of 12 representative genomes. The genomic plasticity of A. fabrum is apparent in its putative sequences for mobile genetic elements (MGEs), limited horizontal gene transfer barriers, and potentially horizontally transferred genes. The evolutionary constraints and functional enrichment in the pan-chromosome were measured by the Clusters of Orthologous Groups (COG) categories using eggNOG-mapper software, and the nonsynonymous/synonymous rate ratio (dN/dS) was determined using HYPHY software. Comparative analysis revealed significant differences in the functional enrichment and the degree of purifying selection between the core genome and non-core genome. We demonstrate that the core gene families undergo stronger purifying selection but have a significant bias to contain one or more positively selected sites. Furthermore, although they shared similar genetic diversity, we observed significant differences between chromosome 1 (Chr I) and the chromid in their functional features and evolutionary constraints. We demonstrate that putative genetic elements responsible for plant infection, ecological adaptation, and speciation represent the core genome, highlighting their importance in the adaptation of A. fabrum to plant-related niches. Our pan-chromosome analysis of A. fabrum provides comprehensive insights into the genetic properties, evolutionary patterns, and niche adaptation of the species. IMPORTANCE Agrobacterium spp. live in diverse plant-associated niches such as soil, the rhizosphere, and vegetation, which are challenged by multiple stressors such as diverse energy sources, plant defenses, and microbial competition. They have evolved the ability to utilize diverse resources, escape plant defenses, and defeat competitors. However, the underlying genetic diversity and evolutionary dynamics of Agrobacterium spp. remain unexplored. We examined the phylogeny and pan-genome of A. fabrum to define intraspecies evolutionary relationships. Our results indicate an open pan-genome and numerous MGEs and horizontally transferred genes among A. fabrum genomes, reflecting the flexibility of the chromosomes and the potential for genetic exchange. Furthermore, we observed significant differences in the functional features and evolutionary constraints between the core and accessory genomes and between Chr I and the chromid, respectively.

Machine Learning Prediction and Phyloanatomic Modeling of Viral Neuroadaptive Signatures in the Macaque Model of HIV-Mediated Neuropathology

In human immunodeficiency virus (HIV) infection, virus replication in and adaptation to the central nervous system (CNS) can result in neurocognitive deficits in approximately 25% of patients with unsuppressed viremia. While no single viral mutation can be agreed upon as distinguishing the neuroadapted population, earlier studies have demonstrated that a machine learning (ML) approach could be applied to identify a collection of mutational signatures within the virus envelope glycoprotein (Gp120) predictive of disease. The S[imian]IV-infected macaque is a widely used animal model of HIV neuropathology, allowing in-depth tissue sampling infeasible for human patients. Yet, translational impact of the ML approach within the context of the macaque model has not been tested, much less the capacity for early prediction in other, noninvasive tissues. We applied the previously described ML approach to prediction of SIV-mediated encephalitis (SIVE) using gp120 sequences obtained from the CNS of animals with and without SIVE with 97% accuracy. The presence of SIVE signatures at earlier time points of infection in non-CNS tissues indicated these signatures cannot be used in a clinical setting; however, combined with protein structural mapping and statistical phylogenetic inference, results revealed common denominators associated with these signatures, including 2-acetamido-2-deoxy-beta-d-glucopyranose structural interactions and high rate of alveolar macrophage (AM) infection. AMs were also determined to be the phyloanatomic source of cranial virus in SIVE animals, but not in animals that did not develop SIVE, implicating a role for these cells in the evolution of the signatures identified as predictive of both HIV and SIV neuropathology. IMPORTANCE HIV-associated neurocognitive disorders remain prevalent among persons living with HIV (PLWH) owing to our limited understanding of the contributing viral mechanisms and ability to predict disease onset. We have expanded on a machine learning method previously used on HIV genetic sequence data to predict neurocognitive impairment in PLWH to the more extensively sampled SIV-infected macaque model in order to (i) determine the translatability of the animal model and (ii) more accurately characterize the predictive capacity of the method. We identified eight amino acid and/or biochemical signatures in the SIV envelope glycoprotein, the most predominant of which demonstrated the potential for aminoglycan interaction characteristic of previously identified HIV signatures. These signatures were not isolated to specific points in time or to the central nervous system, limiting their use as an accurate clinical predictor of neuropathogenesis; however, statistical phylogenetic and signature pattern analyses implicate the lungs as a key player in the emergence of neuroadapted viruses.

Genomic coevolution of papillomavirus and immune system in placental mammals indicates the role of IFN-γ in the emergence of new variants

Papillomaviruses (PVs) are causative agents for warts and cancers in different parts of the body in the mammalian lineage. Therefore, these viruses are proposed as model organisms to study host immune responses to pathogens causing chronic infections. The virus-associated cancer progression depends on two integral processes namely angiogenesis and immune response (AIR). The angiogenesis process aids in tumour progression through vessel formation and maturation but the host immune response, in contrast, makes every attempt to eliminate pathogens and thereby maintain healthy tissues. However, the evolutionary contribution of individual viral genes and host AIR genes in carcinogenesis is yet to be explored. Here, we applied the evolutionary genomics approach to find correlated evolution between six PV genes and 23 host AIR-related genes. We estimated that IFN-γ is the only host gene evolving in a correlated manner with all six PV genes under study. Furthermore, three papillomavirus genes, L2, E6, and E7, are found to interact with two third of host AIR-related genes. Moreover, a combined differential gene expression analysis and network analysis showed that inflammatory cytokine IFN-γ is the key regulator of hub genes in the PPI network of the differentially expressed genes. Functional enrichment of these hub genes is consistent with their established role in different cancers and viral infections. Overall, we conclude that IFN-γ maintains selective pressure on mammalian PV genes and seems to be a potential biomarker for PV-related cancers. This study demonstrates the evolutionary importance of IFN-γ in deciding the fate of carcinogenic PV variants.

HexSE: Simulating evolution in overlapping reading frames

Gene overlap occurs when two or more genes are encoded by the same nucleotides. This phenomenon is found in all taxonomic domains, but is particularly common in viruses, where it may provide a mechanism to increase the information content of compact genomes. The presence of overlapping reading frames (OvRFs) can skew estimates of selection based on the rates of non-synonymous and synonymous substitutions, since a substitution that is synonymous in one reading frame may be non-synonymous in another and vice versa. To understand the impact of OvRFs on molecular evolution, we implemented a versatile simulation model of nucleotide sequence evolution along a phylogeny with any distribution of open reading frames in linear or circular genomes. We use a custom data structure to track the substitution rates at every nucleotide site, which is determined by the stationary nucleotide frequencies, transition bias and the distribution of selection biases (dN/dS) in the respective reading frames. Our simulation model is implemented in the Python scripting language. All source code is released under the GNU General Public License version 3 and are available at https://github.com/PoonLab/HexSE.

Genomic Epidemiology and Transmission Dynamics of Global Coxsackievirus B4

The aim of this study was to determine the global genetic diversity and transmission dynamics of coxsackievirus B4 (CVB4) and to propose future directions for disease surveillance. Next-generation sequencing was performed to obtain the complete genome sequence of CVB4, and the genetic diversity and transmission dynamics of CVB4 worldwide were analyzed using bioinformatics methods such as phylogenetic analysis, evolutionary dynamics, and phylogeographic analysis. Forty complete genomes of CVB4 were identified from asymptomatic infected individuals and hand, foot, and mouth disease (HFMD) patients. Frequent recombination between CVB4 and EV-B multiple serotypes in the 3D^pol region was found and formed 12 recombinant patterns (A-L). Among these, the CVB4 isolated from asymptomatic infected persons and HFMD patients belonged to lineages H and I, respectively. Transmission dynamics analysis based on the VP1 region revealed that CVB4 epidemics in countries outside China were dominated by the D genotype, whereas the E genotype was dominant in China, and both genotypes evolved at a rate of > 6.50 × 10^-3 substitutions/site/year. CVB4 spreads through the population unseen, with the risk of disease outbreaks persisting as susceptible individuals accumulate. Our findings add to publicly available CVB4 genomic sequence data and deepen our understanding of CVB4 molecular epidemiology.

That H9N2 avian influenza viruses circulating in different regions gather in the same live-poultry market poses a potential threat to public health

H9N2 avian influenza viruses are endemic and persistent in China, but those that are prevalent in different provinces are also causes of wide epidemics, related to the spread of wild birds and the cross-regional trade in live poultry. For the past 4 years, beginning in 2018, we have sampled a live-poultry market in Foshan, Guangdong, in this ongoing study. In addition to the prevalence of H9N2 avian influenza viruses in China during this period, we identified isolates from the same market belonging to clade A and clade B, which diverged in 2012-2013, and clade C, which diverged in 2014-2016, respectively. An analysis of population dynamics revealed that, after a critical divergence period from 2014 to 2016, the genetic diversity of H9N2 viruses peaked in 2017. Our spatiotemporal dynamics analysis found that clade A, B, and C, which maintain high rates of evolution, have different prevalence ranges and transmission paths. Clades A and B were mainly prevalent in East China in the early stage, and then spread to Southern China, becoming epidemic with clade C. Strains from different regions converge at the same live-poultry market to communicate, which may be one reasons the H9N2 viruses are difficult to eradicate and increasingly dominant throughout China. Selection pressure and molecular analysis have demonstrated that single amino acid polymorphisms at key receptor binding sites 156, 160, and 190 under positive selection pressure, suggesting that H9N2 viruses are undergoing mutations to adapt to new hosts. Live-poultry markets are important because people who visit them have frequent contact with poultry, H9N2 viruses from different regions converge at these markets and spread through contact between live birds and humans, generating increased risks of human exposure to these viruses and threatening public health safety. Thus, it is important to reducing the cross-regional trade of live poultry and strengthening the monitoring of avian influenza viruses in live-poultry markets to reduce the spread of avian influenza viruses.

Evidence of natural selection and dominance of SARS-CoV-2 variant Lambda (C.37) over variants of concern in Cusco, Peru

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage C.37 (Lambda) has spread rapidly in Peru and other Latin American countries. However, most studies in Peru have focused on Lima, the capital city, without knowing the dynamics of the spread of the variant in other departments. Cusco, Peru, is one of the most popular departments in the country for tourists, so the introduction of new variants of SARS-CoV-2 might occur despite closure of the borders. Therefore, in this work, we analyzed the variants circulating in Cusco. The aim of this work was to better understand the distribution of SARS-CoV-2 lineages circulating in Cusco and to characterize the genomes of these strains. To this end, 46 SARS-CoV-2 genomes from vaccinated and unvaccinated patients were sequenced in the first half of 2021. The genomes were analyzed using phylogenetic and natural selection methods. Phylogenetic trees from Cusco showed dominance of the Lambda lineage over the variants of concern (VOCs), and there was no clustering of variants by district. Natural selection analysis revealed mutations, mainly in the spike protein, at positions 75, 246, 247, 707, 769, and 1020. In addition, we found that unvaccinated patients accumulated more new mutations than did vaccinated patients, and these included the F101Y mutation in ORF7a, E419A in NSP3, a deletion in S (21,618-22,501), and a deletion in ORF3a (25,437-26,122).

Genomic Analysis of G2P[4] Group A Rotaviruses in Zambia Reveals Positive Selection in Amino Acid Site 7 of Viral Protein 3

The G2P[4] genotype is among the rotavirus strains that circulate commonly in humans. Several countries have reported its immediate upsurge after the introduction of rotavirus vaccination, raising concern about sub-optimal vaccine effectiveness against this genotype in the long term. This study aimed to gain insight into the evolution of post-vaccine Zambian G2P[4] group A rotavirus (RVA) strains and their overall genetic make-up by analysis of sequence alignments at the amino acid (AA) level. Twenty-nine Zambian G2P[4] rotavirus strains were subjected to whole-genome sequencing using the Illumina MiSeq® platform. All the strains exhibited the typical DS-1-like genotype constellation, and the nucleotide sequences of the 11 genome segments showed high nucleotide similarities (>97%). Phylogenetic analyses together with representative global G2P[4] RVA showed that Zambian strains clustered into human lineages IV (for VP2, VP4, VP7, NSP1, and NSP5), V (for VP1, VP3, VP6, NSP2, and NSP3), and XXIII (for NSP4). The AA differences between the lineages where the study strains clustered and lineages of global reference strains were identified and analyzed. Selection pressure analysis revealed that AA site seven in the Viral Protein 3 (VP3) genome segment was under positive selection. This site occurs in the region of intrinsic disorder in the VP3 protein, and Zambian G2P[4] strains could potentially be utilizing this intrinsically disordered region to survive immune pressure. The Zambian G2P[4] strains from 2012 to 2016 comprised the G2P[4] strains that have been circulating globally since the early 2000s, highlighting the epidemiological fitness of these contemporary G2P[4] strains. Continuous whole-genome surveillance of G2P[4] strains remains imperative to understand their evolution during the post-vaccination period.

The deep-rooted origin of disulfide-rich spider venom toxins

Spider venoms are a complex concoction of enzymes, polyamines, inorganic salts, and disulfide-rich peptides (DRPs). Although DRPs are widely distributed and abundant, their bevolutionary origin has remained elusive. This knowledge gap stems from the extensive molecular divergence of DRPs and a lack of sequence and structural data from diverse lineages. By evaluating DRPs under a comprehensive phylogenetic, structural and evolutionary framework, we have not only identified 78 novel spider toxin superfamilies but also provided the first evidence for their common origin. We trace the origin of these toxin superfamilies to a primordial knot - which we name 'Adi Shakti', after the creator of the Universe according to Hindu mythology - 375 MYA in the common ancestor of Araneomorphae and Mygalomorphae. As the lineages under evaluation constitute nearly 60% of extant spiders, our findings provide fascinating insights into the early evolution and diversification of the spider venom arsenal. Reliance on a single molecular toxin scaffold by nearly all spiders is in complete contrast to most other venomous animals that have recruited into their venoms diverse toxins with independent origins. By comparatively evaluating the molecular evolutionary histories of araneomorph and mygalomorph spider venom toxins, we highlight their contrasting evolutionary diversification rates. Our results also suggest that venom deployment (e.g. prey capture or self-defense) influences evolutionary diversification of DRP toxin superfamilies.

Evaluating the Potential Fitness Effects of Chinook Salmon (Oncorhynchus tshawytscha) Aquaculture Using Non-Invasive Population Genomic Analyses of MHC Nucleotide Substitution Spectra

Genetic diversity plays a vital role in the adaptability of salmon to changing environmental conditions that can introduce new selective pressures on populations. Variability among local subpopulations may increase the chance that certain advantageous genes are passed down to future generations to mitigate susceptibility to novel diseases, warming oceans, loss of genetic stocks, and ocean acidification. Class I and II genes of the major histocompatibility complex (MHC) are crucial for the fitness of Chinook salmon due to the role they play in disease and pathogen resistance. The objective of this study was to assess the DNA sequence variability among wild and hatchery populations of Alaskan Chinook salmon at the class I α1 and class II β1 exons of the MHC. We hypothesized that the 96 wild samples taken from the Deshka River would display greater levels of observed heterozygosity (Ho) relative to expected heterozygosity (He) in suggesting that individuals with similar phenotypes mate with one another more frequently than would be expected under random mating patterns. Conversely, since no mate selection occurs in the William Jack Hernandez Sport Fish hatchery, we would not expect to see this discrepancy (He = Ho) in the 96 hatchery fish tested in this study. Alternatively, we hypothesized that post-mating selection is driving higher levels of observed heterozygosity as opposed to mate selection. If this is the case, we will observe higher than expected levels of heterozygosity among hatchery salmon. Both populations displayed higher levels of observed heterozygosity than expected heterozygosity at the Class I and II loci but genetic differentiation between the spatially distinct communities was minimal. Class I sequences showed evidence of balancing selection, despite high rates of non-synonymous substitutions observed, specifically at the peptide binding regions of both MHC genes.

The evolution, characterization and phylogeography of avian influenza H9N2 viruses from India

The low pathogenic avian influenza H9N2 virus is a significant zoonotic agent and contributes genes to highly pathogenic avian influenza (HPAI) viruses. H9N2 viruses are prevalent in India with a reported human case. We elucidate the spatio-temporal origins of the H9N2 viruses from India. A total of 30H9N2 viruses were isolated from poultry and environmental specimens (years 2015-2020). Genome sequences of H9N2 viruses (2003-2020) from India were analyzed, revealing several substitutions. We found five reassortant genotypes. The HA, NA and PB2 genes belonged to the Middle-Eastern B sublineage; NP and M to the classical G1 lineage; PB1, PA and NS showed resemblance to genes from either HPAI-H7N3/H5N1 viruses. Molecular clock and phylogeography revealed that the introduction of all the genes to India took place around the year 2000. This is the first report of the genesis and evolution of the H9N2 viruses from India, and highlights the need for surveillance.

Early Genomic Surveillance and Phylogeographic Analysis of Getah Virus, a Reemerging Arbovirus, in Livestock in China

Getah virus (GETV) mainly causes disease in livestock and may pose an epidemic risk due to its expanding host range and the potential of long-distance dispersal through animal trade. Here, we used metagenomic next-generation sequencing (mNGS) to identify GETV as the pathogen responsible for reemerging swine disease in China and subsequently estimated key epidemiological parameters using phylodynamic and spatially-explicit phylogeographic approaches. The GETV isolates were able to replicate in a variety of cell lines, including human cells, and showed high pathogenicity in a mouse model, suggesting the potential for more mammal hosts. We obtained 16 complete genomes and 79 E2 gene sequences from viral strains collected in China from 2016 to 2021 through large-scale surveillance among livestock, pets, and mosquitoes. Our phylogenetic analysis revealed that three major GETV lineages are responsible for the current epidemic in livestock in China. We identified three potential positively selected sites and mutations of interest in E2, which may impact the transmissibility and pathogenicity of the virus. Phylodynamic inference of the GETV demographic dynamics identified an association between livestock meat consumption and the evolution of viral genetic diversity. Finally, phylogeographic reconstruction of GETV dispersal indicated that the sampled lineages have preferentially circulated within areas associated with relatively higher mean annual temperature and pig population density. Our results highlight the importance of continuous surveillance of GETV among livestock in southern Chinese regions associated with relatively high temperatures. IMPORTANCE Although livestock is known to be the primary reservoir of Getah virus (GETV) in Asian countries, where identification is largely based on serology, the evolutionary history and spatial epidemiology of GETV in these regions remain largely unknown. Through our sequencing efforts, we provided robust support for lineage delineation of GETV and identified three major lineages that are responsible for the current epidemic in livestock in China. We further analyzed genomic and epidemiological data to reconstruct the recent demographic and dispersal history of GETV in domestic animals in China and to explore the impact of environmental factors on its genetic diversity and its diffusion. Notably, except for livestock meat consumption, other pig-related factors such as the evolution of live pig transport and pork production do not show a significant association with the evolution of viral genetic diversity, pointing out that further studies should investigate the potential contribution of other host species to the GETV outbreak. Our analysis of GETV demonstrates the need for wider animal species surveillance and provides a baseline for future studies of the molecular epidemiology and early warning of emerging arboviruses in China.

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

The molecular evolution of genes previously associated with large sizes reveals possible pathways to cetacean gigantism

Cetaceans are a group of aquatic mammals with the largest body sizes among living animals, including giant representatives such as blue and fin whales. To understand the genetic bases of gigantism in cetaceans, we performed molecular evolutionary analyses on five genes (GHSR, IGF2, IGFBP2, IGFBP7, and EGF) from the growth hormone/insulin-like growth factor axis, and four genes (ZFAT, EGF, LCORL, and PLAG1) previously described as related to the size of species evolutionarily close to cetaceans, such as pigs, cows, and sheep. Our dataset comprised 19 species of cetaceans, seven of which are classified as giants because they exceed 10 m in length. Our results revealed signs of positive selection in genes from the growth hormone/insulin-like growth factor axis and also in those related to body increase in cetacean-related species. In addition, pseudogenization of the EGF gene was detected in the lineage of toothless cetaceans, Mysticeti. Our results suggest the action of positive selection on gigantism in genes that act both in body augmentation and in mitigating its consequences, such as cancer suppression when involved in processes such as division, migration, and cell development control.

Evolutionary dissection of monkeypox virus: Positive Darwinian selection drives the adaptation of virus-host interaction proteins

The emerging and ongoing outbreak of human monkeypox (hMPX) in 2022 is a serious global threat. An understanding of the evolution of the monkeypox virus (MPXV) at the single-gene level may provide clues for exploring the unique aspects of the current outbreak: rapidly expanding and sustained human-to-human transmission. For the current investigation, alleles of 156 MPXV coding genes (which account for >95% of the genomic sequence) have been gathered from roughly 1,500 isolates, including those responsible for the previous outbreaks. Using a range of molecular evolution approaches, we demonstrated that intra-species homologous recombination has a negligible effect on MPXV evolution. Despite the fact that the majority of the MPXV genes (64.10%) were subjected to negative selection at the whole gene level, 10 MPXV coding genes (MPXVgp004, 010, 012, 014, 044, 098, 138, 178, 188, and 191) were found to have a total of 15 codons or amino acid sites that are known to evolve under positive Darwinian selection. Except for MPXVgp138, almost all of these genes encode proteins that interact with the host. Of these, five ankyrin proteins (MPXVgp004, 010, 012, 178, and 188) and one Bcl-2-like protein (MPXVgp014) are involved in poxviruses' host range determination. We discovered that the majority (80%) of positive amino acid substitutions emerged several decades ago, indicating that these sites have been under constant selection pressure and that more adaptable alleles have been circulating in the natural reservoir. This finding was also supported by the minimum spanning networks of the gene alleles. The three positive amino acid substitutions (T/A426V in MPXVgp010, A423D in MPXVgp012, and S105L in MPXVgp191) appeared in 2019 or 2022, indicating that they would be crucial for the virus' eventual adaptation to humans. Protein modeling suggests that positive amino acid substitutions may affect protein functions in a variety of ways. Further study should focus on revealing the biological effects of positive amino acid substitutions in the genes for viral adaptation to humans, virulence, transmission, and so on. Our study advances knowledge of MPXV's adaptive mechanism and provides insights for exploring factors that are responsible for the unique aspects of the current outbreak.

Diversifying selection and climatic effects on major histocompatibility complex class II gene diversity in the greater horseshoe bat

Heterogeneous pathogenic stress can shape major histocompatibility complex (MHC) diversity by influencing the functional plasticity of the immune response. Therefore, MHC diversity could reflect environmental stress, demonstrating its importance in uncovering the mechanisms of adaptive genetic variation. In this study, we combined neutral microsatellite loci, an immune-related MHC II-DRB locus, and climatic factors to unravel the mechanisms affecting the diversity and genetic differentiation of MHC genes in the greater horseshoe bat (Rhinolophus ferrumequinum), a species with a wide geographical distribution that has three distinct genetic lineages in China. First, increased genetic differentiation at the MHC locus among populations compared using microsatellites indicated diversifying selection. Second, the genetic differentiation of MHC and microsatellites were significantly correlated, suggesting that demographic processes exist. However, MHC genetic differentiation was significantly correlated with geographical distance among populations, even after controlling for the neutral markers, suggesting a major effect of selection. Third, although the MHC genetic differentiation was larger than that for microsatellites, there was no significant difference in the genetic differentiation between the two markers among genetic lineages, indicating the effect of balancing selection. Fourth, combined with climatic factors, MHC diversity and supertypes showed significant correlations with temperature and precipitation, but not with the phylogeographic structure of R. ferrumequinum, suggesting an effect of local adaptation driven by climate on MHC diversity. Moreover, the number of MHC supertypes varied between populations and lineages, suggesting regional characteristics and support for local adaptation. Taken together, the results of our study provide insights into the adaptive evolutionary driving forces at different geographic scales in R. ferrumequinum. In addition, climate factors may have played a vital role in driving adaptive evolution in this species.

Adaptive Evolution of the OAS Gene Family Provides New Insights into the Antiviral Ability of Laurasiatherian Mammals

Many mammals risk damage from virus invasion due to frequent environmental changes. The oligoadenylate synthesis (OAS) gene family, which is an important component of the immune system, provides an essential response to the antiviral activities of interferons by regulating immune signal pathways. However, little is known about the evolutionary characteristics of OASs in Laurasiatherian mammals. Here, we examined the evolution of the OAS genes in 64 mammals to explore the accompanying molecular mechanisms of the antiviral ability of Laurasiatherian mammals living in different environments. We found that OAS2 and OAS3 were found to be pseudogenes in Odontoceti species. This may be related to the fact that they live in water. Some Antilopinae, Caprinae, and Cervidae species lacked the OASL gene, which may be related to their habitats being at higher altitudes. The OASs had a high number of positive selection sites in Cetartiodactyla, which drove the expression of strong antiviral ability. The OAS gene family evolved in Laurasiatherian mammals at different rates and was highly correlated with the species' antiviral ability. The gene evolution rate in Cetartiodactyla was significantly higher than that in the other orders. Compared to other species of the Carnivora family, the higher selection pressure on the OAS gene and the absence of positive selection sites in Canidae may be responsible for its weak resistance to rabies virus. The OAS gene family was relatively conserved during evolution. Conserved genes are able to provide better maintenance of gene function. The rate of gene evolution and the number of positively selected sites combine to influence the resistance of a species to viruses. The positive selection sites demonstrate the adaptive evolution of the OAS gene family to the environment. Adaptive evolution combined with conserved gene function improves resistance to viruses. Our findings offer insights into the molecular and functional evolution of the antiviral ability of Laurasian mammals.

Wide distribution of D-xylose dehydrogenase in yeasts reveals a new element in the D-xylose metabolism for bioethanol production

D-xylose utilization by yeasts is an essential feature for improving second-generation ethanol production. However, industrial yeast strains are incapable of consuming D-xylose. Previous analyzes of D-xylose-consuming or fermenting yeast species reveal that the genomic features associated with this phenotype are complex and still not fully understood. Here we present a previously neglected yeast enzyme related to D-xylose metabolism, D-xylose dehydrogenase (XylDH), which is found in at least 105 yeast genomes. By analyzing the XylDH gene family, we brought evidence of gene evolution marked by purifying selection on codons and positive selection evidence in D-xylose-consuming and fermenting species, suggesting the importance of XylDH for D-xylose-related phenotypes in yeasts. Furthermore, although we found no putative metabolic pathway for XylDH in yeast genomes, namely the absence of three bacterial known pathways for this enzyme, we also provide its expression profile on D-xylose media following D-xylose reductase for two yeasts with publicly available transcriptomes. Based on these results, we suggest that XylDH plays an important role in D-xylose usage by yeasts, likely being involved in a cofactor regeneration system by reducing cofactor imbalance in the D-xylose reductase pathway.

Cytokine gene polymorphism and parasite susceptibility in free-living rodents: Importance of non-coding variants

Associations between genetic variants and susceptibility to infections have long been studied in free-living hosts so as to infer the contemporary evolutionary forces that shape the genetic polymorphisms of immunity genes. Despite extensive studies of proteins interacting with pathogen-derived ligands, such as MHC (major histocompatilbility complex) or TLR (Toll-like receptors), little is known about the efferent arm of the immune system. Cytokines are signalling molecules that trigger and modulate the immune response, acting as a crucial link between innate and adaptive immunity. In the present study we investigated how genetic variation in cytokines in bank voles Myodes glareolus affects their susceptibility to infection by parasites (nematodes: Aspiculuris tianjensis, Heligmosomum mixtum, Heligmosomoides glareoli) and microparasites (Cryptosporidium sp, Babesia microti, Bartonella sp.). We focused on three cytokines: tumour necrosis factor (TNF), lymphotoxin alpha (LTα), and interferon beta (IFNβ1). Overall, we identified four single nucleotide polymorphisms (SNPs) associated with susceptibility to nematodes: two located in LTα and two in IFNβ1. One of those variants was synonymous, another located in an intron. Each SNP associated with parasite load was located in or next to a codon under selection, three codons displayed signatures of positive selection, and one of purifying selection. Our results indicate that cytokines are prone to parasite-driven selection and that non-coding variants, although commonly disregarded in studies of the genetic background of host-parasite co-evolution, may play a role in susceptibility to infections in wild systems.

The epidemiology and evolutionary dynamics of massive dengue outbreak in China, 2019

Introduction

In 2019, China experienced massive dengue outbreaks with high incidence and expanded outbreak areas. The study aims to depict dengue's epidemiology and evolutionary dynamics in China and explore the possible origin of these outbreaks.

Methods

Records of confirmed dengue cases in 2019 were obtained from the China Notifiable Disease Surveillance System. The sequences of complete envelope gene detected from the outbreak provinces in China in 2019 were retrieved from GenBank. Maximum Likelihood trees were constructed to genotype the viruses. The median-joining network was used to visualize fine-scale genetic relationships. Four methods were used to estimate the selective pressure.

Results

A total of 22,688 dengue cases were reported, 71.4% of which were indigenous cases and 28.6% were imported cases (including from abroad and from other domestic provinces). The abroad cases were predominantly imported from Southeast Asia countries (94.6%), with Cambodia (3,234 cases, 58.9%), and Myanmar (1,097 cases, 20.0%) ranked as the top two. A total of 11 provinces with dengue outbreaks were identified in the central-south of China, of which Yunnan and Guangdong provinces had the highest number of imported and indigenous cases. The primary source of imported cases in Yunnan was from Myanmar, while in the other ten provinces, the majority of imported cases were from Cambodia. Guangdong, Yunnan and Guangxi provinces were China's primary sources of domestically imported cases. Phylogenetic analysis of the viruses in outbreak provinces revealed three genotypes: (I, IV, and V) in DENV 1, Cosmopolitan and Asian I genotypes in DENV 2, and two genotypes (I and III) in DENV 3. Some genotypes concurrently circulated in different outbreak provinces. Most of the viruses were clustered with those from Southeast Asia. Haplotype network analysis showed that Southeast Asia, possibly Cambodia and Thailand, was the respective origin of the viruses in clade 1 and 4 for DENV 1. Positive selection was detected at codon 386 in clade 1.

Conclusion

Dengue importation from abroad, especially from Southeast Asia, resulted in the dengue epidemic in China in 2019. Domestic transmission between provinces and positive selection on virus evolution may contribute to the massive dengue outbreaks.

Effect of population size and selection on Toll-like receptor diversity in populations of Galápagos mockingbirds

The interactions of evolutionary forces are difficult to analyse in free-living populations. However, when properly understood, they provide valuable insights into evolutionary biology and conservation genetics. This is particularly important for the interplay of genetic drift and natural selection in immune genes that confer resistance to disease. The Galápagos Islands are inhabited by four closely related species of mockingbirds (Mimus spp.). We used 12 different-sized populations of Galápagos mockingbirds and one population of their continental relative northern mockingbird (Mimus polyglottos) to study the effects of genetic drift on the molecular evolution of immune genes, the Toll-like receptors (TLRs: TLR1B, TLR4 and TLR15). We found that neutral genetic diversity was positively correlated with island size, indicating an important effect of genetic drift. However, for TLR1B and TLR4, there was little correlation between functional (e.g., protein) diversity and island size, and protein structural properties were largely conserved, indicating only a limited effect of genetic drift on molecular phenotype. By contrast, TLR15 was less conserved and even its putative functional polymorphism correlated with island size. The patterns observed for the three genes suggest that genetic drift does not necessarily dominate selection even in relatively small populations, but that the final outcome depends on the degree of selection constraint that is specific for each TLR locus.

Transposon-derived transcription factors across metazoans

Transposable elements (TE) could serve as sources of new transcription factors (TFs) in plants and some other model species, but such evidence is lacking for most animal lineages. Here, we discovered multiple independent co-options of TEs to generate 788 TFs across Metazoa, including all early-branching animal lineages. Six of ten superfamilies of DNA transposon-derived conserved TF families (ZBED, CENPB, FHY3, HTH-Psq, THAP, and FLYWCH) were identified across nine phyla encompassing the entire metazoan phylogeny. The most extensive convergent domestication of potentially TE-derived TFs occurred in the hydroid polyps, polychaete worms, cephalopods, oysters, and sea slugs. Phylogenetic reconstructions showed species-specific clustering and lineage-specific expansion; none of the identified TE-derived TFs revealed homologs in their closest neighbors. Together, our study established a framework for categorizing TE-derived TFs and informing the origins of novel genes across phyla.

Detecting signals of adaptive evolution in grape plastomes with a focus on the Cretaceous-Palaeogene (K/Pg) transition

BACKGROUND AND AIMS

Although plastid genes are widely used in phylogenetic studies, signals of positive selection have been scarcely investigated in the grape family. The plastomes from 91 accessions of Vitaceae were examined to understand the extent to which positive selection is present and to identify which genes are involved. Moreover, the changes through time of genes under episodic positive selection were investigated and the hypothesis of an adaptive process following the Cretaceous-Palaeogene (K/Pg) transition about 66 million years ago was tested.

METHODS

Different codon-substitution models were used to assess pervasive and episodic positive selection events on 70 candidate plastid genes. Divergence times between lineages were estimated and stochastic character mapping analysis was used to simulate variation over time of the genes found to be under episodic positive selection.

KEY RESULTS

A total of 20 plastid genes (29 %) showed positive selection. Among them, 14 genes showed pervasive signatures of positive selection and nine genes showed episodic signatures of positive selection. In particular, four of the nine genes (psbK, rpl20, rpoB, rps11) exhibited a similar pattern showing an increase in the rate of variation close to the K/Pg transition.

CONCLUSION

Multiple analyses have shown that the grape family has experienced ancient and recent positive selection events and that the targeted genes are involved in essential functions such as photosynthesis, self-replication and metabolism. Our results are consistent with the idea that the K/Pg transition has favoured an increased rate of change in some genes. Intense environmental perturbations have influenced the rapid diversification of certain lineages, and new mutations arising on some plastid genes may have been fixed by natural selection over the course of many generations.

RNA structure-altering mutations underlying positive selection on Spike protein reveal novel putative signatures to trace crossing host-species barriers in Betacoronavirus

Similar to other RNA viruses, the emergence of Betacoronavirus relies on cross-species viral transmission, which requires careful health surveillance monitoring of protein-coding information as well as genome-wide analysis. Although the evolutionary jump from natural reservoirs to humans may be mainly traced-back by studying the effect that hotspot mutations have on viral proteins, it is largely unexplored if other impacts might emerge on the structured RNA genome of Betacoronavirus. In this survey, the protein-coding and viral genome architecture were simultaneously studied to uncover novel insights into cross-species horizontal transmission events. We analysed 1,252,952 viral genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 distributed across the world in bats, intermediate animals, and humans to build a new landscape of changes in the RNA viral genome. Phylogenetic analyses suggest that bat viruses are the most closely related to the time of most recent common ancestor of Betacoronavirus, and missense mutations in viral proteins, mainly in the S protein S1 subunit: SARS-CoV (G > T; A577S); MERS-CoV (C > T; S746R and C > T; N762A); and SARS-CoV-2 (A > G; D614G) appear to have driven viral diversification. We also found that codon sites under positive selection on S protein overlap with non-compensatory mutations that disrupt secondary RNA structures in the RNA genome complement. These findings provide pivotal factors that might be underlying the eventual jumping the species barrier from bats to intermediate hosts. Lastly, we discovered that nearly half of the Betacoronavirus genomes carry highly conserved RNA structures, and more than 90% of these RNA structures show negative selection signals, suggesting essential functions in the biology of Betacoronavirus that have not been investigated to date. Further research is needed on negatively selected RNA structures to scan for emerging functions like the potential of coding virus-derived small RNAs and to develop new candidate antiviral therapeutic strategies.

Genomic analysis of hypoxia inducible factor alpha in ray-finned fishes reveals missing Ohnologs and evidence of widespread positive selection

As aquatic hypoxia worsens on a global scale, fishes will become increasingly challenged by low oxygen, and understanding the molecular basis of their response to hypoxia may help to better define the capacity of fishes to cope with this challenge. The hypoxia inducible factor (HIF) plays a critical role in the molecular response to hypoxia by activating the transcription of genes that serve to improve oxygen delivery to the tissues or enhance the capacity of tissues to function at low oxygen. The current study examines the molecular evolution of genes encoding the oxygen-dependent HIFα subunit (HIFA) in the ray-finned fishes (Actinopterygii). Genomic analyses demonstrate that several lineages retain four paralogs of HIFA predicted from two rounds of genome duplication at the base of vertebrate evolution, broaden the known distribution of teleost-specific HIFA paralogs, and provide evidence for salmonid-specific HIFA duplicates. Evolution of the HIFA gene family is characterized by widespread episodic positive selection at amino acid sites that potentially mediate protein stability, protein-protein interactions, and transcriptional regulation. HIFA transcript abundance depends upon paralog, tissue, and fish lineage. A phylogenetically-informed gene nomenclature is proposed along with avenues for future research on this critical family of transcription factors.

The Evolution and Global Spatiotemporal Dynamics of Senecavirus A

Recurrent outbreaks of senecavirus A (SVA)-associated vesicular disease have led to a large number of infected pigs being culled and has caused considerable economic losses to the swine industry. Although SVA was discovered 2 decades ago, knowledge about the evolutionary and transmission histories of SVA remains unclear. Herein, we performed an integrated analysis of the recombination, phylogeny, selection, and spatiotemporal dynamics of SVA. Phylogenetic analysis demonstrated that SVA diverged into two main branches, clade I (pre-2007 strains) and clade II (post-2007 strains). Importantly, analysis of selective strength showed that clade II was evolving under relaxed selection compared with clade I. Positive selection analysis identified 27 positive selective sites, most of which are located on the outer surface of capsid protomer or on the important functional domains of nonstructure proteins. Bayesian phylodynamics suggested that the estimated time to the most recent common ancestor of SVA was around 1986, and the estimated substitution rate of SVA was 3.3522 × 10-3 nucleotide substitutions/site/year. Demographic history analysis revealed that the effective population size of SVA has experienced a gradually increasing trend with slight fluctuation until 2017 followed by a sharp decline. Notably, Bayesian phylogeographic analysis inferred that Brazil might be the source of SVA's global transmission since 2015. In summary, these data illustrated that the ongoing evolution of SVA drove the lineage-specific innovation and potentially phenotypically important variation. Our study sheds new light on the fundamental understanding of SVA evolution and spread history. IMPORTANCE Recurrent outbreaks and global epidemics of senecavirus A-associated vesicular disease have caused heavy economic losses and have threatened the development of the pig industry. However, the question of where the virus comes from has been one of the biggest puzzles due to the stealthy nature of the virus. Consequently, tracing the source, evolution, and transmission pattern of SVA is a very challenging task. Based on the most comprehensive analysis, we revealed the origin time, rapid evolution, epidemic dynamics, and selection of SVA. We observed two main genetic branches, clade I (pre-2007 strains) and clade II (post-2007 strains), and described the epidemiological patterns of SVA in different countries. We also first identified Brazil as the source of SVA's global transmission since 2015. Findings in this study provide important implications for the control and prevention of the virus.

Molecular Evolution of the Bactericidal/Permeability-Increasing Protein (BPIFA1) Regulating the Innate Immune Responses in Mammals

Bactericidal/permeability-increasing protein, a primary factor of the innate immune system of mammals, participates in natural immune protection against invading bacteria. BPIFA1 actively contributes to host defense via multiple mechanisms, such as antibacterial, surfactant, airway surface liquid control, and immunomodulatory activities. However, the evolutionary history and selection forces on the BPIFA1 gene in mammals during adaptive evolution are poorly understood. This study examined the BPIFA1 gene of humans compared with that of other mammalian species to estimate the selective pressure derived by adaptive evolution. To assess whether or not positive selection occurred, we employed several different possibility tests (M1 vs. M2 and M7 vs. M8). The proportions of positively selected sites were significant, with a likelihood log value of 93.63 for the BPIFA1 protein. The Selecton server was used on the same dataset to reconfirm positive selection for specific sites by employing the Mechanistic-Empirical Combination model, thus providing additional evidence supporting the findings of positive selection. There was convincing evidence for positive selection signals in the BPIFA1 genes of mammalian species, which was more significant for selection signs and creating signals. We performed probability tests comparing various models based on dN/dS ratios to recognize specific codons under positive selection pressure. We identified positively selected sites in the LBP-BPI domain of BPIFA1 proteins in the mammalian genome, including a lipid-binding domain with a very high degree of selectivity for DPPC. BPIFA1 activates the upper airway's innate immune system in response to numerous genetic signals in the mammalian genome. These findings highlight evolutionary advancements in immunoregulatory effects that play a significant role in the antibacterial and antiviral defenses of mammalian species.