Short Paired-end Reads Research Articles

MetaAll: integrative bioinformatics workflow for analysing clinical metagenomic data.

Over the past decade, there have been many improvements in the field of metagenomics, including sequencing technologies, advances in bioinformatics and the development of reference databases, but a one-size-fits-all sequencing and bioinformatics pipeline does not yet seem achievable. In this study, we address the bioinformatics part of the analysis by combining three methods into a three-step workflow that increases the sensitivity and specificity of clinical metagenomics and improves pathogen detection. The individual tools are combined into a user-friendly workflow suitable for analysing short paired-end (PE) and long reads from metagenomics datasets-MetaAll. To demonstrate the applicability of the developed workflow, four complicated clinical cases with different disease presentations and multiple samples collected from different biological sites as well as the CAMI Clinical pathogen detection challenge dataset were used. MetaAll was able to identify putative pathogens in all but one case. In this case, however, traditional microbiological diagnostics were also unsuccessful. In addition, co-infection with Haemophilus influenzae and Human rhinovirus C54 was detected in case 1 and co-infection with SARS-Cov-2 and Influenza A virus (FluA) subtype H3N2 was detected in case 3. In case 2, in which conventional diagnostics could not find a pathogen, mNGS pointed to Klebsiella pneumoniae as the suspected pathogen. Finally, this study demonstrated the importance of combining read classification, contig validation and targeted reference mapping for more reliable detection of infectious agents in clinical metagenome samples.

Easy analysis of bacterial whole-genome sequencing data for clinical microbiologists using open-source Galaxy platform: Characterization of ESBL-producing Enterobacterales from bloodstream infections

ObjectivesClinical microbiologists require easy-to-use open access tools with graphical interfaces to perform bacterial whole-genome sequencing (WGS) in routine practice. This study aimed to build a bioinformatics pipeline on the open-source Galaxy platform, facilitating comprehensive and reproducible analysis of bacterial WGS data in a few steps. We then used it to characterize our local epidemiology of ESBL-producing Enterobacterales isolated from patients with bacteremia. MethodsWe built a bioinformatics pipeline consisting of the following sequential tools: Fastp (input data trimming); FastQC (read quality control); SPAdes (genome assembly); Quast (quality control of genome assembly); Prokka (gene annotation); Staramr (ResFinder database) and ABRicate (CARD database) for antimicrobial resistance (AMR) gene screening and molecular strain typing. Paired-end short read WGS data from all ESBL-producing Enterobacterales strains isolated from patients with bacteremia over one year were analysed. ResultsThe Galaxy platform does not require command line tools. The bioinformatics pipeline was constructed within one hour. It only required uploading fastq files and facilitated systematization of de novo assembly of genomes, MLST typing, and AMR gene screening in one step. Among the 66 ESBL-producing strains analysed, the two most frequent ESBL genes were blaCTX−M-15 (62.1%) and blaCTX−M-27 (13.6%). ConclusionsThe open-access Galaxy platform provides a graphical interface and easy-to-use tools suitable for routine use in clinical microbiology laboratories without bioinformatics specialists. We believe that this platform will facilitate fast and low-cost analysis of bacterial WGS data, especially in resource-limited settings.

CoRAL accurately resolves extrachromosomal DNA genome structures with long-read sequencing.

Extrachromosomal DNA (ecDNA) is a central mechanism for focal oncogene amplification in cancer, occurring in ∼15% of early-stage cancers and ∼30% of late-stage cancers. ecDNAs drive tumor formation, evolution, and drug resistance by dynamically modulating oncogene copy number and rewiring gene-regulatory networks. Elucidating the genomic architecture of ecDNA amplifications is critical for understanding tumor pathology and developing more effective therapies. Paired-end short-read (Illumina) sequencing and mapping have been utilized to represent ecDNA amplifications using a breakpoint graph, in which the inferred architecture of ecDNA is encoded as a cycle in the graph. Traversals of breakpoint graphs have been used to successfully predict ecDNA presence in cancer samples. However, short-read technologies are intrinsically limited in the identification of breakpoints, phasing together complex rearrangements and internal duplications, and deconvolution of cell-to-cell heterogeneity of ecDNA structures. Long-read technologies, such as from Oxford Nanopore Technologies, have the potential to improve inference as the longer reads are better at mapping structural variants and are more likely to span rearranged or duplicated regions. Here, we propose Complete Reconstruction of Amplifications with Long reads (CoRAL) for reconstructing ecDNA architectures using long-read data. CoRAL reconstructs likely cyclic architectures using quadratic programming that simultaneously optimizes parsimony of reconstruction, explained copy number, and consistency of long-read mapping. CoRAL substantially improves reconstructions in extensive simulations and 10 data sets from previously characterized cell lines compared with previous short- and long-read-based tools. As long-read usage becomes widespread, we anticipate that CoRAL will be a valuable tool for profiling the landscape and evolution of focal amplifications in tumors.

