CMG-biotools, a free workbench for basic comparative microbial genomics

Metagenomic exploration of cold-active enzymes for detergent applications: Characterization of a novel, cold-active and alkali-stable GH8 endoglucanase from ikaite columns in SW Greenland

Microbial communities from extreme environments are largely understudied, but are essential as producers of metabolites, including enzymes, for industrial processes. As cultivation of most microorganisms remains a challenge, culture-independent approaches for enzyme discovery in the form of metagenomics to analyse the genetic potential of a community are rapidly becoming the way forward. This study focused on analysing a metagenome from the cold and alkaline ikaite columns in Greenland, identifying 282 open reading frames (ORFs) that encoded putative carbohydrate-modifying enzymes with potential applications in, for example detergents and other processes where activity at low temperature and high pH is desired. Seventeen selected ORFs, representing eight enzyme families were synthesized and expressed in two host organisms, Escherichia coli and Aliivibrio wodanis. Aliivibrio wodanis demonstrated expression of a more diverse range of enzyme classes compared to E. coli, emphasizing the importance of alternative expression systems for enzymes from extremophilic microorganisms. To demonstrate the validity of the screening strategy, we chose a recombinantly expressed cellulolytic enzyme from the metagenome for further characterization. The enzyme, Cel240, exhibited close to 40% of its relative activity at low temperatures (4°C) and demonstrated endoglucanase characteristics, with a preference for cellulose substrates. Despite low sequence similarity with known enzymes, computational analysis and structural modelling confirmed its cellulase-family affiliation. Cel240 displayed activity at low temperatures and good stability at 25°C, activity at alkaline pH and increased activity in the presence of CaCl2, making it a promising candidate for detergent and washing industry applications.

Phylogenomic analysis of the Ponticus clade: strains isolated from the spotted rose snapper (Lutjanus guttatus)

The family Vibrionaceae is classified into many clades based on their phylogenetic relationships. The Ponticus clade is one of its clades and consists of four species, Vibrio panuliri, V. ponticus, V. rhodolitus, and V. taketomensis. Two strains, CAIM 703 and CAIM 1902, were isolated from the diseased spotted rose snapper external lesion (Lutjanus guttatus), they were analyzed to determine their taxonomic position, a phylogenetic analysis was performed based on the 16S rRNA sequences proved that the two strains are members of the genus Vibrio and they belong to the Ponticus clade. Then, a phylogenomic analysis was performed with four type strains and four reference strains isolated from marine organisms and aquatic environments. Multilocus Sequence Analysis (MLSA) of 139 single-copy genes showed that CAIM 703 and CAIM 1902 belong to V. panuliri. The 16S rRNA sequence similarity value between CAIM 703 and CAIM 1902 was 99.61%. The Ponticus clade species showed Average Nucleotide Identity (ANI) values between 78 to 80% against the two strains for ANIb, except V. panuliri LBS2T (99% and 100% similarity). Finally, this analysis represents the first phylogenomic analysis of the Ponticus clade where V. panuliri strains are reported from Mexico.

Probability of outbreaks and cross-border dissemination of the emerging pathogen: a genomic survey of Elizabethkingia meningoseptica

IMPORTANCE

Elizabethkingia meningoseptica is an emerging infectious agent associated with life-threatening infections in immunocompromised individuals. However, there are limited data available on the genomic features of E. meningoseptica. This study aims to characterize the geographical distribution, phylogenetic evolution, pathogenesis, and transmission of this bacterium. A systematic analysis of the E. meningoseptica genome revealed that a common ancestor of this bacterium existed 90 years ago. The evolutionary history showed no significant relationship with the sample source, origin, or region, despite the presence of genetic diversity. Whole genome sequencing data also demonstrated that E. meningoseptica bacteria possess inherent resistance and pathogenicity, enabling them to spread within the same hospital and even across borders. This study highlights the potential for E. meningoseptica to cause severe nosocomial outbreaks and horizontal transmission between countries worldwide. The available evidence is crucial for the development of evidence-based public health policies to prevent global outbreaks caused by emerging pathogens.

Phylogenomic Insights on the Xanthomonas translucens Complex, and Development of a TaqMan Real-Time Assay for Specific Detection of pv. translucens on Barley

The reemergence and spread of Xanthomonas translucens, the causal agent of bacterial leaf streak in cereal crops and wilt in turfgrass and forage species, is a concern to growers in the United States and Canada. The pathogen is seedborne and listed as an A2 quarantine organism by EPPO, making it a major constraint to international trade and exchange of germplasm. The pathovar concept of the X. translucens group is confusing due to overlapping of plant host ranges and specificity. Here, comparative genomics, phylogenomics, and 81 up-to-date bacterial core gene set (ubcg2) were used to assign the pathovars of X. translucens into three genetically and taxonomically distinct clusters. The study also showed that whole genome-based digital DNA-DNA hybridization unambiguously can differentiate the pvs. translucens and undulosa. Orthologous gene and proteome matrix analyses suggest that the cluster consisting of graminis, poae, arrhenatheri, phlei, and phleipratensis is very divergent. Whole-genome data were exploited to develop the first pathovar-specific TaqMan real-time PCR tool for detection of pv. translucens on barley. Specificity of the TaqMan assay was validated using 62 Xanthomonas and non-Xanthomonas strains as well as growth chamber-inoculated and naturally infected barley leaves. Sensitivity levels of 0.1 pg (purified DNA) and 23 CFUs per reaction (direct culture) compared favorably with other previously reported real-time PCR assays. The phylogenomics data reported here suggest that the clusters could constitute novel taxonomic units or new species. Finally, the pathovar-specific diagnostic tool will have significant benefits to growers and facilitate international exchange of barley germplasm and trade.

Genomic and phenotypic characterization of Pseudomonas sp. GOM7, a novel marine bacterial species with antimicrobial activity against multidrug-resistant Staphylococcus aureus

Antimicrobial resistance (AMR) represents a serious threat to global health. The development of new drugs to combat infections caused by bacteria resistant to multiple or even all available antibiotics is urgent. Most antibiotics used up to date have been identified from soil microorganisms. The marine environment represents an alternative source with great potential for the identification of microorganisms that produce bioactive molecules, including antibiotics. In this study, we analyzed the antibacterial activity of a collection of 82 bacterial strains isolated from marine water and sediment samples collected from the Southwestern Gulf of Mexico. Eight of the marine isolates inhibited the growth of different pathogenic bacteria, seven of which were identified as presumptive Pseudomonas aeruginosa. Interestingly, genome sequencing and phylogenetic analysis revealed that the remaining marine isolate showing antibacterial activity is a novel Pseudomonas species that we denominated Pseudomonas sp. GOM7, which was not pathogenic in the Galleria mellonella infection model in the conditions tested. Notably, Pseudomonas sp. GOM7 inhibited the growth of multidrug and methicillin-resistant strains of the priority pathogen Staphylococcus aureus. Our results show that the anti-S. aureus compound(s) produced by Pseudomonas sp. GOM7 can be extracted from the culture supernatant of this bacterium with the organic solvent ethyl acetate. Annotation of the Pseudomonas sp. GOM7 genome revealed the presence of several biosynthetic gene clusters predicted to code for possible antimicrobial compounds. Our results further highlight the potential of bacteria from the Gulf of Mexico as a source of novel antimicrobials.

Revisiting chloroplast genomic landscape and annotation towards comparative chloroplast genomes of Rhamnaceae

BACKGROUND

Massive parallel sequencing technologies have enabled the elucidation of plant phylogenetic relationships from chloroplast genomes at a high pace. These include members of the family Rhamnaceae. The current Rhamnaceae phylogenetic tree is from 13 out of 24 Rhamnaceae chloroplast genomes, and only one chloroplast genome of the genus Ventilago is available. Hence, the phylogenetic relationships in Rhamnaceae remain incomplete, and more representative species are needed.

RESULTS

The complete chloroplast genome of Ventilago harmandiana Pierre was outlined using a hybrid assembly of long- and short-read technologies. The accuracy and validity of the final genome were confirmed with PCR amplifications and investigation of coverage depth. Sanger sequencing was used to correct for differences in lengths and nucleotide bases between inverted repeats because of the homopolymers. The phylogenetic trees reconstructed using prevalent methods for phylogenetic inference were topologically similar. The clustering based on codon usage was congruent with the molecular phylogenetic tree. The groups of genera in each tribe were in accordance with tribal classification based on molecular markers. We resolved the phylogenetic relationships among six Hovenia species, three Rhamnus species, and two Ventilago species. Our reconstructed tree provides the most complete and reliable low-level taxonomy to date for the family Rhamnaceae. Similar to other higher plants, the RNA editing mostly resulted in converting serine to leucine. Besides, most genes were subjected to purifying selection. Annotation anomalies, including indel calling errors, unaligned open reading frames of the same gene, inconsistent prediction of intergenic regions, and misannotated genes, were identified in the published chloroplast genomes used in this study. These could be a result of the usual imperfections in computational tools, and/or existing errors in reference genomes. Importantly, these are points of concern with regards to utilizing published chloroplast genomes for comparative genomic analysis.

