Differential expression (DE) analysis and gene set enrichment (GSE) analysis are commonly applied in single cell RNA sequencing (scRNA-seq) studies. Here, we develop an integrative and scalable computational method, iDEA, to perform joint DE and GSE analysis through a hierarchical Bayesian framework. By integrating DE and GSE analyses, iDEA can improve the power and consistency of DE analysis and the accuracy of GSE analysis. Importantly, iDEA uses only DE summary statistics as input, enabling effective data modeling through complementing and pairing with various existing DE methods. We illustrate the benefits of iDEA with extensive simulations. We also apply iDEA to analyze three scRNA-seq data sets, where iDEA achieves up to five-fold power gain over existing GSE methods and up to 64% power gain over existing DE methods. The power gain brought by iDEA allows us to identify many pathways that would not be identified by existing approaches in these data. <<<

Shaohua Hu, Hao Xu, Xiaohua Meng, Xiangxiang Bai, Junli Xu, Jinru Ji, Chaoqun Ying, Yunbo Chen, Ping Shen, Yunxiao Zhou, Beiwen Zheng, Yonghong Xiao <<<

is an emerging species and has increasingly been reported to cause life-threatening infections and even outbreaks in humans. Nevertheless, there is little data regarding the geographical distribution, phylogenetic structure, and transmission across the globe, especially in Asia. We utilize whole-genome sequencing (WGS) data to define a global population framework, phylogenetic structure, geographical distribution, and transmission evaluation of pathogens. The geographical distribution diagram revealed the emerging pathogenic bacteria already distributed in various countries worldwide, especially in the USA and China. Strikingly, phylogenetic analysis showed a part of our China original shared the same ancestor with the USA outbreak strain, which implies the possibility of localized outbreaks and global spread. These closer related strains also contained ICEEaI, which might insert into a disrupted DNA repair gene and made the strain more liable to mutation and outbreak infection. BEAST analysis showed that the most recent common ancestor for ICEEaI was dated twelve years ago, and China might be the most likely recent source of this bacteria. Our study sheds light on the potential possibility of causing the large-scale outbreak and rapid global dissemination. Continued genomic surveillance of the dynamics of populations will generate further knowledge for optimizing future prevent global outbreak infections. <<<

Maria Hurskainen, Ivana Mižíková, David P Cook, Noora Andersson, Chanèle Cyr-Depauw, Flore Lesage, Emmi Helle, Laurent Renesme, Robert P Jankov, Markku Heikinheimo, Barbara C Vanderhyden, Bernard Thébaud <<<

During late lung development, alveolar and microvascular development is finalized to enable sufficient gas exchange. Impaired late lung development manifests as bronchopulmonary dysplasia (BPD) in preterm infants. Single-cell RNA sequencing (scRNA-seq) allows for assessment of complex cellular dynamics during biological processes, such as development. Here, we use MULTI-seq to generate scRNA-seq profiles of over 66,000 cells from 36 mice during normal or impaired lung development secondary to hyperoxia with validation of some of the findings in lungs from BPD patients. We observe dynamic populations of cells, including several rare cell types and putative progenitors. Hyperoxia exposure, which mimics the BPD phenotype, alters the composition of all cellular compartments, particularly alveolar epithelium, stromal fibroblasts, capillary endothelium and macrophage populations. Pathway analysis and predicted dynamic cellular crosstalk suggest inflammatory signaling as the main driver of hyperoxia-induced changes. Our data provides a single-cell view of cellular changes associated with late lung development in health and disease. <<<

Joseph M. Replogle , Reuben A. Saunders , Angela N. Pogson , Jeffrey A. Hussmann , Alexander Lenail , Alina Guna , Lauren Mascibroda , Eric J. Wagner , Karen Adelman , Gila Lithwick-Yanai , Nika Iremadze , Florian Oberstrass , Doron Lipson , Jessica L. Bonnar , Marco Jost , Thomas M. Norman , Jonathan S. Weissman <<<