ProcaryaSV: structural variation detection pipeline for bacterial genomes using short-read sequencing

BackgroundStructural variations play an important role in bacterial genomes. They can mediate genome adaptation quickly in response to the external environment and thus can also play a role in antibiotic resistance. The detection of structural variations in bacteria is challenging, and the recognition of even small rearrangements can be important. Even though most detection tools are aimed at and benchmarked on eukaryotic genomes, they can also be used on prokaryotic genomes. The key features of detection are the ability to detect small rearrangements and support haploid genomes. Because of the limiting performance of a single detection tool, combining the detection abilities of multiple tools can lead to more robust results. There are already available workflows for structural variation detection for long-reads technologies and for the detection of single-nucleotide variation and indels, both aimed at bacteria. Yet we are unaware of structural variations detection workflows for the short-reads sequencing platform. Motivated by this gap we created our workflow. Further, we were interested in increasing the detection performance and providing more robust results.ResultsWe developed an open-source bioinformatics pipeline, ProcaryaSV, for the detection of structural variations in bacterial isolates from paired-end short sequencing reads. Multiple tools, starting with quality control and trimming of sequencing data, alignment to the reference genome, and multiple structural variation detection tools, are integrated. All the partial results are then processed and merged with an in-house merging algorithm. Compared with a single detection approach, ProcaryaSV has improved detection performance and is a reproducible easy-to-use tool.ConclusionsThe ProcaryaSV pipeline provides an integrative approach to structural variation detection from paired-end next-generation sequencing of bacterial samples. It can be easily installed and used on Linux machines. It is publicly available on GitHub at https://github.com/robinjugas/ProcaryaSV.

Illumina short read-based shotgun metagenomic data sets of simulated bacterial communities derived from fresh spinach and surface water.

Paired-end short reads of Illumina HiSeq, MiSeq, and NovaSeq of simulated bacterial communities from fresh spinach and surface water were generated in silico at various sequencing depths. Multidrug-resistant Salmonella enterica serotype Indiana was included in the spinach community, while the water community contained multidrug-resistant Pseudomonas aeruginosa.

CAREx: context-aware read extension of paired-end sequencing data

