Genomic Analyses Identify Novel Molecular Signatures Specific for the Caenorhabditis and other Nematode Taxa Providing Novel Means for Genetic and Biochemical Studies

The phylum Nematoda encompasses numerous free-living as well as parasitic members, including the widely used animal model Caenorhabditis elegans, with significant impact on human health, agriculture, and environment. In view of the importance of nematodes, it is of much interest to identify novel molecular characteristics that are distinctive features of this phylum, or specific taxonomic groups/clades within it, thereby providing innovative means for diagnostics as well as genetic and biochemical studies. Using genome sequences for 52 available nematodes, a robust phylogenetic tree was constructed based on concatenated sequences of 17 conserved proteins. The branching of species in this tree provides important insights into the evolutionary relationships among the studied nematode species. In parallel, detailed comparative analyses on protein sequences from nematodes (Caenorhabditis) species reported here have identified 52 novel molecular signatures (or synapomorphies) consisting of conserved signature indels (CSIs) in different proteins, which are uniquely shared by the homologs from either all genome-sequenced Caenorhabditis species or a number of higher taxonomic clades of nematodes encompassing this genus. Of these molecular signatures, 39 CSIs in proteins involved in diverse functions are uniquely present in all Caenorhabditis species providing reliable means for distinguishing this group of nematodes in molecular terms. The remainder of the CSIs are specific for a number of higher clades of nematodes and offer important insights into the evolutionary relationships among these species. The structural locations of some of the nematodes-specific CSIs were also mapped in the structural models of the corresponding proteins. All of the studied CSIs are localized within the surface-exposed loops of the proteins suggesting that they may potentially be involved in mediating novel protein–protein or protein–ligand interactions, which are specific for these groups of nematodes. The identified CSIs, due to their exclusivity for the indicated groups, provide reliable means for the identification of species within these nematodes groups in molecular terms. Further, due to the predicted roles of these CSIs in cellular functions, they provide important tools for genetic and biochemical studies in Caenorhabditis and other nematodes

Novel Molecular Signatures in the PIP4K/PIP5K Family of Proteins Specific for Different Isozymes and Subfamilies Provide Important Insights into the Evolutionary Divergence of this Protein Family

Members of the PIP4K/PIP5K family of proteins, which generate the highly important secondary messenger phosphatidylinositol-4,5-bisphosphate, play central roles in regulating diverse signaling pathways. In eukaryotic organisms, multiple isozymes and subfamilies of PIP4K/PIP5K proteins are found and it is of much interest to understand their evolution and species distribution and what unique molecular and biochemical characteristics distinguish specific isozymes and subfamilies of proteins. We report here the species distribution of different PIP4K/PIP5K family of proteins in eukaryotic organisms and phylogenetic analysis based on their protein sequences. Our results indicate that the distinct homologs of both PIP4K and PIP5K are found in different organisms belonging to the Holozoa clade of eukaryotes, which comprises of various metazoan phyla as well as their close unicellular relatives Choanoflagellates and Filasterea. In contrast, the deeper-branching eukaryotic lineages, as well as plants and fungi, contain only a single homolog of the PIP4K/PIP5K proteins. In parallel, our comparative analyses of PIP4K/PIP5K protein sequences have identified six highly-specific molecular markers consisting of conserved signature indels (CSIs) that are uniquely shared by either the PIP4K or PIP5K proteins, or both, or specific subfamilies of these proteins. Of these molecular markers, 2 CSIs are distinctive characteristics of all PIP4K homologs, 1 CSI distinguishes the PIP4K and PIP5K homologs from the Holozoa clade of species from the ancestral form of PIP4K/PIP5K found in deeper-branching eukaryotic lineages. The remaining three CSIs are specific for the PIP5Kα, PIP5Kβ, and PIP4Kγ subfamilies of proteins from vertebrate species. These molecular markers provide important means for distinguishing different PIP4K/PIP5K isozymes as well as some of their subfamilies. In addition, the distribution patterns of these markers in different isozymes provide important insights into the evolutionary divergence of PIP4K/PIP5K proteins. Our results support the view that the Holozoa clade of eukaryotic organisms shared a common ancestor exclusive of the other eukaryotic lineages and that the initial gene duplication event leading to the divergence of distinct types of PIP4K and PIP5K homologs occurred in a common ancestor of this clade. Based on the results gleaned from different studies presented here, a model for the evolutionary divergence of the PIP4K/PIP5K family of proteins is presented

