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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

High allelic diversity in Arabidopsis NLRs is associated with distinct genomic features

Plants rely on Nucleotide-binding, Leucine-rich repeat Receptors (NLRs) for pathogen recognition. Highly variable NLRs (hvNLRs) show remarkable intraspecies diversity, while their low-variability paralogs (non-hvNLRs) are conserved between ecotypes. At a population level, hvNLRs provide new pathogen-recognition specificities, but the association between allelic diversity and genomic and epigenomic features has not been established. Our investigation of NLRs in Arabidopsis Col-0 has revealed that hvNLRs show higher expression, less gene body cytosine methylation, and closer proximity to transposable elements than non-hvNLRs. hvNLRs show elevated synonymous and nonsynonymous nucleotide diversity and are in chromatin states associated with an increased probability of mutation. Diversifying selection maintains variability at a subset of codons of hvNLRs, while purifying selection maintains conservation at non-hvNLRs. How these features are established and maintained, and whether they contribute to the observed diversity of hvNLRs is key to understanding the evolution of plant innate immune receptors.

Evolutionary characteristics of the mitochondrial NADH dehydrogenase subunit 6 gene in some populations of four sympatric Mustela species (Mustelidae, Mammalia) from central Europe

BACKGROUND

Selection on or reticulate evolution of mtDNA is documented in various mammalian taxa and could lead to misleading phylogenetic conclusions if not recognized. We sequenced the MT-ND6 gene of four sympatric Mustelid species of the genus Mustela from some central European populations. We hypothesised positive selection on MT-ND6, given its functional importance and the different body sizes and life histories of the species, even though climatic differences may be unimportant for adaptation in sympatry.

METHODS AND RESULTS

MT-ND6 genes were sequenced in 187 sympatric specimens of weasels, Mustela nivalis, stoats, M. erminea, polecats, M. putorius, and steppe polecats, M. eversmannii, from eastern Austria and of fourteen allopatric polecats from eastern-central Germany. Median joining networks, neighbour joining and maximum likelihood analyses as well as Bayesian inference grouped all species according to earlier published phylogenetic models. However, polecats and steppe polecats, two very closely related species, shared the same two haplotypes. We found only negative selection within the Mustela sequences, including 131 downloaded ones covering thirteen species. Positive selection was observed on three MT-ND6 codons of other mustelid genera retrieved from GenBank.

CONCLUSIONS

Negative selection for MT-ND6 within the genus Mustela suggests absence of both environmental and species-specific effects of cellular energy metabolism despite large species-specific differences in body size. The presently found shared polymorphism in European polecats and steppe polecats may result from ancestral polymorphism before speciation and historical or recent introgressive hybridization; it may indicate mtDNA capture of steppe polecats by M. putorius in Europe.

Genomic signatures of exceptional longevity and negligible aging in the long-lived red sea urchin

The red sea urchin (Mesocentrotus franciscanus) is one of the Earth's longest-living animals, reported to live more than 100 years with indeterminate growth, life-long reproduction, and no increase in mortality rate with age. To understand the genetic underpinnings of longevity and negligible aging, we constructed a chromosome-level assembly of the red sea urchin genome and compared it to that of short-lived sea urchin species. Genome-wide syntenic alignments identified chromosome rearrangements that distinguish short- and long-lived species. Expanded gene families in long-lived species play a role in innate immunity, sensory nervous system, and genome stability. An integrated network of genes under positive selection in the red sea urchin was involved in genomic regulation, mRNA fidelity, protein homeostasis, and mitochondrial function. Our results implicated known longevity genes in sea urchin longevity but also revealed distinct molecular signatures that may promote long-term maintenance of tissue homeostasis, disease resistance, and negligible aging.

Accurate Detection of Convergent Mutations in Large Protein Alignments With ConDor

Evolutionary convergences are observed at all levels, from phenotype to DNA and protein sequences, and changes at these different levels tend to be correlated. Notably, convergent mutations can lead to convergent changes in phenotype, such as changes in metabolism, drug resistance, and other adaptations to changing environments. We propose a two-component approach to detect mutations subject to convergent evolution in protein alignments. The "Emergence" component selects mutations that emerge more often than expected, while the "Correlation" component selects mutations that correlate with the convergent phenotype under study. With regard to Emergence, a phylogeny deduced from the alignment is provided by the user and is used to simulate the evolution of each alignment position. These simulations allow us to estimate the expected number of mutations in a neutral model, which is compared to the observed number of mutations in the data studied. In Correlation, a comparative phylogenetic approach, is used to measure whether the presence of each of the observed mutations is correlated with the convergent phenotype. Each component can be used on its own, for example Emergence when no phenotype is available. Our method is implemented in a standalone workflow and a webserver, called ConDor. We evaluate the properties of ConDor using simulated data, and we apply it to three real datasets: sedge PEPC proteins, HIV reverse transcriptase, and fish rhodopsin. The results show that the two components of ConDor complement each other, with an overall accuracy that compares favorably to other available tools, especially on large datasets.

