Carbohydrate sulfotransferases (CHSTs) play a vital role in the production of sulfated polysaccharides (SPs) in algae by catalyzing the sulfation of carbohydrate moieties through the transfer of a sulfuryl group from the donor, 3'-phosphoadenosine 5'-phosphosulfate (PAPS). In the present study, putative algal CHSTs with a PF00685, PF03567. PF06990 and PF13469 domain were identified by HMMER search and Protein Basic Local Alignment Search Tool (BLAST) using the well-characterized human CHSTs as queries. Approximately half of the algal CHSTs that contained a PF00685 domain also possessed a PF13469 domain in an overlapping region. These CHSTs were structurally and phylogenetically distinct from algal CHSTs containing PF03567 or/and PF06990 domains. The PF00685/PF13469 domain is commonly found in Chlorophyta, while PF03567 and PF06990 domains are more prevalent in red algae and brown algae, respectively, reflecting the different types of SPs produced by these distinct phyla. Our phylogenetic analyses of algal CHSTs support the hypothesis of a polyphyletic origin, suggesting complex evolutionary histories involving both lineage-specific evolution and significant horizontal gene transfer (HGT) events between algae and organisms from other diverse taxa, including bacteria. In addition, the specificities of algal CHSTs for different carbohydrate moieties and site-specific sulfation patterns were inferred from the phylogenies of human CHSTs and the CHSTs from of algae with known SPs and chemical structures. This approach helps us to bridge the gap in knowledge, as a limited number of algal CHSTs have been biochemically characterized experimentally.

DNA codon mutations involving Stop signals or the amino acid cysteine can be especially damaging. The former can break protein sequences or add extraneous amino acids. The latter can add or subtract disulfide bonds crucial in protein folding. We present a hypothetical scenario where Stop codons were present early in the evolution of the genetic code; and minimizing catastrophic mutations for code networks affected all subsequent amino acid/codon pairings. Predicted features of this "Catastrophic Mutation Minimization Hypothesis" (CMMH) are that: (1) Cysteine is mutationally adjacent to Stop, isolating a contiguous codon 'neighborhood' with high potential for catastrophe. (2) The sequence of amino acid additions order determines codon assignments through minimizing network-wide mutation costs. Overall, codon locations for 16 of the 20 amino acids in the genetic code are consistent with the CMMH, as are multiple other predictions. We propose an antecedent genetic code consisted of 16 doublet codons specifying 13-14 amino acids. Two variations of these networks are less susceptible to catastrophic mutations than 88.2-97.5% of randomly generated ones. Unlike some previous hypotheses, CMMH does not require the total replacement or rearrangement of amino acids at codons, with its disruptive potential for protein synthesis. Finally, the composition of this ancestral doublet genetic code has all the modern code's utility: amino acids from four chemical types; start and stop signals; metal-binding ability; disulfide bridging for creating protein shapes; and possible epigenetic gene regulation. Thus, the modern code likely evolutionarily fine-tuned antecedent capabilities, rather than significantly increasing competence for making complex proteins.

The OrthoXML file is a standard file format for orthology data. It provides a standardized structure for describing orthologous and paralogous relationships while allowing the user to store custom and database-specific data in the same format. Although many orthology databases use it as a way to export data, there is no comprehensive toolkit for working with this format. Here, we introduce the OrthoXML-tools ( https://github.com/DessimozLab/orthoxml-tools ), a comprehensive toolkit for loading, manipulating, and exporting the OrthoXML files to other formats. We show its capabilities and performance on our benchmarks.

The human natural anti-carbohydrate antibodies anti-Gal, anti-Neu5Gc, and anti-Forssman are "living-fossils" that appeared in ancestral apes, monkeys and hominins millions of years ago. These antibodies appeared at various evolutionary periods in few mutated-offspring that lost the ability to synthesize the corresponding dispensable (i.e., nonessential) carbohydrate-antigens, α-gal epitope, Neu5Gc (N-glycolyl neuraminic acid) and Forssman-antigen, respectively. Production of these antibodies is stimulated by environmental antigens such as those of the human microbiota. Elimination of carbohydrate-antigens in the few mutated-offspring was caused by accidental nonsense or missense mutations that inactivated genes encoding enzymes involved in their biosynthesis, while mostindividuals in parental-populations continued synthesizing these carbohydrate-antigens. It has been suggested that dispensable carbohydrate-antigens which are absent in some mammalian species were evolutionary eliminated due to selective pressure by lethal viruses using these carbohydrate-antigens as "docking" receptors. An alternative selective mechanism which is based on the distribution of anti-Gal, anti-Neu5Gc and anti-Forssman in mammals, is presented in this review and is associated with the protective effects of these natural antibodies. It is suggested that epidemics of lethal enveloped-viruses caused the extinction of parental-populations synthesizing the corresponding carbohydrate-antigens of these antibodies, independent of the cell adhesion mechanisms of such viruses. However, the few mutated offspring were protected by these natural antibodies which bound to carbohydrate-antigens synthesized on viruses as a result of their replication in individuals of the parental-populations. Recent studies suggest that these antibodies continue to contribute to the immune protection of humans against zoonotic infections by viruses presenting α-gal, Neu5Gc or Forssman antigens.

Monocyte chemoattractant proteins (MCPs) are involved in monocyte trafficking during severe inflammation by modulating the chemokine-glycosaminoglycan-receptor signaling axis. MCPs comprise a family of four chemokines (CCL2, CCL7, CCL8, and CCL13/12) that exhibit differential expression patterns in mammals, functional diversity, and receptor/glycosaminoglycan (GAG) binding promiscuity. In this context, the evolution-structure-function paradigm of MCP chemokines in mammals was established by assessing phylogeny, functional divergence, selection pressure, and coevolution in correlation with structural and surface characteristics. Comprehensive analyses were performed using an array of evolutionary and structural bioinformatic methods including molecular dynamics simulations. Our findings demonstrate that substitutions in receptor/GAG-interacting residues mediate episodic diversification and functional diversity in MCP chemokines. Additionally, a balanced interplay of selection pressures has driven the functional changes observed among MCP paralogs, with positive selection at various receptor/GAG-binding sites contributing to their promiscuous receptor/GAG interactions. Meanwhile, processes like purifying selection and coevolution maintain the classical chemokine structure and preserve the ancestral functions of MCP chemokines. Overall, this study suggests that selection pressure on sites within the N-terminal region [N-loop and 310-helix] and 40S loop of MCP chemokines alters surface properties to fine-tune the molecular interactions and functional characteristics without altering the overall chemokine structure.

One of the most important and difficult challenges in the research of molecular evolution is modeling the process of amino acid substitutions. Although single-matrix models, such as the LG model, are popular, their capability to properly capture the heterogeneity of the substitution process across sites is still questioned. Several mixture models with multiple matrices have been introduced and shown to offer advantages over single-matrix models. Current general mixture models assume the reversibility of the evolutionary process, implying that substitution rates between any two amino acids are equal in both forward and backward directions. This assumption is not based on biological properties but rather on computational simplicity. The well-known hypothesis is that more realistic models can yield more accurate evolutionary inferences; therefore, our aim is to estimate more biologically realistic models. To this end, we relax the assumption of reversibility and introduce two new general non-reversible 4-matrix mixture models, called nT4M and nT4X. Using alignments from HSSP and TreeBASE databases as data, our newly estimated models outperformed all single-matrix models and almost all reversible mixture models. Moreover, the new non-reversible mixture models enable us to infer rooted trees.

The major histocompatibility complex (MHC) is a cluster of functionally related genes encoding proteins which, among other functions, mediate immune system activation. While the MHC of many vertebrates has been extensively studied, less is known about the amphibian MHC. This represents an important knowledge gap because amphibians mark the evolutionary transition from an aquatic to a terrestrial lifestyle and often maintain a biphasic lifestyle. Hence, they tend to be exposed to both aquatic and terrestrial pathogen communities, providing opportunities to gain fundamental insights into how the immune system responds to different environmental challenges. Moreover, amphibians are globally threatened by invasive pathogens and the MHC may play a role in combating population decline. In this review, we summarize the current state of knowledge regarding the amphibian MHC and identify the major differences with other vertebrates. We also review how the number of MHC gene copies varies across amphibian groups and how MHC-based variation relates to amphibian ontogeny, behaviour, disease, and phylogeography. We conclude by identifying knowledge gaps and proposing priorities for future research.

Major evolutionary transitions, such as the shift of cetaceans from terrestrial to marine life, can put pressure on sensory systems to adapt to a new set of relevant stimuli. Relatively little is known about the role of smell in the evolution of mysticetes (baleen whales). While their toothed cousins, the odontocetes, lack the anatomical features to smell, it is less clear whether baleen whales have retained this sense, and if so, when the pressure on olfaction diverged in the cetacean evolutionary lineage. We examined eight genes encoding olfactory signal transduction pathway components and key chaperones for signs of inactivating mutations and selective pressures. All of the genes we examined were intact in all eight mysticete genomes examined, despite inactivating mutations in odontocete homologs in multiple genes. We also tested several models representing various hypotheses regarding the evolutionary history of olfaction in cetaceans. Our results support a model where olfactory ability is specifically reduced in the odontocete lineage following their split from stem cetaceans and serve to clarify the evolutionary history of olfaction in cetaceans.

The highly dynamic nature of the Cotton leaf curl virus (CLCuV) complex (causing Cotton leaf curl disease, a significant global threat to cotton) presents a formidable challenge in unraveling precise molecular mechanisms governing viral-host interactions. To address this challenge, the present study investigated the molecular interactions of 6 viral proteins (Rep, TrAP, C4, C5, V2, and βC1) with 18 cotton Transcriptional Gene Silencing (TGS) proteins. Protein-protein dockings conducted for different viral-host protein pairs using Clustered Protein Docking (ClusPro) and Global RAnge Molecular Matching (GRAMM) (216 docking runs), revealed variable binding energies. The interacting pairs with the highest binding affinities were further scrutinized using bioCOmplexes COntact MAPS (COCOMAPS) server, which revealed robust binding of three viral proteins- TrAP, C4, and C5 with 14 TGS proteins, identifying several novel interactions (not reported yet by earlier studies), such as TrAP targeting DCL3, HDA6, and SUVH6; C4 targeting RAV2, CMT2, and DMT1; and C5 targeting CLSY1, RDR1, RDR2, AGO4, SAMS, and SAHH. Visualizing these interactions in PyMol provided a detailed insight into interacting regions. Further assessment of the impact of 18 variants of the C4 protein on interaction with CMT2 revealed no correlation between sequence variation and docking energies. However, conserved residues in the C4 binding regions emerged as potential targets for disrupting viral integrity. Hence, this study provides valuable insights into the viral-host interplay, advancing our understanding of Cotton leaf curl Multan virus pathogenicity and opening novel avenues for devising various antiviral strategies by targeting the host-viral interacting regions after experimental validation.