BackgroundCommonly used next generation sequencing machines typically produce large amounts of short reads of a few hundred base-pairs in length. However, many downstream applications would generally benefit from longer reads.ResultsWe present CAREx—an algorithm for the generation of pseudo-long reads from paired-end short-read Illumina data based on the concept of repeatedly computing multiple-sequence-alignments to extend a read until its partner is found. Our performance evaluation on both simulated data and real data shows that CAREx is able to connect significantly more read pairs (up to 99%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$99\\%$$\\end{document} for simulated data) and to produce more error-free pseudo-long reads than previous approaches. When used prior to assembly it can achieve superior de novo assembly results. Furthermore, the GPU-accelerated version of CAREx exhibits the fastest execution times among all tested tools.ConclusionCAREx is a new MSA-based algorithm and software for producing pseudo-long reads from paired-end short read data. It outperforms other state-of-the-art programs in terms of (i) percentage of connected read pairs, (ii) reduction of error rates of filled gaps, (iii) runtime, and (iv) downstream analysis using de novo assembly. CAREx is open-source software written in C++ (CPU version) and in CUDA/C++ (GPU version). It is licensed under GPLv3 and can be downloaded at (https://github.com/fkallen/CAREx).

Genetic variation of the Plasmodium falciparum circumsporozoite protein in parasite isolates from Homabay County in Kenya

The Plasmodium falciparum Circumsporozoite Protein (PfCSP) has been used in developing the RTS,S, and R21 malaria vaccines. However, genetic polymorphisms within Pfcsp compromise the effectiveness of the vaccine. Thus, it is essential to continuously assess the genetic diversity of Pfcsp, especially when deploying it across different geographical regions. In this study, we assessed the genetic diversity of the Pfcsp on isolates from Homabay County, a malaria-endemic region in western Kenya, and compared it against other isolates from Kenya. We extracted DNA from 27 microscopically confirmed P. falciparum positive samples and conducted Illumina sequencing to generate paired-end short reads. The sequences were then mapped to the Pf3D7 reference genome, and genetic variation was analyzed using bcftools. Additionally, we retrieved isolates from two other malaria-endemic regions in Kenya, Kisumu (n=58) and Kilifi (n=596), from MalariaGEN version 7 and compared their genetic diversity and natural selection. We also evaluated the predicted binding affinities for HLA class I and II supertype alleles for the identified haplotypes using NetMHCpan and NetMHCIIpan. Our results show that the N-terminal of PfCSP was relatively conserved with a notable mutation at A98G across all isolates. The number of NANP repeats varied across the three Kenyan sites within the central repeat region. Furthermore, the C-terminal region showed polymorphism within the Th2R and Th3R regions. Haplotype network analysis of the Kenyan isolates revealed 69 haplotypes, with the 3D7 reference being found in the most prevalent haplotype. When assessing the predicted binding affinities between supertypes in HLA class I and II with the identified haplotypes, we observed stronger predicted binding affinities to multiple haplotypes except for those containing the 3D7 reference. The results suggest the need to take into account the existing changes occurring in Pfcsp while developing malaria vaccines.

Ocrelizumab B cell depletion has no effect on HERV RNA expression in PBMC in MS patients

BackgroundEpstein barr virus (EBV) infection of B cells is now understood to be one of the triggering events for the development of Multiple Sclerosis (MS), a progressive immune-mediated disease of the central nervous system. EBV infection is also linked to expression of human endogenous retroviruses (HERVs) of the HERV-W group, a further risk factor for the development of MS. Ocrelizumab is a high-potency disease-modifying treatment (DMT) for MS, which depletes B cells by targeting CD20. ObjectivesWe studied the effects of ocrelizumab on gene expression in peripheral blood mononuclear cells (PBMC) from paired samples from 20 patients taken prior to and 6 months after beginning ocrelizumab therapy. We hypothesised that EBV and HERV-W loads would be lower in post-treatment samples. MethodsSamples were collected in Paxgene tubes, subject to RNA extraction and Illumina paired end short read mRNA sequencing with mapping of sequence reads to the human genome using Salmon and differential gene expression compared with DeSeq2. Mapping was also performed separately to the HERV-D database of HERV sequences and the EBV reference sequence. ResultsPatient samples were more strongly clustered by individual rather than disease type (relapsing/remitting or primary progressive), treatment (pre and post), age, or sex. Fourteen genes, all clearly linked to B cell function were significantly down regulated in the post treatment samples. Interestingly only one pre-treatment sample had detectable EBV RNA and there were no significant differences in HERV expression (of any group) between pre- and post-treatment samples. ConclusionsWhile EBV and HERV expression are clearly linked to triggering MS pathogenesis, it does not appear that high level expression of these viruses is a part of the ongoing disease process or that changes in virus load are associated with ocrelizumab treatment.

Improved metagenome assemblies through selective enrichment of bacterial genomic DNA from eukaryotic host genomic DNA using ATAC-seq.

Genomics can be used to study the complex relationships between hosts and their microbiota. Many bacteria cannot be cultured in the laboratory, making it difficult to obtain adequate amounts of bacterial DNA and to limit host DNA contamination for the construction of metagenome-assembled genomes (MAGs). For example, Wolbachia is a genus of exclusively obligate intracellular bacteria that live in a wide range of arthropods and some nematodes. While Wolbachia endosymbionts are frequently described as facultative reproductive parasites in arthropods, the bacteria are obligate mutualistic endosymbionts of filarial worms. Here, we achieve 50-fold enrichment of bacterial sequences using ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) with Brugia malayi nematodes, containing Wolbachia (wBm). ATAC-seq uses the Tn5 transposase to cut and attach Illumina sequencing adapters to accessible DNA lacking histones, typically thought to be open chromatin. Bacterial and mitochondrial DNA in the lysates are also cut preferentially since they lack histones, leading to the enrichment of these sequences. The benefits of this include minimal tissue input (<1 mg of tissue), a quick protocol (<4 h), low sequencing costs, less bias, correct assembly of lateral gene transfers and no prior sequence knowledge required. We assembled the wBm genome with as few as 1 million Illumina short paired-end reads with >97% coverage of the published genome, compared to only 12% coverage with the standard gDNA libraries. We found significant bacterial sequence enrichment that facilitated genome assembly in previously published ATAC-seq data sets from human cells infected with Mycobacterium tuberculosis and C. elegans contaminated with their food source, the OP50 strain of E. coli. These results demonstrate the feasibility and benefits of using ATAC-seq to easily obtain bacterial genomes to aid in symbiosis, infectious disease, and microbiome research.

