PePPER: a webserver for prediction of prokaryote promoter elements and regulons

Advances in human papillomavirus detection for cervical cancer screening and diagnosis: challenges of conventional methods and opportunities for emergent tools

Human papillomavirus (HPV) infection is the main cause of cervical cancer and other cancers such as anogenital and oropharyngeal cancers. The prevention screening and treatment of cervical cancer has remained one of the top priorities of the World Health Organization (WHO). In 2020, the WHO came up with the 90-70-90 strategy aimed at eliminating cervical cancers as a public health problem by the year 2030. One of the key priorities of this strategy is the recommendation for countries to ensure that 70% of their women are screened using a high-performance test by the age of 35, and again by the age of 45. Over the years, several traditional methods (notably, Pap smear and nucleic acid-based techniques) have been used for the detection of cervical cancer. While these methods have significantly reduced the incidence of cervical cancer and death, they still come short of excellence for the total eradication of HPV infection. The challenges include low sensitivity, low specificity, poor reproducibility, the need for high-level specialists, and the high cost of access to the facilities, to mention a few. Interestingly, however, several efforts are being made today to mitigate these challenges. In this review, we discussed the pros and cons of the traditional screening and testing of HPV infections, the efforts being made to improve their performances, and the emergent tools (especially, the electrochemical methods) that promise to revolutionize the screening and testing of HPV infections. The main aim of the review is to provide some novel clues to researchers that would allow for the development of high-performance, affordable, and triage-suitable electrochemical-based diagnostic tools for HPV and cervical cancer.

Characterization of two lytic bacteriophages infecting carbapenem-resistant clinical Klebsiella pneumoniae in Dhaka, Bangladesh

Bacteriophages or bacteria infecting viruses are genetically diverse. Due to the emergence of antimicrobial-resistant bacteria, lytic bacteriophages are gaining enormous attention for treating superbug infections. Klebsiella pneumoniae is one of the eight most significant nosocomial pathogens and is addressed as a critical priority pathogen by WHO, requiring alternative treatment options. We reported two highly lytic bacteriophages, Klebsiella phage Kpn BM7 and the novel Klebsiella phage Kpn BU9, isolated from hospital wastewater and exhibiting lytic activity against different clinical isolates. Whole-genome analysis revealed that phages BM7 and BU9 belong to class Caudoviricetes. Phage BM7, with a genome length of 170,558 bp, is a member of the genus Marfavirus and the species Marfavirus F48. While phage BU9, with a genome length of 60,450 bp, remains unclassified. Neither phage harbors any lysogenic, toxin, or antimicrobial resistance genes. Both phages can steadily survive up to 40 °C and at pH 5-7. The optimal MOI was 0.1 for BM7 and 1 for BU9, with short latent periods of 10 and 25 min and burst sizes of 85 PFU/cell and 12 PFU/cell, respectively. This is the first carbapenem-resistant K. pneumoniae targeting lytic phages to be reported from Bangladesh. This study suggests that BM7 and BU9 are potential candidates for targeting carbapenem-resistant K. pneumoniae.

From straight to curved: A historical perspective of DNA shape

DNA is the macromolecule responsible for storing the genetic information of a cell and it has intrinsic properties such as deformability, stability and curvature. DNA Curvature plays an important role in gene transcription and, consequently, in the subsequent production of proteins, a fundamental process of cells. With recent advances in bioinformatics and theoretical biology, it became possible to analyze and understand the involvement of DNA Curvature as a discriminatory characteristic of gene-promoting regions. These regions act as sites where RNAp (ribonucleic acid-polymerase) binds to initiate transcription. This review aims to describe the formation of Curvature, as well as highlight its importance in predicting promoters. Furthermore, this article provides the potential of DNA Curvature as a distinguishing feature for promoter prediction tools, as well as outlining the calculation procedures that have been described by other researchers. This work may support further studies directed towards the enhancement of promoter prediction software.

Characterization of four novel bacteriophages targeting multi-drug resistant Klebsiella pneumoniae strains of sequence type 147 and 307

The global dissemination of multi-drug resistant (MDR) pathogenic bacteria requires the rapid research and development of alternative therapies that can support or replace conventional antibiotics. Among MDR pathogens, carbapenem-resistant Klebsiella pneumoniae (CR-Kp) are of particular concern due to their extensive resistance profiles, global dissemination in hospital environments, and their major role in some life-threatening infections. Phages, or some of their components, are recognized as one of the potential alternatives that might be helpful to treat bacterial infections. In this study, we have isolated and characterized four lytic bacteriophages targeting K. pneumoniae strains of Sequence Type (ST) 307 or ST147, two predominant high-risk clones of CR-Kp. Phages, designated vB_KpS_GP-1, vB_KpP_GP-2, vB_KpP_GP-4, and vB_KpP_GP-5, were isolated from sewage wastewater samples. The vB_KpS_GP-1 phage was a siphovirus unable to establish lysogeny with its host, while the other three were podoviruses. While 85.7% of K. pneumoniae strains of ST307 were selectively lysed by the phages vB_KpS_GP-1 or vB_KpP_GP-5, the other two phages were able to lyse all the tested strains of ST147 (n = 12). Phages were stable over a broad pH and temperature range and were characterized by burst sizes of 10-100 plaque forming units and latency periods of 10-50 minutes. Genome sequencing confirmed the absence of antibiotic resistance genes, virulence factors or toxins and revealed that two phages were likely members of new genera. Given their strictly lytic nature and high selectivity towards two of the major high-risk clones of K. pneumoniae, cocktails of these phages could represent promising candidates for further evaluation in in vivo experimental models of K. pneumoniae infection.

The TolC and Lipopolysaccharide-Specific Escherichia coli Bacteriophage TLS-the Tlsvirus Archetype Virus

Introduction

TLS is a virulent bacteriophage of Escherichia coli that utilizes TolC and lipopolysaccharide as its cell surface receptors.

Methods

The genome was reannotated using the latest online resources and compared to other T1-like phages.

Results

The TLS genome consists of 49,902 base pairs, encoding 86 coding sequences that display considerable sequence similarity with the T1 phage genome. It also contains 18 intergenic 21-base long repeats, each of them upstream of a predicted start codon and in the direction of transcription. Data revealed that DNA packaging occurs through the pac site-mediated headful mechanism.

Conclusions

Based on sequence analysis of its genome, TLS belongs to the Drexlerviridae family and represents the type member of the Tlsvirus genus.

An insight into the role of branched-chain α-keto acid dehydrogenase (BKD) complex in branched-chain fatty acid biosynthesis and virulence of Listeria monocytogenes

Listeria monocytogenes is a foodborne bacterial pathogen that causes listeriosis. Positive regulatory factor A (PrfA) is a pleiotropic master activator of virulence genes of L. monocytogenes that becomes active upon the entry of the bacterium into the cytosol of infected cells. L. monocytogenes can survive and multiply at low temperatures; this is accomplished through the maintenance of appropriate membrane fluidity via branched-chain fatty acid (BCFA) synthesis. Branched-chain α-keto acid dehydrogenase (BKD), which is composed of four polypeptides encoded by lpd, bkdA1, bkdA2, and bkdB, is known to play a vital role in BCFA biosynthesis. Here, we constructed BKD-deficient Listeria strains by in-frame deletion of lpd, bkdA1, bkdA2, and bkdB genes. To determine the role in in vivo and in vitro, mouse model challenges, plaque assay in murine L2 fibroblast, and intracellular replication in J744A.1 macrophage were conducted. BKD-deficient strains exhibited defects in BCFA composition, virulence, and PrfA-regulon function within the host cells. Transcriptomics analysis revealed that the transcript level of the PrfA-regulon was lower in ΔbkdA1 strain than those in the wild-type. This study demonstrates that L. monocytogenes strains lacking BKD complex components were defective in PrfA-regulon function, and full activation of wild-type prfA may not occur within host cells in the absence of BKD. Further study will investigate the consequences of BKD deletion on PrfA function through altering BCFA catabolism.IMPORTANCEListeria monocytogenes is the causative agent of listeriosis, a disease with a high mortality rate. In this study, we have shown that the deletion of BKD can impact the function of PrfA and the PrfA-regulon. The production of virulence proteins within host cells is necessary for L. monocytogenes to promote its intracellular survival and is likely dependent on membrane integrity. We thus report a link between L. monocytogenes membrane integrity and the function of PrfA. This knowledge will increase our understanding of L. monocytogenes pathogenesis, which may provide insight into the development of antimicrobial agents.

From mec cassette to rdhA: a key Dehalobacter genomic neighborhood in a chloroform and dichloromethane-transforming microbial consortium

Chloroform (CF) and dichloromethane (DCM) are groundwater contaminants of concern due to their high toxicity and inhibition of important biogeochemical processes such as methanogenesis. Anaerobic biotransformation of CF and DCM has been well documented but typically independently of one another. CF is the electron acceptor for certain organohalide-respiring bacteria that use reductive dehalogenases (RDases) to dechlorinate CF to DCM. In contrast, known DCM degraders use DCM as their electron donor, which is oxidized using a series of methyltransferases and associated proteins encoded by the mec cassette to facilitate the entry of DCM to the Wood-Ljungdahl pathway. The SC05 culture is an enrichment culture sold commercially for bioaugmentation, which transforms CF via DCM to CO2. This culture has the unique ability to dechlorinate CF to DCM using electron equivalents provided by the oxidation of DCM to CO2. Here, we use metagenomic and metaproteomic analyses to identify the functional genes involved in each of these transformations. Though 91 metagenome-assembled genomes were assembled, the genes for an RDase-named acdA-and a complete mec cassette were found to be encoded on a single contig belonging to Dehalobacter. AcdA and critical Mec proteins were also highly expressed by the culture. Heterologously expressed AcdA dechlorinated CF and other chloroalkanes but had 100-fold lower activity on DCM. Overall, the high expression of Mec proteins and the activity of AcdA suggest a Dehalobacter capable of dechlorination of CF to DCM and subsequent mineralization of DCM using the mec cassette.

IMPORTANCE

Chloroform (CF) and dichloromethane (DCM) are regulated groundwater contaminants. A cost-effective approach to remove these pollutants from contaminated groundwater is to employ microbes that transform CF and DCM as part of their metabolism, thus depleting the contamination as the microbes continue to grow. In this work, we investigate bioaugmentation culture SC05, a mixed microbial consortium that effectively and simultaneously degrades both CF and DCM coupled to the growth of Dehalobacter. We identified the functional genes responsible for the transformation of CF and DCM in SC05. These genetic biomarkers provide a means to monitor the remediation process in the field.

The Streptococcus pyogenes stand-alone regulator RofA exhibits characteristics of a PRD-containing virulence regulator

Streptococcus pyogenes [group A streptococcus (GAS)] is a human pathogen capable of infecting diverse tissues. To successfully infect these sites, GAS must detect available nutrients and adapt accordingly. The phosphoenolpyruvate transferase system (PTS) mediates carbohydrate uptake and metabolic gene regulation to adapt to the nutritional environment. Regulation by the PTS can occur through phosphorylation of transcriptional regulators at conserved PTS-regulatory domains (PRDs). GAS has several PRD-containing stand-alone regulators with regulons encoding both metabolic genes and virulence factors [PRD-containing virulence regulators (PCVRs)]. One is RofA, which regulates the expression of virulence genes in multiple GAS serotypes. It was hypothesized that RofA is phosphorylated by the PTS in response to carbohydrate levels to coordinate virulence gene expression. In this study, the RofA regulon of M1T1 strain 5448 was determined using RNA sequencing. Two operons were consistently differentially expressed across growth in the absence of RofA; the pilus operon was downregulated, and the capsule operon was upregulated. This correlated with increased capsule production and decreased adherence to keratinocytes. Purified RofA-His was phosphorylated in vitro by PTS proteins EI and HPr, and phosphorylated RofA-FLAG was detected in vivo when GAS was grown in low-glucose C medium. Phosphorylated RofA was not observed when C medium was supplemented 10-fold with glucose. Mutations of select histidine residues within the putative PRDs contributed to the in vivo phosphorylation of RofA, although phosphorylation of RofA was still observed, suggesting other phosphorylation sites exist in the protein. Together, these findings support the hypothesis that RofA is a PCVR that may couple sugar metabolism with virulence regulation.

Diamide-based screening method for the isolation of improved oxidative stress tolerance phenotypes in Bacillus mutant libraries

IMPORTANCE

During their life cycle, bacteria are exposed to a range of different stresses that need to be managed appropriately in order to ensure their growth and viability. This applies not only to bacteria in their natural habitats but also to bacteria employed in biotechnological production processes. Oxidative stress is one of these stresses that may originate either from bacterial metabolism or external factors. In biotechnological settings, it is of critical importance that production strains are resistant to oxidative stresses. Accordingly, this also applies to the major industrial cell factory Bacillus subtilis. In the present study, we, therefore, developed a screen for B. subtilis strains with enhanced oxidative stress tolerance. The results show that our approach is feasible and time-, space-, and resource-efficient. We, therefore, anticipate that it will enhance the development of more robust industrial production strains with improved robustness under conditions of oxidative stress.

Classification of Promoter Sequences from Human Genome

We have developed a new method for promoter sequence classification based on a genetic algorithm and the MAHDS sequence alignment method. We have created four classes of human promoters, combining 17,310 sequences out of the 29,598 present in the EPD database. We searched the human genome for potential promoter sequences (PPSs) using dynamic programming and position weight matrices representing each of the promoter sequence classes. A total of 3,065,317 potential promoter sequences were found. Only 1,241,206 of them were located in unannotated parts of the human genome. Every other PPS found intersected with either true promoters, transposable elements, or interspersed repeats. We found a strong intersection between PPSs and Alu elements as well as transcript start sites. The number of false positive PPSs is estimated to be 3 × 10-8 per nucleotide, which is several orders of magnitude lower than for any other promoter prediction method. The developed method can be used to search for PPSs in various eukaryotic genomes.

SFPQ and Its Isoform as Potential Biomarker for Non-Small-Cell Lung Cancer

Cancer markers are measurable molecules in the blood or tissue that are produced by tumor cells or immune cells in response to cancer progression. They play an important role in clinical diagnosis, prognosis, and anti-drug monitoring. Although DNA, RNA, and even physical images have been used, proteins continue to be the most common marker. There are currently no specific markers for lung cancer. Metastatic lung cancer, particularly non-small-cell lung cancer (NSCLC), is one of the most common causes of death. SFPQ, YY1, RTN4, RICTOR, LARP6, and HELLS are expressed at higher levels in cells from NSCLC than in control or cells from inflammatory diseases. SFPQ shows the most difference between the three cell types. Furthermore, the cytoplasmic isoform of SFPQ is only found in advanced cancers. We have developed ELISAs to detect SFPQ and the long and short isoforms. Evidence has shown that the short isoform exists primarily in cancers. Furthermore, immunocytometry studies and IHC analysis have revealed that SFPQ levels are consistent with ELISA results. In addition, enhanced DNA methylation in the SFPQ gene may facilitate the SFPQ expression differences between control and cancer cells. Considering this, elevated SFPQ level and the isoform location could serve as a cancer diagnostic and prognostic marker.

From genotype to phenotype: computational approaches for inferring microbial traits relevant to the food industry

When selecting microbial strains for the production of fermented foods, various microbial phenotypes need to be taken into account to achieve target product characteristics, such as biosafety, flavor, texture, and health-promoting effects. Through continuous advances in sequencing technologies, microbial whole-genome sequences of increasing quality can now be obtained both cheaper and faster, which increases the relevance of genome-based characterization of microbial phenotypes. Prediction of microbial phenotypes from genome sequences makes it possible to quickly screen large strain collections in silico to identify candidates with desirable traits. Several microbial phenotypes relevant to the production of fermented foods can be predicted using knowledge-based approaches, leveraging our existing understanding of the genetic and molecular mechanisms underlying those phenotypes. In the absence of this knowledge, data-driven approaches can be applied to estimate genotype-phenotype relationships based on large experimental datasets. Here, we review computational methods that implement knowledge- and data-driven approaches for phenotype prediction, as well as methods that combine elements from both approaches. Furthermore, we provide examples of how these methods have been applied in industrial biotechnology, with special focus on the fermented food industry.

