Prediction of protein architectures involved in the signaling-pathway initiating sporulation in Firmicutes
BMC Research Notes volume 12, Article number: 686 (2019)
Like many other proteins, those belonging to the signal transduction cascade initiating sporulation (Spo0 pathway) have conserved protein domains (Capra and Laub in Annu Rev Microbiol 66:325–47, 2012). Improvements in bioinformatics applications to discover proteins involved in the initiation of the sporulating cascade in newly sequenced genomes is an important task that requires rigorous comparative genomic methods and manual curation to identify endospore-forming bacteria. This note aims to present a collection of predicted proteins involved in the Spo0 pathway found in the proteomes of fully sequenced and manually curated endospore-forming Firmicutes species. This collection may serve as a guide to conduct future experiments in endospore formers in genomic and metagenomic projects.
Similar to the report of Davidson et al. (PLoS Genet 14:1–33, 2018), we used Pfam profiles (El-Gebali et al. in Nucleic Acids Res 47:D427–32, 2019) defining each protein and the genomic context surrounding the query gene to predict probable orthologs of the Spo0 pathway in Firmicutes. We present in this note a collection of 325 Firmicutes species organized by phylogenetic class and classified as spore formers, non-spore formers or unknown spore phenotype based on published literature, for which we predicted probable orthologs defining the signal transduction pathway initiating sporulation.
Comparative genomic analysis is a valuable tool to explore new genomes and metagenomes to search for functional processes, phylogenetic relationships, and evolutionary traits among organisms. An example of this comparative analysis is finding that proteins involved in the Spo0 pathway constitute of at least a sensor kinase, the phosphotransferases Spo0F and Spo0B, and the response regulator Spo0A , were suggested to be ancestral since it was found not only in Bacilli but also in some Clostridia . While an original article from our group that describes the distribution and prevalence of the Spo0 pathway and also confirms that this pathway is ancestral was under review, the results mentioned above were published by Davidson and collaborators using a set 84 Firmicutes . Therefore, in this note, we provide a more extensive list of manually curated Firmicutes possessing different forms of the Spo0 pathway initiating sporulation in Firmicutes.
The data are a collection of predicted protein architectures defining the proteins shaping the Spo0 pathway in Firmicutes. As a guide, we used the proteins shaping the sporulation cascades in the model organisms shown in Data File 1 . Using the well-curated profiles of the PfamA database , we constructed the architectures shown in Table 1, Data File 1 . The architectures were then used to inspect the genomes of 325 Firmicutes extracted from the Kyoto Encyclopedia of Genes and Genomes (KEGG) genome database [5,6,7], using the hmmscan program from the HMMER suite . We preserved those hits with an expectation value > 0.001. We discarded sequences that do not preserve the same kind of domains and with a length not longer than twice the length of the model. We used the genome neighborhood as a parameter of selection to discriminate homologous proteins not belonging to the Spo0 pathway. To this end, we inspected three genes upstream and downstream from the query gene. Frequently, neighbor genes encode proteins that were organized in a cluster of orthologous genes (COG) . COGs were assigned by finding homologous proteins for each query and neighbors using a hidden Markov model (HMM) search using the hmmsearch program . This HMM search process employs a previously constructed model set that represents each of the 4873 COGs [9, 10]. The genomes tested were annotated using Operon Mapper, which is able to classify genes into COGs [10, 11]. The list of the COGs assigned to the Spo0 proteins and their neighbors is available in Table 1, Data File 2 .
To detect orphan histidine kinases (HKo), defined as kinases not having as a neighboring response regulator , we performed a hmmscan using the profiles of the PfamA shown in Table 1, Data File 1  and kept those hits with an expectation value > 0.001 that fulfill the architectures described in Table 1, Data File 1 . As an additional condition, we discarded those sequences that did not have the same number of domains and with a length no more than twice the length of the model. As the architectures were constructed using proteins experimentally proven to participate in sporulation (Table 1, Data File 1) , we considered for this analysis the architecture found in Clostridium thermocellum, which is composed of a HisKA, an HATPase, and a response regulator domain (Table 1, Data File 1) .
Spo0 pathways found for each endospore former are shown in Table 1, Data File 3 . These results show the architectures describing twelve HKos, three architectures representing transferases, and two architectures describing the response regulator Spo0A. The final table includes 185 Bacilli, 134 Clostridia, 4 Negativicutes, and 2 Erysipelotrichia species.
The data presented in this note were not published as a research article since the main findings were published while our manuscript was under review. Nonetheless, the data are still important and useful since they are a collection of well-curated Firmicutes species for which a Spo0 pathway was present or absent.
Experimental evidence of a sporulation phenotype is not available for a vast number of strains since the conditions have not been identified.
New protein profiles should be constructed to identify novel orphan kinases.
A study should be performed in newly sequenced genomes to increase the collection; nonetheless, the groups presented in this note are probably enough to confirm the main conclusions found by both groups.
hidden Markov model
cluster of orthologous genes
Capra EJ, Laub MT. Evolution of two-component signal transduction systems. Annu Rev Microbiol. 2012;66:325–47.
Davidson P, Eutsey R, Redler B, Hiller NL, Laub MT, Durand D. Flexibility and constraint: evolutionary remodeling of the sporulation initiation pathway in firmicutes. PLoS Genet. 2018;14:1–33.
Data file 1. Proteins model architectures of the signal transduction cascade in endospore formers. 2019. https://doi.org/10.6084/m9.figshare.9701522.
El-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, et al. The Pfam protein families database in 2019. Nucleic Acids Res. 2019;47:D427–32.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases, and drugs. Nucleic Acids Res. 2017;45:D353–61.
Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016;44:D457–62.
KEGG Organisms: complete genomes. https://www.genome.jp/kegg/catalog/org_list.html.
Finn RD, Clements J, Arndt W, Miller BL, Wheeler TJ, Schreiber F, et al. HMMER web server: 2015 update. Nucleic Acids Res. 2015;43:W30–8.
Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, et al. The COG database: an updated version includes eukaryotes. BMC Bioinform. 2003;4:1–14.
Taboada B, Verde C, Merino E. High accuracy operon prediction method based on STRING database scores. Nucleic Acids Res. 2010;38:1–10.
Taboada B, Estrada K, Ciria R, Merino E. Operon-mapper: a web server for precise operon identification in bacterial and archaeal genomes. Bioinformatics. 2018;34:4118–20.
Data file 2. COG defining the Spo0 pathway. 2019. https://doi.org/10.6084/m9.figshare.9701630.
Dürre P. Ancestral sporulation initiation. Mol Microbiol. 2011;80:584–7.
Data file 3. Spo0 pathwaySpo0 predicted pathways and curated Firmicutes. 2019. https://doi.org/10.6084/m9.figshare.9876683.
We thank Ricardo Ciria, Walter Santos, Juan Manuel Hurtado, David Santiago Castañeda Carreón, and Shirley Ainsworth for their technical and bibliographical support.
All data and materials described herein are available in the links provided in the main text.
The design of the study and data collection was supported by the grant IN202118 from PAPIIT-UNAM awarded to RMGR. A CONACyT scholarship supported Paola Martinez-Amador, Lizeth Soto-Avila, and Nori Castañeda-Gomez. Antonio Loza was supported by a postdoctoral scholarship of the Consorcio de Investigación del Golfo de México (CIGoM), which is supported by the National Council of Science and Technology of Mexico—Mexican Ministry of Energy—Hydrocarbon Trust, project 201441.
Ethics approval and consent to participate
Consent for publication
The author declares that she has no competing interests.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Martinez-Amador, P., Castañeda, N., Loza, A. et al. Prediction of protein architectures involved in the signaling-pathway initiating sporulation in Firmicutes. BMC Res Notes 12, 686 (2019). https://doi.org/10.1186/s13104-019-4712-3