A Telomere-to-Telomere Genome Assembly Resource of Bipolaris sorokiniana, the Fungal Pathogen Causing Spot Blotch and Common Root Rot Diseases in Barley and Wheat

Bipolaris sorokiniana ( Cochliobolus sativus) is a fungal pathogen that causes spot blotch, common root rot, and kernel blight in barley and wheat. Four pathotypes (0, 1, 2, and 7) of the fungus were previously identified based on their virulence on three barley differential lines (Bowman, ND 5883, and NDB 112). Although several genome assemblies have been reported for B. sorokiniana, a telomere-to-telomere genome assembly is still lacking for this fungus. In this study, we assembled the genome of a pathotype 0 isolate (ND93-1) of B. sorokiniana using both PacBio HiFi reads and ultralong Oxford Nanopore Technology reads corrected with Illumina paired-end 100-bp short reads. The combined genome assembly of ND93-1 has an estimated size of 35.7 Mb and consists of 16 scaffolds, each having two telomeres. A total of 11,564 protein-coding genes were predicted, including 1,633 genes encoding for secretory proteins and 473 genes for effectors. This telomere-to-telomere genome assembly provides an important resource for comparative genomics and understanding molecular biology of B. sorokiniana and related Bipolaris species. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

De novo transcriptome analysis of Dysoxylum binectariferum to unravel the biosynthesis of pharmaceutically relevant specialized metabolites.

The tropical tree, D. binectariferum, is a prominent source of chromone alkaloid rohitukine, which is used in the semi-syntheses of anticancer molecules such as flavopiridol and P-276-00. The biosynthetic pathway of rohitukine or its derivatives is currently unknown in plants. Here, we explored chromone alkaloid biosynthesis in D. binectariferum through targeted transcriptome sequencing. Illumina sequencing of leaves and roots of a year-old D. binectariferum seedling generated, 42.43 and 38.74 million paired-end short reads, respectively. Quality filtering and de novo assembly of the transcriptome generated 274,970 contigs and 126,788 unigenes with an N50 contig length of 1560 bp. The assembly generated 117,619 translated unigene protein sequences and 51,598 non-redundant sequences. Nearly 80% of these non-redundant sequences were annotated to publicly available protein and nucleotide databases, suggesting the completeness and effectiveness of the transcriptome assembly. Using the assembly, we identified a chalcone synthase (CHS) and three type III polyketide synthases (PKS-III; non-CHS type) that are likely to be involved in the biosynthesis of chromone ring/noreugenin moiety of rohitukine. We also identified key enzymes like lysine decarboxylase in the piperidine pathway that make the piperidine moiety of rohitukine. Besides these, the upstream enzymes in flavonoid biosynthesis like phenylalanine ammonia-lyase (PAL), trans-cinnamate 4-hydroxylase (C4H),4-coumarate-CoA ligase (4CL), and chalcone isomerase (CHI) have also been identified. Also, terpene synthases that are likely to be involved in the biosynthesis of various terpenoid scaffolds have been identified. Together, the D. binectariferum transcriptome resource forms a basis for further exploration of biosynthetic pathways of these valuable compounds through functional validation of the candidate genes and metabolic engineering in heterologous hosts. Additionally, the transcriptome dataset generated will serve as an important resource for research on functional genomics and enzyme discovery in D. binectariferum and comparative analysis with other Meliaceae family members.

Genome sequencing and application of Taiwanese macaque Macaca cyclopis

