Previously we showed that a massively parallel reporter assay, mSTARR-seq, could be used to simultaneously test for both enhancer-like activity and DNA methylation-dependent enhancer activity for millions of loci in a single experiment (Lea et al., 2018). Here we apply mSTARR-seq to query nearly the entire human genome, including almost all CpG sites profiled either on the commonly used Illumina Infinium MethylationEPIC array or via reduced representation bisulfite sequencing. We show that fragments containing these sites are enriched for regulatory capacity, and that methylation-dependent regulatory activity is in turn sensitive to the cellular environment. In particular, regulatory responses to interferon alpha (IFNA) stimulation are strongly attenuated by methyl marks, indicating widespread DNA methylation-environment interactions. In agreement, methylation-dependent responses to IFNA identified via mSTARR-seq predict methylation-dependent transcriptional responses to challenge with influenza virus in human macrophages. Our observations support the idea that pre-existing DNA methylation patterns can influence the response to subsequent environmental exposures-one of the tenets of biological embedding. However, we also find that, on average, sites previously associated with early life adversity are not more likely to functionally influence gene regulation than expected by chance.

The tsunami of new multiplexed spatial profiling technologies has opened a range of computational challenges focused on leveraging these powerful data for biological discovery. A key challenge underlying computation is a suitable representation for features of cellular niches. Here, we develop the covariance environment (COVET), a representation that can capture the rich, continuous multivariate nature of cellular niches by capturing the gene-gene covariate structure across cells in the niche, which can reflect the cell-cell communication between them. We define a principled optimal transport-based distance metric between COVET niches and develop a computationally efficient approximation to this metric that can scale to millions of cells. Using COVET to encode spatial context, we develop environmental variational inference (ENVI), a conditional variational autoencoder that jointly embeds spatial and single-cell RNA-seq data into a latent space. Two distinct decoders either impute gene expression across spatial modality, or project spatial information onto dissociated single-cell data. We show that ENVI is not only superior in the imputation of gene expression but is also able to infer spatial context to disassociated single-cell genomics data.

The phase separation of intrinsically disordered proteins is emerging as an important mechanism for cellular organization. However, efforts to connect protein sequences to the physical properties of condensates, i.e., the molecular grammar, are hampered by a lack of effective approaches for probing high-resolution structural details. Using a combination of multiscale simulations and fluorescence lifetime imaging microscopy experiments, we systematically explored a series of systems consisting of diblock elastin-like polypeptides (ELP). The simulations succeeded in reproducing the variation of condensate stability upon amino acid substitution and revealed different microenvi-ronments within a single condensate, which we verified with environmentally sensitive fluorophores. The interspersion of hydrophilic and hydrophobic residues and a lack of secondary structure formation result in a frustrated microphase separation, which explains both the strong correlation between ELP condensate stability and interfacial hydrophobicity scales, as well as the prevalence of protein-water hydrogen bonds. Our study uncovers new mechanisms for condensate stability and organization that may be broadly applicable.

Shifts in microbiome community composition can have large effects on host health. It is therefore important to understand how perturbations, like those caused by the introduction of exogenous chemicals, modulate microbiome community composition. In poison frogs within the family Dendrobatidae, the skin microbiome is exposed to the alkaloids that the frogs sequester from their diet and use for defense. Given the demonstrated antimicrobial effects of these poison frog alkaloids, these compounds may be structuring the skin microbial community. To test this, we first characterized microbial communities from chemically defended and closely related non-defended frogs from Ecuador. Then we conducted a laboratory experiment to monitor the effect of the alkaloid decahydroquinoline (DHQ) on the microbiome of a single frog species. In both the field and lab experiments, we found that alkaloid-exposed microbiomes are more species rich and phylogenetically diverse, with an increase in rare taxa. To better understand the strain-specific behavior in response to alkaloids, we cultured microbial strains from poison frog skin and found the majority of strains exhibited either enhanced growth or were not impacted by the addition of DHQ. Additionally, stable isotope tracing coupled to nanoSIMS suggests that some of these strains are able to metabolize DHQ. Taken together, these data suggest that poison frog chemical defenses open new niches for skin-associated microbes with specific adaptations, including the likely metabolism of alkaloids, that enable their survival in this toxic environment. This work helps expand our understanding of how exposure to exogenous compounds like alkaloids can impact host microbiomes.