Computational identification of promoters in Klebsiella aerogenes by using support vector machine

Promoters are the basic functional cis-elements to which RNA polymerase binds to initiate the process of gene transcription. Comprehensive understanding gene expression and regulation depends on the precise identification of promoters, as they are the most important component of gene expression. This study aimed to develop a machine learning-based model to predict promoters in Klebsiella aerogenes (K. aerogenes). In the prediction model, the promoter sequences in K. aerogenes genome were encoded by pseudo k-tuple nucleotide composition (PseKNC) and position-correlation scoring function (PCSF). Numerical features were obtained and then optimized using mRMR by combining with support vector machine (SVM) and 5-fold cross-validation (CV). Subsequently, these optimized features were inputted into SVM-based classifier to discriminate promoter sequences from non-promoter sequences in K. aerogenes. Results of 10-fold CV showed that the model could yield the overall accuracy of 96.0% and the area under the ROC curve (AUC) of 0.990. We hope that this model will provide help for the study of promoter and gene regulation in K. aerogenes.

Computational prediction of promotors in Agrobacterium tumefaciens strain C58 by using the machine learning technique

Promotors are those genomic regions on the upstream of genes, which are bound by RNA polymerase for starting gene transcription. Because it is the most critical element of gene expression, the recognition of promoters is crucial to understand the regulation of gene expression. This study aimed to develop a machine learning-based model to predict promotors in Agrobacterium tumefaciens (A. tumefaciens) strain C58. In the model, promotor sequences were encoded by three different kinds of feature descriptors, namely, accumulated nucleotide frequency, k-mer nucleotide composition, and binary encodings. The obtained features were optimized by using correlation and the mRMR-based algorithm. These optimized features were inputted into a random forest (RF) classifier to discriminate promotor sequences from non-promotor sequences in A. tumefaciens strain C58. The examination of 10-fold cross-validation showed that the proposed model could yield an overall accuracy of 0.837. This model will provide help for the study of promoters in A. tumefaciens C58 strain.

Evidence of a Set of Core-Function Genes in 16 Bacillus Podoviral Genomes with Considerable Genomic Diversity

Bacteriophage genomes represent an enormous level of genetic diversity and provide considerable potential to acquire new insights about viral genome evolution. In this study, the genome sequences of sixteen Bacillus-infecting bacteriophages were explored through comparative genomics approaches to reveal shared and unique characteristics. These bacteriophages are in the Salasmaviridae family with small (18,548-27,206 bp) double-stranded DNA genomes encoding 25-46 predicted open reading frames. We observe extensive nucleotide and amino acid sequence divergence among a set of core-function genes that present clear synteny. We identify two examples of sequence directed recombination within essential genes, as well as explore the expansion of gene content in these genomes through the introduction of novel open reading frames. Together, these findings highlight the complex evolutionary relationships of phage genomes that include old, common origins as well as new components introduced through mosaicism.

Tuning of Gene Expression in Clostridium phytofermentans Using Synthetic Promoters and CRISPRi

Control of gene expression is fundamental to cell engineering. Here we demonstrate a set of approaches to tune gene expression in Clostridia using the model Clostridium phytofermentans. Initially, we develop a simple benchtop electroporation method that we use to identify a set of replicating plasmids and resistance markers that can be cotransformed into C. phytofermentans. We define a series of promoters spanning a >100-fold expression range by testing a promoter library driving the expression of a luminescent reporter. By insertion of tet operator sites upstream of the reporter, its expression can be quantitatively altered using the Tet repressor and anhydrotetracycline (aTc). We integrate these methods into an aTc-regulated dCas12a system with which we show in vivo CRISPRi-mediated repression of reporter and fermentation genes in C. phytofermentans. Together, these approaches advance genetic transformation and experimental control of gene expression in Clostridia.

Heterologous Gene Regulation in Clostridia: Rationally Designed Gene Regulation for Industrial and Medical Applications

Several species from the Clostridium genus show promise as industrial solvent producers and cancer therapeutic delivery vehicles. Previous development of shuttle plasmids and genome editing tools has aided the study of these species and enabled their exploitation in industrial and medical applications. Nevertheless, the precise control of gene expression is still hindered by the limited range of characterized promoters. To address this, libraries of promoters (native and synthetic), 5' UTRs, and alternative start codons were constructed. These constructs were tested in Escherichia coli K-12, Clostridium sporogenes NCIMB 10696, and Clostridium butyricum DSM 10702, using β-glucuronidase (gusA) as a gene reporter. Promoter activity was corroborated using a second gene reporter, nitroreductase (nmeNTR) from Neisseria meningitides. A strong correlation was observed between the two reporters. In C. sporogenes and C. butyricum, respectively, changes in GusA activity between the weakest and strongest expressing levels were 129-fold and 78-fold. Similar results were obtained with the nmeNTR. Using the GusA reporter, translation initiation from six alternative (non-AUG) start codons was measured in E. coli, C. sporogenes, and C. butyricum. Clearly, species-specific differences between clostridia and E. coli in translation initiation were observed, and the performance of the start codons was influenced by the upstream 5' UTR sequence. These results highlight a new opportunity for gene control in recombinant clostridia. To demonstrate the value of these results, expression of the sacB gene from Bacillus subtilis was optimized for use as a novel negative selection marker in C. butyricum. In summary, these results indicate improvements in the understanding of heterologous gene regulation in Clostridium species and E. coli cloning strains. This new knowledge can be utilized for rationally designed gene regulation in Clostridium-mediated industrial and medical applications, as well as fundamental research into the biology of Clostridium species.

Biochemical and Regulatory Analyses of Xylanolytic Regulons in Caldicellulosiruptor bescii Reveal Genus-Wide Features of Hemicellulose Utilization

Caldicellulosiruptor species scavenge carbohydrates from runoff containing plant biomass that enters hot springs and from grasses that grow in more moderate parts of thermal features. While only a few Caldicellulosiruptor species can degrade cellulose, all known species are hemicellulolytic. The most well-characterized species, Caldicellulosiruptor bescii, decentralizes its hemicellulase inventory across five different genomic loci and two isolated genes. Transcriptomic analyses, comparative genomics, and enzymatic characterization were utilized to assign functional roles and determine the relative importance of its six putative endoxylanases (five glycoside hydrolase family 10 [GH10] enzymes and one GH11 enzyme) and two putative exoxylanases (one GH39 and one GH3) in C. bescii. Two genus-wide conserved xylanases, C. bescii XynA (GH10) and C. bescii Xyl3A (GH3), had the highest levels of sugar release on oat spelt xylan, were in the top 10% of all genes transcribed by C. bescii, and were highly induced on xylan compared to cellulose. This indicates that a minimal set of enzymes are used to drive xylan degradation in the genus Caldicellulosiruptor, complemented by hemicellulolytic inventories that are tuned to specific forms of hemicellulose in available plant biomasses. To this point, synergism studies revealed that the pairing of specific GH family proteins (GH3, -11, and -39) with C. bescii GH10 proteins released more sugar in vitro than mixtures containing five different GH10 proteins. Overall, this work demonstrates the essential requirements for Caldicellulosiruptor to degrade various forms of xylan and the differences in species genomic inventories that are tuned for survival in unique biotopes with variable lignocellulosic substrates. IMPORTANCE Microbial deconstruction of lignocellulose for the production of biofuels and chemicals requires the hydrolysis of heterogeneous hemicelluloses to access the microcrystalline cellulose portion. This work extends previous in vivo and in vitro efforts to characterize hemicellulose utilization by integrating genomic reconstruction, transcriptomic data, operon structures, and biochemical characteristics of key enzymes to understand the deployment and functionality of hemicellulases by the extreme thermophile Caldicellulosiruptor bescii. Furthermore, comparative genomics of the genus revealed both conserved and divergent mechanisms for hemicellulose utilization across the 15 sequenced species, thereby paving the way to connecting functional enzyme characterization with metabolic engineering efforts to enhance lignocellulose conversion.

SPREAD: An ensemble predictor based on DNA autoencoder framework for discriminating promoters in Pseudomonas aeruginosa

Regulatory elements in DNA sequences, such as promoters, enhancers, terminators and so on, are essential for gene expression in physiological and pathological processes. A promoter is the specific DNA sequence that is located upstream of the coding gene and acts as the "switch" for gene transcriptional regulation. Lots of promoter predictors have been developed for different bacterial species, but only a few are designed for Pseudomonas aeruginosa, a widespread Gram-negative conditional pathogen in nature. In this work, an ensemble model named SPREAD is proposed for the recognition of promoters in Pseudomonas aeruginosa. In SPREAD, the DNA sequence autoencoder model LSTM is employed to extract potential sequence information, and the mean output probability value of CNN and RF is applied as the final prediction. Compared with G4PromFinder, the only state-of-the-art classifier for promoters in Pseudomonas aeruginosa, SPREAD improves the prediction performance significantly, with an accuracy of 0.98, recall of 0.98, precision of 0.98, specificity of 0.97 and F1-score of 0.98.

A pair of isoleucyl-tRNA synthetases in Bacilli fulfills complementary roles to keep fast translation and provide antibiotic resistance

Isoleucyl-tRNA synthetase (IleRS) is an essential enzyme that covalently couples isoleucine to the corresponding tRNA. Bacterial IleRSs group in two clades, ileS1 and ileS2, the latter bringing resistance to the natural antibiotic mupirocin. Generally, bacteria rely on either ileS1 or ileS2 as a standalone housekeeping gene. However, we have found an exception by noticing that Bacillus species with genomic ileS2 consistently also keep ileS1, which appears mandatory in the family Bacillaceae. Taking Priestia (Bacillus) megaterium as a model organism, we showed that PmIleRS1 is constitutively expressed, while PmIleRS2 is stress-induced. Both enzymes share the same level of the aminoacylation accuracy. Yet, PmIleRS1 exhibited a two-fold faster aminoacylation turnover (kcat ) than PmIleRS2 and permitted a notably faster cell-free translation. At the same time, PmIleRS2 displayed a 104 -fold increase in its Ki for mupirocin, arguing that the aminoacylation turnover in IleRS2 could have been traded-off for antibiotic resistance. As expected, a P. megaterium strain deleted for ileS2 was mupirocin-sensitive. Interestingly, an attempt to construct a mupirocin-resistant strain lacking ileS1, a solution not found among species of the family Bacillaceae in nature, led to a viable but compromised strain. Our data suggest that PmIleRS1 is kept to promote fast translation, whereas PmIleRS2 is maintained to provide antibiotic resistance when needed. This is consistent with an emerging picture in which fast-growing organisms predominantly use IleRS1 for competitive survival.

Database of Potential Promoter Sequences in the Capsicum annuum Genome

In this study, we used a mathematical method for the multiple alignment of highly divergent sequences (MAHDS) to create a database of potential promoter sequences (PPSs) in the Capsicum annuum genome. To search for PPSs, 20 statistically significant classes of sequences located in the range from −499 to +100 nucleotides near the annotated genes were calculated. For each class, a position–weight matrix (PWM) was computed and then used to identify PPSs in the C. annuum genome. In total, 825,136 PPSs were detected, with a false positive rate of 0.13%. The PPSs obtained with the MAHDS method were tested using TSSFinder, which detects transcription start sites. The databank of the found PPSs provides their coordinates in chromosomes, the alignment of each PPS with the PWM, and the level of statistical significance as a normal distribution argument, and can be used in genetic engineering and biotechnology.

Escherichia coli transcription factors of unknown function: sequence features and possible evolutionary relationships

Organisms need mechanisms to perceive the environment and respond accordingly to environmental changes or the presence of hazards. Transcription factors (TFs) are required for cells to respond to the environment by controlling the expression of genes needed. Escherichia coli has been the model bacterium for many decades, and still, there are features embedded in its genome that remain unstudied. To date, 58 TFs remain poorly characterized, although their binding sites have been experimentally determined. This study showed that these TFs have sequence variation at the third codon position G+C content but maintain the same Codon Adaptation Index (CAI) trend as annotated functional transcription factors. Most of these transcription factors are in areas of the genome where abundant repetitive and mobile elements are present. Sequence divergence points to groups with distinctive sequence signatures but maintaining the same type of DNA binding domain. Finally, the analysis of the promoter sequences of the 58 TFs showed A+T rich regions that agree with the features of horizontally transferred genes. The findings reported here pave the way for future research of these TFs that may uncover their role as spare factors in case of lose-of-function mutations in core TFs and trace back their evolutionary history.

Global Transcriptomic Analysis of Bacteriophage-Host Interactions between a Kayvirus Therapeutic Phage and Staphylococcus aureus

Kayviruses are polyvalent broad host range staphylococcal phages with a potential to combat staphylococcal infections. However, the implementation of rational phage therapy in medicine requires a thorough understanding of the interactions between bacteriophages and pathogens at omics level. To evaluate the effect of a phage used in therapy on its host bacterium, we performed differential transcriptomic analysis by RNA-Seq from bacteriophage K of genus Kayvirus infecting two Staphylococcus aureus strains, prophage-less strain SH1000 and quadruple lysogenic strain Newman. The temporal transcriptional profile of phage K was comparable in both strains except for a few loci encoding hypothetical proteins. Stranded sequencing revealed transcription of phage noncoding RNAs that may play a role in the regulation of phage and host gene expression. The transcriptional response of S. aureus to phage K infection resembles a general stress response with differential expression of genes involved in a DNA damage response. The host transcriptional changes involved upregulation of nucleotide, amino acid and energy synthesis and transporter genes and downregulation of host transcription factors. The interaction of phage K with variable genetic elements of the host showed slight upregulation of gene expression of prophage integrases and antirepressors. The virulence genes involved in adhesion and immune evasion were only marginally affected, making phage K suitable for therapy. IMPORTANCE Bacterium Staphylococcus aureus is a common human and veterinary pathogen that causes mild to life-threatening infections. As strains of S. aureus are becoming increasingly resistant to multiple antibiotics, the need to search for new therapeutics is urgent. A promising alternative to antibiotic treatment of staphylococcal infections is a phage therapy using lytic phages from the genus Kayvirus. Here, we present a comprehensive view on the phage-bacterium interactions on transcriptomic level that improves the knowledge of molecular mechanisms underlying the Kayvirus lytic action. The results will ensure safer usage of the phage therapeutics and may also serve as a basis for the development of new antibacterial strategies.

Precision Methylome and In Vivo Methylation Kinetics Characterization of Klebsiella pneumoniae

Klebsiella pneumoniae (K. pneumoniae) is an important pathogen that can cause severe hospital- and community-acquired infections. To systematically investigate its methylation features, we determined the whole-genome sequences of 14 K. pneumoniae strains covering varying serotypes, multilocus sequence types, clonal groups, viscosity/virulence, and drug resistance. Their methylomes were further characterized using Pacific Biosciences single-molecule real-time and bisulfite technologies. We identified 15 methylation motifs [13 N6-methyladenine (6mA) and two 5-methylcytosine (5mC) motifs], among which eight were novel. Their corresponding DNA methyltransferases were also validated. Additionally, we analyzed the genomic distribution of GATC and CCWGG methylation motifs shared by all strains, and identified differential distribution patterns of some hemi-/un-methylated GATC motifs, which tend to be located within intergenic regions (IGRs). Specifically, we characterized the in vivo methylation kinetics at single-base resolution on a genome-wide scale by simulating the dynamic processes of replication-mediated passive demethylation and MTase-catalyzed re-methylation. The slow methylation of the GATC motifs in the replication origin (oriC) regions and IGRs implicates the epigenetic regulation of replication initiation and transcription. Our findings illustrate the first comprehensive dynamic methylome map of K. pneumoniae at single-base resolution, and provide a useful reference to better understand epigenetic regulation in this and other bacterial species.

Characterisation of Bacteriophage vB_SmaM_Ps15 Infective to Stenotrophomonas maltophilia Clinical Ocular Isolates