Conserved molecular signatures in the spike protein provide evidence indicating the origin of SARS-CoV-2 and a Pangolin-CoV (MP789) by recombination(s) between specific lineages of Sarbecoviruses

Both SARS-CoV-2 and SARS coronaviruses (CoVs) are members of the subgenus Sarbecovirus. To understand the origin of SARS-CoV-2, sequences for the spike and nucleocapsid proteins from sarbecoviruses were analyzed to identify molecular markers consisting of conserved inserts or deletions (termed CSIs) that are specific for either a particular clade of Sarbecovirus or are commonly shared by two or more clades of these viruses. Three novel CSIs in the N-terminal domain (NTD) of the spike protein S1-subunit (S1-NTD) are uniquely shared by SARS-CoV-2, Bat-CoV-RaTG13 and most pangolin CoVs (SARS-CoV-2r clade). Three other sarbecoviruses viz. bat-CoVZXC21, -CoVZC45 and -PrC31 (forming CoVZC/PrC31 clade), and a pangolin-CoV_MP789 also contain related CSIs in the same positions. In contrast to the S1-NTD, both SARS and SARS-CoV-2r viruses contain two large CSIs in the S1-C-terminal domain (S1-CTD) that are absent in the CoVZC/PrC31 clade. One of these CSIs, consisting of a 12 aa insert, is also present in the RShSTT clade (Cambodia-CoV strains). Sequence similarity studies show that the S1-NTD of SARS-CoV-2r viruses is most similar to the CoVZC/PrC31 clade, whereas their S1-CTD exhibits highest similarity to the RShSTT- (and the SARS-related) CoVs. Results from the shared presence of CSIs and sequence similarity studies on different CoV lineages support the inference that the SARS-CoV-2r cluster of viruses has originated by a genetic recombination between the S1-NTD of the CoVZC/PrC31 clade of CoVs and the S1-CTD of RShSTT/SARS viruses, respectively. We also present compelling evidence, based on the shared presence of CSIs and sequence similarity studies, that the pangolin-CoV_MP789, whose receptor-binding domain is most similar to the SARS-CoV-2 virus, has resulted from another independent recombination event involving the S1-NTD of the CoVZC/PrC31 CoVs and the S1-CTD of an unidentified SARS-CoV-2r related virus. The SARS-CoV-2 virus involved in this latter recombination event is postulated to be most similar to the SARS-CoV-2. Several other CSIs reported here are specific for other clusters of sarbecoviruses including a clade consisting of bat-SARS-CoVs (BM48-31/BGR/2008 and SARS_BtKY72). Structural mapping studies show that the identified CSIs form distinct loops/patches on the surface of the spike protein. It is hypothesized that these novel loops/patches on the spike protein, through their interactions with other host components, should play important roles in the biology/pathology of SARS-CoV-2 virus. Lastly, the CSIs specific for different clades of sarbecoviruses including SARS-CoV-2r clade provide novel means for the identification of these viruses and other potential applications

Conserved Molecular Signatures in the Spike, Nucleocapsid, and Polymerase Proteins Specific for the Genus Betacoronavirus and Its Different Subgenera

The genus Betacoronavirus, consisting of four main subgenera (Embecovirus, Merbecovirus, Nobecovirus, and Sarbecovirus), encompasses all clinically significant coronaviruses (CoVs), including SARS, MERS, and the SARS-CoV-2 virus responsible for current COVID-19 pandemic. Very few molecular characteristics are known that are specific for the genus Betacoronavirus or its different subgenera. In this study, our analyses of the sequences of four essential proteins of CoVs, viz., spike, nucleocapsid, envelope, and RNA-dependent RNA polymerase (RdRp), identified ten novel molecular signatures consisting of conserved signature indels (CSIs) in these proteins which are specific for the genus Betacoronavirus or its subgenera. Of these CSIs, two 14-aa-conserved deletions found within the heptad repeat motifs 1 and 2 of the spike protein are specific for all betacoronaviruses, except for their shared presence in the highly infectious avian coronavirus. Six additional CSIs present in the nucleocapsid protein and one CSI in the RdRp protein are distinctive characteristics of either the Merbecovirus, Nobecovirus, or Sarbecovirus subgenera. In addition, a 4-aa insert is present in the spike protein, which is uniquely shared by all viruses from the subgenera Merbecovirus, Nobecovirus, and Sarbecovirus, but absent in Embecovirus and all other genera of CoVs. This molecular signature provides evidence that viruses from the three subgenera sharing this CSI are more closely related to each other, and they evolved after the divergence of embecoviruses and other CoVs. As all CSIs specific for different groups of CoVs are flanked by conserved regions, their sequences provide novel means for identifying the above groups of CoVs and for developing novel diagnostic tests. Furthermore, our analyses of the structures of the spike and nucleocapsid proteins show that all identified CSIs are localized in the surface-exposed loops of these protein. It is postulated that these surface loops, through their interactions with other cellular proteins/ligands, play important roles in the biology/pathology of these viruses