A central goal of genetics is to define the relationships between genotypes and phenotypes. High-content phenotypic screens such as Perturb-seq (CRISPR-based screens with single-cell RNA-sequencing readouts) enable massively parallel functional genomic mapping but, to date, have been used at limited scales. Here, we perform genome-scale Perturb-seq targeting all expressed genes with CRISPR interference (CRISPRi) across >2.5 million human cells. We use transcriptional phenotypes to predict the function of poorly characterized genes, uncovering new regulators of ribosome biogenesis (including CCDC86, ZNF236, and SPATA5L1), transcription (C7orf26), and mitochondrial respiration (TMEM242). In addition to assigning gene function, single-cell transcriptional phenotypes allow for in-depth dissection of complex cellular phenomena—from RNA processing to differentiation. We leverage this ability to systematically identify genetic drivers and consequences of aneuploidy and to discover an unanticipated layer of stress-specific regulation of the mitochondrial genome. Our information-rich genotype-phenotype map reveals a multidimensional portrait of gene and cellular function. <<<

Arash Jamshidi, Minetta C Liu, Eric A Klein, Oliver Venn, Earl Hubbell, John F Beausang, Samuel Gross, Collin Melton, Alexander P Fields, Qinwen Liu, Nan Zhang, Eric T Fung, Kathryn N Kurtzman, Hamed Amini, Craig Betts, Daniel Civello, Peter Freese, Robert Calef, Konstantin Davydov, Saniya Fayzullina, Chenlu Hou, Roger Jiang, Byoungsok Jung, Susan Tang, Vasiliki Demas, Joshua Newman, Onur Sakarya, Eric Scott, Archana Shenoy, Seyedmehdi Shojaee, Kristan K Steffen, Virgil Nicula, Tom C Chien, Siddhartha Bagaria, Nathan Hunkapiller, Mohini Desai, Zhao Dong, Donald A Richards, Timothy J Yeatman, Allen L Cohn, David D Thiel, Donald A Berry, Mohan K Tummala, Kristi McIntyre, Mikkael A Sekeres, Alan Bryce, Alexander M Aravanis, Michael V Seiden, Charles Swanton <<<

In the Circulating Cell-free Genome Atlas (NCT02889978) substudy 1, we evaluate several approaches for a circulating cell-free DNA (cfDNA)-based multi-cancer early detection (MCED) test by defining clinical limit of detection (LOD) based on circulating tumor allele fraction (cTAF), enabling performance comparisons. Among 10 machine-learning classifiers trained on the same samples and independently validated, when evaluated at 98% specificity, those using whole-genome (WG) methylation, single nucleotide variants with paired white blood cell background removal, and combined scores from classifiers evaluated in this study show the highest cancer signal detection sensitivities. Compared with clinical stage and tumor type, cTAF is a more significant predictor of classifier performance and may more closely reflect tumor biology. Clinical LODs mirror relative sensitivities for all approaches. The WG methylation feature best predicts cancer signal origin. WG methylation is the most promising technology for MCED and informs development of a targeted methylation MCED test. <<<

While spontaneous yawning is common across all vertebrate classes, contagious yawning is less common and has been observed only in a few species of social animals. Interspecific contagious yawning in response to yawning by humans has been observed only by chimpanzees and dogs. After confirming additional occurrences of intraspecific contagious yawning in a group of captive African elephants previously studied, we further investigated the potential for the same group of elephants to engage in interspecific contagious yawning with familiar human handlers. Ten captive African elephants, most of whom had been previously studied, were observed over 13 nights for evidence of intraspecific contagious yawning. Seven of these elephants were also involved in trials where familiar handlers performed staged yawns, as well as trials with staged non-yawning gapes, or trials with no yawns or gapes. Incorporating previously collected contagious yawning data, we describe nine instances of intraspecific contagious yawning in the elephants. Three of the seven elephants yawned contagiously in response to humans during the interspecific yawning trials. This is the first report of interspecific contagious yawning by elephants in response to yawns by familiar humans. <<<

Face processing supports our ability to recognize friend from foe, form tribes and understand the emotional implications of changes in facial musculature. This skill relies on a distributed network of brain regions, but how these regions interact is poorly understood. Here we integrate anatomical and functional connectivity measurements with behavioural assays to create a global model of the face connectome. We dissect key features, such as the network topology and fibre composition. We propose a neurocognitive model with three core streams; face processing along these streams occurs in a parallel and reciprocal manner. Although long-range fibre paths are important, the face network is dominated by short-range fibres. Finally, we provide evidence that the well-known right lateralization of face processing arises from imbalanced intra- and interhemispheric connections. In summary, the face network relies on dynamic communication across highly structured fibre tracts, enabling coherent face processing that underpins behaviour and cognition. <<<

