Peroxide turnover and signalling are involved in many biological phenomena relevant to human diseases. Yet, all the players and mechanisms involved in peroxide perception are not known. Elucidating very remote evolutionary relationships between proteins is an approach that allows the discovery of novel protein functions. Here, we start with three human proteins, SRPX, SRPX2 and CCDC80, involved in tumor suppression and progression, which possess a conserved region of similarity. Structure and function prediction allowed the definition of P-DUDES, a phylogenetically widespread, possibly ancient protein structural domain, common to vertebrates and many bacterial species.

Evolutionary advantages over cousin cells in bacterial pathogens may decide about the success of a specific cell in its environment. Bacteria use a plethora of methods to defend against other cells and many devices to attack their opponents when competing for resources. Bacteriocins are antibacterial proteins that are used to eliminate competition. We report the discovery of a putative membrane-bound bacteriocin encoded by the Bacillus anthracis pathogenic pXO1 plasmid. We analyze the genomic structure of the bacteriocin operon. The proposed mechanisms of action predestine this operon as a potent competitive advantage over cohabitants of the same niche.

The zinc-dependent metalloproteases with His-Glu-x-x-His (HExxH) active site motif, zincins, are a broad group of proteins involved in many metabolic and regulatory functions, and found in all forms of life. Human genome contains more than 100 genes encoding proteins with known zincin-like domains. A survey of all proteins containing the HExxH motif shows that approximately 52% of HExxH occurrences fall within known protein structural domains (as defined in the Pfam database). Domain families with majority of members possessing a conserved HExxH motif include, not surprisingly, many known and putative metalloproteases. Furthermore, several HExxH-containing protein domains thus identified can be confidently predicted to be putative peptidases of zincin fold. Thus, we predict zincin-like fold for eight uncharacterised Pfam families. Besides the domains with the HExxH motif strictly conserved, and those with sporadic occurrences, intermediate families are identified that contain some members with a conserved HExxH motif, but also many homologues with substitutions at the conserved positions. Such substitutions can be evolutionarily conserved and non-random, yet functional roles of these inactive zincins are not known. The CLCAs are a novel zincin-like protease family with many cases of substituted active sites. We show that this allegedly metazoan family has a number of bacterial and archaeal members. An extremely patchy phylogenetic distribution of CLCAs in prokaryotes and their conserved protein domain composition strongly suggests an evolutionary scenario of horizontal gene transfer (HGT) from multicellular eukaryotes to bacteria, providing an example of eukaryote-derived xenologues in bacterial genomes. Additionally, in a protein family identified here as closely homologous to CLCA, the CLCA_X (CLCA-like) family, a number of proteins is found in phages and plasmids, supporting the HGT scenario.

Background. Members of the Runx gene family encode transcription factors that bind to DNA in a sequence-specific manner. Among the three Runx proteins, Runx2 comprises 607 amino acid (aa) residues, is expressed in bone, and plays crucial roles in osteoblast differentiation and bone development. We examined whether the Runx2 gene is also expressed in testes. Methods. Murine testes from 1-, 2-, 3-, 4-, and 10-week-old male mice of the C57BL/6J strain and W∕W (v) strain were used throughout the study. Northern Blot Analyses were performed using extracts form the murine testes. Sequencing of cDNA clones and 5'-rapid amplification of cDNA ends were performed to determine the full length of the transcripts, which revealed that the testicular Runx2 comprises 106 aa residues coding novel protein. Generating an antiserum using the amino-terminal 15 aa of Runx2 (Met(1) to Gly(15)) as an antigen, immunoblot analyses were performed to detect the predicted polypeptide of 106 aa residues with the initiating Met(1). With the affinity-purified anti-Runx2 antibody, immunohistochemical analyses were performed to elucidate the localization of the protein. Furthermore, bioinformatic analyses were performed to predict the function of the protein. Results. A Runx2 transcript was detected in testes and was specifically expressed in germ cells. Determination of the transcript structure indicated that the testicular Runx2 is a splice isoform. The predicted testicular Runx2 polypeptide is composed of only 106 aa residues, lacks a Runt domain, and appears to be a basic protein with a predominantly alpha-helical conformation. Immunoblot analyses with an anti-Runx2 antibody revealed that Met(1) in the deduced open reading frame of Runx2 is used as the initiation codon to express an 11 kDa protein. Furthermore, immunohistochemical analyses revealed that the Runx2 polypeptide was located in the nuclei, and was detected in spermatocytes at the stages of late pachytene, diplotene and second meiotic cells as well as in round spermatids. Bioinformatic analyses suggested that the testicular Runx2 is a histone-like protein. Discussion. A variant of Runx2 that differs from the bone isoform in its splicing is expressed in pachytene spermatocytes and round spermatids in testes, and encodes a histone-like, nuclear protein of 106 aa residues. Considering its nuclear localization and differentiation stage-dependent expression, Runx2 may function as a chromatin-remodeling factor during spermatogenesis. We thus conclude that a single Runx2 gene can encode two different types of nuclear proteins, a previously defined transcription factor in bone and cartilage and a short testicular variant that lacks a Runt domain.

