Optimized Illumina PCR-free library preparation for bacterial whole genome sequencing and analysis of factors influencing de novo assembly

Acinetobacter corruptisaponis sp. nov., Isolated from a Spoiled Bath Lotion

Strain DM2021935T representing a novel Acinetobacter species was isolated from a spoiled bath lotion in Guangdong, China. Based on 16S rRNA gene phylogenetic analysis, strain DM2021935T was closely related to 'Acinetobacter thutiue' VNH17T, Acinetobacter junii CIP 64.5 T, and Acinetobacter tibetensis Y-23 T. Cells of strain DM2021935T were Gram-stain-negative, non-spore-forming, strictly aerobic, catalase-positive, oxidase-negative, α-hemolytic, and non-motile. Strain DM2021935T exhibited growth in 1-3% (w/v) NaCl at temperatures ranging from 4 to 37 °C and tolerated pH levels from 6.0 to 8.0. The predominant fatty acids in strain DM2021935T are C12:0, C16:0, C18:1 ω9c, and summed feature 3. Polar lipid profiles included glycolipids, phospholipids, phosphatidylethanolamine, and phosphatidyl-N-methylethanolamine. The identified respiratory quinones were ubiquinone Q-8 and Q-9. The genomic size of DM2021935T comprised 4.15 Mb, consisting of one chromosome (3,827,633 bp) and two plasmids (241,357 and 83,010 bp). The G + C content was 41.8%. The average nucleotide identity, average amino acid identity, and digital DNA-DNA hybridization values between strain DM2021935T and phylogenetically related type strains were below the species delineation thresholds (72.2-95.4, 53.1-87.0, and 20.4-66.4%, respectively). AntiSMASH analysis identified four gene clusters: non-ribosomal peptide synthetase, non-alpha poly-amino group acids, YcaO cyclodehydratase, and aryl polyene biosynthesis. Based on genotypic data, strain DM2021935T represents a novel species within the genus Acinetobacter. The proposed name for the novel species is Acinetobacter corruptisaponis sp. nov. (type strain DM2021935T = KCTC 92772 T = GDMCC 1.3703 T).

An overlooked phenomenon: complex interactions of potential error sources on the quality of bacterial de novo genome assemblies

BACKGROUND

Parameters adversely affecting the contiguity and accuracy of the assemblies from Illumina next-generation sequencing (NGS) are well described. However, past studies generally focused on their additive effects, overlooking their potential interactions possibly exacerbating one another's effects in a multiplicative manner. To investigate whether or not they act interactively on de novo genome assembly quality, we simulated sequencing data for 13 bacterial reference genomes, with varying levels of error rate, sequencing depth, PCR and optical duplicate ratios.

RESULTS

We assessed the quality of assemblies from the simulated sequencing data with a number of contiguity and accuracy metrics, which we used to quantify both additive and multiplicative effects of the four parameters. We found that the tested parameters are engaged in complex interactions, exerting multiplicative, rather than additive, effects on assembly quality. Also, the ratio of non-repeated regions and GC% of the original genomes can shape how the four parameters affect assembly quality.

CONCLUSIONS

We provide a framework for consideration in future studies using de novo genome assembly of bacterial genomes, e.g. in choosing the optimal sequencing depth, balancing between its positive effect on contiguity and negative effect on accuracy due to its interaction with error rate. Furthermore, the properties of the genomes to be sequenced also should be taken into account, as they might influence the effects of error sources themselves.

Clinical application of amplification-based versus amplification-free metagenomic next-generation sequencing test in infectious diseases

Background

Recently, metagenomic next-generation sequencing (mNGS) has been used in the diagnosis of infectious diseases (IDs) as an emerging and powerful tool. However, whether the complicated methodological variation in mNGS detections makes a difference in their clinical performance is still unknown. Here we conducted a method study on the clinical application of mNGS tests in the DNA detection of IDs.

Methods

We analyzed the effect of several potential factors in the whole process of mNGS for DNA detection on microorganism identification in 98 samples of suspected ID patients by amplification-based mNGS. The amplification-based and amplification-free mNGS tests were successfully performed in 41 samples. Then we compared the clinical application of the two mNGS methods in the DNA detection of IDs.

Results

We found that a higher concentration of extracted nucleic acid was more conducive to detecting microorganisms. Other potential factors, such as read depth and proportion of human reads, might not be attributed to microorganism identification. The concordance rate of amplification-based and amplification-free mNGS results was 80.5% (33/41) in the patients with suspected IDs. Amplification-based mNGS showed approximately 16.7% higher sensitivity than amplification-free mNGS. However, 4 cases with causative pathogens only detected by amplification-based mNGS were finally proved false-positive. In addition, empirical antibiotic treatments were adjusted in 18 patients following mNGS testing with unexpected pathogens.

Conclusions

Amplification-based and amplification-free mNGS tests showed their specific advantages and disadvantages in the diagnosis of IDs. The clinical application of mNGS still needs more exploration from a methodological perspective. With advanced technology and standardized procedure, mNGS will play a promising role in the diagnosis of IDs and help guide the use of antibiotics.

Fundicoccus culcitae sp. nov., a novel potential bacteriocin producing bacterium isolated from a spoiled eye mask

A Gram-positive, facultatively anaerobic, oval beaded-shape, oxidase-negative, and non-motile bacterium designated DM20194951T was isolated from a spoiled eye mask obtained from Guangdong, China. Based on the 16S rRNA gene sequence, phylogenetic analysis indicated that strain DM20194951T showed the highest sequence similarity (95.8%) to Fundicoccus ignavus WS4937T. Meanwhile, strain DM20194951T could be distinguished from the type strains in the genus Fundicoccus by distinct phenotypic and genotypic traits. Strain DM20194951T grew variably with 1-2% (w/v) NaCl and tolerated pH 6.0-10.0. Growth was observed from 28 to 37 °C. The diagnostic diamino acids in the cell-wall peptidoglycan consisted of aspartic and glutamic acids as well as alanine. The predominant fatty acids were C18:1 ω9c, C16:0, and C16:1 ω9c. In the polar lipid profile, two glycolipids, three phospholipids, one phosphatidylglycerol, and one diphosphatidylglycerol were found. No respiratory quinones were detected. The DM20194951T genome is 3.2 Mb in size and contains a G + C content of 38.1%. A gene cluster for lactococcin 972 family bacteriocin production was found in the DM20194951T genome. Based on morphological, genotypic, and phylogenetic data, strain DM20194951T should be considered to represent a novel species in the genus Fundicoccus, for which the name Fundicoccus culcitae sp. nov. is proposed with the type strain DM20194951T (= KCTC 43472T = GDMCC 1.3614T).

FASTdRNA: a workflow for the analysis of ONT direct RNA sequencing

Motivation

Direct RNA-seq (dRNA-seq) using Oxford Nanopore Technology (ONT) has revolutionized transcript mapping by offering enhanced precision due to its long-read length. Unlike traditional techniques, dRNA-seq eliminates the need for PCR amplification, reducing the impact of GC bias, and preserving valuable base physical information, such as RNA modification and poly(A) length estimation. However, the rapid advancement of ONT devices has set higher standards for analytical software, resulting in potential challenges of software incompatibility and reduced efficiency.

Results

We present a novel workflow, called FASTdRNA, to manipulate dRNA-seq data efficiently. This workflow comprises two modules: a data preprocessing module and a data analysis module. The preprocessing data module, dRNAmain, encompasses basecalling, mapping, and transcript counting, which are essential for subsequent analyses. The data analysis module consists of a range of downstream analyses that facilitate the estimation of poly(A) length, prediction of RNA modifications, and assessment of alternative splicing events across different conditions with duplication. The FASTdRNA workflow is designed for the Snakemake framework and can be efficiently executed locally or in the cloud. Comparative experiments have demonstrated its superior performance compared to previous methods. This innovative workflow enhances the research capabilities of dRNA-seq data analysis pipelines by optimizing existing processes and expanding the scope of analysis.

Availability and implementation

The workflow is freely available at https://github.com/Tomcxf/FASTdRNA under an MIT license. Detailed install and usage guidance can be found in the GitHub repository.

Bombella pluederhausensis sp. nov., Bombella pollinis sp. nov., Bombella saccharophila sp. nov. and Bombella dulcis sp. nov., four Bombella species isolated from the environment of the western honey bee Apis mellifera

