Efficient Bioinformatics Pipeline Research Articles

A Customized Human Mitochondrial DNA Database (hMITO DB v1.0) for Rapid Sequence Analysis, Haplotyping and Geo-Mapping.

The field of mitochondrial genomics has advanced rapidly and has revolutionized disciplines such as molecular anthropology, population genetics, and medical genetics/oncogenetics. However, mtDNA next-generation sequencing (NGS) analysis for matrilineal haplotyping and phylogeographic inference remains hindered by the lack of a consolidated mitogenome database and an efficient bioinformatics pipeline. To address this, we developed a customized human mitogenome database (hMITO DB) embedded in a CLC Genomics workflow for read mapping, variant analysis, haplotyping, and geo-mapping. The database was constructed from 4286 mitogenomes. The macro-haplogroup (A to Z) distribution and representative phylogenetic tree were found to be consistent with published literature. The hMITO DB automated workflow was tested using mtDNA-NGS sequences derived from Pap smears and cervical cancer cell lines. The auto-generated read mapping, variants track, and table of haplotypes and geo-origins were completed in 15 min for 47 samples. The mtDNA workflow proved to be a rapid, efficient, and accurate means of sequence analysis for translational mitogenomics.

A Customized Monkeypox Virus Genomic Database (MPXV DB v1.0) for Rapid Sequence Analysis and Phylogenomic Discoveries in CLC Microbial Genomics.

Monkeypox has been a neglected, zoonotic tropical disease for over 50 years. Since the 2022 global outbreak, hundreds of human clinical samples have been subjected to next-generation sequencing (NGS) worldwide with raw data deposited in public repositories. However, sequence analysis for in-depth investigation of viral evolution remains hindered by the lack of a curated, whole genome Monkeypox virus (MPXV) database (DB) and efficient bioinformatics pipelines. To address this, we developed a customized MPXV DB for integration with "ready-to-use" workflows in the CLC Microbial Genomics Module for whole genomic and metagenomic analysis. After database construction (218 MPXV genomes), whole genome alignment, pairwise comparison, and evolutionary analysis of all genomes were analyzed to autogenerate tabular outputs and visual displays (collective runtime: 16 min). The clinical utility of the MPXV DB was demonstrated by using a Chimpanzee fecal, hybrid-capture NGS dataset (publicly available) for metagenomic, phylogenomic, and viral/host integration analysis. The clinically relevant MPXV DB embedded in CLC workflows proved to be a rapid method of sequence analysis useful for phylogenomic exploration and a wide range of applications in translational science.

Deciphering complex genome rearrangements in C. elegans using short-read whole genome sequencing

Genomic rearrangements cause congenital disorders, cancer, and complex diseases in human. Yet, they are still understudied in rare diseases because their detection is challenging, despite the advent of whole genome sequencing (WGS) technologies. Short-read (srWGS) and long-read WGS approaches are regularly compared, and the latter is commonly recommended in studies focusing on genomic rearrangements. However, srWGS is currently the most economical, accurate, and widely supported technology. In Caenorhabditis elegans (C. elegans), such variants, induced by various mutagenesis processes, have been used for decades to balance large genomic regions by preventing chromosomal crossover events and allowing the maintenance of lethal mutations. Interestingly, those chromosomal rearrangements have rarely been characterized on a molecular level. To evaluate the ability of srWGS to detect various types of complex genomic rearrangements, we sequenced three balancer strains using short-read Illumina technology. As we experimentally validated the breakpoints uncovered by srWGS, we showed that, by combining several types of analyses, srWGS enables the detection of a reciprocal translocation (eT1), a free duplication (sDp3), a large deletion (sC4), and chromoanagenesis events. Thus, applying srWGS to decipher real complex genomic rearrangements in model organisms may help designing efficient bioinformatics pipelines with systematic detection of complex rearrangements in human genomes.

Fast-GBS v2.0: an analysis toolkit for genotyping-by-sequencing data.

Genotyping-by-sequencing (GBS) is a rapid, flexible, low-cost, and robust genotyping method that simultaneously discovers variants and calls genotypes within a broad range of samples. These characteristics make GBS an excellent tool for many applications and research questions from conservation biology to functional genomics in both model and non-model species. Continued improvement of GBS relies on a more comprehensive understanding of data analysis, development of fast and efficient bioinformatics pipelines, accurate missing data imputation, and active post-release support. Here, we present the second generation of Fast-GBS (v2.0) that offers several new options (e.g., processing paired-end reads and imputation of missing data) and features (e.g., summary statistics of genotypes) to improve the GBS data analysis process. The performance assessment analysis showed that Fast-GBS v2.0 outperformed other available analytical pipelines, such as GBS-SNP-CROP and Gb-eaSy. Fast-GBS v2.0 provides an analysis platform that can be run with different types of sequencing data, modest computational resources, and allows for missing-data imputation for various species in different contexts.

CSA: A high-throughput chromosome-scale assembly pipeline for vertebrate genomes.