A number of bacterial virulence factors have been observed to adopt structures similar to that of aerolysin, the principal toxin of Aeromonas species. However, a comprehensive description of architecture and structure of the aerolysin-like superfamily has not been determined. In this study, we define a more compact aerolysin-like domain--or aerolysin fold--and show that this domain is far more widely spread than anticipated since it can be found throughout kingdoms. The aerolysin-fold could be found in very diverse domain and functional contexts, although a toxic function could often be assigned. Due to this diversity, the borders of the superfamily could not be set on a sequence level. As a border-defining member, we therefore chose pXO2-60--a protein from the pathogenic pXO2 plasmid of Bacillus anthracis. This fascinating protein, which harbors a unique ubiquitin-like fold domain at the C-terminus of the aerolysin-domain, nicely illustrates the diversity of the superfamily. Its putative role in the virulence of B. anthracis and its three dimensional model are discussed.

The rare human disorder chorea-acanthocytosis (ChAc) is caused by mutations in hVPS13A gene. The hVps13A protein interacts with actin and regulates the level of phosphatidylinositol 4-phosphate (PI4P) in the membranes of neuronal cells. Yeast Vps13 is involved in vacuolar protein transport and, like hVps13A, participates in PI4P metabolism. Vps13 proteins are conserved in eukaryotes, but their molecular function remains unknown. One of the mutations found in ChAc patients causes amino acids substitution I2771R which affects the localization of hVps13A in skeletal muscles. To dissect the mechanism of pathogenesis of I2771R, we created and analyzed a yeast strain carrying the equivalent mutation. Here we show that in yeast, substitution I2749R causes dysfunction of Vps13 protein in endocytosis and vacuolar transport, although the level of the protein is not affected, suggesting loss of function. We also show that Vps13, like hVps13A, influences actin cytoskeleton organization and binds actin in immunoprecipitation experiments. Vps13-I2749R binds actin, but does not function in the actin cytoskeleton organization. Moreover, we show that Vps13 binds phospholipids, especially phosphatidylinositol 3-phosphate (PI3P), via its SHR_BD and APT1 domains. Substitution I2749R attenuates this ability. Finally, the localization of Vps13-GFP is altered when cellular levels of PI3P are decreased indicating its trafficking within the endosomal membrane system. These results suggest that PI3P regulates the functioning of Vps13, both in protein trafficking and actin cytoskeleton organization. Attenuation of PI3P-binding ability in the mutant hVps13A protein may be one of the reasons for its mislocalization and disrupted function in cells of patients suffering from ChAc.

