Direct Selection Pressure Research Articles

The significance of environmentally bioavailable antimicrobials in driving antimicrobial resistance in soils

Antimicrobial resistance (AMR) stands as an escalating public health crisis fueled by antimicrobial residues in the environment, particularly in soil, which acts as a reservoir for antimicrobial resistance genes (ARGs). Merely quantifying the total extractable concentration of antimicrobials, instead of bioavailable fractions, may substantially underestimate their minimal selection concentration for propagating ARGs. To shed light on the role of bioavailability in ARG abundance within soil, a systematic bioavailability assessment method was established for accurately quantifying the partitioning of multi-class antimicrobials in representative Chinese soils. Microcosm studies unveiled that antimicrobials persisting in the bioavailable fraction could potentially prolong their selection pressure duration to trigger AMR. Notably, the co-occurrence of pesticide or steroid hormone influenced the development trends of ARG subtypes, with fluoroquinolone resistance genes (RGs) being particularly susceptible. Partial least squares path model (PLS-PM) analysis uncovered potentially distinct induction mechanisms of antimicrobials: observable results suggested that extractable residual concentration may exert a direct selection pressure on the development of ARGs, while bioavailable concentration could potentially play a stepwise role in affecting the abundance of mobile genetic elements and initiating ARG dissemination. Such unprecedented scrutinization of the interplay between bioavailable antimicrobials in soils and ARG abundance provides valuable insights into strategizing regulatory policy or guidelines for soil remediation.

Evolution of hierarchy and irreversibility in theoretical cell differentiation model.

The process of cell differentiation in multicellular organisms is characterized by hierarchy and irreversibility in many cases. However, the conditions and selection pressures that give rise to these characteristics remain poorly understood. By using a mathematical model, here we show that the network of differentiation potency (differentiation diagram) becomes necessarily hierarchical and irreversible by increasing the number of terminally differentiated states under certain conditions. The mechanisms generating these characteristics are clarified using geometry in the cell state space. The results demonstrate that the hierarchical organization and irreversibility can manifest independently of direct selection pressures associated with these characteristics, instead they appear to evolve as byproducts of selective forces favoring a diversity of differentiated cell types. The study also provides a new perspective on the structure of gene regulatory networks that produce hierarchical and irreversible differentiation diagrams. These results indicate some constraints on cell differentiation, which are expected to provide a starting point for theoretical discussion of the implicit limits and directions of evolution in multicellular organisms.

Homozygosity mapping in the Kazakh national dog breed Tazy

The Tazy is a breed of sighthound common in Kazakhstan. The identification of runs of homozygosity (ROH) is an informative approach to assessing the history and possible patterns of directional selection pressure. To our knowledge, the present study is the first to provide an overview of the ROH pattern in the Tazy dogs from a genome-wide perspective. The ROH of the Tazy was found to be mainly composed of shorter segments (1–2 Mb), accounting for approximately 67% of the total ROH. The estimated ROH-based inbreeding coefficients (FROH) ranged from 0.028 to 0.058 with a mean of 0.057. Five genomic regions under positive selection were identified on chromosomes 18, 22, and 25. The regions on chromosomes 18 and 22 may be breed specific, while the region on chromosome 22 overlaps with regions of hunting traits in other hunting dog breeds. Among the 12 candidate genes located in these regions, the gene CAB39L may be a candidate that affects running speed and endurance of the Tazy dog. Eight genes could belong to an evolutionarily conserved complex as they were clustered in a large protein network with strong linkages. The results may enable effective interventions when incorporated into conservation planning and selection of the Tazy breed.

Ancistroceruswasps (Hymenoptera, Vespidae) from the Centre of Europe: Phylogeny, cryptic species, neutral and non‐neutral markers

