Review

An overview of the antimicrobial resistance mechanisms of bacteria

  • Received: 18 April 2018 Accepted: 13 June 2018 Published: 26 June 2018
  • Resistance to antimicrobial agents has become a major source of morbidity and mortality worldwide. When antibiotics were first introduced in the 1900’s, it was thought that we had won the war against microorganisms. It was soon discovered however, that the microorganisms were capable of developing resistance to any of the drugs that were used. Apparently most pathogenic microorganisms have the capability of developing resistance to at least some antimicrobial agents. The main mechanisms of resistance are: limiting uptake of a drug, modification of a drug target, inactivation of a drug, and active efflux of a drug. These mechanisms may be native to the microorganisms, or acquired from other microorganisms. Understanding more about these mechanisms should hopefully lead to better treatment options for infective diseases, and development of antimicrobial drugs that can withstand the microorganisms attempts to become resistant.

    Citation: Wanda C Reygaert. An overview of the antimicrobial resistance mechanisms of bacteria[J]. AIMS Microbiology, 2018, 4(3): 482-501. doi: 10.3934/microbiol.2018.3.482

    Related Papers:

    [1] Fadoua El Moustaid, Amina Eladdadi, Lafras Uys . Modeling bacterial attachment to surfaces as an early stage of biofilm development. Mathematical Biosciences and Engineering, 2013, 10(3): 821-842. doi: 10.3934/mbe.2013.10.821
    [2] Yousef Rohanizadegan, Stefanie Sonner, Hermann J. Eberl . Discrete attachment to a cellulolytic biofilm modeled by an Itô stochastic differential equation. Mathematical Biosciences and Engineering, 2020, 17(3): 2236-2271. doi: 10.3934/mbe.2020119
    [3] Paolo Fergola, Marianna Cerasuolo, Edoardo Beretta . An allelopathic competition model with quorum sensing and delayed toxicant production. Mathematical Biosciences and Engineering, 2006, 3(1): 37-50. doi: 10.3934/mbe.2006.3.37
    [4] Nikodem J. Poplawski, Abbas Shirinifard, Maciej Swat, James A. Glazier . Simulation of single-species bacterial-biofilm growth using the Glazier-Graner-Hogeweg model and the CompuCell3D modeling environment. Mathematical Biosciences and Engineering, 2008, 5(2): 355-388. doi: 10.3934/mbe.2008.5.355
    [5] Ting Kang, Yanyan Du, Ming Ye, Qimin Zhang . Approximation of invariant measure for a stochastic population model with Markov chain and diffusion in a polluted environment. Mathematical Biosciences and Engineering, 2020, 17(6): 6702-6719. doi: 10.3934/mbe.2020349
    [6] Ana I. Muñoz, José Ignacio Tello . Mathematical analysis and numerical simulation of a model of morphogenesis. Mathematical Biosciences and Engineering, 2011, 8(4): 1035-1059. doi: 10.3934/mbe.2011.8.1035
    [7] Xiaomei Bao, Canrong Tian . Turing patterns in a networked vegetation model. Mathematical Biosciences and Engineering, 2024, 21(11): 7601-7620. doi: 10.3934/mbe.2024334
    [8] Chang-Yuan Cheng, Shyan-Shiou Chen, Xingfu Zou . On an age structured population model with density-dependent dispersals between two patches. Mathematical Biosciences and Engineering, 2019, 16(5): 4976-4998. doi: 10.3934/mbe.2019251
    [9] Jacques A. L. Silva, Flávia T. Giordani . Density-dependent dispersal in multiple species metapopulations. Mathematical Biosciences and Engineering, 2008, 5(4): 843-857. doi: 10.3934/mbe.2008.5.843
    [10] Yue Xia, Lijuan Chen, Vaibhava Srivastava, Rana D. Parshad . Stability and bifurcation analysis of a two-patch model with the Allee effect and dispersal. Mathematical Biosciences and Engineering, 2023, 20(11): 19781-19807. doi: 10.3934/mbe.2023876
  • Resistance to antimicrobial agents has become a major source of morbidity and mortality worldwide. When antibiotics were first introduced in the 1900’s, it was thought that we had won the war against microorganisms. It was soon discovered however, that the microorganisms were capable of developing resistance to any of the drugs that were used. Apparently most pathogenic microorganisms have the capability of developing resistance to at least some antimicrobial agents. The main mechanisms of resistance are: limiting uptake of a drug, modification of a drug target, inactivation of a drug, and active efflux of a drug. These mechanisms may be native to the microorganisms, or acquired from other microorganisms. Understanding more about these mechanisms should hopefully lead to better treatment options for infective diseases, and development of antimicrobial drugs that can withstand the microorganisms attempts to become resistant.


    DNase I hypersensitive sites (DHSs) are a specific genomic region in the chromatin, which is of hypersensitivity to cleavage by the DNase I enzyme [1]. DHSs untie its condensed structure, which makes the DNA exposed and accessible to the regulatory proteins. DHSs are functionally associated with the cis-regulatory elements such as promoters, enhancers, suppressors, insulators, as well as locus control regions [2]. Thus, mapping DHSs is becoming one of the most effective methods to precisely identify the location of many different regulatory elements in specific, well-studied genes [3]. Genetic variations in DHSs were found to be implicated in a wide spectrum of common diseases and traits, including Alzheimer's disease [4,5,6,7,8]. For example, DHSs were identified as driver distal regulatory elements in breast cancer, and were responsible for the aberrant expression of neighboring genes [9].

    Identifying DHSs is of great interest to cis-regulatory element annotation. With advances in next-generation sequencing, many high-throughput techniques have been developed to detect DHSs in the past decades [5,10,11,12], such as the Southern blot approach [13] and DNase-seq [14]. Zhang et al. [15] developed a DNase-seq procedure for genome-wide mapping of DHSs in Arabidopsis thaliana and rice, while Wang et al. proposed a modified DNase-seq for genome-wide identification of DHSs in plants [16]. These experimentally verified DHSs were collected to be deposited in several public databases for further exploration [1].

    Although these high-throughput techniques have contributed much to the discovery of thousands of DHSs, they have two inherent limitations: they are expensive and laborious, which make them insufficient to complete the challenging task of detecting DHSs from tremendous volumes of genomes. Computational identification is another routine to detect DHSs. The computational identification is defined as computational models or functions which are able to predict DHSs after they are trained by known DHSs. Computational identification is extremely cheaper and faster than high-throughput techniques for DHSs detection, and thus it is becoming an alternative to identify DHSs. Computational identification based on machine learning, especially deep learning, has extensively been applied to predict transcription factor binding sites [17,18,19,20,21,22,23] and to mine DNA/RNA motif [24]. For example, Wang et al. [17] created a hybrid convolutional recurrent neural network (RNN) for predicting transcription factor binding sites which obtained the state-of-the-art performance on 66 in vitro datasets. Zhang et al. [18] presented a deep learning-based method for transcription factor-DNA binding signal prediction that was able to deal with up to four transcription factor binding-related tasks. Wang et al. [19] employed fully convolutional neural networks (CNNs), along with gated recurrent units, to localize transcription factor binding sites. Following these successful practices, no less than 10 computational models or methods have been created for DHSs detection over the recent decade. These models or methods can be grouped into traditional machine learning-based methods [25,26,27,28], ensemble learning-based methods [29,30,31,32,33,34], and deep learning-based methods [35,36]. To the best of our knowledge, the first computational predictor for DHSs was the support vector machine-based method, which was proposed by Noble et al. [25] in 2005. This method used nucleotide composition as a representation of DNA sequences. Evidently, the nucleotide composition is unable to sufficiently represent DNA sequences because it drops out information about position and order. Feng et al. [37] used pseudo nucleotide composition [38,39,40] to integrate local and global sequence-order effects of DNA sequences. The pseudo nucleotide composition is similar to the pseudo amino acid composition which is a popular and effective representation for protein sequences. Liu et al. [30] computed nucleotide composition, reverse nucleotide composition, and pseudo dinucleotide composition to build three respective random forest-based classifiers. Three single random forest-based classifiers were fused as an ensemble classifier for DHSs prediction. Zhang et al. [41] employed reverse complement k-mer and dinucleotide-based auto covariance to represent DNA sequences. Zhang et al. [29] stacked multiple traditional machine learning algorithms to build an ensemble classifier for DHSs prediction. Zhang et al. [29] also employed the LASSO to reduce the dimension of representations and the SMOTE-Tomek to overcome imbalance between positive and negative samples. Zheng et al. [32] extracted composition information and physicochemical properties and used a boosting algorithm to optimize informative representations. Zhang et al. [33] proposed a dinucleotide-based spatial autocorrelation to represent DNA sequences. The aforementioned methods heavily rely on representations since it determines, to a great extent, whether methods are performed well or not. High effective representations are generally difficult to obtain in practice. Deep learning is merging as a representative of next-generation artificial intelligence, exhibiting vast potential to solve challenges unsolved in the past. Deep learning has been applied in a wide range of fields, including academic and industrial communities. Lyu et al. [36] developed a deep learning method for DHSs prediction which employed CNNs, along with the gate mechanism, to extract representation. To deal with variable-length inputs, Lyu et al. [36] used the spatial pyramid pooling [42], which was initially proposed to deal with variable-size images. The CNNs are able to capture local properties, and thus has extensively been used in the fields of image analysis and natural language processing. However, the CNNs are insufficient to capture dependencies between local parts, such as words. In the text sequences, dependencies between words are vital because they determine whether one understands it or not. Dao et al. [43] combined the CNNs and long-short term memory (LSTM), which is a special RNN for DHSs prediction. Dao et al. [43] stacked the CNNs and the LSTM in order, i.e., they first used the CNNs to capture local characteristics, and then used the LSTM to catch dependencies between local characteristics. The CNNs and the LSTM characterized different properties of sequences. Linearly stacking the CNN and the LSTM would lose their respective merits. In this paper, we stacked the CNNs and the LSTM in a parallel manner, which absorbed two respective representations learned by the CNN and the LSTM. In addition, we used feed-forward attention to improve representations by the LSTM.

    We downloaded DHS datasets from 14 different tissues and 4 developmental stages in mouse which are available at the following website: http://lin-group.cn/server/iDHS-Deep/. These DHSs were collected according to the atlas of the DHSs created by Breeze et al. [44]. Dao et al. [43] further preprocessed these datasets for training of the iDHS-Deep model, including choosing the DHSs of 50 to 300 bp as positive samples, selecting specific DNA fragments as negative samples, removing or reducing the homology between sequences by using CD-Hit [45,46], which is a sequence-clustering tool, and dividing these datasets into the training set and the independent set at the ratio of 7 to 3. The numbers of positive and negative samples were not equal for the Stomach tissues, which were respectively 1062 and 2125 in the training set, and which were respectively 456 and 911 in the independent set. Except for Stomach tissues, the numbers of positive and negative samples were identical for each tissue or each developmental stage. The number of positive samples were respectively 7114, 10,299, 5766, 6519, 7424, 30,929, 6316, 4978, 1612, 2515, 3511, 2877, 1224, 7512, 52,418, 16,172 and 21,247 for 13 tissues (i.e., Forebrain, Midbrain, Hindbrain, Liver, Lung, Heart, Kidney, Limb, Thymus, Craniofacial, Retina, Muller retina, and Neural tube) and 4 developmental stages (i.e., ESC, Early-Fetal, Late-Fetal and Adult) in the training set, while they were respectively 3049, 4414, 2472, 2795, 3183, 13,256, 2708, 2134,692, 1078, 1506, 1234,525, 3224, 22,466, 6933, 9106 in the independent set.

    As shown in Figure 1, the architecture of the proposed LangMoDHS mainly comprised the embedding, the CNNs, the Bi-LSTM followed by the feed-forward attention, the dropout, the fully-connected layer and the output layer. Unlike the iDHS-Deep [43], the LangMoDHS stacked the CNNs and the Bi-LSTM layer in a parallel manner. DNA sequences were first translated into integer sequences which were actually immediate input to the LangMoDHS. Each character in the DNA sequence was mapped into an integer as follows: A into 1, C into 2, G into 6 and T into 17 [43]. The same character-encoding scheme was adopted by the iDHS-Deep [43]. The integer sequences were embedded into continuous vectors, which were further characterized by the Bi-LSTM, CNNs and feed-forward attention. The output layer consisted of only a neuron that represented the probability of the input sequence containing the DHSs. The LangMoDHS is similar to the Deep6mAPred [47] in terms of the architecture, which is a deep learning method for predicting DNA N6-methyladenosine sites except that the former replaced one-hot encoding with an embedding layer and used two layers of CNN.

    Figure 1.  Flowchart of the LangMoDHS. The numbers in parentheses represent the parameters of each layer of the network, $ {h}_{1} $, $ {h}_{2} $, $ \cdots $ $ {h}_{T} $ represent the vector in the sequence of hidden states, $ a\left({h}_{T}\right) $ is a learnable function, $ {a}_{1} $, $ {a}_{2} $, $ \cdots $ $ {a}_{T} $ is the probability vector, the vector c is computed as a weighted average of $ {h}_{T} $, $ \oplus $ represents element-wise sum and $ \otimes $ represents the element-wise product. Dense denotes the fully-connected layer.

    The embedding layer was intended to bridge the sequences of discrete variables and sequences of continuous vectors. In the deep neural network, the embedding layer is actually a specific neural network analogue to the Word2vec [48,49]. The embedding layer overcame the conventional issues, such as sparsity and no correlation between words. We used embedding in the Keras as the first layer in the LangMoDHS. The Keras is an open-source and extensively used deep learning toolkit.

    The CNN is increasingly becoming one of the best popular neural networks, which was initially pioneered by Fukushima et al. [50] as an extension of the concept of receptive fields [51]. Since LeCun et al. [52] introduced the gradient back propagation algorithm to train the CNNs, the CNNs have attracted more and more attention, especially from deep learning communities. The CNNs contain two basic operations: convolution and pooling. The convolution is to multiply the input by a fixed convolutional kernel in the same layer. The convolutional kernel is like a filter in the digital signal field, which is shared by the same input. In addition, the CNNs use the pooling to reduce overfitting or computation. Actually, the pooling is a down-sampling method, including average pooling and max pooling. The CNNs are divided into 1D, 2D and 3D CNNs. The 2D CNN and the 3D CNN are generally applied for image data analysis, while the 1D CNN is applied to the field of text analysis. In this study, we used two 1D CNNs, each followed by the max pooling layer.

    LSTM [53] is a special RNN that is different from the CNN. One of the main characteristics of LSTM is to share weights at all the time steps. LSTM is capable of preserving previous semantics by the cell state, which is controlled by the gate mechanism. For example, it uses an input gate to determine how much information is updated, and it uses a forget gate to decide what information is removed in the cell state. Therefore, LSTM is able to capture dependency between words in a sequence and thus is especially suitable to deal with sequence analysis. We used two LSTMs, i.e., Bi-LSTM, to capture the semantic relationship between words.

    The attention mechanism is a newly developed technique of deep learning, and it has extensively been applied in the field of computer vision, natural language processing and bioinformatics. All of the deep learning-based language models, such as transformer and Bert employed attention mechanisms. Even Vaswani et al. declared that attention was all you need [54]. The attention mechanism is actually a scheme to allocate different weights to different parts. In the recent five years, many attention mechanisms have been proposed, including feed-forward attention [55] and self-attention [54]. Here, we used feed-forward attention to compensate for the deficiency of Bi-LSTM. The feed-forward attention is computed by

    $ {a}_{t} = \frac{exp\left({e}_{t}\right)}{{\sum }_{k = 1}^{T}exp\left({e}_{k}\right)} , $ (1)

    where

    $ {e}_{t} = \delta \left({h}_{t}\right) . $ (2)

    $ \delta $ is the learnable parameter and $ {h}_{t} $ is the hidden state at the time step of t in a Bi-LSTM. The output is a sum of the attentions multiplying corresponding hidden states, which is computed by

    $ c = {\sum }_{t = 1}^{T}{\theta }_{t}{h}_{t} . $ (3)

    We employed the Receiver Operator Characteristic (ROC) curve, Precision-Recall (PR) curve and F1-score to measure performance. The ROC curve plots the true positive rate (TPR) on the y-axis against the false positive rate (FPR) on the x-axis under the various thresholds. The PR curve plots the Precision on the y-axis against Recall on the x-axis. The F1-score is the summed mean of the Precision and Recall. These metrics are respectively defined by

    $ FPR = \frac{{F}_{P}}{{F}_{P}+{T}_{N}} , $ (4)
    $ TPR = Recall = \frac{{T}_{P}}{{T}_{P}+{F}_{N}} , $ (5)
    $ Precision = \frac{{T}_{p}}{{T}_{P}+{F}_{P}} , $ (6)
    $ F1-score = \frac{2\times Precision\times Recall}{Precision+Recall} . $ (7)

    The more the ROC curve trends to the upper left, the better performance. The more upper right the PR curve trends to upper right, the better performance. The area under the ROC curve (AUROC) and the area under the PR curve (AUPRC) were used to quantitatively assess performance. The AUROC and AUPRC range from 0 to 1, where 1 means a perfect prediction, 0.5 is a random prediction, and 0 is a completely reverse prediction.

    We performed a 5-fold cross-validation and an independent test to check the performance of the proposed method. In the 5-fold cross-validation, the training set was randomly divided into 5 parts of equal or approximate size. 4 parts were used to train the model, and the remaining part was used to test the trained model. The process was repeated five times. The independent test was to use the independent set to test the model trained by the training set. Figure 2A, B showed the AUROC values for the 5-fold cross-validations for the training sets from the 14 tissues and the 4 developmental stages, respectively. Obviously, all of the standard deviations of AUROC values over 5-fold cross-validation were less than 0.058, indicating that the LangMoDHS performed stably. Figure 3A, B showed the ROC curves for the independent set. The LangMoDHS achieved the best performance for the Heart tissue (AUROC = 0.960) and performed the worst for the Thymus tissue (AUROC = 0.770). The range of AUROC values for the 14 tissues was 0.19, implying that the LangMoDHS performed differently for different tissues. The LangMoDHS performed stably across the 4 stages. The highest AUROC value was 0.952, the lowest value was 0.910, and the range was 0.042, which was far smaller than that for the 14 tissues. Figures 4 and 5 respectively showed the PR curves and F1-scores for the independent set. The similar phenomenon was observed. For example, the LangMoDHS reached the best AUPRC value, as well as the best F1-score for the Heart tissue, and it achieved the best AUPRC value, as well as the best F1-score for the Early-Fetal stage.

    Figure 2.  The box plot of AUROC values of 5-fold cross-validations for the training sets from the (A) 14 tissues and (B) from 4 developmental stages.
    Figure 3.  ROC curves of independent tests for the (A) 14 different tissues and (B) 4 different developmental stages.
    Figure 4.  PR curves of independent tests for (A) 14 different tissues and (B) 4 different developmental stages.
    Figure 5.  Bar charts of F1-score results for the independent tests for (A) 14 different tissues and (B) 4 different developmental stages.

    We compared the LangMoDHS with the iDHS-Deep [43], which is a newly developed method for predicting DHSs. Table 1 listed the AUROC values of 5-fold cross-validation and independent test for the 14 tissues. The iDHS-Deep outperformed the LangMoDHS in both 5-fold cross-validation and independent test over 3 tissues: Stomach, Thymus, and Neural tube, while LangMoDHS completely outperformed iDHS-Deep over 4 tissues: Hindbrain, Liver, Lung and Heart. Over 3 tissues: Limb, Craniofacial, and Retina, iDHS-Deep performed better in 5-fold cross-validation, while LangMoDHS performed better in independent test. For the Kidney and Midbrain tissues, LangMoDHS was equal to iDHS-Deep in terms of the 5-fold cross-validation, while LangMoDHS performed better than iDHS-Deep in the independent test. Over the Forebrain and the Muller retina tissues, both methods is equivalent in the independent test. Table 2 listed all of the AUROC values in the 5-fold cross-validations and independent tests over 4 developmental stages. Although LangMoDHS performed worse than iDHS-Deep over two developmental stages: ESC and Late-Fetal in the 5-fold cross-validations, the former completely outperformed the latter over all 4 developmental stages in the independent tests.

    Table 1.  AUROC values of 5-fold cross-validations and independent tests for 14 tissues.
    DATASETS (TISSUES) METHOD
    IDHS-Deep LangMoDHS IDHS-Deep LangMoDHS
    Training datasets Independent datasets
    Forebrain 0.934 0.938 0.939 0.939
    Midbrain 0.931 0.931 0.920 0.932
    Hindbrain 0.911 0.915 0.914 0.926
    Liver 0.927 0.932 0.924 0.935
    Lung 0.906 0.920 0.885 0.919
    Heart 0.955 0.957 0.949 0.960
    Kidney 0.934 0.934 0.923 0.938
    Limb 0.909 0.907 0.908 0.918
    Stomach 0.877 0.848 0.931 0.836
    Thymus 0.921 0.738 0.896 0.770
    Craniofacial 0.908 0.871 0.894 0.901
    Retina 0.902 0.900 0.894 0.911
    Muller retina 0.904 0.882 0.901 0.901
    Neural tube 0.896 0.763 0.900 0.804

     | Show Table
    DownLoad: CSV
    Table 2.  AUROC values in the 5-fold cross-validation and independent tests for the 4 developmental stages.
    DATASETS (STAGES) METHOD
    IDHS-Deep LangMoDHS IDHS-Deep LangMoDHS
    Training datasets Independent datasets
    ESC 0.923 0.920 0.899 0.921
    Early-Fetal 0.949 0.950 0.940 0.952
    Late-Fetal 0.923 0.907 0.901 0.910
    Adult 0.916 0.925 0.905 0.930

     | Show Table
    DownLoad: CSV

    We further tested LangMoDHS for the ability to predict DHSs across tissues (developmental stages). That is to say, the LangMoDHS trained by the dataset from A tissue (developmental stage), was used to predict DHSs from B tissue (developmental stage). Tables 3 and 4 listed the AUROC value of independent tests across tissues and developmental stages, respectively. Except for seven tissues: Heart, Kidney, Stomach, Thymus, Craniofacial, Muller retina, and Neural tube, LangMoDHS exhibited better performance over other tissues different from itself. For example, the LangMoDHS trained by the Forebrain training set achieved an AUROC value of 0.939 over the independent set from Forebrain, but obtained a better AUROC value (0.950) over the independent set from Heart tissue. This indicated that there existed a potential for LangMoDHS to predict DHSs across tissues. However, not all of the cross-tissue performance of LangMoDHS was better. For example, the LangMoDHS trained by Craniofacial training set reached an AUROC value of 0.901 over the Craniofacial independent set, which was better than those over all of the independent sets from other tissues. The similar phenomenon was observed in Table 4.

    Table 3.  AUROC values of independent tests across tissues in mouse genomes.
    Training datasets Forebr-ain Midbr-ain Hindb-rain Liver Lung Heart Kidn-ey Limb Stom-ach Thym-us Cranio-facial Retina Muller retina Neural tube
    Independent datasets
    Forebrain 0.939 0.939 0.922 0.923 0.926 0.944 0.904 0.930 0.789 0.721 0.796 0.914 0.708 0.791
    Midbrain 0.923 0.932 0.918 0.917 0.918 0.929 0.901 0.922 0.759 0.725 0.786 0.908 0.730 0.764
    Hindbrain 0.913 0.918 0.926 0.910 0.910 0.919 0.895 0.914 0.723 0.698 0.765 0.908 0.734 0.732
    Liver 0.901 0.907 0.907 0.935 0.924 0.907 0.922 0.910 0.671 0.658 0.728 0.901 0.747 0.672
    Lung 0.885 0.884 0.885 0.896 0.919 0.884 0.892 0.885 0.662 0.655 0.722 0.882 0.732 0.659
    Heart 0.950 0.953 0.938 0.938 0.939 0.960 0.920 0.947 0.797 0.710 0.789 0.934 0.700 0.792
    Kidney 0.889 0.896 0.901 0.919 0.913 0.884 0.938 0.894 0.594 0.608 0.673 0.903 0.755 0.624
    Limb 0.911 0.916 0.910 0.910 0.906 0.920 0.890 0.918 0.744 0.704 0.780 0.898 0.739 0.742
    Stomach 0.935 0.945 0.912 0.920 0.920 0.950 0.893 0.917 0.836 0.735 0.812 0.895 0.689 0.801
    Thymus 0.899 0.909 0.895 0.902 0.907 0.900 0.883 0.899 0.719 0.770 0.774 0.882 0.743 0.724
    Craniofacial 0.892 0.903 0.904 0.897 0.904 0.895 0.890 0.903 0.710 0.712 0.901 0.888 0.754 0.715
    Retina 0.888 0.895 0.897 0.899 0.892 0.887 0.894 0.893 0.673 0.674 0.742 0.911 0.747 0.680
    Muller retina 0.843 0.855 0.861 0.875 0.866 0.824 0.886 0.855 0.571 0.627 0.680 0.860 0.901 0.613
    Neural tube 0.931 0.933 0.913 0.915 0.918 0.933 0.895 0.926 0.750 0.711 0.788 0.910 0.747 0.804

     | Show Table
    DownLoad: CSV
    Table 4.  AUROC values of independent tests across developmental stages in mouse genomes.
    Training datasets ESC Early-Fetal Late-Fetal Adult
    Independent datasets
    ESC 0.921 0.919 0.908 0.911
    Early-Fetal 0.940 0.952 0.937 0.934
    Late-Fetal 0.890 0.891 0.910 0.903
    Adult 0.905 0.902 0.908 0.930

     | Show Table
    DownLoad: CSV

    LangMoDHS consisted mainly of three components: CNN, Bi-LSTM, and feed-forward attention. We investigated further how much the CNN, Bi-LSTM, and feed-forward attention contributed to the recognition of DHSs. Figure 6AC showed the ROC curves of the independent tests by respectively removing the CNN, Bi-LSTM and attention from LangMoDHS. In contrast, Figure 6D showed the ROC curves of LangMoDHS over the independent tests. It was easy to observe that removing any one of 3 components caused the AUROC values to descend, implying that each contributed substantially to the recognition of DHSs. However, the contributions varied with the tissue (developmental stage) and component. That is to say, some components contributed more for some tissues or developmental stages than for other tissues or stages. For example, the AUROC value after removal of Bi-LSTM was higher than those after respective removal of the CNN and attention for Neural tube tissue, indicating that the CNN and attention contributed more than Bi-LSTM for Neural tube tissue. The AUROC value after removal of attention was more than those after the respective removal of the CNN and Bi-LSTM for Liver tissue, indicating that the CNN and Bi-LSTM contributed more than attention for the Liver tissue. The CNN contributed more than both the Bi-LSTM and attention for 4 developmental stages, the Bi-LSTM contributed more than the attention for 3 developmental stages: ESC, Late-Fetal, and Adult, while attention contributed more than Bi-LSTM for the Early-Fetal developmental stage.

    Figure 6.  ROC curves of the independent tests by removing the CNN, Bi-LSTM and feed-forward attention. (A) ROC curves of only Bi-LSTM. (B) ROC curves of Bi-LSTM + feed-forward attention. (C) ROC curves of CNN + BiLSTM. (D) ROC curves of the LangMoDHS.

    We employed information entropy to analyze the motif of DHS sequences. The position-specific nucleotide matrix is defined by

    $ Z = \left(\begin{array}{cc}\begin{array}{cc}{z}_{11}& \begin{array}{cc}{z}_{12}& \cdots \end{array}
    \begin{array}{cc}\begin{array}{cc}{z}_{11}& \begin{array}{cc}{z}_{12}& \cdots \end{array}
    \end{array}& {z}_{1n}\\ \begin{array}{cc}\begin{array}{c}\begin{array}{c}{z}_{21}\\ {z}_{31}\end{array}
    \begin{array}{cc}\begin{array}{c}\begin{array}{c}{z}_{21}\\ {z}_{31}\end{array}
    \\ {z}_{41}\end{array}& \begin{array}{c}\begin{array}{cc}\begin{array}{c}{z}_{22}\\ {z}_{32}\end{array}
    \begin{array}{c}\begin{array}{cc}\begin{array}{c}{z}_{22}\\ {z}_{32}\end{array}
    &
    \end{array}\\ z42
    \end{array}\end{array}& \begin{array}{c}\begin{array}{c}{z}_{2n}\\ {z}_{3n}\end{array}
    \\ {z}_{4n}\end{array}\end{array}\right) , $
    (8)

    where $ {Z}_{in} $ denotes the occurring probability of the nucleotide i at the position j, and n is the length of the sequence. The position-specific nucleotide matrix is estimated by computing the position-specific nucleotide frequencies over all DHS sequences in the benchmark dataset. The nucleotide information entropy and the position information entropy are respectively defined by

    $ N{P}^{i} = {\sum }_{j = 1}^{n}-{Z}_{ij}\mathrm{l}\mathrm{o}\mathrm{g}\left({Z}_{ij}\right) $ (9)

    and

    $ P{P}^{j} = {\sum }_{i = 1}^{4}{-Z}_{ij}\mathrm{l}\mathrm{o}\mathrm{g}\left({Z}_{ij}\right) . $ (10)

    The lower the information entropy, the more certain the distribution of the nucleotides in the sequences. Tables 5 and 6 showed the nucleotide information entropy for all 14 the tissues and 4 developmental stages, respectively. The nucleotide information entropy of DHSs was lower than those of non-DHSs, indicating that the distribution of nucleotides of DHSs was more certain than those of non-DHSs. The nucleotide information entropy was tissue-specific and stage-specific. For example, the nucleotide information entropy for the Muller retina tissue was lower than those for the Kidney tissue, and the nucleotide information entropy for the Early-Fetal stage was lower than those for the Late-Fetal stage. This implied that the distribution of nucleotides was tissue-specific and stage-specific. Figures 7 and 8 showed the positions of entropy values for all 14 the tissues and all 4 stages, respectively. It can be obviously observed that the position information entropy in the range of 230 to 300 was far less than 2, implying that the distribution of nucleotides in these regions were not completely stochastic. In addition, at the nearby position of 100, the position information entropy for tissues was lower than 2, indicating that the distribution of nucleotides was more certain in these regions.

    Table 5.  Nucleotide information entropy in different tissues of the mouse genome.
    DATASETS (TISSUES) Information entropy
    A T C G
    POS NEG ALL POS NEG ALL POS NEG ALL POS NEG ALL
    Forebrain 7.782 7.997 7.922 7.793 8.001 7.914 7.782 7.999 7.909 7.789 7.997 7.925
    Midbrain 7.729 7.998 7.912 7.736 8.001 7.901 7.728 7.997 7.895 7.732 8.000 7.913
    Hindbrain 7.796 7.997 7.927 7.804 7.994 7.918 7.794 7.997 7.917 7.795 7.994 7.923
    Liver 7.775 8.001 7.918 7.769 7.998 7.908 7.758 7.998 7.904 7.768 8.000 7.915
    Lung 7.770 7.993 7.910 7.761 7.993 7.902 7.767 7.993 7.904 7.763 7.990 7.908
    Heart 7.842 7.992 7.935 7.837 7.994 7.924 7.831 7.992 7.919 7.844 7.993 7.936
    Kidney 7.732 7.990 7.895 7.724 7.994 7.894 7.725 7.988 7.890 7.727 7.988 7.891
    Limb 7.781 8.002 7.926 7.787 8.001 7.916 7.777 7.996 7.911 7.786 7.999 7.925
    Stomach 7.788 7.999 7.957 7.767 7.992 7.934 7.765 7.996 7.937 7.772 8.000 7.956
    Thymus 7.592 7.979 7.867 7.613 7.980 7.857 7.602 7.984 7.859 7.612 7.987 7.881
    Craniofacial 7.620 7.994 7.880 7.634 7.992 7.870 7.632 7.994 7.871 7.626 7.991 7.878
    Retina 7.714 7.984 7.891 7.705 7.986 7.882 7.718 7.986 7.886 7.717 7.982 7.891
    Muller retina 7.571 7.990 7.862 7.560 7.995 7.873 7.581 7.994 7.876 7.555 7.988 7.858
    Neural tube 7.770 8.001 7.908 7.695 7.997 7.890 7.694 7.994 7.888 7.690 7.994 7.901

     | Show Table
    DownLoad: CSV
    Table 6.  Nucleotide information entropy in different developmental stages of the mouse genome.
    DATASETS (STAGES) Information Entropy
    A T C G
    POS NEG ALL POS NEG ALL POS NEG ALL POS NEG ALL
    ESC 7.717 7.992 7.904 7.708 8.000 7.888 7.710 7.993 7.883 7.720 7.997 7.908
    Early-Fetal 7.830 7.995 7.934 7.828 7.998 7.924 7.821 7.996 7.920 7.832 7.996 7.935
    Late-Fetal 7.650 7.998 7.888 7.656 8.000 7.885 7.652 7.998 7.881 7.659 7.998 7.890
    Adult 7.763 7.992 7.908 7.754 7.994 7.900 7.757 7.993 7.900 7.755 7.992 7.905

     | Show Table
    DownLoad: CSV
    Figure 7.  Position information entropy for all tissues.
    Figure 8.  Position information entropy for all developmental stages.

    To facilitate better using DHS sites and non-DHS sites, we have provided a useful web server at http:/www.biolscience.cn/LangMoDHS/. The web server interface is illustrated in Figure 9. Users first select the organization or developmental stage that needs to be predicted. Then, users input the sequence into the inputting box in the FASTA format. Alternatively, users upload a sequence file in the FASTA format. Finally, by clicking the submit button, users will get the predictive result in a certain amount of time which is determined by the number of inputted sequences.

    Figure 9.  Web server page of the LangMoDHS.

    DHSs play a key role in the cellular process. Sequence motif of DHSs is complicated, and thus, identifying DHSs is a challenging task at present. We have presented a deep language model for detecting DHSs in mouse genome. Extensive experiments showed that the LangMoDHS is an effective and efficient method for detecting DHSs. However, LangMoDHS performed differently with tissues and developmental stages. The LangMoDHS performed best for the Heart tissue, where the AUROC, AUPRC, and F1-score values were 0.960, 0.966, and 0.875, respectively, while it performed worse for 2 tissues, i.e., Thymus and Stomach, where the minimum AUROC and the minimum AUPRC values were 0.770 and 0.637, respectively. The range between the maximum and the minimum AUPRC value was up to 0.329, indicating that the sequence motif of DHSs would vary with tissue. This analysis also indicated that LangMoDHS is tissue-specific and stage-specific to predict DHSs.

    It is desirable to develop a universal method which is able to detect DHSs in all tissues or species. However, there is a difference between tissues and species, so it is very difficult to develop such a universal method in practice. Like the iDHS-Deep [43], the LangMoDHS exhibited a certain ability to detect DHSs across tissues or developmental stages. The LangMoDHS achieved better or competitive performance across tissues and developmental stages, indicating that these tissues and stages would be of a similar mechanism of DHSs.

    As mentioned previously, there are many computational approaches to detect DHSs. Compared with the methods [25,26,27,28,29,30,31,33,34,35,37], the LangMoDHS is an end-to-end method which requires no artificial design of features. The iDHS-Deep [43] is a newly developed deep-learning-based method to predict DHSs. The iDHS-Deep consisted mainly of two CNN layers and LSTM. The LSTM was attached at the end of the second CNN layer. The main difference between the iDHS-Deep and the LangMoDHS is that the latter used CNNs and Bi-LSTM in a parallel manner. The CNNs and the Bi-LSTM capture different characterization of the DHSs sequences respectively. Therefore, using CNNs and Bi-LSTM in a parallel manner would be more helpful to accumulate different characterization than stacking CNNs and Bi-LSTM in order. This might be a reason why the LangMoDHS performed better than the iDHS-Deep for most tissues and stages. Another difference is that the LangMoDHS to use feed-forward attention to improve the representations captured by the Bi-LSTM. Although the LangMoDHS exhibited competitive performance, the interpretability needs improving.

    Due to limitations of the methods or techniques, precisely and high effectively identifying DHSs remains challenging. We have presented a deep learning-based language model for computationally predicting DHSs in mouse genomes. Our method achieved competitive performance with the state-of-the-art methods. We developed a user-friendly web server to facilitate the identification of DHSs. The LangMoDHS presented has certain ability to predict DHSs across tissues and across stages, and it is tissue-specific and stage-specific. The nucleotide distributions of DHSs in some regions, such as nearby the position of 100 and the range from 230 to 300, is more certain.

    This work was supported by the National Natural Science Foundation of China (62272310, 62162025), Hunan Provincial Natural Science Foundation of China (2022JJ50177, 2020JJ4034), Scientific Research Fund of Hunan Provincial Education Department (21A0466, 19A215), and Shaoyang University Innovation Foundation for Postgraduate (CX2022SY046).

    The authors declare there is no conflict of interest.

    [1] World Health Organization (2014) World Health Statistics 2014.
    [2] World Health Organization (2015) Global action plan on antimicrobial resistance.
    [3] Griffith M, Postelnick M, Scheetz M (2012) Antimicrobial stewardship programs: methods of operation and suggested outcomes. Expert Rev Anti-Infe 10: 63–73. doi: 10.1586/eri.11.153
    [4] Yu VL (2011) Guidelines for hospital-acquired pneumonia and health-care-associated pneumonia: a vulnerability, a pitfall, and a fatal flaw. Lancet Infect Dis 11: 248–252. doi: 10.1016/S1473-3099(11)70005-6
    [5] Goossens H (2009) Antibiotic consumption and link to resistance. Clin Microbiol Infec 15 3:12–15.
    [6] Pakyz AL, MacDougall C, Oinonen M, et al. (2008) Trends in antibacterial use in US academic health centers: 2002 to 2006. Arch Intern Med 168: 2254–2260. doi: 10.1001/archinte.168.20.2254
    [7] Tacconelli E (2009) Antimicrobial use: risk driver of multidrug resistant microorganisms in healthcare settings. Curr Opin Infect Dis 22: 352–358. doi: 10.1097/QCO.0b013e32832d52e0
    [8] Landers TF, Cohen B, Wittum TE, et al. (2012) A review of antibiotic use in food animals: perspective, policy, and potential. Public Health Rep 127: 4–22. doi: 10.1177/003335491212700103
    [9] Wegener HC (2012) Antibiotic resistance-Linking human and animal health, In: Improving food safety through a One Health approach, Washington: National Academy of Sciences, 331–349.
    [10] Centers for Disease Control and Prevention (CDC) (2013) Antibiotic resistance threats in the United States, 2013, U.S, Department of Health and Human Services, CS239559-B.
    [11] Maragakis LL, Perencevich EN, Cosgrove SE (2008) Clinical and economic burden of antimicrobial resistance. Expert Rev Anti-Infe 6: 751–763. doi: 10.1586/14787210.6.5.751
    [12] Filice GA, Nyman JA, Lexau C, et al. (2010) Excess costs and utilization associated with methicillin resistance for patients with Staphylococcus aureus infection. Infect Cont Hosp Ep 31: 365–373. doi: 10.1086/651094
    [13] Hübner C, Hübner NO, Hopert K, et al. (2014) Analysis of MRSA-attributed costs of hospitalized patients in Germany. Eur J Clin Microbiol 33: 1817–1822. doi: 10.1007/s10096-014-2131-x
    [14] Macedo-Viñas M, De Angelis G, Rohner P, et al. (2013) Burden of methicillin-resistant Staphylococcus aureus infections at a Swiss University hospital: excess length of stay and costs. J Hosp Infect 84: 132–137. doi: 10.1016/j.jhin.2013.02.015
    [15] Pakyz A, Powell JP, Harpe SE, et al. (2008) Diversity of antimicrobial use and resistance in 42 hospitals in the United States. Pharmacotherapy 28: 906–912. doi: 10.1592/phco.28.7.906
    [16] Sandiumenge A, Diaz E, Rodriguez A, et al. (2006) Impact of diversity of antibiotic use on the development of antimicrobial resistance. J Antimicrob Chemoth 57: 1197–1204. doi: 10.1093/jac/dkl097
    [17] Wood TK, Knabel SJ, Kwan BW (2013) Bacterial persister cell formation and dormancy. Appl Environ Microbiol 79: 7116–7121. doi: 10.1128/AEM.02636-13
    [18] Keren I, Kaldalu N, Spoering A, et al. (2004) Persister cells and tolerance to antimicrobials. FEMS Microbiol Lett 230: 13–18. doi: 10.1016/S0378-1097(03)00856-5
    [19] Coculescu BI (2009) Antimicrobial resistance induced by genetic changes. J Med Life 2: 114–123.
    [20] Martinez JL (2014) General principles of antibiotic resistance in bacteria. Drug Discov Today 11: 33–39. doi: 10.1016/j.ddtec.2014.02.001
    [21] Cox G, Wright GD (2013) Intrinsic antibiotic resistance: mechanisms, origins, challenges and solutions. Int J Med Microbiol 303: 287–292. doi: 10.1016/j.ijmm.2013.02.009
    [22] Fajardo A, Martinez-Martin N, Mercadillo M, et al. (2008) The neglected intrinsic resistome of bacterial pathogens. PLoS One 3: e1619. doi: 10.1371/journal.pone.0001619
    [23] Davies J, Davies D (2010) Origins and evolution of antibiotic resistance. Microbiol Mol Biol Rev 74: 417–433. doi: 10.1128/MMBR.00016-10
    [24] Reygaert WC (2009) Methicillin-resistant Staphylococcus aureus (MRSA): molecular aspects of antimicrobial resistance and virulence. Clin Lab Sci 22: 115–119.
    [25] Blázquez J, Couce A, Rodríguez-Beltrán J, et al. (2012) Antimicrobials as promoters of genetic variation. Curr Opin Microbiol 15: 561–569. doi: 10.1016/j.mib.2012.07.007
    [26] Chancey ST, Zähner D, Stephens DS (2012) Acquired inducible antimicrobial resistance in Gram-positive bacteria. Future Microbiol 7: 959–978. doi: 10.2217/fmb.12.63
    [27] Mahon CR, Lehman DC, Manuselis G (2014) Antimicrobial agent mechanisms of action and resistance, In: Textbook of Diagnostic Microbiology, St. Louis: Saunders, 254–273.
    [28] Blair JM, Richmond GE, Piddock LJ (2014) Multidrug efflux pumps in Gram-negative bacteria and their role in antibiotic resistance. Future Microbiol 9: 1165–1177. doi: 10.2217/fmb.14.66
    [29] Kumar A, Schweizer HP (2005) Bacterial resistance to antibiotics: active efflux and reduced uptake. Adv Drug Deliver Rev 57: 1486–1513. doi: 10.1016/j.addr.2005.04.004
    [30] Lambert PA (2002) Cellular impermeability and uptake of biocides and antibiotics in gram-positive bacteria and mycobacteria. J Appl Microbiol 92: 46S–54S. doi: 10.1046/j.1365-2672.92.5s1.7.x
    [31] Bébéar CM, Pereyre S (2005) Mechanisms of drug resistance in Mycoplasma pneumoniae. Curr Drug Targets 5: 263–271. doi: 10.2174/1568005054880109
    [32] Miller WR, Munita JM, Arias CA (2014) Mechanisms of antibiotic resistance in enterococci. Expert Rev Anti-Infe 12: 1221–1236. doi: 10.1586/14787210.2014.956092
    [33] Gill MJ, Simjee S, Al-Hattawi K, et al. (1998) Gonococcal resistance to β-lactams and tetracycline involves mutation in loop 3 of the porin encoded at the penB locus. Antimicrob Agents Ch 42: 2799–2803.
    [34] Cornaglia G, Mazzariol A, Fontana R, et al. (1996) Diffusion of carbapenems through the outer membrane of enterobacteriaceae and correlation of their activities with their periplasmic concentrations. Microb Drug Resist 2: 273–276. doi: 10.1089/mdr.1996.2.273
    [35] Chow JW, Shlaes DM (1991) Imipenem resistance associated with the loss of a 40 kDa outer membrane protein in Enterobacter aerogenes. J Antimicrob Chemoth 28: 499–504. doi: 10.1093/jac/28.4.499
    [36] Thiolas A, Bornet C, Davin-Régli A, et al. (2004) Resistance to imipenem, cefepime, and cefpirome associated with mutation in Omp36 osmoporin of Enterobacter aerogenes. Biochem Bioph Res Co 317: 851–856. doi: 10.1016/j.bbrc.2004.03.130
    [37] Mah TF (2012) Biofilm-specific antibiotic resistance. Future Microbiol 7: 1061–1072. doi: 10.2217/fmb.12.76
    [38] Soto SM (2013) Role of efflux pumps in the antibiotic resistance of bacteria embedded in a biofilm. Virulence 4: 223–229. doi: 10.4161/viru.23724
    [39] Van Acker H, Van Dijck P, Coenye T (2014) Molecular mechanisms of antimicrobial tolerance and resistance in bacterial and fungal biofilms. Trends Microbiol 22: 326–333. doi: 10.1016/j.tim.2014.02.001
    [40] Beceiro A, Tomás M, Bou G (2013) Antimicrobial resistance and virulence: a successful or deleterious association in the bacterial world? Clin Microbiol Rev 26: 185–230. doi: 10.1128/CMR.00059-12
    [41] Randall CP, Mariner KR, Chopra I, et al. (2013) The target of daptomycin is absent form Escherichia coli and other gram-negative pathogens. Antimicrob Agents Ch 57: 637–639. doi: 10.1128/AAC.02005-12
    [42] Yang SJ, Kreiswirth BN, Sakoulas G, et al. (2009) Enhanced expression of dltABCD is associated with development of daptomycin nonsusceptibility in a clinical endocarditis isolate of Staphylococcus aureus. J Infect Dis 200: 1916–1920. doi: 10.1086/648473
    [43] Mishra NN, Bayer AS, Weidenmaier C, et al. (2014) Phenotypic and genotypic characterization of daptomycin-resistant methicillin-resistant Staphylococcus aureus strains: relative roles of mprF and dlt operons. PLoS One 9: e107426. doi: 10.1371/journal.pone.0107426
    [44] Stefani S, Campanile F, Santagati M, et al. (2015) Insights and clinical perspectives of daptomycin resistance in Staphylococcus aureus: a review of the available evidence. Int J Antimicrob Agents 46: 278–289. doi: 10.1016/j.ijantimicag.2015.05.008
    [45] Kumar S, Mukherjee MM, Varela MF (2013) Modulation of bacterial multidrug resistance efflux pumps of the major facilitator superfamily. Int J Bacteriol.
    [46] Roberts MC (2003) Tetracycline therapy: update. Clin Infect Dis 36: 462–467. doi: 10.1086/367622
    [47] Roberts MC (2004) Resistance to macrolide, lincosamide, streptogramin, ketolide, and oxazolidinone antibiotics. Mol Biotechnol 28: 47–62. doi: 10.1385/MB:28:1:47
    [48] Hawkey PM (2003) Mechanisms of quinolone action and microbial response. J Antimicrob Chemoth 1: 28–35.
    [49] Redgrave LS, Sutton SB, Webber MA, et al. (2014) Fluoroquinolone resistance: mechanisms, impact on bacteria, and role in evolutionary success. Trends Microbiol 22: 438–445. doi: 10.1016/j.tim.2014.04.007
    [50] Huovinen P, Sundström L, Swedberg G, et al. (1995) Trimethoprim and sulfonamide resistance. Antimicrob Agents Ch 39: 279–289. doi: 10.1128/AAC.39.2.279
    [51] Vedantam G, Guay GG, Austria NE, et al. (1998) Characterization of mutations contributing to sulfathiazole resistance in Escherichia coli. Antimicrob Agents Ch 42: 88–93.
    [52] Blair JM, Webber MA, Baylay AJ, et al. (2015) Molecular mechanisms of antibiotic resistance. Nat Rev Microbiol 13: 42–51. doi: 10.1038/nrmicro3380
    [53] Ramirez MS, Tolmasky ME (2010) Aminoglycoside modifying enzymes. Drug Resist Update 13: 151–171. doi: 10.1016/j.drup.2010.08.003
    [54] Robicsek A, Strahilevitz J, Jacoby GA, et al. (2006) Fluoroquinolone-modifying enzyme: a new adaptation of a common aminoglycoside acetyltransferase. Nat Med 12: 83–88. doi: 10.1038/nm1347
    [55] Schwarz S, Kehrenberg C, Doublet B, et al. (2004) Molecular basis of bacterial resistance to chloramphenicol and florfenicol. FEMS Microbiol Rev 28: 519–542. doi: 10.1016/j.femsre.2004.04.001
    [56] Pfeifer Y, Cullik A, Witte W (2010) Resistance to cephalosporins and carbapenems in Gram-negative bacterial pathogens. Int J Med Microbiol 300: 371–379. doi: 10.1016/j.ijmm.2010.04.005
    [57] Bush K, Bradford PA (2016) β-Lactams and β-lactamase inhibitors: an overview. CSH Perspect Med 6: a02527.
    [58] Bush K, Jacoby GA (2010) Updated functional classification of β-lactamases. Antimicrob Agents Ch 54: 969–976. doi: 10.1128/AAC.01009-09
    [59] Schultsz C, Geerlings S (2012) Plasmid-mediated resistance in Enterobacteriaceae. Drugs 72: 1–16.
    [60] Bush K (2013) Proliferation and significance of clinically relevant β-lactamases. Ann NY Acad Sci 1277: 84–90. doi: 10.1111/nyas.12023
    [61] Reygaert WC (2013) Antimicrobial resistance mechanisms of Staphylococcus aureus, In: Microbial pathogens and strategies for combating them: science, technology and education, Spain: Formatex, 297–310.
    [62] Toth M, Antunes NT, Stewart NK, et al. (2016) Class D β-lactamases do exist in Gram-positive bacteria. Nat Chem Biol 12: 9–14. doi: 10.1038/nchembio.1950
    [63] Jacoby GA (2009) AmpC β-lactamases. Clin Microbiol Rev 22: 161–182. doi: 10.1128/CMR.00036-08
    [64] Thomson KS (2010) Extended-spectrum-β-lactamase, AmpC, and carbapenemase issues. J Clin Microbiol 48: 1019–1025. doi: 10.1128/JCM.00219-10
    [65] Lahlaoui H, Khalifa ABH, Mousa MB (2014) Epidemiology of Enterobacteriaceae producing CTX-M type extended spectrum β-lactamase (ESBL). Med Maladies Infect 44: 400–404. doi: 10.1016/j.medmal.2014.03.010
    [66] Bevan ER, Jones AM, Hawkey PM (2017) Global epidemiology of CTX-M β-lactamases: temporal and geographical shifts in genotype. J Antimicrob Chemoth 72: 2145–2155. doi: 10.1093/jac/dkx146
    [67] Bajaj P, Singh NS, Virdi JS (2016) Escherichia coli β-lactamases: what really matters. Front Microbiol 7: 417.
    [68] Friedman ND, Tomkin E, Carmeli Y (2016) The negative impact of antibiotic resistance. Clin Microbiol Infect 22: 416–422. doi: 10.1016/j.cmi.2015.12.002
    [69] Zhanel GG, Lawson CD, Adam H, et al. (2013) Ceftazidime-Avibactam: a novel cephalosporin/β-lactamase inhibitor combination. Drugs 73: 159–177. doi: 10.1007/s40265-013-0013-7
    [70] Bush K (2018) Game changers: new β-lactamase inhibitor combinations targeting antibiotic resistance in gram-negative bacteria. ACS Infect Dis 4: 84–87. doi: 10.1021/acsinfecdis.7b00243
    [71] Docquier JD, Mangani S (2018) An update on β-lactamase inhibitor discovery and development. Drug Resist Update 36: 13–29. doi: 10.1016/j.drup.2017.11.002
    [72] Villagra NA, Fuentes JA, Jofré MR, et al. (2012) The carbon source influences the efflux pump-mediated antimicrobial resistance in clinically important Gram-negative bacteria. J Antimicrob Chemoth 67: 921–927. doi: 10.1093/jac/dkr573
    [73] Piddock LJ (2006) Clinically relevant chromosomally encoded multidrug resistance efflux pumps in bacteria. Clin Microbiol Rev 19: 382–402. doi: 10.1128/CMR.19.2.382-402.2006
    [74] Poole K (2007) Efflux pumps as antimicrobial resistance mechanisms. Ann Med 39: 162–176. doi: 10.1080/07853890701195262
    [75] Tanabe M, Szakonyi G, Brown KA, et al. (2009) The multidrug resistance efflux complex, EmrAB from Escherichia coli forms a dimer in vitro. Biochem Bioph Res Co 380: 338–342. doi: 10.1016/j.bbrc.2009.01.081
    [76] Jo I, Hong S, Lee M, et al. (2017) Stoichiometry and mechanistic implications of the MacAB-TolC tripartite efflux pump. Biochem Bioph Res Co 494: 668–673. doi: 10.1016/j.bbrc.2017.10.102
    [77] Jonas BM, Murray BE, Weinstock GM (2001) Characterization of emeA, a norA homolog and multidrug resistance efflux pump, in Enterococcus faecalis. Antimicrob Agents Ch 45: 3574–3579. doi: 10.1128/AAC.45.12.3574-3579.2001
    [78] Truong-Bolduc QC, Dunman PM, Strahilevitz J, et al. (2005) MgrA is a multiple regulator of two new efflux pumps in Staphylococcus aureus. J Bacteriol 187: 2395–2405. doi: 10.1128/JB.187.7.2395-2405.2005
    [79] Kourtesi C, Ball AR, Huang YY, et al. (2013) Microbial efflux systems and inhibitors: approaches to drug discovery and the challenge of clinical implementation. Open Microbiol J 7: 34–52. doi: 10.2174/1874285801307010034
    [80] Costa SS, Viveiros M, Amaral L, et al. (2013) Multidrug efflux pumps in Staphylococcus aureus: an update. Open Microbiol J 7: 59–71. doi: 10.2174/1874285801307010059
    [81] Lubelski J, Konings WN, Driessen AJ (2007) Distribution and physiology of ABC-type transporters contributing to multidrug resistance in bacteria. Microbiol Mol Biol Rev 71: 463–476. doi: 10.1128/MMBR.00001-07
    [82] Putman M, van Veen HW, Konings WN (2000) Molecular properties of bacterial multidrug transporters. Microbiol Mol Biol Rev 64: 672–693. doi: 10.1128/MMBR.64.4.672-693.2000
    [83] Kuroda T, Tsuchiya T (2009) Multidrug efflux transporters in the MATE family. BBA-Proteins Proteom 1794: 763–768. doi: 10.1016/j.bbapap.2008.11.012
    [84] Rouquette-Loughlin, C, Dunham SA, Kuhn M, et al. (2003) The NorM efflux pump of Neisseria gonorrhoeae and Neisseria meningitidis recognizes antimicrobial cationic compounds. J Bacteriol 185: 1101–1106. doi: 10.1128/JB.185.3.1101-1106.2003
    [85] Bay DC, Rommens KL, Turner RJ (2008) Small multidrug resistance proteins: a multidrug transporter family that continues to grow. BBA-Biomembranes 1778: 1814–1838. doi: 10.1016/j.bbamem.2007.08.015
    [86] Yerushalmi H, Lebendiker M, Schuldiner S (1995) EmrE, an Escherichia coli 12-kDa multidrug transporter, exchanges toxic cations and H+ and is soluble in organic solvents. J Biol Chem 270: 6856–6863. doi: 10.1074/jbc.270.12.6856
    [87] Collu F, Cascella M (2013) Multidrug resistance and efflux pumps: insights from molecular dynamics simulations. Curr Top Med Chem 13: 3165–3183. doi: 10.2174/15680266113136660224
    [88] Martinez JL, Sánchez MB, Martinez-Solano L, et al. (2009) Functional role of bacterial multidrug efflux pumps in microbial natural ecosystems. FEMS Microbiol Rev 33: 430–449. doi: 10.1111/j.1574-6976.2008.00157.x
    [89] Deak D, Outterson K, Powers JH, et al. (2016) Progress in the fight against multidrug-resistant bacteria? A review of U.S. Food and Drug Administration-approved antibiotics, 2010-2015. Ann Intern Med 165: 363–372.
  • This article has been cited by:

    1. Sarangam Majumdar, Sukla Pal, Information transmission in microbial and fungal communication: from classical to quantum, 2018, 12, 1873-9601, 491, 10.1007/s12079-018-0462-6
    2. H.J. Eberl, E.M. Jalbert, A. Dumitrache, G.M. Wolfaardt, A spatially explicit model of inverse colony formation of cellulolytic biofilms, 2017, 122, 1369703X, 141, 10.1016/j.bej.2017.03.007
    3. B. D’Acunto, L. Frunzo, I. Klapper, M.R. Mattei, P. Stoodley, Mathematical modeling of dispersal phenomenon in biofilms, 2019, 307, 00255564, 70, 10.1016/j.mbs.2018.07.009
    4. M. R. Mattei, L. Frunzo, B. D’Acunto, Y. Pechaud, F. Pirozzi, G. Esposito, Continuum and discrete approach in modeling biofilm development and structure: a review, 2018, 76, 0303-6812, 945, 10.1007/s00285-017-1165-y
    5. Liu Feng, Dai Xiangjuan, Mei Qicheng, Ling Guang, Wang Xinmei, 2018, Bifurcation Analysis and Control of Fractional-Order Quorum Sensing Network, 978-988-15639-5-8, 10197, 10.23919/ChiCC.2018.8483057
    6. Jia Zhao, Qi Wang, Three-Dimensional Numerical Simulations of Biofilm Dynamics with Quorum Sensing in a Flow Cell, 2017, 79, 0092-8240, 884, 10.1007/s11538-017-0259-4
    7. Hana Ueda, Kristina Stephens, Konstantina Trivisa, William E. Bentley, Sang Yup Lee, Bacteria Floc, but Do They Flock? Insights from Population Interaction Models of Quorum Sensing, 2019, 10, 2150-7511, 10.1128/mBio.00972-19
    8. Pavel Zarva, Hermann J. Eberl, 2020, Chapter 17, 978-3-030-50435-9, 228, 10.1007/978-3-030-50436-6_17
    9. Estefanía Garibay‐Valdez, Luis Rafael Martínez‐Córdova, Francisco Vargas‐Albores, Maurício G. C. Emerenciano, Anselmo Miranda‐Baeza, Edilmar Cortés‐Jacinto, Ángel M. Ortiz‐Estrada, Francesco Cicala, Marcel Martínez‐Porchas, The biofouling process: The science behind a valuable phenomenon for aquaculture, 2022, 1753-5123, 10.1111/raq.12770
    10. Christoph Helmer, Ansgar Jüngel, Antoine Zurek, Analysis of a finite-volume scheme for a single-species biofilm model, 2023, 185, 01689274, 386, 10.1016/j.apnum.2022.12.002
    11. Blessing O. Emerenini, Hermann J. Eberl, Reactor scale modeling of quorum sensing induced biofilm dispersal, 2022, 418, 00963003, 126792, 10.1016/j.amc.2021.126792
    12. Christoph Helmer, Ansgar Jüngel, Existence analysis for a reaction-diffusion Cahn–Hilliard-type system with degenerate mobility and singular potential modeling biofilm growth, 2023, 0, 1078-0947, 0, 10.3934/dcds.2023069
  • Reader Comments
  • © 2018 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(53331) PDF downloads(4675) Cited by(1223)

Figures and Tables

Figures(3)  /  Tables(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog