Comparative Protein Sequence Analysis Research Articles

Prospects for a sequence-based taxonomy of influenza A virus subtypes.

Hemagglutinin (HA) and neuraminidase (NA) proteins are the primary antigenic targets of influenza A virus (IAV) infections. IAV infections are generally classified into subtypes of HA and NA proteins, e.g. H3N2. Most of the known subtypes were originally defined by a lack of antibody cross-reactivity. However, genetic sequencing has played an increasingly important role in characterizing the evolving diversity of IAV. Novel subtypes have recently been described solely by their genetic sequences, and IAV infections are routinely subtyped by molecular assays, or the comparison of sequences to references. In this study, I carry out a comparative analysis of all available IAV protein sequences in the Genbank database (over 1.1 million, reduced to 272,292 unique sequences prior to phylogenetic reconstruction) to determine whether the serologically defined subtypes can be reproduced with sequence-based criteria. I show that a robust genetic taxonomy of HA and NA subtypes can be obtained using a simple clustering method, namely, by progressively partitioning the phylogeny on its longest internal branches. However, this taxonomy also requires some amendments to the current nomenclature. For example, two IAV isolates from bats previously characterized as a divergent lineage of H9N2 should be separated into their own subtype. With the exception of these small and highly divergent lineages, the phylogenies relating each of the other six genomic segments do not support partitions into major subtypes.

Chromosome-level genome assembly of Quercus variabilis provides insights into the molecular mechanism of cork thickness

Quercus variabilis is a deciduous woody species with high ecological and economic value, and is a major source of cork in East Asia. Cork from thick softwood sheets have higher commercial value than those from thin sheets. It is extremely difficult to genetically improve Q. variabilis to produce high quality softwood due to the lack of genomic information. Here, we present a high-quality chromosomal genome assembly for Q. variabilis with length of 791,89 Mb and 54,606 predicted genes. Comparative analysis of protein sequences of Q. variabilis with 11 other species revealed that specific and expanded gene families were significantly enriched in the "fatty acid biosynthesis" pathway in Q. variabilis, which may contribute to the formation of its unique cork. Based on weighted correlation network analysis of time-course (i.e., five important developmental ages) gene expression data in thick-cork versus thin-cork genotypes of Q. variabilis, we identified one co-expression gene module associated with the thick-cork trait. Within this co-expression gene module, 10 hub genes were associated with suberin biosynthesis. Furthermore, we identified a total of 198 suberin biosynthesis-related new candidate genes that were up-regulated in trees with a thick cork layer relative to those with a thin cork layer. Also, we found that some genes related to cell expansion and cell division were highly expressed in trees with a thick cork layer. Collectively, our results revealed that two metabolic pathways (i.e., suberin biosynthesis, fatty acid biosynthesis), along with other genes involved in cell expansion, cell division, and transcriptional regulation, were associated with the thick-cork trait in Q. variabilis, providing insights into the molecular basis of cork development and knowledge for informing genetic improvement of cork thickness in Q. variabilis and closely related species.

Comparative analysis of SilA-laccase mediated degradation of ciprofloxacin, norfloxacin and ofloxacin and interpretation of the possible catalytic mechanism

Fluoroquinolones (FQs) are the most commonly used antimicrobial drugs and regardless of their advantages in the healthcare sector, the pollution of these antimicrobial drugs in the environment has big concerns about human and environmental health. The presence of these antibiotic drugs even at the lowest concentrations in the environment has resulted in the emergence and spread of antibiotic resistance. Hence, it is necessary to remediate these pollutants from the environment. Previously alkaline laccase (SilA) from Streptomyces ipomoeae has been demonstrated to show degrading potentials against two of the FQs, Ciprofloxacin (CIP) and Norfloxacin (NOR); however, the molecular mechanism was not elucidated in detail. In this study, we have analyzed the possible molecular catalytic mechanism of FQ degrading SilA-laccase for the degradation of the FQs, CIP, NOR and Ofloxacin (OFL) using three-dimensional protein structure modeling, molecular docking and molecular dynamic (MD) studies. The comparative protein sequence analysis revealed the presence of tetrapeptide conserved catalytic motif, His102-X-His104-Gly105. After evaluating the active site of the enzyme in depth using CDD, COACH and S-site tools, we have identified the catalytic triad composed of three conserved amino acid residues, His102, Val103 and Tyr108 with which ligands interacted during the catalysis process. By analyzing the MD trajectories, it is revealed that the highest degradation potential of SilA is for CIP followed by NOR and OFL. Ultimately, this study provides the possible comparative catalytic mechanism for the degradation of CIP, NOR and OFL by the SilA enzyme. Communicated by Ramaswamy H. Sarma

