Using genotype array data to compare multi- and single-sample variant calls and improve variant call sets from deep coverage whole-genome sequencing data

Time-Series Expression Profile Analysis of Post-Traumatic Joint Contracture in Rats at the Early Stages of the Healing Process

Objective

This study aimed to characterize the gene expression profile at the early stages of the healing process of post-traumatic joint contracture (PTJC).

Methods

Twelve rats were used for PTJC model establishment and were divided into four groups according to the sampling time: S0d, S3d, S7d and S2w. Transcriptome sequencing was performed on fibrotic joint capsule samples in four groups followed by bioinformatics analyses including differentially expressed genes (DEGs) screening, Short Time-series Expression Miner (STEM) analysis, network construction, and pathway analysis. Five important genes were validated by qRT-PCR.

Results

A total of 1171, 1052 and 793 DEGs were screened in S3d vs S0d, S7d vs S0d, and S2w vs S0d comparison groups, respectively. A total of 383 overlapping genes were screened out, which were significantly enriched in some inflammatory functions and pathways. Through STEM analysis, three clusters were identified, including 105, 57 and 57 DEGs, respectively. Then, based on the cluster genes, 10 genes, such as Il6, Timp1, Cxcl1, Cxcr4 and Mmp3, were further selected after PPI and pathway analyses. The expression levels of Il6, Timp1, Cxcl1, Cxcr4 and Mmp3 were validated by qRT-PCR.

Conclusion

The present study screened out several genes with significant changes in expression levels at the early stages of the healing process in PTJC, such as Il6, Timp1, Cxcl1, Cxcr4 and Mmp3. Our study offers a valuable contribution to the understanding pathomechanism of PTJC.

Evaluation of serverless computing for scalable execution of a joint variant calling workflow

Advances in whole-genome sequencing have greatly reduced the cost and time of obtaining raw genetic information, but the computational requirements of analysis remain a challenge. Serverless computing has emerged as an alternative to using dedicated compute resources, but its utility has not been widely evaluated for standardized genomic workflows. In this study, we define and execute a best-practice joint variant calling workflow using the SWEEP workflow management system. We present an analysis of performance and scalability, and discuss the utility of the serverless paradigm for executing workflows in the field of genomics research. The GATK best-practice short germline joint variant calling pipeline was implemented as a SWEEP workflow comprising 18 tasks. The workflow was executed on Illumina paired-end read samples from the European and African super populations of the 1000 Genomes project phase III. Cost and runtime increased linearly with increasing sample size, although runtime was driven primarily by a single task for larger problem sizes. Execution took a minimum of around 3 hours for 2 samples, up to nearly 13 hours for 62 samples, with costs ranging from $2 to $70.

The application of deep learning for the classification of correct and incorrect SNP genotypes from whole-genome DNA sequencing pipelines

A downside of next-generation sequencing technology is the high technical error rate. We built a tool, which uses array-based genotype information to classify next-generation sequencing–based SNPs into the correct and the incorrect calls. The deep learning algorithms were implemented via Keras. Several algorithms were tested: (i) the basic, naïve algorithm, (ii) the naïve algorithm modified by pre-imposing different weights on incorrect and correct SNP class in calculating the loss metric and (iii)–(v) the naïve algorithm modified by random re-sampling (with replacement) of the incorrect SNPs to match 30%/60%/100% of the number of correct SNPs. The training data set was composed of data from three bulls and consisted of 2,227,995 correct (97.94%) and 46,920 incorrect SNPs, while the validation data set consisted of data from one bull with 749,506 correct (98.05%) and 14,908 incorrect SNPs. The results showed that for a rare event classification problem, like incorrect SNP detection in NGS data, the most parsimonious naïve model and a model with the weighting of SNP classes provided the best results for the classification of the validation data set. Both classified 19% of truly incorrect SNPs as incorrect and 99% of truly correct SNPs as correct and resulted in the F1 score of 0.21 — the highest among the compared algorithms. We conclude the basic models were less adapted to the specificity of a training data set and thus resulted in better classification of the independent, validation data set, than the other tested models.