Detecting individual sites subject to episodic diversifying selection

Complete coding region sequence analyses and antigenic characterization of emerging lineage G-IX of foot- and-mouth disease virus serotype Asia1

Foot-and-mouth disease Virus (FMDV) serotype Asia1 is prevalent in the Indian subcontinent, with only G-III and G-VIII reported in India until 2020. However, in 2019, a novel genetic group within serotype Asia1, designated as G-IX, emerged in Bangladesh, followed by its detection in India in 2020. This report presents analyses of the complete coding region sequences of the G-IX lineage isolates. The length of the open reading frame (ORF) of the two G-IX isolates was 6990 nucleotides without any deletion or insertion. The G-IX isolates showed the highest sequence similarity with an isolate of G-III at the ORF, L, P2, and P3 regions, and with an isolate of G-VIII at the P1 region. Phylogenetic analysis based on the capsid region (P1) supports the hypothesis that G-VIII and G-IX originated from a common ancestor, as speculated earlier. Further, VP1 region-based phylogenetic analyses revealed the re-emergence of G-VIII after a gap of 3 years. One isolate of G-VIII collected during 2023 revealed a codon insertion in the G-H loop of VP1. The vaccine matching studies support the suitability of the currently used Indian vaccine strain IND63/1972 to contain outbreaks due to viruses belonging to G-IX.

Genetic analysis of the circumsporozoite gene in Plasmodium falciparum isolates from Cameroon: Implications for efficacy and deployment of RTS,S/AS01 vaccine

Current understanding of genetic polymorphisms and natural selection in Plasmodium falciparum circumsporozoite (PfCSP), the leading malaria vaccine, is crucial for the development of next-generation vaccines, and such data is lacking in Africa. Blood samples were collected among Plasmodium-infected individuals living in four Cameroonian areas (Douala, Maroua, Mayo-Oulo, Pette). DNA samples were amplified using nested PCR protocols, sequenced, and BLASTed. Single nucleotide polymorphisms (SNPs) were analysed in each PfCSP region, and their impact on PfCSP function/structure was predicted in silico. The N-terminal region showed a limited polymorphism with four haplotypes, and three novel SNPs (N68Y, R87W, K93E) were found. Thirty-five haplotypes were identified in the central region, with several variants (e.g., NVNP and KANP). The C-terminal region was also highly diverse, with 25 haplotypes and eight novel SNPs (N290D, N308I, S312G, K317A, V344I, D356E, E357L, D359Y). Most polymorphic codon sites were mainly observed in the Th2R subregion in isolates from Douala and Pette. The codon site 321 was under episodic positive selection. One novel (E357L) and three known (K322I, G349D, D359Y) SNPs show an impact on function/structure. This study showed extensive genetic diversity with geographical patterns and evidence of the selection of Cameroonian PfCSP central and C-terminal regions.

Investigating role of positively selected genes and mutation sites of ERG11 in drug resistance of Candida albicans

The steep increase in acquired drug resistance in Candida isolates has posed a great challenge in the clinical management of candidiasis globally. Information of genes and codon sites that are positively selected during evolution can provide insights into the mechanisms driving antifungal resistance in Candida. This study aimed to create a manually curated list of genes of Candida spp. reported to be associated with antifungal resistance in literature, and further investigate the structure-function implications of positively selected genes and mutation sites. Sequence analysis of antifungal drug resistance associated gene sequences from various species and strains of Candida revealed that ERG11 and MRR1 of C. albicans were positively selected during evolution. Four sites in ERG11 and two sites in MRR1 of C. albicans were positively selected and associated with drug resistance. These four sites (132, 405, 450, and 464) of ERG11 are predictive markers for azole resistance and have evolved over time. A well-characterized crystal structure of sterol-14-α-demethylase (CYP51) encoded by ERG11 is available in PDB. Therefore, the stability of CYP51 in complex with fluconazole was evaluated using MD simulations and molecular docking studies for two mutations (Y132F and Y132H) reported to be associated with azole resistance in literature. These mutations induced high flexibility in functional motifs of CYP51. It was also observed that residues such as I304, G308, and I379 of CYP51 play a critical role in fluconazole binding affinity. The insights gained from this study can further guide drug design strategies addressing antimicrobial resistance.

Comparative selective pressure analysis on mitochondrial protein-coding genes in flying squirrels (Pteromyini) and tree squirrels (Sciurini)

Different animal groups with varying locomotion modes may have unique energy requirements. Mitochondria produce adenosine triphosphate (ATP) and reactive oxygen species via oxidative phosphorylation to support organisms energy requirements. The tribes Pteromyini (flying squirrels) and Sciurini (tree squirrels), two closely related taxa within the family Sciuridae, exhibit distinct locomotion modes, energy requirements, and likely face different selective pressures on mitochondrial protein-coding genes (PCGs). We analysed 13 mitochondrial genome sequences from species belonging to the tribe Pteromyini and 117 from species belonging to the tribe Sciurini. Phylogenetic analysis revealed Pteromyini and Sciurini formed a sister relationship within the family Sciuridae. Among the 13 PCGs, ATP8 exhibited the highest dN/dS values, while COX1 showed the lowest. The background selection ratio (ω2) values for six genes (ND1, ND2, ND4, ATP6, ND5, and COX3) in Pteromyini were lower than the foreground selection ratio (ω0) values observed in Sciurini. A RELAX analysis revealed that CYTB, ND4, ATP6, and COX3 genes experienced intensified in selection strength. BUSTED analysis identified stronger signatures of diversifying selection in CYTB and ATP6, highlighting amino acid changes. MEME identified episodic diversifying selection at specific sites among eight PCGs. These findings revealed distinct selective pressures on PCGs in flying and tree squirrels.

Positive selection, genetic recombination, and intra-host evolution in novel equine coronavirus genomes and other members of the Embecovirus subgenus

There are several examples of coronaviruses in the Betacoronavirus subgenus Embecovirus that have jumped from an animal to the human host. Studying how evolutionary factors shape coronaviruses in non-human hosts may provide insight into the coronavirus host-switching potential. Equids, such as horses and donkeys, are susceptible to equine coronaviruses (ECoVs). With increased testing prevalence, several ECoV genome sequences have become available for molecular evolutionary analyses, especially those from the United States of America (USA). To date, no analyses have been performed to characterize evolution within coding regions of the ECoV genome. Here, we obtain and describe four new ECoV genome sequences from infected equines from across the USA presenting clinical symptoms of ECoV, and infer ECoV-specific and Embecovirus-wide patterns of molecular evolution. Within two of the four data sets analyzed, we find evidence of intra-host evolution within the nucleocapsid (N) gene, suggestive of quasispecies development. We also identify 12 putative genetic recombination events within the ECoV genome, 11 of which fall in ORF1ab. Finally, we infer and compare sites subject to positive selection on the ancestral branch of each major Embecovirus member clade. Specifically, for the two currently identified human coronavirus (HCoV) embecoviruses that have spilled from animals to humans (HCoV-OC43 and HCoV-HKU1), we find that there are 42 and 2 such sites, respectively, perhaps reflective of the more complex ancestral evolutionary history of HCoV-OC43, which involves several different animal hosts.IMPORTANCEThe Betacoronavirus subgenus Embecovirus contains coronaviruses that not only pose a health threat to animals and humans, but also have jumped from animal to human host. Equids, such as horses and donkeys are susceptible to equine coronavirus (ECoV) infections. No studies have systematically examined evolutionary patterns within ECoV genomes. Our study addresses this gap and provides insight into intra-host ECoV evolution from infected horses. Further, we identify and report natural selection pattern differences between two embecoviruses that have jumped from animals to humans [human coronavirus OC43 and HKU1 (HCoV-OC43 and HCoV-HKU1, respectively)], and hypothesize that the differences observed may be due to the different animal host(s) that each virus circulated in prior to its jump into humans. Finally, we contribute four novel, high-quality ECoV genomes to the scientific community.

Functional and regulatory diversification of Period genes responsible for circadian rhythm in vertebrates

The Period genes (Per) play essential roles in modulating the molecular circadian clock timing in a broad range of species, which regulates the physiological and cellular rhythms through the transcription-translation feedback loop. While the Period gene paralogs are widely observed among vertebrates, the evolutionary history and the functional diversification of Per genes across vertebrates are not well known. In this study, we comprehensively investigated the evolution of Per genes at the copy number and sequence levels, including de novo binding motif discovery by comparative genomics. We also determined the lineage-specific transcriptome landscape across tissues and developmental stages and phenotypic effects in public RNA-seq data sets of model species. We observed multiple lineage-specific gain and loss events Per genes, though no simple association was observed between ecological factors and Per gene numbers in each species. Among salmonid fish species, the per3 gene has been lost in the majority, whereas those retaining the per3 gene exhibit not a signature of relaxed selective constraint but rather a signature of intensified selection. We also determined the signature of adaptive diversification of the CRY-binding region in Per1 and Per3, which modulates the circadian rhythm. We also discovered putative regulatory sequences, which are lineage-specific, suggesting that these cis-regulatory elements may have evolved rapidly and divergently across different lineages. Collectively, our findings revealed the evolution of Per genes and their fine-tuned contribution to the plastic and precise regulation of circadian rhythms in various vertebrate taxa.

