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

Visualizing metagenomic and metatranscriptomic data: A comprehensive review

The fields of Metagenomics and Metatranscriptomics involve the examination of complete nucleotide sequences, gene identification, and analysis of potential biological functions within diverse organisms or environmental samples. Despite the vast opportunities for discovery in metagenomics, the sheer volume and complexity of sequence data often present challenges in processing analysis and visualization. This article highlights the critical role of advanced visualization tools in enabling effective exploration, querying, and analysis of these complex datasets. Emphasizing the importance of accessibility, the article categorizes various visualizers based on their intended applications and highlights their utility in empowering bioinformaticians and non-bioinformaticians to interpret and derive insights from meta-omics data effectively.

Sphingobacterium oryzagri sp. nov., isolated from rice paddy soil

An aerobic, Gram-stain-negative and short rod-shaped bacterial strain, designated M6-31T, was isolated from rice paddy soil sampled in Miryang, Republic of Korea. Growth was observed at 4-35 °C (optimum, 28 °C), pH 6.0-9.0 (optimum, pH 7.0-8.0) and in the presence of 0-4 % (w/v) NaCl (optimum, 0 % w/v). Phylogenetic analysis based on 16S rRNA gene sequences grouped strain M6-31T with Sphingobacterium bambusae IBFC2009T, Sphingobacterium griseoflavum SCU-B140T and Sphingobacterium solani MLS-26-JM13-11T in the same clade, with the 16S rRNA gene sequence similarities ranging from 95.8 to 96.6 %. A genome-based phylogenetic tree reconstructed by using all publicly available Sphingobacterium genomes placed strain M6-31T with S. bambusae KACC 22910T, 'Sphingobacterium deserti' ACCC 05744T, S. griseoflavum CGMCC 1.12966T and Sphingobacterium paludis CGMCC 1.12801T. Orthologous average nucleotide identity and digital DNA-DNA hybridization values between strain M6-31T and its closely related strains were lower than 74.6 and 22.0 %, respectively. The respiratory quinone was menaquinone-7, and the major polar lipid was phosphatidylethanolamine. The major fatty acids (>10 %) were C15 : 0 iso, C17 : 0 iso 3OH and summed feature 3. The phenotypic, chemotaxonomic and genotypic data obtained in this study showed that strain M6-31T represents a novel species of the genus Sphingobacterium, for which the name Sphingobacterium oryzagri sp. nov. (type strain M6-31T=KACC 22765T=JCM 35893T) is proposed.

Unraveling metagenomics through long-read sequencing: a comprehensive review

The study of microbial communities has undergone significant advancements, starting from the initial use of 16S rRNA sequencing to the adoption of shotgun metagenomics. However, a new era has emerged with the advent of long-read sequencing (LRS), which offers substantial improvements over its predecessor, short-read sequencing (SRS). LRS produces reads that are several kilobases long, enabling researchers to obtain more complete and contiguous genomic information, characterize structural variations, and study epigenetic modifications. The current leaders in LRS technologies are Pacific Biotechnologies (PacBio) and Oxford Nanopore Technologies (ONT), each offering a distinct set of advantages. This review covers the workflow of long-read metagenomics sequencing, including sample preparation (sample collection, sample extraction, and library preparation), sequencing, processing (quality control, assembly, and binning), and analysis (taxonomic annotation and functional annotation). Each section provides a concise outline of the key concept of the methodology, presenting the original concept as well as how it is challenged or modified in the context of LRS. Additionally, the section introduces a range of tools that are compatible with LRS and can be utilized to execute the LRS process. This review aims to present the workflow of metagenomics, highlight the transformative impact of LRS, and provide researchers with a selection of tools suitable for this task.

Brevibacillus ruminantium sp. nov., isolated from cow faeces