Formosan macaque (Macaca cyclopis) is the only non-human primate in Taiwan Island. We performed de novo hybrid assembly for M. cyclopis using Illumina paired-end short reads, mate-pair reads and Nanopore long reads and obtained 5065 contigs with a N50 of 2.66 megabases. M. cyclopis contigs > = 10 kb were assigned to chromosomes using Indian rhesus macaque (Macaca mulatta mulatta) genome assembly Mmul_10 as reference, resulting in a draft of M. cyclopis genome of 2,846,042,475 bases, distributed in 21 chromosomes. The draft genome contains 23,462 transcriptional origins (genes), capable of expressing 716,231 exons in 59,484 transcripts. Genome-based phylogenetic study using the assembled M. cyclopis genome together with genomes of four other macaque species, human, orangutan and chimpanzee showed similar result as previously reported. However, the M. cyclopis species was found to diverge from Chinese M. mulatta lasiota about 1.8 million years ago. Fossil gene analysis detected the presence of gap and pol endogenous viral elements of simian retrovirus in all macaques tested, including M. fascicularis, M. m. mulatta and M. cyclopis. However, M. cyclopis showed ~ 2 times less in number and more uniform in chromosomal locations. The constrain in foreign genome disturbance, presumably due to geographical isolation, should be able to simplify genomics-related investigations, making M. cyclopis an ideal primate species for medical research.

High-Quality Genome of the Tree Pathogen Phytophthora plurivora—A Novel Resource for Epidemiological Research

Phytophthora plurivora can affect a range of ecologically and silviculturally important tree species, including European beech ( Fagus sylvatica), a common late successional tree species native to Europe. Here, we report on the high-quality genome of P. plurivora strain TJ71 (CBS 124093). We sequenced it using Oxford Nanopore MinION and PacBio Sequel II long-read sequencing with 80× coverage, chromatin conformation capture (Hi-C) sequencing with 400× coverage, and DNBSEQ 150-bp paired-end short reads sequencing with 200× coverage. This complex sequencing approach allowed for assembly of the genome at the chromosome level. Specifically, the P. plurivora genome resulted in 18 scaffolds of 47 Mbp total size with 95% completeness of the eukaryotic gene set as implemented in BUSCO. This is a considerable improvement relative to the previous NCBI reference genome of P. plurivora (NMPK00000000.1) with ∼41 Mbp organized in 1,898 scaffolds with 93.8% eukaryotic BUSCO completeness. This high-quality genome provides a valuable resource for further evolutionary, epidemiological, and population genomic studies. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

The Crown Pearl V2: an improved genome assembly of the European freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758).

Contiguous assemblies are fundamental to deciphering the composition of extant genomes. In molluscs, this is considerably challenging owing to the large size of their genomes, heterozygosity, and widespread repetitive content. Consequently, long-read sequencing technologies are fundamental for high contiguity and quality. The first genome assembly of Margaritifera margaritifera (Linnaeus, 1758) (Mollusca: Bivalvia: Unionida), a culturally relevant, widespread, and highly threatened species of freshwater mussels, was recently generated. However, the resulting genome is highly fragmented since the assembly relied on short-read approaches. Here, an improved reference genome assembly was generated using a combination of PacBio CLR long reads and Illumina paired-end short reads. This genome assembly is 2.4 Gb long, organized into 1,700 scaffolds with a contig N50 length of 3.4 Mbp. The ab initio gene prediction resulted in 48,314 protein-coding genes. Our new assembly is a substantial improvement and an essential resource for studying this species' unique biological and evolutionary features, helping promote its conservation.

Investigating Possible Interspecies Communication of Plasmids Associated with Transfer of Third-Generation Cephalosporin, Quinolone, and Colistin Resistance Between Simultaneously Isolated Escherichia Coli and Klebsiella Pneumoniae.