Rapid evolution leads to extensive genetic diversification of cattle flu Influenza D virus

Influenza D virus (IDV), the cattle flu virus, is a novel multi-host RNA virus, circulating silently worldwide, with widespread seropositivity among US cattle, reaching up to 80% in some areas raising a potential threat of cattle-to-human transmission. Currently, five genetic lineages of IDV have been described, but their evolutionary dynamics have not been studied. Although IDV was first identified in 2011, our comprehensive analysis of all known IDV genomes suggests that the earliest ancestors of IDV likely to have evolved towards the end of the 20th century and D/OK lineage appears to have emerged in 2005. We confirmed a significantly higher substitution rate in IDV than in Influenza C virus, which is consistent with their global distribution and multi-host tropism. We identified multiple sub-populations within the D/OK lineage, highlighting extensive diversification and dissemination. Other findings are evidence for potential reassortment among IDV strains in the USA and transboundary circulation of IDV in Europe with introductions into Danish cattle, some of which potentially originated from France. IDV, an emerging virus with a higher rate of evolution and uncontrolled circulation, could facilitate its adaptation to humans. Our findings underscore the importance of targeted surveillance for IDV in humans and at-risk animal populations.

Squash vein yellowing virus from California emerged in the Middle East via intragenic and intergeneric recombination events in the hypervariable potyvirus P1 and ipomovirus P1a genes

We present the complete sequence of the genomic RNA of an isolate of squash vein yellowing virus (Ipomovirus cucurbitavenaflavi) from California (SqVYV-CA) and show it is a recombinant virus with a highly divergent 5' UTR and proximal P1a gene. The evolution of SqVYV-CA involved an intrageneric event between unknown potyviruses, related to isolates of papaya ringspot virus (Potyvirus papayanuli) from the Old World, and an intergeneric event between this recombinant potyvirus (minor parent) and an isolate of SqVYV from Israel (SqVYV-IL) (major parent). These events occurred in mixed infections and in the potyvirus P1 and ipomovirus P1a recombination hotspots and resulted in SqVYV-CA having a potyvirus 5' UTR and chimeric P1-P1a gene/protein and the remainder of the genome from SqVYV-IL. The SqVYV-CA chimeric P1-P1a gene is under positive selection, and the protein is intrinsically disordered and may localize to the nucleus via nuclear localization signals in the P1 part. The C-terminal SqVYV-IL P1a part also diverged but retained the conserved serine protease motif. Furthermore, substantial divergence in SqVYV isolates from the Middle East was associated with genetic drift and a long evolutionary history in this region. The finding that the host range and symptomatology in cucurbits of SqVYV-CA is similar to those of SqVYV from Florida and SqVYV-IL, indicated that the recombinant part of the genome had no obvious effect on the virus-host interaction. A divergent part of the P1 sequence of the SqVYV-CA P1-P1a gene was used to develop a primer pair and RT-PCR test for specific detection of SqVYV-CA. This test was used to detect spread of SqVYV-CA to a new production area of California in 2021 and 2022. Together, these results demonstrate (i) a high level of genetic diversity exists among isolates of SqVYV and involved intra- and intergeneric recombination and genetic drift (mutation), (ii) evidence that SqVYV originated in the Middle East and that there were independent introductions into the New World and (iii) the remarkable genetic flexibility of the 5' proximal genes of these viruses.

Adaptive evolution and functional significance of the PPARGC1A gene across diverse animal species

Codon-based analyses of the PPARGC1A gene across 38 vertebrate species were deployed to elucidate patterns of evolutionary change. Employing maximum likelihood assessments through MEGA, we scrutinized 447 codon positions addressing the entire coding region, excluding positions mired by gaps or missing data. Distinct codons manifested variance in selection pressures, particularly codons 4, 11, 66, and 123, which exhibited positive dN-dS values suggestive of positive selection. Codon 137 displayed the most pronounced dN-dS value, signifying intensified selective advantage. Meanwhile, codons 30 and 90 portrayed near-neutral scores, indicative of purifying selection. Complementary computational methods (IFEL, REL, FUBAR, and SLAC) confirmed positive selection at specific codon sites, with varying degrees of corroboration. The integration of mixed-effect modeling (MEME) identified episodic diversifying selection, pinpointing codons that underwent selection episodes in certain lineages. Refined codon model selection lent insight into substitution rates, revealing nuanced degrees of evolutionary conservation among different codons. Supporting these genetic insights, the phylogenetic analysis highlighted relationships among the PPARGC1A sequences and domain analysis confirmed conserved features across species, while protein-protein interaction networks suggested a complex web of functional interdependencies. These findings advance our understanding of the PPARGC1A gene's evolutionary trajectory and underscore the gene's potential adaptive significance within diverse vertebrate lineages.

Relaxed selection in evolution of genes regulating limb development gives clue to variation in forelimb morphology of cetaceans and other mammals

Cetaceans have evolved unique limb structures, such as flippers, due to genetic changes during their transition to aquatic life. However, the full understanding of the genetic and evolutionary mechanisms behind these changes is still developing. By examining 25 limb-related protein-coding genes across various mammalian species, we compared genetic changes between aquatic mammals, like whales, and other mammals with unique limb structures such as bats, rodents and elephants. Our findings revealed significant modifications in limb-related genes, including variations in the Hox, GDF5 and Evx genes. Notably, a relaxed selection in several key genes was observed, suggesting a lifting of developmental constraints, which might have facilitated the emergence of morphological innovations in cetacean limb morphology. We also uncovered non-synonymous changes, insertions and deletions in these genes, particularly in the polyalanine tract of HOXD13, which are distinctive to cetaceans or convergent with other aquatic mammals. These genetic variations correlated with the diverse and specialized limb structures observed in cetaceans, indicating a complex interplay of relaxed selection and specific mutations in mammalian limb evolution.

Interaction of chikungunya virus glycoproteins with macrophage factors controls virion production

Despite their role as innate sentinels, macrophages can serve as cellular reservoirs of chikungunya virus (CHIKV), a highly-pathogenic arthropod-borne alphavirus that has caused large outbreaks among human populations. Here, with the use of viral chimeras and evolutionary selection analysis, we define CHIKV glycoproteins E1 and E2 as critical for virion production in THP-1 derived human macrophages. Through proteomic analysis and functional validation, we further identify signal peptidase complex subunit 3 (SPCS3) and eukaryotic translation initiation factor 3 subunit K (eIF3k) as E1-binding host proteins with anti-CHIKV activities. We find that E1 residue V220, which has undergone positive selection, is indispensable for CHIKV production in macrophages, as its mutation attenuates E1 interaction with the host restriction factors SPCS3 and eIF3k. Finally, we show that the antiviral activity of eIF3k is translation-independent, and that CHIKV infection promotes eIF3k translocation from the nucleus to the cytoplasm, where it associates with SPCS3. These functions of CHIKV glycoproteins late in the viral life cycle provide a new example of an intracellular evolutionary arms race with host restriction factors, as well as potential targets for therapeutic intervention.

Evaluation of the Deletion of African Swine Fever Virus E111R Gene from the Georgia Isolate in Virus Replication and Virulence in Domestic Pigs

African swine fever virus (ASFV) is the causative agent of an often lethal disease in domestic pigs, African swine fever (ASF). ASF is currently a pandemic disease challenging pig production in Eurasia. While the ASFV genome encodes for over 160 proteins, the function of most of them are still not characterized. Among those ASF genes with unknown functions is the E111R gene. It has been recently reported that the deletion of the E111R gene from the genome of the virulent Chinese field isolate SY18 strain produced a reduction of virus virulence when pigs were inoculated at relatively low doses. Conversely, we report here that deletion of the ASFV gene E111R in the Georgia 2010 isolate does not alter the virulence of the parental virus in experimentally inoculated pigs. A recombinant virus lacking the E111R gene, ASFV-G-∆E111R was intramuscularly (IM) inoculated in domestic pigs at a dose of 102 HAD50 of ASFV-G-∆E111R and compared with animals that received a similar dose of virulent ASFV-G. Both, animals inoculated with either the recombinant ASFV-G-∆E111R or the parental virus developed a fatal form of the disease and were euthanized around the 6th-7th day post-inoculation (dpi).

Whole genome characteristics of hedgehog coronaviruses from Poland and analysis of the evolution of the Spike protein for its interspecies transmission potential

BACKGROUND

