Abstract:Abstract Background Optimization of an RNA-Sequencing (RNA-Seq) pipeline is critical to maximize power and accuracy to identify genetic variants, including SNPs, which may serve as genetic markers to select for feed efficiency, leading to economic benefits for beef production. This study used RNA-Seq data (GEO Accession ID: PRJEB7696 and PRJEB15314) from muscle and liver tissue, respectively, from 12 Nellore beef steers selected from 585 steers with residual feed intake measures (RFI; n = 6 low-RFI, n = 6 high-RFI). Three RNA-Seq pipelines were compared including multi-sample calling from i) non-merged samples; ii) merged samples by RFI group, iii) merged samples by RFI and tissue group. The RNA-Seq reads were aligned against the UMD3.1 bovine reference genome (release 94) assembly using STAR aligner. Variants were called using BCFtools and variant effect prediction (VeP) and functional annotation (ToppGene) analyses were performed. Results On average, total reads detected for Approach i) non-merged samples for liver and muscle, were 18,362,086.3 and 35,645,898.7, respectively. For Approach ii), merging samples by RFI group, total reads detected for each merged group was 162,030,705, and for Approach iii), merging samples by RFI group and tissues, was 324,061,410, revealing the highest read depth for Approach iii). Additionally, Approach iii) merging samples by RFI group and tissues, revealed the highest read depth per variant coverage (572.59 ± 3993.11) and encompassed the majority of localized positional genes detected by each approach. This suggests Approach iii) had optimized detection power, read depth, and accuracy of SNP calling, therefore increasing confidence of variant detection and reducing false positive detection. Approach iii) was then used to detect unique SNPs fixed within low- (12,145) and high-RFI (14,663) groups. Functional annotation of SNPs revealed positional candidate genes, for each RFI group (2886 for low-RFI, 3075 for high-RFI), which were significantly ( P < 0.05) associated with immune and metabolic pathways. Conclusion The most optimized RNA-Seq pipeline allowed for more accurate identification of SNPs, associated positional candidate genes, and significantly associated metabolic pathways in muscle and liver tissues, providing insight on the underlying genetic architecture of feed efficiency in beef cattle.

Nextflow vs. plain Bash: Different Approaches to the Parallelisation of SNP Calling from the Whole Genome Sequence Data

RNAflow: An Effective and Simple RNA-Seq Differential Gene Expression Pipeline Using Nextflow

Gene Sequence Alignment on a Public Computing Platform

eQTL-Detect: nextflow-based pipeline for eQTL detection in modular format with sharable and parallelizable scripts

GenMPI: Cluster Scalable Variant Calling for Short/Long Reads Sequencing Data

DNAscan: a fast, computationally and memory efficient bioinformatics pipeline for the analysis of DNA next-generation-sequencing data

Comparison of high-throughput single-cell RNA sequencing data processing pipelines

Sentieon DNASeq Variant Calling Workflow Demonstrates Strong Computational Performance and Accuracy

Metapipeline-DNA: A Comprehensive Germline & Somatic Genomics Nextflow Pipeline

Computational performance and accuracy of Sentieon DNASeq variant calling workflow

VC@Scale: Scalable and high-performance variant calling on cluster environments

Genome-Wide SNP Calling from Genotyping by Sequencing (GBS) Data: A Comparison of Seven Pipelines and Two Sequencing Technologies

Accelerating genomic workflows using NVIDIA Parabricks

Scalable and efficient DNA sequencing analysis on different compute infrastructures aiding variant discovery

Scalepopgen: bioinformatics workflow resources implemented in Nextflow for comprehensive population genomic analyses

A graphical, interactive and GPU-enabled workflow to process long-read sequencing data

Benchmarking of bioinformatics tools for the hybrid de novo assembly of human whole-genome sequencing data

Development and comparison of RNA-sequencing pipelines for more accurate SNP identification: practical example of functional SNP detection associated with feed efficiency in Nellore beef cattle

A Novel Multi-Alignment Pipeline for High-Throughput Sequencing Data.

Paean: A Parallel Transcriptome Quantification Tool Combining Gene Expression And Alternative Splicing Events Using Gpu

Modelling Computational Resources for Next Generation Sequencing Bioinformatics Analysis of 16S rRNA Samples