The coinfection process producing multiple species of pathogens provides a specific ecological niche for the exchange of genetic materials between pathogens, in which plasmids play a vital role in horizontal gene transfer, especially for drug resistance, but the underlying transfer pathway remains unclear. Interspecies communication of the plasmids associated with the transfer of third-generation cephalosporins, quinolones, and colistin resistance has been observed in simultaneously isolated Escherichia coli and Klebsiella pneumoniae from abdominal drainage following surgery. The MICs of antimicrobial agents were determined by the broth microdilution method. The complete chromosome and plasmid sequences were obtained by combining Illumina paired-end short reads and MinION long reads. S1-PFGE, southern blot analysis and conjugation assay confirmed the transferability of the mcr-1-harboring plasmid. Both the E. coli isolate EC15255 and K. pneumoniae isolate KP15255 from the same specimen presented multidrug resistance. Each of them harbored one chromosome and three plasmids, and two plasmids and their mediated resistance could be transferred to the recipient by conjugation. Comparison of their genome sequences suggested that several genetic communication events occurred between species, especially among their plasmids, such as whole-plasmid transfer, insertion, deletion, amplification, or inversion. Exchange of plasmids or the genetic elements they harbor plays a critical role in antimicrobial resistance gene transmission and poses a substantial threat to nosocomial infection control, necessitating the continued surveillance of multidrug resistant pathogens, especially during coinfection. IMPORTANCE The genome sequence of bacterial pathogens commonly provides a detailed clue of genetic communication among clones or even distinct species. The intestinal microecological environment is a representative ecological niche for genetic communication. However, it is still difficult to describe the details of horizontal gene transfer or other genetic events within them because the evidence in the genome sequence is incomplete and limited. In this study, the simultaneously isolated Escherichia coli and Klebsiella pneumoniae from a coinfection process provided an excellent example for observation of interspecies communication between the two genomes and the plasmids they harbor. A complete genome sequence acquired by combining the Illumina and MinION sequencing platforms facilitated the understanding of genetic communication events, such as whole-plasmid transfer, insertion, deletion, amplification, or inversion, which contribute to antimicrobial resistance gene transmission and are a substantial threat to nosocomial infection control.

Critical Assessment of Short-Read Assemblers for the Metagenomic Identification of Foodborne and Waterborne Pathogens Using Simulated Bacterial Communities.

Metagenomics offers the highest level of strain discrimination of bacterial pathogens from complex food and water microbiota. With the rapid evolvement of assembly algorithms, defining an optimal assembler based on the performance in the metagenomic identification of foodborne and waterborne pathogens is warranted. We aimed to benchmark short-read assemblers for the metagenomic identification of foodborne and waterborne pathogens using simulated bacterial communities. Bacterial communities on fresh spinach and in surface water were simulated by generating paired-end short reads of Illumina HiSeq, MiSeq, and NovaSeq at different sequencing depths. Multidrug-resistant Salmonella Indiana SI43 and Pseudomonas aeruginosa PAO1 were included in the simulated communities on fresh spinach and in surface water, respectively. ABySS, IDBA-UD, MaSuRCA, MEGAHIT, metaSPAdes, and Ray Meta were benchmarked in terms of assembly quality, identifications of plasmids, virulence genes, Salmonella pathogenicity island, antimicrobial resistance genes, chromosomal point mutations, serotyping, multilocus sequence typing, and whole-genome phylogeny. Overall, MEGHIT, metaSPAdes, and Ray Meta were more effective for metagenomic identification. We did not obtain an optimal assembler when using the extracted reads classified as Salmonella or P. aeruginosa for downstream genomic analyses, but the extracted reads showed consistent phylogenetic topology with the reference genome when they were aligned with Salmonella or P. aeruginosa strains. In most cases, HiSeq, MiSeq, and NovaSeq were comparable at the same sequencing depth, while higher sequencing depths generally led to more accurate results. As assembly algorithms advance and mature, the evaluation of assemblers should be a continuous process.

A long-read and short-read transcriptomics approach provides the first high-quality reference transcriptome and genome annotation for Pseudotsuga menziesii (Douglas-fir).

Douglas-fir (Pseudotsuga menziesii) is native to western North America. It grows in a wide range of environmental conditions and is an important timber tree. Although there are several studies on the gene expression responses of Douglas-fir to abiotic cues, the absence of high-quality transcriptome and genome data is a barrier to further investigation. Like for most conifers, the available transcriptome and genome reference dataset for Douglas-fir remains fragmented and requires refinement. We aimed to generate a highly accurate, and complete reference transcriptome and genome annotation. We deep-sequenced the transcriptome of Douglas-fir needles from seedlings that were grown under nonstress control conditions or a combination of heat and drought stress conditions using long-read (LR) and short-read (SR) sequencing platforms. We used 2 computational approaches, namely de novo and genome-guided LR transcriptome assembly. Using the LR de novo assembly, we identified 1.3X more high-quality transcripts, 1.85X more "complete" genes, and 2.7X more functionally annotated genes compared to the genome-guided assembly approach. We predicted 666 long noncoding RNAs and 12,778 unique protein-coding transcripts including 2,016 putative transcription factors. We leveraged the LR de novo assembled transcriptome with paired-end SR and a published single-end SR transcriptome to generate an improved genome annotation. This was conducted with BRAKER2 and refined based on functional annotation, repetitive content, and transcriptome alignment. This high-quality genome annotation has 51,419 unique gene models derived from 322,631 initial predictions. Overall, our informatics approach provides a new reference Douglas-fir transcriptome assembly and genome annotation with considerably improved completeness and functional annotation.