An aerobic, Gram-stain-positive, rod-shaped, endospore-forming bacterial strain, designated BB3-R1T, was isolated from cow faeces sampled in Daejeon, Republic of Korea. Growth was observed at 25-45 °C (optimum, 35-40 °C) and pH 7.0-9.0 (optimum, pH 8.0), with up to 3 % (w/v) NaCl (optimum, 0 % NaCl). blast analysis of 16S rRNA gene sequences revealed the highest sequence similarity of strain BB3-R1T to Brevibacillus borstelensis NRRL NRS-818T (98.8 %) followed by Brevibacillus panacihumi JCM 15085T (97.5 %). According to 16S rRNA gene and whole-genome based phylogenetic trees, strain BB3-R1T clustered with Brevibacillus composti FJAT-54423T and B. borstelensis NRRL NRS-818T. OrthoANI and dDDH values of strain BB3-R1T with the closely related strains were lower than 77.5 and 26.8 %, respectively. The major menaquinones and polar lipids of the strain were MK-7 and phosphatidylmonomethylethanolamine, diphosphatidylglycerol, phosphatidylglycerol and phosphatidylethanolamine, respectively. The major fatty acids (>10 %) were C14 : 0 iso, C15 : 0 iso, C15 : 0 anteiso and C16 : 1  ω7c alcohol. The cell-wall peptidoglycan contained cross-linked meso-diaminopimelic acid (type A1 gamma). The phenotypic, chemotaxonomic and genotypic data obtained in this study showed that the strain represents a novel species of the genus Brevibacillus, for which the name Brevibacillus ruminantium sp. nov. (type strain BB3-R1T=KACC 22663T=NBRC 115962T) is proposed.

Halobacillus salinarum sp. nov., Halobacillus shinanisalinarum sp. nov. and Halobacillus amylolyticus sp. nov., isolated from saltern soil

Three bacterial strains, designated SSBR10-3T, SSTM10-2T and SSHM10-5T, were isolated from saltern soil sampled in Jeollanam-do, Republic of Korea. Cells were aerobic, Gram-stain-positive, flagellated and rod-shaped. The strains grew optimally at 28°C and at pH 7.0. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strains SSBR10-3T, SSTM10-2T and SSHM10-5T were placed within the genus Halobacillus, showing the highest similarity to Halobacillus alkaliphilus FP5T (98.6 %), 'Halobacillus ihumii' Marseille-Q1234T (98.5 %) and Halobacillus locisalis MSS-155T (98.6 %), respectively. The genomic similarity values between strains SSBR10-3T, SSTM10-2T and SSHM10-5T and their related species were 17.6-22.6 % for digital DNA-DNA hybridization (dDDH) and 69.6-78.5 % for orthologous average nucleotide identity (OrthoANI), which were lower than the thresholds recommended for species delineation. The dDDH and OrthoANI values among the three strains were below 38.3 and 89.4 %, respectively. Besides the differences in genomic features, strains SSBR10-3T, SSTM10-2T and SSHM10-5T were distinct from each other and from members of the genus in terms of phenotypic traits related to substrate assimilation. The cell-wall peptidoglycan contained meso-diaminopimelic acid, the major fatty acids were anteiso-C15 : 0, iso-C16 : 0 and anteiso-C17 : 0, and the predominant menaquinone was MK-7 for all three strains. Diphosphatidylglycerol, phosphatidylglycerol and an unidentified phospholipid were present in their polar lipid profiles. Based on a polyphasic approach incorporating genomic data, strains SSBR10-3T, SSTM10-2T and SSHM10-5T represent novel species, for which the names Halobacillus salinarum sp. nov. (SSBR10-3T=DSM 114353T=KACC 21935T=NBRC 115504T), Halobacillus shinanisalinarum sp. nov. (SSTM10-2T=DSM 114354T=KACC 21936T=NBRC 115505T) and Halobacillus amylolyticus sp. nov. (SSHM10-5T=DSM 114355T= KACC 21937T=NBRC 115506T) are proposed.

Amygdalobacter indicium gen. nov., sp. nov., and Amygdalobacter nucleatus sp. nov., gen. nov.: novel bacteria from the family Oscillospiraceae isolated from the female genital tract

