riboSeed: leveraging prokaryotic genomic architecture to assemble across ribosomal regions

rRNA operon improves species-level classification of bacteria and microbial community analysis compared to 16S rRNA

Precise identification of species is fundamental in microbial genomics and is crucial for understanding the microbial communities. While the 16S rRNA gene, particularly its V3-V4 regions, has been extensively employed for microbial identification, however has limitations in achieving species-level resolution. Advancements in long-read sequencing technologies have highlighted the rRNA operon as a more accurate marker for microbial classification and analysis than the 16S rRNA gene. This study aims to compare the accuracy of species classification and microbial community analysis using the rRNA operon versus the 16S rRNA gene. We evaluated the species classification accuracy of the rRNA operon,16S rRNA gene, and 16S rRNA V3-V4 regions using a BLAST-based method and a k-mer matching-based method with public data available from NCBI. We further performed simulations to model microbial community analysis. We accessed the performance using each marker in community composition estimation and differential abundance analysis. Our findings demonstrate that the rRNA operon offers an advantage over the 16S rRNA gene and its V3-V4 regions for species-level classification within the genus. When applied to microbial community analysis, the rRNA operon enables a more accurate determination of composition. Using the rRNA operon yielded more reliable results in differential abundance analysis as well.

IMPORTANCE

We quantitatively demonstrated that the rRNA operon outperformed the 16S rRNA and its V3-V4 regions in accuracy for both individual species identification and species-level microbial community analysis. Our findings can provide guidelines for selecting appropriate markers in the field of microbial research.

Comparative Genome Analysis of Enterococcus cecorum Reveals Intercontinental Spread of a Lineage of Clinical Poultry Isolates

Enterococcus cecorum is an emerging pathogen responsible for osteomyelitis, spondylitis, and femoral head necrosis causing animal suffering and mortality and requiring antimicrobial use in poultry. Paradoxically, E. cecorum is a common inhabitant of the intestinal microbiota of adult chickens. Despite evidence suggesting the existence of clones with pathogenic potential, the genetic and phenotypic relatedness of disease-associated isolates remains little investigated. Here, we sequenced and analyzed the genomes and characterized the phenotypes of more than 100 isolates, the majority of which were collected over the last 10 years from 16 French broiler farms. Comparative genomics, genome-wide association studies, and the measured susceptibility to serum, biofilm-forming capacity, and adhesion to chicken type II collagen were used to identify features associated with clinical isolates. We found that none of the tested phenotypes could discriminate the origin of the isolates or the phylogenetic group. Instead, we found that most clinical isolates are grouped phylogenetically, and our analyses selected six genes that discriminate 94% of isolates associated with disease from those that are not. Analysis of the resistome and the mobilome revealed that multidrug-resistant clones of E. cecorum cluster into a few clades and that integrative conjugative elements and genomic islands are the main carriers of antimicrobial resistance. This comprehensive genomic analysis shows that disease-associated clones of E. cecorum belong mainly to one phylogenetic clade. IMPORTANCE Enterococcus cecorum is an important pathogen of poultry worldwide. It causes a number of locomotor disorders and septicemia, particularly in fast-growing broilers. Animal suffering, antimicrobial use, and associated economic losses require a better understanding of disease-associated E. cecorum isolates. To address this need, we performed whole-genome sequencing and analysis of a large collection of isolates responsible for outbreaks in France. By providing the first data set on the genetic diversity and resistome of E. cecorum strains circulating in France, we pinpoint an epidemic lineage that is probably also circulating elsewhere that should be targeted preferentially by preventive strategies in order to reduce the burden of E. cecorum-related diseases.

Comprehensive investigation of antibiotic resistance gene content in cfiA-harboring Bacteroides fragilis isolates of human and animal origins by whole genome sequencing

INTRODUCTION

The emergence of multidrug resistance in Bacteroides fragilis, especially the phylogenetic lineage carrying the carbapenemase gene cfiA, represents an increasing threat to human health. However, knowledge on the diversity of the multidrug-resistant strains and the genetic elements carrying the antibiotic resistance genes (ARGs) remains limited.

AIM

The objective of the study was to describe the resistome in cfiA-positive B. fragilis.

METHODS

A collection of cfiA-positive B. fragilis from diverse human (8 bacteremias, 15 wound infections) and animal (2 chickens, 2 pigs, 6 dogs, 3 cats) sources in Hong Kong, 2015-2017 was analysed by whole genome sequencing.

RESULTS

In the 36 isolates, 13 distinct ARGs (total number 83, median 2, range 0-7 per isolate) other than cfiA were detected. ARGs encoding resistance to aminoglycosides, β-lactams, macrolides, sulphonamides and tetracyclines were carried by CTn341-like, CTnHyb-like, Tn5220-like, Tn4555-like and Tn613-like transposons and were detected in phylogenetically diverse isolates of different host sources. Only few ARGs encoding resistance to metronidazole and tetracyclines were localized on plasmids. In two chicken isolates, a novel transposon (designated as Tn6994) was found to be involved in the dissemination of multiple ARGs mediating resistance to multiple antibiotics, including metronidazole and linezolid that are critically important for treatment of anaerobic infections. In mating experiments, Tn6994 and the associated phenotypic resistance could be transferred to Bacteroides nordii recipient.

CONCLUSION