AbstractThe aim of the study was to clarify the phylogenetic relationships among Northern EuropeanAncistrocerusand comparison of the applicability of evolutionarily neutral and non‐neutral markers for reconstruction of phylogeny. We used a 19,400 bp long dataset that included parts of mitochondrial DNA, nuclear rDNA operon, and 10 nuclear protein‐coding genes. Application of molecular barcoding and species delimitation algorithms unveiled a presence of cryptic species,A. balticussp. n., in the trap‐nesting wasp communities of the centre of Europe. We assessed the morphological, biological, and ecological differences of it from the siblingA. trifasciatusand updated the regional identification key. Phylogeny reconstruction using the neutral and the presumably non‐neutral markers resulted in different tree topologies. Evolutionary congruence of the rDNA operon with the other markers was relatively low. Evolutionary rate of the mitochondrial genes was 7–8 times as high as that of the exons of the nuclear genes, therefore, the mitochondrial markers overshadowed the nuclear ones in the phylogeny reconstructions. We assumed that at the speciation level, we might consider two different patterns of phylogeny: one based on evolutionary time and neutral changes, and the other based on adaptive evolutionary pathways under directional selection pressures. We assessed the effect of directional selection on the nuclear protein‐coding genes, applying the Spearman's rank correlation between pairwise phylogenetic distances among species, estimated using exons, and these distances, estimated using introns. One of these markers demonstrated a lack of positive correlation, implying a variable directional selection pressure on the coded protein. The publication has been registered on ZooBank: urn:lsid:zoobank.org:pub:13BD28D0‐736D‐4B2A‐B5CF‐4824BD4CDCFB.

Gene network robustness as a multivariate character

Robustness to genetic or environmental disturbances is often considered as a key property of living systems. Yet, in spite of being discussed since the 1950s, how robustness emerges from the complexity of genetic architectures and how it evolves still remains unclear. In particular, whether or not robustness to various sources of perturbations is independent conditions the range of adaptive scenarios that can be considered. For instance, selection for robustness to heritable mutations is likely to be modest and indirect, and its evolution might result from indirect selection on a pleiotropically-related character (e.g., homeostasis) rather than adaptation. Here, I propose to treat various robustness measurements as quantitative characters, and study theoretically, by individual-based simulations, their propensity to evolve independently. Based on a simple evolutionary model of a gene regulatory network, I showed that different ways to measure the robustness of gene expression to genetic or non-genetic disturbances were substantially correlated. Yet, robustness was evolvable in several dimensions, and robustness components could evolve differentially under direct selection pressure. Therefore, the fact that the sensitivity of gene expression to e.g. mutations and environmental factors rely on the same gene networks does not preclude that robustness components may have distinct evolutionary histories.

The Influencing Factors of Bacterial Resistance Related to Livestock Farm: Sources and Mechanisms

Bacterial resistance is a complex scientific issue. To manage this issue, we need to deeply understand the influencing factors and mechanisms. Based on the background of livestock husbandry, this paper reviews the factors that affect the acquisition of bacterial resistance. Meanwhile, the resistance mechanism is also discussed. “Survival of the fittest” is the result of genetic plasticity of bacterial pathogens, which brings about specific response, such as producing adaptive mutation, gaining genetic material or changing gene expression. To a large extent, bacterial populations acquire resistance genes directly caused by the selective pressure of antibiotics. However, mobile resistance genes may be co-selected by other existing substances (such as heavy metals and biocides) without direct selection pressure from antibiotics. This is because the same mobile genetic elements as antibiotic resistance genes can be co-located by the resistance determinants of some of these compounds. Furthermore, environmental factors are a source of resistance gene acquisition. Here, we describe some of the key measures that should be taken to mitigate the risk of antibiotic resistance. We call on the relevant governments or organizations around the world to formulate and improve the monitoring policies of antibiotic resistance, strengthen the supervision, strengthen the international cooperation and exchange, and curb the emergence and spread of drug-resistant strains.

Adult resistance genes to barley powdery mildew confer basal penetration resistance associated with broad-spectrum resistance.

Powdery mildew isa major disease of barley (Hordeum vulgare L.) for which breeders have traditionally relied on dominant, pathogen race-specific resistance genes for genetic control. Directional selection pressures in extensive monocultures invariably result in such genes being overcome as the pathogen mutates to evade recognition. This has led to a widespread reliance on fungicides and a single broad-spectrum recessive resistance provided by the mlo gene. The range of resistance genes and alleles found in wild crop relatives and landraces has been reduced in agricultural cultivars through an erosion of genetic diversity during domestication and selective breeding. Three novel major-effect adult plant resistance (APR) genes from landraces, designated Resistance to Blumeria graminis f. sp. hordei (Rbgh1 to Rbgh3), were identified in the terminal regions of barley chromosomes 5HL, 7HS, and 1HS, respectively. The phenotype of the new APR genes showed neither pronounced penetration resistance, nor the spontaneous necrosis and mesophyll cell death typical of mlo resistance, nor a whole epidermal cell hypersensitive response, typical of race-specific resistance. Instead, resistance was localized to the site of attempted penetration in an epidermal cell and was associated with cell wall appositions and cytosolic vesicle-like bodies, and lacked strong induction of reactive oxygen species. The APR genes exhibited differences in vesicle-like body sizes, their distribution, and the extent of localized 3,3-diaminobenzidine staining in individual doubled haploid lines. The results revealed a set of unique basal penetration resistance genes that offer opportunities for combining different resistance mechanisms in breeding programs for robust mildew resistance.

