SCCNAInfer: a robust and accurate tool to infer the absolute copy number on scDNA-seq data.

MOTIVATION: Copy number alterations (CNAs) play an important role in disease progression, especially in cancer. Single-cell DNA sequencing (scDNA-seq) facilitates the detection of CNAs of each cell that is sequenced at a shallow and uneven coverage. However, the state-of-the-art CNA detection tools based on scDNA-seq are still subject to genome-wide errors due to the wrong estimation of the ploidy. RESULTS: We developed SCCNAInfer, a computational tool that utilizes the subclonal signal inside the tumor cells to more accurately infer each cell's ploidy and CNAs. Given the segmentation result of an existing CNA detection method, SCCNAInfer clusters the cells, infers the ploidy of each subclone, refines the read count by bin clustering, and accurately infers the CNAs for each cell. Both simulated and real datasets show that SCCNAInfer consistently improves upon the state-of-the-art CNA detection tools such as Aneufinder, Ginkgo, SCOPE and SeCNV. AVAILABILITY AND IMPLEMENTATION: SCCNAInfer is freely available at https://github.com/compbio-mallory/SCCNAInfer. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Laminar pattern of adolescent development changes in working memory neuronal activity.

Adolescent development is characterized by an improvement in cognitive abilities, such as working memory. Neurophysiological recordings in a nonhuman primate model of adolescence have revealed changes in neural activity that mirror improvement in behavior, including higher firing rate during the delay intervals of working memory tasks. The laminar distribution of these changes is unknown. By some accounts, persistent activity is more pronounced in superficial layers, so we sought to determine whether changes are most pronounced there. We therefore analyzed neurophysiological recordings from the young and adult stage of male monkeys, at different cortical depths. Superficial layers exhibited an increased baseline firing rate in the adult stage. Unexpectedly, we also detected substantial increases in delay period activity in the middle layers after adolescence, which was confirmed even after excluding penetrations near sulci. Finally, improved discriminability around the saccade period was most evident in the deeper layers. These results reveal the laminar pattern of neural activity maturation that is associated with cognitive improvement.NEW & NOTEWORTHY Structural brain changes are evident during adolescent development particularly in the cortical thickness of the prefrontal cortex, at a time when working memory ability increases markedly. The depth distribution of neurophysiological changes during adolescence is not known. Here, we show that neurophysiological changes are not confined to superficial layers, which have most often been implicated in the maintenance of working memory. Contrary to expectations, substantial changes were evident in intermediate layers of the prefrontal cortex.

VolcanoSV enables accurate and robust structural variant calling in diploid genomes from single-molecule long read sequencing.

Structural variants (SVs) significantly contribute to human genome diversity and play a crucial role in precision medicine. Although advancements in single-molecule long-read sequencing offer a groundbreaking resource for SV detection, identifying SV breakpoints and sequences accurately and robustly remains challenging. We introduce VolcanoSV, an innovative hybrid SV detection pipeline that utilizes both a reference genome and local de novo assembly to generate a phased diploid assembly. VolcanoSV uses phased SNPs and unique k-mer similarity analysis, enabling precise haplotype-resolved SV discovery. VolcanoSV is adept at constructing comprehensive genetic maps encompassing SNPs, small indels, and all types of SVs, making it well-suited for human genomics studies. Our extensive experiments demonstrate that VolcanoSV surpasses state-of-the-art assembly-based tools in the detection of insertion and deletion SVs, exhibiting superior recall, precision, F1 scores, and genotype accuracy across a diverse range of datasets, including low-coverage (10x) datasets. VolcanoSV outperforms assembly-based tools in the identification of complex SVs, including translocations, duplications, and inversions, in both simulated and real cancer data. Moreover, VolcanoSV is robust to various evaluation parameters and accurately identifies breakpoints and SV sequences.

Parallel signatures of cognitive maturation in primate antisaccade performance and prefrontal activity.