CONCLUSIONS

In summary, we successfully demonstrated the use of comprehensive genomic data, including DNA and amino acid sequences, to build a reliable and high-resolution phylogenetic tree for the family Rhamnaceae. Additionally, our study indicates that the revision of genome annotation before comparative genomic analyses is necessary to prevent the propagation of errors and complications in downstream analysis and interpretation.

Comparison of Monkeypox virus genomes from the 2017 Nigeria outbreak and the 2022 outbreak

AIMS

The current Monkeypox virus (MPX) outbreak is not only the largest known outbreak to date caused by a strain belonging to the West-African clade, but also results in remarkably different clinical and epidemiological features compared to previous outbreaks of this virus. Here, we consider the possibility that mutations in the viral genome may be responsible for its changed characteristics.

METHODS AND RESULTS

Six genome sequences of isolates from the current outbreak were compared to five genomes of isolates from the 2017 outbreak in Nigeria and to two historic genomes, all belonging to the West-African clade. We report differences that are consistently present in the 2022 isolates but not in the others. Although some variation in repeat units was observed, only two were consistently found in the 2022 genomes only, and these were located in intergenic regions. A total of 55 single nucleotide polymorphisms were consistently present in the 2022 isolates compared to the 2017 isolates. Of these, 25 caused an amino acid substitution in a predicted protein.

CONCLUSIONS

The nature of the substitution and the annotation of the affected protein identified potential candidates that might affect the virulence of the virus. These included the viral DNA helicase and transcription factors.

SIGNIFICANCE

This bioinformatic analysis provides guidance for wet-lab research to identify changed properties of the MPX.

First Record of the Rare Species Aeromonas lusitana from Rainbow Trout (Oncorhynchus mykiss, Walbaum): Comparative Analysis with the Existing Strains

The species Aeromonas lusitana was first described in 2016 with five strains recovered from untreated water and vegetables from Portugal. Since then, no further records exist of this species. During a surveillance study on the presence of Aeromonas in fish farms in Mexico, a new strain (ESV-351) of the mentioned species isolated from a rainbow trout was recovered. It was identified because it clustered phylogenetically with the type strain of A. lusitana based on the analysis of the rpoD gene sequences. In the present study, phenotypic characteristics, antimicrobial resistance profiles, and the presence of putative virulence genes of this novel strain (ESV-351) were determined in parallel to the five isolates from the original species description. Phenotypic differential characteristics exhibited by A. lusitana ESV-351 depicted an evident similarity to the characteristics exhibited by the other evaluated strains. However, the novel strain was positive for the production of indole using conventional methods, while the rest of the strains, including the type strain, were negative for its production. Furthermore, intermediate resistance to ampicillin, amoxicillin-clavulanic acid and cephalothin was detected in both the novel and the type strain. Five different virulence-related genes were detected in the novel strain and in the previously described strains, with the type strain exhibiting the highest number of virulence-related genes. In addition to this, the genome of the novel strain (ESV-351) was sequenced and compared with the genomes from the type strain (A. lusitana CECT 7828T) and other Aeromonas spp. The genomic analysis defined Aeromonas tecta as the closest species to A. lusitana with a highly similar number of predicted proteins. The genomic size, the number of protein-encoding genes and the number of different tRNAs, among other characteristics, make it possible to propose that the ESV-351 strain could potentially have the capacity to adapt to different environments. Genome comparison of the ESV-351 strain with the type strain revealed that both possess a similar sequence of the citrate synthase gene. In addition to this finding, the chromosomal region containing the citrate synthase locus of the novel strain exhibits some similarity to the chromosomal region in the genome of the A. hydrophila type strain and other known human pathogens, such as Vibrio cholerae. This could suggest a possible virulence role for the citrate synthase gene in A. lusitana (ESV-351).

Deep Learning for Clothing Style Recognition Using YOLOv5

With the rapid development of artificial intelligence, much more attention has been paid to deep learning. However, as the complexity of learning algorithms increases, the needs of computation power of hardware facilities become more crucial. Instead of the focus being on computing devices like GPU computers, a lightweight learning algorithm could be the answer for this problem. Cross-domain applications of deep learning have attracted great interest amongst researchers in academia and industries. For beginners who do not have enough support with software and hardware, an open-source development environment is very helpful. In this paper, a relatively lightweight algorithm YOLOv5s is addressed, and the Google Colab is used for model training and testing. Based on the developed environment, many state-of-art learning algorithms can be studied for performance comparisons. To highlight the benefits of one-stage object detection algorithms, the recognition of clothing styles is investigated. The image samples are selected from datasets of fashion clothes and the web crawling of online stores. The image data are categorized into five groups: plaid; plain; block; horizontal; and vertical. Average precison, mean average precison, recall, F1-score, model size, and frame per second are the metrics used for performance validations. From the experimental outcomes, it shows that YOLOv5s is better than other learning algorithms in the recognition accuracy and detection speed.

Whole-genome sequence analysis for evaluating the safety and probiotic potential of Lactiplantibacillus pentosus 9D3, a gamma-aminobutyric acid (GABA)-producing strain isolated from Thai pickled weed

Lactiplantibacillus pentosus 9D3, a prominent gamma-aminobutyric acid (GABA)-producing bacteria isolated from Thai pickled weed was characterized for its safety and probiotic properties via whole-genome analysis and in vitro testing. The whole-genome sequence of L. pentosus 9D3 was determined using a hybrid-sequencing approach, combining PacBio and Illumina technologies. A 3.81-Mbp genome of L. pentosus 9D3 consisting of one 3.65-Mbp chromosome and six plasmids (1.9–71.9 Kbp) was identified with an estimated GC content of 46.09% and 3,456 predicted genes. The strain was confirmed to be Lactiplantibacillus pentosus according to the high average nucleotide identity value of >95% and digital DNA–DNA hybridization scores of >70% to the L. pentosus type strain. Comparative genome analysis with other L. pentosus strains showed that the GABA-producing capability was specific to the strain 9D3. Genes related to GABA biosynthesis and transport were identified on a plasmid, pLPE-70K, indicating the acquired nature of this property. The safety of L. pentosus 9D3 was demonstrated through the lack of genes related to the production of toxins, biogenic amines, and antimicrobial drugs. Although the strain exhibited resistance to ampicillin and chloramphenicol, none of the antimicrobial resistance (AMR) genes were associated with mobile elements, i.e., plasmids and prophages. Therefore, the strain is considered to have low risk of transferring the AMR genes to other, potentially pathogenic bacteria. In addition, L. pentosus 9D3 showed good survivability in the gastrointestinal tract environment and was able to adhere to the intestinal cell in vitro. Therefore, L. pentosus 9D3 is concluded to be safe, with the potential to be used as a probiotic, exerting its health benefit through GABA production in the food system. The GABA-producing capability of the strain in vivo is the subject of further investigation.

Global genomic and proteomic analysis indicates co-evolution of Neisseria species and with their human host

Neisseria, a genus from the beta-proteobacteria class, is of potential clinical importance. This genus contains both pathogenic and commensal strains. Gonorrhea and meningitis are two major diseases caused by pathogens belonging to this genus. With the increased use of antimicrobial agents against these pathogens they have evolved the antimicrobial resistance capacity making these diseases nearly untreatable. The set of anti-bacterial resistance genes (resistome) and genes associated with signal processing (secretomes) are crucial for the host-microbial interaction. With the virtue of whole-genome sequences and computational biology, it is now possible to study the genomic and proteomic riddles of Neisseria along with their comprehensive evolutionary and metabolic profiling. We have studied relative synonymous codon usage, amino acid usage, reverse ecology, comparative genomics, evolutionary analysis and pathogen-host (Neisseria-human) interaction through bioinformatics analysis. Our analysis revealed the co-evolution of Neisseria genomes with the human host. Moreover, the co-occurrence of Neisseria and humans has been supported through reverse ecology analysis. A differential pattern of the evolutionary rate of resistomes and secretomes was evident among the pathogenic and commensal strains. Comparative genomics supported the presence of virulent genes in both pathogenic and commensal strains of the select genus. Our analysis also indicated a transition from commensal to pathogenic Neisseria strains through the long run of evolution.