Analyses and Insights into Genetic Reassortment and Natural Selection as Key Drivers of Piscine orthoreovirus Evolution

Piscine orthoreovirus (PRV) is a pathogen that causes heart and skeletal muscle inflammation in Salmo salar and has also been linked to circulatory disorders in other farmed salmonids, such as Oncorhynchus kisutch and Oncorhynchus mykiss. The virus has a segmented, double-stranded RNA genome, which makes it possible to undergo genetic reassortment and increase its genomic diversity through point mutations. In this study, genetic reassortment in PRV was assessed using the full genome sequences available in public databases. This study used full genome sequences that were concatenated and genome-wide reassortment events, and phylogenetic analyses were performed using the recombination/reassortment detection program version 5 (RDP5 V 5.5) software. Additionally, each segment was aligned codon by codon, and overall mean distance and selection was tested using the Molecular Evolutionary Genetics Analysis X software, version 10.2 (MEGA X version 10.2). The results showed that there were 17 significant reassortment events in 12 reassortant sequences, involving genome exchange between low and highly virulent genotypes. PRV sequences from different salmonid host species did not appear to limit the reassortment. This study found that PRV frequently undergoes reassortment events to increase the diversity of its segmented genome, leading to antigenic variation and increased virulence. This study also noted that to date, no reassortment events have been described between PRV-1 and PRV-3 genotypes. However, the number of complete genomic sequences within each genotype is uneven. This is important because PRV-3 induces cross-protection against PRV-1, making it a potential vaccine candidate.

Dynamic co-evolution of transposable elements and the piRNA pathway in African cichlid fishes

East African cichlid fishes have diversified in an explosive fashion, but the (epi)genetic basis of the phenotypic diversity of these fishes remains largely unknown. Although transposable elements (TEs) have been associated with phenotypic variation in cichlids, little is known about their transcriptional activity and epigenetic silencing. Here, we describe dynamic patterns of TE expression in African cichlid gonads and during early development. Orthology inference revealed an expansion of piwil1 genes in Lake Malawi cichlids, likely driven by PiggyBac TEs. The expanded piwil1 copies have signatures of positive selection and retain amino acid residues essential for catalytic activity. Furthermore, the gonads of African cichlids express a Piwi-interacting RNA (piRNA) pathway that target TEs. We define the genomic sites of piRNA production in African cichlids and find divergence in closely related species, in line with fast evolution of piRNA-producing loci. Our findings suggest dynamic co-evolution of TEs and host silencing pathways in the African cichlid radiations. We propose that this co-evolution has contributed to cichlid genomic diversity.

Evidence for Selection on Mitochondrial OXPHOS Genes in the Mediterranean Killifish Aphanius fasciatus Valenciennes, 1821

Mitochondrial oxidative phosphorylation (OXPHOS) genes are a system subject to selection under determined environmental constraints despite a neutral evolution model that has long been hypothesized for the mitochondrial genome. In this study, the sequences of ND1, Cytb, and COI OXPHOS genes were analyzed in six populations of the eurythermal and euryhaline killifish A. fasciatus, to detect non-synonymous mutations leading to amino acid changes and to check whether selection acted on them using tests of recombination and selection. The results indicate a high COI and Cytb gene diversity and a high percentage of private haplotypes in all populations. In the Greek population, non-synonymous nucleotide substitutions were observed in the N-terminal region of COI and Cytb. Positively selected sites were also found. The information we obtained from the mitochondrial DNA sequences of A. fasciatus adds to the growing data on selective pressure acting on mitochondrial DNA in non-model species. These results should be explored from the perspective of the local adaptation of eurythermal and euryhaline species and supported using experimental evidence to better understand the interplay between historical climatic events and local adaptation and how each of them contributes to shaping the genetic structure of this species.

The Effect of Global Spread, Epidemiology, and Control Strategies on the Evolution of the GI-19 Lineage of Infectious Bronchitis Virus

