SequelTools: a suite of tools for working with PacBio Sequel raw sequence data.

BACKGROUND: PacBio sequencing is an incredibly valuable third-generation DNA sequencing method due to very long read lengths, ability to detect methylated bases, and its real-time sequencing methodology. Yet, hitherto no tool was available for analyzing the quality of, subsampling, and filtering PacBio data. RESULTS: Here we present SequelTools, a command-line program containing three tools: Quality Control, Read Subsampling, and Read Filtering. The Quality Control tool quickly processes PacBio Sequel raw sequence data from multiple SMRTcells producing multiple statistics and publication-quality plots describing the quality of the data including N50, read length and count statistics, PSR, and ZOR. The Read Subsampling tool allows the user to subsample reads by one or more of the following criteria: longest subreads per CLR or random CLR selection. The Read Filtering tool provides options for normalizing data by filtering out certain low-quality scraps reads and/or by minimum CLR length. SequelTools is implemented in bash, R, and Python using only standard libraries and packages and is platform independent. CONCLUSIONS: SequelTools is a program that provides the only free, fast, and easy-to-use quality control tool, and the only program providing this kind of read subsampling and read filtering for PacBio Sequel raw sequence data, and is available at https://github.com/ISUgenomics/SequelTools .

Thiol Starvation Induces Redox-Mediated Dysregulation of Escherichia coli Biofilm Components.

A hallmark of bacterial biofilms is the production of an extracellular matrix (ECM) that encases and protects the community from environmental stressors. Biofilm formation is an integral portion of the uropathogenic Escherichia coli (UPEC) life cycle. Approximately 2% of UPEC isolates are cysteine auxotrophs. Here, we investigated how cysteine homeostasis impacted UPEC UTI89 strain biofilm formation and, specifically, the production of the ECM components curli and cellulose. Cysteine auxotrophs produced less cellulose and slightly more curli compared to wild-type (WT) strains, and cysteine auxotrophs formed smooth, nonrugose colonies. Cellulose production was restored in cysteine auxotrophs when YfiR was inactivated. YfiR is a redox-sensitive regulator of the diguanylate cyclase, YfiN. The production of curli, a temperature-regulated appendage, was independent of temperature in UTI89 cysteine auxotrophs. In a screen of UPEC isolates, we found that ∼60% of UPEC cysteine auxotrophs produced curli at 37°C, but only ∼2% of cysteine prototrophic UPEC isolates produced curli at 37°C. Interestingly, sublethal concentrations of amdinocillin and trimethoprim-sulfamethoxazole inhibited curli production, whereas strains auxotrophic for cysteine continued to produce curli even in the presence of amdinocillin and trimethoprim-sulfamethoxazole. The dysregulation of ECM components and resistance to amdinocillin in cysteine auxotrophs may be linked to hyperoxidation, since the addition of exogenous cysteine or glutathione restored WT biofilm phenotypes to mutants unable to produce cysteine and glutathione.IMPORTANCE Uropathogenic Escherichia coli (UPEC) bacteria are the predominant causative agent of urinary tract infections (UTIs). UTIs account for billions of dollars of financial burden annually to the health care industry in the United States. Biofilms are an important aspect of the UPEC pathogenesis cascade and for the establishment of chronic infections. Approximately 2% of UPEC isolates from UTIs are cysteine auxotrophs, yet there is relatively little known about the biofilm formation of UPEC cysteine auxotrophs. Here we show that cysteine auxotrophs have dysregulated biofilm components due to a change in the redox state of the periplasm. Additionally, we show the relationship between cysteine auxotrophs, biofilms, and antibiotics frequently used to treat UTIs.

Consequences of Whole-Genome Triplication as Revealed by Comparative Genomic Analyses of the Wild Radish Raphanus raphanistrum and Three Other Brassicaceae Species.