Whole-Genome Resources and Species-Level Taxonomic Validation of 89 Plant-Pathogenic Xanthomonas Strains Isolated from Various Host Plants

Bacterial spot disease caused by Xanthomonas spp. is a global threat to tomato and pepper plants. A recent classification of these pathogens indicated the need for a diverse dataset of whole-genome resources. We report whole-genome resources of 89 Xanthomonas strains isolated from Canada (n = 44), the United States (n = 29), Argentina (n = 4), Brazil (n = 3), Costa Rica (n = 3), New Zealand (n = 1), Australia (n = 1), Mexico (n = 1), Taiwan (n = 1), Thailand (n = 1), and unknown (n = 1). Of these strains, 48 were previously identified to species-level based on nongenome-based approaches while 41 strains were classified only at the genus level. The average coverage of the sequencing reads was 103×. The draft genome sizes ranged from 4.53 to 5.46 Mbp with a G + C content of 63.53 to 67.78% and comprised 4,233-5,178 protein-coding sequences. Using average nucleotide identity (ANI) and genome-based DNA-DNA hybridization (gDDH) values, the taxonomic classifications were validated for 38 of the 48 strains previously assigned to species level using other methods. Ten strains previously identified as Xanthomonas campestris, X. axonopodis, X. vasicola, and X. arboricola were incorrectly assigned, and new species-level delineations are proposed. Data from ANI, gDDH, and pangenome phylogeny of shared protein families were used to assign the 41 strains, previously identified only to genus level, into five distinct species: X. euvesicatoria (pv. euvesicatoria or pv. perforans), X. hortorum pv. gardneri, X. vesicatoria, X. campestris, and X. arboricola. These 89 whole-genome sequences of Xanthomonas strains, the majority (49.4%) of which are from Canada, could be useful resources in our understanding of the global population structure and evolution of these pathogens.

The phylogeny of the genus Azohydromonas supports its transfer to the family Comamonadaceae

The genus Azohydromonas encompasses five validly described species belonging to the betaproteobacterial class. Recognized for their potential biotechnological uses, they were first described as belonging to the genus Alcaligenes. The phylogeny of the 16S rRNA gene of the original strains as well as newly described species led to a description of these strains within a new bacterial genus, Azohydromonas. However, the phylogenetic position of this genus remains described as part of the family Alcaligenaceae, even those some authors have placed it within the order Burkholderiales. To unravel the precise position of the genus Azohydromonas, a wide phylogenomic analysis was performed. The results of 16S rRNA gene phylogeny, as well as those obtained by the multilocus analysis of homologous proteins and overall genome relatedness indices, support the reclassification of Azohydromonas in the Rubrivivax-Ideonella lineage of the family Comamonadaceae, so the transfer of this genus is proposed.

Comparative genomics of Photobacterium species from terrestrial and marine habitats

Photobacterium (P.) is a genus widely studied in regards to its association with and ubiquitous presence in marine environments. However, certain species (P. phosphoreum, P. carnosum, P. iliopiscarium) have been recently described to colonize and spoil raw meats without a marine link. We have studied 27 strains from meat as well as 26 strains from marine environments in order to probe for intraspecies marine/terrestrial subpopulations and identify distinct genomic features acquired by environmental adaptation. We have conducted phylogenetic analysis (MLSA, ANI, fur, codon usage), search of plasmids (plasmidSPADES), phages (PHASTER), CRISPR-cas operons (CRISPR-finder) and secondary metabolites gene clusters (antiSMASH, BAGEL), in addition to a targeted gene search for specific pathways (e.g. TCA cycle, pentose phosphate, respiratory chain) and elements relevant for growth, adaptation and competition (substrate utilization, motility, bioluminescence, sodium and iron transport). P. carnosum appears as a conserved single clade, with one isolate from MAP fish clustering apart that doesn't, however, show distinct features that could indicate different adaptation. The species harbors genes for a wide carbon source utilization (glycogen/starch, maltose, pullulan, fucose) for colonization of diverse niches in its genome. P. phosphoreum is represented by two different clades on the phylogenetic analyses not correlating to their origin or distribution of other features analyzed that can be divided into two novel subspecies based on genome-wide values. A more diverse antimicrobial activity (sactipeptides, microcins), production of secondary metabolites (siderophores and arylpolyenes), stress response and adaptation (bioluminescence, sodium transporters, catalase, high affinity for oxygen cytochrome cbb3 oxidase, DMSO reductase and proton translocating NADH dehydrogenase) is predicted compared to the other species. P. iliopiscarium was divided into two clades based on source of isolation correlating with phylogeny and distribution of several traits. The species shows traits common to the other two species, similar carbon utilization/transport gene conservation as P. carnosum for the meat-isolated strains, and predicted utilization of marine-common DMSO and flagellar cluster for the sea-isolated strains. Results additionally suggest that photobacteria are highly prone to horizontal acquisition/loss of genetic material and genetic transduction, and that it might be a strategy for increasing the frequency of strain- or species-specific features that offers a growth/competition advantage.

Genomics and phylogeny of the proposed phylum 'Candidatus Poribacteria' associated with the excavating sponge Thoosa mismalolli

Members of the proposed phylum 'Candidatus Poribacteria' are among the most abundant microorganisms in the highly diverse microbiome of the sponge mesohyl. Genomic and phylogenetic characteristics of this proposed phylum are barely known. In this study, we analyzed metagenome-assembled genomes (MAGs) obtained from the coral reef excavating sponge Thoosa mismalolli from the Mexican Pacific Ocean. Two MAGs were extracted and analyzed together with 32 MAGs and single-amplified genomes (SAGs) obtained from NCBI. The phylogenetic tree based on the sequences of 139 single-copy genes (SCG) showed two clades. Clade A (23 genomes) represented 67.7% of the total of the genomes, while clade B (11 genomes) comprised 32.3% of the genomes. The Average Nucleotide Identity (ANI) showed values between 66 and 99% for the genomes of the proposed phylum, and the pangenome of genomes revealed a total of 37,234 genes that included 1722 core gene. The number of genes used in the phylogenetic analysis increased from 28 (previous studies) to 139 (this study), which allowed a better resolution of the phylogeny of the proposed phylum. The results supported the two previously described classes, 'Candidatus Entoporibacteria' and 'Candidatus Pelagiporibacteria', and the genomes SB0101 and SB0202 obtained in this study belong to two new species of the class 'Candidatus Entoporibacteria'. This is the first comparative study that includes MAGs from a non-sponge host (Porites lutea) to elucidate the taxonomy of the poorly known Candidatus phylum in a polyphasic approach. Finally, our study also contributes to the sponge microbiome project by reporting the first MAGs of the proposed phylum 'Candidatus Poribacteria' isolated from the excavating sponge T. mismalolli.

Genomic stability among O3:K6 V. parahaemolyticus pandemic strains isolated between 1996 to 2012 in American countries

BACKGROUND

The V. parahaemolyticus pandemic clone, results in the development of gastrointestinal illness in humans. Toxigenic strains of this species are frequently isolated from aquatic habitats and organisms such as mollusks and crustaceans. Reports on the isolation of the pandemic clone started in 1996, when a new O3:K6 clone was identified in Asia, that rapidly spread worldwide, becoming the predominant clone isolated from clinical cases. In this study whole genome sequencing was accomplished with an Illumina MiniSeq platform, upon six novel V. parahaemolyticus strains, that have been isolated in Mexico since 1998 and three representative genomes of strains that were isolated from reported outbreaks in other American countries, and were deposited in the GenBank. These nine genomes were compared against the reference sequence of the O3:K6 pandemic strain (RIMD 2210633), which was isolated in 1996, to determine sequence differences within American isolates and between years of isolation.

RESULTS

The results indicated that strains that were isolated at different times and from different countries, were highly genetically similar, among them as well as to the reference strain RIMD 2210633, indicating a high level of genetic stability among the strains from American countries between 1996 to 2012, without significant genetic changes relative to the reference strain RIMD 2210633, which was isolated in 1996 and was considered to be representative of a novel O3:K6 pandemic strain.

CONCLUSIONS

The genomes of V. parahaemolyticus strains isolated from clinical and environmental sources in Mexico and other American countries, presented common characteristics that have been reported for RIMD 2210633 O3:K6 pandemic strain. The major variations that were registered in this study corresponded to genes non associated to virulence factors, which could be the result of adaptations to different environmental conditions. Nevertheless, results do not show a clear pattern with the year or locality where the strains were isolated, which is an indication of a genomic stability of the studied strains.

Genomic Diversity of Pigeon Pea (Cajanus cajan L. Millsp.) Endosymbionts in India and Selection of Potential Strains for Use as Agricultural Inoculants