Cloning, homology modelling and expression analysis of Oryza sativa WNK gene family

With No Lysine kinases (WNKs) represents a gene family that encodes Ser/Thr kinases, with anomalous disposition of catalytic lysine residue in subdomain I. In plants, WNKs had been linked to circadian rhythm, photoperiodic response and abiotic stress tolerance with mechanism yet undeciphered. In the present study, full-length CDS sequences of rice WNKs (OsWNK1 to 8) were cloned from indica cultivar IR64. A total of six highly conserved kinase subdomains were identified. Comparative analysis of protein sequences from six different species of rice showed varying magnitudes of substitution (76.2 %), deletion (15.4 %), and addition (8.4 %) events. ConSurf analysis coupled with CASTp results identified functional residues that were clustered together in modelled 3-D structures. Among post-translational modifications (PTMs) studied, 87.7 % of phosphorylation sites were predicted. Mined protein-protein interactions (PPIs) depicted OsWNKs to interact notably with other OsWNK members and with key proteins like PRR95 involved in photoperiodic response and protein phosphatase like PP2C involved in ABA signalling. Gene duplication analysis revealed two paralogous duplicated gene pairs: WNK6-WNK9 and WNK7-WNK8. Oryza sativa showed maximum syntenic relationship with Sorghum bicolor among the compared species. OsWNKs showed differential transcript expression profiles on treatment with plant growth regulators indicating its versatile role in plant growth and development.

Collagen Pentablock Copolymers Form Smectic Liquid Crystals as Precursors for Mussel Byssus Fabrication.

Protein-based biological materials are important role models for the design and fabrication of next generation advanced polymers. Marine mussels (Mytilus spp.) fabricate hierarchically structured collagenous fibers known as byssal threads via bottom-up supramolecular assembly of fluid protein precursors. The high degree of structural organization in byssal threads is intimately linked to their exceptional toughness and self-healing capacity. Here, we investigated the hypothesis that multidomain collagen precursor proteins, known as preCols, are stored in secretory vesicles as a colloidal liquid crystal (LC) phase prior to thread self-assembly. Using advanced electron microscopy methods, including scanning TEM and FIB-SEM, we visualized the detailed smectic preCol LC nanostructure in 3D, including various LC defects, confirming this hypothesis and providing quantitative insights into the mesophase structure. In light of these findings, we performed an in-depth comparative analysis of preCol protein sequences from multiple Mytilid species revealing that the smectic organization arises from an evolutionarily conserved ABCBA pentablock copolymer-like primary structure based on demarcations in hydropathy and charge distribution as well as terminal pH-responsive domains that trigger fiber formation. These distilled supramolecular assembly principles provide inspiration and strategies for sustainable assembly of nanostructured polymeric materials for potential applications in engineering and biomedical applications.

A novel β-glucosidase from a hot-spring metagenome shows elevated thermal stability and tolerance to glucose and ethanol