The GI-19 lineage of infectious bronchitis virus (IBV) has emerged as one of the most impactful, particularly in the "Old World". Originating in China several decades ago, it has consistently spread and evolved, often forming independent clades in various areas and countries, each with distinct production systems and control strategies. This study leverages this scenario to explore how different environments may influence virus evolution. Through the analysis of the complete S1 sequence, four datasets were identified, comprising strains of monophyletic clades circulating in different continents or countries (e.g., Asia vs. Europe and China vs. Thailand), indicative of single introduction events and independent evolution. The population dynamics and evolutionary rate variation over time, as well as the presence and intensity of selective pressures, were estimated and compared across these datasets. Since the lineage origin (approximately in the mid-20th century), a more persistent and stable viral population was estimated in Asia and China, while in Europe and Thailand, a sharp increase following the introduction (i.e., 2005 and 2007, respectively) of GI-19 was observed, succeeded by a rapid decline. Although a greater number of sites on the S1 subunit were under diversifying selection in the Asian and Chinese datasets, more focused and stronger pressures were evident in both the European (positions 2, 52, 54, 222, and 379 and Thai (i.e., positions 10, 12, 32, 56, 62, 64, 65, 78, 95, 96, 119, 128, 140, 182, 292, 304, 320, and 323) strains, likely reflecting a more intense and uniform application of vaccines in these regions. This evidence, along with the analysis of control strategies implemented in different areas, suggests a strong link between effective, systematic vaccine implementation and infection control. However, while the overall evolutionary rate was estimated at approximately 10-3 to 10-4, a significant inverse correlation was found between viral population size and the rate of viral evolution over time. Therefore, despite the stronger selective pressure imposed by vaccination, effectively constraining the former through adequate control strategies can efficiently prevent viral evolution and the emergence of vaccine-escaping variants.

Comparative evolutionary analyses of peste des petits ruminants virus genetic lineages

Peste des petits ruminants virus (PPRV) causes a highly infectious disease affecting mainly goats and sheep in large parts of Africa, Asia, and the Middle East and has an important impact on the global economy and food security. Full genome sequencing of PPRV strains has proved to be critical to increasing our understanding of PPR epidemiology and to inform the ongoing global efforts for its eradication. However, the number of full PPRV genomes published is still limited and with a heavy bias towards recent samples and genetic Lineage IV (LIV), which is only one of the four existing PPRV lineages. Here, we generated genome sequences for twenty-five recent (2010-6) and seven historical (1972-99) PPRV samples, focusing mainly on Lineage II (LII) in West Africa. This provided the first opportunity to compare the evolutionary pressures and history between the globally dominant PPRV genetic LIV and LII, which is endemic in West Africa. Phylogenomic analysis showed that the relationship between PPRV LII strains was complex and supported the extensive transboundary circulation of the virus within West Africa. In contrast, LIV sequences were clearly separated per region, with strains from West and Central Africa branched as a sister clade to all other LIV sequences, suggesting that this lineage also has an African origin. Estimates of the time to the most recent common ancestor place the divergence of modern LII and LIV strains in the 1960s-80s, suggesting that this period was particularly important for the diversification and spread of PPRV globally. Phylogenetic relationships among historical samples from LI, LII, and LIII and with more recent samples point towards a high genetic diversity for all these lineages in Africa until the 1970s-80s and possible bottleneck events shaping PPRV's evolution during this period. Molecular evolution analyses show that strains belonging to LII and LIV have evolved under different selection pressures. Differences in codon usage and adaptative selection pressures were observed in all viral genes between the two lineages. Our results confirm that comparative genomic analyses can provide new insights into PPRV's evolutionary history and molecular epidemiology. However, PPRV genome sequencing efforts must be ramped up to increase the resolution of such studies for their use in the development of efficient PPR control and surveillance strategies.

Genome-wide analysis of cellulose synthase gene superfamily in Tectona grandis L.f

This study aimed to explore Cellulose synthase gene superfamily of teak, and its evolutionary relationship with homologous genes of other woody species. The incidence of evolutionary events like gene duplication and gene loss, influence of the selection pressure, and consequent adaptive functional divergence of the duplicated TgCes gene were assessed alongside it's role in wood coloration. This study identified 39 full-length non-redundant proteins belonging to CesA and Csl gene families. TgCesA and TgCsl proteins with Cellulose synthase domain repeats indicated tandem gene duplication and probable genetic variability, enabling local adaptation. Further, multi-domain protein (MYB-like DNA-binding domain and CesA domain) with maximum introns was also identified indicating gene fusion and formation of complex protein with novel functions. Phylogenetic analysis grouped the genes into seven subfamilies (CesA, CslA, CslC, CslD, CslE, CslG, and CslM) with each undergoing gene duplication and loss along their evolutionary history. Post-species gene duplications and probable neofunctionalization were identified in TgCesA and TgCsl gene families. Each subfamily was found to be under strong purifying selection with a few or no sites under positive selection. Functional divergence analysis further revealed site-specific selective constraints in CesA and Csl genes of the teak Cellulose synthase gene family. Furthermore, protein-protein interaction network analysis identified co-expression of Cellulose synthase gene with flavonoid 3',5'-hydroxylase (F3'5'H, CYP75A), involved in the biosynthesis of xylem anthocyanin compounds, probably responsible for wood coloration. This study thus offers a foundation for future research in wood formation and wood property traits specific to teak and its provenances.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03927-6.

Evolutionary analysis of major histocompatibility complex variants in chytrid-resistant and susceptible amphibians

An amphibian emerging infectious disease (EID), chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), originated in Asia but primarily led to declines and extinctions in amphibian populations outside of Asia. Host major histocompatibility complex (MHC) molecules exhibit high polymorphism, and the evolution of MHC can be influenced by recombination and pathogens. Previous studies have indicated that host MHC class II is associated with Bd resistance. In this study, I conducted recombination and selection tests on functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) and an Australasian Bd-susceptible species (Litoria caerulea). Recombination at the same site was identified in both species, supporting the hypothesis that recombination contributes to MHC IIß1 diversity in amphibians. Positive selection was observed in MHC IIß1 alleles in both species. In L. caerulea, at least four amino acid sites were identified under significant positive selection in the MHC IIß1, whereas these sites were either negatively selected or conserved in B. gargarizans. This suggests these sites might be selected for Bd resistance. Hydrophobicity was detected in certain amino acid sites relating to Bd resistance, suggesting this physicochemical property may be a factor selected to counteract Bd infection. These findings of this study provide an evolutionary basis for understanding how amphibian MHC IIß1 may undergo selection in response to chytrid infection.

