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.

Whole-genome sequence of the Halobacillus naozhouensis type strain KACC 21980

We present the whole-genome sequence of Halobacillus naozhouensis Korean Agricultural Culture Collection (KACC) 21980T, isolated from China by Chen et al.. The genome of Halobacillus naozhouensis KACC 21980T comprises a circular chromosome (4.2 Mb) and one plasmid (17 kb). It includes a total of 4,168 predicted coding genes.

Illuminating the dark side of the human transcriptome with long read transcript sequencing

Background

The human transcriptome annotation is regarded as one of the most complete of any eukaryotic species. However, limitations in sequencing technologies have biased the annotation toward multi-exonic protein coding genes. Accurate high-throughput long read transcript sequencing can now provide additional evidence for rare transcripts and genes such as mono-exonic and non-coding genes that were previously either undetectable or impossible to differentiate from sequencing noise.

Results

We developed the Transcriptome Annotation by Modular Algorithms (TAMA) software to leverage the power of long read transcript sequencing and address the issues with current data processing pipelines. TAMA achieved high sensitivity and precision for gene and transcript model predictions in both reference guided and unguided approaches in our benchmark tests using simulated Pacific Biosciences (PacBio) and Nanopore sequencing data and real PacBio datasets. By analyzing PacBio Sequel II Iso-Seq sequencing data of the Universal Human Reference RNA (UHRR) using TAMA and other commonly used tools, we found that the convention of using alignment identity to measure error correction performance does not reflect actual gain in accuracy of predicted transcript models. In addition, inter-read error correction can cause major changes to read mapping, resulting in potentially over 6 K erroneous gene model predictions in the Iso-Seq based human genome annotation. Using TAMA’s genome assembly based error correction and gene feature evidence, we predicted 2566 putative novel non-coding genes and 1557 putative novel protein coding gene models.

Conclusions

Long read transcript sequencing data has the power to identify novel genes within the highly annotated human genome. The use of parameter tuning and extensive output information of the TAMA software package allows for in depth exploration of eukaryotic transcriptomes. We have found long read data based evidence for thousands of unannotated genes within the human genome. More development in sequencing library preparation and data processing are required for differentiating sequencing noise from real genes in long read RNA sequencing data.

Networking and Specificity-Changing DNA Methyltransferases in Helicobacter pylori

Epigenetic DNA base methylation plays important roles in gene expression regulation. We here describe a gene expression regulation network consisting of many DNA methyltransferases each frequently changing its target sequence-specificity. Our object Helicobacter pylori, a bacterium responsible for most incidence of stomach cancer, carries a large and variable repertoire of sequence-specific DNA methyltransferases. By creating a dozen of single-gene knockout strains for the methyltransferases, we revealed that they form a network controlling methylome, transcriptome and adaptive phenotype sets. The methyltransferases interact with each other in a hierarchical way, sometimes regulated positively by one methyltransferase but negatively with another. Motility, oxidative stress tolerance and DNA damage repair are likewise regulated by multiple methyltransferases. Their regulation sometimes involves translation start and stop codons suggesting coupling of methylation, transcription and translation. The methyltransferases frequently change their sequence-specificity through gene conversion of their target recognition domain and switch their target sets to remodel the network. The emerging picture of a metamorphosing gene regulation network, or firework, consisting of epigenetic systems ever-changing their specificity in search for adaptation, provides a new paradigm in understanding global gene regulation and adaptive evolution.

Complete Genomic Data of Burkholderia glumae Strain GX Associated with Bacterial Panicle Blight of Rice in China

Burkholderia glumae is a seedborne pathogen causing bacterial panicle blight of rice. Here, we report the complete genome of B. glumae strain GX, which represents the first whole-genome sequence of an isolate from China. The assembled genome consisted of five contigs, with two circular chromosomes of 3,712,850 and 2,750,046 bp and three plasmids of 201,571, 105,587, and 96,100 bp. This complete genome will provide a valuable resource for further studies on bacterial panicle blight worldwide.

Genome Sequence Resource for Ilyonectria mors-panacis, Causing Rusty Root Rot of Panax notoginseng

Ilyonectria mors-panacis is the cause of a serious disease hampering the production of Panax notoginseng, an important Chinese medicinal herb, widely used for its anti-inflammatory, antifatigue, hepato-protective, and coronary heart disease prevention effects. Here, we report the first Illumina-Pacbio hybrid sequenced draft genome assembly of I. mors-panacis strain G3B and its annotation. The availability of this genome sequence not only represents an important tool toward understanding the genetics behind the infection mechanism of I. mors-panacis strain G3B but also will help illuminate the complexities of the taxonomy of this species.

