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. <<<

We present Momentum Contrast (MoCo) for unsupervised visual representation learning. From a perspective on contrastive learning as dictionary look-up, we build a dynamic dictionary with a queue and a moving-averaged encoder. This enables building a large and consistent dictionary on-the-fly that facilitates contrastive unsupervised learning. MoCo provides competitive results under the common linear protocol on ImageNet classification. More importantly, the representations learned by MoCo transfer well to downstream tasks. MoCo can outperform its supervised pre-training counterpart in 7 detection/segmentation tasks on PASCAL VOC, COCO, and other datasets, sometimes surpassing it by large margins. This suggests that the gap between unsupervised and supervised representation learning has been largely closed in many vision tasks. <<<

The history of human populations in Africa is complex and includes various demographic events that influenced patterns of genetic variation across the continent. Through genetic studies of modern-day, and most recently, ancient African genetic variation, it became evident that deep African history is captured by the relationships among hunter-gatherers. Furthermore, it was shown that agriculture had a large influence on the distribution of current-day Africans. These later population movements changed the demographic face of the continent and descendants of farming groups today form the majority populations across Africa. Ancient DNA methods are continually evolving, and we see evidence of this in how research has advanced in the last decade. With the increased availability of full genomic data from diverse sets of modern-day and prehistoric Africans we now have more power to infer human demography. Future ancient DNA research promises to reveal more detailed stories of human prehistory in Africa. <<<

Anti-PD-1-based immunotherapy has had a major impact on cancer treatment but has only benefited a subset of patients. Among the variables that could contribute to interpatient heterogeneity is differential composition of the patients' microbiome, which has been shown to affect antitumor immunity and immunotherapy efficacy in preclinical mouse models. We analyzed baseline stool samples from metastatic melanoma patients before immunotherapy treatment, through an integration of 16 ribosomal RNA gene sequencing, metagenomic shotgun sequencing, and quantitative polymerase chain reaction for selected bacteria. A significant association was observed between commensal microbial composition and clinical response. Bacterial species more abundant in responders included , , and Reconstitution of germ-free mice with fecal material from responding patients could lead to improved tumor control, augmented T cell responses, and greater efficacy of anti-PD-L1 therapy. Our results suggest that the commensal microbiome may have a mechanistic impact on antitumor immunity in human cancer patients. <<<

Xiangdong Pu, Zhen Li, Ya Tian, Ranran Gao, Lijun Hao, Yating Hu, Chunnian He, Wei Sun, Meimei Xu, Reuben J Peters, Yves Van de Peer, Zhichao Xu, Jingyuan Song <<<

Lonicera japonica is a widespread member of the Caprifoliaceae (honeysuckle) family utilized in traditional medical practices. This twining vine honeysuckle also is a much-sought ornamental, in part due to its dynamic flower coloration, which changes from white to gold during development. The molecular mechanism underlying dynamic flower coloration in L. japonica was elucidated by integrating whole genome sequencing, transcriptomic analysis and biochemical assays. Here, we report a chromosome-level genome assembly of L. japonica, comprising nine pseudochromosomes with a total size of 843.2 Mb. We also provide evidence for a whole-genome duplication event in the lineage leading to L. japonica, which occurred after its divergence from Dipsacales and Asterales. Moreover, gene expression analysis not only revealed correlated expression of the relevant biosynthetic genes with carotenoid accumulation, but also suggested a role for carotenoid degradation in L. japonica's dynamic flower coloration. The variation of flower color is consistent with not only the observed carotenoid accumulation pattern, but also with the release of volatile apocarotenoids that presumably serve as pollinator attractants. Beyond novel insights into the evolution and dynamics of flower coloration, the high-quality L. japonica genome sequence also provides a foundation for molecular breeding to improve desired characteristics. <<<

The aim of this paper is to provide a comprehensive review of statistical challenges in neuroimaging data analysis from neuroimaging techniques to large-scale neuroimaging studies to statistical learning methods. We briefly review eight popular neuroimaging techniques and their potential applications in neuroscience research and clinical translation. We delineate the four common themes of neuroimaging data and review major image processing analysis methods for processing neuroimaging data at the individual level. We briefly review four large-scale neuroimaging-related studies and a consortium on imaging genomics and discuss four common themes of neuroimaging data analysis at the population level. We review nine major population-based statistical analysis methods and their associated statistical challenges and present recent progress in statistical methodology to address these challenges. <<<