Evaluating the genetic diversity of the Plasmodium vivax siap2 locus: A promising candidate for an effective malaria vaccine?

Malaria is the deadliest parasitic disease in the world. Traditional control measures have become less effective; hence, there is a need to explore alternative strategies, such as antimalarial vaccines. However, designing an anti-Plasmodium vivax vaccine is considered a challenge due to the complex parasite biology and the antigens' high genetic diversity. Recently, the sporozoite invasion-associated protein 2 (SIAP2) has been suggested as a potential antigen to be considered in vaccine design due to its significance during hepatocyte invasion. However, its use may be limited by the incomplete understanding of gene/protein diversity. Here, the genetic diversity of pvsiap2 using P. vivax DNA samples from Colombia was assessed. Through PCR amplification and sequencing, we compared the Colombian sequences with available worldwide sequences, revealing that pvsiap2 displays low genetic diversity. Molecular evolutionary analyses showed that pvsiap2 appears to be influenced by directional selection. Moreover, the haplotypes found differ by a few mutational steps and several of them were shared between different geographical areas. On the other hand, several conserved regions within PvSIAP2 were predicted as potential B-cell or T-cell epitopes. Considering these characteristics and its role in hepatocyte invasion, the PvSIAP2 protein emerges as a promising antigen to be considered in a multi-antigen-multi-stage (multivalent) fully effective vaccine against P. vivax malaria.

Rapid evolution of genes with anti-cancer functions during the origins of large bodies and cancer resistance in elephants

Elephants have emerged as a model system to study the evolution of body size and cancer resistance because, despite their immense size, they have a very low prevalence of cancer. Previous studies have found that duplication of tumor suppressors at least partly contributes to the evolution of anti-cancer cellular phenotypes in elephants. Still, many other mechanisms must have contributed to their augmented cancer resistance. Here, we use a suite of codon-based maximum-likelihood methods and a dataset of 13,310 protein-coding gene alignments from 261 Eutherian mammals to identify positively selected and rapidly evolving elephant genes. We found 496 genes (3.73% of alignments tested) with statistically significant evidence for positive selection and 660 genes (4.96% of alignments tested) that likely evolved rapidly in elephants. Positively selected and rapidly evolving genes are statistically enriched in gene ontology terms and biological pathways related to regulated cell death mechanisms, DNA damage repair, cell cycle regulation, epidermal growth factor receptor (EGFR) signaling, and immune functions, particularly neutrophil granules and degranulation. All of these biological factors are plausibly related to the evolution of cancer resistance. Thus, these positively selected and rapidly evolving genes are promising candidates for genes contributing to elephant-specific traits, including the evolution of molecular and cellular characteristics that enhance cancer resistance.

Ancient reindeer mitogenomes reveal island-hopping colonisation of the Arctic archipelagos

Climate warming at the end of the last glacial period had profound effects on the distribution of cold-adapted species. As their range shifted towards northern latitudes, they were able to colonise previously glaciated areas, including remote Arctic islands. However, there is still uncertainty about the routes and timing of colonisation. At the end of the last ice age, reindeer/caribou (Rangifer tarandus) expanded to the Holarctic region and colonised the archipelagos of Svalbard and Franz Josef Land. Earlier studies have proposed two possible colonisation routes, either from the Eurasian mainland or from Canada via Greenland. Here, we used 174 ancient, historical and modern mitogenomes to reconstruct the phylogeny of reindeer across its whole range and to infer the colonisation route of the Arctic islands. Our data shows a close affinity among Svalbard, Franz Josef Land and Novaya Zemlya reindeer. We also found tentative evidence for positive selection in the mitochondrial gene ND4, which is possibly associated with increased heat production. Our results thus support a colonisation of the Eurasian Arctic archipelagos from the Eurasian mainland and provide some insights into the evolutionary history and adaptation of the species to its High Arctic habitat.

Multiple Lines of Evidence Support 199 SARS-CoV-2 Positively Selected Amino Acid Sites