Mobile Oxazolidinone Resistance Genes in Gram-Positive and Gram-Negative Bacteria.

Seven mobile oxazolidinone resistance genes, including cfr, cfr(B), cfr(C), cfr(D), cfr(E), optrA, and poxtA, have been identified to date. The cfr genes code for 23S rRNA methylases, which confer a multiresistance phenotype that includes resistance to phenicols, lincosamides, oxazolidinones, pleuromutilins, and streptogramin A compounds. The optrA and poxtA genes code for ABC-F proteins that protect the bacterial ribosomes from the inhibitory effects of oxazolidinones. The optrA gene confers resistance to oxazolidinones and phenicols, while the poxtA gene confers elevated MICs or resistance to oxazolidinones, phenicols, and tetracycline. These oxazolidinone resistance genes are most frequently found on plasmids, but they are also located on transposons, integrative and conjugative elements (ICEs), genomic islands, and prophages. In these mobile genetic elements (MGEs), insertion sequences (IS) most often flanked the cfr, optrA, and poxtA genes and were able to generate translocatable units (TUs) that comprise the oxazolidinone resistance genes and occasionally also other genes. MGEs and TUs play an important role in the dissemination of oxazolidinone resistance genes across strain, species, and genus boundaries. Most frequently, these MGEs also harbor genes that mediate resistance not only to antimicrobial agents of other classes, but also to metals and biocides. Direct selection pressure by the use of antimicrobial agents to which the oxazolidinone resistance genes confer resistance, but also indirect selection pressure by the use of antimicrobial agents, metals, or biocides (the respective resistance genes against which are colocated on cfr-, optrA-, or poxtA-carrying MGEs) may play a role in the coselection and persistence of oxazolidinone resistance genes.

Epigenomic modifications induced by hatchery rearing persist in germ line cells of adult salmon after their oceanic migration.

Human activities induce direct or indirect selection pressure on natural population and may ultimately affect population's integrity. While numerous conservation programs aimed to minimize human‐induced genomic variation, human‐induced environmental variation may generate epigenomic variation potentially affecting fitness through phenotypic modifications. Major questions remain pertaining to how much epigenomic variation arises from environmental heterogeneity, whether this variation can persist throughout life, and whether it can be transmitted across generations. We performed whole genome bisulfite sequencing (WGBS) on the sperm of genetically indistinguishable hatchery and wild‐born migrating adults of Coho salmon (Oncorhynchus kisutch) from two geographically distant rivers at different epigenome scales. Our results showed that coupling WGBS with fine‐scale analyses (local and chromosomal) allows the detection of parallel early‐life hatchery‐induced epimarks that differentiate wild from hatchery‐reared salmon. Four chromosomes and 183 differentially methylated regions (DMRs) displayed a significant signal of methylation differentiation between hatchery and wild‐born Coho salmon. Moreover, those early‐life epimarks persisted in germ line cells despite about 1.5 year spent in the ocean following release from hatchery, opening the possibility for transgenerational inheritance. Our results strengthen the hypothesis that epigenomic modifications environmentally induced during early‐life development persist in germ cells of adults until reproduction, which could potentially impact their fitness.

Characteristics of Salmonella From Chinese Native Chicken Breeds Fed on Conventional or Antibiotic-Free Diets.

