A benchmark of batch-effect correction methods for single-cell RNA sequencing data.

BACKGROUND: Large-scale single-cell transcriptomic datasets generated using different technologies contain batch-specific systematic variations that present a challenge to batch-effect removal and data integration. With continued growth expected in scRNA-seq data, achieving effective batch integration with available computational resources is crucial. Here, we perform an in-depth benchmark study on available batch correction methods to determine the most suitable method for batch-effect removal. RESULTS: We compare 14 methods in terms of computational runtime, the ability to handle large datasets, and batch-effect correction efficacy while preserving cell type purity. Five scenarios are designed for the study: identical cell types with different technologies, non-identical cell types, multiple batches, big data, and simulated data. Performance is evaluated using four benchmarking metrics including kBET, LISI, ASW, and ARI. We also investigate the use of batch-corrected data to study differential gene expression. CONCLUSION: Based on our results, Harmony, LIGER, and Seurat 3 are the recommended methods for batch integration. Due to its significantly shorter runtime, Harmony is recommended as the first method to try, with the other methods as viable alternatives.

Genetic variation of the Turnip mosaic virus population of Vietnam: a case study of founder, regional and local influences.

Turnip mosaic virus (TuMV) is one of the most important viruses infecting a wide range of plant species, primarily from the family Brassicaceae. Thirty TuMV isolates were collected from Brassica and Raphanus plants in Vietnam during 2006-2008. Host reaction studies showed that many of the isolates belonged to Brassica/Raphanus (BR) host-infecting type. Sequence-based phylogenetic and population genetic analyses were made of the complete polyprotein gene sequences, and of four non-recombinogenic regions of those sequences (i.e. genes of the helper-component proteinase protein, protein 3, nuclear inclusion b protein and coat protein). These were used to assess the subpopulation differentiation and divergence between Vietnamese TuMV populations and those of nearby Asian countries. Nine inter- and intralineage recombination type patterns were identified in the genomes of the Vietnamese isolates, of which seven were novel. All the Vietnamese non-recombinant isolates fell into the world-B group of TuMV and clustered with Chinese isolates. The estimates of genetic differentiation and gene flow reveal that the TuMV populations of Vietnam, China and Japan are genetically linked but have clear local founder effects. This, the first population genetic study of a TuMV population in Southeast Asia, indicates the importance of such studies for providing the scientific basis of control strategies.