Pseudomonas syringae pv. phaseolicola is the causative agent of halo blight in common bean (Phaseolus vulgaris). Similar to other pathogenic gram-negative bacteria, it secrets a set of type III effectors into host cells to subvert defense mechanisms. HopQ1 (for Hrp outer protein Q) is one of these type III effectors contributing to virulence of bacteria. Upon delivery into a plant cell, HopQ1 undergoes phosphorylation, binds host 14-3-3 proteins and suppresses defense-related signaling. Some plants however, evolved systems to recognize HopQ1 and respond to its presence and thus to prevent infection. HopQ1 shows homology to Nucleoside Hydrolases (NHs), but it contains a modified calcium binding motif not found in the canonical enzymes. CLuster ANalysis of Sequences (CLANS) revealed that HopQ1 and alike proteins make a distinct group of putative NHs located distantly from the classical enzymes. The HopQ1 - like protein (HLP) group comprises sequences from plant pathogenic bacteria, fungi, and lower plants. Our data suggest that the evolution of HopQ1 homologs in bacteria, fungi, and algae was independent. The location of moss HopQ1 homologs inside the fungal clade indicates a possibility of horizontal gene transfer (HGT) between those taxa. We identified a HLP in the moss Physcomitrella patens. Our experiments show that this protein (referred to as PpHLP) extended by a TTSS signal of HopQ1 promoted P. syringae growth in bean and was recognized by Nicotiana benthamiana immune system. Thus, despite the low sequence similarity to HopQ1 the engineered PpHLP acted as a bacterial virulence factor and displayed similar to HopQ1 virulence properties.

The extrachromosomal gene pool plays a significant role both in evolution and in the environmental adaptation of bacteria. The L. lactis subsp. lactis IL594 strain contains seven plasmids, named pIL1 to pIL7, and is the parental strain of the plasmid-free L. lactis IL1403, which is one of the best characterized lactococcal strains of LAB. Complete nucleotide sequences of pIL1 (6,382 bp), pIL2 (8,277 bp), pIL3 (19,244 bp), pIL4 (48,979), pIL5 (23,395), pIL6 (28,435 bp) and pIL7 (28,546) were established and deposited in the generally accessible database (GeneBank). Nine highly homologous repB-containing replicons, belonging to the lactococcal theta-type replicons, have been identified on the seven plasmids. Moreover, a putative region involved in conjugative plasmid mobilization was found on four plasmids, through identification of the presence of mob genes and/or oriT sequences. Detailed bioinformatic analysis of the plasmid nucleotide sequences provided new insight into the repertoire of plasmid-encoded functions in L. lactis, and indicated that plasmid genes from IL594 strain can be important for L. lactis adaptation to specific environmental conditions (e.g. genes coding for proteins involved in DNA repair or cold shock response) as well as for technological processes (e.g. genes encoding citrate and lactose utilization, oligopeptide transport, restriction-modification system). Moreover, global gene analysis indicated cooperation between plasmid- and chromosome-encoded metabolic pathways.

Here, we report the extensive bioinformatic and functional analyses of the unusual pLOCK 0919, a plasmid originating from the probiotic Lactobacillus casei LOCK 0919 strain. This plasmid is atypical because it harbors the spaCBA-srtC gene cluster encoding SpaCBA pili. We show that all other spaCBA-srtC sequences of the Lactobacillus genus that have been previously described and deposited in GenBank are present in the chromosomal DNA. Another important observation for pLOCK 0919 is that the spaCBA-srtC gene cluster and its surrounding genes are highly similar to the respective DNA region that is present in the most well-known and active SpaCBA pili producer, the probiotic Lactobacillus rhamnosus GG strain. Our results demonstrate that the spaCBA-srtC clusters of pLOCK 0919 and L. rhamnosus GG are genealogically similar, located in DNA regions that are rich in transposase genes and are poorly conserved among the publicly available sequences of Lactobacillus sp. In contrast to chromosomally localized pilus gene clusters from L. casei and Lactobacillus paracasei, the plasmidic spaC of L. casei LOCK 0919 is expressed and undergoes a slight glucose-induced repression. Moreover, results of series of in vitro tests demonstrate that L. casei LOCK 0919 has an adhesion potential, which is largely determined by the presence of the pLOCK 0919 plasmid. In particular, the plasmid occurrence positively influenced the hydrophobicity and aggregation abilities of L. casei LOCK 0919. Moreover, in vivo studies indicate that among the three Lactobacillus strains used to colonize the gastrointestinal tract of germ-free mice, already after 2 days of colonization, L. casei LOCK 0919 became the dominant strain and persisted there for at least 48 days.

