A Novel marRAB Operon Contributes to the Rifampicin Resistance in Mycobacterium smegmatis

The multiple-antibiotic resistance regulator (MarR) plays an important role in modulating bacterial antibiotic resistance. However, the regulatory model of the marRAB operon in mycobacteria remains to be characterized. Here we report that a MarR, encoded by Ms6508, and its marRAB operon specifically contribute to rifampicin (RIF) resistance in Mycobacterium smegmatis. We show that the MarR recognizes a conserved 21-bp palindromic motif and negatively regulates the expression of two ABC transporters in the operon, encoded by Ms6509–6510. Unlike other known drug efflux pumps, overexpression of these two ABC transporters unexpectedly increased RIF sensitivity and deletion of these two genes increased mycobacterial resistance to the antibiotic. No change can be detected for the sensitivity of recombinant mycobacterial strains to three other anti-TB drugs. Furthermore, HPLC experiments suggested that Ms6509–Ms6510 could pump RIF into the mycobacterial cells. These findings indicated that the mycobacterial MarR functions as a repressor and constitutively inhibits the expression of the marRAB operon, which specifically contributes to RIF resistance in M. smegmatis. Therefore, our data suggest a new regulatory mechanism of RIF resistance and also provide the new insight into the regulatory model of a marRAB operon in mycobacteria.     

An Operon of Three Transcriptional Regulators Controls Horizontal Gene Transfer of the Integrative and Conjugative Element ICEclc in Pseudomonas knackmussii B13

The integrative and conjugative element ICEclc is a mobile genetic element in Pseudomonas knackmussii B13, and an experimental model for a widely distributed group of elements in Proteobacteria. ICEclc is transferred from specialized transfer competent cells, which arise at a frequency of 3-5% in a population at stationary phase. Very little is known about the different factors that control the transfer frequency of this ICE family. Here we report the discovery of a three-gene operon encoded by ICEclc, which exerts global control on transfer initiation. The operon consists of three consecutive regulatory genes, encoding a TetR-type repressor MfsR, a MarR-type regulator and a LysR-type activator TciR. We show that MfsR autoregulates expression of the operon, whereas TciR is a global activator of ICEclc gene expression, but no clear role was yet found for MarR. Deletion of mfsR increases expression of tciR and marR, causing the proportion of transfer competent cells to reach almost 100% and transfer frequencies to approach 1 per donor. mfsR deletion also caused a two orders of magnitude loss in population viability, individual cell growth arrest and loss of ICEclc. This indicates that autoregulation is an important feature maintaining ICE transfer but avoiding fitness loss. Bioinformatic analysis showed that the mfsR-marR-tciR operon is unique for ICEclc and a few highly related ICE, whereas tciR orthologues occur more widely in a large variety of suspected ICE among Proteobacteria.     

Structural mechanism of transcription regulation of the Staphylococcus aureus multidrug efflux operon mepRA by the MarR family repressor MepR

The multidrug efflux pump MepA is a major contributor to multidrug resistance in Staphylococcus aureus. MepR, a member of the multiple antibiotic resistance regulator (MarR) family, represses mepA and its own gene. Here, we report the structure of a MepR–mepR operator complex. Structural comparison of DNA-bound MepR with ‘induced’ apoMepR reveals the large conformational changes needed to allow the DNA-binding winged helix-turn-helix motifs to interact with the consecutive major and minor grooves of the GTTAG signature sequence. Intriguingly, MepR makes no hydrogen bonds to major groove nucleobases. Rather, recognition-helix residues Thr60, Gly61, Pro62 and Thr63 make sequence-specifying van der Waals contacts with the TTAG bases. Removing these contacts dramatically affects MepR–DNA binding activity. The wings insert into the flanking minor grooves, whereby residue Arg87, buttressed by Asp85, interacts with the O2 of T₄ and O4′ ribosyl oxygens of A₂₃ and T₄. Mutating Asp85 and Arg87, both conserved throughout the MarR family, markedly affects MepR repressor activity. The His14′:Arg59 and Arg10′:His35:Phe108 interaction networks stabilize the DNA-binding conformation of MepR thereby contributing significantly to its high affinity binding. A structure-guided model of the MepR–mepA operator complex suggests that MepR dimers do not interact directly and cooperative binding is likely achieved by DNA-mediated allosteric effects.     

