ORFograph: search for novel insecticidal protein genes in genomic and metagenomic assembly graphs

BACKGROUND

Since the prolonged use of insecticidal proteins has led to toxin resistance, it is important to search for novel insecticidal protein genes (IPGs) that are effective in controlling resistant insect populations. IPGs are usually encoded in the genomes of entomopathogenic bacteria, especially in large plasmids in strains of the ubiquitous soil bacteria, Bacillus thuringiensis (Bt). Since there are often multiple similar IPGs encoded by such plasmids, their assemblies are typically fragmented and many IPGs are scattered through multiple contigs. As a result, existing gene prediction tools (that analyze individual contigs) typically predict partial rather than complete IPGs, making it difficult to conduct downstream IPG engineering efforts in agricultural genomics.

METHODS

Although it is difficult to assemble IPGs in a single contig, the structure of the genome assembly graph often provides clues on how to combine multiple contigs into segments encoding a single IPG.

RESULTS

We describe ORFograph, a pipeline for predicting IPGs in assembly graphs, benchmark it on (meta)genomic datasets, and discover nearly a hundred novel IPGs. This work shows that graph-aware gene prediction tools enable the discovery of greater diversity of IPGs from (meta)genomes.

CONCLUSIONS

We demonstrated that analysis of the assembly graphs reveals novel candidate IPGs. ORFograph identified both already known genes "hidden" in assembly graphs and potential novel IPGs that evaded existing tools for IPG identification. As ORFograph is fast, one could imagine a pipeline that processes many (meta)genomic assembly graphs to identify even more novel IPGs for phenotypic testing than would previously be inaccessible by traditional gene-finding methods. While here we demonstrated the results of ORFograph only for IPGs, the proposed approach can be generalized to any class of genes. Video abstract.

Phylogenetic relationship and genomic characterization of Salmonella Typhimurium strains isolated from swine in Brazil

Salmonella Typhimurium has been transmitted between humans and animals. Although, Brazil has been one of the largest pork meat exporters worldwide, there are few studies that characterized epidemiologically S. Typhimurium strains from swine. The aims of this work were to study the phylogenetic relationship of S. Typhimurium genomes isolated from swine in Brazil among themselves and with other genomes isolated from several sources and countries using wgMLST and cgMLST and to perform the search of Salmonella pathogenicity islands (SPIs). In addition, for S. Typhimurium strains from swine to compare the virulence and antimicrobial resistance genes by VFDB and ResFinder, genetic content by BLAST Atlas and orthologous proteins clusters by OrthoVenn. The constructed phylogenetic trees by wgMLST and cgMLST grouped the majority (92.3% and 80.7%, respectively) of the strains isolated from swine in Brazil into the same group. All the isolates contained important SPIs (SPI-1, SPI-2, SPI-3, SPI-5 and SPI-9). A total of 100 and 31 virulence and resistance genes were detected in the S. Typhimurium strains isolated from swine, respectively. The BLAST Atlas and orthologous proteins analysis found regions of phages and differences in metabolic, regulatory and cellular processes among S. Typhimurium LT2 and S. Typhimurium isolates from swine. In conclusion, molecular typing based in the wgMLST and cgMLST suggested that the S. Typhimurium isolates from swine studied were genetically related. The pathogenic potential of the strains studied was corroborated by the presence of important SPIs and virulence genes. The high number of antimicrobial resistance genes detected is worrying and reinforced their potential risk in swine in Brazil. The comparison by BLAST Atlas suggested differences in mobile genetic elements among S. Typhimurium LT2 and S. Typhimurium isolates from swine in Brazil. The orthologous proteins analysis revealed unique genes related to important cellular processes in the strains from swine.

Characterization of two newly isolated bacteriophages PW2 and PW4 and derived endolysins with lysis activity against Bacillus cereus group strains

This study characterized two novel Siphoviridae phages, PW2 and PW4, that can infect 52% and 44% of the tested Bacillus cereus group isolates and display relatively high activity against four cereulide-producing isolates belonging to B. weihenstephanensis and B. paranthracis. The genome sequences of PW2 and PW4 are similar to six known phages infecting B. cereus group isolates, which can be classified into two conserved groups, with the PW2 genome harboring conserved coding sequences (CDSs) from both groups. Two phage-derived endolysins, LysPW2 and LysPW4, which are predicted to encode N-acetylmuramoyl-L-alanine amidase, and their enzymatically active domains (EADs), LysPW2-EAD and LysPW4-EAD, were heterologously expressed. Both LysPW2 and LysPW4, especially the former, show a much wider host range than the phages, albeit still limited to the B. cereus group for the tested bacteria. The optimal temperature and pH for LysPW2 ability is 37 °C and pH 8.0 and for LysPW4 is 50 °C and pH 9.0. Neither LysPW2-EAD nor LysPW4-EAD show any lytic activity against vegetative cells of the tested B. cereus group isolates but can inhibit germination in 66.3% and 65.7% of spores, respectively. In addition, both LysPW2-EAD and LysPW4-EAD exhibit spore-binding capabilities.

Characterization of the complete mitochondrial genome of Ectropis grisescens (Lepidoptera, Geometridae)

Ectropis grisescens (Lepidoptera, Geometridae) is one of the main leaf-eating pests in tea plantations in China. In this study, the complete mitochondrial genome of this species was sequenced and assembled. The total length of the mitochondrial genome of E. grisescens was 15,794 bp (GenBank accession No. MW337302). The base composition was 41.26% for A, 39.49% for T, 7.92% for G, and 11.33% for C. The circular mitogenome contained 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Phylogenetic analysis performed using 13 protein-coding genes of 15 species of Geometridae and an out-group Pieris melete (Lepidoptera, Pieridae) showed that E. grisescens is closely related to species of E. obliqua, and this is consistent with the morphological identification.

High-throughput sequencing yields a complete mitochondrial genome of the Cryptotympana atrata (Insecta: Hemiptera: Cicadidae)

Cryptotympana atrata is a common insect pest found in forest ecosystem throughout East and South Asia. In this study, the complete mitochondrial genome of one individual was determined using high-throughput sequencing. The mitogenome is 15,338 bp in length with an A + T content of 77.9%, and contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and one control region (CR or D-loop). The gene arrangement and composition is similar to other published mitogenomes of Cicadidae. The concatenated PCGs were used to conduct Bayesian phylogenetic analyses together with several related Cicadidae with mitogenome data in GenBank. Phylogenetic analysis shows that two species (C. atrata and C. facialis) and Auritibicen bihamatus were herein corroborated to be the tribe of Cryptotympanini. Our results show the location of genus Cryptotympana in Cicadinae and the location of the subfamily in Cicadidae, and provide data for further study of phylogeny in Hemiptera.

Genomic Characterization of mcr-1.1-Producing Escherichia coli Recovered From Human Infections in São Paulo, Brazil

Polymyxins are one of most important antibiotics available for multidrug-resistant Gram-negative infections. Diverse chromosomal resistance mechanisms have been described, but the polymyxin resistance phenotype is not yet completely understood. The objective of this study was to characterize colistin resistant mcr-1-producing strains isolated from human infections over one year in a hospital setting (Hospital das Clínicas, São Paulo, Brazil). We isolated 490 colistin-resistant Gram-negative rods, of which eight were mcr-1.1-positive Escherichia coli, the only species with this result, indicating a low incidence of the mcr-1 production mechanism among colistin-resistant isolates. All mcr-1.1 positive isolates showed similarly low MICs for colistin and were susceptible to most antibiotics tested. The isolates showed diversity of MLST classification. The eight mcr-1.1-positive E. coli genomes were sequenced. In seven of eight isolates the mcr-1.1 gene is located in a contig that is presumed to be a part of an IncX4 plasmid; in one isolate, it is located in a contig that is presumed to be part of an IncHI2A plasmid. Three different genomic contexts for mcr-1.1 were observed, including a genomic cassette mcr-1.1-pap2 disrupting a DUF2806 domain-containing gene in six isolates. In addition, an IS1-family transposase was found inserted next to the mcr-1.1 cassette in one isolate. An mcr-1.1-pap2 genomic cassette not disrupting any gene was identified in another isolate. Our results suggest that plasmid dissemination of hospital-resident strains took place during the study period and highlight the need for continued genomic surveillance.

