DNA structure at the plasmid origin-of-transfer indicates its potential transfer range

Multiple relaxases contribute to the horizontal transfer of the virulence plasmids from the tumorigenic bacterium Pseudomonas syringae pv. savastanoi NCPPB 3335

Pseudomonas syringae pv. savastanoi NCPPB 3335 is the causal agent of olive knot disease and contains three virulence plasmids: pPsv48A (pA), 80 kb; pPsv48B (pB), 45 kb, and pPsv48C (pC), 42 kb. Here we show that pB contains a complete MPF_T (previously type IVA secretion system) and a functional origin of conjugational transfer adjacent to a relaxase of the MOB_P family; pC also contains a functional oriT-MOB_P array, whereas pA contains an incomplete MPF_I (previously type IVB secretion system), but not a recognizable oriT. Plasmid transfer occurred on solid and in liquid media, and on leaf surfaces of a non-host plant (Phaseolus vulgaris) with high (pB) or moderate frequency (pC); pA was transferred only occasionally after cointegration with pB. We found three plasmid-borne and three chromosomal relaxase genes, although the chromosomal relaxases did not contribute to plasmid dissemination. The MOB_P relaxase genes of pB and pC were functionally interchangeable, although with differing efficiencies. We also identified a functional MOB_Q mobilization region in pC, which could only mobilize this plasmid. Plasmid pB could be efficiently transferred to strains of six phylogroups of P. syringae sensu lato, whereas pC could only be mobilized to two strains of phylogroup 3 (genomospecies 2). In two of the recipient strains, pB was stably maintained after 21 subcultures in liquid medium. The carriage of several relaxases by the native plasmids of P. syringae impacts their transfer frequency and, by providing functional diversity and redundancy, adds robustness to the conjugation system.

Insights into the ancestry evolution of the Mycobacterium tuberculosis complex from analysis of Mycobacterium riyadhense

Current evolutionary scenarios posit the emergence of Mycobacterium tuberculosis from an environmental saprophyte through a cumulative process of genome adaptation. Mycobacterium riyadhense, a related bacillus, is being increasingly isolated from human clinical cases with tuberculosis-like symptoms in various parts of the world. To elucidate the evolutionary relationship between M. riyadhense and other mycobacterial species, including members of the M. tuberculosis complex (MTBC), eight clinical isolates of M. riyadhense were sequenced and analyzed. We show, among other features, that M. riyadhense shares a large number of conserved orthologs with M. tuberculosis and shows the expansion of toxin/antitoxin pairs, PE/PPE family proteins compared with other non-tuberculous mycobacteria. We observed M. riyadhense lacks wecE gene which may result in the absence of lipooligosaccharides (LOS) IV. Comparative transcriptomic analysis of infected macrophages reveals genes encoding inducers of Type I IFN responses, such as cytosolic DNA sensors, were relatively less expressed by macrophages infected with M. riyadhense or M. kansasii, compared to BCG or M. tuberculosis. Overall, our work sheds new light on the evolution of M. riyadhense, its relationship to the MTBC, and its potential as a system for the study of mycobacterial virulence and pathogenesis.

Performance of Regression Models as a Function of Experiment Noise

Background:

A challenge in developing machine learning regression models is that it is difficult to know whether maximal performance has been reached on the test dataset, or whether further model improvement is possible. In biology, this problem is particularly pronounced as sample labels (response variables) are typically obtained through experiments and therefore have experiment noise associated with them. Such label noise puts a fundamental limit to the metrics of performance attainable by regression models on the test dataset.

Results:

We address this challenge by deriving an expected upper bound for the coefficient of determination (R²) for regression models when tested on the holdout dataset. This upper bound depends only on the noise associated with the response variable in a dataset as well as its variance. The upper bound estimate was validated via Monte Carlo simulations and then used as a tool to bootstrap performance of regression models trained on biological datasets, including protein sequence data, transcriptomic data, and genomic data.

Conclusions:

The new method for estimating upper bounds for model performance on test data should aid researchers in developing ML regression models that reach their maximum potential. Although we study biological datasets in this work, the new upper bound estimates will hold true for regression models from any research field or application area where response variables have associated noise.

Learning the Regulatory Code of Gene Expression