β-glucosidase causes hydrolysis of β-1,4-glycosidic bond in glycosides and oligosaccharides. It is an industrially important enzyme owing to its potential in biomass processing applications. In this study, computational screening of an extreme temperature aquatic habitat metagenomic resource was done, leading to the identification of a novel gene, bglM, encoding a β-glucosidase. The comparative protein sequence and homology structure analyses designated it as a GH1 family β-glucosidase. The bglM gene was expressed in a heterologous host, Escherichia coli. The purified protein, BglM, was biochemically characterized for β-glucosidase activity. BglM exhibited noteworthy hydrolytic potential towards cellobiose and lactose. BglM, showed substantial catalytic activity in the pH range of 5.0–7.0 and at the temperature 40 °C–70 °C. The enzyme was found quite stable at 50 °C with a loss of hardly 20% after 40 h of heat exposure. Furthermore, any drastically negative effect was not observed on the enzyme’s activity in the presence of metal ions, non-ionic surfactants, metal chelating, and denaturing agents. A significantly high glucose tolerance, retaining 80% relative activity at 1 M, and 40% at 5 M glucose, and ethanol tolerance, exhibiting 80% relative activity in 10% ethanol, enrolled BglM as a promising enzyme for cellulose saccharification. Furthermore, its ability to catalyze the hydrolysis of daidzin and polydatin ascertained it as an admirably suited biocatalyst for enhancement of nutritional values in soya and wine industries.

SARS-CoV-2 attachment to host cells is possibly mediated via RGD-integrin interaction in a calcium-dependent manner and suggests pulmonary EDTA chelation therapy as a novel treatment for COVID 19

SARS-CoV-2 is a highly contagious virus that has caused serious health crisis world-wide resulting into a pandemic situation. As per the literature, the SARS-CoV-2 is known to exploit humanACE2 receptors (similar toprevious SARS-CoV-1) for gaining entry into the host cell for invasion, infection, multiplication and pathogenesis. However, considering the higher infectivity of SARS-CoV-2 along with the complex etiology and pathophysiological outcomes seen in COVID-19 patients, it seems that there may be an alternate receptor for SARS-CoV-2. I performed comparative protein sequence analysis, database based gene expression profiling, bioinformatics based molecular docking using authentic tools and techniques for unveiling the molecular basis of high infectivity of SARS-CoV-2 as compared to previous known coronaviruses. My study revealed that SARS-CoV-2 (previously known as 2019-nCoV) harbors a RGD motif in its receptor binding domain (RBD) and the motif is absent in all other previously known SARS-CoVs. The RGD motif is well known for its role in cell-attachment and cell-adhesion. My hypothesis is that the SARS-CoV-2 may be (via RGD) exploiting integrins, that have high expression in lungs and all other vital organs, for invading host cells. However, an experimental verification is required. The expression of ACE2, which is a known receptor for SARS-CoV-2, was found to be negligible in lungs. I assume that higher infectivity of SARS-CoV-2 could be due to this RGD-integrin mediated acquired cell-adhesive property. Gene expression profiling revealed that expression of integrins is significantly high in lung cells, in particular αvβ6, α5β1, αvβ8 and an ECM protein, ICAM1. The molecular docking experiment showed the RBD of spike protein binds with integrins precisely at RGD motif in a similar manner as a synthetic RGD peptide binds to integrins as found by other researchers. SARS-CoV-2 spike protein has a number of phosphorylation sites that can induce cAMP, PKC, Tyr signaling pathways. These pathways either activate calcium ion channels or get activated by calcium. In fact, integrins have calcium & metal binding sites that were predicted around and in vicinity of RGD-integrin docking site in our analysis which suggests that RGD-integrins interaction possibly occurs in calcium-dependent manner. The higher expression of integrins in lungs along with their previously known high binding affinity (~KD = 4.0 nM) for virus RGD motif could serve as a possible explanation for high infectivity of SARS-CoV-2. On the contrary, human ACE2 has lower expression in lungs and its high binding affinity (~KD = 15 nM) for spike RBD alone could not manifest significant virus-host attachment. This suggests that besides human ACE2, an additional or alternate receptor for SARS-CoV-2 is likely to exist. A highly relevant evidence never reported earlier which corroborate in favor of RGD-integrins mediated virus-host attachment is an unleashed cytokine storm which causes due to activation of TNF-α and IL-6 activation; and integrins role in their activation is also well established. Altogether, the current study has highlighted possible role of calcium and other divalent ions in RGD-integrins interaction for virus invasion into host cells and suggested that lowering divalent ion in lungs could avert virus-host cells attachment.