The complete mitochondrial genome of Pentatoma rufipes (Hemiptera, Pentatomidae) and its phylogenetic implications

Pentatoma rufipes (Linnaeus, 1758) is an important agroforestry pest widely distributed in the Palaearctic region. In this study, we sequence and annotate the complete mitochondrial genome of P. rufipes and reconstruct the phylogenetic trees for Pentatomoidea using existing data for eight families published in the National Center for Biotechnology Information database. The mitogenome of P. rufipes is 15,887-bp-long, comprising 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a control region, with an A+T content of 77.7%. The genome structure, gene order, nucleotide composition, and codon usage of the mitogenome of P. rufipes were consistent with those of typical Hemiptera insects. Among the protein-coding genes of Pentatomoidea, the evolutionary rate of ATP8 was the fastest, and COX1 was found to be the most conservative gene in the superfamily. Substitution saturation assessment indicated that neither transition nor transversion substitutions were saturated in the analyzed datasets. Phylogenetic analysis using the Bayesian inference method showed that P. rufipes belonged to Pentatomidae. The node support values based on the dataset concatenated from protein-coding and RNA genes were the highest. Our results enrich the mitochondrial genome database of Pentatomoidea and provide a reference for further studies of phylogenetic systematics.

Genome profiling of fluoroquinolone-resistant uropathogenic Escherichia coli isolates from Brazil

Urinary tract infections (UTIs) are a major public health concern in both community and hospital settings worldwide. Uropathogenic Escherichia coli (UPEC) is the main causative agent of UTI and increasingly associated with antibiotic resistance. Herein, we report the draft genome sequence of 9 fluoroquinolone-resistant UPEC isolates from Brazil and examine selected major phenotypic features, such as antimicrobial resistance profile, phylogroup, serotype, sequence type (ST), virulence genes, and resistance marks. Besides the quinolone resistance, beta-lactams, ESBL production, aminoglycosides, and tetracycline resistance were observed. High prevalence of 20 virulence genes was detected in all isolates, such as those encoding type 1 fimbriae, acid tolerance system, and hemolysin E, particularly within E. coli B2 phylogroup, as ST131 and ST1193 strains, among other genomic analyses as genomic islands, resistance plasmids, and integron identification.

Effect of treatment of pneumonia and otitis media with tildipirosin or florfenicol + flunixin meglumine on health and upper respiratory tract microbiota of preweaned Holstein dairy heifers

The objective of this randomized clinical study was to compare the effect of 2 antimicrobial interventions, tildipirosin or florfenicol + flunixin meglumine, used for treatment of pneumonia and extralabel treatment for otitis on health parameters and upper respiratory tract (URT) microbiota of preweaned Holstein calves. Housed preweaned Holstein heifers diagnosed with either otitis or pneumonia were assigned into 1 of 2 treatment groups, receiving a single subcutaneous injection of either 4 mg/kg of tildipirosin (TLD; n = 444) or 40 mg/kg of florfenicol combined with 2.2 mg/kg of a nonsteroidal anti-inflammatory, flunixin meglumine (FLF; n = 442). Calves were enrolled and treated on the day of diagnosis of the first case of pneumonia or otitis. If a calf had a recurrent case, the opposite drug was administered, respecting an interval of 5 d between drug injections. Blood samples for leukocyte counts were collected at 0, 2, 4, and 6 d after treatment, and rectal temperature was measured daily during the 5 d after treatment. Ear scores were observed from calves with otitis. Additionally, swabs of the URT were collected from a subset of 20 calves in each treatment group at d 0, 3, 6, 9, and 11 following enrollment for analysis of URT microbiota through next-generation sequencing of the 16S rRNA gene and quantitative PCR. Swabs were also collected from a comparative group of 20 healthy calves that did not receive any drug. No differences were observed between groups for recurrence risk of either pneumonia (TLD = 32.4%; FLF = 29.7%) or otitis (TLD = 72.7%; FLF = 73.6%). Similarly, no differences were observed for the total number of treatments for pneumonia (TLD = 1.45; FLF = 1.42) or otitis (TLD = 2.96; FLF = 3.07). On the other hand, both drugs reduced rectal temperature, ear scores, and leukocyte counts, with FLF calves having a greater reduction in rectal temperature within 4 d after treatment. Both TLD and FLF reduced the total bacterial load when compared with healthy untreated calves, but no differences were observed between treatment groups. Furthermore, compared with the untreated group, treated calves had lower mean relative abundances (MRA) of the genera Mannheimia, Moraxella, and Pasteurella within 11, 9, and 3 d after treatment, respectively; however, no significant differences were observed between TLD and FLF. On the other hand, MRA of Mycoplasma was not decreased by both treatments compared to untreated animals, and a higher MRA was observed in the TLD group during 11 d after treatment in comparison to FLF and untreated calves. Based on this data, we concluded that both drugs used in the study were effective in reducing rectal temperature, ear scores, leukocyte counts, and MRA of the genera Mannheimia, Pasteurella, and Moraxella in the URT, and calves treated with FLF had a greater reduction in rectal temperature.

Isolation and Genomic Analysis of the Phage vB_PaeP_fHoPae04 Infecting Pseudomonas aeruginosa

Here, we report the genomic sequence of Pseudomonas aeruginosa phage vB_PaeP_fHoPae04, isolated from hospital wastewater in Helsinki, Finland. The phage genome is 45,491 bp long, has a G+C content of 52.2%, and contains 70 protein-coding genes and 3 tRNA genes.

Characterization of Three Complete Mitogenomes of Flatidae (Hemiptera: Fulgoroidea) and Compositional Heterogeneity Analysis in the Planthoppers' Mitochondrial Phylogenomics

Although sequences of mitogenomes have been widely used for investigating phylogenetic relationship, population genetics, and biogeography in many members of Fulgoroidea, only one complete mitogenome of a member of Flatidae has been sequenced. Here, the complete mitogenomes of Cerynia lineola, Cromna sinensis, and Zecheuna tonkinensis are sequenced. The gene arrangements of the three new mitogenomes are consistent with ancestral insect mitogenomes. The strategy of using mitogenomes in phylogenetics remains in dispute due to the heterogeneity in base composition and the possible variation in evolutionary rates. In this study, we found compositional heterogeneity and variable evolutionary rates among planthopper mitogenomes. Phylogenetic analysis based on site-homogeneous models showed that the families (Delphacidae and Derbidae) with high values of Ka/Ks and A + T content tended to fall together at a basal position on the trees. Using a site-heterogeneous mixture CAT + GTR model implemented in PhyloBayes yielded almost the same topology. Our results recovered the monophyly of Fulgoroidea. In this study, we apply the heterogeneous mixture model to the planthoppers’ phylogenetic analysis for the first time. Our study is based on a large sample and provides a methodological reference for future phylogenetic studies of Fulgoroidea.

Extensive Comparative Genomic Analysis of Enterococcus faecalis and Enterococcus faecium Reveals a Direct Association between the Absence of CRISPR-Cas Systems, the Presence of Anti-Endonuclease (ardA) and the Acquisition of Vancomycin Resistance in E. faecium