Goal-directed tasks involve acquiring an internal model, known as a predictive map, of relevant stimuli and associated outcomes to guide behavior. Here, we identified neural signatures of a predictive map of task behavior in perirhinal cortex (Prh). Mice learned to perform a tactile working memory task by classifying sequential whisker stimuli over multiple training stages. Chemogenetic inactivation demonstrated that Prh is involved in task learning. Chronic two-photon calcium imaging, population analysis, and computational modeling revealed that Prh encodes stimulus features as sensory prediction errors. Prh forms stable stimulus-outcome associations that expand in a retrospective manner and generalize as animals learn new contingencies. Stimulus-outcome associations are linked to prospective network activity encoding possible expected outcomes. This link is mediated by cholinergic signaling to guide task performance, demonstrated by acetylcholine imaging and perturbation. We propose that Prh combines error-driven and map-like properties to acquire a predictive map of learned task behavior.

The number and escape levels of genes that escape X chromosome inactivation (XCI) in female somatic cells vary among tissues and cell types, potentially contributing to specific sex differences. Here we investigate the role of CTCF, a master chromatin conformation regulator, in regulating escape from XCI. CTCF binding profiles and epigenetic features were systematically examined at constitutive and facultative escape genes using mouse allelic systems to distinguish the inactive X (Xi) and active X (Xa) chromosomes.

Dopamine in the nucleus accumbens helps motivate behavior based on expectations of future reward ("values"). These values need to be updated by experience: after receiving reward, the choices that led to reward should be assigned greater value. There are multiple theoretical proposals for how this credit assignment could be achieved, but the specific algorithms that generate updated dopamine signals remain uncertain. We monitored accumbens dopamine as freely behaving rats foraged for rewards in a complex, changing environment. We observed brief pulses of dopamine both when rats received reward (scaling with prediction error), and when they encountered novel path opportunities. Furthermore, dopamine ramped up as rats ran towards reward ports, in proportion to the value at each location. By examining the evolution of these dopamine place-value signals, we found evidence for two distinct update processes: progressive propagation along taken paths, as in temporal-difference learning, and inference of value throughout the maze, using internal models. Our results demonstrate that within rich, naturalistic environments dopamine conveys place values that are updated via multiple, complementary learning algorithms.

Pacinian corpuscle neurons are specialized low-threshold mechanoreceptors (LTMRs) that are tuned to detect high-frequency vibration (~40-2000 Hz), however it is unclear how Pacinians and other LTMRs encode mechanical forces encountered during naturalistic behavior. Here, we developed methods to record LTMRs in awake, freely moving mice. We find that Pacinians, but not other LTMRs, encode subtle vibrations of surfaces encountered by the animal, including low-amplitude vibrations initiated over two meters away. Strikingly, Pacinians are also highly active during a wide variety of natural behaviors, including walking, grooming, digging, and climbing. Pacinians in the hindlimb are sensitive enough to be activated by forelimb- or upper-body-dominant behaviors. Finally, we find that Pacinian LTMRs have diverse tuning and sensitivity. Our findings suggest a Pacinian population code for the representation of vibro-tactile features generated by self-initiated movements and low-amplitude environmental vibrations emanating from distant locations.

AGAP1 is an Arf1 GAP that regulates endolysosomal trafficking. Damaging variants have been linked to cerebral palsy and autism. We report 3 new individuals with microdeletion variants in AGAP1 . Affected individuals have intellectual disability (3/3), autism (3/3), dystonia with axial hypotonia (1/3), abnormalities of brain maturation (1/3), growth impairment (2/3) and facial dysmorphism (2/3). We investigated mechanisms potentially underlying AGAP1 neurodevelopmental impairments using the Drosophila ortholog, CenG1a . We discovered reduced axon terminal size, increased neuronal endosome abundance, and elevated autophagy at baseline. Given potential incomplete penetrance, we assessed gene-environment interactions. We found basal elevation in phosphorylation of the integrated stress-response protein eIF2α and inability to further increase eIF2α-P with subsequent cytotoxic stressors. CenG1a -mutant flies have increased lethality from exposure to environmental insults. We propose a model wherein disruption of AGAP1 function impairs endolysosomal trafficking, chronically activating the integrated stress response, and leaving AGAP1-deficient cells susceptible to a variety of second hit cytotoxic stressors. This model may have broader applicability beyond AGAP1 in instances where both genetic and environmental insults co-occur in individuals with neurodevelopmental disorders.