Four strains of members of the genus Bombella were isolated from samples associated with the western honey bee Apis mellifera, which could not be assigned to a species with a validly published name. Strains TMW 2.2543^T, TMW 2.2556^T, TMW 2.2558^T and TMW 2.2559^T exhibit in silico DNA-DNA hybridisation (isDDH) and orthologous average nucleotide identity (orthoANI) values below species delineation thresholds compared with all described species of the genus Bombella and with each other. TMW 2.2556^T and TMW 2.2558^T form their own clade within the genus. The major respiratory quinone of all strains was Q-10. The composition of cellular fatty acids was diverse between strains. All strains stained Gram-negative, were rod-shaped, strictly aerobic, pellicle-forming, catalase-positive, oxidase-negative, mesophilic and grew over a wide pH range; they were halosensitive but glucose-tolerant. Unlike the other studied strains, TMW 2.2558^T was non-motile. Phylogenetic, chemotaxonomic and physiological analyses revealed a clear distinction between all the strains and species with validly published names. All the data support the proposition of four novel species within the genus Bombella, namely Bombella pluederhausensis sp. nov., Bombella pollinis sp. nov., Bombella saccharophila sp. nov. and Bombella dulcis sp. nov., with the respective type strains Bombella pluederhausensis sp. nov. TMW 2.2543^T (= DSM 114872^T, = LMG 32791^T), Bombella pollinis sp. nov. TMW 2.2556^T (= DSM 114874^T, = LMG 32792^T), Bombella saccharophila sp. nov. TMW 2.2558^T (= DSM 114875^T, = LMG 32793^T) and Bombella dulcis sp. nov. TMW 2.2559^T (= DSM 114877^T, = LMG 32794^T). Moreover, three genomes available in the NCBI database that have not yet been described as species with validly published names could be assigned to the proposed species. Bombella sp. ESL0378 and Bombella sp. ESL0385 to Bombella pollinis sp. nov. and Bombella sp. AS1 to Bombella saccharophila sp. nov.

DREAMS: deep read-level error model for sequencing data applied to low-frequency variant calling and circulating tumor DNA detection

Circulating tumor DNA detection using next-generation sequencing (NGS) data of plasma DNA is promising for cancer identification and characterization. However, the tumor signal in the blood is often low and difficult to distinguish from errors. We present DREAMS (Deep Read-level Modelling of Sequencing-errors) for estimating error rates of individual read positions. Using DREAMS, we develop statistical methods for variant calling (DREAMS-vc) and cancer detection (DREAMS-cc). For evaluation, we generate deep targeted NGS data of matching tumor and plasma DNA from 85 colorectal cancer patients. The DREAMS approach performs better than state-of-the-art methods for variant calling and cancer detection.

Sequencing and analysis of the complete mitochondrial genomes of Toona sinensis and Toona ciliata reveal evolutionary features of Toona

BACKGROUND

Toona is a critical genus in the Meliaceae, and the plants of this group are an asset for both restorative and restorative purposes, the most flexible of which are Toona sinensis and Toona ciliata. To concentrate on the advancement of mitochondrial(Mt) genome variety in T.sinensis and T.ciliata, the Mt genomes of the two species were sequenced in high throughput independently, after de novo assembly and annotation to construct a Mt genome map for comparison in genome structure. Find their repetitive sequences and analyze them in comparison with the chloroplast genome, along with Maximum-likelihood(ML) phylogenetic analysis with 16 other relatives.

RESULTS

(1) T. sinensis and T.ciliata are both circular structures with lengths of 683482 bp and 68300 bp, respectively. They share a high degree of similarity in encoding genes and have AT preferences. All of them have the largest Phe concentration and are the most frequently used codons. (2) Both of their Mt genome are highly preserved in terms of structural and functional genes, while the main variability is reflected in the length of tRNA, the number of genes, and the value of RSCU. (3) T. siniensis and T. ciliata were detected to have 94 and 87 SSRs, respectively, of which mononucleotides accounted for the absolute proportion. Besides, the vast majority of their SSRs were found to be poly-A or poly-T. (4)10 and 11 migrating fragments were identified in the comparison with the chloroplast genome, respectively. (5) In the ML evolutionary tree, T.sinensis and T.ciliata clustered individually into a small branch with 100% support, reflecting two species of Toona are very similarly related to each other.

CONCLUSIONS

This research provides a basis for the exploitation of T.sinensis and T.ciliata in terms of medicinal, edible, and timber resources to avoid confusion; at the same time, it can explore the evolutionary relationship between the Toona and related species, which does not only have an important practical value, but also provides a theoretical basis for future hybrid breeding of forest trees, molecular markers, and evolutionary aspects of plants, which has great scientific significance.

Heteroplasmic Mutant Load Differences in Mitochondrial DNA-Associated Leigh Syndrome

BACKGROUND

Mitochondrial DNA (mtDNA)-associated Leigh syndrome is influenced by mutant pathogenicity and corresponding heteroplasmic loads; however, the manner in which heteroplasmic mutant load affects patient phenotypes and the relationship between mutant types and heteroplasmic mutant loads remain unknown. We aimed to investigate the distribution of the mutant load of different mtDNA mutations in a single-center cohort.

METHODS

We used next-generation sequencing to confirm mtDNA mutations in 31 patients with Leigh syndrome. Subsequently, we counted the number of mtDNA reads to quantitatively analyze the heteroplasmic mutant load and categorize the patients according to the mtDNA mutations they harbored. Confirmed cases of mtDNA-associated Leigh syndrome were classified according to the mutations observed in six genes and 10 nucleotides.

RESULTS

Of the 31 patients with Leigh syndrome, 27 harbored known pathogenic mutations. We discovered that MT-ATP6 was the most commonly mutated gene (n = 13 patients), followed by MT-ND3 (n = 7) and MT-ND5 (n = 4). MT-ATP6 had a significantly higher mutant load than MT-ND3 and MT-ND5 (P < 0.001, each). By contrast, MT-ND5 had a significantly lower mutant load than MT-ND3 (P = 0.007). Notably, the mutation loads varied significantly among patients carrying the MT-ATP6, MT-ND3, and MT-ND5 mutations.

CONCLUSIONS

Our study illustrated the heteroplasmic diversity and phenotypic expression threshold of mutated mitochondrial genes in mtDNA-associated Leigh syndrome. The results provide promising insights into the genotype-phenotype correlation in mtDNA-associated Leigh syndrome that are expected to guide the development of tailored treatments for Leigh syndrome.

Enhancement of hexavalent chromium reduction by Shewanella oneidensis MR-1 in presence of copper nanoparticles via stimulating bacterial extracellular electron transfer and environmental adaptability

The bioreduction of hexavalent chromium (Cr(VI)) depends highly on bacterial activity, while the release of copper nanoparticles (Cu NPs) poses threats to microorganisms in the environment. This work demonstrated that Cr(VI) reduction efficiency of Shewanella oneidensis MR-1 was remarkably enhanced by 83.7% under 20 mg/L Cu NPs exposure. Cu NPs improved the electron migration capacity of Shewanella oneidensis MR-1 by enhancing bioelectrochemical performance and flavin mononucleotide secretion. Moreover, key genes related to extracellular electron transfer pathways, including direct electron transfer through outer-membrane proteins, flavin-mediated electron transfer, and conductive flagella, were generally upregulated under Cu NPs exposure. In addition, environmental adaptability of Shewanella oneidensis MR-1 was enhanced under Cu NPs exposure by improving environmental information processing and energy and reducing power production, promoting Cr(VI) reduction by Shewanella oneidensis MR-1. This work indicated that Cu NPs could enhance Cr(VI) reduction by Shewanella oneidensis MR-1 through regulating extracellular electron transfer and environmental adaptability.

Subspecies Classification and Comparative Genomic Analysis of Lactobacillus kefiranofaciens HL1 and M1 for Potential Niche-Specific Genes and Pathways

(1) Background: Strains HL1 and M1, isolated from kefir grains, have been tentatively identified, based on their partial 16S rRNA gene sequences, as Lactobacillus kefiranofaciens. The two strains demonstrated different health benefits. Therefore, not only the genetic factors exerting diverse functionalities in different L. kefiranofaciens strains, but also the potential niche-specific genes and pathways among the L. kefiranofaciens strains, should be identified. (2) Methods: Phenotypic and genotypic approaches were employed to identify strains HL1 and M1 at the subspecies level. For the further characterization of the probiotic properties of both strains, comparative genomic analyses were used. (3) Results: Both strains were identified as L. kefiranofaciens subsp. kefirgranum. According to the COG function category, dTDP-rhamnose and rhamnose-containing glycans were specifically detected in the L. kefiranofaciens subsp. Kefirgranum genomes. Three unique genes (epsI, epsJ, and epsK) encoding glycosyltransferase in the EPS gene cluster, and the ImpB/MucB/SamB family protein encoding gene were found in HL1 and M1. The specific ability to degrade arginine via the ADI pathway was found in HL1. The presence of the complete glycogen metabolism (glg) operon in the L. kefiranofaciens strains suggested the importance of glycogen synthesis to enable colonization in kefir grains and extend survival under environmental stresses. (4) Conclusions: The obtained novel information on the potential genes and pathways for polysaccharide synthesis and other functionalities in our HL1 and M1 strains could be applied for further functionality predictions for potential probiotic screening.