Tao Zhong , Jingkuan Wei , Kunhua Wu , Liangjun Chen , Fenqiang Zhao , Yuchen Pei , Ya Wang , Hongjiang Zhang , Zhengwang Wu , Ying Huang , Tengfei Li , Li Wang , Yongchang Chen , Weizhi Ji , Yu Zhang , Gang Li , Yuyu Niu <<<

Longitudinal brain imaging atlases with densely sampled time-points and ancillary anatomical information are of fundamental importance in studying early developmental characteristics of human and non-human primate brains during infancy, which feature extremely dynamic imaging appearance, brain shape and size. However, for non-human primates, which are highly valuable animal models for understanding human brains, the existing brain atlases are mainly developed based on adults or adolescents, denoting a notable lack of temporally densely-sampled atlases covering the dynamic early brain development. To fill this critical gap, in this paper, we construct a comprehensive set of longitudinal brain atlases and associated tissue probability maps (gray matter, white matter, and cerebrospinal fluid) with totally 12 time-points from birth to 4 years of age (i.e., 1, 2, 3, 4, 5, 6, 9, 12, 18, 24, 36, and 48 months of age) based on 175 longitudinal structural MRI scans from 39 typically-developing cynomolgus macaques, by leveraging state-of-the-art computational techniques tailored for early developing brains. Furthermore, to facilitate region-based analysis using our atlases, we also provide two popular hierarchy parcellations, i.e., cortical hierarchy maps (6 levels) and subcortical hierarchy maps (6 levels), on our longitudinal macaque brain atlases. These early developing atlases, which have the densest time-points during infancy (to the best of our knowledge), will greatly facilitate the studies of macaque brain development. <<<

The hallmarks of cancer conceptualization is a heuristic tool for distilling the vast complexity of cancer phenotypes and genotypes into a provisional set of underlying principles. As knowledge of cancer mechanisms has progressed, other facets of the disease have emerged as potential refinements. Herein, the prospect is raised that phenotypic plasticity and disrupted differentiation is a discrete hallmark capability, and that nonmutational epigenetic reprogramming and polymorphic microbiomes both constitute distinctive enabling characteristics that facilitate the acquisition of hallmark capabilities. Additionally, senescent cells, of varying origins, may be added to the roster of functionally important cell types in the tumor microenvironment. SIGNIFICANCE: Cancer is daunting in the breadth and scope of its diversity, spanning genetics, cell and tissue biology, pathology, and response to therapy. Ever more powerful experimental and computational tools and technologies are providing an avalanche of "big data" about the myriad manifestations of the diseases that cancer encompasses. The integrative concept embodied in the hallmarks of cancer is helping to distill this complexity into an increasingly logical science, and the provisional new dimensions presented in this perspective may add value to that endeavor, to more fully understand mechanisms of cancer development and malignant progression, and apply that knowledge to cancer medicine. <<<

Faces of different people elicit distinct fMRI patterns in several face-selective regions of the human brain. Here we used representational similarity analysis to investigate what type of identity-distinguishing information is encoded in three face-selective regions: fusiform face area (FFA), occipital face area (OFA), and posterior superior temporal sulcus (pSTS). In a sample of 30 human participants (22 females, 8 males), we used fMRI to measure brain activity patterns elicited by naturalistic videos of famous face identities, and compared their representational distances in each region with models of the differences between identities. We built diverse candidate models, ranging from low-level image-computable properties (pixel-wise, GIST, and Gabor-Jet dissimilarities), through higher-level image-computable descriptions (OpenFace deep neural network, trained to cluster faces by identity), to complex human-rated properties (perceived similarity, social traits, and gender). We found marked differences in the information represented by the FFA and OFA. Dissimilarities between face identities in FFA were accounted for by differences in perceived similarity, Social Traits, Gender, and by the OpenFace network. In contrast, representational distances in OFA were mainly driven by differences in low-level image-based properties (pixel-wise and Gabor-Jet dissimilarities). Our results suggest that, although FFA and OFA can both discriminate between identities, the FFA representation is further removed from the image, encoding higher-level perceptual and social face information. Recent studies using fMRI have shown that several face-responsive brain regions can distinguish between different face identities. It is however unclear whether these different face-responsive regions distinguish between identities in similar or different ways. We used representational similarity analysis to investigate the computations within three brain regions in response to naturalistically varying videos of face identities. Our results revealed that two regions, the fusiform face area and the occipital face area, encode distinct identity information about faces. Although identity can be decoded from both regions, identity representations in fusiform face area primarily contained information about social traits, gender, and high-level visual features, whereas occipital face area primarily represented lower-level image features. <<<