Rapid multi-locus sequence typing direct from uncorrected long reads using Krocus

Genome sequencing is rapidly being adopted in reference labs and hospitals for bacterial outbreak investigation and diagnostics where time is critical. Seven gene multi-locus sequence typing is a standard tool for broadly classifying samples into sequence types (STs), allowing, in many cases, to rule a sample out of an outbreak, or allowing for general characteristics about a bacterial strain to be inferred. Long-read sequencing technologies, such as from Oxford Nanopore, can produce read data within minutes of an experiment starting, unlike short-read sequencing technologies which require many hours/days. However, the error rates of raw uncorrected long read data are very high. We present Krocus which can predict a ST directly from uncorrected long reads, and which was designed to consume read data as it is produced, providing results in minutes. It is the only tool which can do this from uncorrected long reads. We tested Krocus on over 700 isolates sequenced using long-read sequencing technologies from Pacific Biosciences and Oxford Nanopore. It provides STs for isolates on average within 90 s, with a sensitivity of 94% and specificity of 97% on real sample data, directly from uncorrected raw sequence reads. The software is written in Python and is available under the open source license GNU GPL version 3.

FMLRC: Hybrid long read error correction using an FM-index

Background

Long read sequencing is changing the landscape of genomic research, especially de novo assembly. Despite the high error rate inherent to long read technologies, increased read lengths dramatically improve the continuity and accuracy of genome assemblies. However, the cost and throughput of these technologies limits their application to complex genomes. One solution is to decrease the cost and time to assemble novel genomes by leveraging “hybrid” assemblies that use long reads for scaffolding and short reads for accuracy.

Results

We describe a novel method leveraging a multi-string Burrows-Wheeler Transform with auxiliary FM-index to correct errors in long read sequences using a set of complementary short reads. We demonstrate that our method efficiently produces significantly more high quality corrected sequence than existing hybrid error-correction methods. We also show that our method produces more contiguous assemblies, in many cases, than existing state-of-the-art hybrid and long-read only de novo assembly methods.

Conclusion

Our method accurately corrects long read sequence data using complementary short reads. We demonstrate higher total throughput of corrected long reads and a corresponding increase in contiguity of the resulting de novo assemblies. Improved throughput and computational efficiency than existing methods will help better economically utilize emerging long read sequencing technologies.

Comparative Whole-Genomic Analysis of an Ancient L2 Lineage Mycobacterium tuberculosis Reveals a Novel Phylogenetic Clade and Common Genetic Determinants of Hypervirulent Strains

Background: Development of improved therapeutics against tuberculosis (TB) is hindered by an inadequate understanding of the relationship between disease severity and genetic diversity of its causative agent, Mycobacterium tuberculosis. We previously isolated a hypervirulent M. tuberculosis strain H112 from an HIV-negative patient with an aggressive disease progression from pulmonary TB to tuberculous meningitis—the most severe manifestation of tuberculosis. Human macrophage challenge experiment demonstrated that the strain H112 exhibited significantly better intracellular survivability and induced lower level of TNF-α than the reference virulent strain H37Rv and other 123 clinical isolates.

Aim: The present study aimed to identify the potential genetic determinants of mycobacterial virulence that were common to strain H112 and hypervirulent M. tuberculosis strains of the same phylogenetic clade isolated in other global regions.

Methods: A low-virulent M. tuberculosis strain H54 which belonged to the same phylogenetic lineage (L2) as strain H112 was selected from a collection of 115 clinical isolates. Both H112 and H54 were whole-genome-sequenced using PacBio sequencing technology. A comparative genomics approach was adopted to identify mutations present in strain H112 but absent in strain H54. Subsequently, an extensive phylogenetic analysis was conducted by including all publically available M. tuberculosis genomes. Single-nucleotide-polymorphisms (SNPs) and structural variations (SVs) common to hypervirulent strains in the global collection of genomes were considered as potential genetic determinants of hypervirulence.

Results:Sequencing data revealed that both H112 and H54 were identified as members of the same sub-lineage L2.2.1. After excluding the lineage-related mutations shared between H112 and H54, we analyzed the phylogenetic relatedness of H112 with global collection of M. tuberculosis genomes (n = 4,338), and identified a novel phylogenetic clade in which four hypervirulent strains isolated from geographically diverse regions were clustered together. All hypervirulent strains in the clade shared 12 SNPs and 5 SVs with H112, including those affecting key virulence-associated loci, notably, a deleterious SNP (rv0178 p. D150E) within mce1 operon and an intergenic deletion (854259_ 854261delCC) in close-proximity to phoP.

