Bacillus velezensis strains are applied as ecologically safe biopesticides, plant growth promoting rhizobacteria (PGPR), and in veterinary probiotics. They are abundant in various environments including soil, plants, marine habitats, the intestinal micro-flora, etc. The mechanisms underlying this adaptive plasticity and bioactivity are not well understood, nor is it clear why several strains outperform other same species isolates by their bioactivities. The main objective of this work was to demonstrate versatility of bioactivities and lifestyle strategies of the selected B. velezensis strains suitable to serve as model organisms in future studies. Here, we performed a comparative study of newly sequenced genomes of four B. velezensis isolates with distinct phenotypes and isolation origin, which were assessed by RNA sequencing under the effect of root exudate stimuli and profiled by epigenetic modifications of chromosomal DNA. Among the selected strains, UCMB5044 is an oligotrophic PGPR strain adapted to nutrient poor desert soils. UCMB5113 and At1 are endophytes that colonize plants and require nutrient rich media. In contrast, the probiotic strain, UCMB5007, is a copiotroph, which shows no propensity to colonize plants. PacBio and Illumina sequencing approaches were used to generate complete genome assemblies, tracing epigenetic modifications, and determine gene expression profiles. All sequence data was deposited at NCBI. The strains, UCMB5113 and At1, show 99% sequence identity and similar phenotypes despite being isolated from geographically distant regions. UCMB5007 and UCMB5044 represent another group of organisms with almost identical genomes but dissimilar phenotypes and plant colonization propensity. The two plant associated strains, UCMB5044 and UCMB5113, share 398 genes putatively associated with root colonization, which are activated by exposure to maize root exudates. In contrast, UCMB5007 did not respond to root exudate stimuli. It was hypothesized that alterations in the global methylation pattern and some other epigenetic modifications enable adaptation of strains to different habitats and therefore may be of importance in terms of the biotechnological applicability of these bacteria. Contrary, the ability to grow on root exudates as a sole source of nutrients or a strong antagonism against phytopathogens showed by the strains in vitro cannot be considered as good predictors of PGPR activities.

A chronic inflammatory state to a large extent explains sickle cell disease (SCD) pathophysiology. Nonetheless, the principal dysregulated factors affecting this major pathway and their mechanisms of action still have to be fully identified and elucidated. Integrating gene expression and genome-wide association study (GWAS) data analysis represents a novel approach to refining the identification of key mediators and functions in complex diseases. Here, we performed gene expression meta-analysis of five independent publicly available microarray datasets related to homozygous SS patients with SCD to identify a consensus SCD transcriptomic profile. The meta-analysis conducted using the MetaDE R package based on combining p values (maxP approach) identified 335 differentially expressed genes (DEGs; 224 upregulated and 111 downregulated). Functional gene set enrichment revealed the importance of several metabolic pathways, of innate immune responses, erythrocyte development, and hemostasis pathways. Advanced analyses of GWAS data generated within the framework of this study by means of the atSNP R package and SIFT tool identified 60 regulatory single-nucleotide polymorphisms (rSNPs) occurring in the promoter of 20 DEGs and a deleterious SNP, affecting CAMKK2 protein function. This novel database of candidate genes, transcription factors, and rSNPs associated with SCD provides new markers that may help to identify new therapeutic targets.

Background: Sickle cell disease (SCD) is a blood disorder caused by a point mutation on the beta globin gene resulting in the synthesis of abnormal hemoglobin. Fetal hemoglobin (HbF) reduces disease severity, but the levels vary from one individual to another. Most research has focused on common genetic variants which differ across populations and hence do not fully account for HbF variation. Methods: We investigated rare and common genetic variants that influence HbF levels in 14 SCD patients to elucidate variants and pathways in SCD patients with extreme HbF levels (≥7.7% for high HbF) and (≤2.5% for low HbF) in Tanzania. We performed targeted next generation sequencing (Illumina_Miseq) covering exonic and other significant fetal hemoglobin-associated loci, including BCL11A, MYB, HOXA9, HBB, HBG1, HBG2, CHD4, KLF1, MBD3, ZBTB7A and PGLYRP1. Results: Results revealed a range of genetic variants, including bi-allelic and multi-allelic SNPs, frameshift insertions and deletions, some of which have functional importance. Notably, there were significantly more deletions in individuals with high HbF levels (11% vs 0.9%). We identified frameshift deletions in individuals with high HbF levels and frameshift insertions in individuals with low HbF. CHD4 and MBD3 genes, interacting in the same sub-network, were identified to have a significant number of pathogenic or non-synonymous mutations in individuals with low HbF levels, suggesting an important role of epigenetic pathways in the regulation of HbF synthesis. Conclusions: This study provides new insights in selecting essential variants and identifying potential biological pathways associated with extreme HbF levels in SCD interrogating multiple genomic variants associated with HbF in SCD.

Despite successful use of Plant Growth Promoting Rhizobacteria (PGPR) in agriculture, little is known about specific mechanisms of gene regulation facilitating the effective communication between bacteria and plants during plant colonization. Active PGPR strain B. atrophaeus UCMB-5137 was studied in this research. RNA sequencing profiles were generated in experiments where root exudate stimulations were used to mimic interactions between bacteria and plants. It was found that the gene regulation in B. atrophaeus UCMB-5137 in response to the root exudate stimuli differed from the reported gene regulation at similar conditions in B. amyloliquefaciens FZB42, which was considered as a paradigm PGPR. This difference was explained by hypersensitivity of UCMB-5137 to the root exudate stimuli impelling it to a sessile root colonization behavior through the CcpA-CodY-AbrB regulation. It was found that the transcriptional factor DegU also could play an important role in gene regulations during plant colonization. A significant stress caused by the root exudates on in vitro cultivated B. atrophaeus UCMB-5137 was noticed and discussed. Multiple cases of conflicted gene regulations showed scantiness of our knowledge on the regulatory network in Bacillus. Some of these conflicted regulations could be explained by interference of non-coding RNA (ncRNA). Search through differential expressed intergenic regions revealed 49 putative loci of ncRNA regulated by the root exudate stimuli. Possible target mRNA were predicted and a general regulatory network of B. atrophaeus UCMB-5137 genome was designed.1 Highlights  Plant colonizing PGPR Bacillus responded differently to the root exudate stimuli;  In UCMB-5137 the CcpA-CodY-AbrB regulation caused fast cell immobilization;  DegU regulon is important for plant colonization behavior in PGPR Bacillus;  ncRNA involved in regulation of plant colonization were identified;  A comprehensive model of gene regulation in UCMB-5137 was developed;