The minimal set of proteins necessary to maintain a vertebrate cell forms an interesting core of cellular machinery. The known proteome of human red blood cell consists of about 1400 proteins. We treated this protein complement of one of the simplest human cells as a model and asked the questions on its function and origins. The proteome was mapped onto phylogenetic profiles, i.e. vectors of species possessing homologues of human proteins. A novel clustering approach was devised, utilising similarity in the phylogenetic spread of homologues as distance measure. The clustering based on phylogenetic profiles yielded several distinct protein classes differing in phylogenetic taxonomic spread, presumed evolutionary history and functional properties. Notably, small clusters of proteins common to vertebrates or Metazoa and other multicellular eukaryotes involve biological functions specific to multicellular organisms, such as apoptosis or cell-cell signaling, respectively. Also, a eukaryote-specific cluster is identified, featuring GTP-ase signalling and ubiquitination. Another cluster, made up of proteins found in most organisms, including bacteria and archaea, involves basic molecular functions such as oxidation-reduction and glycolysis. Approximately one third of erythrocyte proteins do not fall in any of the clusters, reflecting the complexity of protein evolution in comparison to our simple model. Basically, the clustering obtained divides the proteome into old and new parts, the former originating from bacterial ancestors, the latter from inventions within multicellular eukaryotes. Thus, the model human cell proteome appears to be made up of protein sets distinct in their history and biological roles. The current work shows that phylogenetic profiles concept allows protein clustering in a way relevant both to biological function and evolutionary history.

In this work, we report likely recurrent horizontal (lateral) gene transfer events of genes encoding pore-forming toxins of the aerolysin family between species belonging to different kingdoms of life. Clustering based on pairwise similarity and phylogenetic analysis revealed several distinct aerolysin sequence groups, each containing proteins from multiple kingdoms of life. These results strongly support at least six independent transfer events between distantly related phyla in the evolutionary history of one protein family and discount selective retention of ancestral genes as a plausible explanation for this patchy phylogenetic distribution. We discuss the possible roles of these proteins and show evidence for a convergent new function in two extant species. We hypothesize that certain gene families are more likely to be maintained following horizontal gene transfer from commensal or pathogenic organism to its host if they 1) can function alone; and 2) are immediately beneficial for the ecology of the organism, as in the case of pore-forming toxins which can be utilized in multicellular organisms for defense and predation.

The correct folding and dimerization of tubulins, before their addition to the microtubular structure, needs a group of conserved proteins called cofactors A to E. The biochemical analysis of cofactors gave an insight to their general functions, however not much is known about the domain structure and detailed, molecular function of these proteins.

The rapid spread of the COVID-19 demands immediate response from the scientific communities. Appropriate countermeasures mean thoughtful and educated choice of viral targets (epitopes). There are several articles that discuss such choices in the SARS-CoV-2 proteome, other focus on phylogenetic traits and history of the Coronaviridae genome/proteome. However none consider viral protein low complexity regions (LCRs). Recently we created the first methods that are able to compare such fragments.

The widespread occurrence of repetitive stretches of DNA in genomes of organisms across the tree of life imposes fundamental challenges for sequencing, genome assembly, and automated annotation of genes and proteins. This multi-level problem can lead to errors in genome and protein databases that are often not recognized or acknowledged. As a consequence, end users working with sequences with repetitive regions are faced with 'ready-to-use' deposited data whose trustworthiness is difficult to determine, let alone to quantify. Here, we provide a review of the problems associated with tandem repeat sequences that originate from different stages during the sequencing-assembly-annotation-deposition workflow, and that may proliferate in public database repositories affecting all downstream analyses. As a case study, we provide examples of the Atlantic cod genome, whose sequencing and assembly were hindered by a particularly high prevalence of tandem repeats. We complement this case study with examples from other species, where mis-annotations and sequencing errors have propagated into protein databases. With this review, we aim to raise the awareness level within the community of database users, and alert scientists working in the underlying workflow of database creation that the data they omit or improperly assemble may well contain important biological information valuable to others.