Whole-genome sequencing and gene sharing network analysis powered by machine learning identifies antibiotic resistance sharing between animals,humans and environment in livestock farming

Anthropogenic environments such as those created by intensive farming of livestock, have been proposed to provide ideal selection pressure for the emergence of antimicrobial-resistant Escherichia coli bacteria and antimicrobial resistance genes (ARGs) and spread to humans. Here, we performed a longitudinal study in a large-scale commercial poultry farm in China, collecting E. coli isolates from both farm and slaughterhouse; targeting animals, carcasses, workers and their households and environment. By using whole-genome phylogenetic analysis and network analysis based on single nucleotide polymorphisms (SNPs), we found highly interrelated non-pathogenic and pathogenic E. coli strains with phylogenetic intermixing, and a high prevalence of shared multidrug resistance profiles amongst livestock, human and environment. Through an original data processing pipeline which combines omics, machine learning, gene sharing network and mobile genetic elements analysis, we investigated the resistance to 26 different antimicrobials and identified 361 genes associated to antimicrobial resistance (AMR) phenotypes; 58 of these were known AMR-associated genes and 35 were associated to multidrug resistance. We uncovered an extensive network of genes, correlated to AMR phenotypes, shared among livestock, humans, farm and slaughterhouse environments. We also found several human, livestock and environmental isolates sharing closely related mobile genetic elements carrying ARGs across host species and environments. In a scenario where no consensus exists on how antibiotic use in the livestock may affect antibiotic resistance in the human population, our findings provide novel insights into the broader epidemiology of antimicrobial resistance in livestock farming. Moreover, our original data analysis method has the potential to uncover AMR transmission pathways when applied to the study of other pathogens active in other anthropogenic environments characterised by complex interconnections between host species.     

Mechanism for Attenuation of DNA Binding by MarR Family Transcriptional Regulators by Small Molecule Ligands

Members of the multiple antibiotic resistance regulator (MarR) family control gene expression in a variety of metabolic processes in bacteria and archaea. Hypothetical uricase regulator (HucR), which belongs to the ligand-responsive branch of the MarR family, regulates uricase expression in Deinococcus radiodurans by binding a shared promoter region between uricase and HucR genes. We show here that HucR responds only to urate and, to a lesser extent, to xanthine by attenuated DNA binding, compared to other intermediates of purine degradation. Using molecular-dynamics-guided mutational analysis, we identified the ligand-binding site in HucR. Electrophoretic mobility shift assays and intrinsic Trp fluorescence have identified W20 from the N-terminal helix and R80 from helix 3, which serves as a scaffold for the DNA recognition helix, as being essential for ligand binding. Using structural data combined with in silico and in vitro analyses, we propose a mechanism for the attenuation of DNA binding in which a conformational change initiated by charge repulsion due to a bound ligand propagates to DNA recognition helices. This mechanism may apply generally to MarR homologs that bind anionic phenolic ligands.     

Chloroplast genome characteristics and phylogenetic analysis of the medicinal plant Blumea balsamifera (L.) DC

Blumea balsamifera (L.) DC., a medicinal plant with high economic value in the Asteraceae family, is widely distributed in China and Southeast Asia. However, studies on the population structure or phylogenetic relationships with other related species are rare owing to the lack of genome information. In this study, through high-throughput sequencing, we found that the chloroplast genome of B. balsamifera was 151,170 bp in length, with a pair of inverted repeat regions (IRa and IRb) comprising 24,982 bp, a large single-copy (LSC) region comprising 82,740 bp, and a small single-copy (SSC) region comprising 18,466 bp. A total of 130 genes were identified in the chloroplast genome of B. balsamifera, including 85 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes; furthermore, sequence analysis identified 53 simple sequence repeats. Whole chloroplast genome comparison indicated that the inverted regions (IR) were more conserved than large single-copy and SSC regions. Phylogenetic analysis showed that B. balsamifera is closely related to Pluchea indica. Conclusively, the chloroplast genome of B. balsamifera was helpful for species identification and analysis of the genetic diversity and evolution in the genus Blumea and family Asteraceae.     