Data-driven machine learning is the method of choice for predicting molecular phenotypes from nucleotide sequence, modeling gene expression events including protein-DNA binding, chromatin states as well as mRNA and protein levels. Deep neural networks automatically learn informative sequence representations and interpreting them enables us to improve our understanding of the regulatory code governing gene expression. Here, we review the latest developments that apply shallow or deep learning to quantify molecular phenotypes and decode the cis-regulatory grammar from prokaryotic and eukaryotic sequencing data. Our approach is to build from the ground up, first focusing on the initiating protein-DNA interactions, then specific coding and non-coding regions, and finally on advances that combine multiple parts of the gene and mRNA regulatory structures, achieving unprecedented performance. We thus provide a quantitative view of gene expression regulation from nucleotide sequence, concluding with an information-centric overview of the central dogma of molecular biology.

Swine growth promotion with antibiotics or alternatives can increase antibiotic resistance gene mobility potential

Even though the use of antibiotics for food-producing animals may contribute to the emergence of antimicrobial resistance, antibiotics are still used as growth promoters. Due to consumer and regulatory pressures, the use of alternatives to antibiotics as growth promoters is increasing, thus more information is needed on their capability to disseminate antimicrobial resistance compared to antibiotics. We investigated the impacts of carbadox (antibiotic), copper sulfate and zinc oxide (metals) and mushroom powder (natural product) on the pig fecal resistome and microbiome. Antibiotic resistance gene (ARG) and mobile genetic element (MGE) abundances were measured using a high-throughput qPCR array with 382 primer pairs. Bacterial community composition was determined by 16S rRNA gene sequencing. More ARGs co-occurred with MGEs in the growth promoter group samples than in the control group samples. Community composition could not be linked to resistome in the growth promoter group samples, indicating a potential decoupling of ARGs and phylogeny. Additionally, machine-learning methods aided in defining the community and resistome differences in response to treatments. Since increased ARG mobility potential was the primary response to the dietary additives used in this study, we suggest that ARG mobility should be considered when designing antimicrobial use policies and antimicrobial resistance surveillances.

Deep learning suggests that gene expression is encoded in all parts of a co-evolving interacting gene regulatory structure

Understanding the genetic regulatory code governing gene expression is an important challenge in molecular biology. However, how individual coding and non-coding regions of the gene regulatory structure interact and contribute to mRNA expression levels remains unclear. Here we apply deep learning on over 20,000 mRNA datasets to examine the genetic regulatory code controlling mRNA abundance in 7 model organisms ranging from bacteria to Human. In all organisms, we can predict mRNA abundance directly from DNA sequence, with up to 82% of the variation of transcript levels encoded in the gene regulatory structure. By searching for DNA regulatory motifs across the gene regulatory structure, we discover that motif interactions could explain the whole dynamic range of mRNA levels. Co-evolution across coding and non-coding regions suggests that it is not single motifs or regions, but the entire gene regulatory structure and specific combination of regulatory elements that define gene expression levels.

Multiple plasmid origin-of-transfer regions might aid the spread of antimicrobial resistance to human pathogens

Antimicrobial resistance poses a great danger to humanity, in part due to the widespread horizontal gene transfer of plasmids via conjugation. Modeling of plasmid transfer is essential to uncovering the fundamentals of resistance transfer and for the development of predictive measures to limit the spread of resistance. However, a major limitation in the current understanding of plasmids is the incomplete characterization of the conjugative DNA transfer mechanisms, which conceals the actual potential for plasmid transfer in nature. Here, we consider that the plasmid-borne origin-of-transfer substrates encode specific DNA structural properties that can facilitate finding these regions in large datasets and develop a DNA structure-based alignment procedure for typing the transfer substrates that outperforms sequence-based approaches. Thousands of putative DNA transfer substrates are identified, showing that plasmid mobility can be twofold higher and span almost twofold more host species than is currently known. Over half of all putative mobile plasmids contain the means for mobilization by conjugation systems belonging to different mobility groups, which can hypothetically link previously confined host ranges across ecological habitats into a robust plasmid transfer network. This hypothetical network is found to facilitate the transfer of antimicrobial resistance from environmental genetic reservoirs to human pathogens, which might be an important driver of the observed rapid resistance development in humans and thus an important point of focus for future prevention measures.

Plasmid Design for Tunable Two-Enzyme Co-Expression Promotes Whole-Cell Production of Cellobiose