Pigeon pea (Cajanus cajan L. Millsp. ) is a legume crop resilient to climate change due to its tolerance to drought. It is grown by millions of resource-poor farmers in semiarid and tropical subregions of Asia and Africa and is a major contributor to their nutritional food security. Pigeon pea is the sixth most important legume in the world, with India contributing more than 70% of the total production and harbouring a wide variety of cultivars. Nevertheless, the low yield of pigeon pea grown under dry land conditions and its yield instability need to be improved. This may be done by enhancing crop nodulation and, hence, biological nitrogen fixation (BNF) by supplying effective symbiotic rhizobia through the application of "elite" inoculants. Therefore, the main aim in this study was the isolation and genomic analysis of effective rhizobial strains potentially adapted to drought conditions. Accordingly, pigeon pea endosymbionts were isolated from different soil types in Southern, Central, and Northern India. After functional characterisation of the isolated strains in terms of their ability to nodulate and promote the growth of pigeon pea, 19 were selected for full genome sequencing, along with eight commercial inoculant strains obtained from the ICRISAT culture collection. The phylogenomic analysis [Average nucleotide identity MUMmer (ANIm)] revealed that the pigeon pea endosymbionts were members of the genera Bradyrhizobium and Ensifer. Based on nodC phylogeny and nod cluster synteny, Bradyrhizobium yuanmingense was revealed as the most common endosymbiont, harbouring nod genes similar to those of Bradyrhizobium cajani and Bradyrhizobium zhanjiangense. This symbiont type (e.g., strain BRP05 from Madhya Pradesh) also outperformed all other strains tested on pigeon pea, with the notable exception of an Ensifer alkalisoli strain from North India (NBAIM29). The results provide the basis for the development of pigeon pea inoculants to increase the yield of this legume through the use of effective nitrogen-fixing rhizobia, tailored for the different agroclimatic regions of India.

Veronia nyctiphanis gen. nov., sp. nov., Isolated from the Stomach of the Euphausiid Nyctiphanes simplex (Hansen, 1911) in the Gulf of California, and Reclassification of Enterovibrio pacificus as Veronia pacifica comb. nov

The bacterial strain 42Xb2 ^T was isolated from a female adult krill Nyctiphanes simplex infected with the apostome parasitoid ciliate Pseudocollinia brintoni in January 2007 in the Gulf of California. The strain has the morphological, phenotypic, and molecular characteristics of the bacteria of the family Vibrionaceae. The 16S rRNA gene sequence has a similarity of 97.7% with Enterovibrio pacificus SW014 ^T and 96.1% similarity with Enterovibrio norvegicus LMG 19839 ^T. A phylogenomic and a multilocus sequence analyses placed this strain close to the genera Enterovibrio, Grimontia, and Salinivibrio, but clearly forming a separate branch from these bacterial genera. Genomic analyses presented further support this result. A novel genus Veronia gen. nov. and a species Veronia nyctiphanis sp. nov. is here described with CAIM 600 ^T (= DSM 24592 ^T = CECT 7578 ^T) as the type strain. Morphological, physiological, and genetic evidence presented here support the unification of Enterovibrio pacificus and Veronia nyctiphanis in the new genus Veronia. Enterovibrio pacificus is reclassified as Veronia pacifica. V. pacifica is assigned as the type species of the new genus Veronia.Genome Sequencing Data The GenBank/EMBL/DDBJ accession numbers for the genome sequence of Veronia nyctiphanis CAIM 600 ^T is PEIB01 and of Enterovibrio pacificus CAIM 1920 ^T is LYBM01. The 16S rRNA gene sequence of V. nyctiphanis CAIM 600 ^T is JX129353.

Genomic and Phenotypic Biology of Novel Strains of Dickeya zeae Isolated From Pineapple and Taro in Hawaii: Insights Into Genome Plasticity, Pathogenicity, and Virulence Determinants

Dickeya zeae, a bacterial plant pathogen of the family Pectobacteriaceae, is responsible for a wide range of diseases on potato, maize, rice, banana, pineapple, taro, and ornamentals and significantly reduces crop production. D. zeae causes the soft rot of taro (Colocasia esculenta) and the heart rot of pineapple (Ananas comosus). In this study, we used Pacific Biosciences single-molecule real-time (SMRT) sequencing to sequence two high-quality complete genomes of novel strains of D. zeae: PL65 (size: 4.74997 MB; depth: 701x; GC: 53.6%) and A5410 (size: 4.7792 MB; depth: 558x; GC: 53.5%) isolated from economically important Hawaiian crops, taro, and pineapple, respectively. Additional complete genomes of D. zeae representing three additional hosts (philodendron, rice, and banana) and other species used for a taxonomic comparison were retrieved from the NCBI GenBank genome database. Genomic analyses indicated the truncated type III and IV secretion systems (T3SS and T4SS) in the taro strain, which only harbored one and two genes of T3SS and T4SS, respectively, and showed high heterogeneity in the type VI secretion system (T6SS). Unlike strain EC1, which was isolated from rice and recently reclassified as D. oryzae, neither the genome PL65 nor A5410 harbors the zeamine biosynthesis gene cluster, which plays a key role in virulence of other Dickeya species. The percentages of average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) between the two genomes were 94.47 and 57.00, respectively. In this study, we compared the major virulence factors [plant cell wall-degrading extracellular enzymes and protease (Prt)] produced by D. zeae strains and evaluated the virulence on taro corms and pineapple leaves. Both strains produced Prts, pectate lyases (Pels), and cellulases but no significant quantitative differences were observed (p > 0.05) between the strains. All the strains produced symptoms on taro corms and pineapple leaves, but the strain PL65 produced symptoms more rapidly than others. Our study highlights the genetic constituents of pathogenicity determinants and genomic heterogeneity that will help to understand the virulence mechanisms and aggressiveness of this plant pathogen.

Core and Accessory Genome Analysis of Vibrio mimicus

Vibrio mimicus is an emerging pathogen, mainly associated with contaminated seafood consumption. However, little is known about its evolution, biodiversity, and pathogenic potential. This study analyzes the pan-, core, and accessory genomes of nine V. mimicus strains. The core genome yielded 2424 genes in chromosome I (ChI) and 822 genes in chromosome II (ChII), with an accessory genome comprising an average of 10.9% of the whole genome for ChI and 29% for ChII. Core genome phylogenetic trees were obtained, and V. mimicus ATCC-33654 strain was the closest to the outgroup in both chromosomes. Additionally, a phylogenetic study of eight conserved genes (ftsZ, gapA, gyrB, topA, rpoA, recA, mreB, and pyrH), including Vibrio cholerae, Vibrio parilis, Vibrio metoecus, and Vibrio caribbenthicus, clearly showed clade differentiation. The main virulence genes found in ChI corresponded with type I secretion proteins, extracellular components, flagellar proteins, and potential regulators, while, in ChII, the main categories were type-I secretion proteins, chemotaxis proteins, and antibiotic resistance proteins. The accessory genome was characterized by the presence of mobile elements and toxin encoding genes in both chromosomes. Based on the genome atlas, it was possible to characterize differential regions between strains. The pan-genome of V. mimicus encompassed 3539 genes for ChI and 2355 genes for ChII. These results give us an insight into the virulence and gene content of V. mimicus, as well as constitute the first approach to its diversity.

The identification of novel immunogenic antigens as potential Shigella vaccine components

BACKGROUND

Shigella is a major diarrheal pathogen for which there is presently no vaccine. Whole genome sequencing provides the ability to predict and derive novel antigens for use as vaccines. Here, we aimed to identify novel immunogenic Shigella antigens that could serve as Shigella vaccine candidates, either alone, or when conjugated to Shigella O-antigen.

METHODS

Using a reverse vaccinology approach, where genomic analysis informed the Shigella immunome via an antigen microarray, we aimed to identify novel immunogenic Shigella antigens. A core genome analysis of Shigella species, pathogenic and non-pathogenic Escherichia coli, led to the selection of 234 predicted immunogenic Shigella antigens. These antigens were expressed and probed with acute and convalescent serum from microbiologically confirmed Shigella infections.

RESULTS

Several Shigella antigens displayed IgG and IgA seroconversion, with no difference in sero-reactivity across by sex or age. IgG sero-reactivity to key Shigella antigens was observed at birth, indicating transplacental antibody transfer. Six antigens (FepA, EmrK, FhuA, MdtA, NlpB, and CjrA) were identified in in vivo testing as capable of producing binding IgG and complement-mediated bactericidal antibody.

CONCLUSIONS

These findings provide six novel immunogenic Shigella proteins that could serve as candidate vaccine antigens, species-specific carrier proteins, or targeted adjuvants.