The pediatric urobiome in genitourinary conditions: a narrative review

The microbial ecosystem within the bladder that can be measured within the urine, or urobiome, is an emerging field of study with little published data regarding children. However, investigations into urobiome research have the potential to significantly impact the understanding of the pathophysiology of genitourinary conditions, as well as potentially identify novel therapeutics. Therefore, both researchers and clinicians should be aware of pediatric urobiome research. The purpose of this review is to highlight the literature around urobiome research in urinary tract infections, nephrolithiasis, and neurogenic bladder; comment on pediatric-specific considerations when reading and interpreting the urobiome literature; and to identify new potential areas of research.

Comparative Genomics of Acetic Acid Bacteria within the Genus Bombella in Light of Beehive Habitat Adaptation

It is known that the bacterial microbiota in beehives is essential for keeping bees healthy. Acetic acid bacteria of the genus Bombella colonize several niches in beehives and are associated with larvae protection against microbial pathogens. We have analyzed the genomes of 22 Bombella strains of different species isolated in eight different countries for taxonomic affiliation, central metabolism, prophages, bacteriocins and tetracycline resistance to further elucidate the symbiotic lifestyle and to identify typical traits of acetic acid bacteria. The genomes can be assigned to four different species. Three genomes show ANIb values and DDH values below species demarcation values to any validly described species, which identifies them as two potentially new species. All Bombella spp. lack genes in the Embden–Meyerhof–Parnas pathway and the tricarboxylic acid cycle, indicating a focus of intracellular carbohydrate metabolism on the pentose phosphate pathway or the Entner–Doudoroff pathway for which all genes were identified within the genomes. Five membrane-bound dehydrogenases were identified that catalyze oxidative fermentation reactions in the periplasm, yielding oxidative energy. Several complete prophages, but no bacteriocins, were identified. Resistance to tetracycline, used to prevent bacterial infections in beehives, was only found in Bombella apis MRM1T. Bombella strains exhibit increased osmotolerance in high glucose concentrations compared to Gluconobacter oxydans, indicating adaption to high sugar environments such as beehives.

Microbial communities in swamps of four mangrove reserves driven by interactions between physicochemical properties and microbe in the North Beibu Gulf, China

Mangroves are distributed in coastal and estuarine regions and are characterized as a sink for terrestrial pollution. It is believed that complex interactions between environmental factors and microbial communities exist in mangrove swamps. However, little is known about environment-microbe interactions. There is a need to clarify some important environmental factors shaping microbial communities and how environmental factors interact with microbial assemblages in mangrove swamps. In the present study, physicochemical and microbial characteristics in four mangrove reserves (named ZZW, Qin, Bei, and GQ) in the North Beibu Gulf were determined. The interactions between environmental factors and microbial assemblages were analyzed with statistical methods in addition to CCA and RDA. Higher concentrations of sulfate (SO₄^2--S) and Fe but lower concentrations of total phosphorus (TP) and NO₃^--N were detected in ZZW and Qin. Nutrient elements (NO₃^--N, NH₄⁺-N, organic matter (OM), SO₄^2--S, Fe, and TP) were more important than heavy metals for determining the microbial assemblages, and NO₃^--N was the most important factor. NO₃^--N, SO₄^2--S, TP, and Fe formed a significant co-occurrence network in conjunction with some bacterial taxa, most of which were Proteobacteria. Notably, comparatively elevated amounts of sulfate-reducing bacteria (Desulfatibacillum, Desulfomonile, and Desulfatiglans) and sulfur-oxidizing bacteria (Thioprofundum and Thiohalophilus) were found in ZZW and Qin. The co-occurrence network suggested that some bacteria involved in sulfate reduction and sulfur oxidation drive the transformation of P and N, resulting in the reduction of P and N in mangrove swamps. Through the additional utilization of multivariate regression tree (MRT) and co-occurrence network analysis, our research provides a new perspective for understanding the interactions between environmental factors and microbial communities in mangroves.

Genotype-phenotype analysis of MT-ATP6-associated Leigh syndrome

OBJECTIVES

Mitochondrial DNA (mtDNA)-associated Leigh syndrome (LS) is characterized by maternal inheritance, and the heteroplasmic mutant load of mtDNA pathogenic variants is known to affect clinical phenotypes. Among mtDNA pathogenic variants, variants of the MT-ATP6 gene account for most of reported cases. In this report, we aimed to describe the clinical and genetic findings of MT-ATP6-associated LS patients diagnosed at a single tertiary institution in Korea.

METHODS

Thirteen patients with genetically confirmed MT-ATP6-associated LS were selected. We reviewed each patient's clinical findings, including general characteristics, biochemical parameters, brain MR images, muscle biopsy results, and heteroplasmic mutant load over a long-term follow-up period.

RESULTS

MT-ATP6-associated LS was of predominantly early onset (age <2 years), although we identified 2 late-onset (>60 months) LS patients. The heteroplasmic mutant load estimated by next-generation sequencing was 96%-100% in all nucleotide change groups. Compared with other forms of MT-ATP6-associated LS, the m.8993T>G point mutation elicited a significantly higher rate of symptom onset before 2 years of age. Brain MRI showed bilateral basal ganglia involvement in all patients, followed by cerebral atrophy, brainstem and thalamus involvement, and cerebellar atrophy. After follow-up (median 7.2 years, range 1.4 to 11.5 years), LS with m.8993T>G point mutations had a slightly more severe clinical progression compared with other forms of MT-ATP6-associated LS.

CONCLUSIONS

MT-ATP6-associated LS patients presented with a broad spectrum of clinical diagnoses and had a very high heteroplasmic mutant load. This study provides valuable data on MT-ATP6-associated LS that will inform subsequent studies on LS.

Optimization of enzymatic fragmentation is crucial to maximize genome coverage: a comparison of library preparation methods for Illumina sequencing

Background

Novel commercial kits for whole genome library preparation for next-generation sequencing on Illumina platforms promise shorter workflows, lower inputs and cost savings. Time savings are achieved by employing enzymatic DNA fragmentation and by combining end-repair and tailing reactions. Fewer cleanup steps also allow greater DNA input flexibility (1 ng-1 μg), PCR-free options from 100 ng DNA, and lower price as compared to the well-established sonication and tagmentation-based DNA library preparation kits.

Results

We compared the performance of four enzymatic fragmentation-based DNA library preparation kits (from New England Biolabs, Roche, Swift Biosciences and Quantabio) to a tagmentation-based kit (Illumina) using low input DNA amounts (10 ng) and PCR-free reactions with 100 ng DNA. With four technical replicates of each input amount and kit, we compared the kits’ fragmentation sequence-bias as well as performance parameters such as sequence coverage and the clinically relevant detection of single nucleotide and indel variants. While all kits produced high quality sequence data and demonstrated similar performance, several enzymatic fragmentation methods produced library insert sizes which deviated from those intended. Libraries with longer insert lengths performed better in terms of coverage, SNV and indel detection. Lower performance of shorter-insert libraries could be explained by loss of sequence coverage to overlapping paired-end reads, exacerbated by the preferential sequencing of shorter fragments on Illumina sequencers. We also observed that libraries prepared with minimal or no PCR performed best with regard to indel detection.

Conclusions

The enzymatic fragmentation-based DNA library preparation kits from NEB, Roche, Swift and Quantabio are good alternatives to the tagmentation based Nextera DNA flex kit from Illumina, offering reproducible results using flexible DNA inputs, quick workflows and lower prices. Libraries with insert DNA fragments longer than the cumulative sum of both read lengths avoid read overlap, thus produce more informative data that leads to strongly improved genome coverage and consequently also increased sensitivity and precision of SNP and indel detection. In order to best utilize such enzymatic fragmentation reagents, researchers should be prepared to invest time to optimize fragmentation conditions for their particular samples.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08316-y.

Sediment-associated microbial community profiling: sample pre-processing through sequential membrane filtration for 16S rRNA amplicon sequencing

BACKGROUND