Here, we performed a comparative genomic analysis of all available genomes of E. faecalis (n = 1591) and E. faecium (n = 1981) and investigated the association between the presence or absence of CRISPR-Cas systems, endonuclease/anti-endonuclease systems and the acquisition of antimicrobial resistance, especially vancomycin resistance genes. Most of the analysed Enterococci were isolated from humans and less than 14% of them were from foods and animals. We analysed and detected CRISPR-Cas systems in 75.36% of E. faecalis genomes and only 4.89% of E. faecium genomes with a significant difference (p-value < 10-5). We found a negative correlation between the number of CRISPR-Cas systems and genome size (r = -0.397, p-value < 10-5) and a positive correlation between the genome %GC content and the number of CRISPR-Cas systems (r = 0.215, p-value < 10-5). Our findings showed that the presence of the anti-endonuclease ardA gene may explain the decrease in the number of CRISPR-Cas systems in E. faecium, known to deactivate the endonucleases' protective activities and enable the E. faecium genome to be versatile in acquiring mobile genetic elements, including carriers of antimicrobial resistance genes, especially vanB. Most importantly, we observed that there was a direct association between the absence of CRISPR-Cas, the presence of the anti-CRISPR ardA gene and the acquisition of vancomycin resistance genes.

A novel temperate phage, vB_PstS-pAN, induced from the naphthalene-degrading bacterium Pseudomonas stutzeri AN10

A novel temperate phage named vB_PstS-pAN was induced by mitomycin C treatment from the naphthalene-degrading bacterium Pseudomonas stutzeri AN10. The phage particles have icosahedral heads and long non-contractile tails, and vB_PstS-pAN can therefore be morphologically classified as a member of the family Siphoviridae. The whole genome of vB_PstS-pAN is 39,466 bp in length, with an 11-nt 3' overhang cohesive end. There are 53 genes in the vB_PstS-pAN genome, including genes responsible for phage integration, replication, morphogenesis, and bacterial lysis. The vB_PstS-pAN genome has low similarity to other phage genomes in the GenBank database, suggesting that vB_PstS-pAN is a novel member of the family Siphoviridae.

Evolution of virulence in a novel family of transmissible mega-plasmids

Some Serratia entomophila isolates have been successfully exploited in biopesticides due to their ability to cause amber disease in larvae of the Aotearoa (New Zealand) endemic pasture pest, Costelytra giveni. Anti-feeding prophage and ABC toxin complex virulence determinants are encoded by a 153-kb single-copy conjugative plasmid (pADAP; amber disease-associated plasmid). Despite growing understanding of the S. entomophila pADAP model plasmid, little is known about the wider plasmid family. Here, we sequence and analyse mega-plasmids from 50 Serratia isolates that induce variable disease phenotypes in the C. giveni insect host. Mega-plasmids are highly conserved within S. entomophila, but show considerable divergence in Serratia proteamaculans with other variants in S. liquefaciens and S. marcescens, likely reflecting niche adaption. In this study to reconstruct ancestral relationships for a complex mega-plasmid system, strong co-evolution between Serratia species and their plasmids were found. We identify 12 distinct mega-plasmid genotypes, all sharing a conserved gene backbone, but encoding highly variable accessory regions including virulence factors, secondary metabolite biosynthesis, Nitrogen fixation genes and toxin-antitoxin systems. We show that the variable pathogenicity of Serratia isolates is largely caused by presence/absence of virulence clusters on the mega-plasmids, but notably, is augmented by external chromosomally encoded factors.

Characterization of the complete chloroplast genome of basswood (Tilia mongolica maxim.)

Herein, we report the complete chloroplast genome of Tilia mongolica Maxim. from Tiliaceae. The chloroplast genome of T. mongolica is 162,804 bp, with a large single copy region of 91,255 bp, small single copy region of 20,355 bp, and two inverted-repeat regions of 25,597 bp. The chloroplast genome contains 130 genes, including 85 protein-coding, 8 rRNA, and 37 tRNA. The total GC content is 36.46%. The phylogenetic analysis of T. mongolica showed a relatively close relationship with Tilia taishanensis.

Integrated Genomic and Functional Characterization of the Anti-diabetic Potential of Arthrobacter sp. SW1

For decades, bacterial natural products have served as valuable resources for developing novel drugs to treat several human diseases. Recent advancements in the integrative approach of using genomic and functional tools have proved beneficial in obtaining a comprehensive understanding of these biomolecules. This study presents an in-depth characterization of the anti-diabetic activity exhibited by a bacterial isolate SW1, isolated from an effluent treatment plant. As a primary screening, we assessed the isolate for its potential to inhibit alpha-amylase and alpha-glucosidase enzymes. Upon confirmation, we further utilized LC-MS, ESI-MS/MS, and NMR spectroscopy to identify and characterize the biomolecule. These efforts were coupled with the genomic assessment of the biosynthetic gene cluster involved in the anti-diabetic compound production. Our investigation discovered that the isolate SW1 inhibited both α-amylase and α-glucosidase activity. The chemical analysis suggested the production of acarbose, an anti-diabetic biomolecule, which was further confirmed by the presence of biosynthetic gene cluster "acb" in the genome. Our in-depth chemical characterization and genome mining approach revealed the potential of bacteria from an unconventional niche, an effluent treatment plant. To the best of our knowledge, it is one of the first few reports of acarbose production from the genus Arthrobacter.

Cell division in the archaeon Haloferax volcanii relies on two FtsZ proteins with distinct functions in division ring assembly and constriction

In bacteria, the tubulin homologue FtsZ assembles a cytokinetic ring, termed the Z ring, and plays a key role in the machinery that constricts to divide the cells. Many archaea encode two FtsZ proteins from distinct families, FtsZ1 and FtsZ2, with previously unclear functions. Here, we show that Haloferax volcanii cannot divide properly without either or both FtsZ proteins, but DNA replication continues and cells proliferate in alternative ways, such as blebbing and fragmentation, via remarkable envelope plasticity. FtsZ1 and FtsZ2 colocalize to form the dynamic division ring. However, FtsZ1 can assemble rings independent of FtsZ2, and stabilizes FtsZ2 in the ring, whereas FtsZ2 functions primarily in the constriction mechanism. FtsZ1 also influenced cell shape, suggesting it forms a hub-like platform at midcell for the assembly of shape-related systems too. Both FtsZ1 and FtsZ2 are widespread in archaea with a single S-layer envelope, but archaea with a pseudomurein wall and division septum only have FtsZ1. FtsZ1 is therefore likely to provide a fundamental recruitment role in diverse archaea, and FtsZ2 is required for constriction of a flexible S-layer envelope, where an internal constriction force might dominate the division mechanism, in contrast with the single-FtsZ bacteria and archaea that divide primarily by wall ingrowth.

Screening of Bacteriophage Encoded Toxic Proteins with a Next Generation Sequencing-Based Assay

Bacteriophage vB_EcoM_fHy-Eco03 (fHy-Eco03 for short) was isolated from a sewage sample based on its ability to infect an Escherichia coli clinical blood culture isolate. Altogether, 32 genes encoding hypothetical proteins of unknown function (HPUFs) were identified from the genomic sequence of fHy-Eco03. The HPUFs were screened for toxic properties (toxHPUFs) with a novel, Next Generation Sequencing (NGS)-based approach. This approach identifies toxHPUF-encoding genes through comparison of gene-specific read coverages in DNA from pooled ligation mixtures before electroporation and pooled transformants after electroporation. The performance and reliability of the NGS screening assay was compared with a plating efficiency-based method, and both methods identified the fHy-Eco03 gene g05 product as toxic. While the outcomes of the two screenings were highly similar, the NGS screening assay outperformed the plating efficiency assay in both reliability and efficiency. The NGS screening assay can be used as a high throughput method in the search for new phage-inspired antimicrobial molecules.

Differential regulation of undecylprodigiosin biosynthesis in the yeast-scavenging Streptomyces strain MBK6