The hedgehogs have been recently identified as possible reservoir of Middle East respiratory syndrome coronavirus like (MERS-CoV-like). These viruses were classified as a distinct Betacoronavirus erinacei (BCoV-Eri) species within the MerBCoV-Eriirus subgenus. As coronaviruses are known for their ability to jump between different hosts, including humans, this can pose a particular threat to people in direct contact with hedgehogs, such as those working at animal asylums. Our previous studies have shown the presence of BCoV-Eri strains in animals collected in the wildlife rehabilitation centre. This study aimed to investigate the presence of CoV in subsequent hedgehogs collected from the urban area of Poland and their molecular characteristics.

RESULTS

Monitoring for the presence of coronavirus infection in hedgehogs revealed five positive individuals. The presence of BCoV-Eri was found in a total of 20% of animals tested. Our analyses revealed no correlation between CoVs positivity and animal health conditions but a higher probability of such infection in juveniles and females. The whole genome of two Polish Hedgehog coronavirus 1 strains were sequenced and compared with available counterparts from European and Asian countries. Phylogenetic analysis showed that both CoV strains formed common cluster with other similar MerBCoV-Eriirus, but they were also found to be genetically variable and most changes in the S protein were identified. Our analysis revealed that some S protein sites of the Hedgehog coronavirus 1 strains evolved under positive selection pressure and of five such sites, three are in the S1 region while the other two in the S2 region of the Spike.

CONCLUSIONS

BCoV-Eri is to some extent prevalent in wildlife asylums in Poland. Given that the S protein of BCoVs-Eri is highly variable and that some sites of this protein evolve under positive selection pressure, these strains could potentially acquire a favourable feature for cross-species transmission. Consequently, the threat to humans working in such asylums is particularly high. Adequate biosecurity safeguards, but also human awareness of such risks, are therefore essential.

Incorporation of nitrogen in antinutritional Solanum alkaloid biosynthesis

Steroidal glycoalkaloids (SGAs) are specialized metabolites produced by hundreds of Solanum species including food crops, such as tomato, potato and eggplant. Unlike true alkaloids, nitrogen is introduced at a late stage of SGA biosynthesis through an unknown transamination reaction. Here, we reveal the mechanism by which GLYCOALKALOID METABOLISM12 (GAME12) directs the biosynthesis of nitrogen-containing steroidal alkaloid aglycone in Solanum. We report that GAME12, a neofunctionalized γ-aminobutyric acid (GABA) transaminase, undergoes changes in both active site specificity and subcellular localization to switch from its renown and generic activity in core metabolism to function in a specialized metabolic pathway. Moreover, overexpression of GAME12 alone in engineered S. nigrum leaves is sufficient for de novo production of nitrogen-containing SGAs. Our results highlight how hijacking a core metabolism GABA shunt enzyme is crucial in numerous Solanum species for incorporating a nitrogen to a steroidal-specialized metabolite backbone and form defensive alkaloids.

Independent repeated mutations within the alphaviruses Ross River virus and Barmah Forest virus indicates convergent evolution and past positive selection in ancestral populations despite ongoing purifying selection

Ross River virus (RRV) and Barmah Forest virus (BFV) are arthritogenic arthropod-borne viruses (arboviruses) that exhibit generalist host associations and share distributions in Australia and Papua New Guinea (PNG). Using stochastic mapping and discrete-trait phylogenetic analyses, we profiled the independent evolution of RRV and BFV signature mutations. Analysis of 186 RRV and 88 BFV genomes demonstrated their viral evolution trajectories have involved repeated selection of mutations, particularly in the nonstructural protein 1 (nsP1) and envelope 3 (E3) genes suggesting convergent evolution. Convergent mutations in the nsP1 genes of RRV (residues 248 and 441) and BFV (residues 297 and 447) may be involved with catalytic enzyme mechanisms and host membrane interactions during viral RNA replication and capping. Convergent E3 mutations (RRV site 59 and BFV site 57) may be associated with enzymatic furin activity and cleavage of E3 from protein precursors assisting viral maturation and infectivity. Given their requirement to replicate in disparate insect and vertebrate hosts, convergent evolution in RRV and BFV may represent a dynamic link between their requirement to selectively 'fine-tune' intracellular host interactions and viral replicative enzymatic processes. Despite evidence of evolutionary convergence, selection pressure analyses did not reveal any RRV or BFV amino acid sites under strong positive selection and only weak positive selection for nonstructural protein sites. These findings may indicate that their alphavirus ancestors were subject to positive selection events which predisposed ongoing pervasive convergent evolution, and this largely supports continued purifying selection in RRV and BFV populations during their replication in mosquito and vertebrate hosts.

The classification, origin, and evolutionary dynamics of severe fever with thrombocytopenia syndrome virus circulating in East Asia

The classification of severe fever with thrombocytopenia syndrome virus (SFTSV) lacked consistency due to limited virus sequences used across previous studies, and the origin and transmission dynamics of the SFTSV remains not fully understood. In this study, we analyzed the diversity and phylodynamics of SFTSV using the most comprehensive and largest dataset publicly available for a better understanding of SFTSV classification and transmission. A total of 1267 L segments, 1289 M segments, and 1438 S segments collected from China, South Korea, and Japan were included in this study. Maximum likelihood trees were reconstructed to classify the lineages. Discrete phylogeographic analysis was conducted to infer the phylodynamics of SFTSV. We found that the L, M, and S segments were highly conserved, with mean pairwise nucleotide distances of 2.80, 3.36, and 3.35% and could be separated into 16, 13, and 15 lineages, respectively. The evolutionary rate for L, M, and the S segment was 0.61 × 10-4 (95% HPD: 0.48-0.73 × 10-4), 1.31 × 10-4 (95% HPD: 0.77-1.77 × 10-4) and 1.27 × 10-4 (95% HPD: 0.65-1.85 × 10-4) subs/site/year. The SFTSV most likely originated from South Korea around the year of 1617.6 (95% HPD: 1513.1-1724.3), 1700.4 (95% HPD: 1493.7-1814.0), and 1790.1 (95% HPD: 1605.4-1887.2) for L, M, and S segments, respectively. Hubei Province in China played a critical role in the geographical expansion of the SFTSV. The effective population size of SFTSV peaked around 2010 to 2013. We also identified several codons under positive selection in the RdRp, Gn-Gc, and NS genes. By leveraging the largest dataset of SFTSV, our analysis could provide new insights into the evolution and dispersal of SFTSV, which may be beneficial for the control and prevention of severe fever with thrombocytopenia syndrome.

Evolution of Key Oxygen-Sensing Genes Is Associated with Hypoxia Tolerance in Fishes

Low dissolved oxygen (hypoxia) is recognized as a major threat to aquatic ecosystems worldwide. Because oxygen is paramount for the energy metabolism of animals, understanding the functional and genetic drivers of whole-animal hypoxia tolerance is critical to predicting the impacts of aquatic hypoxia. In this study, we investigate the molecular evolution of key genes involved in the detection of and response to hypoxia in ray-finned fishes: the prolyl hydroxylase domain (PHD)-hypoxia-inducible factor (HIF) oxygen-sensing system, also known as the EGLN (egg-laying nine)-HIF oxygen-sensing system. We searched fish genomes for HIFA and EGLN genes, discovered new paralogs from both gene families, and analyzed protein-coding sites under positive selection. The physicochemical properties of these positively selected amino acid sites were summarized using linear discriminants for each gene. We employed phylogenetic generalized least squares to assess the relationship between these linear discriminants for each HIFA and EGLN and hypoxia tolerance as reflected by the critical oxygen tension (Pcrit) of the corresponding species. Our results demonstrate that Pcrit in ray-finned fishes correlates with the physicochemical variation of positively selected sites in specific HIFA and EGLN genes. For HIF2A, two linear discriminants captured more than 90% of the physicochemical variation of these sites and explained between 20% and 39% of the variation in Pcrit. Thus, variation in HIF2A among fishes may contribute to their capacity to cope with aquatic hypoxia, similar to its proposed role in conferring tolerance to high-altitude hypoxia in certain lineages of terrestrial vertebrates.

Positive Selection Drives the Evolution of the Structural Maintenance of Chromosomes (SMC) Complexes

Structural Maintenance of Chromosomes (SMC) complexes are an evolutionary conserved protein family. In most eukaryotes, three SMC complexes have been characterized, as follows: cohesin, condensin, and SMC5/6 complexes. These complexes are involved in a plethora of functions, and defects in SMC genes can lead to an increased risk of chromosomal abnormalities, infertility, and cancer. To investigate the evolution of SMC complex genes in mammals, we analyzed their selective patterns in an extended phylogeny. Signals of positive selection were identified for condensin NCAPG, for two SMC5/6 complex genes (SMC5 and NSMCE4A), and for all cohesin genes with almost exclusive meiotic expression (RAD21L1, REC8, SMC1B, and STAG3). For the latter, evolutionary rates correlate with expression during female meiosis, and most positively selected sites fall in intrinsically disordered regions (IDRs). Our results support growing evidence that IDRs are fast evolving, and that they most likely contribute to adaptation through modulation of phase separation. We suggest that the natural selection signals identified in SMC complexes may be the result of different selective pressures: a host-pathogen arms race in the condensin and SMC5/6 complexes, and an intragenomic conflict for meiotic cohesin genes that is similar to that described for centromeres and telomeres.