Novel molecular, structural and evolutionary characteristics of the phosphoketolases from bifidobacteria and Coriobacteriales.

Members from the order Bifidobacteriales, which include many species exhibiting health promoting effects, differ from all other organisms in using a unique pathway for carbohydrate metabolism, known as the "bifid shunt", which utilizes the enzyme phosphoketolase (PK) to carry out the phosphorolysis of both fructose-6-phosphate (F6P) and xylulose-5-phosphate (X5P). In contrast to bifidobacteria, the PKs found in other organisms (referred to XPK) are able to metabolize primarily X5P and show very little activity towards F6P. Presently, very little is known about the molecular or biochemical basis of the differences in the two forms of PKs. Comparative analyses of PK sequences from different organisms reported here have identified multiple high-specific sequence features in the forms of conserved signature inserts and deletions (CSIs) in the PK sequences that clearly distinguish the X5P/F6P phosphoketolases (XFPK) of bifidobacteria from the XPK homologs found in most other organisms. Interestingly, most of the molecular signatures that are specific for the XFPK from bifidobacteria are also shared by the PK homologs from the Coriobacteriales order of Actinobacteria. Similarly to the Bifidobacteriales, the order Coriobacteriales is also made up of commensal organisms, that are saccharolytic and able to metabolize wide variety of carbohydrates, producing lactate and other metabolites. Phylogenetic studies provide evidence that the XFPK from bifidobacteria are specifically related to those found in the Coriobacteriales and suggest that the gene for PK (XFPK) was horizontally transferred between these two groups. A number of the identified CSIs in the XFPK sequence, which serve to distinguish the XFPK homologs from XPK homologs, are located at the subunit interface in the structure of the XFPK dimer protein. The results of protein modelling and subunit docking studies indicate that these CSIs are involved in the formation/stabilization of the protein dimer. The significance of these observations regarding the differences in the activities of the XFPK and XPK homologs are discussed. Additionally, this work also discusses the significance of the XFPK-like homologs, similar to those found in bifidobacteria, in the order Coriobacteriales

Ribonucleotide Reductases from Bifidobacteria Contain Multiple Conserved Indels Distinguishing Them from All Other Organisms: In Silico Analysis of the Possible Role of a 43 aa Bifidobacteria-Specific Insert in the Class III RNR Homolog

Bifidobacteria comprises an important group/order of bacteria whose members have widespread usage in the food and health industry due to their health-promoting activity in the human gastrointestinal tract. However, little is known about the underlying molecular properties that are responsible for the probiotic effects of these bacteria. The enzyme ribonucleotide reductase (RNR) plays a key role in all organisms by reducing nucleoside di- or tri- phosphates into corresponding deoxyribose derivatives required for DNA synthesis, and RNR homologs belonging to classes I and III are present in either most or all Bifidobacteriales. Comparative analyses of these RNR homologs have identified several novel sequence features in the forms of conserved signature indels (CSIs) that are exclusively found in bifidobacterial RNRs. Specifically, in the large subunit of the aerobic class Ib RNR, three CSIs have been identified that are uniquely found in the Bifidobacteriales homologs. Similarly, the large subunit of the anaerobic class III RNR contains five CSIs that are also distinctive characteristics of bifidobacteria. Phylogenetic analyses indicate that these CSIs were introduced in a common ancestor of the Bifidobacteriales and retained by all descendants, likely due to their conferring advantageous functional roles. The identified CSIs in the bifidobacterial RNR homologs provide useful tools for further exploration of the novel functional aspects of these important enzymes that are exclusive to these bacteria. We also report here the results of homology modeling studies, which indicate that most of the bifidobacteria-specific CSIs are located within the surface loops of the RNRs, and of these, a large 43 amino acid insert in the class III RNR homolog forms an extension of the allosteric regulatory site known to be essential for protein function. Preliminary docking studies suggest that this large CSI may be playing a role in enhancing the stability of the RNR dimer complex. The possible significances of the identified CSIs, as well as the distribution of RNR homologs in the Bifidobacteriales, are discussed