Catalyst development for biochemical cascade reactions often follows a "whole-cell-approach" in which a single microbial cell is made to express all required enzyme activities. Although attractive in principle, the approach can encounter limitations when efficient overall flux necessitates precise balancing between activities. This study shows an effective integration of major design strategies from synthetic biology to a coherent development of plasmid vectors, enabling tunable two-enzyme co-expression in E. coli, for whole-cell-production of cellobiose. An efficient transformation of sucrose and glucose into cellobiose by a parallel (countercurrent) cascade of disaccharide phosphorylases requires the enzyme co-expression to cope with large differences in specific activity of cellobiose phosphorylase (14 U mg^-1 ) and sucrose phosphorylase (122 U mg^-1 ). Mono- and bicistronic co-expression strategies controlling transcription, transcription-translation coupling or plasmid replication are analyzed for effect on activity and stable producibility of the whole-cell-catalyst. A key role of bom (basis of mobility) for plasmid stability dependent on the ori is reported and the importance of RBS (ribosome binding site) strength is demonstrated. Whole cell catalysts show high specific rates (460 µmol cellobiose min^-1  g^-1 dry cells) and performance metrics (30 g L^-1 ; ∼82% yield; 3.8 g L^-1 h^-1 overall productivity) promising for cellobiose production.

Understanding Horizontal Gene Transfer network in human gut microbiota

Background

Horizontal Gene Transfer (HGT) is the process of transferring genetic materials between species. Through sharing genetic materials, microorganisms in the human microbiota form a network. The network can provide insights into understanding the microbiota. Here, we constructed the HGT networks from the gut microbiota sequencing data and performed network analysis to characterize the HGT networks of gut microbiota.

Results

We constructed the HGT network and perform the network analysis to two typical gut microbiota datasets, a 283-sample dataset of Mother-to-Child and a 148-sample dataset of longitudinal inflammatory bowel disease (IBD) metagenome. The results indicated that (1) the HGT networks are scale-free. (2) The networks expand their complexities, sizes, and edge numbers, accompanying the early stage of lives; and microbiota established in children shared high similarity as their mother (p-value = 0.0138), supporting the transmission of microbiota from mother to child. (3) Groups harbor group-specific network edges, and network communities, which can potentially serve as biomarkers. For instances, IBD patient group harbors highly abundant communities of Proteobacteria (p-value = 0.0194) and Actinobacteria (p-value = 0.0316); children host highly abundant communities of Proteobacteria (p-value = 2.8785\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^{-5}$$\end{document}e-5) and Actinobacteria (p-value = 0.0015), and the mothers host highly abundant communities of Firmicutes (p-value = 8.0091\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^{-7}$$\end{document}e-7). IBD patient networks contain more HGT edges in pathogenic genus, including Mycobacterium, Sutterella, and Pseudomonas. Children’s networks contain more edges from Bifidobacterium and Escherichia.

Conclusion

Hence, we proposed the HGT network constructions from the gut microbiota sequencing data. The HGT networks capture the host state and the response of microbiota to the environmental and host changes, and they are essential to understand the human microbiota.

DNA processing by the MOBH family relaxase TraI encoded within the gonococcal genetic island

Relaxases of the MOB_H family are often found on large plasmids, genetic islands and integrative conjugative elements. Many members of this family contain an N-terminal relaxase domain (TraI_2) followed by a disordered middle part and a C-terminal domain of unknown function (TraI_2_C). The TraI_2 domain contains two putative metal-binding motifs, an HD domain motif and an alternative 3H motif. TraI, encoded within the gonococcal genetic island of Neisseria gonorrhoeae, is the prototype of the MOB_H family. SAXS experiments showed that TraI_2 and TraI_2_C form globular structures separated by an extended middle domain. The TraI_2 domain cleaves oriT-ssDNA in a site-specific Mn²⁺ or Co²⁺ dependent manner. The minimal oriT encompasses 50 nucleotides, requires an inverted repeat 3′ of the nic-site and several nucleotides around nic for efficient cleavage. Surprisingly, no stable covalent relaxase-DNA intermediate was observed. Mutagenesis of conserved tyrosines showed that cleavage was abolished in the Y212A mutant, whereas the Y212F and Y212H mutants retained residual activity. The HD and the alternative 3H motifs were essential for cleavage and the HD domain residues D162 and D267 for metal ion binding. We propose that the active site binds two metal ions, one in a high-affinity and one in a low-affinity site.