Salmonella is a common food-borne Gram-negative pathogen with multiple serotypes. Pullorum disease, caused by Salmonella Pullorum, seriously threatens the poultry industry. Many previous studies were focused on the epidemiological characteristics of Salmonella infections in conventional antibiotic use poultry. However, little is known about Salmonella infections in chicken flocks fed on antibiotic-free diets. Herein, we investigated and compared Salmonella infections in three Chinese native breeders fed on antibiotic-free diets, including the Luhua, Langya, and Qingjiaoma chickens, and one conventional breeder, the Bairi chicken, via analyzing 360 dead embryos in 2019. The results showed that the main Salmonella serotypes detected in a total of 155 isolates were S. Pullorum (82.6%) and S. Enteritidis (17.4%). Coinfection with two serotypes of Salmonella was specifically found in Bairi chicken. The sequence type (ST) in S. Pullorum was ST92 (n = 96) and ST2151 (n = 32), whereas only ST11 (n = 27) was found in S. Enteritidis. The Salmonella isolates from three breeder flocks fed on antibiotic-free diets exhibited phenotypic heterogeneity with a great variety of drug resistance spectrum. Most of the isolates among three chicken breeds Luhua (64.9%, 50/77), Langya (60%, 12/20) and Qingjiaoma (58.3%, 7/12) fed on antibiotic-free diets were resistant to only one antibiotic (erythromycin), whereas the rate of resistance to one antibiotic in conventional Bairi chicken isolates was only 4.3% (2/46). The multidrug-resistance rate in Salmonella isolates from layer flocks fed on antibiotic-free diets (20.2%, 22/109) was significantly (P < 0.0001) lower than that from chickens fed on conventional diets (93.5%, 43/46). However, high rate of resistance to erythromycin (97.4%~100%) and streptomycin (26%~41.7%) were also found among three breeder flocks fed on antibiotic-free diets, indicating resistance to these antibiotics likely spread before antibiotic-free feeding in poultry farms. The findings of this study supplement the epidemiological data of salmonellosis and provide an example of the characteristics of Salmonella in the chicken flocks without direct antibiotic selective pressure.

Long-term changes of plumage between urban and rural populations of white-crowned sparrows (Zonotrichia leucophrys)

Abstract Urbanization is one of the most extreme forms of land transformation and results in changes to ecosystems and species compositions. As a result, there are strong directional selection pressures compared to nearby rural areas. Despite a surge in research on the different selection pressures on acoustic communication in urban and rural areas, there has been comparatively little investigation into traits involved with visual communication. We measured the plumage of museum specimens of white-crowned sparrows (Zonotrichia leucophrys) from urban and adjacent rural habitats in San Francisco, CA, to assess the effects of divergent habitats on plumage. We found significant differences in dorsal plumage, but not crown plumage, between urban and rural populations that have been diverging over the past 100 years. Urban birds have increasingly darker and duller dorsal plumage, whereas rural birds in adjacent areas have plumage with richer hues and more color complexity. Our findings suggest a newly observed adaptation to urban environments by native species and suggest that many traits, in addition to acoustic signals, may be changing in response to urban selection pressures. Additional collections in urban areas are needed to explore likely divergences in plumage coloration between urban and rural environments.

Lack of phenological shift leads to increased camouflage mismatch in mountain hares.

Understanding whether organisms will be able to adapt to human-induced stressors currently endangering their existence is an urgent priority. Globally, multiple species moult from a dark summer to white winter coat to maintain camouflage against snowy landscapes. Decreasing snow cover duration owing to climate change is increasing mismatch in seasonal camouflage. To directly test for adaptive responses to recent changes in snow cover, we repeated historical (1950s) field studies of moult phenology in mountain hares (Lepus timidus) in Scotland. We found little evidence that population moult phenology has shifted to align seasonal coat colour with shorter snow seasons, or that phenotypic plasticity prevented increases in camouflage mismatch. The lack of responses resulted in 35 additional days of mismatch between 1950 and 2016. We emphasize the potential role of weak directional selection pressure and low genetic variability in shaping the scope for adaptive responses to anthropogenic stressors.

Do Non-Human Primates Really Represent Others' Beliefs?

Over two decades of research have produced compelling evidence that non-human primates understand some psychological states in other individuals but are unable to represent others' beliefs. Recently, three studies employing anticipatory looking (AL) paradigms reported that non-human primates do show hints of implicitly understanding the beliefs of others. However, measures of AL have been increasingly scrutinized in the human literature owing to extensive replication problems. We argue that new reports of belief representation in non-human primates using AL should be interpreted cautiously because of methodological and theoretical challenges paralleling trends in the human literature. We explore how future work can address these challenges, and conclude by identifying new evolutionary questions raised by the prospect that non-human primates implicitly represent others' beliefs without an explicit belief representation system that guides fitness-relevant behavior.

Loss of consumers constrains phenotypic evolution in the resulting food web.