Recent acknowledgment that multidrug resistant Stenotrophomonas maltophilia strains can cause severe infections has led to increasing global interest in addressing its pathogenicity. While being primarily associated with hospital-acquired respiratory tract infections, this bacterial species is also relevant to ophthalmology, particularly to contact lens-related diseases. In the current study, the capacity of Stenotrophomonas phage vB_SmaM_Ps15 to infect ocular S. maltophilia strains was investigated to explore its future potential as a phage therapeutic. The phage proved to be lytic to a range of clinical isolates collected in Australia from eye swabs, contact lenses and contact lens cases that had previously shown to be resistant to several antibiotics and multipurpose contact lenses disinfectant solutions. Morphological analysis by transmission electron microscopy placed the phage into the Myoviridae family. Its genome size was 161,350 bp with a G + C content of 54.2%, containing 276 putative protein-encoding genes and 24 tRNAs. A detailed comparative genomic analysis positioned vB_SmaM_Ps15 as a new species of the Menderavirus genus, which currently contains six very similar globally distributed members. It was confirmed as a virulent phage, free of known lysogenic and pathogenicity determinants, which supports its potential use for the treatment of S. maltophilia eye infections.

Promotech: a general tool for bacterial promoter recognition

Promoters are genomic regions where the transcription machinery binds to initiate the transcription of specific genes. Computational tools for identifying bacterial promoters have been around for decades. However, most of these tools were designed to recognize promoters in one or few bacterial species. Here, we present Promotech, a machine-learning-based method for promoter recognition in a wide range of bacterial species. We compare Promotech's performance with the performance of five other promoter prediction methods. Promotech outperforms these other programs in terms of area under the precision-recall curve (AUPRC) or precision at the same level of recall. Promotech is available at https://github.com/BioinformaticsLabAtMUN/PromoTech .

Development of a miRNA Sensor by an Inducible CRISPR-Cas9 Construct in Ciona Embryogenesis

MicroRNAs (miRNAs) regulate multicellular processes and diverse signaling pathways in organisms. The detection of the spatiotemporal expression of miRNA in vivo is crucial for uncovering the function of miRNA. However, most of the current detecting techniques cannot reflect the dynamics of miRNA sensitively in vivo. Here, we constructed a miRNA-induced CRISPR-Cas9 platform (MICR) used in marine chordate Ciona. The key component of MICR is a pre-single guide RNA (sgRNA) flanked by miRNA-binding sites that can be released by RNA-induced silencing complex (RISC) cleavage to form functional sgRNA in the presence of complementary miRNA. By using the miRNA-inducible CRISPR-on system (MICR-ON), we successfully detected the dynamic expression of a miRNA csa-miR-4018a during development of Ciona embryo. The detected patterns were validated to be consistent with the results by in situ hybridization. It is worth noting that the expression of csa-miR-4018a was examined by MICR-ON to be present in additional tissues, where no obvious signaling was detected by in situ hybridization, suggesting that the MICR-ON might be a more sensitive approach to detect miRNA signal in living animal. Thus, MICR-ON was demonstrated to be a sensitive and highly efficient approach for monitoring the dynamics of expression of miRNA in vivo and will facilitate the exploration of miRNA functions in biological systems.

Mathematical Algorithm for Identification of Eukaryotic Promoter Sequences

Identification of promoter sequences in the eukaryotic genome, by computer methods, is an important task of bioinformatics. However, this problem has not been solved since the best algorithms have a false positive probability of 10−3–10−4 per nucleotide. As a result of full genome analysis, there may be more false positives than annotated gene promoters. The probability of a false positive should be reduced to 10−6–10−8 to reduce the number of false positives and increase the reliability of the prediction. The method for multi alignment of the promoter sequences was developed. Then, mathematical methods were developed for calculation of the statistically important classes of the promoter sequences. Five promoter classes, from the rice genome, were created. We developed promoter classes to search for potential promoter sequences in the rice genome with a false positive number less than 10−8 per nucleotide. Five classes of promoter sequences contain 1740, 222, 199, 167 and 130 promoters, respectively. A total of 145,277 potential promoter sequences (PPSs) were identified. Of these, 18,563 are promoters of known genes, 87,233 PPSs intersect with transposable elements, and 37,390 PPSs were found in previously unannotated sequences. The number of false positives for a randomly mixed rice genome is less than 10−8 per nucleotide. The method developed for detecting PPSs was compared with some previously used approaches. The developed mathematical method can be used to search for genes, transposable elements, and transcript start sites in eukaryotic genomes.

Autorepressor PsrA is required for optimal Legionella pneumophila growth in Acanthamoeba castellanii protozoa

Legionella pneumophila possesses a unique intracellular lifecycle featuring distinct morphological stages that include replicative forms and transmissive cyst forms. Expression of genes associated with virulence traits and cyst morphogenesis is concomitant, and governed by a complex stringent response based-regulatory network and the stationary phase sigma factor RpoS. In Pseudomonas spp., rpoS expression is controlled by the autorepressor PsrA, and orthologs of PsrA and RpoS are required for cyst formation in Azotobacter. Here we report that the L. pneumophila psrA ortholog, expressed as a leaderless monocistronic transcript, is also an autorepressor, but is not a regulator of rpoS expression. Further, the binding site sequence recognized by L. pneumophila PsrA is different from that of Pseudomonas PsrA, suggesting a repertoire of target genes unique to L. pneumophila. While PsrA was dispensable for growth in human U937-derived macrophages, lack of PsrA affected bacterial intracellular growth in Acanthamoeba castellanii protozoa, but also increased the quantity of poly-3-hydroxybutyrate (PHB) inclusions in matured transmissive cysts. Interestingly, overexpression of PsrA increased the size and bacterial load of the replicative vacuole in both host cell types. Taken together, we report that PsrA is a host-specific requirement for optimal temporal progression of L. pneumophila intracellular lifecycle in A. castellanii.

Use of Mathematical Methods for the Biosafety Assessment of Agricultural Crops

In Russia and around the world, there are important questions regarding the potential threats to national and biological safety created by genetic technologies and the need to improve or introduce new, justified, and adequate measures for their control, regulation, and prevention. The article shows that a significant volume of the global market is occupied by five major transgenic crops, and producers are ready to switch to crops with an edited genome that has been approved in the United States, Argentina, and other countries. We propose a qualitatively new approach to the risk assessment of edited plants, "Safe Design," and we have also developed an extremely important, fundamentally new approach to the development of methods that combine next-generation sequencing (NGS) and Bioinformatics for the assessment of the crop import biosafety. The proposed mathematical approach provides a detailed analysis of the possible insertions of DNA fragments into the genome of edited crops and a clarification of their biological significance. The developed method can be used in the rapid screening of plants for the presence of potentially dangerous genes, viral sequences, and nonspecific promoter sequences.

Multiple Alignment of Promoter Sequences from the Arabidopsis thaliana L. Genome

In this study, we developed a new mathematical method for performing multiple alignment of highly divergent sequences (MAHDS), i.e., sequences that have on average more than 2.5 substitutions per position (x). We generated sets of artificial DNA sequences with x ranging from 0 to 4.4 and applied MAHDS as well as currently used multiple sequence alignment algorithms, including ClustalW, MAFFT, T-Coffee, Kalign, and Muscle to these sets. The results indicated that most of the existing methods could produce statistically significant alignments only for the sets with x < 2.5, whereas MAHDS could operate on sequences with x = 4.4. We also used MAHDS to analyze a set of promoter sequences from the Arabidopsis thaliana genome and discovered many conserved regions upstream of the transcription initiation site (from -499 to +1 bp); a part of the downstream region (from +1 to +70 bp) also significantly contributed to the obtained alignments. The possibilities of applying the newly developed method for the identification of promoter sequences in any genome are discussed. A server for multiple alignment of nucleotide sequences has been created.

Functional screening of a Caatinga goat (Capra hircus) rumen metagenomic library reveals a novel GH3 β-xylosidase