Conclusion: The present study identified common genetic factors in a novel phylogenetic clade of hypervirulent M. tuberculosis. The causative role of these mutations in mycobacterial virulence should be validated in future study.

Complete genome sequence of Stenotrophomonas sp. KACC 91585, an efficient bacterium for unsaturated fatty acid hydration

Hydroxy fatty acids (HFAs) such as 10-hydroxystearic acid (10-HSA) and 10-hydroxy-12(Z)-octadecenoic acid (10-HOD), which are similar to ricinoleic acid, are important starting materials and intermediates for the industrial manufacture of many commodities. Stenotrophomonas sp. KACC 91585, which was isolated from lake sediment, is an efficient bacterium for unsaturated fatty acid hydration that produces 10-HSA and 10-HOD from oleic acid and linoleic acid, respectively, with high conversion rates. The complete genome of this strain is 4,541,729bp with 63.83% GC content and devoid of plasmids. Sets of genes involved in the fatty acid biosynthesis and modification as well as modified lipids were identified in the genome, and these genes were concerned with HFA production. This genome sequence provides molecular information and elucidation for HFA production, and will be used as an efficient biocatalyst source for the biotechnological production of HFA.

Whole genome sequence of Pantoea ananatis R100, an antagonistic bacterium isolated from rice seed

Pantoea ananatis is a group of bacteria, which was first reported as plant pathogen. Recently, several papers also described its biocontrol ability. In 2003, P. ananatis R100, which showed strong antagonism against several plant pathogens, was isolated from rice seeds. In this study, whole genome sequence of this strain was determined by SMRT Cell technology. The total genome size of R100 is 4,857,861bp with 4659 coding genes (CDS), 82 tRNAs and 22 rRNAs. The genome sequence of R100 may shed a light on the research of antagonism P. ananatis.

Complete genome sequence of Lactobacillus heilongjiangensis DSM 28069(T): Insight into its probiotic potential

Lactobacillus heilongjiangensis DSM 28069(T) is a potential probiotic isolated from traditional Chinese pickle. Here we report the complete genome sequence of this strain. The complete genome is 2,790,548bp with the GC content of 37.5% and devoid of plasmids. Sets of genes involved in the biosynthesis of riboflavin and folate were identified in the genome, which revealed its potential application in biotechnological industry. The genome sequence of L. heilongjiangensis DSM 28069(T) now provides the fundamental information for future studies.

Whole genome sequence of Pseudomonas aeruginosa F9676, an antagonistic bacterium isolated from rice seed

Pseudomonas aeruginosa is a group of bacteria, which can be isolated from diverse ecological niches. P. aeruginosa strain F9676 was first isolated from a rice seed sample in 2003. It showed strong antagonism against several plant pathogens. In this study, whole genome sequencing was carried out. The total genome size of F9676 is 6368,008bp with 5586 coding genes (CDS), 67 tRNAs and 3 rRNAs. The genome sequence of F9676 may shed a light on antagonism P. aeruginosa.

Full genome sequence of Brevibacillus laterosporus strain B9, a biological control strain isolated from Zhejiang, China

Brevibacillus laterosporus was newly classified from Bacillus laterosporus, which has ability to be used as a biological control agent in crop field. B. laterosporus strain B9 is an aerobic, motile, Gram-positive, spore-forming rod that was isolated from a field of Oryza sativa in Zhejiang, China in 2011. This bacterium has been confirmed to be a strong antagonist against bacterial brown strip of rice caused by Acidovorex avenae subsp. avenae. Here we describe the features of B. laterosporus strain B9, together with the complete genome sequence and its annotation. The 5,272,435bp genome contains 4804 protein-coding genes and 227 RNA-only encoding genes with 2 plasmids.

Genome sequence of Xanthomonas sacchari R1, a biocontrol bacterium isolated from the rice seed

Xanthomonas sacchari, was first identified as a pathogenic bacterium isolated from diseased sugarcane in Guadeloupe. In this study, R1 was first isolated from rice seed samples from Philippines in 2002. The antagonistic ability against several rice pathogens raises our attention. The genomic feature of this strain was described in this paper. The total genome size of X. sacchari R1 is 5,000,479 bp with 4315 coding sequences (CDS), 59 tRNAs, 2rRNAs and one plasmid.