There is a growing awareness that tumor-adjacent normal tissues used as control samples in cancer studies do not represent fully healthy tissues. Instead, they are intermediates between healthy tissues and tumors. The factors that contribute to the deviation of such control samples from healthy state include exposure to the tumor-promoting factors, tumor-related immune response, and other aspects of tumor microenvironment. Characterizing the relation between gene expression of tumor-adjacent control samples and tumors is fundamental for understanding roles of microenvironment in tumor initiation and progression, as well as for identification of diagnostic and prognostic biomarkers for cancers. To address the demand, we developed and validated TranNet, a computational approach that utilizes gene expression in matched control and tumor samples to study the relation between their gene expression profiles. TranNet infers a sparse weighted bipartite graph from gene expression profiles of matched control samples to tumors. The results allow us to identify predictors (potential regulators) of this transition. To our knowledge, TranNet is the first computational method to infer such regulation. We applied TranNet to the data of several cancer types and their matched control samples from The Cancer Genome Atlas (TCGA). Many predictors identified by TranNet are genes associated with regulation by the tumor microenvironment as they are enriched in G-protein coupled receptor signaling, cell-to-cell communication, immune processes, and cell adhesion. Correspondingly, targets of inferred predictors are enriched in pathways related to tissue remodelling (including the epithelial-mesenchymal Transition (EMT)), immune response, and cell proliferation. This implies that the predictors are markers and potential stromal facilitators of tumor progression. Our results provide new insights for the relationships between tumor adjacent control sample, tumor and the tumor environment. Moreover, the set of predictors identified by TranNet will provide a valuable resource for future investigations. The TranNet method was implemented in python, source codes and the data sets used for and generated during this study are available at the Github site https://github.com/ncbi/TranNet .

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LTP2 greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs, and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Genetically heterogenous cell lines from laboratory mice are promising tools for population-based screening as they offer power for genetic mapping, and potentially, predictive value for in vivo experimentation in genetically matched individuals. To explore this further, we derived a panel of fibroblast lines from a genetic reference population of laboratory mice (the Diversity Outbred, DO). We then used high-content imaging to capture hundreds of cell morphology traits in cells exposed to the oxidative stress-inducing arsenic metabolite monomethylarsonous acid (MMAIII). We employed dose-response modeling to capture latent parameters of response and we then used these parameters to identify several hundred cell morphology quantitative trait loci (cmQTL). Response cmQTL encompass genes with established associations with cellular responses to arsenic exposure, including Abcc4 and Txnrd1, as well as novel gene candidates like Xrcc2. Moreover, baseline trait cmQTL highlight the influence of natural variation on fundamental aspects of nuclear morphology. We show that the natural variants influencing response include both coding and non-coding variation, and that cmQTL haplotypes can be used to predict response in orthogonal cell lines. Our study sheds light on the major molecular initiating events of oxidative stress that are under genetic regulation, including the NRF2-mediated antioxidant response, cellular detoxification pathways, DNA damage repair response, and cell death trajectories.

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation (CNV) provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene over-expression (OE), revealing that the fitness costs of CNV can vary substantially with genetic background in a common-garden environment. But the interplay between CNV tolerance and environment remains unexplored on a genomic scale. Here we measured the tolerance to gene OE in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride (NaCl) stress. OE genes that are commonly deleterious during NaCl stress recapitulated those commonly deleterious under standard conditions. However, NaCl stress uncovered novel differences in strain responses to gene OE. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to NaCl stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to gene OE. Although most genes were deleterious, hundreds were beneficial when overexpressed - remarkably, most of these effects were strain specific. Few beneficial genes were shared between the NaCl-sensitive isolates, implicating mechanistic differences behind their NaCl sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural CNV of a sodium export pump that likely contributes to strain-specific responses to OE of other genes. Our results reveal extensive strain-by-environment interaction in the response to gene CNV, raising important implications for the accessibility of CNV-dependent evolutionary routes under times of stress.

Targeting deep brain structures during electrophysiology and injections requires intensive training and expertise. Even with experience, researchers often can't be certain that a probe is placed precisely in a target location and this complexity scales with the number of simultaneous probes used in an experiment. Here, we present Pinpoint, open-source software that allows for interactive exploration of stereotaxic insertion plans. Once an insertion plan is created, Pinpoint allows users to save these online and share them with collaborators. 3D modeling tools allow users to explore their insertions alongside rig and implant hardware and ensure plans are physically possible. Probes in Pinpoint can be linked to electronic micro-manipulators allowing real-time visualization of current brain region targets alongside neural data. In addition, Pinpoint can control manipulators to automate and parallelize the insertion process. Compared to previously available software, Pinpoint's easy access through web browsers, extensive features, and real-time experiment integration enable more efficient and reproducible recordings.