Determining the Virulence Properties of Escherichia coli ST131 Containing Bacteriocin-Encoding Plasmids Using Short- and Long-Read Sequencing and Comparing Them with Those of Other E. coli Lineages

Escherichia coli ST131 is a clinical challenge due to its multidrug resistant profile and successful global spread. They are often associated with complicated infections, particularly urinary tract infections (UTIs). Bacteriocins play an important role to outcompete other microorganisms present in the human gut. Here, we characterized bacteriocin-encoding plasmids found in ST131 isolates of patients suffering from a UTI using both short- and long-read sequencing. Colicins Ia, Ib and E1, and microcin V, were identified among plasmids that also contained resistance and virulence genes. To investigate if the potential transmission range of the colicin E1 plasmid is influenced by the presence of a resistance gene, we constructed a strain containing a plasmid which had both the colicin E1 and blaCMY-2 genes. No difference in transmission range was found between transformant and wild-type strains. However, a statistically significantly difference was found in adhesion and invasion ability. Bacteriocin-producing isolates from both ST131 and non-ST131 lineages were able to inhibit the growth of other E. coli isolates, including other ST131. In summary, plasmids harboring bacteriocins give additional advantages for highly virulent and resistant ST131 isolates, improving the ability of these isolates to compete with other microbiota for a niche and thereby increasing the risk of infection.

The secret life of conjugative relaxases

Conjugative relaxases are well-characterized proteins responsible for the site- and strand-specific endonucleolytic cleavage and strand transfer reactions taking place at the start and end of the conjugative DNA transfer process. Most of the relaxases characterized biochemically and structurally belong to the HUH family of endonucleases. However, an increasing number of new families of relaxases are revealing a variety of protein folds and catalytic alternatives to accomplish conjugative DNA processing. Relaxases show high specificity for their cognate target DNA sequences, but several recent reports underscore the importance of their activity on secondary targets, leading to widespread mobilization of plasmids containing an oriT-like sequence. Some relaxases perform other functions associated with their nicking and strand transfer ability, such as catalyzing site-specific recombination or initiation of plasmid replication. They perform these roles in the absence of conjugation, and the validation of these functions in several systems strongly suggest that they are not mere artifactual laboratory observations. Other unexpected roles recently assigned to relaxases include controlling plasmid copy number and promoting retrotransposition. Their capacity to mediate promiscuous mobilization and genetic reorganizations can be exploited for a number of imaginative biotechnological applications. Overall, there is increasing evidence that conjugative relaxases are not only key enzymes for horizontal gene transfer, but may have been adapted to perform other roles which contribute to prokaryotic genetic plasticity. Relaxed target specificity may be key to this versatility.

Evidence of Illegitimate Recombination Between Two Pasteurellaceae Plasmids Resulting in a Novel Multi-Resistance Replicon, pM3362MDR, in Actinobacillus pleuropneumoniae

Evidence of plasmids carrying the tetracycline resistance gene, tet(B), was found in the previously reported whole genome sequences of 14 United Kingdom, and 4 Brazilian, isolates of Actinobacillus pleuropneumoniae. Isolation and sequencing of selected plasmids, combined with comparative sequence analysis, indicated that the four Brazilian isolates all harbor plasmids that are nearly identical to pB1001, a plasmid previously found in Pasteurella multocida isolates from Spain. Of the United Kingdom isolates, 13/14 harbor plasmids that are (almost) identical to pTetHS016 from Haemophilus parasuis. The remaining United Kingdom isolate, MIDG3362, harbors a 12666 bp plasmid that shares extensive regions of similarity with pOV from P. multocida (which carries bla_ROB-1 , sul2, and strAB genes), as well as with pTetHS016. The newly identified multi-resistance plasmid, pM3362MDR, appears to have arisen through illegitimate recombination of pTetHS016 into the stop codon of the truncated strB gene in a pOV-like plasmid. All of the tet(B)-carrying plasmids studied were capable of replicating in Escherichia coli, and predicted origins of replication were identified. A putative origin of transfer (oriT) sequence with similar secondary structure and a nic-site almost identical to that of RP4 was also identified in these plasmids, however, attempts to mobilize them from an RP4-encoding E. coli donor strain were not successful, indicating that specific conjugation machinery may be required.