Sequential membrane filtration as a pre-processing step for capturing sediment-associated microorganisms could provide good quality and integrity DNA that can be preserved and kept at ambient temperatures before community profiling through culture-independent molecular techniques. However, the effects of sample pre-processing via filtration on DNA-based profiling of sediment-associated microbial community diversity and composition are poorly understood. Specifically, the influences of pre-processing on the quality and quantity of extracted DNA, high-throughput DNA sequencing reads, and detected microbial taxa need further evaluation.

RESULTS

We assessed the impact of pre-processing freshwater sediment samples by sequential membrane filtration (from 10, 5 to 0.22 μm pore size) for 16S rRNA-based community profiling of sediment-associated microorganisms. Specifically, we examined if there would be method-driven differences between non- and pre-processed sediment samples regarding the quality and quantity of extracted DNA, PCR amplicon, resulting high-throughput sequencing reads, microbial diversity, and community composition. We found no significant difference in the qualities and quantities of extracted DNA and PCR amplicons, and the read abundance after bioinformatics processing (i.e., denoising and chimeric-read filtering steps) between the two methods. Although the non- and pre-processed sediment samples had more unique than shared amplicon sequence variants (ASVs), we report that their shared ASVs accounted for 74% of both methods' absolute read abundance. More so, at the genus level, the final collection filter identified most of the genera (95% of the reads) captured from the non-processed samples, with a total of 51 false-negative (2%) and 59 false-positive genera (3%). We demonstrate that while there were differences in shared and unique taxa, both methods revealed comparable microbial diversity and community composition.

CONCLUSIONS

Our observations highlight the feasibility of pre-processing sediment samples for community analysis and the need to further assess sampling strategies to help conceptualize appropriate study designs for sediment-associated microbial community profiling.

Hackflex: low-cost, high-throughput, Illumina Nextera Flex library construction

We developed a low-cost method for the production of Illumina-compatible sequencing libraries that allows up to 14 times more libraries for high-throughput Illumina sequencing to be generated for the same cost. We call this new method Hackflex. The quality of library preparation was tested by constructing libraries from Escherichia coli MG1655 genomic DNA using either Hackflex, standard Nextera Flex (recently renamed as Illumina DNA Prep) or a variation of standard Nextera Flex in which the bead-linked transposase is diluted prior to use. In order to test the library quality for genomes with a higher and a lower G+C content, library construction methods were also tested on Pseudomonas aeruginosa PAO1 and Staphylococcus aureus ATCC 25923, respectively. We demonstrated that Hackflex can produce high-quality libraries and yields a highly uniform coverage, equivalent to the standard Nextera Flex kit. We show that strongly size-selected libraries produce sufficient yield and complexity to support de novo microbial genome assembly, and that assemblies of the large-insert libraries can be much more contiguous than standard libraries without strong size selection. We introduce a new set of sample barcodes that are distinct from standard Illumina barcodes, enabling Hackflex samples to be multiplexed with samples barcoded using standard Illumina kits. Using Hackflex, we were able to achieve a per-sample reagent cost for library prep of A$7.22 (Australian dollars) (US $5.60; UK £3.87, £1=A$1.87), which is 9.87 times lower than the standard Nextera Flex protocol at advertised retail price. An additional simple modification and further simplification of the protocol by omitting the wash step enables a further price reduction to reach an overall 14-fold cost saving. This method will allow researchers to construct more libraries within a given budget, thereby yielding more data and facilitating research programmes where sequencing large numbers of libraries is beneficial.

Co-occurrence of cfr-mediated linezolid-resistance in ST398 LA-MRSA and non-aureus staphylococci isolated from a pig farm

Linezolid resistance, mediated by the cfr gene, which confers resistant phenotypes to phenicols, lincosamides, oxazolidinones, pleuromutilins, and streptogramin A antimicrobials, has emerged in S. aureus and non-aureus staphylococci (NAS). Moreover, due to the transferable potential via plasmids, the spread of cfr among staphylococci is of great concern. In the present study, we investigated the prevalence of cfr-mediated linezolid resistance in ST398 methicillin-resistant S. aureus (MRSA) and NAS strains isolated from a pig farm. Among the 26 staphylococci isolates collected from a pig farm, 14 cfr-harboring ST398 MRSA and NAS (S. epidermidis, S. pasteuri, S. cohnii, and S. rostri) strains were resistant to linezolid and also carried the fexA gene. Comparative genome analysis of cfr-carrying linezolid-resistant ST398 MRSA and NAS (S. pasteuri, S. cohnii, and S. epidermidis) strains revealed that the segments harboring cfr in different staphylococcal strains showed ≥ 99 % sequence identity and the corresponding region containing the cfr, fexA, and Tn558 elements were located in a 38-kb plasmid, designated pSA12 of ST398 MRSA. These observations indicate that the cfr-carrying plasmids and/or fragments may be disseminated among staphylococci in a pig farm and possibly transmitted to staphylococci of human origin, subsequently posing a threat to public health. This is the first report of the co-existence of cfr in linezolid-resistant ST398 MRSA and NAS isolated from a pig farm in South Korea.

Association Between Epilepsy and Leigh Syndrome With MT-ND3 Mutation, Particularly the m.10191T>C Point Mutation

Background and Purpose: Recent advances in molecular genetic testing have led to a rapid increase in the understanding of the genetics of Leigh syndrome. Several studies have suggested that Leigh syndrome with MT-ND3 mutation is strongly associated with epilepsy. This study focused on the epilepsy-related characteristics of Leigh syndrome with MT-ND3 mutation identified in a single tertiary hospital in South Korea. Methods: We selected 31 patients with mitochondrial DNA (mtDNA) mutations who were genetically diagnosed with mtDNA-associated Leigh syndrome. Among them, seven patients with MT-ND3 mutations were detected. We reviewed various clinical findings such as laboratory findings, brain images, electroencephalography data, seizure types, seizure frequency, antiepileptic drug use history, and current seizure status. Results: The nucleotide changes in the seven patients with the Leigh syndrome with MT-ND3 mutation were divided into two groups: m.10191T>C and m.10158T>C. Six of the seven patients were found to have the m.10191T>C mutations. The median value of the mutant load was 82.5%, ranging from 57.9 to 93.6%. No particular tendency was observed for the first symptom or seizure onset or mutant load. The six patients with the m.10191T>C mutation were diagnosed with epilepsy. Three of these patients were diagnosed with Lennox-Gastaut syndrome (LGS). Conclusion: We reported a very strong association between epilepsy and MT-ND3 mutation in Leigh syndrome, particularly the m.10191T>C mutation. The possibility of an association between the epilepsy phenotype of the m.10191T>C mutation and LGS was noted.

Evaluating whole-genome sequencing quality metrics for enteric pathogen outbreaks

Background

Whole genome sequencing (WGS) has gained increasing importance in responses to enteric bacterial outbreaks. Common analysis procedures for WGS, single nucleotide polymorphisms (SNPs) and genome assembly, are highly dependent upon WGS data quality.

Methods

Raw, unprocessed WGS reads from Escherichia coli, Salmonella enterica, and Shigella sonnei outbreak clusters were characterized for four quality metrics: PHRED score, read length, library insert size, and ambiguous nucleotide composition. PHRED scores were strongly correlated with improved SNPs analysis results in E. coli and S. enterica clusters.

Results

Assembly quality showed only moderate correlations with PHRED scores and library insert size, and then only for Salmonella. To improve SNP analyses and assemblies, we compared seven read-healing pipelines to improve these four quality metrics and to see how well they improved SNP analysis and genome assembly. The most effective read healing pipelines for SNPs analysis incorporated quality-based trimming, fixed-width trimming, or both. The Lyve-SET SNPs pipeline showed a more marked improvement than the CFSAN SNP Pipeline, but the latter performed better on raw, unhealed reads. For genome assembly, SPAdes enabled significant improvements in healed E. coli reads only, while Skesa yielded no significant improvements on healed reads.

Conclusions

PHRED scores will continue to be a crucial quality metric albeit not of equal impact across all types of analyses for all enteric bacteria. While trimming-based read healing performed well for SNPs analyses, different read healing approaches are likely needed for genome assembly or other, emerging WGS analysis methodologies.

Genomic variation in captive deer mouse (Peromyscus maniculatus) populations

BACKGROUND

Deer mice (genus Peromyscus) are the most common rodents in North America. Despite the availability of reference genomes for some species, a comprehensive database of polymorphisms, especially in those maintained as living stocks and distributed to academic investigators, is missing. In the present study we surveyed two populations of P. maniculatus that are maintained at the Peromyscus Genetic Stock Center (PGSC) for polymorphisms across their 2.5 × 109 bp genome.

RESULTS