Four obligately anaerobic Gram-positive bacteria representing one novel genus and two novel species were isolated from the female genital tract. Both novel species, designated UPII 610-JT and KA00274T, and an additional isolate of each species were characterized utilizing biochemical, genotypic and phylogenetic analyses. All strains were non-motile and non-spore forming, asaccharolytic, non-cellulolytic and indole-negative coccobacilli. Fatty acid methyl ester analysis for UPII 610-JT and KA00274T and additional isolates revealed C16 : 0, C18 : 0, C18:1ω9c and C18:2ω6,9c to be the major fatty acids for both species. UPII 610-JT had a 16S rRNA gene sequence similarity of 99.4 % to an uncultured clone sequence (AY724740) designated as Bacterial Vaginosis Associated Bacterium 2 (BVAB2). KA00274T had a 16S rRNA gene sequence similarity of 96.5 % to UPII 610-JT. Whole genomic DNA mol% G+C content was 42.2 and 39.3 % for UPII 610-JT and KA00274T, respectively. Phylogenetic analyses indicate these isolates represent a novel genus and two novel species within the Oscillospiraceae family. We propose the names Amygdalobacter indicium gen. nov., sp. nov., for UPII 610-JT representing the type strain of this species (=DSM 112989T, =ATCC TSD-274T) and Amygdalobacter nucleatus gen. nov., sp. nov., for KA00274T representing the type strain of this species (=DSM 112988T, =ATCC TSD-275T).

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai, Drosophila takahashii, Drosophila bipectinata, and Drosophila ananassae

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma. However, the assembled portion of the Muller F element (orthologous to the fourth chromosome in Drosophila melanogaster) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to >18 Mb. Here, we present chromosome-level long-read genome assemblies for 4 Drosophila species with expanded F elements ranging in size from 2.3 to 20.5 Mb. Each Muller element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai , Drosophila takahashii , Drosophila bipectinata , and Drosophila ananassae

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma . However, the assembled portion of the Muller F Element (orthologous to the fourth chromosome in Drosophila melanogaster ) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to > 18 Mb. Here, we present chromosome-level long read genome assemblies for four Drosophila species with expanded F Elements ranging in size from 2.3 Mb to 20.5 Mb. Each Muller Element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

Clinical Diagnostics of Bacterial Infections and Their Resistance to Antibiotics-Current State and Whole Genome Sequencing Implementation Perspectives

Antimicrobial resistance (AMR), defined as the ability of microorganisms to withstand antimicrobial treatment, is responsible for millions of deaths annually. The rapid spread of AMR across continents warrants systematic changes in healthcare routines and protocols. One of the fundamental issues with AMR spread is the lack of rapid diagnostic tools for pathogen identification and AMR detection. Resistance profile identification often depends on pathogen culturing and thus may last up to several days. This contributes to the misuse of antibiotics for viral infection, the use of inappropriate antibiotics, the overuse of broad-spectrum antibiotics, or delayed infection treatment. Current DNA sequencing technologies offer the potential to develop rapid infection and AMR diagnostic tools that can provide information in a few hours rather than days. However, these techniques commonly require advanced bioinformatics knowledge and, at present, are not suited for routine lab use. In this review, we give an overview of the AMR burden on healthcare, describe current pathogen identification and AMR screening methods, and provide perspectives on how DNA sequencing may be used for rapid diagnostics. Additionally, we discuss the common steps used for DNA data analysis, currently available pipelines, and tools for analysis. Direct, culture-independent sequencing has the potential to complement current culture-based methods in routine clinical settings. However, there is a need for a minimum set of standards in terms of evaluating the results generated. Additionally, we discuss the use of machine learning algorithms regarding pathogen phenotype detection (resistance/susceptibility to an antibiotic).

Bioinformatic Tools for NGS-Based Metagenomics to Improve the Clinical Diagnosis of Emerging, Re-Emerging and New Viruses