The loss of biodiversity is altering the structure of ecological networks; however, we are currently in a poor position to predict how these altered communities will affect the evolution of remaining populations. Theory on fitness landscapes provides a framework for predicting how selection alters the evolutionary trajectory and adaptive potential of populations, but often treats the network of interacting populations as a “black box.” Here, we integrate ecological networks and fitness landscapes to examine how changes in food‐web structure shape phenotypic evolution. We conducted a field experiment that removed a guild of larval parasitoids that imposed direct and indirect selection pressures on an insect herbivore. We then measured herbivore survival as a function of three key phenotypic traits to estimate directional, quadratic, and correlational selection gradients in each treatment. We used these selection gradients to characterize the slope and curvature of the fitness landscape to understand the direct and indirect effects of consumer loss on phenotypic evolution. We found that the number of traits under directional selection increased with the removal of larval parasitoids, indicating evolution was more constrained toward a specific combination of traits. Similarly, we found that the removal of larval parasitoids altered the curvature of the fitness landscape in such a way that tended to decrease the evolvability of the traits we measured in the next generation. Our results suggest that the loss of trophic interactions can impose greater constraints on phenotypic evolution. This indicates that the simplification of ecological communities may constrain the adaptive potential of remaining populations to future environmental change.

Introgression of Two Quantitative Trait Loci for Stripe Rust Resistance into Three Chinese Wheat Cultivars

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most devastating diseases in wheat. Due to the large-scale and widely-distributed planting pattern of wheat, the directional selection pressure of the pathogen is very strong. Therefore, it is urgent to pyramid more stripe rust resistance genes in wheat cultivars to enhance resistance durability and ensure wheat production safety. In this study, two quantitative trait loci (QTL) for adult plant resistance (APR) to stripe rust, QYr.nafu-2BL and QYr.nafu-3BS, were validated and introgressed from wheat line P9897 into three Chinese elite wheat cultivars, Chuanmai 42, Xiangmai 25, and Zhengmai 9023, through marker validation. The three Chinese elite varieties were used as the female parent to cross with wheat line P9897, and they were selfed to the F6 generation. A total of 114 lines were then selected based on field agronomic traits and stripe rust resistance. Four markers (Xcfd73, Xgwm120, Xbarc87 and Xbarc133) linked with the QTL’s regions were employed to screen the 114 F6 lines. Subsequently, 27 lines combining two target QTL from P9897 were selected. The combination of agronomic traits and disease resistance results showed that 13 of these selected lines had favorable application prospects. The promising lines selected in this study could enrich the genetic resources of wheat stripe rust resistance genes, as well as provide material support and a theoretical basis for the prevention and control of wheat stripe rust in China.

Antibiotic Resistance Genes in Antibiotic-Free Chicken Farms.

Rising concern about the use of antibiotics in food production has resulted in many studies on the occurrence of antibiotic resistance genes (ARGs) in animal-associated bacterial communities. There are few baseline data on the abundance of ARGs on farms where chickens are intensively raised with little or no use of antibiotics. This study used a high-throughput quantitative PCR array to survey two antibiotic-free chicken farms for the occurrence of ARGs and mobile genetic elements known to enhance the spread of ARGs. No antibiotics had been used on the study farms for five years prior to this study. The results provide a baseline for the occurrence of resistance genes in the chicken production system without direct selective pressure.

Directed Evolution of the Stringency of the LuxR Vibrio fischeri Quorum Sensor without OFF-State Selection.

Stringency (low leak) is one of the most important specifications required for genetic circuits and induction systems, but it is challenging to evolve without sacrificing the maximum output level. This problem also comes from the absence of truly tunable negative selection methods. This paper reports that stringently switching variants can sometimes emerge with surprising frequency upon mutations. We randomly mutated the previously generated leaky variants of LuxR, the quorum-sensing transcription activator from Vibrio fischeri, to restore the stringency. We found as much as 10-20% of the entire population exhibited significantly improved signal-to-noise ratios compared with their parents. This indicated that these mutants arose by the loss of folding capability by accumulating destabilizing mutations, not by introducing rare adaptive mutations, thereby becoming AHL-dependent folders. Only four rounds of mutagenesis and ON-state selection resulted in the domination of the entire population by the improved variants with low leak, without direct selection pressure for stringency. With this surprising frequency, conversion into the "ligand-addicted folders" should be one of the prevailing modes of evolving stringency both in the laboratory and in nature, and the workflow described here provides a rapid and versatile method of improving the signal-to-noise ratio of various genetic switches.

Negative frequency dependent selection on plasmid carriage and low fitness costs maintain extended spectrum \u03b2-lactamases in Escherichia coli