Primary sequence of phosphoketolase protein from <i>Bifidobacterium longum</i> depicting the location of different CSIs.

<p>Secondary structure elements are displayed above the primary sequence, with helices shown as cylinders and sheets shown as arrows. The secondary structure information was obtained directly from the solved structure of <i>B</i>. <i>longum</i> PK from Protein Data Bank (PDB ID: 3AI7). NCBI CD-search webserver was used to demarcate the three domains (yellow boxes): N-terminal PP-domain (PP-D), middle PYR-domain (PYR-D), and the C-terminal domain (CT-D). The eight identified CSIs in bifidobacteria PKs are indicated by bold and underlined residues, insertions are shown in red and the positions where deletions are present are shown in blue. The arrows on the top of CSIs indicate the residues contributing in subunit-subunit interactions, as determined via the PISA webserver.</p

Novel molecular, structural and evolutionary characteristics of the phosphoketolases from bifidobacteria and <i>Coriobacteriales</i>

<div><p>Members from the order <i>Bifidobacteriales</i>, which include many species exhibiting health promoting effects, differ from all other organisms in using a unique pathway for carbohydrate metabolism, known as the “bifid shunt”, which utilizes the enzyme phosphoketolase (PK) to carry out the phosphorolysis of both fructose-6-phosphate (F6P) and xylulose-5-phosphate (X5P). In contrast to bifidobacteria, the PKs found in other organisms (referred to XPK) are able to metabolize primarily X5P and show very little activity towards F6P. Presently, very little is known about the molecular or biochemical basis of the differences in the two forms of PKs. Comparative analyses of PK sequences from different organisms reported here have identified multiple high-specific sequence features in the forms of conserved signature inserts and deletions (CSIs) in the PK sequences that clearly distinguish the X5P/F6P phosphoketolases (XFPK) of bifidobacteria from the XPK homologs found in most other organisms. Interestingly, most of the molecular signatures that are specific for the XFPK from bifidobacteria are also shared by the PK homologs from the <i>Coriobacteriales</i> order of Actinobacteria. Similarly to the <i>Bifidobacteriales</i>, the order <i>Coriobacteriales</i> is also made up of commensal organisms, that are saccharolytic and able to metabolize wide variety of carbohydrates, producing lactate and other metabolites. Phylogenetic studies provide evidence that the XFPK from bifidobacteria are specifically related to those found in the <i>Coriobacteriales</i> and suggest that the gene for PK (XFPK) was horizontally transferred between these two groups. A number of the identified CSIs in the XFPK sequence, which serve to distinguish the XFPK homologs from XPK homologs, are located at the subunit interface in the structure of the XFPK dimer protein. The results of protein modelling and subunit docking studies indicate that these CSIs are involved in the formation/stabilization of the protein dimer. The significance of these observations regarding the differences in the activities of the XFPK and XPK homologs are discussed. Additionally, this work also discusses the significance of the XFPK-like homologs, similar to those found in bifidobacteria, in the order <i>Coriobacteriales</i>.</p></div

A maximum likelihood distance tree based on PKs sequences for members of the phylum Actinobacteria are representative outgroup species from the phylum <i>Firmicutes</i>.

<p>The numbers on the nodes indicate bootstrap scores for the group of species represented by different nodes. In this tree, members of the order <i>Bifidobacteriales</i> branch with the <i>Coriobacteriales</i> and the clade comprising of these two orders is separated from all other Actinobacteria/bacteria by a long branch. Based on the species distributions of different CSIs, the evolutionary stages where genetic changes giving rise to different CSIs have likely occurred are marked.</p