Genomic Insights of Cryobacterium Isolated From Ice Core Reveal Genome Dynamics for Adaptation in Glacier

Glacier is the dominant cold habitat in terrestrial environments, providing a model ecosystem to explore extremophilic strategies and study early lives on Earth. The dominant form of life in glaciers is bacteria. However, little is known about past evolutionary processes that bacteria underwent during adaptation to the cryosphere and the connection of their genomic traits to environmental stressors. Aiming to test the hypothesis that bacterial genomic content and dynamics are driven by glacial environmental stressors, we compared genomes of 21 psychrophilic Cryobacterium strains, including 14 that we isolated from three Tibetan ice cores, to their mesophilic counterparts from the same family Microbacteriaceae of Actinobacteria. The results show that psychrophilic Cryobacterium underwent more dynamic changes in genome content, and their genomes have a significantly higher number of genes involved in stress response, motility, and chemotaxis than their mesophilic counterparts (P < 0.05). The phylogenetic birth-and-death model imposed on the phylogenomic tree indicates a vast surge in recent common ancestor of psychrophilic Cryobacterium (gained the greatest number of genes by 1,168) after the division of the mesophilic strain Cryobacterium mesophilum. The expansion in genome content brought in key genes primarily of the categories “cofactors, vitamins, prosthetic groups, pigments,” “monosaccharides metabolism,” and “membrane transport.” The amino acid substitution rates of psychrophilic Cryobacterium strains are two orders of magnitude lower than those in mesophilic strains. However, no significantly higher number of cold shock genes was found in psychrophilic Cryobacterium strains, indicating that multi-copy is not a key factor for cold adaptation in the family Microbacteriaceae, although cold shock genes are indispensable for psychrophiles. Extensive gene acquisition and low amino acid substitution rate might be the strategies of psychrophilic Cryobacterium to resist low temperature, oligotrophy, and high UV radiation on glaciers. The exploration of genome evolution and survival strategies of psychrophilic Cryobacterium deepens our understanding of bacterial cold adaptation.

Isolation and characterisation of novel phages infecting Lactobacillus plantarum and proposal of a new genus, "Silenusvirus"

Bacteria of Lactobacillus sp. are very useful to humans. However, the biology and genomic diversity of their (bacterio)phage enemies remains understudied. Knowledge on Lactobacillus phage diversity should broaden to develop efficient phage control strategies. To this end, organic waste samples were screened for phages against two wine-related Lactobacillus plantarum strains. Isolates were shotgun sequenced and compared against the phage database and each other by phylogenetics and comparative genomics. The new isolates had only three distant relatives from the database, but displayed a high overall degree of genomic similarity amongst them. The latter allowed for the use of one isolate as a representative to conduct transmission electron microscopy and structural protein sequencing, and to study phage adsorption and growth kinetics. The microscopy and proteomics tests confirmed the observed diversity of the new isolates and supported their classification to the family Siphoviridae and the proposal of the new phage genus "Silenusvirus".

EpitoCore: Mining Conserved Epitope Vaccine Candidates in the Core Proteome of Multiple Bacteria Strains

In reverse vaccinology approaches, complete proteomes of bacteria are submitted to multiple computational prediction steps in order to filter proteins that are possible vaccine candidates. Most available tools perform such analysis only in a single strain, or a very limited number of strains. But the vast amount of genomic data had shown that most bacteria contain pangenomes, i.e., their genomic information contains core, conserved genes, and random accessory genes specific to each strain. Therefore, in reverse vaccinology methods it is of the utmost importance to define core proteins and core epitopes. EpitoCore is a decision-tree pipeline developed to fulfill that need. It provides surfaceome prediction of proteins from related strains, defines core proteins within those, calculate their immunogenicity, predicts epitopes for a given set of MHC alleles defined by the user, and then reports if epitopes are located extracellularly and if they are conserved among the core homologs. Pipeline performance is illustrated by mining peptide vaccine candidates in Mycobacterium avium hominissuis strains. From a total proteome of ~4,800 proteins per strain, EpitoCore predicted 103 highly immunogenic core homologs located at cell surface, many of those related to virulence and drug resistance. Conserved epitopes identified among these homologs allows the users to define sets of peptides with potential to immunize the largest coverage of tested HLA alleles using peptide-based vaccines. Therefore, EpitoCore is able to provide automated identification of conserved epitopes in bacterial pangenomic datasets.

Genomic and Metabolomic Analysis of Antarctic Bacteria Revealed Culture and Elicitation Conditions for the Production of Antimicrobial Compounds

Concern about finding new antibiotics against drug-resistant pathogens is increasing every year. Antarctic bacteria have been proposed as an unexplored source of bioactive metabolites; however, most biosynthetic gene clusters (BGCs) producing secondary metabolites remain silent under common culture conditions. Our work aimed to characterize elicitation conditions for the production of antibacterial secondary metabolites from 34 Antarctic bacterial strains based on MS/MS metabolomics and genome mining approaches. Bacterial strains were cultivated under different nutrient and elicitation conditions, including the addition of lipopolysaccharide (LPS), sodium nitroprusside (SNP), and coculture. Metabolomes were obtained by HPLC-QTOF-MS/MS and analyzed through molecular networking. Antibacterial activity was determined, and seven strains were selected for genome sequencing and analysis. Biosynthesis pathways were activated by all the elicitation treatments, which varies among strains and dependents of culture media. Increased antibacterial activity was observed for a few strains and addition of LPS was related with inhibition of Gram-negative pathogens. Antibiotic BGCs were found for all selected strains and the expressions of putative actinomycin, carotenoids, and bacillibactin were characterized by comparison of genomic and metabolomic data. This work established the use of promising new elicitors for bioprospection of Antarctic bacteria and highlights the importance of new "-omics" comparative approaches for drug discovery.

Genomics of Aerobic Photoheterotrophs in Wheat Phyllosphere Reveals Divergent Evolutionary Patterns of Photosynthetic Genes in Methylobacterium spp

Phyllosphere is a habitat to a variety of viruses, bacteria, fungi, and other microorganisms, which play a fundamental role in maintaining the health of plants and mediating the interaction between plants and ambient environments. A recent addition to this catalogue of microbial diversity was the aerobic anoxygenic phototrophs (AAPs), a group of widespread bacteria that absorb light through bacteriochlorophyll α (BChl a) to produce energy without fixing carbon or producing molecular oxygen. However, culture representatives of AAPs from phyllosphere and their genome information are lacking, limiting our capability to assess their potential ecological roles in this unique niche. In this study, we investigated the presence of AAPs in the phyllosphere of a winter wheat (Triticum aestivum L.) in Denmark by employing bacterial colony based infrared imaging and MALDI-TOF mass spectrometry (MS) techniques. A total of ∼4,480 colonies were screened for the presence of cellular BChl a, resulting in 129 AAP isolates that were further clustered into 21 groups based on MALDI-TOF MS profiling, representatives of which were sequenced using the Illumina NextSeq and Oxford Nanopore MinION platforms. Seventeen draft and four complete genomes of AAPs were assembled belonging in Methylobacterium, Rhizobium, Roseomonas, and a novel Alsobacter. We observed a diverging pattern in the evolutionary rates of photosynthesis genes among the highly homogenous AAP strains of Methylobacterium (Alphaproteobacteria), highlighting an ongoing genomic innovation at the gene cluster level.

Comparative Genomics of Rumen Butyrivibrio spp. Uncovers a Continuum of Polysaccharide-Degrading Capabilities

Feeding a global population of 8 billion people and climate change are the primary challenges facing agriculture today. Ruminant livestock are important food-producing animals, and maximizing their productivity requires an understanding of their digestive systems and the roles played by rumen microbes in plant polysaccharide degradation. Members of the genera Butyrivibrio and Pseudobutyrivibrio are a phylogenetically diverse group of bacteria and are commonly found in the rumen, where they are a substantial source of polysaccharide-degrading enzymes for the depolymerization of lignocellulosic material. Our findings have highlighted the immense enzymatic machinery of Butyrivibrio and Pseudobutyrivibrio species for the degradation of plant fiber, suggesting that these bacteria occupy similar niches but apply different degradation strategies in order to coexist in the competitive rumen environment.