SARS-CoV-2 amino acid variants that contribute to an increased transmissibility or to host immune system escape are likely to increase in frequency due to positive selection and may be identified using different methods, such as codeML, FEL, FUBAR, and MEME. Nevertheless, when using different methods, the results do not always agree. The sampling scheme used in different studies may partially explain the differences that are found, but there is also the possibility that some of the identified positively selected amino acid sites are false positives. This is especially important in the context of very large-scale projects where hundreds of analyses have been performed for the same protein-coding gene. To account for these issues, in this work, we have identified positively selected amino acid sites in SARS-CoV-2 and 15 other coronavirus species, using both codeML and FUBAR, and compared the location of such sites in the different species. Moreover, we also compared our results to those that are available in the COV2Var database and the frequency of the 10 most frequent variants and predicted protein location to identify those sites that are supported by multiple lines of evidence. Amino acid changes observed at these sites should always be of concern. The information reported for SARS-CoV-2 can also be used to identify variants of concern in other coronaviruses.

Diversity of selected toll-like receptor genes in cheetahs (Acinonyx jubatus) and African leopards (Panthera pardus pardus)

The anthropogenic impact on wildlife is ever increasing. With shrinking habitats, wild populations are being pushed to co-exist in proximity to humans leading to an increased threat of infectious diseases. Therefore, understanding the immune system of a species is key to assess its resilience in a changing environment. The innate immune system (IIS) is the body's first line of defense against pathogens. High variability in IIS genes, like toll-like receptor (TLR) genes, appears to be associated with resistance to infectious diseases. However, few studies have investigated diversity in TLR genes in vulnerable species for conservation. Large predators are threatened globally including leopards and cheetahs, both listed as 'vulnerable' by IUCN. To examine IIS diversity in these sympatric species, we used next-generation-sequencing to compare selected TLR genes in African leopards and cheetahs. Despite differences, both species show some TLR haplotype similarity. Historic cheetahs from all subspecies exhibit greater genetic diversity than modern Southern African cheetahs. The diversity in investigated TLR genes is lower in modern Southern African cheetahs than in African leopards. Compared to historic cheetah data and other subspecies, a more recent population decline might explain the observed genetic impoverishment of TLR genes in modern Southern African cheetahs. However, this may not yet impact the health of this cheetah subspecies.

Genomic Analyses of Major SARS-CoV-2 Variants Predicting Multiple Regions of Pathogenic and Transmissive Importance

The rapid evolution of SARS-CoV-2 has fueled its global proliferation since its discovery in 2019, with several notable variants having been responsible for increases in cases of coronavirus disease 2019 (COVID-19). Analyses of codon bias and usage in these variants between phylogenetic clades or lineages may grant insights into the evolution of SARS-CoV-2 and identify target codons indicative of evolutionary or mutative trends that may prove useful in tracking or defending oneself against emerging strains. We processed a cohort of 120 SARS-CoV-2 genome sequences through a statistical and bioinformatic pipeline to identify codons presenting evidence of selective pressure as well as codon coevolution. We report the identification of two codon sites in the orf8 and N genes demonstrating such evidence with real-world impacts on pathogenicity and transmissivity.

Genotype F of Echovirus 25 with multiple recombination pattern have been persistently and extensively circulating in Chinese mainland

Echovirus 25 (E25), a member of the Enterovirus B (EV-B) species, can cause aseptic meningitis (AM), viral meningitis (VM), and acute flaccid paralysis (AFP). However, systematic studies on the molecular epidemiology of E25, especially those concerning its evolution and recombination, are lacking. In this study, 18 strains of E25, isolated from seven provinces of China between 2009 and 2018, were collected based on the Chinese hand, foot, and mouth disease (HFMD) surveillance network, and 95 sequences downloaded from GenBank were also screened. Based on the phylogenetic analysis of 113 full-length VP1 sequences worldwide, globally occurring E25 strains were classified into 9 genotypes (A-I), and genotype F was the dominant genotype in the Chinese mainland. The average nucleotide substitution rate of E25 was 6.08 × 10-3 substitutions/site/year, and six important transmission routes were identified worldwide. Seventeen recombination patterns were determined, of which genotype F can be divided into 9 recombination patterns. A positive selector site was found in the capsid protein region of genotype F. Recombination analysis and pressure selection analysis for genotype F showed multiple recombination patterns and evolution characteristics, which may be responsible for it being the dominant genotype in the Chinese mainland. This study provides a theoretical basis for the subsequent prevention and control of E25.

Genetic diversity and natural selection of apical membrane antigen-1 (ama-1) in Cameroonian Plasmodium falciparum isolates