Streptomyces are efficient chemists with a capacity to generate diverse and potent chemical scaffolds. The secondary metabolism of these soil-dwelling prokaryotes is stimulated upon interaction with other microbes in their complex ecosystem. We observed such an interaction when a Streptomyces isolate was cultivated in a media supplemented with dead yeast cells. Whole-genome analysis revealed that Streptomyces sp. MBK6 harbors the red cluster that is cryptic under normal environmental conditions. An interactive culture of MBK6 with dead yeast triggered the production of the red pigments metacycloprodigiosin and undecylprodigiosin. Streptomyces sp. MBK6 scavenges dead-yeast cells and preferentially grows in aggregates of sequestered yeasts within its mycelial network. We identified that the activation depends on the cluster-situated regulator, mbkZ, which may act as a cross-regulator. Cloning of this master regulator mbkZ in S. coelicolor with a constitutive promoter and promoter-deprived conditions generated different production levels of the red pigments. These surprising results were further validated by DNA–protein binding assays. The presence of the red cluster in Streptomyces sp. MBK6 provides a vivid example of horizontal gene transfer of an entire metabolic pathway followed by differential adaptation to a new environment through mutations in the receiver domain of the key regulatory protein MbkZ.

Association between the blaCTX-M-14-harboring Escherichia coli Isolated from Weasels and Domestic Animals Reared on a University Campus

Antimicrobial-resistant (AMR) bacteria affect human and animal health worldwide. Here, CTX-M-14-producing Escherichia coli isolates were isolated from Siberian weasels (Mustela sibirica) that were captured on a veterinary campus. To clarify the source of bacteria in the weasels, we examined the domestic animals reared in seven facilities on the campus. Extended-spectrum β-lactamase (ESBL)-producing E. coli were isolated on deoxycholate hydrogen sulfide lactose agar, containing cephalexin (50 μg/mL) or cefotaxime (2 μg/mL), and were characterized with antimicrobial susceptibility testing, pulsed-field gel electrophoresis (PFGE), replicon typing, and β-lactamase typing analyses. Next-generation sequencing of the ESBL-encoding plasmids was also performed. CTX-M-14 producers isolated from both domestic animals and weasels were classified into six clusters with seven PFGE profiles. The PFGE and antimicrobial resistance profiles were characterized by the animal facility. All CTX-M-14 plasmids belonged to the IncI1 type with a similar size (98.9-99.3 kb), except for one plasmid that was 105.5 kb in length. The unweighted pair group method with arithmetic mean (UPGMA) revealed that the CTX-M-14 plasmid in the weasel isolates might have the same origin as the CTX-M-14 plasmid in the domestic animals. Our findings shed further light on the association of antimicrobial resistance between wild and domestic animals.

Complete nucleotide sequence analysis and identification of 7-cyano-7-deazaguanine (PreQ0) biosynthesis-related genes in the novel Bacillus subtilis-infecting Siphoviridae family phage BSP7

In this study, bacteriophage BSP7, a novel Bacillus subtilis-infecting member of the family Siphoviridae, was isolated from a Korean soybean-based fermented food, Deonjang, using B. subtilis ATCC 21336 as a host. The genome is 55,455 bp long with 39.92% G+C content. A total of 70 ORFs with no tRNA were detected in the genome. A distinct feature of the BSP7 genome among B. subtilis-infecting Siphoviridae family phages is the presence of putative ORFs related to biosynthesis of 7-cyano-7-deazaguanine (PreQ₀), a precursor of queuosine and archaeosine biosynthesis. Bioinformatic analysis revealed that the genome of BSP7 does not exhibit any significant similarities to other phages with sequences in the NCBI database. A comparative genomic analysis also confirmed the uniqueness of BSP7 within the family Siphoviridae.

Paenibacillus glycinis sp. nov., an Endophytic Bacterium Isolated from the Nodules of Soybean (Glycine max (L.) Merr)

A Gram-negative-staining, endospore-forming, rod-shaped bacterium, designated T1T, was isolated from root nodules of soybean (Glycine max (L.) Merr) in Heilongjiang Province of China. The isolate was identified as a member of the genus Paenibacillus based on phenotypic and phylogenetic characteristics. The 16S rRNA sequence was closely related to that of Paenibacillus sacheonensis SY01T with a similarity of 98.4%. Average nucleotide identity and in silico DNA-DNA hybridization values between strain T1T and P. sacheonensis DSM 23054 T were 81.4% and 25.4%, respectively. The DNA G + C content of strain T1T was 58.2 mol%. meso-Diaminopimelic acid was detected in the cell-wall peptidoglycan. The major cellular fatty acids were anteiso-C15:0, iso-C16:0 and iso-C15:0. The predominant respiratory quinone was MK-7. The polar lipids were diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, five unidentified phospholipids, four unidentified aminophospholipids, an unidentified glycolipid and an unidentified lipid. Based on these results, T1T is considered to represent a novel species of the genus Paenibacillus, for which the name Paenibacillus glycinis sp. nov. is proposed. The type strain is T1T (= CGMCC 1.18563 = KCTC43227).

Streptomyces typhae sp. nov., a novel endophytic actinomycete with antifungal activity isolated the root of cattail (Typha angustifolia L.)

A novel endophytic actinomycete with antagonistic activity against various phytopathogenic fungi, designated strain p1417^T, was isolated from the root of cattail (Typha angustifolia L.) collected from Yunnan Province, Southwest China. A polyphasic taxonomic study was carried out to establish the taxonomic status of this strain. Strain p1417^T was found to have morphological and chemotaxonomic characteristics typical of the genus Streptomyces. The diamino acid present in the cell wall was LL-diaminopimelic acid. Xylose and arabinose occurred in whole cell hydrolysates. The menaquinones were identified as MK-9(H₈), MK-9(H₆) and MK-9(H₄). The polar lipid profile was found to contain diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylinositol and phosphatidylinositol mannoside. The major fatty acids were found to be iso-C_16:0, anteiso-C_15:0, iso-C_15:0 and C_16:0. The genomic DNA G + C content of strain p1417^T based on the genome sequence was 72.0 mol%. Based on 16 S rRNA gene, five housekeeping genes and whole genome sequences analysis, strain p1417^T was most closely related to Streptomyces flavofungini JCM 4753^T (99.4% 16 S rRNA gene sequence similarity), Streptomyces alboflavus JCM 4615^T (98.8%) and Streptomyces aureoverticillatus JCM 4347^T (98.2%). However, the average nucleotide identity values, the digital DNA-DNA hybridization values and the multilocus sequence analysis evolutionary distances between this strain and its closely related strains showed that it belonged to one distinct species. In addition, these results were also supported by differences in the phenotypic and chemotaxonomic characteristics between strain p1417^T and three closely related type strains. Therefore, it is concluded that strain p1417^T represents a novel species of the genus of Streptomyces, for which the name Streptomyces typhae sp. nov. is proposed. The type strain is p1417^T (= CCTCC AA 2019091^T = DSM 110636^T).

Evolutionary genomic and bacteria GWAS analysis of Mycobacterium avium subsp. paratuberculosis and dairy cattle Johne's disease phenotypes

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle. The Regional Dairy Quality Management Alliance (RDQMA) project is a multistate research program involving MAP isolates taken from three intensively studied commercial dairy farms in the northeastern United States, which emphasized longitudinal data collection of both MAP isolates and animal health in three regional dairy herds for a period of about 7 years. This paper reports the results of a pan-GWAS analysis involving 318 MAP isolates and dairy cow Johne's disease phenotypes, taken from these three farms. Based on our highly curated accessory gene count the pan-GWAS analysis identified several MAP genes associated with bovine Johne's disease phenotypes scored from these three farms, with some of the genes having functions suggestive of possible cause/effect relationships to these phenotypes. This paper reports a pan-genomic comparative analysis between MAP and Mycobacterium tuberculosis, assessing functional Gene Ontology category enrichments between these taxa. Finally, we also provide a population genomic perspective on the effectiveness of herd isolation, involving closed dairy farms, in preventing MAP inter-farm cross infection on a micro-geographic scale.IMPORTANCE Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle, and enormous economic consequences for the dairy industry. Understanding which genes in this bacterium are correlated with key disease phenotypes can lead to functional experiments targeting these genes and ultimately lead to improved control strategies. This study represents a rare example of a prolonged longitudinal study of dairy cattle where the disease was measured and the bacteria were isolated from the same cows. The genome sequences of over 300 MAP isolates were analyzed for genes that were correlated with a wide range of Johne's disease phenotypes. A number of genes were identified that were significantly associated with several aspects of the disease and suggestive of further experimental follow-up.