Plant polysaccharide breakdown by microbes in the rumen is fundamental to digestion in ruminant livestock. Bacterial species belonging to the rumen genera Butyrivibrio and Pseudobutyrivibrio are important degraders and utilizers of lignocellulosic plant material. These bacteria degrade polysaccharides and ferment the released monosaccharides to yield short-chain fatty acids that are used by the ruminant for growth and the production of meat, milk, and fiber products. Although rumen Butyrivibrio and Pseudobutyrivibrio species are regarded as common rumen inhabitants, their polysaccharide-degrading and carbohydrate-utilizing enzymes are not well understood. In this study, we analyzed the genomes of 40 Butyrivibrio and 6 Pseudobutyrivibrio strains isolated from the plant-adherent fraction of New Zealand dairy cows to explore the polysaccharide-degrading potential of these important rumen bacteria. Comparative genome analyses combined with phylogenetic analysis of their 16S rRNA genes and short-chain fatty acid production patterns provide insight into the genomic diversity and physiology of these bacteria and divide Butyrivibrio into 3 species clusters. Rumen Butyrivibrio bacteria were found to encode a large and diverse spectrum of degradative carbohydrate-active enzymes (CAZymes) and binding proteins. In total, 4,421 glycoside hydrolases (GHs), 1,283 carbohydrate esterases (CEs), 110 polysaccharide lyases (PLs), 3,605 glycosyltransferases (GTs), and 1,706 carbohydrate-binding protein modules (CBM) with predicted activities involved in the depolymerization and transport of the insoluble plant polysaccharides were identified. Butyrivibrio genomes had similar patterns of CAZyme families but varied greatly in the number of genes within each category in the Carbohydrate-Active Enzymes database (CAZy), suggesting some level of functional redundancy. These results suggest that rumen Butyrivibrio species occupy similar niches but apply different degradation strategies to be able to coexist in the rumen.

IMPORTANCE Feeding a global population of 8 billion people and climate change are the primary challenges facing agriculture today. Ruminant livestock are important food-producing animals, and maximizing their productivity requires an understanding of their digestive systems and the roles played by rumen microbes in plant polysaccharide degradation. Members of the genera Butyrivibrio and Pseudobutyrivibrio are a phylogenetically diverse group of bacteria and are commonly found in the rumen, where they are a substantial source of polysaccharide-degrading enzymes for the depolymerization of lignocellulosic material. Our findings have highlighted the immense enzymatic machinery of Butyrivibrio and Pseudobutyrivibrio species for the degradation of plant fiber, suggesting that these bacteria occupy similar niches but apply different degradation strategies in order to coexist in the competitive rumen environment.

Genome-Wide Analyses Revealed Remarkable Heterogeneity in Pathogenicity Determinants, Antimicrobial Compounds, and CRISPR-Cas Systems of Complex Phytopathogenic Genus Pectobacterium

The Pectobacterium genus comprises pectolytic enterobacteria defined as the causal agents of soft rot, blackleg, and aerial stem rot diseases of potato and economically important crops. In this study, we undertook extensive genome-wide comparative analyses of twelve species that conform the Pectobacterium genus. Bioinformatics approaches outlined a low nucleotide identity of P. parmentieri and P. wasabiae with other species, while P. carotovorum subsp. odoriferum was shown to harbor numerous pseudogenes, which suggests low coding capacity and genomic degradation. The genome atlases allowed for distinguishing distinct DNA structures and highlighted suspicious high transcription zones. The analyses unveiled a noteworthy heterogeneity in the pathogenicity determinants. Specifically, phytotoxins, polysaccharides, iron uptake systems, and the type secretion systems III-V were observed in just some species. Likewise, a comparison of gene clusters encoding antimicrobial compounds put in evidence for high conservation of carotovoricin, whereas a few species possessed the phenazine, carbapenem, and carocins. Moreover, three clustered regularly interspaced short palindromic repeats-Cas (CRISPR-Cas) systems: I-E, I-F, and III-A were identified. Surrounding some CRISPR-Cas regions, different toxin and antitoxin systems were found, which suggests bacterial suicide in the case of an immune system failure. Multiple whole-genome alignments shed light on to the presence of a novel cellobiose phosphotransferase system (PTS) exclusive to P. parmenteri, and an unreported T5SS conserved in almost all species. Several regions that were associated with virulence, microbe antagonism, and adaptive immune systems were predicted within genomic islands, which underscored the essential role that horizontal gene transfer has imparted in the dynamic evolution and speciation of Pectobacterium species. Overall, the results decipher the different strategies that each species has developed to infect their hosts, outcompete for food resources, and defend against bacteriophages. Our investigation provides novel genetic insights that will assist in understanding the pathogenic lifestyle of Pectobacterium, a genus that jeopardizes the agriculture sustainability of important crops worldwide.

Taxogenomics and Systematics of the Genus Pantoea

Members of the genus Pantoea are Gram-negative bacteria isolated from various environments. Taxonomic affiliation based on multilocus sequence analysis (MLSA) is used routinely for inferring accurate phylogeny and identification of bacterial species and genera. Partial sequences of five housekeeping genes (fusA, gyrB, leuS, rpoB, and pyrG) were extracted from 206 draft or complete genomes of Pantoea strains publicly available in databases and analyzed together with the representative sequences of the 25 validly published Pantoea type strains to verify and assess their phylogenetic assignations. Of a total of 159 strains assigned to species level, 11.3% of the non-type strains were incorrectly assigned within suitable Pantoea species. The highest proportion of misidentified strains was recorded in Pantoea vagans, 8 out of 15 (53.3%) inaccurate assignations at the species level. One probable reason for this incorrect classification could be the method previously used for strain identification. Forty-seven (22.8%) genome sequences were from strains identified at the genus level only (Pantoea sp.). A combination of MLSA, average nucleotide identities [ANI and MuMmer-based ANI (ANIm)], tetranucleotide usage pattern (TETRA), and genome-based DNA-DNA hybridization (gDDH) data was used to accurately assign 25 of the 47 strains to validly published Pantoea species, while 17 strains could be assigned as putative novel species within the genus Pantoea. Four genomes designed as Pantoea sp. were identified as Mixta calida. Positive and significant correlation coefficients were computed between MLSA and all the indices derived from whole-genome sequences being proposed for species delimitation. gDDH exhibited the best correlation with MLSA while TETRA was the worst. Accurate species-level identification is key to a better understanding of bacterial diversity and evolution. The MLSA scheme used here could be instrumental to determine the correct taxonomic status of new whole-genome sequenced Pantoea strains, especially non-type strains, before depositing into public databases.

Adaptive Evolution of Sphingobium hydrophobicum C1T in Electronic Waste Contaminated River Sediment

Electronic waste (e-waste) has caused a severe worldwide pollution problem. Despite increasing isolation of degradative microorganisms from e-waste contaminated environments, the mechanisms underlying their adaptive evolution in such habitats remain unclear. Sphingomonads generally have xenobiotic-degrading ability and may play important roles in bioremediation. Sphingobium hydrophobicum C1^T, characterized with superior cell surface hydrophobicity, was recently isolated from e-waste contaminated river sediment. To dissect the mechanisms driving its adaptive evolution, we evaluated its stress resistance, sequenced its genome and performed comparative genomic analysis with 19 other Sphingobium strains. Strain C1^T can feed on several kinds of e-waste-derived xenobiotics, exhibits a great resistance to heavy metals and possesses a high colonization ability. It harbors abundant genes involved in environmental adaptation, some of which are intrinsic prior to experiencing e-waste contamination. The extensive genomic variations between strain C1^T and other Sphingobium strains, numerous C1^T-unique genes, massive mobile elements and frequent genome rearrangements reflect a high genome plasticity. Positive selection, gene duplication, and especially horizontal gene transfer drive the adaptive evolution of strain C1^T. Moreover, presence of type IV secretion systems may allow strain C1^T to be a source of beneficial genes for surrounding microorganisms. This study provides new insights into the adaptive evolution of sphingomonads, and potentially guides bioremediation strategies.

Expanding the Diversity of Myoviridae Phages Infecting Lactobacillus plantarum-A Novel Lineage of Lactobacillus Phages Comprising Five New Members

Lactobacillus plantarum is a bacterium with probiotic properties and promising applications in the food industry and agriculture. So far, bacteriophages of this bacterium have been moderately addressed. We examined the diversity of five new L. plantarum phages via whole genome shotgun sequencing and in silico protein predictions. Moreover, we looked into their phylogeny and their potential genomic similarities to other complete phage genome records through extensive nucleotide and protein comparisons. These analyses revealed a high degree of similarity among the five phages, which extended to the vast majority of predicted virion-associated proteins. Based on these, we selected one of the phages as a representative and performed transmission electron microscopy and structural protein sequencing tests. Overall, the results suggested that the five phages belong to the family Myoviridae, they have a long genome of 137,973-141,344 bp, a G/C content of 36.3-36.6% that is quite distinct from their host's, and surprisingly, 7 to 15 tRNAs. Only an average 41/174 of their predicted genes were assigned a function. The comparative analyses unraveled considerable genetic diversity for the five L. plantarum phages in this study. Hence, the new genus "Semelevirus" was proposed, comprising exclusively of the five phages. This novel lineage of Lactobacillus phages provides further insight into the genetic heterogeneity of phages infecting Lactobacillus sp. The five new Lactobacillus phages have potential value for the development of more robust starters through, for example, the selection of mutants insensitive to phage infections. The five phages could also form part of phage cocktails, which producers would apply in different stages of L. plantarum fermentations in order to create a range of organoleptic outputs.