Complete chloroplast genome of the medicinal plant Evolvulus alsinoides: comparative analysis,identification of mutational hotspots and evolutionary dynamics with species of Solanales

Evolvulus alsinoides, belonging to the family Convolvulaceae, is an important medicinal plant widely used as a nootropic in the Indian traditional medicine system. In the genus Evolvulus, no research on the chloroplast genome has been published. Hence, the present study focuses on annotation, characterization, identification of mutational hotspots, and phylogenetic analysis in the complete chloroplast genome (cp) of E. alsinoides. Genome comparison and evolutionary dynamics were performed with the species of Solanales. The cp genome has 114 genes (80 protein-coding genes, 30 transfer RNA, and 4 ribosomal RNA genes) that were unique with total genome size of 157,015 bp. The cp genome possesses 69 RNA editing sites and 44 simple sequence repeats (SSRs). Predicted SSRs were randomly selected and validated experimentally. Six divergent hotspots such as trnQ-UUG, trnF-GAA, psaI, clpP, ndhF, and ycf1 were discovered from the cp genome. These microsatellites and divergent hot spot sequences of the Taxa ‘Evolvulus’ could be employed as molecular markers for species identification and genetic divergence investigations. The LSC area was found to be more conserved than the SSC and IR region in genome comparison. The IR contraction and expansion studies show that nine genes rpl2, rpl23, ycf1, ycf2, ycf1, ndhF, ndhA, matK, and psbK were present in the IR-LSC and IR-SSC boundaries of the cp genome. Fifty-four protein-coding genes in the cp genome were under negative selection pressure, indicating that they were well conserved and were undergoing purifying selection. The phylogenetic analysis reveals that E. alsinoides is closely related to the genus Cressa with some divergence from the genus Ipomoea. This is the first time the chloroplast genome of the genus Evolvulus has been published. The findings of the present study and chloroplast genome data could be a valuable resource for future studies in population genetics, genetic diversity, and evolutionary relationship of the family Convolvulaceae.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12298-021-01051-w.     

Overview of chitin metabolism enzymes in Manduca sexta: Identification,domain organization,phylogenetic analysis and gene expression

Chitin is one of the most abundant biomaterials in nature. The biosynthesis and degradation of chitin in insects are complex and dynamically regulated to cope with insect growth and development. Chitin metabolism in insects is known to involve numerous enzymes, including chitin synthases (synthesis of chitin), chitin deacetylases (modification of chitin by deacetylation) and chitinases (degradation of chitin by hydrolysis). In this study, we conducted a genome-wide search and analysis of genes encoding these chitin metabolism enzymes in Manduca sexta. Our analysis confirmed that only two chitin synthases are present in M. sexta as in most other arthropods. Eleven chitin deacetylases (encoded by nine genes) were identified, with at least one representative in each of the five phylogenetic groups that have been described for chitin deacetylases to date. Eleven genes encoding for family 18 chitinases (GH18) were found in the M. sexta genome. Based on the presence of conserved sequence motifs in the catalytic sequences and phylogenetic relationships, two of the M. sexta chitinases did not cluster with any of the current eight phylogenetic groups of chitinases: two new groups were created (groups IX and X) and their characteristics are described. The result of the analysis of the Lepidoptera-specific chitinase-h (group h) is consistent with its proposed bacterial origin. By analyzing chitinases from fourteen species that belong to seven different phylogenetic groups, we reveal that the chitinase genes appear to have evolved sequentially in the arthropod lineage to achieve the current high level of diversity observed in M. sexta. Based on the sequence conservation of the catalytic domains and on their developmental stage- and tissue-specific expression, we propose putative functions for each group in each category of enzymes.