BackgroundEasy-to-use and fast bioinformatics pipelines for long-read assembly that go beyond the contig level to generate highly continuous chromosome-scale genomes from raw data remain scarce.ResultChromosome-Scale Assembler (CSA) is a novel computationally highly efficient bioinformatics pipeline that fills this gap. CSA integrates information from scaffolded assemblies (e.g., Hi-C or 10X Genomics) or even from diverged reference genomes into the assembly process. As CSA performs automated assembly of chromosome-sized scaffolds, we benchmark its performance against state-of-the-art reference genomes, i.e., conventionally built in a laborious fashion using multiple separate assembly tools and manual curation. CSA increases the contig lengths using scaffolding, local re-assembly, and gap closing. On certain datasets, initial contig N50 may be increased up to 4.5-fold. For smaller vertebrate genomes, chromosome-scale assemblies can be achieved within 12 h using low-cost, high-end desktop computers. Mammalian genomes can be processed within 16 h on compute-servers. Using diverged reference genomes for fish, birds, and mammals, we demonstrate that CSA calculates chromosome-scale assemblies from long-read data and genome comparisons alone. Even contig-level draft assemblies of diverged genomes are helpful for reconstructing chromosome-scale sequences. CSA is also capable of assembling ultra-long reads.ConclusionsCSA can speed up and simplify chromosome-level assembly and significantly lower costs of large-scale family-level vertebrate genome projects.

Improving Copy Number Variant Detection from Sequencing Data with a Combination of Programs and a Predictive Model

Bioinformatics tools for analyzing copy number variants (CNVs) from massively parallel sequencing (MPS) data are less well developed compared with other variant types. We present an efficient bioinformatics pipeline for CNV detection from gene panel MPS data in neuromuscular disorders. CNVs were generated in silico into samples sequenced with a previously published MPS gene panel. The in silico CNVs from these samples were analyzed with four programs having complementary CNV detection ranges: CoNIFER, XHMM, ExomeDepth, and CODEX. A logistic regression model was trained with the obtained set of in silico CNV detections to predict true-positive CNV detections among all CNV detections from samples. This model was validated using 66 control samples with a verified true-positive (n=58) or false-positive (n=8) CNV detection. Applying all four programs together provided more sensitive detection results with in silico CNVs than other program combinations or any program alone. Furthermore, a model with CNV detection-specific scores from all four programs as variables performed overall best in the validation. No single program could detect all CNV sizes and types equally or with enough accuracy. Therefore, a combination of carefully selected programs should be used to maximize detection accuracy. In addition, the detected CNVs should be reviewed with a statistical model to streamline and standardize the filtering of the detections for annotation.

DNAscan: personal computer compatible NGS analysis, annotation and visualisation

BackgroundNext Generation Sequencing (NGS) is a commonly used technology for studying the genetic basis of biological processes and it underpins the aspirations of precision medicine. However, there are significant challenges when dealing with NGS data. Firstly, a huge number of bioinformatics tools for a wide range of uses exist, therefore it is challenging to design an analysis pipeline. Secondly, NGS analysis is computationally intensive, requiring expensive infrastructure, and many medical and research centres do not have adequate high performance computing facilities and cloud computing is not always an option due to privacy and ownership issues. Finally, the interpretation of the results is not trivial and most available pipelines lack the utilities to favour this crucial step.ResultsWe have therefore developed a fast and efficient bioinformatics pipeline that allows for the analysis of DNA sequencing data, while requiring little computational effort and memory usage. DNAscan can analyse a whole exome sequencing sample in 1 h and a 40x whole genome sequencing sample in 13 h, on a midrange computer. The pipeline can look for single nucleotide variants, small indels, structural variants, repeat expansions and viral genetic material (or any other organism). Its results are annotated using a customisable variety of databases and are available for an on-the-fly visualisation with a local deployment of the gene.iobio platform. DNAscan is implemented in Python. Its code and documentation are available on GitHub: https://github.com/KHP-Informatics/DNAscan. Instructions for an easy and fast deployment with Docker and Singularity are also provided on GitHub.ConclusionsDNAscan is an extremely fast and computationally efficient pipeline for analysis, visualization and interpretation of NGS data. It is designed to provide a powerful and easy-to-use tool for applications in biomedical research and diagnostic medicine, at minimal computational cost. Its comprehensive approach will maximise the potential audience of users, bringing such analyses within the reach of non-specialist laboratories, and those from centres with limited funding available.

A new versatile primer set targeting a short fragment of the mitochondrial COI region for metabarcoding metazoan diversity: application for characterizing coral reef fish gut contents