Molecular evolution of Hokkaido virus, a genotype of Orthohantavirus puumalaense, among Myodes rodents

Viruses in the genus Orthohantavirus within the family Hantaviridae cause human hantavirus infections and represent a threat to public health. Hokkaido virus (HOKV), a genotype of Orthohantavirus puumalaense (Puumala virus; PUUV), was first identified in Tobetsu, Hokkaido, Japan. Although it is genetically related to the prototype of PUUV, the evolutionary pathway of HOKV is unclear. We conducted a field survey in a forest in Tobetsu in 2022 and captured 44 rodents. Complete coding genome sequences of HOKVs were obtained from five viral-RNA-positive rodents (four Myodes rufocanus bedfordiae and one Apodemus speciosus). Phylogenetic analysis revealed a close relationship between the phylogenies and geographical origins of M. rufocanus-related orthohantaviruses. Comparison of the phylogenetic trees of the S segments of orthohantaviruses and the cytochrome b genes of Myodes species suggested that Myodes-related orthohantaviruses evolved in Myodes rodent species as a result of genetic isolation and host switching.

Detection and Molecular Characterization of Porcine Teschoviruses in India: Identification of New Genotypes

Porcine Teschoviruses (PTVs) are ubiquitous enteric viral pathogens that infect pigs and wild boars worldwide. PTVs have been responsible for causing the severe clinical disease (Teschen disease) to asymptomatic infections. However, to date, limited information is available on large-scale epidemiological data and molecular characterization of PTVs in several countries. In this study, we report epidemiological data on PTVs based on screening of 534 porcine fecal samples from different states of India and a RT-PCR based detection of PTVs shows a percent positivity of 8.24% (44/534). The PTV prevalence varied among different regions of the country with the highest detection rates observed in the state of Karnataka (38.1%). Phylogenetic analysis based on VP1 gene reveals the presence of PTV genotype 6 and 13 along with some unassigned novel genotypes which did not cluster with any of the established PTV genotypes (PTV 1-PTV 13). Indian PTV 6 strains are genetically closest to the Spanish strains (85.7-94.4%) whereas PTV 13 and novel genotype strains were found to be more similar to the Chinese strains (88.1-99.1%). Using recombination detection software, no Indian PTVs found to be recombinant on VP1 gene and selection pressure analysis revealed the purifying selection in the several sites of the VP1 gene of PTVs. The Bayesian analysis of Indian PTVs shows 1.16 × 10-4 substitution/site/year as the mean evolutionary rate. Further, isolation of the novel PTV strains from India and more detailed investigation much needed to know the evolutionary history of PTV strains circulating in porcine populations in India.

Characterizing human respiratory syncytial virus among children admitted with acute respiratory tract infections from 2019 to 2022

Respiratory syncytial virus is a major causative agent of lower respiratory tract infection in children, especially infants with substantial morbidity and mortality implications. The virus undergoes continuous evolution documented by accumulation of mutations in the glycoprotein gene necessitating vigilant surveillance to provide essential data to epidemiologists and researchers involved in development of vaccines. This study was aimed to perform molecular characterization of respiratory syncytial virus (RSV) among children ≤ 5 years admitted in hospital. In the current study we observed RSV-A (2019 (n = 95) and 2021 (n = 61) seasons) and RSV-B (2022 season (n = 68)). Phylogenetic analysis revealed all RSV-A strains (n = 47) to be GA.2.3.5 and RSV-B (n = 22) were classified as GB.5.0.5a. Selection pressure analysis identified one positive (P274L/V) and one negative site (P230T) in RSV-A, while in RSV-B there was only one negatively selected site (S295). This study spanning over three seasons contributes to RSV evolutionary dynamics in India emphasizing the importance of on-going surveillance to inform effective public health strategies and vaccine development efforts.

Genetic substructure and host-specific natural selection trend across vaccine-candidate ORF-2 capsid protein of hepatitis-E virus

Hepatitis E virus is a primary cause of acute hepatitis worldwide. The present study attempts to assess the genetic variability and evolutionary divergence among HEV genotypes. A vaccine promising capsid-protein coding ORF-2 gene sequences of HEV was evaluated using phylogenetics, model-based population genetic methods and principal component analysis. The analyses unveiled nine distinct clusters as subpopulations for six HEV genotypes. HEV-3 genotype samples stratified into four different subgroups, while HEV-4 stratified into three additional subclusters. Rabbit-infectious HEV-3ra samples constitute a distinct cluster. Pairwise analysis identified marked genetic distinction of HEV-4c and HEV-4i subgenotypes compared to other genotypes. Numerous admixed, inter and intragenotype recombinant strains were detected. The MEME method identified several ORF-2 codon sites under positive selection. Some selection signatures lead to amino acid substitutions within ORF-2, resulting in altered physicochemical features. Moreover, a pattern of host-specific adaptive signatures was identified among HEV genotypes. The analyses conclusively depict that recombination and episodic positive selection events have shaped the observed genetic diversity among different HEV genotypes. The significant genetic diversity and stratification of HEV-3 and HEV-4 genotypes into subgroups, as identified in the current study, are noteworthy and may have implications for the efficacy of anti-HEV vaccines.

Differences between phytophagous and predatory species in Pentatomidae based on the mitochondrial genome

Pentatomidae includes many species of significant economic value as plant pests and biological control agents. The feeding habits of Pentatomidae are closely related to their energy metabolism and ecological adaptations. In this study, we sequenced the mitochondrial genomes of 12 Asopinae species using the next-generation sequencing to explore the effect of dietary changes on mitochondrial genome evolution. Notably, all sequences were double-stranded circular DNA molecules containing 37 genes and one control region. We then compared and analyzed the mitochondrial genome characteristics of phytophagous and predatory bugs. Notably, no significant difference was observed in the length of the mitochondrial genomes between the predatory and phytophagous bugs. However, the AT content was higher in the mitochondrial genomes of phytophagous bugs than that of predatory bugs. Moreover, phytophagous bugs prefer codon usage patterns ending in A/T compared with predatory bugs. The evolution rate of predatory bugs was lower than that of phytophagous bugs. The phylogenetic relationships across phytophagous bugs' lineages were largely consistent at depth nodes based on different datasets and tree-reconstructing methods, and strongly supported the monophyly of predatory bugs. Additionally, the estimated divergence times indicated that Pentatomidae explosively radiated in the Early Cretaceous. Subsequently, the subfamily Asopinae and the genus Menida diverged in the Late Cretaceous. Our research results provide data supporting for the evolutionary patterns and classification of Pentatomidae.

Chromosome-level genome assembly of a deep-sea Venus flytrap sea anemone sheds light upon adaptations to an extremely oligotrophic environment

The Venus flytrap sea anemone Actinoscyphia liui inhabits the nutrient-limited deep ocean in the tropical western Pacific. Compared with most other sea anemones, it has undergone a distinct modification of body shape similar to that of the botanic flytrap. However, the molecular mechanism by which such a peculiar sea anemone adapts to a deep-sea oligotrophic environment is unknown. Here, we report the chromosomal-level genome of A. liui constructed from PacBio and Hi-C data. The assembled genome is 522 Mb in size and exhibits a continuous scaffold N50 of 58.4 Mb. Different from most other sea anemones, which typically possess 14-18 chromosomes per haplotype, A. liui has only 11. The reduced number of chromosomes is associated with chromosome fusion, which likely represents an adaptive strategy to economize energy in oligotrophic deep-sea environments. Comparative analysis with other deep-sea sea anemones revealed adaptive evolution in genes related to cellular autophagy (TMBIM6, SESN1, SCOCB and RPTOR) and mitochondrial energy metabolism (MDH1B and KAD2), which may aid in A. liui coping with severe food scarcity. Meanwhile, the genome has undergone at least two rounds of expansion in gene families associated with fast synaptic transmission, facilitating rapid responses to water currents and prey. Positive selection was detected on putative phosphorylation sites of muscle contraction-related proteins, possibly further improving feeding efficiency. Overall, the present study provides insights into the molecular adaptation to deep-sea oligotrophic environments and sheds light upon potential effects of a novel morphology on the evolution of Cnidaria.

Viroporin-like activity of the hairpin transmembrane domain of African swine fever virus B169L protein