This study illustrates the importance of transposons in the dissemination of ARGs in the cfiA-positive division of B. fragilis. One Health approach is necessary to track the dissemination of ARGs.

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.

Ecological Divergence Within the Enterobacterial Genus Sodalis: From Insect Symbionts to Inhabitants of Decomposing Deadwood

The bacterial genus Sodalis is represented by insect endosymbionts as well as free-living species. While the former have been studied frequently, the distribution of the latter is not yet clear. Here, we present a description of a free-living strain, Sodalis ligni sp. nov., originating from decomposing deadwood. The favored occurrence of S. ligni in deadwood is confirmed by both 16S rRNA gene distribution and metagenome data. Pangenome analysis of available Sodalis genomes shows at least three groups within the Sodalis genus: deadwood-associated strains, tsetse fly endosymbionts and endosymbionts of other insects. This differentiation is consistent in terms of the gene frequency level, genome similarity and carbohydrate-active enzyme composition of the genomes. Deadwood-associated strains contain genes for active decomposition of biopolymers of plant and fungal origin and can utilize more diverse carbon sources than their symbiotic relatives. Deadwood-associated strains, but not other Sodalis strains, have the genetic potential to fix N2, and the corresponding genes are expressed in deadwood. Nitrogenase genes are located within the genomes of Sodalis, including S. ligni, at multiple loci represented by more gene variants. We show decomposing wood to be a previously undescribed habitat of the genus Sodalis that appears to show striking ecological divergence.

Ornithobacterium rhinotracheale: MALDI-TOF MS and Whole Genome Sequencing Confirm That Serotypes K, L and M Deviate from Well-Known Reference Strains and Numerous Field Isolates

Ornithobacterium rhinotracheale is one of the most important bacterial agents of respiratory diseases in poultry. For correct identification and characterization of this fastidious bacterium, reliable diagnostic tools are essential. Still, phenotypic tests are used to identify O. rhinotracheale and serotyping is the most common method for characterization, despite known drawbacks and disadvantages such as divergent results, cross-reactivity between strains, or the non-typeability of strains. The intention of the present study was to evaluate MALDI-TOF MS and whole genome sequencing for the identification and characterization of O. rhinotracheale. For this purpose, a selection of 59 well-defined reference strains and 47 field strains derived from outbreaks on Austrian turkey farms were investigated by MALDI-TOF MS. The field strains originated from different geographical areas in Austria with some of the isolates derived from multiple outbreaks on farms within a year, or recurrent outbreaks over several years. MALDI-TOF MS proved a suitable method for identification of O. rhinotracheale to genus or species level except for 3 strains representing serotypes M, K and F. Phylogenetic analysis showed that most strains grouped within one cluster even though they were comprised of different serotypes, while serotypes F, K, and M clearly formed a different cluster. All field isolates from turkey farms clustered together, independent of the origin of the isolates, e.g., geographical area, multiple outbreaks within a year or recurrent outbreaks over several years. Whole genome sequencing of serotype M, K and F strains confirmed the extraordinary status and deviation from known fully-sequenced strains due to a lack of sequence similarity. This was further confirmed by alignments of single genes (16S-RNA and rpoB) and multilocus sequence typing although the demarcation was less obvious. Altogether, the results indicate that these three serotypes belong to a different species than O. rhinotracheale, and might even be members of multiple new species.

Uncovering patterns of the evolution of genomic sequence entropy and complexity

The lack of consensus concerning the biological meaning of entropy and complexity of genomes and the different ways to assess these data hamper conclusions concerning what are the causes of genomic entropy variation among species. This study aims to evaluate the entropy and complexity of genomic sequences of several species without using homologies to assess relationships among these variables and non-molecular data (e.g., the number of individuals) to seek a trigger of interspecific genomic entropy variation. The results indicate a relationship among genomic entropy, genome size, genomic complexity, and the number of individuals: species with a small number of individuals harbors large genome, and hence, low entropy but a higher complexity. We defined that the complexity of a genome relies on the entropy of each DNA segment within genome. Then, the entropy and complexity of a genome reflects its organization solely. Exons of vertebrates harbor smaller entropies than non-exon regions (likely by the repeats that accumulated from duplications), whereas other taxonomic groups do not present this pattern. Our findings suggest that small initial population might have defined current genomic entropy and complexity: actual genomes are less complex than ancestral ones. Besides, our data disagree with the relationship between phenotype and genomic entropies previously established. Finally, by establishing the relationship between genomic entropy/complexity with the number of individuals and genome size, under an evolutive perspective, ideas concerning the genomic variability may emerge.

Re-sequencing and optical mapping reveals misassemblies and real inversions on Corynebacterium pseudotuberculosis genomes

The number of draft genomes deposited in Genbank from the National Center for Biotechnology Information (NCBI) is higher than the complete ones. Draft genomes are assemblies that contain fragments of misassembled regions (gaps). Such draft genomes present a hindrance to the complete understanding of the biology and evolution of the organism since they lack genomic information. To overcome this problem, strategies to improve the assembly process are developed continuously. Also, the greatest challenge to the assembly progress is the presence of repetitive DNA regions. This article highlights the use of optical mapping, to detect and correct assembly errors in Corynebacterium pseudotuberculosis. We also demonstrate that choosing a reference genome should be done with caution to avoid assembly errors and loss of genetic information.