IntroductionThe PCR-based analysis of homologous genes has become one of the most powerful approaches for species detection and identification, particularly with the recent availability of Next Generation Sequencing platforms (NGS) making it possible to identify species composition from a broad range of environmental samples. Identifying species from these samples relies on the ability to match sequences with reference barcodes for taxonomic identification. Unfortunately, most studies of environmental samples have targeted ribosomal markers, despite the fact that the mitochondrial Cytochrome c Oxidase subunit I gene (COI) is by far the most widely available sequence region in public reference libraries. This is largely because the available versatile (“universal”) COI primers target the 658 barcoding region, whose size is considered too large for many NGS applications. Moreover, traditional barcoding primers are known to be poorly conserved across some taxonomic groups.ResultsWe first design a new PCR primer within the highly variable mitochondrial COI region, the “mlCOIintF” primer. We then show that this newly designed forward primer combined with the “jgHCO2198” reverse primer to target a 313 bp fragment performs well across metazoan diversity, with higher success rates than versatile primer sets traditionally used for DNA barcoding (i.e. LCO1490/HCO2198). Finally, we demonstrate how the shorter COI fragment coupled with an efficient bioinformatics pipeline can be used to characterize species diversity from environmental samples by pyrosequencing. We examine the gut contents of three species of planktivorous and benthivorous coral reef fish (family: Apogonidae and Holocentridae). After the removal of dubious COI sequences, we obtained a total of 334 prey Operational Taxonomic Units (OTUs) belonging to 14 phyla from 16 fish guts. Of these, 52.5% matched a reference barcode (>98% sequence similarity) and an additional 32% could be assigned to a higher taxonomic level using Bayesian assignment.ConclusionsThe molecular analysis of gut contents targeting the 313 COI fragment using the newly designed mlCOIintF primer in combination with the jgHCO2198 primer offers enormous promise for metazoan metabarcoding studies. We believe that this primer set will be a valuable asset for a range of applications from large-scale biodiversity assessments to food web studies.

Library-free methylation sequencing with bisulfite padlock probes.

We previously developed Bisulfite Padlock Probes (BSPPs) for the specific and parallel digital quantification of DNA methylation1. Here we report the second-generation of BSPPs with a design algorithm to generate more efficient padlock probes, a library-free protocol that dramatically reduces the time and cost of sample preparation and is compatible with automation, and an efficient bioinformatics pipeline that accurately obtains both methylation levels and genotypes from bisulfite sequencing data.

CARMEN - Comparative Analysis and in silico Reconstruction of organism-specific MEtabolic Networks

New sequencing technologies provide ultra-fast access to novel microbial genome data. For their interpretation, an efficient bioinformatics pipeline that facilitates in silico reconstruction of metabolic networks is highly desirable. The software tool CARMEN performs in silico reconstruction of metabolic networks to interpret genome data in a functional context. CARMEN supports the visualization of automatically derived metabolic networks based on pathway information from the KEGG database or from user-defined SBML templates; this software also enables comparative genomics. The reconstructed networks are stored in standardized SBML format. We demonstrated the functionality of CARMEN with a major application example focusing on the reconstruction of glycolysis and related metabolic reactions of Xanthomonas campestris pv. campestris B100. The curation of such pathways facilitates enhanced visualization of experimental results, simulations and comparative genomics. A second application of this software was performed on a set of corynebacteria to compare and to visualize their carbohydrate metabolism. In conclusion, using CARMEN, we developed highly automated data analysis software that rapidly converts sequence data into new knowledge, replacing the time-consuming manual reconstruction of metabolic networks. This tool is particularly useful for obtaining an overview of newly sequenced genomes and their metabolic blueprints and for comparative genome analysis. The generated pathways provide automated access to modeling and simulation tools that are compliant with the SBML standard. A user-friendly web interface of CARMEN is available at http://carmen.cebitec.uni-bielefeld.de.

Evaluation of the bacterial diversity in the feces of cattle using 16S rDNA bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP)

BackgroundThe microbiota of an animal's intestinal tract plays important roles in the animal's overall health, productivity and well-being. There is still a scarcity of information on the microbial diversity in the gut of livestock species such as cattle. The primary reason for this lack of data relates to the expense of methods needed to generate such data. Here we have utilized a bacterial tag-encoded FLX 16s rDNA amplicon pyrosequencing (bTEFAP) approach that is able to perform diversity analyses of gastrointestinal populations. bTEFAP is relatively inexpensive in terms of both time and labor due to the implementation of a novel tag priming method and an efficient bioinformatics pipeline. We have evaluated the microbiome from the feces of 20 commercial, lactating dairy cows.ResultsUbiquitous bacteria detected from the cattle feces included Clostridium, Bacteroides, Porpyhyromonas, Ruminococcus, Alistipes, Lachnospiraceae, Prevotella, Lachnospira, Enterococcus, Oscillospira, Cytophage, Anaerotruncus, and Acidaminococcus spp. Foodborne pathogenic bacteria were detected in several of the cattle, a total of 4 cows were found to be positive for Salmonella spp (tentative enterica) and 6 cows were positive for Campylobacter spp. (tentative lanienae).ConclusionUsing bTEFAP we have examined the microbiota in the feces of cattle. As these methods continue to mature we will better understand the ecology of the major populations of bacteria the lower intestinal tract. This in turn will allow for a better understanding of ways in which the intestinal microbiome contributes to animal health, productivity and wellbeing.