African swine fever virus (ASFV) is the causative agent of a contagious disease affecting wild and domestic swine. The function of B169L protein, as a potential integral structural membrane protein, remains to be experimentally characterized. Using state-of-the-art bioinformatics tools, we confirm here earlier predictions indicating the presence of an integral membrane helical hairpin, and further suggest anchoring of this protein to the ER membrane, with both terminal ends facing the lumen of the organelle. Our evolutionary analysis confirmed the importance of purifying selection in the preservation of the identified domains during the evolution of B169L in nature. Also, we address the possible function of this hairpin transmembrane domain (HTMD) as a class IIA viroporin. Expression of GFP fusion proteins in the absence of a signal peptide supported B169L insertion into the ER as a Type III membrane protein and the formation of oligomers therein. Overlapping peptides that spanned the B169L HTMD were reconstituted into ER-like membranes and the adopted structures analyzed by infrared spectroscopy. Consistent with the predictions, B169L transmembrane sequences adopted α-helical conformations in lipid bilayers. Moreover, single vesicle permeability assays demonstrated the assembly of lytic pores in ER-like membranes by B169L transmembrane helices, a capacity confirmed by ion-channel activity measurements in planar bilayers. Emphasizing the relevance of these observations, pore-forming activities were not observed in the case of transmembrane helices derived from EP84R, another ASFV protein predicted to anchor to membranes through a α-helical HTMD. Overall, our results support predictions of viroporin-like function for the B169L HTMD.IMPORTANCEAfrican swine fever (ASF), a devastating disease affecting domestic swine, is widely spread in Eurasia, producing significant economic problems in the pork industry. Approaches to prevent/cure the disease are mainly restricted to the limited information concerning the role of most of the genes encoded by the large (160-170 kba) virus genome. In this report, we present the experimental data on the functional characterization of the African swine fever virus (ASFV) gene B169L. Data presented here indicates that the B169L gene encodes for an essential membrane-associated protein with a viroporin function.

Enhanced Diversifying Selection on Polymerase Genes in H5N1 Clade 2.3.4.4b: A Key Driver of Altered Species Tropism and Host Range Expansion

Highly pathogenic avian influenza H5N1 clade 2.3.4.4b viruses have shown unprecedented host range and pathogenicity, including infections in cattle, previously not susceptible to H5N1. We investigated whether selection pressures on clade 2.3.4.4b viral genes could shed light on their unique epidemiological features. Our analysis revealed that while the gene products of clade 2.3.4.4b H5N1 primarily undergo purifying selection, there are notable instances of episodic diversifying selection. Specifically, the polymerase genes PB2, PB1, and PA exhibit significantly greater selection pressures in clade 2.3.4.4b than all earlier H5N1 virus clades. Polymerases play critical roles in influenza virus adaptation, including viral fitness, interspecies transmission, and virulence. Our findings provide evidence that significant selection pressures have shaped the evolution of the H5N1 clade 2.3.4.4b viruses, facilitating their expanded host tropism and the potential for further adaptation to mammalian hosts. We discuss how exogenous factors, such as altered bird migration patterns and increased host susceptibility, may have contributed to the expanded host range. As H5N1 viruses continue to infect new hosts, there is a greater risk of emergent novel variants with increased pathogenicity in humans and animals. Thus, comprehensive One Health surveillance is critical to monitor transmission among avian and mammalian hosts.

Evolutionary and functional analyses reveal a role for the RHIM in tuning RIPK3 activity across vertebrates

Receptor interacting protein kinases (RIPK) RIPK1 and RIPK3 play important roles in diverse innate immune pathways. Despite this, some RIPK1/3-associated proteins are absent in specific vertebrate lineages, suggesting that some RIPK1/3 functions are conserved while others are more evolutionarily labile. Here, we perform comparative evolutionary analyses of RIPK1-5 and associated proteins in vertebrates to identify lineage-specific rapid evolution of RIPK3 and RIPK1 and recurrent loss of RIPK3-associated proteins. Despite this, diverse vertebrate RIPK3 proteins are able to activate NF-κB and cell death in human cells. Additional analyses revealed a striking conservation of the RIP homotypic interaction motif (RHIM) in RIPK3, as well as other human RHIM-containing proteins. Interestingly, diversity in the RIPK3 RHIM can tune activation of NF-κB while retaining the ability to activate cell death. Altogether, these data suggest that NF-κB activation is a core, conserved function of RIPK3, and the RHIM can tailor RIPK3 function to specific needs within and between species.

Molecular Evolutionary Analyses of the Fusion Genes in Human Parainfluenza Virus Type 4

The human parainfluenza virus type 4 (HPIV4) can be classified into two distinct subtypes, 4a and 4b. The full lengths of the fusion gene (F gene) of 48 HPIV4 strains collected during the period of 1966-2022 were analyzed. Based on these gene sequences, the time-scaled evolutionary tree was constructed using Bayesian Markov chain Monte Carlo methods. A phylogenetic tree showed that the first division of the two subtypes occurred around 1823, and the most recent common ancestors of each type, 4a and 4b, existed until about 1940 and 1939, respectively. Although the mean genetic distances of all strains were relatively wide, the distances in each subtype were not wide, indicating that this gene was conserved in each subtype. The evolutionary rates of the genes were relatively low (4.41 × 10-4 substitutions/site/year). Moreover, conformational B-cell epitopes were predicted in the apex of the trimer fusion protein. These results suggest that HPIV4 subtypes diverged 200 years ago and the progenies further diverged and evolved.

Molecular Footprints of Potato Virus Y Isolate Infecting Potatoes (Solanum tuberosum) in Kenya

Potato virus Y (PVY) is a highly diverse and genetically variable virus with various strains. Differential evolutionary routes have been reported in the genus Potyvirus, caused by natural selection pressure, mutation, and recombination, with their virulence being dependent on different environmental conditions. Despite its significance and economic impact on Solanaceous species, the understanding of PVY's phylogeography in Kenya remains limited and inadequately documented. The study centers on the molecular characterization of a Kenyan PVY isolate, GenBank accession number PP069009. In-depth phylogenetic analysis unveiled a strong evolutionary association between the Kenyan isolate and isolate [JQ924287] from the United States of America, supported by a robust 92% probability. Recombinant analyses exposed a mosaic-like genetic architecture within the Kenyan isolate, indicating multiple gene recombination events. Selection pressure scrutiny identified specific sites under selective pressure, with evidence of positive/diversifying and negative/purifying selection. Population genetics analysis revealed a calculated nucleotide diversity (π) of 0.00354881, while analysis of molecular variance (AMOVA) unveiled a structured genetic landscape with an øST value of 0.45224. The extensive haplotype network depicted the possibility of diverse PVY strains occurring across continents. This analysis provides valuable insights into the genetic diversity and distribution of PVY globally, highlighting the importance of understanding evolutionary dynamics for effective management and control strategies of PVY on a global scale.

Evolution of Virus-like Features and Intrinsically Disordered Regions in Retrotransposon-derived Mammalian Genes

Several mammalian genes have originated from the domestication of retrotransposons, selfish mobile elements related to retroviruses. Some of the proteins encoded by these genes have maintained virus-like features; including self-processing, capsid structure formation, and the generation of different isoforms through -1 programmed ribosomal frameshifting. Using quantitative approaches in molecular evolution and biophysical analyses, we studied 28 retrotransposon-derived genes, with a focus on the evolution of virus-like features. By analyzing the rate of synonymous substitutions, we show that the -1 programmed ribosomal frameshifting mechanism in three of these genes (PEG10, PNMA3, and PNMA5) is conserved across mammals and originates alternative proteins. These genes were targets of positive selection in primates, and one of the positively selected sites affects a B-cell epitope on the spike domain of the PNMA5 capsid, a finding reminiscent of observations in infectious viruses. More generally, we found that retrotransposon-derived proteins vary in their intrinsically disordered region content and this is directly associated with their evolutionary rates. Most positively selected sites in these proteins are located in intrinsically disordered regions and some of them impact protein posttranslational modifications, such as autocleavage and phosphorylation. Detailed analyses of the biophysical properties of intrinsically disordered regions showed that positive selection preferentially targeted regions with lower conformational entropy. Furthermore, positive selection introduces variation in binary sequence patterns across orthologues, as well as in chain compaction. Our results shed light on the evolutionary trajectories of a unique class of mammalian genes and suggest a novel approach to study how intrinsically disordered region biophysical characteristics are affected by evolution.

Diversifying selection identified in immune epitopes of bovine coronavirus isolates from Irish cattle

Bovine betacoronavirus (BoCoV) is a pneumoenteric pathogen of cattle that is closely related to human coronavirus OC43. Vaccines are administered to protect against diseases caused by BoCoV, but knowledge gaps exist with regard to correlates of protection and the effect of immune evasion on driving evolution. In this study, immune epitopes were mapped onto BoCoV structural proteins, including spike and haemagglutinin esterase (HE), and then supported with targeted gene sequencing of Irish clinical isolates and selective pressure analysis. Increased prevalence of diversifying selection and amino acid changes in some mapped immune epitopes suggests that immune escape is selecting for non-synonymous mutations arising in these regions. Selection analysis and sequencing provided increased support for neutralising antibody (nAb) epitopes compared to others, suggesting that nAbs are an important arm of the immune response to BoCoV. Phylogenetic analysis of spike and HE sequences showed that Irish isolates from this study were in the European clade, except for one HE sequence that sat in the Asian/American clade, while the spike gene of this sample was in the European clade. Recombination between a European and an Asian/American isolate would give rise to such a sequence. This study has gathered evidence suggesting that pressure to evade the nAb response is contributing to BoCoV evolution.