Draft Genome Sequences of Dysgonomonas sp. Strains GY75 and GY617, Isolated from the Hindgut of Reticulitermes flavipes

Dysgonomonas species are facultative heterotrophs capable of growth on lignocellulose-derived polysaccharides. Dysgonomonas species harbor myriad genes involved in glycan modification and are well suited to the lignocellulose-rich conditions within the termite hindgut. Here, we report draft genome sequences for Dysgonomonas sp. strains GY75 and GY617, isolated from the hindgut of Reticulitermes flavipes.

Dysgonomonas species are facultative heterotrophs capable of growth on lignocellulose-derived polysaccharides. Dysgonomonas species harbor myriad genes involved in glycan modification and are well suited to the lignocellulose-rich conditions within the termite hindgut. Here, we report draft genome sequences for Dysgonomonas sp. strains GY75 and GY617, isolated from the hindgut of Reticulitermes flavipes.

Draft Genome Sequence of Enterococcus faecalis AS003, a Strain Possessing All Three Type II-a CRISPR Loci

Enterococcus faecalis is a clinically significant member of the human microbiome. Three CRISPR-Cas loci are located in conserved locations. Previous studies provide evidence that E. faecalis strains with functional CRISPR-Cas genes are negatively correlated with antibiotic resistance. Here, we report the genome sequence of an unusual strain possessing all three CRISPR-Cas loci.

Whole-Genome Sequencing and Comparative Genomics of Three Helicobacter pylori Strains Isolated from the Stomach of a Patient with Adenocarcinoma

Helicobacter pylori is a common pathogen associated with several severe digestive diseases. Although multiple virulence factors have been described, it is still unclear the role of virulence factors on H. pylori pathogenesis and disease progression. Whole genome sequencing could help to find genetic markers of virulence strains. In this work, we analyzed three complete genomes from isolates obtained at the same point in time from a stomach of a patient with adenocarcinoma, using multiple available bioinformatics tools. The genome analysis of the strains B508A-S1, B508A-T2A and B508A-T4 revealed that they were cagA, babA and sabB/hopO negative. The differences among the three genomes were mainly related to outer membrane proteins, methylases, restriction modification systems and flagellar biosynthesis proteins. The strain B508A-T2A was the only one presenting the genotype vacA s1, and had the most distinct genome as it exhibited fewer shared genes, higher number of unique genes, and more polymorphisms were found in this genome. With all the accumulated information, no significant differences were found among the isolates regarding virulence and origin of the isolates. Nevertheless, some B508A-T2A genome characteristics could be linked to the pathogenicity of H. pylori.

The complete chloroplast genome of Chimonobambusa purpurea (Bambusatae)

Chimonobambusa purpurea is one of the important bamboo species in southwest of China. We studied the complete chloroplast (cp) genome of C. purpurea in this study. The cp genome of C. purpurea (GenBank accession: MW030500) was 139,574 bp in length, including a large single-copy (LSC) region of 83,171 bp, a small single-copy (SSC) region of 12,811 bp, and a pair of inverted repeated (IR) regions of 21,796 bp. And the genome contained 133 genes, including 86 protein-coding genes, 39 tRNA genes, and 8 rRNA genes. Based on 30 cp genomes, we used the phylogenetic analysis to build phylogenetic tree, indicating that C. purpurea is closely related to C. tumidissinoda.

The complete mitochondrial genome of the spiny red gurnard Chelidonichthys spinosus McClelland, 1844 (Scorpaeniformes: Triglidae)

The spiny red gurnard, Chelidonichthys spinosus is a common marine economic fish species along the coast of China. In the present study, the complete mitochondrial genome of spiny red gurnard collected from the Yangtze Estuary was determined by next-generation sequencing (NGS). The mitogenome is a circular nucleotide 16,466 bp in length and has the typical vertebrate genome structure of 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and two control regions (Origin of Heavy chain and Origin of Light chain). The termination associated sequence (TAS), the central conserved sequence block (central CSB), and the conserved sequence block (CSB) were detected in the control region. Phylogenetic analysis of C. spinosus places in a fully supported clade with C. kumu in a sister position to Eutrigla gurnardus. The complete mitochondrial genome sequence of the spiny red gurnard provides baseline genetic information for future studies on the molecular systematics and phylogeny of bony fishes in the Triglidae.

The complete chloroplast genome of Schizonepeta tenuifolia (Benth.) Briq., a traditional Chinese herb

Schizonepeta tenuifolia (Benth.) Briq. is a traditional Chinese medicinal herb. The complete chloroplast genome sequence of S. tenuifolia was obtained by high-throughput sequencing platform. The chloroplast genome of S. tenuifolia is a circular form of 151,254 bp in length, with an average GC content of 37.85%. The genome contains a set of 132 genes, including 87 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis based on complete chloroplast genome sequences indicates that S. tenuifolia has a close relationship with Dracocephalum palmatum. This study provides a molecular basis for the classification of S. tenuifolia.

Elevation of Lactococcus lactis subsp. cremoris to the species level as Lactococcus cremoris sp. nov. and transfer of Lactococcus lactis subsp. tructae to Lactococcus cremoris as Lactococcus cremoris subsp. tructae comb. nov

Currently, Lactococcus lactis contains four subspecies: L. lactis subsp. lactis, L. lactis subsp. hordniae, L. lactis subsp. cremoris and L. lactis subsp. tructae. In the study of Pérez et al., these four subspecies could be clearly divided into two groups based on recA sequence analysis: L. lactis subsp. lactis and L. lactis subsp. hordniae; L. lactis subsp. cremoris and L. lactis subsp. tructae. The two groups had a relatively low DNA-DNA hybridization value (about 60 %). In the present study, the taxonomic position of L. lactis subsp. cremoris and L. lactis subsp. tructae was re-examined based on sequence analyses of 16S rRNA, rpoB, recA and pheS genes, average nucleotide identity (ANI) values and digital DNA-DNA hybridization (dDDH) values. The result of 16S rRNA gene sequence analysis indicated that L. lactis subsp. cremoris NCDO 607^T and L. lactis subsp. tructae L105^T were phylogenetically related to the type strains of L. lactis subsp. hordniae, L. lactis subsp. lactis, Lactococcus taiwanensis, Lactococcus kimchii, Lactococcus allomyrinae, Lactococcus protaetiae, Lactococcus hircilactis, Lactococcus fujiensis and Lactococcus nasutitermitis. The 16S rRNA gene, rpoB, recA, pheS and concatenated rpoB, recA and pheS sequence similarities, ANI values, and dDDH values between the type strains of L. lactis subsp. cremoris, L. lactis subsp. tructae and phylogenetically related species were 93.5-99.4 %, 83.3-97.6 %, 80.6-92.4 %, 79.7-92.7 %, 83.5-94.3 %, 72.4-86.9 % and 21.4-32.5 %, respectively. Lower than 95-96 % ANI values and lower than 70 % dDDH values confirmed that the type strains of L. lactis subsp. cremoris and L. lactis subsp. tructae represent a novel species in the genus Lactococcus. Because L. lactis subsp. cremoris was proposed and validated before L. lactis subsp. tructae, L. lactis subsp. cremoris is elevated to the species level as Lactococcus cremoris sp. nov. and L. lactis subsp. tructae is transferred to L. cremoris as L. cremoris subsp. tructae comb. nov. The type strain of L. cremoris sp. nov. is NCDO 607^T (=ATCC 19257^T=DSM 20069^T=JCM 16167^T=LMG 6897^T=NBRC 100676^T). The type strain of L. cremoris subsp. tructae comb. nov. is L105^T (=NBRC 110453^T=DSM 21502^T=JCM 31125^T=LMG 24662^T).