High density of variation was identified, corresponding to one SNP every 55 bp for the high altitude stock (SM2) or 207 bp for the low altitude stock (BW) using snpEff (v4.3). Indels were detected every 1157 bp for BW or 311 bp for SM2. The average Watterson estimator for the BW and SM2 populations is 248813.70388 and 869071.7671 respectively. Some differences in the distribution of missense, nonsense and silent mutations were identified between the stocks, as well as polymorphisms in genes associated with inflammation (NFATC2), hypoxia (HIF1a) and cholesterol metabolism (INSIG1) and may possess value in modeling pathology.

CONCLUSIONS

This genomic resource, in combination with the availability of P. maniculatus from the PGSC, is expected to promote genetic and genomic studies with this animal model.

Facklamia lactis sp. nov., isolated from raw milk

Two strains of a Gram-staining-positive species were isolated from German bulk tank milk. On the basis of their 16S rRNA sequences they were affiliated to the genus Facklamia but could not be assigned to any species with a validly published name. Facklamia miroungae ATCC BAA-466^T (97.3 % 16S rRNA sequence similarity), Facklamia languida CCUG 37842^T (96.9 %), and Facklamia hominis CCUG 36813^T (96.6 %) are the closest relatives. In the 16S rRNA phylogeny and in the core-genome phylogeny strains WS 5301^T and WS 5302 form a well-supported, separate lineage. Pairwise average nucleotide identity calculated using MUMmer (ANIm) between WS 5301^T and type strains of other Facklamia species is well below the species cut-off (95 %) and ranges from 83.4 to 87.7 %. The DNA G+C content of the type strain is 36.4 mol% and the assembly size of the genome is 2.2 Mb. Cells of WS 5301^T are non-motile, non-endospore-forming, oxidase-negative, catalase-negative and facultatively anaerobic cocci. The fastidious species grows at 10-40 °C and with up to 7.0 % (w/v) NaCl in BHI supplemented with 5 g l^-1 yeast extract. Major polar lipids are phosphatidylglycerol, diphosphatidylglycerol and two glycolipids. Predominant fatty acids are C_16 : 1ω9c and C_18 : 1ω9c. On the basis of their genomic, physiological and chemotaxonomic characteristics the strains examined in this study represent the same, hitherto unknown species. We propose the name Facklamia lactis sp. nov. for which WS 5301^T (=DSM 111018^T=LMG 31861^T) is the type strain and WS 5302 (=DSM 111019=LMG 31862) is an additional strain of this novel species.

Evaluating coverage bias in next-generation sequencing of Escherichia coli

Whole-genome sequencing is essential to many facets of infectious disease research. However, technical limitations such as bias in coverage and tagmentation, and difficulties characterising genomic regions with extreme GC content have created significant obstacles in its use. Illumina has claimed that the recently released DNA Prep library preparation kit, formerly known as Nextera Flex, overcomes some of these limitations. This study aimed to assess bias in coverage, tagmentation, GC content, average fragment size distribution, and de novo assembly quality using both the Nextera XT and DNA Prep kits from Illumina. When performing whole-genome sequencing on Escherichia coli and where coverage bias is the main concern, the DNA Prep kit may provide higher quality results; though de novo assembly quality, tagmentation bias and GC content related bias are unlikely to improve. Based on these results, laboratories with existing workflows based on Nextera XT would see minor benefits in transitioning to the DNA Prep kit if they were primarily studying organisms with neutral GC content.

Draft Genome Sequence of Mobilitalea sibirica Strain P3M-3T, the Sole Representative of the Genus Mobilitalea

Mobilitalea sibirica strain P3M-3^T is a strictly anaerobic, halotolerant, organotrophic bacterium of the family Lachnospiraceae that can utilize various plant-derived polysaccharides as its carbon source. The strain was originally isolated by from a microbial mat in western Siberia (Russia). In this study, we present the draft genome sequence of M. sibirica P3M-3^T based on Illumina paired-end sequencing.

Mobilitalea sibirica strain P3M-3^T is a strictly anaerobic, halotolerant, organotrophic bacterium of the family Lachnospiraceae that can utilize various plant-derived polysaccharides as its carbon source. The strain was originally isolated from a microbial mat in western Siberia (Russia). In this study, we present the draft genome sequence of M. sibirica P3M-3^T based on Illumina paired-end sequencing.

Bombella favorum sp. nov. and Bombella mellum sp. nov., two novel species isolated from the honeycombs of Apis mellifera

As part of a study investigating the microbiome of bee hives and honey, two novel strains (TMW 2.1880^T and TMW 2.1889^T) of acetic acid bacteria were isolated and subsequently taxonomically characterized by a polyphasic approach, which revealed that they cannot be assigned to known species. The isolates are Gram-stain-negative, aerobic, pellicle-forming, catalase-positive and oxidase-negative. Cells of TMW 2.1880^T are non-motile, thin/short rods, and cells of TMW 2.1889^T are motile and occur as rods and long filaments. Morphological, physiological and phylogenetic analyses revealed a distinct lineage within the genus Bombella. Strain TMW 2.1880^T is most closely related to the type strain of Bombella intestini with a 16S rRNA gene sequence similarity of 99.5 %, and ANIb and in silico DDH values of 94.16 and 56.3 %, respectively. The genome of TMW 2.1880^T has a size of 1.98 Mb and a G+C content of 55.3 mol%. Strain TMW 2.1889^T is most closely related to the type strain of Bombella apis with a 16S rRNA gene sequence similarity of 99.5 %, and ANIb and in silico DDH values of 85.12 and 29.5 %, respectively. The genome of TMW 2.1889^T has a size of 2.07 Mb and a G+C content of 60.4 mol%. Ubiquinone analysis revealed that both strains contained Q-10 as the main respiratory quinone. Major fatty acids for both strains were C_16 : 0, C_19 : 0 cyclo ω8c and summed feature 8, respectively, and additionally C_14 : 0 2-OH only for TMW 2.1880^T and C_14 : 0 only for TMW 2.1889^T. Based on polyphasic evidence, the two isolates from honeycombs of Apis mellifera represent two novel species of the genus Bombella, for which the names Bombella favorum sp. nov and Bombella mellum sp. nov. are proposed. The designated respective type strains are TMW 2.1880^T (=LMG 31882^T=CECT 30114^T) and TMW 2.1889^T (=LMG 31883^T=CECT 30113^T).

Genomic Information on Linezolid-Resistant Sequence-Type 398 Livestock-Associated Methicillin-Resistant Staphylococcus aureus Isolated from a Pig

The frequent occurrence of sequence-type 398 (ST398) livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) in pigs has become a major public health concern owing to the increased zoonotic potential of the pathogen. Recently, a novel oxazolidinone resistance gene, chloramphenicol-florfenicol resistant (cfr), conferring multiresistance phenotypes to phenicols, lincosamides, oxazolidinones, pleuromutilins, and streptogramin A (PhLOPS_A), has been found among ST398 LA-MRSA strains isolated from pigs. In this study, we report the first in silico genome analysis of a linezolid-resistant ST398 LA-MRSA strain, designated PJFA-521M, recovered from a pig in Korea. Genomic analyses revealed that the presence of the cfr gene was responsible for the observed linezolid resistance in the PJFA-521M strain. Moreover, newer antimicrobial resistance genes, such as the dfrG, aadE, spw, lsa(E), lnu(B), and fexA genes, were found in the PJFA-521M strain. In addition to the genetic elements for antimicrobial resistance, the carriage of various virulence genes for adherence, invasion, and immunomodulation was identified in the genome, especially within several mobile genetic elements (MGEs). The presence of multiple antimicrobial resistance genes and virulence genes on MGEs in the genome of a linezolid-resistant ST398 LA-MRSA should raise awareness regarding the use of other antimicrobial agents in pig farms and may also provide selective pressure for the prevalence of the cfr gene and the associated multidrug-resistant phenotype.

A collection of bacterial isolates from the pig intestine reveals functional and taxonomic diversity

Our knowledge about the gut microbiota of pigs is still scarce, despite the importance of these animals for biomedical research and agriculture. Here, we present a collection of cultured bacteria from the pig gut, including 110 species across 40 families and nine phyla. We provide taxonomic descriptions for 22 novel species and 16 genera. Meta-analysis of 16S rRNA amplicon sequence data and metagenome-assembled genomes reveal prevalent and pig-specific species within Lactobacillus, Streptococcus, Clostridium, Desulfovibrio, Enterococcus, Fusobacterium, and several new genera described in this study. Potentially interesting functions discovered in these organisms include a fucosyltransferase encoded in the genome of the novel species Clostridium porci, and prevalent gene clusters for biosynthesis of sactipeptide-like peptides. Many strains deconjugate primary bile acids in in vitro assays, and a Clostridium scindens strain produces secondary bile acids via dehydroxylation. In addition, cells of the novel species Bullifex porci are coccoidal or spherical under the culture conditions tested, in contrast with the usual helical shape of other members of the family Spirochaetaceae. The strain collection, called ‘Pig intestinal bacterial collection’ (PiBAC), is publicly available at www.dsmz.de/pibac and opens new avenues for functional studies of the pig gut microbiota.