Functional screening of metagenomic libraries is an effective approach for identification of novel enzymes. A Caatinga biome goat rumen metagenomic library was screened using esculin as a substrate, and a gene from an unknown bacterium encoding a novel GH3 enzyme, BGL11, was identified. None of the BGL11 closely related genes have been previously characterized. Recombinant BGL11 was obtained and kinetically characterized. Substrate specificity of the purified protein was assessed using seven synthetic aryl substrates. Activity towards nitrophenyl-β-D-glucopyranoside (pNPG), 4-nitrophenyl-β-D-xylopyranoside (pNPX) and 4-nitrophenyl-β-D-cellobioside (pNPC) suggested that BGL11 is a multifunctional enzyme with β-glucosidase, β-xylosidase, and cellobiohydrolase activities. However, further testing with five natural substrates revealed that, although BGL11 has multiple substrate specificity, it is most active towards xylobiose. Thus, in its native goat rumen environment, BGL11 most likely functions as an extracellular β-xylosidase acting on hemicellulose. Biochemical characterization of BGL11 showed an optimal pH of 5.6, and an optimal temperature of 50°C. Enzyme stability, an important parameter for industrial application, was also investigated. At 40°C purified BGL11 remained active for more than 15 hours without reduction in activity, and at 50°C, after 7 hours of incubation, BGL11 remained 60% active. The enzyme kinetic parameters of Km and Vmax using xylobiose were determined to be 3.88 mM and 38.53 μmol.min-1.mg-1, respectively, and the Kcat was 57.79 s-1. In contrast to BLG11, most β-xylosidases kinetically studied belong to the GH43 family and have been characterized only using synthetic substrates. In industry, β-xylosidases can be used for plant biomass deconstruction, and the released sugars can be fermented into valuable bio-products, ranging from the biofuel ethanol to the sugar substitute xylitol.

A survey of non-coding RNAs in the social and predatory myxobacterium Myxococcus xanthus DK1622

Prokaryotic ncRNAs are important regulators of gene expression, and can be involved in complex signalling networks. The myxobacteria are model organisms for studies into multicellular development and microbial predation, being particularly renowned for their large genomes and exceptionally sophisticated signalling networks. However, apart from two specific examples, little is known about their regulatory ncRNAs. Here, we integrate bioinformatic predictions and transcriptome sequence data to provide a comprehensive survey of the ncRNAs made by the exemplar myxobacterium M. xanthus DK1622. M. xanthus RNA-seq data from four experimental conditions was interrogated to identify transcripts mapping outside coding sequences and to known ncRNAs. The resulting 37 ncRNAs were clustered on the genome and most (30/37) were conserved across the myxobacteria. A majority of ncRNAs (22/37) were intergenic, while 13 were at least partially antisense to protein-coding genes. Predicted promoter and terminator sequences explained the start/stop sites of 18 ncRNAs. mRNA targets for the ncRNAs were predicted, including plausible candidates for a known regulatory ncRNA. 22 ncRNAs were differentially expressed by nutrient availability and expression of 25 predicted targets was found to correlate strongly with that of their regulatory ncRNAs. Sharing of predicted mRNA targets by multiple ncRNAs suggests that some ncRNAs might regulate each other within signalling networks. This genomic survey of M. xanthus ncRNA biology provides a starting point for further studies of myxobacterial ncRNAs, which are likely to have important functions in these industrially important and sophisticated organisms.

A novel method SEProm for prokaryotic promoter prediction based on DNA structure and energetics

MOTIVATION

Despite conservation in general architecture of promoters and protein-DNA interaction interface of RNA polymerases among various prokaryotes, identification of promoter regions in the whole genome sequences remains a daunting challenge. The available tools for promoter prediction do not seem to address the problem satisfactorily, apparently because the biochemical nature of promoter signals is yet to be understood fully. Using 28 structural and 3 energetic parameters, we found that prokaryotic promoter regions have a unique structural and energy state, quite distinct from that of coding regions and the information for this signature state is in-built in their sequences. We developed a novel promoter prediction tool from these 31 parameters using various statistical techniques.

RESULTS

Here, we introduce SEProm, a novel tool that is developed by studying and utilizing the in-built structural and energy information of DNA sequences, which is applicable to all prokaryotes including archaea. Compared to five most recent, diverged and current best available tools, SEProm performs much better, predicting promoters with an 'F-value' of 82.04 and 'Precision' of 81.08. The next best 'F-value' was obtained with PromPredict (72.14) followed by BProm (68.37). On the basis of 'Precision' value, the next best 'Precision' was observed for Pepper (75.39) followed by PromPredict (72.01). SEProm maintained the lead even when comparison was done on two test organisms (not involved in training for SEProm).

AVAILABILITY AND IMPLEMENTATION

The software is freely available with easy to follow instructions (www.scfbio-iitd.res.in/software/TSS_Predict.jsp).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Differential Effects of Varying Concentrations of Phosphorus, Iron, and Nitrogen in N2-Fixing Cyanobacteria

Diazotrophs or N₂-fixers are one of the most ecologically significant groups in marine ecosystems (pelagic and benthic). Inorganic phosphorus (PO₄^3–) and iron (Fe) can limit the growth and N₂-fixing capacities of cyanobacteria. However, studies investigating co-limitation of these factors are lacking. Here, we added different concentrations of PO₄^3– and Fe in two cyanobacterial species whose relatives can be found in seagrass habitats: the unicellular Halothece sp. (PCC 7418) and the filamentous Fischerella muscicola (PCC 73103), grown under different nitrate (NO₃^–) concentrations and under N₂ as sole N source, respectively. Their growth, pigment content, N₂-fixation rates, oxidative stress responses, and morphological and cellular changes were investigated. Our results show a serial limitation of NO₃^– and PO₄^3– (with NO₃^– as the primary limiting nutrient) for Halothece sp. Simultaneous co-limitation of PO₄^3– and Fe was found for both species tested, and high levels of Fe (especially when added with high PO₄^3– levels) inhibited the growth of Halothece sp. Nutrient limitation (PO₄^3–, Fe, and/or NO₃^–) enhanced oxidative stress responses, morphological changes, and apoptosis. Furthermore, an extensive bio-informatic analysis describing the predicted Pho, Fur, and NtcA regulons (involved in the survival of cells to P, Fe, and N limitation) was made using the complete genome of Halothece sp. as a model, showing the potential of this strain to adapt to different nutrient regimes (P, Fe, or N).

Systematic Analysis of REBASE Identifies Numerous Type I Restriction-Modification Systems with Duplicated, Distinct hsdS Specificity Genes That Can Switch System Specificity by Recombination

Many bacterial species contain DNA methyltransferases that have random on/off switching of expression. These systems, called phasevarions (phase-variable regulons), control the expression of multiple genes by global methylation changes. In every previously characterized phasevarion, genes involved in pathobiology, antibiotic resistance, and potential vaccine candidates are randomly varied in their expression, commensurate with methyltransferase switching. Our systematic study to determine the extent of phasevarions controlled by invertible Type I R-M systems will provide valuable information for understanding how bacteria regulate genes and is key to the study of physiology, virulence, and vaccine development; therefore, it is critical to identify and characterize phase-variable methyltransferases controlling phasevarions.

N⁶-Adenine DNA methyltransferases associated with some Type I and Type III restriction-modification (R-M) systems are able to undergo phase variation, randomly switching expression ON or OFF by varying the length of locus-encoded simple sequence repeats (SSRs). This variation of methyltransferase expression results in genome-wide methylation differences and global changes in gene expression. These epigenetic regulatory systems are called phasevarions, phase-variable regulons, and are widespread in bacteria. A distinct switching system has also been described in Type I R-M systems, based on recombination-driven changes in hsdS genes, which dictate the DNA target site. In order to determine the prevalence of recombination-driven phasevarions, we generated a program called RecombinationRepeatSearch to interrogate REBASE and identify the presence and number of inverted repeats of hsdS downstream of Type I R-M loci. We report that 3.9% of Type I R-M systems have duplicated variable hsdS genes containing inverted repeats capable of phase variation. We report the presence of these systems in the major pathogens Enterococcus faecalis and Listeria monocytogenes, which could have important implications for pathogenesis and vaccine development. These data suggest that in addition to SSR-driven phasevarions, many bacteria have independently evolved phase-variable Type I R-M systems via recombination between multiple, variable hsdS genes.