Transmission, adaptation and geographical spread of the Pseudomonas aeruginosa Liverpool epidemic strain

The Liverpool epidemic strain (LES) is an important transmissible clonal lineage of Pseudomonas aeruginosa that chronically infects the lungs of people with cystic fibrosis (CF). Previous studies have focused on the genomics of the LES in a limited number of isolates, mostly from one CF centre in the UK, and from studies highlighting identification of the LES in Canada. Here we significantly extend the current LES genome database by genome sequencing 91 isolates from multiple CF centres across the UK, and we describe the comparative genomics of this large collection of LES isolates from the UK and Canada. Phylogenetic analysis revealed that the 145 LES genomes analysed formed a distinct clonal lineage when compared with the wider P. aeruginosa population. Notably, the isolates formed two clades: one associated with isolates from Canada, and the other associated with UK isolates. Further analysis of the UK LES isolates revealed clustering by clinic geography. Where isolates clustered closely together, the association was often supported by clinical data linking isolates or patients. When compared with the earliest known isolate, LESB58 (from 1988), many UK LES isolates shared common loss-of-function mutations, such as in genes gltR and fleR. Other loss-of-function mutations identified in previous studies as common adaptations during CF chronic lung infections were also identified in multiple LES isolates. Analysis of the LES accessory genome (including genomic islands and prophages) revealed variations in the carriage of large genomic regions, with some evidence for shared genomic island/prophage complement according to clinic location. Our study reveals divergence and adaptation during the spread of the LES, within the UK and between continents.

A novel vieuvirus from multidrug-resistant Acinetobacter baumannii

Bacteriophages are considered the most abundant biological entities on earth, and they are able to modulate the populations of their bacterial hosts. Although the potential of bacteriophages has been accepted as an alternative strategy to combat multidrug-resistant pathogenic bacteria, there still exists a considerable knowledge gap regarding their genetic diversity, which hinders their use as antimicrobial agents. In this study, we undertook a genomic and phylogenetic characterization of the phage Ab11510-phi, which was isolated from a multidrug-resistant Acinetobacter baumannii strain (Ab11510). We found that Ab11510-phi has a narrow host range and belongs to a small group of transposable phages of the genus Vieuvirus that have only been reported to infect Acinetobacter bacteria. Finally, we showed that Ab11510-phi (as well as other vieuvirus phages) has a high level of mosaicism. On a broader level, we demonstrate that comparative genomics and phylogenetic analysis are necessary tools for the proper characterization of phage diversity.

Antimicrobial resistance in Enterobacteriaceae isolated from arthropods in Gifu City, Japan