Antigenic variation associated with genetic diversity in global Plasmodium falciparum apical membrane antigen-1 (PfAMA-1) is a major impediment to designing an effective malaria vaccine. Here, we report the first study on genetic diversity and natural selection of the Pfama-1 gene in P. falciparum isolates from Cameroon. A total of 328 P. falciparum positive samples collected during 2016 and 2019 from five localities of Cameroon were analysed. The ectodomain coding fragment of Pfama-1 gene was amplified for polymorphism profiling and natural selection analysis. A total of 108 distinct haplotypes were found in 203 P. falciparum isolates with considerable nucleotide diversity (π = 0.016) and haplotype diversity (Hd = 0.976). Most amino acid substitutions detected were scattered in ectodomain-I and few specific mutations viz P145L, K148Q, K462I, L463F, N471K, S482L, E537G, K546R and I547F were seen only in Cameroonian isolates. A tendency of natural selection towards positive diversifying selection was observed (Taj-D = 2.058). Five positively selected codon sites (P145L, S283L, Q308E/K, P330S and I547F) were identified, which overlapped with predicted B-cell epitopes and red blood cell (RBC) binding sites, suggesting their potential implication in host immune pressure and parasite-RBC binding complex modulation. The Cameroonian P. falciparum populations indicated a moderate level of genetic differentiation when compared with global sequences, with few exceptions from Vietnam and Venezuela. Our findings provide baseline data on existing Pfama-1 gene polymorphisms in Cameroonian field isolates, which will be useful information for malaria vaccine design.

Evolutionary Analysis of Cnidaria Small Cysteine-Rich Proteins (SCRiPs), an Enigmatic Neurotoxin Family from Stony Corals and Sea Anemones (Anthozoa: Hexacorallia)

Cnidarians (corals, sea anemones, and jellyfish) produce toxins that play central roles in key ecological processes, including predation, defense, and competition, being the oldest extant venomous animal lineage. Cnidaria small cysteine-rich proteins (SCRiPs) were the first family of neurotoxins detected in stony corals, one of the ocean's most crucial foundation species. Yet, their molecular evolution remains poorly understood. Moreover, the lack of a clear classification system has hindered the establishment of an accurate and phylogenetically informed nomenclature. In this study, we extensively surveyed 117 genomes and 103 transcriptomes of cnidarians to identify orthologous SCRiP gene sequences. We annotated a total of 168 novel putative SCRiPs from over 36 species of stony corals and 12 species of sea anemones. Phylogenetic reconstruction identified four distinct SCRiP subfamilies, according to strict discrimination criteria based on well-supported monophyly with a high percentage of nucleotide and amino acids' identity. Although there is a high prevalence of purifying selection for most SCRiP subfamilies, with few positively selected sites detected, a subset of Acroporidae sequences is influenced by diversifying positive selection, suggesting potential neofunctionalizations related to the fine-tuning of toxin potency. We propose a new nomenclature classification system relying on the phylogenetic distribution and evolution of SCRiPs across Anthozoa, which will further assist future proteomic and functional research efforts.

Cellular receptors for mammalian viruses

The interaction of viral surface components with cellular receptors and other entry factors determines key features of viral infection such as host range, tropism and virulence. Despite intensive research, our understanding of these interactions remains limited. Here, we report a systematic analysis of published work on mammalian virus receptors and attachment factors. We build a dataset twice the size of those available to date and specify the role of each factor in virus entry. We identify cellular proteins that are preferentially used as virus receptors, which tend to be plasma membrane proteins with a high propensity to interact with other proteins. Using machine learning, we assign cell surface proteins a score that predicts their ability to function as virus receptors. Our results also reveal common patterns of receptor usage among viruses and suggest that enveloped viruses tend to use a broader repertoire of alternative receptors than non-enveloped viruses, a feature that might confer them with higher interspecies transmissibility.

Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid?

Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?

Patterns of evolution in MHC class II DQA and DQB exon 2 genes of Alpine mountain hares, Lepus timidus varronis, and sympatric and parapatric brown hares, L. europaeus, from Switzerland

In natural populations, hybridization is known to occur between a wide range of species. However, its evolutionary significance is less clear. Genes involved in fighting pathogens are considered excellent candidates for studying adaptive introgression, although both introgression and balancing selection can generate similar patterns of diversity and differentiation. Here, we compared DQA and DQB MHC class II and microsatellite allelic diversity of sympatric and parapatric mountain (Lepus timidus) and brown hare (L. europaeus) populations from Switzerland. We detected higher genetic diversity in brown hares compared to mountain hares at both MHC and microsatellite loci. We consider the observed patterns of microsatellite diversity both for L. europaeus and L. timidus as result of stochastic demographic processes while the pattern of MHC polymorphism of the studied hare populations can be explained by pathogen-driven selection. Rare bidirectional gene flow between both hare species seems to occur specifically for MHC alleles. However, the high number of shared alleles showing similar high frequency in both species suggests that reciprocally exchanged MHC alleles are being maintained via balancing selection. Adaptation to similar pathogen communities can also lead to parallel selection of MHC alleles. Positive selection, recombination and mutations have played different roles in shaping the patterns of MHC allelic diversity in and differentiation between both species. Results for the latter evolutionary forces do not show a better matching between the sympatric populations compared to the parapatric ones, suggesting a minor role of introgression for the observed evolutionary patterns of the studied hare species.

Babesia bovis RON2 binds to bovine erythrocytes through a highly conserved epitope