The ability to suppress inappropriate actions and respond rapidly to appropriate ones matures late in life, after puberty. We investigated the development of this capability in monkeys trained to look away from a lone, bright stimulus (antisaccade task). We evaluated behavioral performance and recorded neural activity in the prefrontal cortex both before and after the transition from puberty to adulthood. Compared to when young, adult monkeys processed the stimulus more rapidly, resisted more effectively the involuntary urge to look at it, and adhered to the task rule more consistently. The spatially selective visuomotor neurons in the prefrontal cortex provided neural correlates of these behavioral changes indicative of a faster transition from stimulus-driven (exogenous) to goal-driven (endogenous) control within the time course of each trial. The results reveal parallel signatures of cognitive maturation in behavior and prefrontal activity that are consistent with improvements in attentional allocation after adolescence.

Tradeoffs in alignment and assembly-based methods for structural variant detection with long-read sequencing data.

Long-read sequencing offers long contiguous DNA fragments, facilitating diploid genome assembly and structural variant (SV) detection. Efficient and robust algorithms for SV identification are crucial with increasing data availability. Alignment-based methods, favored for their computational efficiency and lower coverage requirements, are prominent. Alternative approaches, relying solely on available reads for de novo genome assembly and employing assembly-based tools for SV detection via comparison to a reference genome, demand significantly more computational resources. However, the lack of comprehensive benchmarking constrains our comprehension and hampers further algorithm development. Here we systematically compare 14 read alignment-based SV calling methods (including 4 deep learning-based methods and 1 hybrid method), and 4 assembly-based SV calling methods, alongside 4 upstream aligners and 7 assemblers. Assembly-based tools excel in detecting large SVs, especially insertions, and exhibit robustness to evaluation parameter changes and coverage fluctuations. Conversely, alignment-based tools demonstrate superior genotyping accuracy at low sequencing coverage (5-10×) and excel in detecting complex SVs, like translocations, inversions, and duplications. Our evaluation provides performance insights, highlighting the absence of a universally superior tool. We furnish guidelines across 31 criteria combinations, aiding users in selecting the most suitable tools for diverse scenarios and offering directions for further method development.

Leveraging cross-source heterogeneity to improve the performance of bulk gene expression deconvolution.

A main limitation of bulk transcriptomic technologies is that individual measurements normally contain contributions from multiple cell populations, impeding the identification of cellular heterogeneity within diseased tissues. To extract cellular insights from existing large cohorts of bulk transcriptomic data, we present CSsingle, a novel method designed to accurately deconvolve bulk data into a predefined set of cell types using a scRNA-seq reference. Through comprehensive benchmark evaluations and analyses using diverse real data sets, we reveal the systematic bias inherent in existing methods, stemming from differences in cell size or library size. Our extensive experiments demonstrate that CSsingle exhibits superior accuracy and robustness compared to leading methods, particularly when dealing with bulk mixtures originating from cell types of markedly different cell sizes, as well as when handling bulk and single-cell reference data obtained from diverse sources. Our work provides an efficient and robust methodology for the integrated analysis of bulk and scRNA-seq data, facilitating various biological and clinical studies.

Benchmarking clustering, alignment, and integration methods for spatial transcriptomics.

BACKGROUND: Spatial transcriptomics (ST) is advancing our understanding of complex tissues and organisms. However, building a robust clustering algorithm to define spatially coherent regions in a single tissue slice and aligning or integrating multiple tissue slices originating from diverse sources for essential downstream analyses remains challenging. Numerous clustering, alignment, and integration methods have been specifically designed for ST data by leveraging its spatial information. The absence of comprehensive benchmark studies complicates the selection of methods and future method development. RESULTS: In this study, we systematically benchmark a variety of state-of-the-art algorithms with a wide range of real and simulated datasets of varying sizes, technologies, species, and complexity. We analyze the strengths and weaknesses of each method using diverse quantitative and qualitative metrics and analyses, including eight metrics for spatial clustering accuracy and contiguity, uniform manifold approximation and projection visualization, layer-wise and spot-to-spot alignment accuracy, and 3D reconstruction, which are designed to assess method performance as well as data quality. The code used for evaluation is available on our GitHub. Additionally, we provide online notebook tutorials and documentation to facilitate the reproduction of all benchmarking results and to support the study of new methods and new datasets. CONCLUSIONS: Our analyses lead to comprehensive recommendations that cover multiple aspects, helping users to select optimal tools for their specific needs and guide future method development.