The authors present a public collection of 117 bacterial isolates from the pig gut, including the description of 38 novel taxa. Interesting functions discovered in these organisms include a new fucosyltransferease and sactipeptide-like molecules encoded by biosynthetic gene clusters.

Pseudomonas cremoris sp. nov., a novel proteolytic species isolated from cream

During a study investigating the microbiota of raw milk and its semi-finished products, strains WS 5106^T and WS 5096 were isolated from cream and skimmed milk concentrate. They could be assigned to the genus Pseudomonas by their 16S rRNA sequences, but not to any validly named species. In this work, a polyphasic approach was used to characterize the novel strains and to investigate their taxonomic status. Examinations based on the topology of core genome phylogenomy as well as average nucleotide identity (ANIm) comparisons suggested a novel Pseudomonas species within the Pseudomonas fluorescens subgroup. With pairwise ANIm values of 90.1 and 89.8 %, WS 5106^T was most closely related to Pseudomonas nabeulensis CECT 9765^T and Pseudomonas kairouanensis CECT 9766^T. The G+C content of strain WS 5106^T was 60.1 mol%. Morphologic analyses revealed Gram-stain-negative, aerobic, catalase and oxidase positive, rod-shaped and motile cells. Proteolysis on skimmed milk agar as well as lipolysis on tributyrin agar occurred at both 28 and 6 °C. Tolerated growth conditions were temperatures between 4 and 34 °C, pH values between 6.0 and 8.0, and salt concentrations of up to 5 %. Fatty acid profiles showed a pattern typical for Pseudomonas, with C_16 : 0 as the dominant component. The major cellular polar lipids were phosphatidylethanolamine, phosphatidylglycerol and diphosphatidylglycerol and the dominating quinone was Q-9. Based on these results, it is proposed to classify the strains as a novel species, Pseudomonas cremoris sp. nov., with WS 5106^T (=DSM 111143^T=LMG 31863^T) as type strain and WS 5096 (=DSM 111129=LMG 31864) as an additional strain.

Lactococcus carnosus sp. nov. and Lactococcus paracarnosus sp. nov., two novel species isolated from modified-atmosphere packaged beef steaks

As part of a study investigating the spoilage microbiome of modified-atmosphere packaged beef from Germany, four novel strains of lactic acid bacteria were isolated and subsequently taxonomically characterized by a polyphasic approach, which revealed that they could not be assigned to known species. The isolates were Gram-staining-positive, coccoid, facultatively anaerobic, non-motile, catalase-negative and oxidase-negative. Morphological, physiological and phylogenetic analysis revealed a distinct lineage within the genus Lactococcus , with Lactococcus piscium and Lactococcus plantarum as closest relatives. Results indicated that they represented two different novel species with two strains, (TMW 2.1612T/TMW 2.1613 and TMW 2.1615T/TMW 2.1614), respectively. The two strains of both novel species shared identical 16S rRNA gene sequences but a MLSA allowed their intraspecies differentiation. The 16S rRNA gene sequences of TMW 2.1612T and TMW 2.1615T had a similarity of 99.85 % to each other and a similarity of 99.85 and 99.78 % the most closely related type strain of Lactococcus piscium , respectively. However, the ANIb value between the respective type strains TMW 2.1612T and TMW 2.1615T, and the type strain of Lactococcus piscium was only 94.3 and 92.0 %, respectively, and 92.9 % between TMW 2.1612T and TMW 2.1615T. The in silico DDH estimate value between the respective type strain TMW 2.1612T and TMW 2.1615T and the most closely related type strain of Lactococcus piscium was only 59.9 and 48.9 %, respectively, and 51.1 % between TMW 2.1612T and TMW 2.1615T. Peptidoglycan type of strain TMW 2.1612T is Lys–Thr–Ala and major fatty acids are summed feature 8 and C16 : 0. Peptidoglycan type of strain TMW 2.1615T is Lys–Ala and major fatty acids are C16 : 0, C19 : 0cyclo ω8c and summed feature 8. On the basis of polyphasic evidence, the meat isolates represent two novel species of the genus Lactococcus , for which the names Lactococcus carnosus sp. nov and Lactococcus paracarnosus sp. nov are proposed. The designated respective type strains are TMW 2.1612T (=DSM 111016T =CECT 30115T) and TWM 2.1615T (=DSM 111017T =CECT 30116T).

Fundicoccus ignavus gen. nov., sp. nov., a novel genus of the family Aerococcaceae isolated from bulk tank milk

Three strains of a Gram-stain-positive, catalase-negative, facultative anaerobic, and coccoid species were isolated from German bulk tank milk. Phylogenetic analyses based on the 16S rRNA gene sequences indicated that the three strains (WS4937^T, WS4759 and WS5303) constitute an independent phylogenetic lineage within the family Aerococcaceae with Facklamia hominis CCUG 36813^T (93.7-94.1 %) and Eremococcus coleocola M1831/95/2^T (93.5 %) as most closely related type species. The unclassified strains demonstrated variable growth with 6.5 % (w/v) NaCl and tolerated pH 6.5-9.5. Growth was observed from 12 to 39 °C. Their cell-wall peptidoglycan belongs to the A1α type (l-Lys-direct) consisting of alanine, glutamic acid and lysine. The predominant fatty acids were C_16 : 1 ω9c, C_16 : 0 and C_18 : 1 ω9c and in the polar lipids profile three glycolipids, a phospholipid, phosphatidylglycerol, phosphoglycolipid and diphosphatidylglycerol were found. The G+C content of strain WS4937^T was 37.4 mol% with a genome size of ~3.0 Mb. Based on phylogenetic, phylogenomic and biochemical characterizations, the isolates can be demarcated from all other genera of the family Aerococcaceae and, therefore, the novel genus Fundicoccus gen. nov. is proposed. The type species of the novel genus is Fundicoccus ignavus gen. nov., sp. nov. WS4937^T (=DSM 109652^T=LMG 31441^T).

Importance of Defluviitalea raffinosedens for Hydrolytic Biomass Degradation in Co-Culture with Hungateiclostridium thermocellum

Bacterial hydrolysis of polysaccharides is an important step for the production of sustainable energy, for example during the conversion of plant biomass to methane-rich biogas. Previously, Hungateiclostridium thermocellum was identified as cellulolytic key player in thermophilic biogas microbiomes with a great frequency as an accompanying organism. The aim of this study was to physiologically characterize a recently isolated co-culture of H. thermocellum and the saccharolytic bacterium Defluviitalea raffinosedens from a laboratory-scale biogas fermenter. The characterization focused on cellulose breakdown by applying the measurement of cellulose hydrolysis, production of metabolites, and the activity of secreted enzymes. Substrate degradation and the production of volatile metabolites was considerably enhanced when both organisms acted synergistically. The metabolic properties of H. thermocellum have been studied well in the past. To predict the role of D. raffinosedens in this bacterial duet, the genome of D. raffinosedens was sequenced for the first time. Concomitantly, to deduce the prevalence of D. raffinosedens in anaerobic digestion, taxonomic composition and transcriptional activity of different biogas microbiomes were analyzed in detail. Defluviitalea was abundant and metabolically active in reactor operating at highly efficient process conditions, supporting the importance of this organism for the hydrolysis of the raw substrate.

Genetic Organization of the aprX-lipA2 Operon Affects the Proteolytic Potential of Pseudomonas Species in Milk