The relationships between polypeptide composition, sequence, structure and function have been puzzling biologists ever since first protein sequences were determined. Here, we study the statistics of occurrence of all possible pentapeptide sequences in known proteins. To compensate for the non-uniform distribution of individual amino acid residues in protein sequences, we investigate separately all possible permutations of every given amino acid composition. For the majority of permutation groups we find that pentapeptide occurrences deviate strongly from the expected binomial distributions, and that the observed distributions are also characterized by high numbers of outlier sequences. An analysis of identified outliers shows they often contain known motifs and rare amino acids, suggesting that they represent important functional elements. We further compare the pentapeptide composition of regions known to correspond to protein domains with that of non-domain regions. We find that a substantial number of pentapeptides is clearly strongly favored in protein domains. Finally, we show that over-represented pentapeptides are significantly related to known functional motifs and to predicted ancient structural peptides.

Low complexity regions are fragments of protein sequences composed of only a few types of amino acids. These regions frequently occur in proteins and can play an important role in their functions. However, scientists are mainly focused on regions characterized by high diversity of amino acid composition. Similarity between regions of protein sequences frequently reflect functional similarity between them. In this article, we discuss strengths and weaknesses of the similarity analysis of low complexity regions using BLAST, HHblits and CD-HIT. These methods are considered to be the gold standard in protein similarity analysis and were designed for comparison of high complexity regions. However, we lack specialized methods that could be used to compare the similarity of low complexity regions. Therefore, we investigated the existing methods in order to understand how they can be applied to compare such regions. Our results are supported by exploratory study, discussion of amino acid composition and biological roles of selected examples. We show that existing methods need improvements to efficiently search for similar low complexity regions. We suggest features that have to be re-designed specifically for comparing low complexity regions: scoring matrix, multiple sequence alignment, e-value, local alignment and clustering based on a set of representative sequences. Results of this analysis can either be used to improve existing methods or to create new methods for the similarity analysis of low complexity regions.

Low complexity regions (LCRs) in protein sequences are characterized by a less diverse amino acid composition compared to typically observed sequence diversity. Recent studies have shown that LCRs may co-occur with intrinsically disordered regions, are highly conserved in many organisms, and often play important roles in protein functions and in diseases. In previous decades, several methods have been developed to identify regions with LCRs or amino acid bias, but most of them as stand-alone applications and currently there is no web-based tool which allows users to explore LCRs in protein sequences with additional functional annotations. We aim to fill this gap by providing PlaToLoCo - PLAtform of TOols for LOw COmplexity-a meta-server that integrates and collects the output of five different state-of-the-art tools for discovering LCRs and provides functional annotations such as domain detection, transmembrane segment prediction, and calculation of amino acid frequencies. In addition, the union or intersection of the results of the search on a query sequence can be obtained. By developing the PlaToLoCo meta-server, we provide the community with a fast and easily accessible tool for the analysis of LCRs with additional information included to aid the interpretation of the results. The PlaToLoCo platform is available at: http://platoloco.aei.polsl.pl/.

Selenoproteins serve important functions in many organisms, usually providing essential oxidoreductase enzymatic activity, often for defense against toxic xenobiotic substances. Most eukaryotic genomes possess a small number of these proteins, usually not more than 20. Selenoproteins belong to various structural classes, often related to oxidoreductase function, yet a few of them are completely uncharacterised.Here, the structural and functional prediction for the uncharacterised selenoprotein O (SELO) is presented. Using bioinformatics tools, we predict that SELO protein adopts a three-dimensional fold similar to protein kinases. Furthermore, we argue that despite the lack of conservation of the "classic" catalytic aspartate residue of the archetypical His-Arg-Asp motif, SELO kinases might have retained catalytic phosphotransferase activity, albeit with an atypical active site. Lastly, the role of the selenocysteine residue is considered and the possibility of an oxidoreductase-regulated kinase function for SELO is discussed.The novel kinase prediction is discussed in the context of functional data on SELO orthologues in model organisms, FMP40 a.k.a.YPL222W (yeast), and ydiU (bacteria). Expression data from bacteria and yeast suggest a role in oxidative stress response. Analysis of genomic neighbourhoods of SELO homologues in the three domains of life points toward a role in regulation of ABC transport, in oxidative stress response, or in basic metabolism regulation. Among bacteria possessing SELO homologues, there is a significant over-representation of aquatic organisms, also of aerobic ones. The selenocysteine residue in SELO proteins occurs only in few members of this protein family, including proteins from Metazoa, and few small eukaryotes (Ostreococcus, stramenopiles). It is also demonstrated that enterobacterial mchC proteins involved in maturation of bactericidal antibiotics, microcins, form a distant subfamily of the SELO proteins.The new protein structural domain, with a putative kinase function assigned, expands the known kinome and deserves experimental determination of its biological role within the cell-signaling network.