Complete Genome of Bacillus velezensis CMT-6 and Comparative Genome Analysis Reveals Lipopeptide Diversity

The complete genome sequence of Bacillus velezensis type strain CMT-6 is presented for the first time. A comparative analysis between the genome sequences of CMT-6 with the genome of Bacillus amyloliquefaciens DSM7^T, B. velezensis FZB42, and Bacillus subtilis 168 revealed major differences in the lipopeptide synthesis genes. Of the above, only the CMT-6 strain possessed an integrated synthetase gene for synthesizing surfactin, iturin, and fengycin. However, CMT-6 shared 14, 12, and 10 other lipopeptide-producing genes with FZB42, DSM7^T, and 168 respectively. The largest numbers of non-synonymous mutations were detected in 205 gene sequences that produced these three lipopeptides in CMT-6 and 168. Comparing CMT-6 with DSM7^T, 58 non-synonymous mutations were detected in gene sequences that contributed to produce lipopeptides. In addition, InDels were identified in yczE and glnR genes. CMT-6 and FZB42 had the lowest number of non-synonymous mutations with 8 lipopeptide-related gene sequences. And InDels were identified in only yczE. The numbers of core genes, InDels, and non-synonymous mutations in genes were the main reasons for the differences in yield and variety of lipopeptides. These results will enrich the genomic resources available for B. velezensis and provide fundamental information to construct strains that can produce specific lipopeptides.

Iron-associated protein interaction networks reveal the key functional modules related to survival and virulence of Pasteurella multocida

Pasteurella multocida causes respiratory infectious diseases in a multitude of birds and mammals. A number of virulence-associated genes were reported across different strains of P. multocida, including those involved in the iron transport and metabolism. Comparative iron-associated genes of P. multocida among different animal hosts towards their interaction networks have not been fully revealed. Therefore, this study aimed to identify the iron-associated genes from core- and pan-genomes of fourteen P. multocida strains and to construct iron-associated protein interaction networks using genome-scale network analysis which might be associated with the virulence. Results showed that these fourteen strains had 1587 genes in the core-genome and 3400 genes constituting their pan-genome. Out of these, 2651 genes associated with iron transport and metabolism were selected to construct the protein interaction networks and 361 genes were incorporated into the iron-associated protein interaction network (iPIN) consisting of nine different iron-associated functional modules. After comparing with the virulence factor database (VFDB), 21 virulence-associated proteins were determined and 11 of these belonged to the heme biosynthesis module. From this study, the core heme biosynthesis module and the core outer membrane hemoglobin receptor HgbA were proposed as candidate targets to design novel antibiotics and vaccines for preventing pasteurellosis across the serotypes or animal hosts for enhanced precision agriculture to ensure sustainability in food security.

Contrasted evolutionary constraints on carbohydrate active enzymes (CAZymes) in selected Frankia strains

Carbohydrate active enzymes (CAZymes) are capable of breaking complex polysaccharides into simpler form. In plant-host-associated microorganisms CAZymes are known to be involved in plant cell wall degradation. However, the biology and evolution of Frankia CAZymes are largely unknown. In the present study, we took a genomic approach to evaluate the presence and putative roles of CAZymes in Frankia. The CAZymes were found to be potentially highly expressed (PHX) proteins and contained more aromatic amino acids, which increased their biosynthetic energy cost. These energy rich amino acids were present in the active sites of CAZymes aiding in their carbohydrate binding capacity. Phylogenetic and evolutionary analyses showed that, in Frankia strains with the capacity to nodulate host plants, CAZymes were evolving slower than the other PHX genes, whereas similar genes from non-nodulating (or ineffectively nodulating) Frankia strains showed little variation in their evolutionary constraints compared to other PHX genes. Thus, the present study revealed the persistence of a strong purifying selection on CAZymes of Frankia indicating their crucial role.

Comparative Transcriptome Profiling Reveals a Potential Role of Type VI Secretion System and Fimbriae in Virulence of Non-O157 Shiga Toxin-Producing Escherichia coli

Shiga toxin-producing Escherichia coli (STEC) cause both sporadic infections and outbreaks of enteric disease in humans, with symptoms ranging from asymptomatic carriage to severe disease like haemolytic uremic syndrome (HUS). Bacterial virulence factors like subtypes of the Shiga toxin (Stx) and the locus of enterocyte effacement (LEE) pathogenicity island, as well as host factors like young age, are strongly associated with development of HUS. However, these factors alone do not accurately differentiate between strains that cause HUS and those that do not cause severe disease, which is important in the context of diagnosis, treatment, as well as infection control. We have used RNA sequencing to compare transcriptomes of 30 stx2a and eae positive STEC strains of non-O157 serogroups isolated from children <5 years of age. The strains were from children with HUS (HUS group, n = 15), and children with asymptomatic or mild disease (non-HUS group, n = 15), either induced with mitomycin C or non-induced, to reveal potential differences in gene expression levels between groups. When the HUS and non-HUS group were compared for differential expression of protein-encoding gene families, 399 of 6,119 gene families were differentially expressed (log2 fold change ≥ 1, FDR < 0.05) in the non-induced condition, whereas only one gene family was differentially expressed in the induced condition. Gene ontology and cluster analysis showed that several fimbrial operons, as well as a putative type VI secretion system (T6SS) were more highly expressed in the HUS group than in the non-HUS group, indicating a role of these in the virulence of STEC strains causing severe disease.

Integrative and Conjugative Elements (ICEs) in Pasteurellaceae Species and Their Detection by Multiplex PCR

Strains of the Pasteurellaceae bacteria Pasteurella multocida and Mannheimia haemolytica are major etiological agents of bovine respiratory disease (BRD). Treatment of BRD with antimicrobials is becoming more challenging due to the increasing occurrence of resistance in infecting strains. In Pasteurellaceae strains exhibiting resistance to multiple antimicrobials including aminoglycosides, beta-lactams, macrolides and sulfonamides, the resistance determinants are often chromosomally encoded within integrative and conjugative elements (ICEs). To gain a more comprehensive picture of ICE structures, we sequenced the genomes of six strains of P. multocida and four strains of M. haemolytica; all strains were independent isolates and eight of them were multiple-resistant. ICE sequences varied in size from 49 to 79 kb, and were comprised of an array of conserved genes within a core region and varieties of resistance genes within accessory regions. These latter regions mainly account for the variation in the overall ICE sizes. From the sequence data, we developed a multiplex PCR assay targeting four conserved core genes required for integration and maintenance of ICE structures. Application of this assay on 75 isolates of P. multocida and M. haemolytica reveals how the presence and structures of ICEs are related to their antibiotic resistance phenotypes. The assay is also applicable to other members of the Pasteurellaceae family including Histophilus somni and indicates how clustering and dissemination of the resistance genes came about.

Subtractive proteomics to identify novel drug targets and reverse vaccinology for the development of chimeric vaccine against Acinetobacter baumannii

The emergence of drug-resistant Acinetobacter baumannii is the global health problem associated with high mortality and morbidity. Therefore it is high time to find a suitable therapeutics for this pathogen. In the present study, subtractive proteomics along with reverse vaccinology approaches were used to predict suitable therapeutics against A. baumannii. Using subtractive proteomics, we have identified promiscuous antigenic membrane proteins that contain the virulence factors, resistance factors and essentiality factor for this pathogenic bacteria. Selected promiscuous targeted membrane proteins were used for the design of chimeric-subunit vaccine with the help of reverse vaccinology. Available best tools and servers were used for the identification of MHC class I, II and B cell epitopes. All selected epitopes were further shortlisted computationally to know their immunogenicity, antigenicity, allergenicity, conservancy and toxicity potentials. Immunogenic predicted promiscuous peptides used for the development of chimeric subunit vaccine with immune-modulating adjuvants, linkers, and PADRE (Pan HLA-DR epitopes) amino acid sequence. Designed vaccine construct V4 also interact with the MHC, and TLR4/MD2 complex as confirm by docking and molecular dynamics simulation studies. Therefore designed vaccine construct V4 can be developed to control the host-pathogen interaction or infection caused by A. baumannii.

Genome analysis of two novel Pseudomonas strains exhibiting differential hypersensitivity reactions on tobacco seedlings reveals differences in nonflagellar T3SS organization and predicted effector proteins