People living with HIV co-infected with the Kaposi Sarcoma-associated Herpes Virus have a distinct HIV Tat profile and higher rates of antiretroviral virologic failure, more evident among those with Kaposi's sarcoma

Kaposi sarcoma (KS) is a neoplasm of vascular origin that promotes angiogenesis and the growth of endothelial cells triggered by the Kaposi Sarcoma-associated Herpes Virus (KSHV). When associated with HIV, KSHV becomes more aggressive and rapidly evolves. The HIV-1 TAT protein can be essential in developing AIDS-associated KS by promoting angiogenesis and increasing KSHV replication. Therefore, we evaluated the genetic profile of the first exon of tat gene among groups of people living with HIV (PLHIV) with (case group, n = 36) or without KS, this later with (positive control group, n = 46) and without KSHV infection (negative control group, n = 24); all individuals under antiretroviral therapy. The genetic diversity, the DN/DS ratio, and the genetic entropy of the first exon of tat were higher in the case group, followed by the positive control group, which was higher than the negative control group. The number of tat codons under positive selection was seven in the case group, six in the positive control group, and one in the negative control group. The prevalence of HIV viral loads below the detection limit was equal in the case and positive control groups, which were lower than in the negative control group. The mean CD4+ T cell counts were higher in the negative control group, followed by the positive control group, and followed by the case group. These results emphasize the negative influence of KSHV in antiretroviral treatment, as well as the HIV-specific TAT profile among PLHIV who developed KS.

Genetic structure of apical membrane antigen-1 in Plasmodium falciparum isolates from Pakistan

Plasmodium falciparum apical membrane antigen-1 (PfAMA-1) is a major candidate for the blood-stage malaria vaccine. Genetic polymorphisms of global pfama-1suggest that the genetic diversity of the gene can disturb effective vaccine development targeting this antigen. This study was conducted to explore the genetic diversity and gene structure of pfama-1 among P. falciparum isolates collected in the Khyber Pakhtunkhwa (KP) province of Pakistan. A total of 19 full-length pfama-1 sequences were obtained from KP-Pakistan P. falciparum isolates, and genetic polymorphism and natural selection were investigated. KP-Pakistan pfama-1 exhibited genetic diversity, wherein 58 amino acid changes were identified, most of which were located in ectodomains, and domains I, II, and III. The amino acid changes commonly found in the ectodomain of global pfama-1 were also detected in KP-Pakistan pfama-1. Interestingly, 13 novel amino acid changes not reported in the global population were identified in KP-Pakistan pfama-1. KP-Pakistan pfama-1 shared similar levels of genetic diversity with global pfama-1. Evidence of natural selection and recombination events were also detected in KP-Pakistan pfama-1.

Evidence for gene flow and trait reversal during radiation of Mexican Goodeid fish

Understanding the phylogeographic history of a group and identifying the factors contributing to speciation is an important challenge in evolutionary biology. The Goodeinae are a group of live-bearing fishes endemic to Mexico. Here, we develop genomic resources for species within the Goodeinae and use phylogenomic approaches to characterise their evolutionary history. We sequenced, assembled and annotated the genomes of four Goodeinae species, including Ataeniobius toweri, the only matrotrophic live-bearing fish without a trophotaenia in the group. We estimated timings of species divergence and examined the extent and timing of introgression between the species to assess if this may have occurred during an early radiation, or in more recent episodes of secondary contact. We used branch-site models to detect genome-wide positive selection across Goodeinae, and we specifically asked whether this differs in A. toweri, where loss of placental viviparity has recently occurred. We found evidence of gene flow between geographically isolated species, suggesting vicariant speciation was supplemented by limited post-speciation gene flow, and gene flow may explain previous uncertainties about Goodeid phylogeny. Genes under positive selection in the group are likely to be associated with the switch to live-bearing. Overall, our studies suggest that both volcanism-driven vicariance and changes in reproductive mode influenced radiation in the Goodeinae.

Positive selection analyses identify a single WWE domain residue that shapes ZAP into a more potent restriction factor against alphaviruses

The host interferon pathway upregulates intrinsic restriction factors in response to viral infection. Many of them block a diverse range of viruses, suggesting that their antiviral functions might have been shaped by multiple viral families during evolution. Host-virus conflicts have led to the rapid adaptation of host and viral proteins at their interaction hotspots. Hence, we can use evolutionary genetic analyses to elucidate antiviral mechanisms and domain functions of restriction factors. Zinc finger antiviral protein (ZAP) is a restriction factor against RNA viruses such as alphaviruses, in addition to other RNA, retro-, and DNA viruses, yet its precise antiviral mechanism is not fully characterized. Previously, an analysis of 13 primate ZAP orthologs identified three positively selected residues in the poly(ADP-ribose) polymerase-like domain. However, selective pressure from ancient alphaviruses and others likely drove ZAP adaptation in a wider representation of mammals. We performed positive selection analyses in 261 mammalian ZAP using more robust methods with complementary strengths and identified seven positively selected sites in all domains of the protein. We generated ZAP inducible cell lines in which the positively selected residues of ZAP are mutated and tested their effects on alphavirus replication and known ZAP activities. Interestingly, the mutant in the second WWE domain of ZAP (N658A) is dramatically better than wild-type ZAP at blocking replication of Sindbis virus and other ZAP-sensitive alphaviruses due to enhanced viral translation inhibition. The N658A mutant is adjacent to the previously reported poly(ADP-ribose) (PAR) binding pocket, but surprisingly has reduced binding to PAR. In summary, the second WWE domain is critical for engineering a more potent ZAP and fluctuations in PAR binding modulate ZAP antiviral activity. Our study has the potential to unravel the role of ADP-ribosylation in the host innate immune defense and viral evolutionary strategies that antagonize this post-translational modification.

Investigating the emergence of a zoonotic virus: phylogenetic analysis of European bat lyssavirus 1 in the UK

European bat lyssavirus 1 (EBLV-1, Lyssavirus hamburg) is predominantly detected in serotine bats (Eptesicus serotinus) and is responsible for the majority of bat rabies cases in mainland Europe. A passive bat rabies surveillance scheme detected the virus in a serotine bat in the UK for the first time in October 2018. As of May 2024, 34 cases have been reported, 20 of which involved contact with an animal and 5 reported human contact. We investigated the emergence of EBLV-1 by undertaking comprehensive sequence analysis and Bayesian phylogenetics, based on complete virus genomes of 33 UK sequences and 108 sequences covering six countries in mainland Europe (1968-2023), including 21 French EBLV-1-positive RNA samples sequenced for this study. Sequence analysis revealed extreme similarity among UK EBLV-1 sequences (99.9%-100%), implying a single source of introduction rather than multiple independent introductions. Bayesian analysis revealed that the UK EBLV-1 sequences shared their most recent common ancestor with an EBLV-1 sequence from a serotine bat detected in Brittany, France, in 2001, with an estimated date of divergence of 1997. Within the UK sequences, the earliest divergence was estimated to occur in 2007. This study provides valuable insights into the molecular epidemiology of an emerging zoonotic pathogen and improved understanding of the risks posed to public and animal health.

Epidemiological and evolutionary analysis of canine circovirus from 1996 to 2023

BACKGROUND

Canine circovirus (CanineCV), a non-enveloped virus with a circular DNA genome, has been identified in various avian and mammalian species, including domestic and wild canids. This study aimed to comprehensively analyze the prevalence of CanineCV across diverse animal species in 11 provinces of China.

RESULTS

A total of 1,666 serum samples were collected, revealing a 5.82% prevalence of CanineCV in dogs, with the highest rates being observed in southern and eastern China. Phylogenetic analysis of 266 global CanineCV genomes sourced from the NCBI identified six distinct genotypes, elucidating the complex dynamics of their evolution. Evidence suggested a potential bat origin for CanineCV, with positive selection and high rates of evolution being observed. Recombination analysis revealed dynamic genetic exchange, highlighting the intricate nature of CanineCV evolution. Mutational analysis identified key amino acid substitutions likely to influence the virus's adaptation. Additionally, glycosylation, palmitoylation, and SUMOylation sites were predicted, shedding light on crucial functional properties of the virus.

CONCLUSIONS

This study provides a global perspective on the origin, genetic diversity, and evolutionary dynamics of CanineCV. Understanding these factors is crucial for elucidating its epidemiology and potential health risks.

A positive selection at binding site 501 in the B.1 lineage might have triggered the highly infectious sub-lineages of SARS-CoV-2