Epidemics and pandemics have occurred since the beginning of time, resulting in millions of deaths. Many such disease outbreaks are caused by viruses. Some viruses, particularly RNA viruses, are characterized by their high genetic variability, and this can affect certain phenotypic features: tropism, antigenicity, and susceptibility to antiviral drugs, vaccines, and the host immune response. The best strategy to face the emergence of new infectious genomes is prompt identification. However, currently available diagnostic tests are often limited for detecting new agents. High-throughput next-generation sequencing technologies based on metagenomics may be the solution to detect new infectious genomes and properly diagnose certain diseases. Metagenomic techniques enable the identification and characterization of disease-causing agents, but they require a large amount of genetic material and involve complex bioinformatic analyses. A wide variety of analytical tools can be used in the quality control and pre-processing of metagenomic data, filtering of untargeted sequences, assembly and quality control of reads, and taxonomic profiling of sequences to identify new viruses and ones that have been sequenced and uploaded to dedicated databases. Although there have been huge advances in the field of metagenomics, there is still a lack of consensus about which of the various approaches should be used for specific data analysis tasks. In this review, we provide some background on the study of viral infections, describe the contribution of metagenomics to this field, and place special emphasis on the bioinformatic tools (with their capabilities and limitations) available for use in metagenomic analyses of viral pathogens.

Complete mitochondrial genome of the Spanish toothcarp, Aphanius iberus (Valenciennes, 1846) (Actinopterygii, Aphaniidae) and its phylogenetic position within the Cyprinodontiformes order

BACKGROUND

The Spanish toothcarp (Aphanius iberus Valenciennes, 1846) is a small fish endemic to the eastern coastline of the Iberian Peninsula and is currently listed as "Endangered" (category IUCN: EN). It mainly inhabits brackish waters which can exhibit large fluctuations in temperature and salinity throughout the year. The genetics of A. iberus are not well-known since most studies have only evaluated the genetic structure of the species under a conservation framework in order to identify its potential conservation units. Different phylogenetic relationships of Aphanius have been published based on some particular genes. In the present study, the entire mitochondrial genome of A. iberus was obtained for the first time in the context of an A. iberus reference genome and a hypothesis regarding its phylogenetic position was considered.

METHODS AND RESULTS

The mitogenome (a circular doble-stranded DNA sequence of 16,708 bp) was reconstructed and aligned against 83 Cyprinodontiformes and two outgroup taxa to identify the phylogenetic position of A. iberus. PartitionFinder was first used to test for the best evolutionary model and the phylogenetic analyses were performed using two methods: Maximun-Likelihood Approximation (IQ-Tree) and Bayesian inference (MrBayes). Our results show that A. iberus forms a sister group with Orestias ascotanensis, a cyprinodontiform species native to South America.

CONCLUSIONS

The results were congruent with the traditional morphometric reconstructed trees and with a geological vicariant hypothesis involving Cyprinodontiformes where Aphaniidae is shown as a monophyletic family separated from the family Cyprinodontidae. The information gathered from this study is not only valuable for improving our understanding of the evolutionary history of A. iberus, but for future genomic studies involving the species.

DNA Methyltransferase Regulates Nitric Oxide Homeostasis and Virulence in a Chronically Adapted Pseudomonas aeruginosa Strain