B. bovis invasion of bovine erythrocytes requires tight junction formation involving AMA-1/RON2 complex interaction. RON2 has been considered a vaccine candidate since antibodies targeting the protein can inhibit parasite invasion of target cells; however, the mechanism controlling B. bovis RON2 interaction with red blood cells is not yet fully understood. This study was thus aimed at identifying B. bovis RON2 protein regions associated with interaction with bovine erythrocytes. Natural selection analysis of the ron2 gene identified predominantly negative selection signals in the C-terminal region. Interestingly, protein-cell and competition assays highlighted the RON2-C region's role in peptide 42918-mediated erythrocyte binding, probably to a sialoglycoprotein receptor. This peptide (1218SFIMVKPPALHCVLKPVETL1237) lies within an intrinsically disordered region of the RON2 secondary structure flanked by two helical residues. The study provides, for the first time, valuable insights into RON2's role in interaction with its target cells. Future studies are required for studying the peptide's potential as an anti-B. bovis vaccine component.

Adaptive evolution of mitochondrial genomes in Triplophysa cavefishes

Extreme conditions in caves pose survival challenges for cave dwellers, who gradually develop adaptive survival features. Cavefishes are one of the most successful animals among cave dwellers. Triplophysa cavefishes are an important group of cavefishes, and they show remarkable adaptability to the extreme environments of caves. However, there is a limited understanding of their adaptation mechanisms. In this study, eight complete mitochondrial genomes of Triplophysa cavefishes were newly obtained, and their genomic characteristics, including the base composition, base bias, and codon usage, were analyzed. Phylogenetic analysis was carried out based on 13 mitochondrial protein-coding genes from 44 Nemacheilidae species. This showed that Triplophysa cavefishes and non-cavefishes separate into two reciprocally monophyletic clades, suggesting a single origin of the cave phenotype. Positive selection analysis strongly suggested that the selection pressure in cavefishes is higher than that in non-cavefishes. Furthermore, the ND5 gene in cavefishes showed evidence of positive selection, which suggests that the gene may play an important role in the adaptation of cavefishes to the cave environment. Protein structure analysis of the ND5 subunit implied that the sites of positive selection in cavefishes might allow them to acquire lower ND5 protein stability, compared to that in non-cavefishes, which might help the accumulation of nonsynonymous (mildly deleterious) mutations. Together, our study revealed the genetic signatures of cave adaptation in Triplophysa cavefishes from the perspective of energy metabolism.

Decoupled evolution of the Sex Peptide gene family and Sex Peptide Receptor in Drosophilidae

Across internally fertilising species, males transfer ejaculate proteins that trigger wide-ranging changes in female behaviour and physiology. Much theory has been developed to explore the drivers of ejaculate protein evolution. The accelerating availability of high-quality genomes now allows us to test how these proteins are evolving at fine taxonomic scales. Here, we use genomes from 264 species to chart the evolutionary history of Sex Peptide (SP), a potent regulator of female post-mating responses in Drosophila melanogaster. We infer that SP first evolved in the Drosophilinae subfamily and has since followed markedly different evolutionary trajectories in different lineages. Outside of the Sophophora-Lordiphosa, SP exists largely as a single-copy gene with independent losses in several lineages. Within the Sophophora-Lordiphosa, the SP gene family has repeatedly and independently expanded. Up to seven copies, collectively displaying extensive sequence variation, are present in some species. Despite these changes, SP expression remains restricted to the male reproductive tract. Alongside, we document considerable interspecific variation in the presence and morphology of seminal microcarriers that, despite the critical role SP plays in microcarrier assembly in D. melanogaster, appears to be independent of changes in the presence/absence or sequence of SP. We end by providing evidence that SP's evolution is decoupled from that of its receptor, Sex Peptide Receptor, in which we detect no evidence of correlated diversifying selection. Collectively, our work describes the divergent evolutionary trajectories that a novel gene has taken following its origin and finds a surprisingly weak coevolutionary signal between a supposedly sexually antagonistic protein and its receptor.

Equine herpesvirus 4 infected domestic horses associated with Sintashta spoke-wheeled chariots around 4,000 years ago

Equine viral outbreaks have disrupted the socio-economic life of past human societies up until the late 19th century and continue to be of major concern to the horse industry today. With a seroprevalence of 60-80 per cent, equine herpesvirus 4 (EHV-4) is the most common horse pathogen on the planet. Yet, its evolutionary history remains understudied. Here, we screen the sequenced data of 264 archaeological horse remains to detect the presence of EHV-4. We recover the first ancient EHV-4 genome with 4.2× average depth-of-coverage from a specimen excavated in the Southeastern Urals and dated to the Early Bronze Age period, approximately 3,900 years ago. The recovery of an EHV-4 virus outside the upper respiratory tract not only points to an animal particularly infected but also highlights the importance of post-cranial bones in pathogen characterisation. Bayesian phylogenetic reconstruction provides a minimal time estimate for EHV-4 diversification to around 4,000 years ago, a time when modern domestic horses spread across the Central Asian steppes together with spoke-wheeled Sintashta chariots, or earlier. The analyses also considerably revise the diversification time of the two EHV-4 subclades from the 16th century based solely on modern data to nearly a thousand years ago. Our study paves the way for a robust reconstruction of the history of non-human pathogens and their impact on animal health.