The descendants of the B lineage are the most predominant variants among the SARS-CoV-2 virus due to the incorporation of new mutations augmenting the infectivity of the virus. There is a substantial increase in the transition transversion bias, nucleotide diversity and purifying selection on the spike protein in the descendants of the B lineage of the SARS-CoV-2 virus on a temporal scale. A strong bias for C-to-U substitutions is found in the genes encoding spike protein in this lineage. The positive selection has operated on the spike gene of B lineages and its sub-lineages. The B.1 lineage has undergone positive selection on site 501 of the receptor binding domain ultimately reflected in a key substitution N501Y in its three descendant lineages namely B.1.1.7, B.1.351 and P.1. The intensity of purifying selection on the multiple sites of the spike gene has increased substantially in the sub-lineages of B.1 in a timescale. The binding site 501 on the spike protein in B lineage is found to coevolve with other amino acid sites. This study sheds light on the evolutionary trajectory of the B lineage into highly infectious descendants in the recent past under the influence of positive and purifying selection exerted by natural immunity and vaccination of the host.

Molecular characterization of canine circovirus based on the Capsid gene in Thailand

BACKGROUND

Canine circovirus (CanineCV) is a single-stranded circular DNA virus that infects domestic and wild canids in many countries. CanineCV is associated with gastroenteritis and diarrhea, respiratory disease, and generalized vasculitis leading to a fatal event. The Capsid protein (Cap) is a structural protein of the virus which has high genetic variability and plays a role in the canine immune response. In this study, we cloned the full-length CanineCV Capsid gene (Cap). In-silico analyses were used to explore the genomic and amino acid variability and natural selection acting on the Cap gene. The immune relevance for T-cell and B-cell epitopes was predicted by the immunoinformatic approach.

RESULTS

According to the Cap gene, our results showed that CanineCV was separated into five phylogenetic groups. The obtained CanineCV strain from this study was grouped with the previously discovered Thai strain (MG737385), as supported by a haplotype network. Entropy analyses revealed high nucleotide and amino acid variability of the Capsid region. Selection pressure analysis revealed four codons at positions 24, 50, 103, and 111 in the Cap protein evolved under diversifying selection. Prediction of B-cell epitopes exhibited four consensus sequences based on physiochemical properties, and eleven peptide sequences were predicted as T-cell epitopes. In addition, the positive selection sites were located within T-cell and B-cell epitopes, suggesting the role of the host immune system as a driving force in virus evolution.

CONCLUSIONS

Our study provides knowledge of CanineCV genetic diversity, virus evolution, and potential epitopes for host cell immune response.

Exploring the genetic diversity pattern of PvEBP/DBP2: A promising candidate for an effective Plasmodium vivax vaccine

Malaria remains a public health challenge. Since many control strategies have proven ineffective in eradicating this disease, new strategies are required, among which the design of a multivalent vaccine stands out. However, the effectiveness of this strategy has been hindered, among other reasons, by the genetic diversity observed in parasite antigens. In Plasmodium vivax, the Erythrocyte Binding Protein (PvEBP, also known as DBP2) is an alternate ligand to Duffy Binding Protein (DBP); given its structural resemblance to DBP, EBP/DBP2 is proposed as a promising antigen for inclusion in vaccine design. However, the extent of genetic diversity within the locus encoding this protein has not been comprehensively assessed. Thus, this study aimed to characterize the genetic diversity of the locus encoding the P. vivax EBP/DBP2 protein and to determine the evolutionary mechanisms modulating this diversity. Several intrapopulation genetic variation parameters were estimated using 36 gene sequences of PvEBP/DBP2 from Colombian P. vivax clinical isolates and 186 sequences available in databases. The study then evaluated the worldwide genetic structure and the evolutionary forces that may influence the observed patterns of genetic variation. It was found that the PvEBP/DBP2 gene exhibits one of the lowest levels of genetic diversity compared to other vaccine-candidate antigens. Four major haplotypes were shared worldwide. Analysis of the protein's 3D structure and epitope prediction identified five regions with potential antigenic properties. The results suggest that the PvEBP/DBP2 protein possesses ideal characteristics to be considered when designing a multivalent effective antimalarial vaccine against P. vivax.

The evolution of mammalian Rem2: unraveling the impact of purifying selection and coevolution on protein function, and implications for human disorders

Rad And Gem-Like GTP-Binding Protein 2 (Rem2), a member of the RGK family of Ras-like GTPases, is implicated in Huntington's disease and Long QT Syndrome and is highly expressed in the brain and endocrine cells. We examine the evolutionary history of Rem2 identified in various mammalian species, focusing on the role of purifying selection and coevolution in shaping its sequence and protein structural constraints. Our analysis of Rem2 sequences across 175 mammalian species found evidence for strong purifying selection in 70% of non-invariant codon sites which is characteristic of essential proteins that play critical roles in biological processes and is consistent with Rem2's role in the regulation of neuronal development and function. We inferred epistatic effects in 50 pairs of codon sites in Rem2, some of which are predicted to have deleterious effects on human health. Additionally, we reconstructed the ancestral evolutionary history of mammalian Rem2 using protein structure prediction of extinct and extant sequences which revealed the dynamics of how substitutions that change the gene sequence of Rem2 can impact protein structure in variable regions while maintaining core functional mechanisms. By understanding the selective pressures, protein- and gene - interactions that have shaped the sequence and structure of the Rem2 protein, we gain a stronger understanding of its biological and functional constraints.

D for dominant: porcine circovirus 2d (PCV-2d) prevalence over other genotypes in wild boars and higher viral flows from domestic pigs in Italy

Introduction

Porcine circovirus 2 (PCV-2) is a key pathogen for the swine industry at a global level. Nine genotypes, differing in epidemiology and potentially virulence, emerged over time, with PCV-2a, -2b, and -2d being the most widespread and clinically relevant. Conversely, the distribution of minor genotypes appears geographically and temporally restricted, suggesting lower virulence and different epidemiological drivers. In 2022, PCV-2e, the most genetically and phenotypically divergent genotype, was identified in multiple rural farms in North-eastern Italy. Since rural pigs often have access to outdoor environment, the introduction from wild boars was investigated.

Methods

Through a molecular and spatial approach, this study investigated the epidemiology and genetic diversity of PCV-2 in 122 wild boars across different provinces of North-eastern Italy.

Results

Molecular analysis revealed a high PCV-2 frequency (81.1%, 99/122), and classified the majority of strains as PCV-2d (96.3%, 78/81), with sporadic occurrences of PCV-2a (1.2%, 1/81) and PCV-2b (2.5%, 2/81) genotypes. A viral flow directed primarily from domestic pigs to wild boars was estimated by phylogenetic and phylodynamic analyses.

Discussion

These findings attested that the genotype replacement so far described only in the Italian domestic swine sector occurred also in wild boars. and suggested that the current heterogeneity of PCV-2d strains in Italian wild boars likely depends more on different introduction events from the domestic population rather than the presence of independent evolutionary pressures. While this might suggest PCV-2 circulation in wild boars having a marginal impact in the industrial sector, the sharing of PCV-2d strains across distinct wild populations, in absence of a consistent geographical pattern, suggests a complex interplay between domestic and wild pig populations, emphasizing the importance of improved biosecurity measures to mitigate the risk of pathogen transmission.

AOC: Analysis of Orthologous Collections - an application for the characterization of natural selection in protein-coding sequences

Motivation

Modern molecular sequence analysis increasingly relies on automated and robust software tools for interpretation, annotation, and biological insight. The Analysis of Orthologous Collections (AOC) application automates the identification of genomic sites and species/lineages influenced by natural selection in coding sequence analysis. AOC quantifies different types of selection: negative, diversifying or directional positive, or differential selection between groups of branches. We include all steps necessary to go from unaligned homologous sequences to complete results and interactive visualizations that are designed to aid in the useful interpretation and contextualization.

Results

We are motivated by a desire to make evolutionary analyses as simple as possible, and to close the disparity in the literature between genes which draw a significant amount of interest and those that are largely overlooked and underexplored. We believe that such underappreciated and understudied genetic datasets can hold rich biological information and offer substantial insights into the diverse patterns and processes of evolution, especially if domain experts are able to perform the analyses themselves.

Availability and implementation

A Snakemake [Mölder et al., 2021] application implementation is publicly available on GitHub at https://github.com/aglucaci/AnalysisOfOrthologousCollections and is accompanied by software documentation and a tutorial.

Binding Evolution of the Dengue Virus Envelope Against DC-SIGN: A Combined Approach of Phylogenetics and Molecular Dynamics Analyses Over 30 Years of Dengue Virus in Brazil

The Red Queen Hypothesis (RQH), derived from Lewis Carroll's "Through the Looking-Glass", postulates that organisms must continually adapt in response to each other to maintain relative fitness. Within the context of host-pathogen interactions, the RQH implies an evolutionary arms race, wherein viruses evolve to exploit hosts and hosts evolve to resist viral invasion. This study delves into the dynamics of the RQH in the context of virus-cell interactions, specifically focusing on virus receptors and cell receptors. We observed multiple virus-host systems and noted patterns of co-evolution. As viruses evolved receptor-binding proteins to effectively engage with cell receptors, cells countered by altering their receptor genes. This ongoing mutual adaptation cycle has influenced the molecular intricacies of receptor-ligand interactions. Our data supports the RQH as a driving force behind the diversification and specialization of both viral and host cell receptors. Understanding this co-evolutionary dance offers insights into the unpredictability of emerging viral diseases and potential therapeutic interventions. Future research is crucial to dissect the nuanced molecular changes and the broader ecological consequences of this ever-evolving battle. Here, we combine phylogenetic inferences, structural modeling, and molecular dynamics analyses to describe the epidemiological characteristics of major Brazilian DENV strains that circulated from 1990 to 2022 from a combined perspective, thus providing us with a more detailed picture on the dynamics of such interactions over time.