Nazli Dizman, Luis Meza, Paulo Bergerot, Marice Alcantara, Tanya Dorff, Yung Lyou, Paul Frankel, Yujie Cui, Valerie Mira, Marian Llamas, Joann Hsu, Zeynep Zengin, Nicholas Salgia, Sabrina Salgia, Jasnoor Malhotra, Neal Chawla, Alex Chehrazi-Raffle, Ramya Muddasani, John Gillece, Lauren Reining, Jeff Trent, Motomichi Takahashi, Kentaro Oka, Seiya Higashi, Marcin Kortylewski, Sarah K Highlander, Sumanta K Pal <<<

Previous studies have suggested that the gut microbiome influences the response to checkpoint inhibitors (CPIs) in patients with cancer. CBM588 is a bifidogenic live bacterial product that we postulated could augment CPI response through modulation of the gut microbiome. In this open-label, single-center study (NCT03829111), 30 treatment-naive patients with metastatic renal cell carcinoma with clear cell and/or sarcomatoid histology and intermediate- or poor-risk disease were randomized 2:1 to receive nivolumab and ipilimumab with or without daily oral CBM588, respectively. Stool metagenomic sequencing was performed at multiple timepoints. The primary endpoint to compare the relative abundance of Bifidobacterium spp. at baseline and at 12 weeks was not met, and no significant differences in Bifidobacterium spp. or Shannon index associated with the addition of CBM588 to nivolumab-ipilimumab were detected. Secondary endpoints included response rate, progression-free survival (PFS) and toxicity. PFS was significantly longer in patients receiving nivolumab-ipilimumab with CBM588 than without (12.7 months versus 2.5 months, hazard ratio 0.15, 95% confidence interval 0.05-0.47, P = 0.001). Although not statistically significant, the response rate was also higher in patients receiving CBM588 (58% versus 20%, P = 0.06). No significant difference in toxicity was observed between the study arms. The data suggest that CBM588 appears to enhance the clinical outcome in patients with metastatic renal cell carcinoma treated with nivolumab-ipilimumab. Larger studies are warranted to confirm this clinical observation and elucidate the mechanism of action and the effects on microbiome and immune compartments. <<<

Andreas Schuh , Antonios Makropoulos , Emma C. Robinson , Lucilio Cordero-Grande , Emer Hughes , Jana Hutter , Anthony N. Price , Maria Murgasova , Rui Pedro A. G. Teixeira , Nora Tusor , Johannes K. Steinweg , Suresh Victor , Mary A. Rutherford , Joseph V. Hajnal , A. David Edwards , Daniel Rueckert <<<

AbstractPremature birth increases the risk of developing neurocognitive and neurobe-havioural disorders. The mechanisms of altered brain development causing these disorders are yet unknown. Studying the morphology and function of the brain during maturation provides us not only with a better understanding of normal development, but may help us to identify causes of abnormal development and their consequences. A particular difficulty is to distinguish abnormal patterns of neurodevelopment from normal variation. The Developing Human Connectome Project (dHCP) seeks to create a detailed four-dimensional (4D) connectome of early life. This connectome may provide insights into normal as well as abnormal patterns of brain development. As part of this project, more than a thousand healthy fetal and neonatal brains will be scanned in vivo. This requires computational methods which scale well to larger data sets. We propose a novel groupwise method for the construction of a spatio-temporal model of mean morphology from cross-sectional brain scans at different gestational ages. This model scales linearly with the number of images and thus improves upon methods used to build existing public neonatal atlases, which derive correspondence between all pairs of images. By jointly estimating mean shape and longitudinal change, the atlas created with our method overcomes temporal inconsistencies, which are encountered when mean shape and intensity images are constructed separately for each time point. Using this approach, we have constructed a spatio-temporal atlas from 275 healthy neonates between 35 and 44 weeks post-menstrual age (PMA). The resulting atlas qualitatively preserves cortical details significantly better than publicly available atlases. This is moreover confirmed by a number of quantitative measures of the quality of the spatial normalisation and sharpness of the resulting template brain images. <<<