Multilocus sequence analysis (MLSA) of two new biological control strains (S1E40 and S3E12) of Pseudomonas was performed to assess their taxonomic position relative to close lineages, and comparative genomics employed to investigate whether these strains differ in key genetic features involved in hypersensitivity responses (HRs). Strain S3E12, at high concentration, incites HRs on tobacco and corn plantlets while S1E40 does not. Phylogenies based on individual genes and 16S rRNA-gyrB-rpoB-rpoD concatenated sequence data show strains S1E40 and S3E12 clustering in distinct groups. Strain S3E12 consistently clustered with Pseudomonas marginalis, a bacterium causing soft rots on plant tissues. MLSA data suggest that strains S1E40 and S3E12 are novel genotypes. This is consistent with the data of genome-based DNA-DNA homology values that are below the proposed cutoff species boundary. Comparative genomics analysis of the two strains revealed major differences in the type III secretion systems (T3SS) as well as the predicted T3SS secreted effector proteins (T3Es). One nonflagellar (NF-T3SS) and two flagellar T3SSs (F-T3SS) clusters were identified in both strains. While F-T3SS clusters in both strains were relatively conserved, the NF-T3SS clusters differed in the number of core components present. The predicted T3Es also differed in the type and number of CDSs with both strains having unique predicted protease-related effectors. In addition, the T1SS organization of the S3E12 genome has protein-coding sequences (CDSs) encoding for key factors such as T1SS secreted agglutinin repeats-toxins (a group of cytolysins and cytotoxins), a membrane fusion protein (LapC), a T1SS ATPase of LssB family (LapB), and T1SS-associated transglutaminase-like cysteine proteinase (LapP). In contrast, strain S1E40 has all CDSs for the seven-gene operon (pelA-pelG) required for Pel biosynthesis but not S3E12, suggesting that biofilm formation in these strains is modulated differently. The data presented here provide an insight of the genome organization of these two phytobacterial strains.

Comparative genomics of Mycobacterium reveals evolutionary trends of M. avium complex

Mycobacterium is gram positive, slow growing, disease causing Actinobacteria. Beside potential pathogenic species, Mycobacterium also contains opportunistic pathogens as well as free living non-pathogenic species. Disease related various analyses on Mycobacterium tuberculosis are very widespread. However, genomic study of overall Mycobacterium species for understanding the selection pressure on genes as well as evolution of the organism is still illusive. MLSA and 16s rDNA based analysis has been generated for 241 Mycobacterium strains and a detailed analysis of codon and amino acid usage bias of mycobacterial genes, their functional analysis have been done. Further the evolutionary features of M. avium complex also have been revealed. Mycobacterial genes are moderately GC rich showed higher expression level in PPs and significant negative correlation with biosynthetic cost of proteins. Translational selection pressure was observed in mycobacterial genes. MAC showed close relationship with NPs and higher evolutionary rate in MAC revealed their constant evolving nature.

Genomic Diversity in the Endosymbiotic Bacterium Rhizobium leguminosarum

Rhizobium leguminosarum bv. viciae is a soil α-proteobacterium that establishes a diazotrophic symbiosis with different legumes of the Fabeae tribe. The number of genome sequences from rhizobial strains available in public databases is constantly increasing, although complete, fully annotated genome structures from rhizobial genomes are scarce. In this work, we report and analyse the complete genome of R. leguminosarum bv. viciae UPM791. Whole genome sequencing can provide new insights into the genetic features contributing to symbiotically relevant processes such as bacterial adaptation to the rhizosphere, mechanisms for efficient competition with other bacteria, and the ability to establish a complex signalling dialogue with legumes, to enter the root without triggering plant defenses, and, ultimately, to fix nitrogen within the host. Comparison of the complete genome sequences of two strains of R. leguminosarum bv. viciae, 3841 and UPM791, highlights the existence of different symbiotic plasmids and a common core chromosome. Specific genomic traits, such as plasmid content or a distinctive regulation, define differential physiological capabilities of these endosymbionts. Among them, strain UPM791 presents unique adaptations for recycling the hydrogen generated in the nitrogen fixation process.

Genome sequence of "Leucobacter massiliensis" sp. nov. isolated from human pharynx after travel to the 2014 Hajj

"Leucobacter massiliensis" strain 122RC15T sp. nov. is a new species within the genus Leucobacter. The genome of this strain is described here. It was isolated from the pharynx of a 76-year-old Algerian female after travelling from the 2014 Hajj. "Leucobacter massiliensis" is a Gram-positive, aerobic bacillus. Here we describe the features including complete genome and annotation of this strain. The 3 136 406-bp long genome contains 2797 protein-coding genes and 49 RNA genes.

The diversity and specificity of the extracellular proteome in the cellulolytic bacterium Caldicellulosiruptor bescii is driven by the nature of the cellulosic growth substrate

BACKGROUND

Caldicellulosiruptor bescii is a thermophilic cellulolytic bacterium that efficiently deconstructs lignocellulosic biomass into sugars, which subsequently can be fermented into alcohols, such as ethanol, and other products. Deconstruction of complex substrates by C. bescii involves a myriad of highly abundant, substrate-specific extracellular solute binding proteins (ESBPs) and carbohydrate-active enzymes (CAZymes) containing carbohydrate-binding modules (CBMs). Mass spectrometry-based proteomics was employed to investigate how these substrate recognition proteins and enzymes vary as a function of lignocellulosic substrates.

RESULTS

Proteomic analysis revealed several key extracellular proteins that respond specifically to either C5 or C6 mono- and polysaccharides. These include proteins of unknown functions (PUFs), ESBPs, and CAZymes. ESBPs that were previously shown to interact more efficiently with hemicellulose and pectin were detected in high abundance during growth on complex C5 substrates, such as switchgrass and xylan. Some proteins, such as Athe_0614 and Athe_2368, whose functions are not well defined were predicted to be involved in xylan utilization and ABC transport and were significantly more abundant in complex and C5 substrates, respectively. The proteins encoded by the entire glucan degradation locus (GDL; Athe_1857, 1859, 1860, 1865, 1867, and 1866) were highly abundant under all growth conditions, particularly when C. bescii was grown on cellobiose, switchgrass, or xylan. In contrast, the glycoside hydrolases Athe_0609 (Pullulanase) and 0610, which both possess CBM20 and a starch binding domain, appear preferential to C5/complex substrate deconstruction. Some PUFs, such as Athe_2463 and 2464, were detected as highly abundant when grown on C5 substrates (xylan and xylose), also suggesting C5-substrate specificity.

CONCLUSIONS

This study reveals the protein membership of the C. bescii secretome and demonstrates its plasticity based on the complexity (mono-/disaccharides vs. polysaccharides) and type of carbon (C5 vs. C6) available to the microorganism. The presence or increased abundance of extracellular proteins as a response to specific substrates helps to further elucidate C. bescii's utilization and conversion of lignocellulosic biomass to biofuel and other valuable products. This includes improved characterization of extracellular proteins that lack discrete functional roles and are poorly/not annotated.

Genomic analysis of the marine fish pathogen Photobacterium damselae subsp. piscicida: Insertion sequences proliferation is associated with chromosomal reorganisations and rampant gene decay

Photobacterium damselae subsp. piscicida (Pdp) is an intracellular fish pathogen that causes photobacteriosis, a disease proven deadly in farmed fish worldwide. This work focuses on the analysis of genome sequences, chromosomes structure and gene contents of two strains from Sparus aurata (DI21) and Solea senegalensis (L091106-03H), isolated on the Spanish Atlantic coast. The comparative genomic analysis revealed that DI21 and L091106-03H share 98% of their genomes, including two virulence plasmids: pPHDP70 encoding siderophore piscibactin synthesis and pPHDP10 encoding the apoptotic toxin AIP56. Both genomes harbour a surprisingly large number of IS elements accounting for 12-17% of the total genome, representing an IS density of 0.15 elements per kb, one of the highest IS density values in a bacterial pathogen. This massive proliferation of ISs is responsible for the generation of a high number of pseudogenes that caused extensive loss of biological functions. Pseudogene formation is one of the main features of Pdp genome that explains most of the ecological and phenotypic differences with respect to its sibling subspecies P. damselae subsp. damselae and to other Vibrionaceae. Evidence was also found proving the existence of two chromosomal configurations depending on the origin of the strains: an European and an Asian/American types of genome organisation, reinforcing the idea of the existence of two geographically-linked clonal lineages in Pdp. In short, our study suggests that the host-dependent lifestyle of Pdp allowed massive IS proliferation and gene decay processes, which are major evolutionary forces in the shaping of the Pdp genome.