Positive Selection in Bifidobacterium Genes Drives Species-Specific Host-Bacteria Communication

paPAML: An Improved Computational Tool to Explore Selection Pressure on Protein-Coding Sequences

Evolution is change over time. Although neutral changes promoted by drift effects are most reliable for phylogenetic reconstructions, selection-relevant changes are of only limited use to reconstruct phylogenies. On the other hand, comparative analyses of neutral and selected changes of protein-coding DNA sequences (CDS) retrospectively tell us about episodic constrained, relaxed, and adaptive incidences. The ratio of sites with nonsynonymous (amino acid altering) versus synonymous (not altering) mutations directly measures selection pressure and can be analysed by using the Phylogenetic Analysis by Maximum Likelihood (PAML) software package. We developed a CDS extractor for compiling protein-coding sequences (CDS-extractor) and parallel PAML (paPAML) to simplify, amplify, and accelerate selection analyses via parallel processing, including detection of negatively selected sites. paPAML compiles results of site, branch-site, and branch models and detects site-specific negative selection with the output of a codon list labelling significance values. The tool simplifies selection analyses for casual and inexperienced users and accelerates computing speeds up to the number of allocated computer threads. We then applied paPAML to examine the evolutionary impact on a new GINS Complex Subunit 3 exon, and neutrophil-associated as well as lysin and apolipoprotein genes. Compared with codeml (PAML version 4.9j) and HyPhy (HyPhy FEL version 2.5.26), all paPAML test runs performed with 10 computing threads led to identical selection pressure results, whereas the total selection analysis via paPAML, including all model comparisons, was about 3 to 5 times faster than the longest running codeml model and about 7 to 15 times faster than the entire processing time of these codeml runs.

The Role of the PFNA Operon of Bifidobacteria in the Recognition of Host's Immune Signals: Prospects for the Use of the FN3 Protein in the Treatment of COVID-19

Bifidobacteria are some of the major agents that shaped the immune system of many members of the animal kingdom during their evolution. Over recent years, the question of concrete mechanisms underlying the immunomodulatory properties of bifidobacteria has been addressed in both animal and human studies. A possible candidate for this role has been discovered recently. The PFNA cluster, consisting of five core genes, pkb2, fn3, aaa-atp, duf58, tgm, has been found in all gut-dwelling autochthonous bifidobacterial species of humans. The sensory region of the species-specific serine-threonine protein kinase (PKB2), the transmembrane region of the microbial transglutaminase (TGM), and the type-III fibronectin domain-containing protein (FN3) encoded by the I gene imply that the PFNA cluster might be implicated in the interaction between bacteria and the host immune system. Moreover, the FN3 protein encoded by one of the genes making up the PFNA cluster, contains domains and motifs of cytokine receptors capable of selectively binding TNF-α. The PFNA cluster could play an important role for sensing signals of the immune system. Among the practical implications of this finding is the creation of anti-inflammatory drugs aimed at alleviating cytokine storms, one of the dire consequences resulting from SARS-CoV-2 infection.

FN3 protein fragment containing two type III fibronectin domains from B. longum GT15 binds to human tumor necrosis factor alpha in vitro

Most species of the genus Bifidobacterium contain the gene cluster PFNA, which is presumably involved in the species-specific communication between bacteria and their hosts. The gene cluster PFNA consists of five genes including fn3, which codes for a protein containing two fibronectin type III domains. Each fibronectin domain contains sites similar to cytokine-binding sites of human receptors. Based on this finding we assumed that this protein would bind specifically to human cytokines in vitro. We cloned a fragment of the fn3 gene (1503 bp; 501 aa) containing two fibronectin domains, from the strain B. longum subsp. longum GT15. After cloning the fragment into the expression vector pET16b and expressing it in E. coli, the protein product was purified to a homogenous state for further analysis. Using the immunoferment method, we tested the purified fragment's ability to bind the following human cytokines: IL-1β, IL-6, IL-10, TNFα. We developed a sandwich ELISA system to detect any specific interactions between the purified protein and any of the studied cytokines. We found that the purified protein fragment only binds to TNFα.

Phylogenetic characterization of two novel species of the genus Bifidobacterium: Bifidobacterium saimiriisciurei sp. nov. and Bifidobacterium platyrrhinorum sp. nov

Three bifidobacterial Gram-stain-positive, non-spore forming and fructose-6-phosphate phosphoketolase-positive strains, SMA1^T, SMB2 and SMA15^T were isolated from the faeces of two adult males of the squirrel monkey (Saimiri sciureus). On the basis of 16S rRNA gene sequence similarities, the type strain of Bifidobacterium primatium DSM 100687^T (99.3%; similarity) was the closest neighbour to strains SMA1^T and SMB2, whereas the type strain of Bifidobacterium stellenboschense DSM 23968^T (96.5%) was the closest neighbour to strain SMA15^T. The average nucleotide identity (ANI) values of SMA1^T and SAM15^T with the closely related type strains were 93.7% and 88.1%, respectively. The in silico DNA‒DNA hybridization values with the closest neighbours were 53.1% and 36.9%, respectively. GC contents of strains SMA1^T and SMA15^T were 63.6 and 66.4 mol%, respectively. Based on the phylogenetic, genotypic and phenotypic data obtained, the strains SMA1^T and SMA15^T clearly represent two novel taxa within the genus Bifidobacterium for which the names Bifidobacterium saimiriisciurei sp. nov. (type strain SMA1^T = BCRC 81223^T = NBRC 114049^T = DSM 106020^T) and Bifidobacterium platyrrhinorum sp. nov. (type strain SMA15^T = BCRC 81224^T = NBRC 114051^T = DSM 106029^T) are proposed.