Polyploidization events are frequent among flowering plants, and the duplicate genes produced via such events contribute significantly to plant evolution. We sequenced the genome of wild radish (Raphanus raphanistrum), a Brassicaceae species that experienced a whole-genome triplication event prior to diverging from Brassica rapa. Despite substantial gene gains in these two species compared with Arabidopsis thaliana and Arabidopsis lyrata, ∼70% of the orthologous groups experienced gene losses in R. raphanistrum and B. rapa, with most of the losses occurring prior to their divergence. The retained duplicates show substantial divergence in sequence and expression. Based on comparison of A. thaliana and R. raphanistrum ortholog floral expression levels, retained radish duplicates diverged primarily via maintenance of ancestral expression level in one copy and reduction of expression level in others. In addition, retained duplicates differed significantly from genes that reverted to singleton state in function, sequence composition, expression patterns, network connectivity, and rates of evolution. Using these properties, we established a statistical learning model for predicting whether a duplicate would be retained postpolyploidization. Overall, our study provides new insights into the processes of plant duplicate loss, retention, and functional divergence and highlights the need for further understanding factors controlling duplicate gene fate.

The Catabolite Repressor Protein-Cyclic AMP Complex Regulates csgD and Biofilm Formation in Uropathogenic Escherichia coli.

The extracellular matrix protects Escherichia coli from immune cells, oxidative stress, predation, and other environmental stresses. Production of the E. coli extracellular matrix is regulated by transcription factors that are tuned to environmental conditions. The biofilm master regulator protein CsgD upregulates curli and cellulose, the two major polymers in the extracellular matrix of uropathogenic E. coli (UPEC) biofilms. We found that cyclic AMP (cAMP) regulates curli, cellulose, and UPEC biofilms through csgD The alarmone cAMP is produced by adenylate cyclase (CyaA), and deletion of cyaA resulted in reduced extracellular matrix production and biofilm formation. The catabolite repressor protein (CRP) positively regulated csgD transcription, leading to curli and cellulose production in the UPEC isolate, UTI89. Glucose, a known inhibitor of CyaA activity, blocked extracellular matrix formation when added to the growth medium. The mutant strains ΔcyaA and Δcrp did not produce rugose biofilms, pellicles, curli, cellulose, or CsgD. Three putative CRP binding sites were identified within the csgD-csgB intergenic region, and purified CRP could gel shift the csgD-csgB intergenic region. Additionally, we found that CRP binded upstream of kpsMT, which encodes machinery for K1 capsule production. Together our work shows that cAMP and CRP influence E. coli biofilms through transcriptional regulation of csgD IMPORTANCE The catabolite repressor protein (CRP)-cyclic AMP (cAMP) complex influences the transcription of ∼7% of genes on the Escherichia coli chromosome (D. Zheng, C. Constantinidou, J. L. Hobman, and S. D. Minchin, Nucleic Acids Res 32:5874-5893, 2004, https://dx.doi.org/10.1093/nar/gkh908). Glucose inhibits E. coli biofilm formation, and ΔcyaA and Δcrp mutants show impaired biofilm formation (D. W. Jackson, J.W. Simecka, and T. Romeo, J Bacteriol 184:3406-3410, 2002, https://dx.doi.org/10.1128/JB.184.12.3406-3410.2002). We determined that the cAMP-CRP complex regulates curli and cellulose production and the formation of rugose and pellicle biofilms through csgD Additionally, we propose that cAMP may work as a signaling compound for uropathogenic E. coli (UPEC) to transition from the bladder lumen to inside epithelial cells for intracellular bacterial community formation through K1 capsule regulation.

Molecular Evidence for Functional Divergence and Decay of a Transcription Factor Derived from Whole-Genome Duplication in Arabidopsis thaliana.

Functional divergence between duplicate transcription factors (TFs) has been linked to critical events in the evolution of land plants and can result from changes in patterns of expression, binding site divergence, and/or interactions with other proteins. Although plant TFs tend to be retained post polyploidization, many are lost within tens to hundreds of million years. Thus, it can be hypothesized that some TFs in plant genomes are in the process of becoming pseudogenes. Here, we use a pair of salt tolerance-conferring transcription factors, DWARF AND DELAYED FLOWERING1 (DDF1) and DDF2, that duplicated through paleopolyploidy 50 to 65 million years ago, as examples to illustrate potential mechanisms leading to duplicate retention and loss. We found that the expression patterns of Arabidopsis thaliana (At)DDF1 and AtDDF2 have diverged in a highly asymmetric manner, and AtDDF2 has lost most inferred ancestral stress responses. Consistent with promoter disablement, the AtDDF2 promoter has fewer predicted cis-elements and a methylated repetitive element. Through comparisons of AtDDF1, AtDDF2, and their Arabidopsis lyrata orthologs, we identified significant differences in binding affinities and binding site preference. In particular, an AtDDF2-specific substitution within the DNA-binding domain significantly reduces binding affinity. Cross-species analyses indicate that both AtDDF1 and AtDDF2 are under selective constraint, but among A. thaliana accessions, AtDDF2 has a higher level of nonsynonymous nucleotide diversity compared with AtDDF1. This may be the result of selection in different environments or may point toward the possibility of ongoing functional decay despite retention for millions of years after gene duplication.

Iron induces bimodal population development by Escherichia coli.

Bacterial biofilm formation is a complex developmental process involving cellular differentiation and the formation of intricate 3D structures. Here we demonstrate that exposure to ferric chloride triggers rugose biofilm formation by the uropathogenic Escherichia coli strain UTI89 and by enteric bacteria Citrobacter koseri and Salmonella enterica serovar typhimurium. Two unique and separable cellular populations emerge in iron-triggered, rugose biofilms. Bacteria at the air-biofilm interface express high levels of the biofilm regulator csgD, the cellulose activator adrA, and the curli subunit operon csgBAC. Bacteria in the interior of rugose biofilms express low levels of csgD and undetectable levels of matrix components curli and cellulose. Iron activation of rugose biofilms is linked to oxidative stress. Superoxide generation, either through addition of phenazine methosulfate or by deletion of sodA and sodB, stimulates rugose biofilm formation in the absence of high iron. Additionally, overexpression of Mn-superoxide dismutase, which can mitigate iron-derived reactive oxygen stress, decreases biofilm formation in a WT strain upon iron exposure. Not only does reactive oxygen stress promote rugose biofilm formation, but bacteria in the rugose biofilms display increased resistance to H(2)O(2) toxicity. Altogether, we demonstrate that iron and superoxide stress trigger rugose biofilm formation in UTI89. Rugose biofilm development involves the elaboration of two distinct bacterial populations and increased resistance to oxidative stress.

Retained duplicate genes in green alga Chlamydomonas reinhardtii tend to be stress responsive and experience frequent response gains.

BACKGROUND: Green algae belong to a group of photosynthetic organisms that occupy diverse habitats, are closely related to land plants, and have been studied as sources of food and biofuel. Although multiple green algal genomes are available, a global comparative study of algal gene families has not been carried out. To investigate how gene families and gene expression have evolved, particularly in the context of stress response that have been shown to correlate with gene family expansion in multiple eukaryotes, we characterized the expansion patterns of gene families in nine green algal species, and examined evolution of stress response among gene duplicates in Chlamydomonas reinhardtii. RESULTS: Substantial variation in domain family sizes exists among green algal species. Lineage-specific expansion of families occurred throughout the green algal lineage but inferred gene losses occurred more often than gene gains, suggesting a continuous reduction of algal gene repertoire. Retained duplicates tend to be involved in stress response, similar to land plant species. However, stress responsive genes tend to be pseudogenized as well. When comparing ancestral and extant gene stress response state, we found that response gains occur in 13% of duplicate gene branches, much higher than 6% in Arabidopsis thaliana. CONCLUSION: The frequent gains of stress response among green algal duplicates potentially reflect a high rate of innovation, resulting in a species-specific gene repertoire that contributed to adaptive response to stress. This could be further explored towards deciphering the mechanism of stress response, and identifying suitable green algal species for oil production.

MAKER-P: a tool kit for the rapid creation, management, and quality control of plant genome annotations.

We have optimized and extended the widely used annotation engine MAKER in order to better support plant genome annotation efforts. New features include better parallelization for large repeat-rich plant genomes, noncoding RNA annotation capabilities, and support for pseudogene identification. We have benchmarked the resulting software tool kit, MAKER-P, using the Arabidopsis (Arabidopsis thaliana) and maize (Zea mays) genomes. Here, we demonstrate the ability of the MAKER-P tool kit to automatically update, extend, and revise the Arabidopsis annotations in light of newly available data and to annotate pseudogenes and noncoding RNAs absent from The Arabidopsis Informatics Resource 10 build. Our results demonstrate that MAKER-P can be used to manage and improve the annotations of even Arabidopsis, perhaps the best-annotated plant genome. We have also installed and benchmarked MAKER-P on the Texas Advanced Computing Center. We show that this public resource can de novo annotate the entire Arabidopsis and maize genomes in less than 3 h and produce annotations of comparable quality to those of the current The Arabidopsis Information Resource 10 and maize V2 annotation builds.

The disulfide bonding system suppresses CsgD-independent cellulose production in Escherichia coli.

The bacterial extracellular matrix encases cells and protects them from host-related and environmental insults. The Escherichia coli master biofilm regulator CsgD is required for the production of the matrix components curli and cellulose. CsgD activates the diguanylate cyclase AdrA, which in turn stimulates cellulose production through cyclic di-GMP (c-di-GMP). Here, we identified and characterized a CsgD- and AdrA-independent cellulose production pathway that was maximally active when cultures were grown under reducing conditions or when the disulfide bonding system (DSB) was compromised. The CsgD-independent cellulose activation pathway was dependent on a second diguanylate cyclase, called YfiN. c-di-GMP production by YfiN was repressed by the periplasmic protein YfiR, and deletion of yfiR promoted CsgD-independent cellulose production. Conversely, when YfiR was overexpressed, cellulose production was decreased. Finally, we found that YfiR was oxidized by DsbA and that intraprotein YfiR disulfide bonds stabilized YfiR in the periplasm. Altogether, we showed that reducing conditions and mutations in the DSB system caused hyperactivation of YfiN and subsequent CsgD-independent cellulose production.

Evolution of TPS20-related terpene synthases influences chemical diversity in the glandular trichomes of the wild tomato relative Solanum habrochaites.

A systematic screen of volatile terpene production in the glandular trichomes of 79 accessions of Solanum habrochaites was conducted and revealed the presence of 21 mono- and sesquiterpenes that exhibit a range of qualitative and quantitative variation. Hierarchical clustering identified distinct terpene phenotypic modules with shared patterns of terpene accumulation across accessions. Several terpene modules could be assigned to previously identified terpene synthase (TPS) activities that included members of the TPS-e/f subfamily that utilize the unusual cis-prenyl diphosphate substrates neryl diphosphate and 2z,6z-farnesyl diphosphate. DNA sequencing and in vitro enzyme activity analysis of TPS-e/f members from S. habrochaites identified three previously unassigned enzyme activities that utilize these cisoid substrates. These produce either the monoterpenes α-pinene and limonene, or the sesquiterpene 7-epizingiberene, with the in vitro analyses that recapitulated the trichome chemistry found in planta. Comparison of the distribution of S. habrochaites accessions with terpene content revealed a strong preference for the presence of particular TPS20 alleles at distinct geographic locations. This study reveals that the unusually high intra-specific variation of volatile terpene synthesis in glandular trichomes of S. habrochaites is due at least in part to evolution at the TPS20 locus.

Pre-existing heterogeneity facilitates development of heteroresistance upon gene acquisition.

Antibiotic resistance causes 1.27 million global deaths annually and is predicted to worsen. Heteroresistance is a form of resistance in which only a minor and unstable subpopulation of cells of a bacterial isolate are resistant to a given antibiotic, and are therefore often undetected by clinical diagnostics. These infrequent and undetected resistant cells can be selected during antibiotic therapy, expand in number, and cause unexplained treatment failures. A major question is how heteroresistance evolves. Here, studying the antibiotic fosfomycin, we report that heteroresistance can develop from a pre-existing state of phenotypic heterogeneity in which an isolate harbors a subpopulation with increased minimum inhibitory concentration (MIC), but below the clinical resistance breakpoint. We call this phenomenon heterosusceptibility and demonstrate that acquisition of a resistance gene, fosA, increases the MIC of the subpopulation beyond the breakpoint, making the isolate heteroresistant. Conversely, deletion of fosA from a heteroresistant isolate led to reduction of the MIC of the resistant subpopulation without a loss of heterogeneity, thus generating heterosusceptibility. A survey of 103 carbapenem-resistant Enterobacterales (CRE) revealed that the Escherichia sp. isolates lacked the fosA gene and uniformly exhibited fosfomycin heterosusceptibility, whereas the Klebsiella and Enterobacter encoded the fosA gene and were almost exclusively heteroresistant. Furthermore, some isolates exhibited heterosusceptibility to other antibiotics, demonstrating that this is a widespread phenomenon. These results highlight a mechanism for the evolution of heteroresistance and suggest that surveillance for heterosusceptibility may facilitate the prediction of impending heteroresistance before it evolves.

Swine-to-Ferret Transmission of Antigenically Drifted Contemporary Swine H3N2 Influenza A Virus Is an Indicator of Zoonotic Risk to Humans.

Human-to-swine transmission of influenza A (H3N2) virus occurs repeatedly and plays a critical role in swine influenza A virus (IAV) evolution and diversity. Human seasonal H3 IAVs were introduced from human-to-swine in the 1990s in the United States and classified as 1990.1 and 1990.4 lineages; the 1990.4 lineage diversified into 1990.4.A-F clades. Additional introductions occurred in the 2010s, establishing the 2010.1 and 2010.2 lineages. Human zoonotic cases with swine IAV, known as variant viruses, have occurred from the 1990.4 and 2010.1 lineages, highlighting a public health concern. If a variant virus is antigenically drifted from current human seasonal vaccine (HuVac) strains, it may be chosen as a candidate virus vaccine (CVV) for pandemic preparedness purposes. We assessed the zoonotic risk of US swine H3N2 strains by performing phylogenetic analyses of recent swine H3 strains to identify the major contemporary circulating genetic clades. Representatives were tested in hemagglutination inhibition assays with ferret post-infection antisera raised against existing CVVs or HuVac viruses. The 1990.1, 1990.4.A, and 1990.4.B.2 clade viruses displayed significant loss in cross-reactivity to CVV and HuVac antisera, and interspecies transmission potential was subsequently investigated in a pig-to-ferret transmission study. Strains from the three lineages were transmitted from pigs to ferrets via respiratory droplets, but there were differential shedding profiles. These data suggest that existing CVVs may offer limited protection against swine H3N2 infection, and that contemporary 1990.4.A viruses represent a specific concern given their widespread circulation among swine in the United States and association with multiple zoonotic cases.

Characterizing a century of genetic diversity and contemporary antigenic diversity of N1 neuraminidase in influenza A virus from North American swine.

Influenza A viruses (IAVs) of the H1N1 classical swine lineage became endemic in North American swine following the 1918 pandemic. Additional human-to-swine transmission events after 1918, and a spillover of H1 viruses from wild birds in Europe, potentiated a rapid increase in genomic diversity via reassortment between introductions and the endemic classical swine lineage. To determine mechanisms affecting reassortment and evolution, we conducted a phylogenetic analysis of N1 and paired HA swine IAV genes in North America between 1930 and 2020. We described fourteen N1 clades within the N1 Eurasian avian lineage (including the N1 pandemic clade), the N1 classical swine lineage, and the N1 human seasonal lineage. Seven N1 genetic clades had evidence for contemporary circulation. To assess antigenic drift associated with N1 genetic diversity, we generated a panel of representative swine N1 antisera and quantified the antigenic distance between wild-type viruses using enzyme-linked lectin assays and antigenic cartography. Within the N1 genes, the antigenic similarity was variable and reflected shared evolutionary history. Sustained circulation and evolution of N1 genes in swine had resulted in a significant antigenic distance between the N1 pandemic clade and the classical swine lineage. Between 2010 and 2020, N1 clades and N1-HA pairings fluctuated in detection frequency across North America, with hotspots of diversity generally appearing and disappearing within 2 years. We also identified frequent N1-HA reassortment events (n = 36), which were rarely sustained (n = 6) and sometimes also concomitant with the emergence of new N1 genetic clades (n = 3). These data form a baseline from which we can identify N1 clades that expand in range or genetic diversity that may impact viral phenotypes or vaccine immunity and subsequently the health of North American swine.

A Chromosome-Level Genome of the Camphor Tree and the Underlying Genetic and Climatic Factors for Its Top-Geoherbalism.

Camphor tree [Cinnamomum camphora (L.) J. Presl], a species in the magnoliid family Lauraceae, is known for its rich volatile oils and is used as a medical cardiotonic and as a scent in many perfumed hygiene products. Here, we present a high-quality chromosome-scale genome of C. camphora with a scaffold N50 of 64.34 Mb and an assembled genome size of 755.41 Mb. Phylogenetic inference revealed that the magnoliids are a sister group to the clade of eudicots and monocots. Comparative genomic analyses identified two rounds of ancient whole-genome duplication (WGD). Tandem duplicated genes exhibited a higher evolutionary rate, a more recent evolutionary history and a more clustered distribution on chromosomes, contributing to the production of secondary metabolites, especially monoterpenes and sesquiterpenes, which are the principal essential oil components. Three-dimensional analyses of the volatile metabolites, gene expression and climate data of samples with the same genotype grown in different locations showed that low temperature and low precipitation during the cold season modulate the expression of genes in the terpenoid biosynthesis pathways, especially TPS genes, which facilitates the accumulation of volatile compounds. Our study lays a theoretical foundation for policy-making regarding the agroforestry applications of camphor tree.

Interspecies Transmission from Pigs to Ferrets of Antigenically Distinct Swine H1 Influenza A Viruses with Reduced Reactivity to Candidate Vaccine Virus Antisera as Measures of Relative Zoonotic Risk.

During the last decade, endemic swine H1 influenza A viruses (IAV) from six different genetic clades of the hemagglutinin gene caused zoonotic infections in humans. The majority of zoonotic events with swine IAV were restricted to a single case with no subsequent transmission. However, repeated introduction of human-seasonal H1N1, continual reassortment between endemic swine IAV, and subsequent drift in the swine host resulted in highly diverse swine IAV with human-origin genes that may become a risk to the human population. To prepare for the potential of a future swine-origin IAV pandemic in humans, public health laboratories selected candidate vaccine viruses (CVV) for use as vaccine seed strains. To assess the pandemic risk of contemporary US swine H1N1 or H1N2 strains, we quantified the genetic diversity of swine H1 HA genes, and identified representative strains from each circulating clade. We then characterized the representative swine IAV against human seasonal vaccine and CVV strains using ferret antisera in hemagglutination inhibition assays (HI). HI assays revealed that 1A.3.3.2 (pdm09) and 1B.2.1 (delta-2) demonstrated strong cross reactivity to human seasonal vaccines or CVVs. However, swine IAV from three clades that represent more than 50% of the detected swine IAVs in the USA showed significant reduction in cross-reactivity compared to the closest CVV virus: 1A.1.1.3 (alpha-deletion), 1A.3.3.3-clade 3 (gamma), and 1B.2.2.1 (delta-1a). Representative viruses from these three clades were further characterized in a pig-to-ferret transmission model and shown to exhibit variable transmission efficiency. Our data prioritize specific genotypes of swine H1N1 and H1N2 to further investigate in the risk they pose to the human population.

1,3-Dioxane-Linked Novel Bacterial Topoisomerase Inhibitors: Expanding Structural Diversity and the Antibacterial Spectrum.

Antibacterial resistance continues its devastation of available therapies. Novel bacterial topoisomerase inhibitors (NBTIs) offer one solution to this critical issue. Two series of amine NBTIs bearing tricyclic DNA-binding moieties as well as amide NBTIs with a bicyclic DNA-binding moiety were synthesized and evaluated against methicillin-resistant Staphylococcus aureus (MRSA). Additionally, these compounds and a series of bicyclic amine analogues displayed high activity against susceptible and drug-resistant Neisseria gonorrhoeae, expanding the spectrum of these dioxane-linked NBTIs.

octoFLUshow: an Interactive Tool Describing Spatial and Temporal Trends in the Genetic Diversity of Influenza A Virus in U.S. Swine.

Influenza A virus (IAV) is passively surveilled in swine in the United States through a U.S. Department of Agriculture administered surveillance system. We present an interactive Web tool to visualize and explore trends in the genetic and geographic diversity of IAV derived from the surveillance system.

Antibiotic-Selected Gene Amplification Heightens Metal Resistance.

The increasing frequency of antibiotic resistance poses myriad challenges to modern medicine. Environmental survival of multidrug-resistant bacteria in health care facilities, including hospitals, creates reservoirs for transmission of these difficult to treat pathogens. To prevent bacterial colonization, these facilities deploy an array of infection control measures, including bactericidal metals on surfaces, as well as implanted devices. Although antibiotics are routinely used in these health care environments, it is unknown whether and how antibiotic exposure affects metal resistance. We identified a multidrug-resistant Enterobacter clinical isolate that displayed heteroresistance to the antibiotic colistin, where only a minor fraction of cells within the population resist the drug. When this isolate was grown in the presence of colistin, a 9-kb DNA region was duplicated in the surviving resistant subpopulation, but surprisingly, was not required for colistin heteroresistance. Instead, the amplified region included a three-gene locus (ncrABC) that conferred resistance to the bactericidal metal, nickel. ncrABC expression alone was sufficient to confer nickel resistance to E. coli K-12. Due to its selection for the colistin-resistant subpopulation harboring the duplicated 9-kb region that includes ncrABC, colistin treatment led to enhanced nickel resistance. Taken together, these data suggest that the use of antibiotics may inadvertently promote enhanced resistance to antimicrobial metals, with potentially profound implications for bacterial colonization and transmission in the health care environment.IMPORTANCE To inhibit bacterial transmission and infection, health care facilities use bactericidal metal coatings to prevent colonization of surfaces and implanted devices. In these environments, antibiotics are commonly used, but their effect on metal resistance is unclear. The data described here reveal that exposure of a human isolate of Enterobacter cloacae to a last-line antibiotic, colistin, resulted in a DNA amplification that does not confer antibiotic resistance but instead facilitates resistance to the toxic metal nickel. This highlights a novel aspect of antibiotic and metal interplay. Concerningly, these data suggest the use of antibiotics could in some cases promote bacterial survival and colonization in the health care environment and ultimately increase transmission and infection of patients.

Colistin Heteroresistance Is Largely Undetected among Carbapenem-Resistant Enterobacterales in the United States.

Heteroresistance is a form of antibiotic resistance where a bacterial strain is comprised of a minor resistant subpopulation and a majority susceptible subpopulation. We showed previously that colistin heteroresistance can mediate the failure of colistin therapy in an in vivo infection model, even for isolates designated susceptible by clinical diagnostics. We sought to characterize the extent of colistin heteroresistance among the highly drug-resistant carbapenem-resistant Enterobacterales (CRE). We screened 408 isolates for colistin heteroresistance. These isolates were collected between 2012 and 2015 in eight U.S. states as part of active surveillance for CRE. Colistin heteroresistance was detected in 10.1% (41/408) of isolates, and it was more common than conventional homogenous resistance (7.1%, 29/408). Most (93.2%, 38/41) of these heteroresistant isolates were classified as colistin susceptible by standard clinical diagnostic testing. The frequency of colistin heteroresistance was greatest in 2015, the last year of the study. This was especially true among Enterobacter isolates, of which specific species had the highest rates of heteroresistance. Among Klebsiella pneumoniae isolates, which were the majority of isolates tested, there was a closely related cluster of colistin-heteroresistant ST-258 isolates found mostly in Georgia. However, cladistic analysis revealed that, overall, there was significant diversity in the genetic backgrounds of heteroresistant K. pneumoniae isolates. These findings suggest that due to being largely undetected in the clinic, colistin heteroresistance among CRE is underappreciated in the United States.IMPORTANCE Heteroresistance is an underappreciated phenomenon that may be the cause of some unexplained antibiotic treatment failures. Misclassification of heteroresistant isolates as susceptible may lead to inappropriate therapy. Heteroresistance to colistin was more common than conventional resistance and was overwhelmingly misclassified as susceptibility by clinical diagnostic testing. Higher proportions of colistin heteroresistance observed in certain Enterobacter species and clustering among heteroresistant Klebsiella pneumoniae strains may inform colistin treatment recommendations. Overall, the rate of colistin nonsusceptibility was more than double the level detected by clinical diagnostics, suggesting that the prevalence of colistin nonsusceptibility among CRE may be higher than currently appreciated in the United States.

De novo assembly, annotation, and comparative analysis of 26 diverse maize genomes.

We report de novo genome assemblies, transcriptomes, annotations, and methylomes for the 26 inbreds that serve as the founders for the maize nested association mapping population. The number of pan-genes in these diverse genomes exceeds 103,000, with approximately a third found across all genotypes. The results demonstrate that the ancient tetraploid character of maize continues to degrade by fractionation to the present day. Excellent contiguity over repeat arrays and complete annotation of centromeres revealed additional variation in major cytological landmarks. We show that combining structural variation with single-nucleotide polymorphisms can improve the power of quantitative mapping studies. We also document variation at the level of DNA methylation and demonstrate that unmethylated regions are enriched for cis-regulatory elements that contribute to phenotypic variation.