Deep generative models have been an upsurge in the deep learning community since they were proposed. These models are designed for generating new synthetic data including images, videos and texts by fitting the data approximate distributions. In the last few years, deep generative models have shown superior performance in drug discovery especially de novo molecular design. In this study, deep generative models are reviewed to witness the recent advances of de novo molecular design for drug discovery. In addition, we divide those models into two categories based on molecular representations in silico. Then these two classical types of models are reported in detail and discussed about both pros and cons. We also indicate the current challenges in deep generative models for de novo molecular design. De novo molecular design automatically is promising but a long road to be explored. <<<

: The ability to engineer mammalian genomes in a quick and cost-effective way has led to rapid adaptation of CRISPR technology in biomedical research. CRISPR-based engineering has the potential to accelerate drug discovery, to support the reduction of high attrition rate in drug development and to enhance development of cell and gene-based therapies.: How CRISPR technology is transforming drug discovery is discussed in this review. From target identification to target validation in both and models, CRISPR technology is positively impacting the early stages of drug development by providing a straightforward way to genome engineering. This property also attracted attention for CRISPR application in the cell and gene therapy area.: CRISPR technology is rapidly becoming the preferred tool for genome engineering and nowadays it is hard to imagine the drug discovery pipeline without this technology. With the years to come, CRISPR technology will undoubtedly be further refined and will flourish into a mature technology that will play a key role in supporting genome engineering requirements in the drug discovery pipeline as well as in cell and gene therapy development. <<<

Deformable medical image registration is widely used in medical image processing with the invertible and one-to-one mapping between images. While state-of-the-art image registration methods are based on convolutional neural networks, few attempts have been made with Transformers which show impressive performance on computer vision tasks. Existing models neglect to employ attention mechanisms to handle the long-range cross-image relevance in embedding learning, limiting such approaches to identify the semantically meaningful correspondence of anatomical structures. These methods also ignore the topology preservation and invertibility of the transformation although they achieve fast image registration. In this paper, we propose a novel, symmetric unsupervised learning network Swin-VoxelMorph based on the Swin Transformer which minimizes the dissimilarity between images and estimates both forward and inverse transformations simultaneously. Specifically, we propose 3D Swin-UNet, which applies hierarchical Swin Transformer with shifted windows as the encoder to extract context features. And a symmetric Swin Transformer-based decoder with patch expanding layer is designed to perform the up-sampling operation to estimate the registration fields. Besides, our objective loss functions can guarantee substantial diffeomorphic properties of the predicted transformations. We verify our method on two datasets including ADNI and PPMI, and it achieves state-of-the-art registration accuracy while maintaining desirable diffeomorphic properties. <<<

Fengling Hu, Sarah M Weinstein, Erica B Baller, Alessandra M Valcarcel, Azeez Adebimpe, Armin Raznahan, David R Roalf, Timothy E Robert-Fitzgerald, Virgilio Gonzenbach, Ruben C Gur, Raquel E Gur, Simon Vandekar, John A Detre, Kristin A Linn, Aaron Alexander-Bloch, Theodore D Satterthwaite, Russell T Shinohara <<<

When individual subjects are imaged with multiple modalities, biological information is present not only within each modality, but also between modalities - that is, in how modalities covary at the voxel level. Previous studies have shown that local covariance structures between modalities, or intermodal coupling (IMCo), can be summarized for two modalities, and that two-modality IMCo reveals otherwise undiscovered patterns in neurodevelopment and certain diseases. However, previous IMCo methods are based on the slopes of local weighted linear regression lines, which are inherently asymmetric and limited to the two-modality setting. Here, we present a generalization of IMCo estimation which uses local covariance decompositions to define a symmetric, voxel-wise coupling coefficient that is valid for two or more modalities. We use this method to study coupling between cerebral blood flow, amplitude of low frequency fluctuations, and local connectivity in 803 subjects ages 8 through 22. We demonstrate that coupling is spatially heterogeneous, varies with respect to age and sex in neurodevelopment, and reveals patterns that are not present in individual modalities. As availability of multi-modal data continues to increase, principal-component-based IMCo (pIMCo) offers a powerful approach for summarizing relationships between multiple aspects of brain structure and function. An R package for estimating pIMCo is available at: https://github.com/hufengling/pIMCo. <<<