MarkerMAG: linking metagenome-assembled genomes (MAGs) with 16S rRNA marker genes using paired-end short reads.

Metagenome-assembled genomes (MAGs) have substantially extended our understanding of microbial functionality. However, 16S rRNA genes, which are commonly used in phylogenetic analysis and environmental surveys, are often missing from MAGs. Here, we developed MarkerMAG, a pipeline that links 16S rRNA genes to MAGs using paired-end sequencing reads. Assessment of MarkerMAG on three benchmarking metagenomic datasets with various degrees of complexity shows substantial increases in the number of MAGs with 16S rRNA genes and a 100% assignment accuracy. MarkerMAG also estimates the copy number of 16S rRNA genes in MAGs with high accuracy. Assessments on three real metagenomic datasets demonstrate 1.1- to 14.2-fold increases in the number of MAGs with 16S rRNA genes. We also show that MarkerMAG-improved MAGs increase the accuracy of functional prediction from 16S rRNA gene amplicon data. MarkerMAG is helpful in connecting information in MAG databases with those in 16S rRNA databases and surveys and hence contributes to our increasing understanding of microbial diversity, function and phylogeny. MarkerMAG is implemented in Python3 and freely available at https://github.com/songweizhi/MarkerMAG. Supplementary data are available at Bioinformatics online.

Genome Sequencing of Amomum tsao-ko Provides Novel Insight Into Its Volatile Component Biosynthesis.

As an important economic and medicinal crop, Amomum tsao-ko is rich in volatile oils and widely used in food additives, essential oils, and traditional Chinese medicine. However, the lack of the genome remains a limiting factor for understanding its medicinal properties at the molecular level. Here, based on 288.72 Gb of PacBio long reads and 105.45 Gb of Illumina paired-end short reads, we assembled a draft genome for A. tsao-ko (2.70 Gb in size, contig N50 of 2.45 Mb). Approximately 90.07% of the predicted genes were annotated in public databases. Based on comparative genomic analysis, genes involved in secondary metabolite biosynthesis, flavonoid metabolism, and terpenoid biosynthesis showed significant expansion. Notably, the DXS, GGPPS, and CYP450 genes, which participate in rate-limiting steps for terpenoid backbone biosynthesis and modification, may form the genetic basis for essential oil formation in A. tsao-ko. The assembled A. tsao-ko draft genome provides a valuable genetic resource for understanding the unique features of this plant and for further evolutionary and agronomic studies of Zingiberaceae species.

Identification of Protein Isoforms Using Reference Databases Built from Long and Short Read RNA-Sequencing.

Alternative splicing can lead to distinct protein isoforms. These can have different functions in specific cells and tissues or in different developmental stages. In this study, we explored whether transcripts assembled from long read, nanopore-based, direct RNA-sequencing (RNA-seq) could improve the identification of protein isoforms in human K562 cells. By comparing with Illumina-based short read RNA-seq, we showed that a large proportion of Ensembl transcripts (5949/14,326) and genes expressing alternatively spliced transcripts (486/2981) identified with long direct reads were missed by short paired-end reads. By co-analyzing proteomic and transcriptomic data, we also showed that some peptides (826/35,976), proteins (262/3215), and protein isoforms arising from distinct transcript variants (574/1212) identified with isoform-specific peptides via custom long-read-based databases were missed in Illumina-derived databases. Finally, we generated unequivocal peptide evidence for a set of protein isoforms and showed that long read, direct RNA-seq allows the discovery of novel protein isoforms not already in reference databases or custom databases built from short read RNA-seq data. Our analysis highlights the benefits of long read RNA-seq data in the generation of reference databases to increase tandem mass spectrometry (MS/MS) identification of protein isoforms.