The wide occurrence of antimicrobial-resistant (AMR) bacteria in various environments is of great concern. Here, we examined the prevalence and antimicrobial susceptibility of Enterobacteriaceae isolated from 88 wild arthropods, collected in Gifu city, Japan. In total, 168 isolates of Enterobacteriaceae were obtained from 61 arthropods. All isolates were susceptible to all the antimicrobial agents tested, except colistin (31 isolates) and kanamycin (one isolate). The aph(3')-Ia gene, responsible for kanamycin resistance, was detected in Klebsiella oxytoca. Although synanthropic arthropods (houseflies and cockroaches) serve as vectors for AMR Enterobacteriaceae, other wild arthropods are not crucial carriers of Enterobacteriaceae resistant to antimicrobial agents.

Proposal of Enterococcus xinjiangensis Ren et al. 2020 as a later heterotypic synonym of Enterococcus lactis Morandi et al. 2012

Enterococcus lactis was published in the International Journal of Systematic and Evolutionary Microbiology in June 2012. Enterococcus xinjiangensis was published in Current Microbiology in 2016 and validated in November 2020. In the present study, the relationship between E. lactis and E. xinjiangensis was re-evaluated. In the 16S rRNA gene phylogeny, E. xinjiangensis was closely related to E. lactis and Enterococcus faecium. The type strains of E. xinjiangensis and E. lactis shared 99.8 % 16S rRNA gene sequence similarity, 98.5 % pheS sequence similarity, 99.9 % rpoA sequence similarity, 98.7 % average nucleotide identity (ANI) value and 88.9 % digital DNA-DNA hybridization (dDDH) value, indicating that they represent the same species. Meanwhile, although the type strains of E. xinjiangensis and E. faecium also shared 99.8 % 16S rRNA gene sequence similarity, 97.1 % pheS sequence similarity and 99.6 % rpoA sequence similarity, 94.7 % ANI and 59.1 % dDDH values indicated that they represent two different species. On the basis of the results present here, we propose Enterococcus xinjiangensis Ren et al. 2020 as a later heterotypic synonym of Enterococcus lactis Morandi et al. 2012.

Identification and complete genome of lytic "Kp34likevirus" phage vB_KpnP_Bp5 and therapeutic potency in the treatment of lethal Klebsiella pneumoniae infections in mice

Klebsiella pneumoniae (K. pneumoniae) infection exist widely in the farming and medical. The treatment of K. pneumoniae infection is primarily based on antibiotics, which not only leads to a large economic burden but also the development of antibiotic resistance. Bacteriophages therapy present a promising alternative. The object of this study was identifying comprehensively a lytic lethal K. pneumoniae phage vB_KpnP_Bp5, and evaluating the phage as an anti-infective agent in an experimental K. pneumoniae infection murine model. The phage Bp5 had the following characteristics: the optimal number of infections was 0.001, the latent period was 5 min, the outbreak period was 40 min, the burst size was 24 plaque-forming unit (PFU)/cell, the phage could withstand 50 °C temperature and the optimal pH value was 4.0-10.0. According to electron microscopy and whole-genome sequence analysis, the phage should be classified as a member of order Caudovirales, family Podoviridae, subfamily Autographiviridae. Meantime, phylogenetic analysis showed high conservation of gene arrangement and gene content. We demonstrated that administration of phage Bp5 significantly reduced colony formation by K. pneumoniae and alleviated damage to lung tissue. In addition, different therapy time point was closely related to body health and the degree of tissue damage. Once treated promptly, it will greatly reduce mortality and alveolar inflammatory exudation and injury.

Genome-Wide Association Studies for the Detection of Genetic Variants Associated With Daptomycin and Ceftaroline Resistance in Staphylococcus aureus

Background

As next generation sequencing (NGS) technologies have experienced a rapid development over the last decade, the investigation of the bacterial genetic architecture reveals a high potential to dissect causal loci of antibiotic resistance phenotypes. Although genome-wide association studies (GWAS) have been successfully applied for investigating the basis of resistance traits, complex resistance phenotypes have been omitted so far. For S. aureus this especially refers to antibiotics of last resort like daptomycin and ceftaroline. Therefore, we aimed to perform GWAS for the identification of genetic variants associated with DAP and CPT resistance in clinical S. aureus isolates.

Materials/methods

To conduct microbial GWAS, we selected cases and controls according to their clonal background, date of isolation, and geographical origin. Association testing was performed with PLINK and SEER analysis. By using in silico analysis, we also searched for rare genetic variants in candidate loci that have previously been described to be involved in the development of corresponding resistance phenotypes.

Results

GWAS revealed MprF P314L and L826F to be significantly associated with DAP resistance. These mutations were found to be homogenously distributed among clonal lineages suggesting convergent evolution. Additionally, rare and yet undescribed single nucleotide polymorphisms could be identified within mprF and putative candidate genes. Finally, we could show that each DAP resistant isolate exhibited at least one amino acid substitution within the open reading frame of mprF. Due to the presence of strong population stratification, no genetic variants could be associated with CPT resistance. However, the investigation of the staphylococcal cassette chromosome mec (SCCmec) revealed various mecA SNPs to be putatively linked with CPT resistance. Additionally, some CPT resistant isolates revealed no mecA mutations, supporting the hypothesis that further and still unknown resistance determinants are crucial for the development of CPT resistance in S. aureus.

Conclusion

We hereby confirmed the potential of GWAS to identify genetic variants that are associated with antibiotic resistance traits in S. aureus. However, precautions need to be taken to prevent the detection of spurious associations. In addition, the implementation of different approaches is still essential to detect multiple forms of variations and mutations that occur with a low frequency.

T4-like Bacteriophages Isolated from Pig Stools Infect Yersinia pseudotuberculosis and Yersinia pestis Using LPS and OmpF as Receptors

The Yersinia bacteriophages fPS-2, fPS-65, and fPS-90, isolated from pig stools, have long contractile tails and elongated heads, and they belong to genus Tequatroviruses in the order Caudovirales. The phages exhibited relatively wide host ranges among Yersinia pseudotuberculosis and related species. One-step growth curve experiments revealed that the phages have latent periods of 50-80 min with burst sizes of 44-65 virions per infected cell. The phage genomes consist of circularly permuted dsDNA of 169,060, 167,058, and 167,132 bp in size, respectively, with a G + C content 35.3%. The number of predicted genes range from 267 to 271. The phage genomes are 84-92% identical to each other and ca 85% identical to phage T4. The phage receptors were identified by whole genome sequencing of spontaneous phage-resistant mutants. The phage-resistant strains had mutations in the ompF, galU, hldD, or hldE genes. OmpF is a porin, and the other genes encode lipopolysaccharide (LPS) biosynthetic enzymes. The ompF, galU, and hldE mutants were successfully complemented in trans with respective wild-type genes. The host recognition was assigned to long tail fiber tip protein Gp38, analogous to that of T-even phages such as Salmonella phage S16, specifically to the distal β-helices connecting loops.

Genomic signatures of adaptation to natural settings in non-typhoidal Salmonella enterica Serovars Saintpaul, Thompson and Weltevreden

Salmonella enterica is a pathogenic bacterium responsible for intestinal illness and systemic diseases such as typhoid and paratyphoid fevers. Among clinical manifestation classification, non-typhoidal Salmonella is mainly known as foodborne pathogen associated with the consumption of fecal contaminated food and water. Even though Salmonella hosts include humans and warm-blooded animals, it has been found in non-host environments as river water where the bacteria use different strategies to fitness the environment persisting and establishment. Now with the availability of WGS and bioinformatics tools, we can explore bacterial genomes with higher resolution to increase our understanding of specific genetic signatures among environmental and clinical isolates, being the goal of this work. Pangenome construction allowed the detection of specific environmental and clinical gene clusters related to metabolism and secretion systems as the main signature respectively. Specifically, D-galactonate degradation pathway was observed mainly in environmental genomes while T3SS and flagellum genes were detected for all clinical but not for all environmental isolates. Gene duplication and pseudogenes accumulation were detected as the main adaptation strategy for environmental isolates; thus, isolation source may play an important role in genome plasticity, conferring a selective advantage to survive and persist for environmental Salmonella isolates. Intact prophage sequences with cargo genes were observable for both isolation sources playing an important role in virulence contribution.

A Large Tn7-like Transposon Confers Hyper-Resistance to Copper in Pseudomonas syringae pv. syringae

Copper resistance mechanisms provide an important adaptive advantage to plant pathogenic bacteria under exposure to copper treatments. Copper resistance determinants have been described in Pseudomonas syringae pv. syringae (Pss) strains isolated from mango intimately associated with 62 kb plasmids belonging to the pPT23A family (PFP). It has been previously described that the indiscriminate use of copper-based compounds promotes the selection of copper resistant bacterial strains and constitutes a selective pressure in the evolution of copper resistance determinants. Hence, we have explored in this study the copper resistance evolution and the distribution of specific genetic determinants in two different Pss mango populations isolated from the same geographical regions, mainly from southern Spain with an average of 20 years of difference. The total content of plasmids, in particular the 62 kb plasmids, and the number of copper resistant Pss strains were maintained at similar levels over the time. Interestingly, the phylogenetic analysis indicated the presence of a phylogenetic subgroup (PSG) in the Pss mango phylotype, mostly composed of the recent Pss population analyzed in this study that was strongly associated with a hyper-resistant phenotype to copper. Genome sequencing of two selected Pss strains from this PSG revealed the presence of a large Tn7-like transposon of chromosomal location, which harbored putative copper and arsenic resistance genes (COARS Tn7-like). Transformation of the copper sensitive Pss UMAF0158 strain with some putative copper resistance genes and RT-qPCR experiments brought into light the role of COARS Tn7-like transposon in the hyper-resistant phenotype to copper in Pss.IMPORTANCECopper compounds have traditionally been used as standard bactericides in agriculture in the past few decades. However, the extensive use of copper has fostered the evolution of bacterial copper resistance mechanisms. Pseudomonas syringae is a plant pathogenic bacterium used worldwide as a model to study plant-pathogen interactions. The adaption of P. syringae to plant surface environment is the most important step prior to an infection. In this scenario, copper resistance mechanisms could play a key role in improving its epiphytic survival. In this work, a novel Tn7-like transposon of chromosomal location was detected in P. syringae pv. syringae strains isolated from mango. This transposon conferred the highest resistance to copper sulfate described to date for this bacterial phytopathogen. Understanding in depth the copper resistance mechanisms and their evolution are important steps to the agricultural industry to get a better improvement of disease management strategies.

Analysis of complete mitochondrial genome sequence of bar-tailed Treecreeper certhia himalayana (psittaciformes: Certhiidae)

Bar-tailed Treecreeper Certhia himalayana usually lives in coniferous or mixed broadleaf-conifer forests, often crawling along the trunk. In this study, we first sequenced and described the complete mitochondrial genome and phylogeny of C. himalayana. The whole genome of C. himalayana was 16,852 bp in length, and contained 13 protein-coding genes, 22 transfer RNA genes, 2 ribosome RNA genes, and 1 non-coding control regions. The overall base composition of the mitochondrial DNA was 25.1% for A, 29.2% for T, 14.5% for C, 31.2% for G, with a GC content of 45.7%. A phylogenetic tree strongly supported that C. himalayana closely related with Family Troglodytidae by highly probability.

Complete mitochondrial genome of Promethis valgipes valgipes (Marseul) (Insecta: Coleoptera: Tenebrionidae)

Promethis valgipes valgipes (Marseul) is one of the important fungus-eating beetle distributed in central China, Korea, and Japan. Beetles were obtained from Pizhou City and the species’ mitochondrial genome was characterized (GenBank accession number MW201671). The mitogenome consists of a circular DNA molecule of 15,801 bp, with 68.51% AT content. It comprises 13 protein-coding genes (PCGs), 22 tRNA genes, and two rRNA genes. The PCGs have typical ATN (Met) start codons, except nad1 (TTG as start codon), and are terminated by typical TAN stop codons.

The mitochondrial genome of Ophiostoma himal-ulmi and comparison with other fungi causing Dutch elm disease

The mitochondrial genome of Ophiostoma himal-ulmi, a species endemic to the Western Himalayas and one of the fungi that cause Dutch elm disease, has been sequenced and characterized. The mitochondrial genome was compared with other available genomes for members of the Ophiostomatales, including other agents of Dutch elm disease (Ophiostoma ulmi, Ophiostoma novo-ulmi subspecies novo-ulmi, and Ophiostoma novo-ulmi subspecies americana), and it was observed that gene synteny is highly conserved, and variability among members of the fungi that cause Dutch-elm disease is primarily due to the number of intron insertions. Among the fungi that cause Dutch elm disease that we examined, O. himal-ulmi has the largest mitochondrial genomes (ranging from 94 934 to 111 712 bp), owing to the expansion of the number of introns.

The Mitogenomes of Ophiostoma minus and Ophiostoma piliferum and Comparisons With Other Members of the Ophiostomatales

Fungi assigned to the Ophiostomatales are of economic concern as many are blue-stain fungi and some are plant pathogens. The mitogenomes of two blue-stain fungi, Ophiostoma minus and Ophiostoma piliferum, were sequenced and compared with currently available mitogenomes for other members of the Ophiostomatales. Species representing various genera within the Ophiostomatales have been examined for gene content, gene order, phylogenetic relationships, and the distribution of mobile elements. Gene synteny is conserved among the Ophiostomatales but some members were missing the atp9 gene. A genome wide intron landscape has been prepared to demonstrate the distribution of the mobile genetic elements (group I and II introns and homing endonucleases) and to provide insight into the evolutionary dynamics of introns among members of this group of fungi. Examples of complex introns or nested introns composed of two or three intron modules have been observed in some species. The size variation among the mitogenomes (from 23.7 kb to about 150 kb) is mostly due to the presence and absence of introns. Members of the genus Sporothrix sensu stricto appear to have the smallest mitogenomes due to loss of introns. The taxonomy of the Ophiostomatales has recently undergone considerable revisions; however, some lineages remain unresolved. The data showed that genera such as Raffaelea appear to be polyphyletic and the separation of Sporothrix sensu stricto from Ophiostoma is justified.

Endophytic Lifestyle of Global Clones of Extended-Spectrum β-Lactamase-Producing Priority Pathogens in Fresh Vegetables: a Trojan Horse Strategy Favoring Human Colonization?

The global spread of antibiotic-resistant bacteria and their resistance genes is a critical issue that is no longer restricted to hospital settings, but also represents a growing problem involving environmental and food safety. In this study, we have performed a microbiological and genomic investigation of critical priority pathogens resistant to broad-spectrum cephalosporins and showing endophytic lifestyles in fresh vegetables sold in a country with high endemicity of extended-spectrum β-lactamases (ESBLs). We report the isolation of international high-risk clones of CTX-M-15-producing Escherichia coli, belonging to clonal complexes CC38 and CC648, and Klebsiella pneumoniae of complex CC307 from macerated tissue of surface-sterilized leaves of spinach, cabbage, arugula, and lettuce. Regardless of species, all ESBL-positive isolates were able to endophytically colonize common bean (Phaseolus vulgaris) seedlings, showed resistance to acid pH, and had a multidrug-resistant (MDR) profile to clinically relevant antibiotics (i.e., broad-spectrum cephalosporins, aminoglycosides, and fluoroquinolones). Genomic analysis of CTX-M-producing endophytic Enterobacterales revealed a wide resistome (antibiotics, biocides, disinfectants, and pesticides) and virulome, and genes for endophytic fitness and for withstanding acidic conditions. Transferable IncFIB and IncHI2A plasmids carried bla _CTX-M-15 genes and, additionally, an IncFIB plasmid (named pKP301cro) also harbored genes encoding resistance to heavy metals. These data support the hypothesis that fresh vegetables marketed for consumption can act as a figurative Trojan horse for the hidden spread of international clones of critical WHO priority pathogens producing ESBLs, and/or their resistance genes, to humans and other animals, which is a critical issue within a food safety and broader public and environmental health perspective.IMPORTANCE Extended-spectrum β-lactamases (ESBL)-producing Enterobacterales are a leading cause of human and animal infections, being classified as critical priority pathogens by the World Health Organization. Epidemiological studies have shown that spread of ESBL-producing bacteria is not a problem restricted to hospitals, but also represents a growing problem involving environmental and food safety. In this regard, CTX-M-type β-lactamases have become the most widely distributed and clinically relevant ESBLs worldwide. Here, we have investigated the occurrence and genomic features of ESBL-producing Enterobacterales in surface-sterilized fresh vegetables. We have uncovered that international high-risk clones of CTX-M-15-producing Escherichia coli and Klebsiella pneumoniae harboring a wide resistome and virulome, carry additional genes for endophytic fitness and resistance to acidic conditions. Furthermore, we have demonstrated that these CTX-M-15-positive isolates are able to endophytically colonize plant tissues. Therefore, we believe that fresh vegetables can act as a figurative Trojan horse for the hidden spread of critical priority pathogens exhibiting endophytic lifestyles.

Draft Genome Sequence of Comamonas aquatilis Strain LK (= CSUR P6418 = CECT 9772), Isolated from the Planarian Schmidtea mediterranea

Comamonas aquatilis was defined as a new Comamonas species based on its 16S rRNA sequence, but the genome from the type strain SB30-Cr27-3^T (= CIP 111491^T = CCM 8815^T) is not available. We have cultivated from the planarian Schmidtea mediterranea a Comamonas aquatilis strain, LK (= CSUR P6418 = CECT 9772), that exhibits 100% 16S rRNA sequence similarity to strain SB30-Cr27-3^T We have sequenced the genome of strain LK and obtained a chromosome of 4,899,818 bp, with a G+C content of 61.75%, assembled into two contigs.

The complete mitochondrial genome of ribbed gunnel Dictyosoma burgeri (van der Hoeven, 1855)

The complete mitochondrial genome of Dictyosoma burgeri collected from Yellow and Bohai Seas was determined by next-generation sequencing. The mitogenome is a circular molecule 16,513 bp in length, including the typical structure of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. The TAS, central CSB, and CSB were detected in the control region. The gene contents of the mitogenome are identical to those observed in most bony fishes.

The complete chloroplast genome of Isodon rubescens, a traditional Chinese herb

Isodon rubescens (Hemsley) H. Hara is a kind of traditional medical herb which can be used for cancer treatment. In this study, the complete chloroplast genome sequence of I. rubescens was assembled. Its complete circular chloroplast DNA length was 152,761 bp. The genome was made up of a large single-copy region of 83,655 bp, a small single-copy region of 17,660 bp, and a pair of inverted repeat regions of 25,723 bp. The genome totally encoded 129 genes, containing 85 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenetic analysis indicated that I. rubescens had a close relationship with basil (Ocimum basilicum).

Draft Genome Sequences of Dysgonomonas sp. Strains BGC7 and HGC4, Isolated from the Hindgut of a Lower Termite

Dysgonomonas spp. are facultative heterotrophs which colonize diverse environments, including the hindgut of the lower termite Reticulitermes flavipes. Dysgonomonas genomes are enriched for genes involving oligo- and polysaccharide utilization, enabling modification of a wide array of complex glycans. Here, we report draft genome sequences for Dysgonomonas sp. strains BGC7 and HGC4.

Dysgonomonas spp. are facultative heterotrophs which colonize diverse environments, including the hindgut of the lower termite Reticulitermes flavipes. Dysgonomonas genomes are enriched for genes involving oligo- and polysaccharide utilization, enabling modification of a wide array of complex glycans. Here, we report draft genome sequences for Dysgonomonas sp. strains BGC7 and HGC4.

The complete chloroplast genome of Chinese medicine cultivar species of Rehmannia glutinosa (Orobanchaceae)

The cultivar of Rehmannia glutinosa (Orobanchaceae) is one of the four famous ‘Huai’ medicine cultivar species endemic to Henan Province in central China. In this study, we report the complete chloroplast (cp) genome of R. glutinosa cultivar Wen 85-5. The cp genome of R. glutinosa cultivar Wen 85-5 was 155,499 bp in length and contained a pair of inverted repeat regions (IR, 25,748 bp) separated by a small single copy (SSC, 17,600 bp) and a large single copy (LSC, 84,403 bp) region. Chloroplast genome sequences of two cultivar of R. glutinosa (Wen 85-5 and Jiwang 1) are identical to each other. The sequence of cp genome of R. glutinosa cultivar Wen 85-5 was 99.70% similar to the wild population of R. glutinosa. Some distinctive insert and deletion in R. glutinosa cultivar Wen 85-5 by comparison with wild population were reported. The maximum-likelihood phylogenetic analysis revealed that R. glutinosa cultivar Wen 85-5.was sister to the R. glutinosa cultivar Jiwang 1 (BS = 100%), and further clustering with R. glutinosa (BS = 100%). This result will be helpful for the conservation and breeding programs of the cultivar of R. glutinosa.