Structure-Based Virtual Screening to Discover Potential Lead Molecules for the SARS-CoV-2 Main Protease.

The COVID-19 disease is caused by a new strain of the coronavirus family (SARS-CoV-2), and it has affected at present millions of people all over the world. The indispensable role of the main protease (Mpro) in viral replication and gene expression makes this enzyme an attractive drug target. Therefore, inhibition of SARS-CoV-2 Mpro as a proposition to halt virus ingression is being pursued by scientists globally. Here we carried out a study with two objectives: the first being to perform comparative protein sequence and 3D structural analysis to understand the effect of 12 point mutations on the active site. Among these, two mutations, viz., Ser46 and Phe134, were found to cause a significant change at the active sites of SARS-CoV-2. The Ser46 mutation present at the entrance of the S5 subpocket of SARS-CoV-2 increases the contribution of other two hydrophilic residues, while the Phe134 mutation, present in the catalytic cysteine loop, can cause an increase in catalytic efficiency of Mpro by facilitating fast proton transfer from the Cys145 to His41 residue. It was observed that active site remained conserved among Mpro of both SARS-CoVs, except at the entrance of the S5 subpocket, suggesting sustenance of substrate specificity. The second objective was to screen the inhibitory effects of three different data sets (natural products, coronaviruses main protease inhibitors, and FDA-approved drugs) using a structure-based virtual screening approach. A total of 73 hits had a combo score >2.0. Eight different structural scaffold classes were identified, such as one/two tetrahydropyran ring(s), dipeptide/tripeptide/oligopeptide, large (approximately 20 atoms) cyclic peptide, and miscellaneous. The screened hits showed key interactions with subpockets of the active site. Further, molecular dynamics studies of selected screened compounds confirmed their perfect fitting into the subpockets of the active site. This study suggests promising structures that can fit into the SARS-CoV-2 Mpro active site and also offers direction for further lead optimization and rational drug design.

Phylogenetic analysis and molecular signatures specific to the Ralstonia solanacearum species complex

Ralstonia solanacearum species complex is a highly destructive bacterial plant pathogen recognized as a group of genotypically-diverse strains with broad geographic distribution and wide host range. Strains are divided into four phylotypes (I, II, III and IV) relating to their geographical origins. Based on the polyphasic approach, phylotypes I and III have been reclassified as R. pseudosolanacearum, phylotype IV as R. syzygii species and phylotype II as R. solanacearum. The large data sets of whole genome sequences for the members of R. solanacearum allow us to discover novel molecular markers to evaluate the phylogenetic relationships within this species. In the phylogenetic tree based on concatenated amino acid sequences of 57 housekeeping genes, members of R. solanacearum species complex formed strongly supported clades that clearly distinguished as R. solanacearum, R. pseudosolanacearum and R. syzygii species. In this study, the conserved signature indels (CSIs) which are specifically present in protein sequences and can be used to differentiate various groups of Ralstonia species are reported. Comparative analysis of whole protein sequences from the three species has identified seven CSIs that are specific to R. solanacearum (phylotype II), 6 CSIs specific to R. pseudosolanacearum strains (phylotype I and III) and five CSIs that are specifically shared by members of R. syzygii (phylotype IV). These CSIs are flanked by conserved sequence regions, providing highly specific tools for identification of previously known or novel strains of these species.

Zebrafish ubiquitin-specific peptidase 5 (USP5) activates interferon resistance to the virus by increase the expression of RIG-I.

The ubiquitin specific peptidase (USP) family is involved in many life processes, of which antiviral is also an important basic function. One of the more important ways is to activate interferon. In this study, we reported the antiviral function of the ubiquitin specific peptidase 5(USP5) gene in zebrafish. Evolutionary and comparative protein sequence analysis of the USP5 was performed. The localization of USP5 in FHM cells cytoplasm was determined. Overexpression of USP5 significantly evoked higher expression of mRNA that encode IFNφ1 and ISGs, the promoteractivities of IFNφ1 and IFNstimulated response element (ISRE) were augmented likewise. USP5 was also able to enhance the expression of RIG-I and activate higher levels of IFNφ1 stimulated by Poly (I: C). Viral infection and interference tests demonstrated that USP5 inhibited the replication of SVCV in vitro. In summary, this study reveals that USP5 is able to activate higher levels of interferon by increasing RIG-I protein levels, and thus implement antivirus functions.