Opportunistic pathogens such as Pseudomonas aeruginosa adapt their genomes rapidly during chronic infections. Understanding their epigenetic regulation may provide biomarkers for diagnosis and reveal novel regulatory mechanisms. We performed single-molecule real-time sequencing (SMRT-seq) to characterize the methylome of a chronically adapted P. aeruginosa clinical strain, TBCF10839. Two N⁶-methyladenine (6mA) methylation recognition motifs (RCCANNNNNNNTGAR and TRGANNNNNNTGC [modification sites are in bold]) were identified and predicted as new type I methylation sites using REBASE analysis. We confirmed that the motif TRGANNNNNNTGC was methylated by the methyltransferase (MTase) M.PaeTBCFII, according to methylation sensitivity assays in vivo and vitro. Transcriptomic analysis showed that a ΔpaeTBCFIIM knockout mutant significantly downregulated nitric oxide reductase (NOR) regulation and expression of coding genes such as nosR and norB, which contain methylated motifs in their promoters or coding regions. The ΔpaeTBCFIIM strain exhibited reduced intercellular survival capacity in NO-producing RAW264.7 macrophages and attenuated virulence in a Galleria mellonella infection model; the complemented strain recovered these defective phenotypes. Further phylogenetic analysis demonstrated that homologs of M.PaeTBCFII occur frequently in P. aeruginosa as well as other bacterial species. Our work therefore provided new insights into the relationship between DNA methylation, NO detoxification, and bacterial virulence, laying a foundation for further exploring the molecular mechanism of DNA methyltransferase in regulating the pathogenicity of P. aeruginosa. IMPORTANCE Pseudomonas aeruginosa is an opportunistic pathogen which causes acute and chronic infections that are difficult to treat. Comparative genomic analysis has showed broad genome diversity among P. aeruginosa clinical strains and revealed their different regulatory traits compared to the laboratory strains. While current investigation of the epigenetics of P. aeruginosa is still lacking, understanding epigenetic regulation may provide biomarkers for diagnosis and facilitate development of novel therapies. Denitrification capability is critical for microbial versatility in response to different environmental stress conditions, including the bacterial infection process, where nitric oxide (NO) can be generated by phagocytic cells. The denitrification regulation mechanisms have been studied intensively at genetic and biochemical levels. However, there is very little evidence about the epigenetic regulation of bacterial denitrification mechanism. P. aeruginosa TBCF10839 is a chronically host-adapted strain isolated from a cystic fibrosis (CF) patient with special antiphagocytosis characteristics. Here, we investigated the regulatory effect of an orphan DNA MTase, M.PaeTBCFII, in P. aeruginosa TBCF10839. We demonstrated that the DNA MTase regulates the transcription of denitrification genes represented by NOR and affects antiphagocytic ability in bacteria. In silico analysis suggested that DNA methylation modification may enhance gene expression by affecting the binding of transacting factors such as DNR and RpoN. Our findings not only deepen the understanding of the role of DNA MTase in transcriptional regulation in P. aeruginosa but also provide a theoretical foundation for the in-depth study of the molecular mechanism of the epigenetic regulation on denitrification, virulence, and host-pathogen interaction.

A De Novo Chromosome-Level Genome Assembly of the White-Tailed Deer, Odocoileus Virginianus

Cervids are distinguished by the shedding and regrowth of antlers. Furthermore, they provide insights into prion and other diseases. Genomic resources can facilitate studies of the genetic underpinnings of deer phenotypes, behavior, and disease resistance. Widely distributed in North America, the white-tailed deer (Odocoileus virginianus) has recreational, commercial, and food source value for many households. We present a genome generated using DNA from a single Illinois white-tailed sequenced on the PacBio Sequel II platform and assembled using Wtdbg2. Omni-C chromatin conformation capture sequencing was used to scaffold the genome contigs. The final assembly was 2.42 Gb, consisting of 508 scaffolds with a contig N50 of 21.7 Mb, a scaffold N50 of 52.4 Mb, and a BUSCO complete score of 93.1%. Thirty-six chromosome pseudomolecules comprised 93% of the entire sequenced genome length. A total of 20 651 predicted genes using the BRAKER pipeline were validated using InterProScan. Chromosome length assembly sequences were aligned to the genomes of related species to reveal corresponding chromosomes.

Comprehensive transcriptome characterization of Grus japonensis using PacBio SMRT and Illumina sequencing

The red-crowned crane (Grus japonensis) is an endangered species distributed across southeast Russia, northeast China, Korea, and Japan. Here, we sequenced for the first time the full-length unreferenced transcriptome of red-crowned crane mixed samples using a PacBio Sequel platform. A total of 359,136 circular consensus sequences (CCS) were obtained via clustering to remove redundancy. A total of 303,544 full-length non-chimeric sequences were identified by judging whether CCS contained 5' and 3' adapters, and the poly(A) tail. Eight samples were sequenced using Illumina, and PacBio sequencing data were corrected according to the collected Illumina data to obtain more accurate full-length transcripts. A total of 4,100 long non-coding RNAs, 13,115 simple sequences repeat loci and 29 transcription factor families were identified. The expression of lncRNAs and TFs in pancreas was lowest comparing with other tissues. Many enriched immune-related transmission pathways (MHC and IL receptors) were identified in the spleen. This study will contribute to a better understanding of the gene structure and post-transcriptional regulatory network, and provide references for future studies on red-crowned cranes.