Plasmids may maintain antibiotic resistance genes in bacterial populations through conjugation, in the absence of direct selection pressure. However, the costs and benefits of conjugation for plasmid and bacterial fitness are not well understood. Using invasion and competition experiments with plasmid mutants we explicitly tested how conjugation contributes to the maintenance of a plasmid bearing a single extended-spectrum ß-lactamase (ESBL) gene (blaCTX-M-14). Surprisingly, conjugation had little impact on overall frequencies, although it imposed a substantial fitness cost. Instead, stability resulted from the plasmid conferring fitness benefits when rare. Frequency dependent fitness did not require a functional blaCTX-M-14 gene, and was independent of culture media. Fitness benefits when rare are associated with the core plasmid backbone but are able to drive up frequencies of antibiotic resistance because fitness burden of the blaCTX-M-14 gene is very low. Negative frequency dependent fitness can contribute to maintaining a stable frequency of resistance genes in the absence of selection pressure from antimicrobials. In addition, persistent, low cost resistance has broad implications for antimicrobial stewardship.

Antibiotic resistomes of healthy pig faecal metagenomes.

Antibiotic resistance reservoirs within food-producing animals are thought to be a risk to animal and human health. This study describes the minimum natural resistome of pig faeces as the bacteria are under no direct antibiotic selective pressure. The faecal resistome of 257 different genes comprised 56 core and 201 accessory resistance genes. The genes present at the highest relative abundances across all samples were tetW, tetQ, tet44, tet37, tet40, mefA, aadE, ant(9)−1, ermB and cfxA2. This study characterized the baseline resistome, the microbiome composition and the metabolic components described by the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in healthy pig faeces, without antibiotic selective pressures. The microbiome hierarchical analysis resulted in a cluster tree with a highly similar pattern to that of the accessory resistome cluster tree. Functional capacity profiling identified genes associated with horizontal gene transfer. We identified a statistically significant positive correlation between the total antibiotic resistome and suggested indicator genes, which agree with using these genes as indicators of the total resistomes. The correlation between total resistome and total microbiome in this study was positive and statistically significant. Therefore, the microbiome composition influenced the resistome composition. This study identified a core and accessory resistome present in a cohort of healthy pigs, in the same conditions without antibiotics. It highlights the presence of antibiotic resistance in the absence of antibiotic selective pressure and the variability between animals even under the same housing, food and living conditions. Antibiotic resistance will remain in the healthy pig gut even when antibiotics are not used. Therefore, the risk of antibiotic resistance transfer from animal faeces to human pathogens or the environment will remain in the absence of antibiotics.

Association of colistin residues and manure treatment with the abundance of mcr-1 gene in swine feedlots

BackgroundThe extensive use of colistin in swine production may have contributed to the recent emergence of corresponding mobile resistance gene mcr-1. The use of colistin as a feed additive was banned in China in April 2017. ObjectivesTo examine the occurrence of colistin and dissemination of mcr-1 in swine feedlots before and after the colistin ban and effects of different manure treatments. MethodsEnvironmental samples were collected from swine feedlots before (December 2016) and after (December 2017) the colistin ban. Colistin concentrations were determined by ultra-high performance liquid chromatography coupled to tandem mass spectrometry. The prevalence of mcr-1 were determined by quantitative PCR analysis, while bacterial community composition was investigated by 16S rRNA sequencing. ResultsIn 2016, colistin was detected in feed and fresh manure samples at 67 mg/kg and 17 mg/kg, respectively, but was absent from all samples in 2017. In 2016, the relative abundance of mcr-1 in fresh manure was lower than that in solid samples after natural drying, while a higher relative abundance was detected in fresh manure samples compared with biogas slurry samples. A strong correlation between colistin concentration and relative abundance of mcr-1 was observed in fresh manure. The samples collected in 2017 showed a lower relative abundance of mcr-1 compared with those collected in 2016. Bacterial community analysis showed that the abundance of Enterobacteriaceae, which act as a vehicle and reservoir of mcr-1, increased with natural dying but decreased with anaerobic digestion. ConclusionsThe presence of colistin exerts direct selection pressure for the accumulation of mcr-1 in manure, while the ban on colistin likely halted the dissemination of mcr-1 on pig farms. Anaerobic digestion is an effective waste treatment process for removing mcr-1, which might be mainly driven by the shift in bacterial community structure.