The role of canines in transmitting antibiotic resistant bacteria to humans in the urban environment is poorly understood. To elucidate this role, we utilized genomic sequencing and phylogenetics to characterize the burden and transmission dynamics of antibiotic resistant Escherichia coli (ABR-Ec) cultured from canine and human feces present on urban sidewalks in San Francisco, California. We collected a total of fifty-nine ABR-Ec from human (n=12) and canine (n=47) fecal samples from the Tenderloin and South of Market (SoMa) neighborhoods of San Francisco. We then analyzed phenotypic and genotypic antibiotic resistance (ABR) of the isolates, as well as clonal relationships based on cgMLST and single nucleotide polymorphisms (SNPs) of the core genomes. Using Bayesian inference, we reconstructed the transmission dynamics between humans and canines from multiple local outbreak clusters using the marginal structured coalescent approximation (MASCOT). Overall, we found human and canine samples to carry similar amounts and profiles of ABR genes. Our results provide evidence for multiple transmission events of ABR-Ec between humans and canines. In particular, we found one instance of likely transmission from canines to humans as well as an additional local outbreak cluster consisting of one canine and one human sample. Based on this analysis, it appears that canine feces act as an important reservoir of clinically relevant ABR-Ec within the urban environment. Our findings support that public health measures should continue to emphasize proper canine feces disposal practices, access to public toilets and sidewalk and street cleaning. Importance: Antibiotic resistance in E. coli is a growing public health concern with global attributable deaths projected to reach millions annually. Current research has focused heavily on clinical routes of antibiotic resistance transmission to design interventions while the role of alternative reservoirs such as domesticated animals remain less well understood. Our results suggest canines are part of the transmission network that disseminates high-risk multidrug resistance in E. coli within the urban San Francisco community. As such, this study highlights the need to consider canines, and potentially domesticated animals more broadly, when designing interventions to reduce the prevalence of antibiotic resistance in the community. Additionally, it showcases the utility of genomic epidemiology to reconstruct the pathways by which antimicrobial resistance spreads.

The effects of context on the subjective experience of serotonergic psychedelics have not been fully examined in human neuroimaging studies, partly due to limitations of the imaging environment. Here, we administered saline or psilocybin to mice in their home cage or an enriched environment, immunofluorescently-labeled brain-wide c-Fos, and imaged cleared tissue with light sheet microscopy to examine the impact of context on psilocybin-elicited neural activity at cellular resolution. Voxel-wise analysis of c-Fos-immunofluorescence revealed differential neural activity, which we validated with c-Fos + cell density measurements. Psilocybin increased c-Fos expression in the neocortex, caudoputamen, central amygdala, and parasubthalamic nucleus and decreased c-Fos in the hypothalamus, cortical amygdala, striatum, and pallidum. Main effects of context and psilocybin-treatment were robust, widespread, and spatially distinct, whereas interactions were surprisingly sparse.

Understanding the relationship between the composition of the human gut microbiota and the ecological forces shaping it is of high importance as progress towards therapeutic modulation of the microbiota advances. However, given the inaccessibility of the gastrointestinal tract, our knowledge of the biogeographical and ecological relationships between physically interacting taxa has been limited to date. It has been suggested that interbacterial antagonism plays an important role in gut community dynamics, but in practice the conditions under which antagonistic behavior is favored or disfavored by selection in the gut environment are not well known. Here, using phylogenomics of bacterial isolate genomes and analysis of infant and adult fecal metagenomes, we show that the contact-dependent type VI secretion system (T6SS) is repeatedly lost from the genomes of Bacteroides fragilis in adults compare to infants. Although this result implies a significant fitness cost to the T6SS, but we could not identify in vitro conditions under which such a cost manifests. Strikingly, however, experiments in mice illustrated that the B. fragilis T6SS can be favored or disfavored in the gut environment, depending on the strains and species in the surrounding community and their susceptibility to T6SS antagonism. We use a variety of ecological modeling techniques to explore the possible local community structuring conditions that could underlie the results of our larger scale phylogenomic and mouse gut experimental approaches. The models illustrate robustly that the pattern of local community structuring in space can modulate the extent of interactions between T6SS-producing, sensitive, and resistant bacteria, which in turn control the balance of fitness costs and benefits of performing contact-dependent antagonistic behavior. Taken together, our genomic analyses, in vivo studies, and ecological theory point toward new integrative models for interrogating the evolutionary dynamics of type VI secretion and other predominant modes of antagonistic interaction in diverse microbiomes.