High prevalence and genetic diversity of Treponema paraluisleporidarum isolates in European lagomorphs

IMPORTANCE

Syphilis is an ancient disease of humans and lagomorphs caused by two distinct but genetically closely related bacteria (>98% sequence identity based on the whole genome) of the genus Treponema. While human syphilis is well studied, little is known about the disease in the lagomorph host. Yet, comparative studies are needed to understand mechanisms in host-pathogen coevolution in treponematoses. Importantly, Treponema paraluisleporidarum-infected hare populations provide ample opportunity to study the syphilis-causing pathogen in a naturally infected model population without antibiotic treatment, data that cannot be obtained from syphilis infection in humans. We provide data on genetic diversity and are able to highlight various types of repetitions in one of the two hypervariable regions at the tp0548 locus that have not been described in the human syphilis-causing sister bacterium Treponema pallidum subsp. pallidum.

Genetic Characteristics of Wuxiang Virus in Shanxi Province, China

Wuxiang virus (WUXV) is the first sandfly-borne Phlebovirus isolated from Phlebotomus chinensis collected in China and has been established as a consistent viral presence in the local sandfly populations of both Wuxiang County and Yangquan City. However, its distribution in the Shanxi Province remains unclear. In this study, three novel WUXV strains were isolated from sandflies collected from Jiexiu City, Shanxi Province, China, in 2022. Subsequently, whole-genome sequences of these novel strains were generated using next-generation sequencing. The open reading frame (ORF) sequences of the WUXV strains from the three locations were subjected to gene analysis. Phylogenetic analysis revealed that WUXV belongs to two distinct clades with geographical differences. Strains from Wuxiang County and Yangquan City belonged to clade 1, whereas strains from Jiexiu City belonged to clade 2. Reassortment and recombination analyses indicated no gene reassortment or recombination between the two clades. However, four reassortments or recombination events could be detected in clade 1 strains. By aligning the amino acid sequences, eighty-seven mutation sites were identified between the two clades, with seventeen, sixty, nine, and one site(s) in the proteins RdRp, M, NSs, and N, respectively. Additionally, selection pressure analysis identified 17 positively selected sites across the entire genome of WUXV, with two, thirteen, one, and one site(s) in the proteins RdRp, M, NSs, and N, respectively. Notably, sites M-312 and M-340 in the M segment not only represented mutation sites but also showed positive selective pressure effects. These findings highlight the need for continuous nationwide surveillance of WUXV.

Comparative genomics reveals probable adaptations for xylose use in Thermoanaerobacterium saccharolyticum

Second-generation ethanol, a promising biofuel for reducing greenhouse gas emissions, faces challenges due to the inefficient metabolism of xylose, a pentose sugar. Overcoming this hurdle requires exploration of genes, pathways, and organisms capable of fermenting xylose. Thermoanaerobacterium saccharolyticum is an organism capable of naturally fermenting compounds of industrial interest, such as xylose, and understanding evolutionary adaptations may help to bring novel genes and information that can be used for industrial yeast, increasing production of current bio-platforms. This study presents a deep evolutionary study of members of the firmicutes clade, focusing on adaptations in Thermoanaerobacterium saccharolyticum that may be related to overall fermentation metabolism, especially for xylose fermentation. One highlight is the finding of positive selection on a xylose-binding protein of the xylFGH operon, close to the annotated sugar binding site, with this protein already being found to be expressed in xylose fermenting conditions in a previous study. Results from this study can serve as basis for searching for candidate genes to use in industrial strains or to improve Thermoanaerobacterium saccharolyticum as a new microbial cell factory, which may help to solve current problems found in the biofuels' industry.

Studying SARS-CoV-2 interactions using phage-displayed receptor binding domain as a model protein

SARS-CoV-2 receptor binding domain (RBD) mediates viral entry into human cells through its interaction with angiotensin converting enzyme 2 (ACE2). Most neutralizing antibodies elicited by infection or vaccination target this domain. Such a functional relevance, together with large RBD sequence variability arising during viral spreading, point to the need of exploring the complex landscape of interactions between RBD-derived variants, ACE2 and antibodies. The current work was aimed at developing a simple platform to do so. Biologically active and antigenic Wuhan-Hu-1 RBD, as well as mutated RBD variants found in nature, were successfully displayed on filamentous phages. Mutational scanning confirmed the global plasticity of the receptor binding motif within RBD, highlighted residues playing a critical role in receptor binding, and identified mutations strengthening the interaction. The ability of vaccine-induced antibodies to inhibit ACE2 binding of many mutated RBD variants, albeit at different extents, was shown. Amino acid replacements which could compromise such inhibitory potential were underscored. The expansion of our approach could be the starting point for a large-scale phage-based exploration of diversity within RBD of SARS-CoV-2 and related coronaviruses, useful to understand structure-function relationships, to engineer RBD proteins, and to anticipate changes to watch during viral evolution.