IMPORTANCE Many bacterial species contain DNA methyltransferases that have random on/off switching of expression. These systems, called phasevarions (phase-variable regulons), control the expression of multiple genes by global methylation changes. In every previously characterized phasevarion, genes involved in pathobiology, antibiotic resistance, and potential vaccine candidates are randomly varied in their expression, commensurate with methyltransferase switching. Our systematic study to determine the extent of phasevarions controlled by invertible Type I R-M systems will provide valuable information for understanding how bacteria regulate genes and is key to the study of physiology, virulence, and vaccine development; therefore, it is critical to identify and characterize phase-variable methyltransferases controlling phasevarions.

Development of Lactococcus lactis Biosensors for Detection of Diacetyl

Some secondary metabolites of fermentative bacteria are desired compounds for the food industry. Examples of these compounds are diacetyl and acetaldehyde, which are produced by species of the lactic acid bacteria (LAB) family. Diacetyl is an aromatic compound, giving the buttery flavor associated with dairy products, and acetaldehyde is the compound responsible for the yogurt flavor and aroma. The quantification of these compounds in food matrices is a laborious task that involves sample preparation and specific analytical methods. The ability of bacteria to naturally sense metabolites has successfully been exploited to develop biosensors that facilitate the identification and quantification of certain metabolites (Mahr and Frunzke, 2016). The presence of a specific metabolite is sensed by the biosensors, and it is subsequently translated into the expression of one or more reporter genes. In this study we aimed to develop fluorescence-based biosensors to detect diacetyl and acetaldehyde. Since the metabolic pathways for production and degradation of these compounds are present in Lactococcus lactis, the sensing mechanisms in this bacterium are expected. Thus, we identified diacetyl and acetaldehyde responsive promoters by performing transcriptome analyses in L. lactis. The characterization of the biosensors showed their response to the presence of these compounds, and a further analysis of the diacetyl-biosensors (its dynamics and orthogonality) was performed. Moreover, we attempted to produce natural diacetyl from producer strains, namely L. lactis subsp. lactis biovar diacetylactis, to benchmark the performance of our biosensors. The diacetyl-biosensors responded linearly to the amounts of diacetyl obtained in the bacterial supernatants, i.e., the increases in GFP expression were proportional to the amounts of diacetyl present in the supernatants of L. lactis subsp. lactis biovar diacetylactis MR3-T7 strain. The biosensors developed in this study may eventually be used to engineer strains or pathways for increased diacetyl and acetaldehyde production, and may facilitate the detection of these metabolites in complex food matrices.

Siccibacter turicensis from Kangaroo Scats: Possible Implication in Cellulose Digestion

Microbiota in the kangaroo gut degrade cellulose, contributing to the kangaroo’s energy and survival. In this preliminary study, to discover more about the gut microbes that contribute to the survival of kangaroos, cellulose-degrading bacteria were isolated from kangaroo scats by selection on solidified media containing carboxymethyl cellulose as the main carbon source. One frequently occurring aerobic bacterium was Siccibacter turicensis, a microbe previously isolated in fruit powder and from a patient with angular cheilitis. The whole genome sequence of the kangaroo isolate was obtained using the Illumina MiSeq platform. Its sequence shared 97.98% identity of the S. turicensis Type strain, and the ability of the Type strain to degrade cellulose was confirmed. Analysis of the genomic data focused on the cellulose operon. In addition to genes from the operon, we suggest that a gene following the operon may have an important role in regulating cellulose metabolism by signal transduction. This is the first report of S. turicensis found within microbiota of the animal gut. Because of its frequent presence in the kangaroo gut, we suggest that S. turicensis plays a role in cellulose digestion for kangaroos.

Depletion of microbiome-derived molecules in the host using Clostridium genetics

The gut microbiota produce hundreds of molecules that are present at high concentrations in the host circulation. Unraveling the contribution of each molecule to host biology remains difficult. We developed a system for constructing clean deletions in Clostridium spp., the source of many molecules from the gut microbiome. By applying this method to the model commensal organism Clostridium sporogenes, we knocked out genes for 10 C. sporogenes-derived molecules that accumulate in host tissues. In mice colonized by a C. sporogenes for which the production of branched short-chain fatty acids was knocked out, we discovered that these microbial products have immunoglobulin A-modulatory activity.

Francisella novicida Two-Component System Response Regulator BfpR Modulates iglC Gene Expression, Antimicrobial Peptide Resistance, and Biofilm Production

Response regulators are a critical part of the two-component system of gene expression regulation in bacteria, transferring a signal from a sensor kinase into DNA binding activity resulting in alteration of gene expression. In this study, we investigated a previously uncharacterized response regulator in Francisella novicida, FTN_1452 that we have named BfpR (Biofilm-regulating Francisella protein Regulator, FTN_1452). In contrast to another Francisella response regulator, QseB/PmrA, BfpR appears to be a negative regulator of biofilm production, and also a positive regulator of antimicrobial peptide resistance in this bacterium. The protein was crystallized and X-ray crystallography studies produced a 1.8 Å structure of the BfpR N-terminal receiver domain revealing interesting insight into its potential interaction with the sensor kinase. Structural analysis of BfpR places it in the OmpR/PhoP family of bacterial response regulators along with WalR and ResD. Proteomic and transcriptomic analyses suggest that BfpR overexpression affects expression of the critical Francisella virulence factor iglC, as well as other proteins in the bacterium. We demonstrate that mutation of bfpR is associated with an antimicrobial peptide resistance phenotype, a phenotype also associated with other response regulators, for the human cathelicidin peptide LL-37 and a sheep antimicrobial peptide SMAP-29. F. novicida with mutated bfpR replicated better than WT in intracellular infection assays in human-derived macrophages suggesting that the down-regulation of iglC expression in bfpR mutant may enable this intracellular replication to occur. Response regulators have been shown to play important roles in the regulation of bacterial biofilm production. We demonstrate that F. novicida biofilm formation was highly increased in the bfpR mutant, corresponding to altered glycogen synthesis. Waxworm infection experiments suggest a role of BfpR as a negative modulator of iglC expression with de-repression by Mg²⁺. In this study, we find that the response regulator BfpR may be a negative regulator of biofilm formation, and a positive regulator of antimicrobial peptide resistance in F. novicida.

A riboswitch gives rise to multi-generational phenotypic heterogeneity in an auxotrophic bacterium

Auxotrophy, the inability to produce an organic compound essential for growth, is widespread among bacteria. Auxotrophic bacteria rely on transporters to acquire these compounds from their environment. Here, we study the expression of both low- and high-affinity transporters of the costly amino acid methionine in an auxotrophic lactic acid bacterium, Lactococcus lactis. We show that the high-affinity transporter (Met-transporter) is heterogeneously expressed at low methionine concentrations, resulting in two isogenic subpopulations that sequester methionine in different ways: one subpopulation primarily relies on the high-affinity transporter (high expression of the Met-transporter) and the other subpopulation primarily relies on the low-affinity transporter (low expression of the Met-transporter). The phenotypic heterogeneity is remarkably stable, inherited for tens of generations, and apparent at the colony level. This heterogeneity results from a T-box riboswitch in the promoter region of the met operon encoding the high-affinity Met-transporter. We hypothesize that T-box riboswitches, which are commonly found in the Lactobacillales, may play as-yet unexplored roles in the predominantly auxotrophic lifestyle of these bacteria.

Genome-scale exploration of transcriptional regulation in the nisin Z producer Lactococcus lactis subsp. lactis IO-1

Transcription is of the most crucial steps of gene expression in bacteria, whose regulation guarantees the bacteria's ability to adapt to varying environmental conditions. Discovering the molecular basis and genomic principles of the transcriptional regulation is thus one of the most important tasks in cellular and molecular biology. Here, a comprehensive phylogenetic footprinting framework was implemented to predict maximal regulons of Lactococcus lactis subsp. lactis IO-1, a lactic acid bacterium known for its high potentials in nisin Z production as well as efficient xylose consumption which have made it a promising biotechnological strain. A total set of 321 regulons covering more than 90% of all the bacterium's operons have been elucidated and validated according to available data. Multiple novel biologically-relevant members were introduced amongst which arsC, mtlA and mtl operon for BusR, MtlR and XylR regulons can be named, respectively. Moreover, the effect of riboflavin on nisin biosynthesis was assessed in vitro and a negative correlation was observed. It is believed that understandings from such networks not only can be useful for studying transcriptional regulatory potentials of the target organism but also can be implemented in biotechnology to rationally design favorable production conditions.