Transposable elements with long terminal direct repeats (LTR TEs) are one of the best studied groups of mobile elements. They are ubiquitous elements present in almost all eukaryotic genomes. Their number and state of conservation can be a highlight of genome dynamics. We searched all published fungal genomes for LTR-containing retrotransposons, including both complete, functional elements and remnant copies. We identified a total of over 66,000 elements, all of which belong to the Ty1/Copia or Ty3/Gypsy superfamilies. Most of the detected Gypsy elements represent Chromoviridae, i.e. they carry a chromodomain in the pol ORF. We analyzed our data from a genome-ecology perspective, looking at the abundance of various types of LTR TEs in individual genomes and at the highest-copy element from each genome. The TE content is very variable among the analyzed genomes. Some genomes are very scarce in LTR TEs (<50 elements), others demonstrate huge expansions (>8000 elements). The data shows that transposon expansions in fungi usually involve an increase both in the copy number of individual elements and in the number of element types. The majority of the highest-copy TEs from all genomes are Ty3/Gypsy transposons. Phylogenetic analysis of these elements suggests that TE expansions have appeared independently of each other, in distant genomes and at different taxonomical levels. We also analyzed the evolutionary relationships between protein domains encoded by the transposon pol ORF and we found that the protease is the fastest evolving domain whereas reverse transcriptase and RNase H evolve much slower and in correlation with each other.

Paracoccidioides is a fungal pathogen and the cause of paracoccidioidomycosis, a health-threatening human systemic mycosis endemic to Latin America. Infection by Paracoccidioides, a dimorphic fungus in the order Onygenales, is coupled with a thermally regulated transition from a soil-dwelling filamentous form to a yeast-like pathogenic form. To better understand the genetic basis of growth and pathogenicity in Paracoccidioides, we sequenced the genomes of two strains of Paracoccidioides brasiliensis (Pb03 and Pb18) and one strain of Paracoccidioides lutzii (Pb01). These genomes range in size from 29.1 Mb to 32.9 Mb and encode 7,610 to 8,130 genes. To enable genetic studies, we mapped 94% of the P. brasiliensis Pb18 assembly onto five chromosomes. We characterized gene family content across Onygenales and related fungi, and within Paracoccidioides we found expansions of the fungal-specific kinase family FunK1. Additionally, the Onygenales have lost many genes involved in carbohydrate metabolism and fewer genes involved in protein metabolism, resulting in a higher ratio of proteases to carbohydrate active enzymes in the Onygenales than their relatives. To determine if gene content correlated with growth on different substrates, we screened the non-pathogenic onygenale Uncinocarpus reesii, which has orthologs for 91% of Paracoccidioides metabolic genes, for growth on 190 carbon sources. U. reesii showed growth on a limited range of carbohydrates, primarily basic plant sugars and cell wall components; this suggests that Onygenales, including dimorphic fungi, can degrade cellulosic plant material in the soil. In addition, U. reesii grew on gelatin and a wide range of dipeptides and amino acids, indicating a preference for proteinaceous growth substrates over carbohydrates, which may enable these fungi to also degrade animal biomass. These capabilities for degrading plant and animal substrates suggest a duality in lifestyle that could enable pathogenic species of Onygenales to transfer from soil to animal hosts.