Since the 1990s, a novel third approach to vaccination against a broad array of target antigens and diseases has been in development. This technology involves the direct administration of plasmid deoxyribonucleic acid (DNA) containing the gene encoding the immunogen against which an immune response is sought, leading to the in-situ production of the target immunogen(s) in the vaccine recipient. Such vaccines are referred to as “plasmid DNA vaccines” or simply “DNA vaccines” (with both terms used interchangeably throughout these Guidelines). The value and advantages of plasmid DNA products need to be assessed on a case-by-case basis; their utility will depend upon: (a) the nature of the organism being vaccinated against or the targeted disease; (b) the nature of the immunogen or activity of the gene insert; (c) the type of immune response required for effectiveness; and (d) the delivery system and route of administration. The development and application of DNA vaccines continues to progress. Since the WHO Guidelines for assuring the quality and nonclinical safety evaluation of DNA vaccines were adopted by the Expert Committee for Biological Standardization in 2005, many clinical trials of DNA vaccines have taken place and considerable experience in their manufacture and control has accrued. The current revision reflects this experience, especially in relation to data derived from nonclinical and clinical safety testing, which address many of the concerns expressed in previous versions of these Guidelines. The control of DNA vaccines should continue to be approached in a flexible manner to enable further modifications as more experience is gained in their production and use, and as other components or delivery systems are included. The intention of the current document is to provide a scientifically sound basis for the consistent manufacture and control of DNA vaccines for human use to ensure their continued safety and efficacy following licensure. <<<

Betul Akgol Oksuz, Liyan Yang, Sameer Abraham, Sergey V Venev, Nils Krietenstein, Krishna Mohan Parsi, Hakan Ozadam, Marlies E Oomen, Ankita Nand, Hui Mao, Ryan M J Genga, Rene Maehr, Oliver J Rando, Leonid A Mirny, Johan H Gibcus, Job Dekker <<<

Chromosome conformation capture (3C) assays are used to map chromatin interactions genome-wide. Chromatin interaction maps provide insights into the spatial organization of chromosomes and the mechanisms by which they fold. Hi-C and Micro-C are widely used 3C protocols that differ in key experimental parameters including cross-linking chemistry and chromatin fragmentation strategy. To understand how the choice of experimental protocol determines the ability to detect and quantify aspects of chromosome folding we have performed a systematic evaluation of 3C experimental parameters. We identified optimal protocol variants for either loop or compartment detection, optimizing fragment size and cross-linking chemistry. We used this knowledge to develop a greatly improved Hi-C protocol (Hi-C 3.0) that can detect both loops and compartments relatively effectively. In addition to providing benchmarked protocols, this work produced ultra-deep chromatin interaction maps using Micro-C, conventional Hi-C and Hi-C 3.0 for key cell lines used by the 4D Nucleome project. <<<

Resting state functional magnetic resonance imaging (rsfMRI) data exhibits complex but structured patterns. However, the underlying origins are unclear and entangled in rsfMRI data. Here we establish a variational auto-encoder, as a generative model trainable with unsupervised learning, to disentangle the unknown sources of rsfMRI activity. After being trained with large data from the Human Connectome Project, the model has learned to represent and generate patterns of cortical activity and connectivity using latent variables. The latent representation and its trajectory represent the spatiotemporal characteristics of rsfMRI activity. The latent variables reflect the principal gradients of the latent trajectory and drive activity changes in cortical networks. Representational geometry captured as covariance or correlation between latent variables, rather than cortical connectivity, can be used as a more reliable feature to accurately identify subjects from a large group, even if only a short period of data is available in each subject. Our results demonstrate that VAE is a valuable addition to existing tools, particularly suited for unsupervised representation learning of resting state fMRI activity. <<<

SUMMARY: The visualization and analysis of genomic repeats is typically accomplished using dot plots; however, the emergence of telomere-to-telomere assemblies with multi-megabase repeats requires new visualization strategies. Here, we introduce StainedGlass, which can generate publication-quality figures and interactive visualizations that depict the identity and orientation of multi-megabase tandem repeat structures at a genome-wide scale. The tool can rapidly reveal higher-order structures and improve the inference of evolutionary history for some of the most complex regions of genomes.AVAILABILITY AND IMPLEMENTATION: StainedGlass is implemented using Snakemake and available open source under the MIT license at https://mrvollger.github.io/StainedGlass/.SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online. <<<

A free-energy principle has been proposed recently that accounts for action, perception and learning. This Review looks at some key brain theories in the biological (for example, neural Darwinism) and physical (for example, information theory and optimal control theory) sciences from the free-energy perspective. Crucially, one key theme runs through each of these theories - optimization. Furthermore, if we look closely at what is optimized, the same quantity keeps emerging, namely value (expected reward, expected utility) or its complement, surprise (prediction error, expected cost). This is the quantity that is optimized under the free-energy principle, which suggests that several global brain theories might be unified within a free-energy framework. <<<