Variable Surface Glycoprotein from Trypanosoma brucei Undergoes Cleavage by Matrix Metalloproteinases: An in silico Approach.

In order to survive as extracellular parasites in the mammalian host environment, Trypanosoma brucei has developed efficient mechanisms of immune system evasion, which include the abundant expression of a variable surface glycoprotein (VSG) coat. VSGs are anchored in the parasite membrane by covalent C-terminal binding to glycosylphosphatidylinositol and may be periodically removed by a phospholipase C (PLC) and a major surface protein (TbMSP). VSG molecules show extraordinary antigenic diversity and a comparative analysis of protein sequences suggests that conserved elements may be a suitable target against African trypanosomiasis. However, the cleavage mechanisms of these molecules remain unclear. Moreover, in protozoan infections, including those caused by Trypanosoma brucei, it is possible to observe an increased expression of the matrix metalloproteinases (MMPs). To address the cleavage mechanism of VSGs, the PROSPER server was used for the identification of VSG sequence cleavage sites. After data compilation, it was observed that 64 VSG consensus sequences showed a high conservation of hydrophobic residues, such as valine (V), methionine (M), leucine (L) and isoleucine (I) in the fifth position—the exact location of the cleavage site. In addition, the PROSPER server identified conserved cleavage site portions of VSG proteins recognized by three matrix metalloproteases (gelatinases: MMP-2, MMP-3 and MMP-9). However, further biological studies are needed in order to analyze and confirm this prediction.

Comparative analyses of whole-genome protein sequences from multiple organisms

Phylogenies based on entire genomes are a powerful tool for reconstructing the Tree of Life. Several methods have been proposed, most of which employ an alignment-free strategy. Average sequence similarity methods are different than most other whole-genome methods, because they are based on local alignments. However, previous average similarity methods fail to reconstruct a correct phylogeny when compared against other whole-genome trees. In this study, we developed a novel average sequence similarity method. Our method correctly reconstructs the phylogenetic tree of in silico evolved E. coli proteomes. We applied the method to reconstruct a whole-proteome phylogeny of 1,087 species from all three domains of life, Bacteria, Archaea, and Eucarya. Our tree was automatically reconstructed without any human decisions, such as the selection of organisms. The tree exhibits a concentric circle-like structure, indicating that all the organisms have similar total branch lengths from their common ancestor. Branching patterns of the members of each phylum of Bacteria and Archaea are largely consistent with previous reports. The topologies are largely consistent with those reconstructed by other methods. These results strongly suggest that this approach has sufficient taxonomic resolution and reliability to infer phylogeny, from phylum to strain, of a wide range of organisms.

An atypical strictosidine synthase, OsSTRL2, plays key roles in anther development and pollen wall formation in rice

Strictosidine synthase (STR) plays an important role in the biosynthesis of terpenoid indole alkaloids (TIAs) and is expressed in a range of active meristematic tissues of higher plants. STR proteins are involved in different physiological and biochemical pathways. However, the function of STR proteins in rice development remains poorly understood. In this study, we identified 21 possible STR-like (OsSTRL) family members in rice genome and found that only one gene, OsSTRL2, exhibited a pre-emergency specific florescence expression pattern. Tissue-specific expression profile analysis, β-glucuronidase histochemical (GUS) staining and RNA in situ hybridization confirmed that OsSTRL2 was highly expressed in tapetal cells and microspores. Comparative protein sequence analysis indicated that OsSTRL2 lacked the key catalytic residue found in a typical STR (STR1), although it possessed conserved β-propellers and α-helices formed the basic structure of STR1. OsSTRL2 knockout mutant resulted to male sterility because of the defects in anther development and pollen wall formation. Subcellular localization of OsSTRL2-YFP revealed that the OsSTRL2 protein was primarily localized in the endoplasmic reticulum (ER). Therefore, OsSTRL2 is an atypical strictosidine synthase that plays crucial roles in regulating anther development and pollen wall formation in rice.