Adolescents' secrecy is intertwined with perception of parents' behaviors as acts of privacy invasion. It is currently untested, however, how this transactional process operates at the within-person level-where these causal processes take place. Dutch adolescents (n = 244, M = 13.84, 38.50% boys) reported three times on perceived parental privacy invasion and secrecy. Cross-lagged panel models (CLPM) confirmed earlier findings. Privacy invasion predicted increased secrecy, but a reverse effect was found from increased secrecy to increased privacy invasion. Controlling for confounding positive group-level associations with a novel random intercept CLPM, negative within-person associations were found. Higher levels of secrecy predicted lower levels of privacy invasive behaviors at the within-person level. These opposing findings within- versus between-persons illustrate a Simpson's paradox. <<<

Abhinav Achreja, Tao Yu, Anjali Mittal, Srinadh Choppara, Olamide Animasahun, Minal Nenwani, Fulei Wuchu, Noah Meurs, Aradhana Mohan, Jin Heon Jeon, Itisam Sarangi, Anusha Jayaraman, Sarah Owen, Reva Kulkarni, Michele Cusato, Frank Weinberg, Hye Kyong Kweon, Chitra Subramanian, Max S Wicha, Sofia D Merajver, Sunitha Nagrath, Kathleen R Cho, Analisa DiFeo, Xiongbin Lu, Deepak Nagrath <<<

Recurrent loss-of-function deletions cause frequent inactivation of tumour suppressor genes but often also involve the collateral deletion of essential genes in chromosomal proximity, engendering dependence on paralogues that maintain similar function. Although these paralogues are attractive anticancer targets, no methodology exists to uncover such collateral lethal genes. Here we report a framework for collateral lethal gene identification via metabolic fluxes, CLIM, and use it to reveal MTHFD2 as a collateral lethal gene in UQCR11-deleted ovarian tumours. We show that MTHFD2 has a non-canonical oxidative function to provide mitochondrial NAD, and demonstrate the regulation of systemic metabolic activity by the paralogue metabolic pathway maintaining metabolic flux compensation. This UQCR11-MTHFD2 collateral lethality is confirmed in vivo, with MTHFD2 inhibition leading to complete remission of UQCR11-deleted ovarian tumours. Using CLIM's machine learning and genome-scale metabolic flux analysis, we elucidate the broad efficacy of targeting MTHFD2 despite distinct cancer genetic profiles co-occurring with UQCR11 deletion and irrespective of stromal compositions of tumours. <<<

Species of the Enterobacter cloacae complex are widely encountered in nature, but they can act as pathogens. The biochemical and molecular studies on E. cloacae have shown genomic heterogeneity, comprising six species: Enterobacter cloacae, Enterobacter asburiae, Enterobacter hormaechei, Enterobacter kobei, Enterobacter ludwigii and Enterobacter nimipressuralis, E. cloacae and E. hormaechei are the most frequently isolated in human clinical specimens. Phenotypic identification of all species belonging to this taxon is usually difficult and not always reliable; therefore, molecular methods are often used. Although the E. cloacae complex strains are among the most common Enterobacter spp. causing nosocomial bloodstream infections in the last decade, little is known about their virulence-associated properties. By contrast, much has been published on the antibiotic-resistance features of these microorganisms. In fact, they are capable of overproducing AmpC β-lactamases by derepression of a chromosomal gene or by the acquisition of a transferable ampC gene on plasmids conferring the antibiotic resistance. Many other resistance determinants that are able to render ineffective almost all antibiotic families have been recently acquired. Most studies on antimicrobial susceptibility are focused on E. cloacae, E. hormaechei and E. asburiae; these studies reported small variations between the species, and the only significant differences had no discriminating features. <<<