Sheila Shanmugan, Jakob Seidlitz, Zaixu Cui, Azeez Adebimpe, Danielle S Bassett, Maxwell A Bertolero, Christos Davatzikos, Damien A Fair, Raquel E Gur, Ruben C Gur, Bart Larsen, Hongming Li, Adam Pines, Armin Raznahan, David R Roalf, Russell T Shinohara, Jacob Vogel, Daniel H Wolf, Yong Fan, Aaron Alexander-Bloch, Theodore D Satterthwaite <<<

Prior work has shown that there is substantial interindividual variation in the spatial distribution of functional networks across the cerebral cortex, or functional topography. However, it remains unknown whether there are sex differences in the topography of individualized networks in youth. Here, we leveraged an advanced machine learning method (sparsity-regularized non-negative matrix factorization) to define individualized functional networks in 693 youth (ages 8 to 23 y) who underwent functional MRI as part of the Philadelphia Neurodevelopmental Cohort. Multivariate pattern analysis using support vector machines classified participant sex based on functional topography with 82.9% accuracy ( < 0.0001). Brain regions most effective in classifying participant sex belonged to association networks, including the ventral attention, default mode, and frontoparietal networks. Mass univariate analyses using generalized additive models with penalized splines provided convergent results. Furthermore, transcriptomic data from the Allen Human Brain Atlas revealed that sex differences in multivariate patterns of functional topography were spatially correlated with the expression of genes on the X chromosome. These results highlight the role of sex as a biological variable in shaping functional topography. <<<

Causal discovery and causal reasoning are classically treated as separate and consecutive tasks: one first infers the causal graph, and then uses it to estimate causal effects of interventions. However, such a two-stage approach is uneconomical, especially in terms of actively collected interventional data, since the causal query of interest may not require a fully-specified causal model. From a Bayesian perspective, it is also unnatural, since a causal query (e.g., the causal graph or some causal effect) can be viewed as a latent quantity subject to posterior inference -- other unobserved quantities that are not of direct interest (e.g., the full causal model) ought to be marginalized out in this process and contribute to our epistemic uncertainty. In this work, we propose Active Bayesian Causal Inference (ABCI), a fully-Bayesian active learning framework for integrated causal discovery and reasoning, which jointly infers a posterior over causal models and queries of interest. In our approach to ABCI, we focus on the class of causally-sufficient, nonlinear additive noise models, which we model using Gaussian processes. We sequentially design experiments that are maximally informative about our target causal query, collect the corresponding interventional data, and update our beliefs to choose the next experiment. Through simulations, we demonstrate that our approach is more data-efficient than several baselines that only focus on learning the full causal graph. This allows us to accurately learn downstream causal queries from fewer samples while providing well-calibrated uncertainty estimates for the quantities of interest. <<<

Bayesian structure learning allows one to capture uncertainty over the causal directed acyclic graph (DAG) responsible for generating given data. In this work, we present Tractable Uncertainty for STructure learning (TRUST), a framework for approximate posterior inference that relies on probabilistic circuits as the representation of our posterior belief. In contrast to sample-based posterior approximations, our representation can capture a much richer space of DAGs, while also being able to tractably reason about the uncertainty through a range of useful inference queries. We empirically show how probabilistic circuits can be used as an augmented representation for structure learning methods, leading to improvement in both the quality of inferred structures and posterior uncertainty. Experimental results on conditional query answering further demonstrate the practical utility of the representational capacity of TRUST. <<<

This paper describes a path integral formulation of the free energy principle. The ensuing account expresses the paths or trajectories that a particle takes as it evolves over time. The main results are a method or principle of least action that can be used to emulate the behaviour of particles in open exchange with their external milieu. Particles are defined by a particular partition, in which internal states are individuated from external states by active and sensory blanket states. The variational principle at hand allows one to interpret internal dynamics - of certain kinds of particles - as inferring external states that are hidden behind blanket states. We consider different kinds of particles, and to what extent they can be imbued with an elementary form of inference or sentience. Specifically, we consider the distinction between dissipative and conservative particles, inert and active particles and, finally, ordinary and strange particles. Strange particles (look as if they) infer their own actions, endowing them with apparent autonomy or agency. In short - of the kinds of particles afforded by a particular partition - strange kinds may be apt for describing sentient behaviour. <<<