Molecular characterization and expression analysis of heat shock protein 70 and 90 from Hermetia illucens reared in a food waste bioconversion pilot plant

Two full-length cDNAs of heat shock protein (HSP) genes (Hihsp70 and Hihsp90) were cloned from the black soldier fly (BSF) Hermetia illucens larvae reared in a food waste bioconversion pilot plant. The Hihsp70 and Hihsp90 transcripts were 2243 and 2507bp long, contained 1923 and 2166bp open reading frames encoding proteins of 640 and 721 amino acids with a molecular mass of 69.8 and 83kDa, respectively. Comparative analysis of protein sequences revealed the presence of the conserved HSP motifs in both proteins, showing high homology to their counterparts in other insect species from six different orders. Hihsp70 and Hihsp90 transcriptional expression profiles during two key developmental stages in the bioconversion process were evaluated by quantitative real time PCR showing that both genes were modulated during larval development. HiHsp70 mRNA expression levels during the II instar larvae was higher in respect to the V instar larvae. A similar difference in mRNA expression levels, but in a less extent, was found for the Hihsp90. Moreover, a diverse transcript level between the two genes at the V larval stage was observed where Hihsp90 was up-regulated compared to Hihsp70.These results suggested the involvement of Hsp70 and Hsp90 in H. illucens development and provide further evidences on the ecological and evolutionary importance of HSPs in the insect developmental processes together with valuable information on molecular features of adaptability to peculiar rearing conditions during food waste bioconversion.

Comparative Analysis of Immune Checkpoint Molecules and Their Potential Role in the Transmissible Tasmanian Devil Facial Tumor Disease.

Immune checkpoint molecules function as a system of checks and balances that enhance or inhibit immune responses to infectious agents, foreign tissues, and cancerous cells. Immunotherapies that target immune checkpoint molecules, particularly the inhibitory molecules programmed cell death 1 and cytotoxic T-lymphocyte-associated protein 4 (CTLA-4), have revolutionized human oncology in recent years, yet little is known about these key immune signaling molecules in species other than primates and rodents. The Tasmanian devil facial tumor disease is caused by transmissible cancers that have resulted in a massive decline in the wild Tasmanian devil population. We have recently demonstrated that the inhibitory checkpoint molecule PD-L1 is upregulated on Tasmanian devil (Sarcophilus harrisii) facial tumor cells in response to the interferon-gamma cytokine. As this could play a role in immune evasion by tumor cells, we performed a thorough comparative analysis of checkpoint molecule protein sequences among Tasmanian devils and eight other species. We report that many of the key signaling motifs and ligand-binding sites in the checkpoint molecules are highly conserved across the estimated 162 million years of evolution since the last common ancestor of placental and non-placental mammals. Specifically, we discovered that the CTLA-4 (MYPPPY) ligand-binding motif and the CTLA-4 (GVYVKM) inhibitory domain are completely conserved across all nine species used in our comparative analysis, suggesting that the function of CTLA-4 is likely conserved in these species. We also found that cysteine residues for intra- and intermolecular disulfide bonds were also highly conserved. For instance, all 20 cysteine residues involved in disulfide bonds in the human 4-1BB molecule were also present in devil 4-1BB. Although many key sequences were conserved, we have also identified immunoreceptor tyrosine-based inhibitory motifs (ITIMs) and immunoreceptor tyrosine-based switch motifs (ITSMs) in genes and protein domains that have not been previously reported in any species. This checkpoint molecule analysis and review of salient features for each of the molecules presented here can serve as road map for the development of a Tasmanian devil facial tumor disease immunotherapy. Finally, the strategies can be used as a guide for veterinarians, ecologists, and other researchers willing to venture into the nascent field of wild immunology.