Psychrotolerant Pseudomonas species are a main cause of proteolytic spoilage of ultra-high temperature (UHT) milk products due to the secretion of the heat-resistant metallopeptidase AprX, which is encoded by the first gene of the aprX-lipA2 operon. While the proteolytic property has been characterized for many different Pseudomonas isolates, the underlying aprX-lipA2 gene organization was only described for a few strains so far. In this study, the phylogenomic analysis of 185 Pseudomonas type strains revealed that the presence of aprX is strongly associated to a monophylum composed of 81 species, of which 83% carried the aprX locus. Furthermore, almost all type strains of known milk-relevant species were shown to be members of the three monophyletic groups P. fluorescens, P. gessardii, and P. fragi. In total, 22 different types of aprX-lipA2 genetic organizations were identified in the genus, whereby 31% of the species tested carried the type 1 operon structure consisting of eight genes (aprXIDEF prtAB lipA2). Other genetic structures differed from type 1 mainly in the presence and location of genes coding for two lipases (lipA1 and lipA2) and putative autotransporters (prtA and prtB). The peptidase activity of 129 strains, as determined on skim milk agar and in UHT-milk, correlated largely with different aprX-lipA2 gene compositions. Particularly, isolates harboring the type 1 operon were highly proteolytic, while strains with other operon types, especially ones lacking prtA and prtB, exhibited significantly lower peptidase activities. In conclusion, the phylogenomic position and the aprX-lipA2 gene organization specify the proteolytic potential of Pseudomonas isolates. In addition, however, an interplay of several environmental factors and intrinsic traits influences production and activity of AprX, leading to strain-specific proteolytic phenotypes.

Application of Next-Generation Sequencing in Neurodegenerative Diseases: Opportunities and Challenges

Genetic factors (gene mutations) lead to various rare and prevalent neurological diseases. Identification of underlying mutations in neurodegenerative diseases is of paramount importance due to the heterogeneous nature of the genome and different clinical manifestations. An early and accurate molecular diagnosis are cardinal for neurodegenerative patients to undergo proper therapeutic regimens. The next-generation sequencing (NGS) method examines up to millions of sequences at a time. As a result, the rare molecular diagnoses, previously presented with "unknown causes", are now possible in a short time. This method generates a large amount of data that can be utilized in patient management. Since each person has a unique genome, the NGS has transformed diagnostic and therapeutic strategies into sequencing and individual genomic mapping. However, this method has disadvantages like other diagnostic methods. Therefore, in this review, we aimed to briefly summarize the NGS method and correlated studies to unravel the genetic causes of neurodegenerative diseases including Alzheimer's disease, Parkinson's disease, epilepsy, and MS. Finally, we discuss the NGS challenges and opportunities in neurodegenerative diseases.

Bifidobacterium tibiigranuli sp. nov. isolated from homemade water kefir

Two Bifidobacterium strains, TMW 2.2057^T and TMW 2.1764 were isolated from two different homemade water kefirs from Germany. Both strains were oxidase- and catalase-negative and Gram-staining-positive. Cells were non-motile, irregular rods that were aerotolerant anaerobes. On basis of fructose 6-phosphate phosphoketolase activity, they were assigned to the family Bifidobacteriaceae. Comparative analysis of 16S rRNA and concatenated housekeeping genes (clpC, dnaB, dnaG, dnaJ, hsp60 and rpoB) demonstrated that both strains represented a member of the genus Bifidobacterium, with Bifidobacterium subtile DSM 20096^T as the closest phylogenetic relative (98.35 % identity). Both strains can be distinguished using randomly amplified polymorphic DNA fingerprinting. Analysis of concatenated marker gene sequences as well as average nucleotide identity by blast (ANIb) and in silico DNA-DNA hybridization (isDDH) calculations of their genome sequences confirmed Bifidobacterium subtile DSM 20096^T as the closest relative (87.91 and 35.80 % respectively). All phylogenetic analyses allow differentiation of strains TMW 2.2057^T and TMW 2.1764 from all hitherto described species of the genus Bifidobacterium with validly published names. We therefore propose a novel species with the name Bifidobacterium tibiigranuli, for which TMW 2.2057^T (=DSM 108414^T=LMG 31086^T) is the type strain.

Pseudomonas haemolytica sp. nov., isolated from raw milk and skimmed milk concentrate

Two strains, WS 5063T and WS 5067, isolated from raw cow’s milk and skimmed milk concentrate, could be affiliated as members of the same, hitherto unknown, Pseudomonas species by 16S rRNA and rpoD gene sequences. Multilocus sequence and average nucleotide identity (ANIm) analyses based on draft genome sequences confirmed the discovery of a novel Pseudomonas species. It was most closely related to Pseudomonas synxantha DSM 18928T with an ANIm of 91.4 %. The DNA G+C content of WS 5063T was 60.0 mol %. Phenotypic characterizations showed that the isolates are rod-shaped, motile, catalase- and oxidase-positive, and aerobic. Growth occurred at 4–34 °C and at pH values of pH 5.5–8.0. Both strains showed strong β-haemolysis on blood agar. The major cellular polar lipids were phosphatidylethanolamine, phosphatidylglycerol and diphosphatidylglycerol. The dominant quinone was Q-9 (90 %), but noticeable amounts of Q-8 (9 %) and traces of Q-7 were also detected. Fatty acid profiles were typical for Pseudomonas species and exhibited C16 : 0 as a major component. Based on these results, we conclude that both strains belong to a novel species, for which the name Pseudomonas haemolytica sp. nov. is proposed. The type strain is WS 5063T (=DSM 108987T=LMG 31232T) and an additional strain is WS 5067 (=DSM 108988=LMG 31233).

Brevilactibacter flavus gen. nov., sp. nov., a novel bacterium of the family Propionibacteriaceae isolated from raw milk and dairy products and reclassification of Propioniciclava sinopodophylli as Brevilactibacter sinopodophylli comb. nov

Eight facultatively anaerobic rod-shaped bacteria were isolated from raw milk and two other dairy products. Results of phylogenetic analyses based on 16S rRNA gene sequences showed that the isolates are placed in a distinct lineage within the family Propionibacteriaceae with Propioniciclava sinopodophylli and Propioniciclava tarda as the closest relatives (94.6 and 93.5 % similarity, respectively). The cell-wall peptidoglycan contained meso-diaminopimelic acid, alanine and glutamic acid and was of the A1γ type (meso-DAP-direct). The major cellular fatty acid was anteiso-C15 : 0 and the major polar lipids were diphosphatidylglycerol, phosphatidyglycerol and three unidentified glycolipids. The quinone system contained predominantly menaquinone MK-9(H4). The G+C content of the genomic DNA of strain VG341T was 67.7 mol%. The whole-cell sugar pattern contained ribose, rhamnose, arabinose and galactose. On the basis of phenotypic and genetic data, eight strains (VG341T, WS4684, WS4769, WS 4882, WS4883, WS4901, WS4902 and WS4904) are proposed to be classified as members of a novel species in a new genus of the family Propionibacteriaceae , for which the name Brevilactibacter flavus gen. nov., sp. nov. is proposed. The type strain is VG341T (=WS4900T=DSM 100885T=LMG 29089T) and seven additional strains are WS4684, WS4769, WS4882, WS4883, WS4901, WS4902 and WS4904. Furthermore, we propose the reclassification of P. sinopodophylli as Brevilactibacter sinopodophylli comb. nov.

Draft Genome Sequence of Paenibacillus polymyxa DSM 292, a Gram-Positive, Spore-Forming Soil Bacterium with High Biotechnological Potential

Paenibacillus polymyxa DSM 292 was originally isolated from soil in 1947 due to its ability to produce antibiotics. The low proteolytic properties of strain DSM 292 warrant its examination as a host for heterologous protein production. Here, we report the draft genome sequence of DSM 292 as established by Illumina MiSeq paired-end sequencing.

Plasmid Reconstruction from Next-Gen Data: A Detailed Protocol for the Use of PLACNETw for the Reconstruction of Plasmids from WGS Datasets

Mobile Genetic Elements (MGE) play essential roles in adaptive bacterial evolution, facilitating genetic exchange for extrachromosomal DNA, especially antibiotic resistance genes and virulence factors. For this reason, high-throughput next-generation sequencing of bacteria is of great relevance, especially for clinical pathogenic bacteria. Accurate identification of MGE from whole-genome sequencing (WGS) datasets is one of the major challenges, still hindered by methodological limitations and high sequencing costs.This chapter encompasses the protocol used for plasmid reconstruction by applying the PLACNETw methodology, from raw reads to assembled plasmids and chromosome. PLACNETw is a graphical user-friendly interface to visualize and reconstruct MGE from short-read WGS datasets. No bioinformatic background or sophisticated computational resources are required and high precision and sensitivity are achieved.

Monitoring of Lactobacillus sanfranciscensis strains during wheat and rye sourdough fermentations by CRISPR locus length polymorphism PCR