Dissecting positive selection events and immunological drives during the evolution of adeno-associated virus lineages

Adeno-associated virus (AAV) serotypes from primates are being developed and clinically used as vectors for human gene therapy. However, the evolutionary mechanism of AAV variants is far from being understood, except that genetic recombination plays an important role. Furthermore, little is known about the interaction between AAV and its natural hosts, human and nonhuman primates. In this study, natural AAV capsid genes were subjected to systemic evolutionary analysis with a focus on selection drives during the diversification of AAV lineages. A number of positively selected sites were identified from these AAV lineages with functional relevance implied by their localization on the AAV structures. The selection drives of the two AAV2 capsid sites were further investigated in a series of biological experiments. These observations did not support the evolution of the site 410 of the AAV2 capsid driven by selection pressure from the human CD4+ T-cell response. However, positive selection on site 548 of the AAV2 capsid was directly related to host humoral immunity because of the profound effects of mutations at this site on the immune evasion of AAV variants from human neutralizing antibodies at both the individual and population levels. Overall, this work provides a novel interpretation of the genetic diversity and evolution of AAV lineages in their natural hosts, which may contribute to their further engineering and application in human gene therapy.

Mapping the Evolutionary Space of SARS-CoV-2 Variants to Anticipate Emergence of Subvariants Resistant to COVID-19 Therapeutics

New sublineages of SARS-CoV-2 variants-of-concern (VOCs) continuously emerge with mutations in the spike glycoprotein. In most cases, the sublineage-defining mutations vary between the VOCs. It is unclear whether these differences reflect lineage-specific likelihoods for mutations at each spike position or the stochastic nature of their appearance. Here we show that SARS-CoV-2 lineages have distinct evolutionary spaces (a probabilistic definition of the sequence states that can be occupied by expanding virus subpopulations). This space can be accurately inferred from the patterns of amino acid variability at the whole-protein level. Robust networks of co-variable sites identify the highest-likelihood mutations in new VOC sublineages and predict remarkably well the emergence of subvariants with resistance mutations to COVID-19 therapeutics. Our studies reveal the contribution of low frequency variant patterns at heterologous sites across the protein to accurate prediction of the changes at each position of interest.

The molecular evolution of cancer associated genes in mammals

Cancer is a disease that many multicellular organisms have faced for millions of years, and species have evolved various tumour suppression mechanisms to control oncogenesis. Although cancer occurs across the tree of life, cancer related mortality risks vary across mammalian orders, with Carnivorans particularly affected. Evolutionary theory predicts different selection pressures on genes associated with cancer progression and suppression, including oncogenes, tumour suppressor genes and immune genes. Therefore, we investigated the evolutionary history of cancer associated gene sequences across 384 mammalian taxa, to detect signatures of selection across categories of oncogenes (GRB2, FGL2 and CDC42), tumour suppressors (LITAF, Casp8 and BRCA2) and immune genes (IL2, CD274 and B2M). This approach allowed us to conduct a fine scale analysis of gene wide and site-specific signatures of selection across mammalian lineages under the lens of cancer susceptibility. Phylogenetic analyses revealed that for most species the evolution of cancer associated genes follows the species' evolution. The gene wide selection analyses revealed oncogenes being the most conserved, tumour suppressor and immune genes having similar amounts of episodic diversifying selection. Despite BRCA2's status as a key caretaker gene, episodic diversifying selection was detected across mammals. The site-specific selection analyses revealed that the two apoptosis associated domains of the Casp8 gene of bats (Chiroptera) are under opposing forces of selection (positive and negative respectively), highlighting the importance of site-specific selection analyses to understand the evolution of highly complex gene families. Our results highlighted the need to critically assess different types of selection pressure on cancer associated genes when investigating evolutionary adaptations to cancer across the tree of life. This study provides an extensive assessment of cancer associated genes in mammals with highly representative, and substantially large sample size for a comparative genomic analysis in the field and identifies various avenues for future research into the mechanisms of cancer resistance and susceptibility in mammals.

Reconstruction of Avian Reovirus History and Dispersal Patterns: A Phylodynamic Study

Avian reovirus (ARV) infection can cause significant losses to the poultry industry. Disease control has traditionally been attempted mainly through vaccination. However, the increase in clinical outbreaks in the last decades demonstrated the poor effectiveness of current vaccination approaches. The present study reconstructs the evolution and molecular epidemiology of different ARV genotypes using a phylodynamic approach, benefiting from a collection of more than one thousand sigma C (σC) sequences sampled over time at a worldwide level. ARVs' origin was estimated to occur several centuries ago, largely predating the first clinical reports. The origins of all genotypes were inferred at least one century ago, and their emergence and rise reflect the intensification of the poultry industry. The introduction of vaccinations had only limited and transitory effects on viral circulation and further expansion was observed, particularly after the 1990s, likely because of the limited immunity and the suboptimal and patchy vaccination application. In parallel, strong selective pressures acted with different strengths and directionalities among genotypes, leading to the emergence of new variants. While preventing the spread of new variants with different phenotypic features would be pivotal, a phylogeographic analysis revealed an intricate network of viral migrations occurring even over long distances and reflecting well-established socio-economic relationships.

Complete mitochondrial genome of critically endangered catfish Hemibagrus punctatus (Jerdon, 1849) and comparative analysis for insights into the phylogeny of hemibagrids through mitogenomic approach

BACKGROUND

Hemibagrus punctatus (Jerdon, 1849) is a critically endangered bagrid catfish endemic to the Western Ghats of India, whose population is declining due to anthropogenic activities. The current study aims to compare the mitogenome of H. punctatus with that of other Bagrid catfishes and provide insights into their evolutionary relationships.

METHODS AND RESULTS

Samples were collected from Hemmige Karnataka, India. In the present study, the mitogenome of H. punctatus was successfully assembled, and its phylogenetic relationships with other Bagridae species were studied. The total genomic DNA of samples was extracted following the phenol-chloroform isoamyl alcohol method. Samples were sequenced, and the Illumina paired-end reads were assembled to a contig length of 16,517 bp. The mitochondrial genome was annotated using MitoFish and MitoAnnotator (Iwasaki et al., 2013). A robust phylogenetic analysis employing NJ (Maximum composite likelihood) and ASAP methods supports the classification of H. punctatus within the Bagridae family, which validates the taxonomic status of this species. In conclusion, this research enriches our understanding of H. punctatus mitogenome, shedding light on its evolutionary dynamics within the Bagridae family and contributing to the broader knowledge of mitochondrial genes in the context of evolutionary biology.

CONCLUSIONS

The study's findings contribute to a better understanding of the mitogenome of H. punctatus and provide insights into the evolutionary relationships within other Hemibagrids.

A survey of the Sli gene in wild and cultivated potato

Inbred-hybrid breeding of diploid potatoes necessitates breeding lines that are self-compatible. One way of incorporating self-compatibility into incompatible cultivated potato (Solanum tuberosum) germplasm is to introduce the S-locus inhibitor gene (Sli), which functions as a dominant inhibitor of gametophytic self-incompatibility. To learn more about Sli diversity and function in wild species relatives of cultivated potato, we obtained Sli gene sequences that extended from the 5'UTR to the 3'UTR from 133 individuals from 22 wild species relatives of potato and eight diverse cultivated potato clones. DNA sequence alignment and phylogenetic trees based on genomic and protein sequences show that there are two highly conserved groups of Sli sequences. DNA sequences in one group contain the 533 bp insertion upstream of the start codon identified previously in self-compatible potato. The second group lacks the insertion. Three diploid and four polyploid individuals of wild species collected from geographically disjointed localities contained Sli with the 533 bp insertion. For most of the wild species clones examined, however, Sli did not have the insertion. Phylogenetic analysis indicated that Sli sequences with the insertion, in wild species and in cultivated clones, trace back to a single origin. Some diploid wild potatoes that have Sli with the insertion were self-incompatible and some wild potatoes that lack the insertion were self-compatible. Although there is evidence of positive selection for some codon positions in Sli, there is no evidence of diversifying selection at the gene level. In silico analysis of Sli protein structure did not support the hypothesis that amino acid changes from wild-type (no insertion) to insertion-type account for changes in protein function. Our study demonstrated that genetic factors besides the Sli gene must be important for conditioning a switch in the mating system from self-incompatible to self-compatible in wild potatoes.