Phylogenomic Analyses and Comparative Studies on Genomes of the Bifidobacteriales: Identification of Molecular Signatures Specific for the Order Bifidobacteriales and Its Different Subclades.

The order Bifidobacteriales comprises a diverse variety of species found in the gastrointestinal tract of humans and other animals, some of which are opportunistic pathogens, whereas a number of others exhibit health-promoting effects. However, currently very few biochemical or molecular characteristics are known which are specific for the order Bifidobacteriales, or specific clades within this order, which distinguish them from other bacteria. This study reports the results of detailed comparative genomic and phylogenetic studies on 62 genome-sequenced species/strains from the order Bifidobacteriales. In a robust phylogenetic tree for the Bifidobacteriales constructed based on 614 core proteins, a number of well-resolved clades were observed including a clade separating the Scarodvia-related genera (Scardovia clade) from the genera Bifidobacterium and Gardnerella, as well as a number of previously reported clusters of Bifidobacterium spp. In parallel, our comparative analyses of protein sequences from the Bifidobacteriales genomes have identified numerous molecular markers that are specific for this group of bacteria. Of these markers, 32 conserved signature indels (CSIs) in widely distributed proteins and 10 signature proteins are distinctive characteristics of all sequenced Bifidobacteriales species and provide novel and highly specific means for distinguishing these bacteria. In addition, multiple other molecular signatures are specific for the following clades of Bifidobacteriales: (i) 5 CSIs specific for a clade comprising of the Scardovia-related genera; (ii) 3 CSIs and 2 CSPs specific for a clade consisting of the Bifidobacterium and Gardnerella spp.; (iii) multiple other signatures demarcating a number of clusters of the B. asteroides-and B. longum- related species. The described molecular markers provide novel and reliable means for distinguishing the Bifidobacteriales and a number of their clades in molecular terms and for the classification of these bacteria. The Bifidobacteriales-specific CSIs, found in important proteins, are predicted to play important roles in modifying the cellular functions of the affected proteins. Hence, biochemical studies on the cellular functions of these CSIs could lead to discovery of novel characteristics of either all Bifidobacteriales, or specific groups of bacteria within this order. Some of the functions affected/modified by these genetic changes could also be important for the probiotic/pathogenic activities of the bifidobacteria.

Evolutionary Genomics Suggests That CheV Is an Additional Adaptor for Accommodating Specific Chemoreceptors within the Chemotaxis Signaling Complex

Escherichia coli and Salmonella enterica are models for many experiments in molecular biology including chemotaxis, and most of the results obtained with one organism have been generalized to another. While most components of the chemotaxis pathway are strongly conserved between the two species, Salmonella genomes contain some chemoreceptors and an additional protein, CheV, that are not found in E. coli. The role of CheV was examined in distantly related species Bacillus subtilis and Helicobacter pylori, but its role in bacterial chemotaxis is still not well understood. We tested a hypothesis that in enterobacteria CheV functions as an additional adaptor linking the CheA kinase to certain types of chemoreceptors that cannot be effectively accommodated by the universal adaptor CheW. Phylogenetic profiling, genomic context and comparative protein sequence analyses suggested that CheV interacts with specific domains of CheA and chemoreceptors from an orthologous group exemplified by the Salmonella McpC protein. Structural consideration of the conservation patterns suggests that CheV and CheW share the same binding spot on the chemoreceptor structure, but have some affinity bias towards chemoreceptors from different orthologous groups. Finally, published experimental results and data newly obtained via comparative genomics support the idea that CheV functions as a “phosphate sink” possibly to off-set the over-stimulation of the kinase by certain types of chemoreceptors. Overall, our results strongly suggest that CheV is an additional adaptor for accommodating specific chemoreceptors within the chemotaxis signaling complex.

HistoneDB 2.0: a histone database with variants—an integrated resource to explore histones and their variants