A review of computational tools for generating metagenome-assembled genomes from metagenomic sequencing data

Metagenomic sequencing provides a culture-independent avenue to investigate the complex microbial communities by constructing metagenome-assembled genomes (MAGs). A MAG represents a microbial genome by a group of sequences from genome assembly with similar characteristics. It enables us to identify novel species and understand their potential functions in a dynamic ecosystem. Many computational tools have been developed to construct and annotate MAGs from metagenomic sequencing, however, there is a prominent gap to comprehensively introduce their background and practical performance. In this paper, we have thoroughly investigated the computational tools designed for both upstream and downstream analyses, including metagenome assembly, metagenome binning, gene prediction, functional annotation, taxonomic classification, and profiling. We have categorized the commonly used tools into unique groups based on their functional background and introduced the underlying core algorithms and associated information to demonstrate a comparative outlook. Furthermore, we have emphasized the computational requisition and offered guidance to the users to select the most efficient tools. Finally, we have indicated current limitations, potential solutions, and future perspectives for further improving the tools of MAG construction and annotation. We believe that our work provides a consolidated resource for the current stage of MAG studies and shed light on the future development of more effective MAG analysis tools on metagenomic sequencing.

PB-Motif-A Method for Identifying Gene/Pseudogene Rearrangements With Long Reads: An Application to CYP21A2 Genotyping

Long read sequencing technologies have the potential to accurately detect and phase variation in genomic regions that are difficult to fully characterize with conventional short read methods. These difficult to sequence regions include several clinically relevant genes with highly homologous pseudogenes, many of which are prone to gene conversions or other types of complex structural rearrangements. We present PB-Motif, a new method for identifying rearrangements between two highly homologous genomic regions using PacBio long reads. PB-Motif leverages clustering and filtering techniques to efficiently report rearrangements in the presence of sequencing errors and other systematic artifacts. Supporting reads for each high-confidence rearrangement can then be used for copy number estimation and phased variant calling. First, we demonstrate PB-Motif's accuracy with simulated sequence rearrangements of PMS2 and its pseudogene PMS2CL using simulated reads sweeping over a range of sequencing error rates. We then apply PB-Motif to 26 clinical samples, characterizing CYP21A2 and its pseudogene CYP21A1P as part of a diagnostic assay for congenital adrenal hyperplasia. We successfully identify damaging variation and patient carrier status concordant with clinical diagnosis obtained from multiplex ligation-dependent amplification (MLPA) and Sanger sequencing. The source code is available at: github.com/zstephens/pb-motif.

The first draft genome of feather grasses using SMRT sequencing and its implications in molecular studies of Stipa

The Eurasian plant Stipa capillata is the most widespread species within feather grasses. Many taxa of the genus are dominants in steppe plant communities and can be used for their classification and in studies related to climate change. Moreover, some species are of economic importance mainly as fodder plants and can be used for soil remediation processes. Although large-scale molecular data has begun to appear, there is still no complete or draft genome for any Stipa species. Thus, here we present a single-molecule long-read sequencing dataset generated using the Pacific Biosciences Sequel System. A draft genome of about 1004 Mb was obtained with a contig N50 length of 351 kb. Importantly, here we report 81,224 annotated protein-coding genes, present 77,614 perfect and 58 unique imperfect SSRs, reveal the putative allopolyploid nature of S. capillata, investigate the evolutionary history of the genus, demonstrate structural heteroplasmy of the chloroplast genome and announce for the first time the mitochondrial genome in Stipa. The assembled nuclear, mitochondrial and chloroplast genomes provide a significant source of genetic data for further works on phylogeny, hybridisation and population studies within Stipa and the grass family Poaceae.