Lactobacillus (L.) sanfranciscensis is a competitive key species in sourdough fermentations. However, the principles involved in establishing the commonly observed phenomenon of strain dominance are unresolved. This has been studied little because the methods for fast and reliable differentiation of strains and their monitoring during fermentation are tedious and cannot be done with large numbers of isolates. In this contribution, we present a strain-specific, PCR-based typing method that uses length heterogeneities of the clustered regularly interspaced short palindromic repeats (CRISPR) loci as they occur in the genomes of different strains. In silico analysis of 21 genomes revealed 14 different CRISPR genotypes. We then designed a primer set to simultaneously detect different strains in a multiplex PCR assay designated CRISPR locus length polymorphism PCR (CLLP-PCR). The usefulness of this method was evaluated in lab-scale sourdough fermentations conducted with rye and wheat flours. First, the flour was mixed with water to a dough yield of 200. Then each dough was inoculated with four different L. sanfranciscensis strains (TMW 1.1150, TMW 1.392, TMW 1.2142, and TMW 1.2138) at levels of 10⁹ cfu/g each. Sourdoughs were propagated at 28 °C for 5 days by back slopping 5% to the flour mass every 24 h. Samples were collected each day; DNA was isolated, and the presence of strains was detected qualitatively in the sourdoughs with PCR. L. sanfranciscensis TMW 1.392 became dominant as early as 2 days into the fermentation and remained the only detectable strain for the rest of the sampling period. CLLP-PCR proved to be useful in investigating the assertiveness of different strains of L. sanfranciscensis in sourdoughs. Therefore, CLLP-PCR may be used as a tool to investigate assertiveness of microorganisms in food fermentations at the strain level.

Application of different DNA extraction procedures, library preparation protocols and sequencing platforms: impact on sequencing results

In this study three DNA extraction procedures, two library preparation protocols and two sequencing platforms were applied to analyse six bacterial cultures and their corresponding DNA obtained as part of a proficiency test. The impact of each variable on sequencing results was assessed using the following parameters: reads quality, assembly and alignment statistics; number of single nucleotide polymorphisms (SNPs), detected applying assembly- and alignment-based strategies; antimicrobial resistance genes (ARGs), identified on de novo assemblies of all sequenced genomes. The investigated nucleic acid extraction procedures, library preparation kits and sequencing platforms do not significantly affect de novo assembly statistics and number of SNPs and ARGs. The only exception was observed for two duplicates, which were associated to one PCR-based library preparation kit. Results from this comparative study can support researchers in the choice toward the available pre-sequencing and sequencing options, and might suggest further comparisons to be performed.

The evolution and characterization of influenza A(H7N9) virus under the selective pressure of peramivir

Human infection with H7N9 virus has provoked global public health concern due to the substantial morbidity and mortality. Neuraminidase inhibitors (NAIs) are used as first-line drugs to treat the infection. However, virus quasispecies can evolve rapidly under drug pressure, which may alter various biological characteristics of virus. Using an in vitro evolution platform and next-generation sequencing, we found the presence of peramivir led to changes to the dominant populations of the virus. Two important amino acid substitutions were identified in NA, I222T and H274Y, which caused reduced susceptibilities to oseltamivir or both oseltamivir and peramivir as confirmed by enzyme- and cell-based assays. The NA-H274Y variant showed decreased replicative fitness at the early stage of infection accompanied with impaired NA function. The quasispecies evolution of H7N9 virus and the potential emergence of these two variants should be closely monitored, which may guide the adjustment of antiviral strategies.

In situ production and characterization of cloud forming dextrans in fruit-juices

Turbidity in beverages is typically achieved by addition of emulsion based cloud systems. Their intrinsic instability necessitates the widespread use of technological measures and use of food additives to prevent emulsion decay. In this work, we explored the possibility to establish a new generation of natural, stable clouding systems based on bacterial dextrans. Lactobacillus hordei TMW 1.1907 originating from water kefir was used to produce dextrans in sucrose supplemented apple or grape juices. By varying the fermentation conditions, two distinct types of dextran molecules could be produced at yields ranging from 2.5 to 8.5 g/L. The dextran-containing fermentates showed an unchanged turbidity after pasteurization at acidic pH and subsequent storage for three months. No sedimentation of particles occurred upon storage. Neutralization of the acidic fruit juices to pH 7 prior to fermentation significantly increased the dextran yields. The molecular weight, rms radii and turbidity of dextrans produced at 20 °C were higher than those produced at 30 °C. Characterization of the isolated dextrans by asymmetric flow field-flow fractionation coupled to multi-angle laser light scattering revealed a random-coil like structure and rms radii ranging from 66.0 to 87.4 nm. The averaged molar masses of the cloud forming dextrans were in the approximate range of 103.1 to 141.6 MDa. In conclusion, our results demonstrate the possibility to ferment fruit juices for in situ production of dextrans exhibiting novel techno-functional properties beyond gelling and thickening.

Detailed Structural Characterization of Glucans Produced by Glucansucrases from Leuconostoc citreum TMW 2.1194

The water kefir organism Leuconostoc citreum TMW 2.1194 forms highly branched dextrans with O3- and O4-bound side chains. To obtain detailed information on the enzymatic synthesis of these polymers, the four glucansucrases encoded by Leuconostoc citreum TMW 2.1194 were cloned, heterologously expressed, and used for polysaccharide production. Molecular and macromolecular structure of the synthesized glucans were analyzed by methylation analysis, two-dimensional NMR spectroscopy, oligosaccharide analysis after partial hydrolysis, and asymmetric flow field-flow fractionation. It was demonstrated that two glucansucrases form insoluble glucans with variously branched dextran sections and varying portions of consecutive, 1,3-linked glucose units. In contrast, the other two glucansucrases synthesized O3- (Lc6255) and O4-branched (Lc1785) soluble dextrans. Analysis, isolation, and characterization of enzymatically liberated oligosaccharides showed that monomeric and elongated side chains are abundant in both polysaccharides. From the structures and size distributions it was concluded that Lc1785 is mainly responsible for synthesis of fermentatively produced soluble dextrans.

Draft Genome Sequence of Anaerosphaera sp. Strain GS7-6-2, a Coccal Bacterium Isolated from a Biogas-Related Environment

Strain GS7-6-2 was isolated from a mesophilically operated biogas fermenter. The 16S rRNA gene sequence (93.27% identity to Anaerosphaera aminiphila WN036^T) indicated that GS7-6-2 represents a putative novel species within the genus Anaerosphaera (family Peptoniphilaceae).

Strain GS7-6-2 was isolated from a mesophilically operated biogas fermenter. The 16S rRNA gene sequence (93.27% identity to Anaerosphaera aminiphila WN036^T) indicated that GS7-6-2 represents a putative novel species within the genus Anaerosphaera (family Peptoniphilaceae). Here, we report the draft genome sequence of GS7-6-2 as established by Illumina paired-end sequencing.

Long fragments achieve lower base quality in Illumina paired-end sequencing

Illumina’s technology provides high quality reads of DNA fragments with error rates below 1/1000 per base. Sequencing runs typically generate millions of reads in which the vast majority of the reads has an average error rate below 1/1000. However, some paired-end sequencing data show the presence of a subpopulation of reads where the second read (R2) has lower average qualities. We show that the fragment length is a major driver of increased error rates in the R2 reads. Fragments above 500 nt tend to yield lower base qualities and higher error rates than shorter fragments. We use publicly available Illumina data to demonstrate that the fragment length dependency of the R2 read qualities exists in various library protocols, in different labs and using different sequencer models. Our finding extends the understanding of the Illumina read quality and has implications on error models for Illumina reads. It also sheds a light on the importance of controlling the fragment size during library preparation.

Sequence and cultivation study of Muribaculaceae reveals novel species, host preference, and functional potential of this yet undescribed family

BACKGROUND

Bacteria within family S24-7 (phylum Bacteroidetes) are dominant in the mouse gut microbiota and detected in the intestine of other animals. Because they had not been cultured until recently and the family classification is still ambiguous, interaction with their host was difficult to study and confusion still exists regarding sequence data annotation.

METHODS

We investigated family S24-7 by combining data from large-scale 16S rRNA gene analysis and from functional and taxonomic studies of metagenomic and cultured species.

RESULTS

A total of 685 species was inferred by full-length 16S rRNA gene sequence clustering. While many species could not be assigned ecological habitats (93,045 samples analyzed), the mouse was the most commonly identified host (average of 20% relative abundance and nine co-occurring species). Shotgun metagenomics allowed reconstruction of 59 molecular species, of which 34 were representative of the 16S rRNA gene-derived species clusters. In addition, cultivation efforts allowed isolating five strains representing three species, including two novel taxa. Genome analysis revealed that S24-7 spp. are functionally distinct from neighboring families and versatile with respect to complex carbohydrate degradation.

CONCLUSIONS

We provide novel data on the diversity, ecology, and description of bacterial family S24-7, for which the name Muribaculaceae is proposed.