Brain structure and activity predicting cognitive maturation in adolescence.

Cognitive abilities of primates, including humans, continue to improve through adolescence 1,2. While a range of changes in brain structure and connectivity have been documented 3,4, how they affect neuronal activity that ultimately determines performance of cognitive functions remains unknown. Here, we conducted a multilevel longitudinal study of monkey adolescent neurocognitive development. The developmental trajectory of neural activity in the prefrontal cortex accounted remarkably well for working memory improvements. While complex aspects of activity changed progressively during adolescence, such as the rotation of stimulus representation in multidimensional neuronal space, which has been implicated in cognitive flexibility, even simpler attributes, such as the baseline firing rate in the period preceding a stimulus appearance had predictive power over behavior. Unexpectedly, decreases in brain volume and thickness, which are widely thought to underlie cognitive changes in humans 5 did not predict well the trajectory of neural activity or cognitive performance changes. Whole brain cortical volume in particular, exhibited an increase and reached a local maximum in late adolescence, at a time of rapid behavioral improvement. Maturation of long-distance white matter tracts linking the frontal lobe with areas of the association cortex and subcortical regions best predicted changes in neuronal activity and behavior. Our results provide evidence that optimization of neural activity depending on widely distributed circuitry effects cognitive development in adolescence.

Laminar pattern of adolescent development changes in working memory neuronal activity.

Adolescent development is characterized by an improvement in cognitive abilities, such as working memory. Neurophysiological recordings in a non-human primate model of adolescence have revealed changes in neural activity that mirror improvement in behavior, including higher firing rate during the delay intervals of working memory tasks. The laminar distribution of these changes is unknown. By some accounts, persistent activity is more pronounced in superficial layers, so we sought to determine whether changes are most pronounced there. We therefore analyzed neurophysiological recordings from neurons recorded in the young and adult stage, at different cortical depths. Superficial layers exhibited increased baseline firing rate in the adult stage. Unexpectedly, changes in persistent activity were most pronounced in the middle layers. Finally, improved discriminability of stimulus location was most evident in the deeper layers. These results reveal the laminar pattern of neural activity maturation that is associated with cognitive improvement. NEW AND NOTEWORTHY: Structural brain changes are evident during adolescent development particularly in the cortical thickness of the prefrontal cortex, at a time when working memory ability increases markedly. The depth distribution of neurophysiological changes during adolescence is not known. Here we show that neurophysiological changes are not confined to superficial layers, which have most often been implicated in the maintenance of working memory. Contrary to expectations, greatest changes were evident in intermediate layers of the prefrontal cortex.

ADEPT: Autoencoder with differentially expressed genes and imputation for robust spatial transcriptomics clustering.

Advancements in spatial transcriptomics (ST) have enabled an in-depth understanding of complex tissues by quantifying gene expression at spatially localized spots. Several notable clustering methods have been introduced to utilize both spatial and transcriptional information in the analysis of ST datasets. However, data quality across different ST sequencing techniques and types of datasets influence the performance of different methods and benchmarks. To harness spatial context and transcriptional profile in ST data, we developed a graph-based, multi-stage framework for robust clustering, called ADEPT. To control and stabilize data quality, ADEPT relies on a graph autoencoder backbone and performs an iterative clustering on imputed, differentially expressed genes-based matrices to minimize the variance of clustering results. ADEPT outperformed other popular methods on ST data generated by different platforms across analyses such as spatial domain identification, visualization, spatial trajectory inference, and data denoising.