Xbp1 splicing and regulated IRE1-dependent RNA decay (RIDD) are two RNase activities of the ER stress sensor IRE1. While Xbp1 splicing has important roles in stress responses and animal physiology, the physiological role(s) of RIDD remain enigmatic. Genetic evidence in C. elegans connects XBP1-independent IRE1 activity to organismal stress adaptation, but whether this is via RIDD, and what are the targets is yet unknown. We show that cytosolic kinase/RNase domain of C. elegans IRE1 is indeed capable of RIDD in human cells, and that sensory neurons use RIDD to signal environmental stress, by degrading mRNA of TGFβ-like growth factor DAF-7. daf-7 was degraded in human cells by both human and worm IRE1 RNAse activity with same efficiency and specificity as Blos1, confirming daf-7 as RIDD substrate. Surprisingly, daf-7 degradation in vivo was triggered by concentrations of ER stressor tunicamycin too low for xbp-1 splicing. Decrease in DAF-7 normally signals food limitation and harsh environment, triggering adaptive changes to promote population survival. Because C. elegans is a bacteriovore, and tunicamycin, like other common ER stressors, is an antibiotic secreted by Streptomyces spp., we asked whether daf-7 degradation by RIDD could signal pending food deprivation. Indeed, pre-emptive tunicamycin exposure increased survival of C. elegans populations under food limiting/high temperature stress, and this protection was abrogated by overexpression of DAF-7. Thus, C. elegans uses stress-inducing metabolites in its environment as danger signals, and employs IRE1's RIDD activity to modulate the neuroendocrine signaling for survival of upcoming environmental challenge.

Environmental influences on brain structure and function during early development have been well-characterized. In pre-registered analyses, we test the theory that socioeconomic status (SES) is associated with differences in trajectories of intrinsic brain network development from birth to three years (n = 261). Prenatal SES is associated with developmental increases in cortical network segregation, with neonates and toddlers from lower-SES backgrounds showing a steeper increase in cortical network segregation with age, consistent with accelerated network development. Associations between SES and cortical network segregation occur at the local scale and conform to a sensorimotor-association hierarchy of cortical organization. SES-associated differences in cortical network segregation are associated with language abilities at two years, such that lower segregation is associated with improved language abilities. These results yield key insight into the timing and directionality of associations between the early environment and trajectories of cortical development.

New approaches for combatting microbial infections are needed. One strategy for disrupting pathogenesis involves developing compounds that interfere with bacterial virulence. A critical molecular determinant of virulence for Gram-negative bacteria are efflux pumps of the resistance-nodulation-division (RND) family, which includes AcrAB-TolC. We previously identified small molecules that bind AcrB, inhibit AcrAB-TolC, and do not appear to damage membranes. These efflux pump modulators (EPMs) were discovered in an in-cell screening platform called SAFIRE (Screen for Anti-infectives using Fluorescence microscopy of IntracellulaR Enterobacteriaceae). SAFIRE identifies compounds that disrupt the growth of a Gram-negative human pathogen, Salmonella enterica serotype Typhimurium (S. Typhimurium) in macrophages. We used medicinal chemistry to iteratively design ~200 EPM35 analogs and test them for activity in SAFIRE, generating compounds with nanomolar potency. Analogs were demonstrated to bind AcrB in a substrate binding pocket by cryo-electron microscopy (cryo-EM). Despite having amphipathic structures, the EPM analogs do not disrupt membrane voltage, as monitored by FtsZ localization to the cell septum. The EPM analogs had little effect on bacterial growth in standard Mueller Hinton Broth. However, under broth conditions that mimic the micro-environment of the macrophage phagosome, acrAB is required for growth, the EPM analogs are bacteriostatic, and increase the potency of antibiotics. These data suggest that under macrophage-like conditions the EPM analogs prevent the export of a toxic bacterial metabolite(s) through AcrAB-TolC. Thus, compounds that bind AcrB could disrupt infection by specifically interfering with the export of bacterial toxic metabolites, host defense factors, and/or antibiotics.