Homo- and heterodimerization is a common feature of the solute carrier family SLC10 members

The solute carrier family SLC10 consists of seven members, including the bile acid transporters Na+/taurocholate co-transporting polypeptide (NTCP) and apical sodium-dependent bile acid transporter (ASBT), the steroid sulfate transporter SOAT as well as four orphan carriers (SLC10A3, SLC10A4, SLC10A5 and SLC10A7). Previously, homodimerization of NTCP, ASBT and SOAT was described and there is increasing evidence that carrier oligomerization is an important regulatory factor for protein sorting and transport function. In the present study, homo- and heterodimerization were systematically analyzed among all SLC10 carriers (except for SLC10A3) using the yeast-two-hybrid membrane protein system. Strong homodimerization occurred for NTCP/NTCP, ASBT/ASBT and SLC10A7/SLC10A7. Heterodimerization was observed for most of the SLC10 carrier combinations. Heterodimerization of NTCP was additionally investigated by co-localization of NTCP-GFP and NTCP-mScarlet with respective SLC10 carrier constructs. NTCP co-localized with SLC10A4, SLC10A5, SOAT and SLC10A7. This co-localization was most pronounced for SLC10A4 and was additionally confirmed by co-immunoprecipitation. Interestingly, SLC10 carrier co-expression decreased the taurocholate transport function of NTCP for most of the analyzed constructs, indicating that SLC10 carrier heterodimerization is of functional relevance. In conclusion, homo- and heterodimerization is a common feature of the SLC10 carriers. The relevance of this finding for regulation and transport function of the SLC10 carriers in vivo needs further investigation.

Characterization and diversity of the complete set of GH family 3 enzymes from Rhodothermus marinus DSM 4253

The genome of Rhodothermus marinus DSM 4253 encodes six glycoside hydrolases (GH) classified under GH family 3 (GH3): RmBgl3A, RmBgl3B, RmBgl3C, RmXyl3A, RmXyl3B and RmNag3. The biochemical function, modelled 3D-structure, gene cluster and evolutionary relationships of each of these enzymes were studied. The six enzymes were clustered into three major evolutionary lineages of GH3: β-N-acetyl-glucosaminidases, β-1,4-glucosidases/β-xylosidases and macrolide β-glucosidases. The RmNag3 with additional β-lactamase domain clustered with the deepest rooted GH3-lineage of β-N-acetyl-glucosaminidases and was active on acetyl-chitooligosaccharides. RmBgl3B displayed β-1,4-glucosidase activity and was the only representative of the lineage clustered with macrolide β-glucosidases from Actinomycetes. The β-xylosidases, RmXyl3A and RmXyl3B, and the β-glucosidases RmBgl3A and RmBgl3C clustered within the major β-glucosidases/β-xylosidases evolutionary lineage. RmXyl3A and RmXyl3B showed β-xylosidase activity with different specificities for para-nitrophenyl (pNP)-linked substrates and xylooligosaccharides. RmBgl3A displayed β-1,4-glucosidase/β-xylosidase activity while RmBgl3C was active on pNP-β-Glc and β-1,3-1,4-linked glucosyl disaccharides. Putative polysaccharide utilization gene clusters were also investigated for both R. marinus DSM 4253 and DSM 4252^T (homolog strain). The analysis showed that in the homolog strain DSM 4252^TRmar_1080 (RmXyl3A) and Rmar_1081 (RmXyl3B) are parts of a putative polysaccharide utilization locus (PUL) for xylan utilization.

Computational Tools for Discovering and Engineering Natural Product Biosynthetic Pathways

Natural products (NPs), also known as secondary metabolites, are produced in bacteria, fungi, and plants. NPs represent a rich source of antibacterial, antifungal, and anticancer agents. Recent advances in DNA sequencing technologies and bioinformatics unveiled nature's great potential for synthesizing numerous NPs that may confer unprecedented structural and biological features. However, discovering novel bioactive NPs by genome mining remains a challenge. Moreover, even with interesting bioactivity, the low productivity of many NPs significantly limits their practical applications. Here we discuss the progress in developing bioinformatics tools for efficient discovery of bioactive NPs. In addition, we highlight computational methods for optimizing the productivity of NPs of pharmaceutical importance.

Interspecies Inhibition of Porphyromonas gingivalis by Yogurt-Derived Lactobacillus delbrueckii Requires Active Pyruvate Oxidase

Despite a growing interest in using probiotic microorganisms to prevent disease, the mechanisms by which probiotics exert their action require further investigation. Porphyromonas gingivalis is an important pathogen implicated in the development of periodontitis. We isolated several strains of Lactobacillus delbrueckii from dairy products and examined their ability to inhibit P. gingivalis growth in vitro We observed strain-specific inhibition of P. gingivalis growth in vitro Whole-genome sequencing of inhibitory and noninhibitory strains of L. delbrueckii revealed significant genetic differences supporting the strain specificity of the interaction. Extracts of the L. delbrueckii STYM1 inhibitory strain contain inhibitory activity that is abolished by treatment with heat, proteinase K, catalase, and sodium sulfite. We purified the inhibitory protein(s) from L. delbrueckii STYM1 extracts using ammonium sulfate precipitation, anion-exchange chromatography, and gel filtration chromatography. Pyruvate oxidase was highly enriched in the purified samples. Lastly, we showed that purified, catalytically active, recombinant pyruvate oxidase is sufficient to inhibit P. gingivalis growth in vitro without the addition of cofactors. Further, using a saturated transposon library, we isolated transposon mutants of P. gingivalis in the feoB2 (PG_1294) gene that are resistant to killing by inhibitory L. delbrueckii, consistent with a mechanism of hydrogen peroxide production by pyruvate oxidase. Our results support the current understanding of the importance of strain selection, not simply species selection, in microbial interactions. Specific L. delbrueckii strains or their products may be effective in the treatment and prevention of P. gingivalis-associated periodontal disease.IMPORTANCE P. gingivalis is implicated in the onset and progression of periodontal disease and associated with some systemic diseases. Probiotic bacteria represent an attractive preventative therapy for periodontal disease. However, the efficacy of probiotic bacteria can be variable between studies. Our data support the known importance of selecting particular strains of bacteria for probiotic use, not simply a single species. Specifically, in the context of probiotic intervention of periodontitis, our data suggest that high-level expression of pyruvate oxidase with hydrogen peroxide production in L. delbrueckii could be an important characteristic for the design of a probiotic supplement or a microbial therapeutic.

Biological and molecular characterization of a bacteriophage infecting Xanthomonas campestris pv. campestris, isolated from brassica fields

Xanthomonas campestris pv. campestris (Xcc) is the causal agent of black rot of crucifers. Here, we report a virus that infects Xcc isolated from brassica fields in Brazil. Morphological, molecular and phylogenetic analysis indicated that the isolated virus is a new member of the genus Pbunavirus, family Myoviridae, and we propose the name "Xanthomonas virus XC 2" for this virus. The isolated virus has a narrow host range, infecting only Xcc isolates, and it did not infect unrelated bacteria. These results indicate that the isolated bacteriophage is highly specific for Xcc and may be a potential agent for biological control.

From biophysics to 'omics and systems biology

Recent decades brought a revolution to biology, driven mainly by exponentially increasing amounts of data coming from "'omics" sciences. To handle these data, bioinformatics often has to combine biologically heterogeneous signals, for which methods from statistics and engineering (e.g. machine learning) are often used. While such an approach is sometimes necessary, it effectively treats the underlying biological processes as a black box. Similarly, systems biology deals with inherently complex systems, characterized by a large number of degrees of freedom, and interactions that are highly non-linear. To deal with this complexity, the underlying physical interactions are often (over)simplified, such as in Boolean modelling of network dynamics. In this review, we argue for the utility of applying a biophysical approach in bioinformatics and systems biology, including discussion of two examples from our research which address sequence analysis and understanding intracellular gene expression dynamics.