Compaction of DNA into chromatin is a characteristic feature of eukaryotic organisms. The core (H2A, H2B, H3, H4) and linker (H1) histone proteins are responsible for this compaction through the formation of nucleosomes and higher order chromatin aggregates. Moreover, histones are intricately involved in chromatin functioning and provide a means for genome dynamic regulation through specific histone variants and histone post-translational modifications. ‘HistoneDB 2.0 – with variants’ is a comprehensive database of histone protein sequences, classified by histone types and variants. All entries in the database are supplemented by rich sequence and structural annotations with many interactive tools to explore and compare sequences of different variants from various organisms. The core of the database is a manually curated set of histone sequences grouped into 30 different variant subsets with variant-specific annotations. The curated set is supplemented by an automatically extracted set of histone sequences from the non-redundant protein database using algorithms trained on the curated set. The interactive web site supports various searching strategies in both datasets: browsing of phylogenetic trees; on-demand generation of multiple sequence alignments with feature annotations; classification of histone-like sequences and browsing of the taxonomic diversity for every histone variant. HistoneDB 2.0 is a resource for the interactive comparative analysis of histone protein sequences and their implications for chromatin function.Database URL: http://www.ncbi.nlm.nih.gov/projects/HistoneDB2.0

First report on the gastropod proapoptotic AIF3 counterpart from disk abalone (Haliotis discus discus) deciphering its transcriptional modulation by induced pathogenic stress.

Apoptosis inducing factor (AIF) is a flavoprotein that is involved in oxidative phosphorylation and induces apoptosis in eukaryotic cells. There are three isozymes of AIF that have been identified to date, designated as AIF1, AIF2, and AIF3; the human AIF3 is also known as an AIF-like protein (AIFL). This study aimed to identify and characterize a homologue of AIF3 from disk abalone (AbAIF3) that belongs to the phylum Mollusca. The open reading frame (ORF) of AbAIF3 is 1749 base pairs (bp) in length and encodes a protein of 583 amino acids, with a predicted molecular mass of 63.14 kDa. Based on our in-silico analysis, the AbAIF3 protein harbored the typical domain architecture as that of the known AIF family proteins, consisting of N-terminal Rieske and pyridine nucleotide-disulphide oxidoreductase domain. Comparative protein sequence analysis confirmed that AbAIF3 is a homolog of AIF3. Moreover, our phylogenetic analysis revealed that AbAIF3 had a close evolutionary relationship with the molluscan counterparts. Interestingly, AbAIF3 was shown to induce apoptosis in HEK293T cells using transfection assays followed by flow cytometric analysis. In addition, we found that AbAIF3 mRNA expression was ubiquitous in physiologically important tissues, and significantly modulated upon experimental immune stimulations in hemocytes. Collectively, our study illustrates the indispensable role of AbAIF3 in inducing apoptosis in disk abalones, which in turn might be involved in hosts' immune defense mechanisms against microbial infections.

A proposed new bacteriophage subfamily: "Jerseyvirinae".

Based on morphology and comparative nucleotide and protein sequence analysis, a new subfamily of the family Siphoviridae is proposed, named "Jerseyvirinae" and consisting of three genera, "Jerseylikevirus", "Sp3unalikevirus" and "K1glikevirus". To date, this subfamily consists of 18 phages for which the genomes have been sequenced. Salmonella phages Jersey, vB_SenS_AG11, vB_SenS-Ent1, vB_SenS-Ent2, vB_SenS-Ent3, FSL SP-101, SETP3, SETP7, SETP13, SE2, SS3e and wksl3 form the proposed genus "Jerseylikevirus". The proposed genus "K1glikevirus" consists of Escherichia phages K1G, K1H, K1ind1, K1ind2 and K1ind3. The proposed genus "Sp3unalikevirus" contains one member so far. Jersey-like phages appear to be widely distributed, as the above phages were isolated in the UK, Canada, the USA and South Korea between 1970 and the present day. The distinguishing features of this subfamily include a distinct siphovirus morphotype, genomes of 40.7-43.6 kb (49.6-51.4 mol % G+C), a syntenic genome organisation, and a high degree of nucleotide sequence identity and shared proteins. All known members of the proposed subfamily are strictly lytic.