title: Classification and genetic targeting of cell types in the primary taste and premotor center of the Drosophila brain

abstract: Neural circuits carry out complex computations that allow animals to evaluate food, select mates, move toward attractive stimuli, and move away from threats. In insects, the subesophageal zone (SEZ) is a brain region that receives gustatory, pheromonal, and mechanosensory inputs and contributes to the control of diverse behaviors, including feeding, grooming, and locomotion. Despite its importance in sensorimotor transformations, the study of SEZ circuits has been hindered by limited knowledge of the underlying diversity of SEZ neurons. Here, we generate a collection of split-GAL4 lines that provides precise genetic targeting of 138 different SEZ cell types in D. melanogaster, comprising approximately one third of all SEZ neurons. We characterize the single cell anatomy of these neurons and find that they cluster by morphology into six supergroups that organize the SEZ into discrete anatomical domains. We find that the majority of local SEZ interneurons are not classically polarized, suggesting rich local processing, whereas SEZ projection neurons tend to be classically polarized, conveying information to a limited number of higher brain regions. This study provides insight into the anatomical organization of the SEZ and generates resources that will facilitate further study of SEZ neurons and their contributions to sensory processing and behavior.


title: Thalamo-hippocampal pathway regulates incidental memory load: insights from sex differences

abstract:
Memory can be challenged by increasing both its required duration and the amount of information to be encoded, namely the memory load. The dorsal hippocampus (dHP) is involved in memory consolidation, which is the stabilization of a trace from short-term (STM) to long-term memory (LTM), as well as in the ability to process high information load. However, how memory load influences memory consolidation, and the underlying neural mechanisms, are yet unknown.

To address these questions, we used male and female mice that, despite having the same STM capacity of 6 objects in our Different Object recognition Task (DOT), spontaneously show differences in the number of objects directly transferred to LTM, when tested over longer delays. Males memorize all 6 encoded objects, while females only remember up to 4, both at 1 and 24 h delays. Interestingly, males activate the dHP more (as measured by c-Fos expression), whereas females activate the nuclei reuniens and rhomboid of the ventral midline thalamus (VMT). Optogenetic inhibition of the VMT-dHP pathway during off-line memory consolidation favors 6-object LTM retention in females by removing inhibitory control over dHP activation, while chemogenetic VMT-activation impairs it in males.

Our data represent a first demonstration of a sub-cortical control of dHP recruitment, that might underlie its sex-dependent activation during incidental memory, with potential also for clinical application.

One Sentence SummaryThe sex-dependent recruitment of a subcortical pathway regulates how many items are spontaneously memorized in mice.

category: neuroscience


title: Identification of a developmental switch in information transfer between whisker S1 and S2 cortex in mice

abstract: The whiskers of rodents are a key sensory organ that provides critical tactile information for animal navigation and object exploration throughout life. Previous work has explored the developmental sensory-driven activation of the primary sensory cortex processing whisker information (wS1), also called barrel cortex. This body of work has shown that the barrel cortex is already activated by sensory stimuli during the first post-natal week. However, it is currently unknown when over the course of development these stimuli begin being processed by higher order cortical areas, such as secondary whisker somatosensory area (wS2). Here we investigate for the first time the developmental engagement of wS2 by sensory stimuli and the emergence of cortico-cortical communication from wS1 to wS2. Using in vivo wide-field imaging and electrophysiological recordings in control and conditional knock-out mice we find that wS1 and wS2 are able to process bottom-up information coming from the thalamus already right after birth. We identify that it is only at the end of the first post-natal week that wS1 begins to provide excitation into wS2, a connection which begins to acquire feed-forward inhibition characteristics after the second post-natal week. Therefore, we have uncovered a developmental window during which excitatory versus inhibitory functional connectivity between wS1 and wS2 takes place.


title: Bispectrum-based Cross-frequency Functional Connectivity: A Study of Alzheimer's Disease

abstract: Alzheimers disease (AD) is a neurodegenerative disorder known to affect functional connectivity (FC) across many brain regions. Linear FC measures have been applied to study the differences in AD by splitting neurophysiological signals, such as electroencephalography (EEG) recordings, into discrete frequency bands and analysing them in isolation from each other. We address this limitation by quantifying cross-frequency FC in addition to the traditional within-band approach. Cross-bispectrum, a higher-order spectral analysis approach, is used to measure the nonlinear FC and is compared with the cross-spectrum, which only measures the linear FC within bands. This work reports the first use of cross-bispectrum to reconstruct a cross-frequency FC network where each frequency band is treated as a layer in a multilayer network with both inter- and intra-layer edges. An increase of within-band FC in AD cases is observed in low-frequency bands using both methods. Cross-bispectrum also detects multiple cross-frequency differences, mainly increased FC in AD cases in{delta} -{theta} coupling. Overall, increased strength of low-frequency coupling and decreased level of high-frequency coupling is observed in AD cases in comparison to healthy controls (HC). Vulnerability analysis reveals that the integration and segregation properties of networks are enabled by different frequency couplings in AD networks compared to HCs. Moreover, the segregation property depends less on the{gamma} interactions in AD cases, suggesting a shift from high to low-frequency connectivity. The results highlight the importance of studying nonlinearity and including cross-frequency FC in characterising AD. Moreover, the results demonstrate the advantages and limitations of using cross-bispectrum to reconstruct FC networks.


title: Cross-Frequency Multilayer Network Analysis with Bispectrum-based Functional Connectivity: A Study of Alzheimer's Disease


title: Thalamocortical and intracortical contributions to stimulus-evoked and oscillatory activity in rodent primary auditory cortex

abstract: Intracortical, horizontal connections seem ideally suited to contribute to cortical processing by spreading information across cortical space and coordinating activity between distant cortical sites. In sensory systems experiments have implicated horizontal connections in the generation of receptive fields and have in turn led to computational models of receptive field generation that rely on the contribution of horizontal connections. Testing the contribution of horizontal connections at the mesoscopic level has been difficult due to the lack of a suitable method to observe the activity of intracortical horizontal connections. Here, we develop such a method based on the analysis of the relative residues of the cortical laminar current source density reconstructions. In the auditory cortex of Mongolian gerbils, the method is then tested by manipulating the contribution of horizontal connections by surgical dissection. Our results indicate that intracortical horizontal connections contribute to the frequency-tuning of mesoscopic cortical patches. Futhermore, we dissociated a type of cortical gamma oscillation based on horizontal connections between mesoscopic patches from gamma oscillations locally generated within mesoscopic patches. The data further imply that global and local coordination of activity during sensory stimulation occur in a low and high gamma frequency band, respectively. Taken together the present data demonstrate that intracortical horizontal connections play an important role in generating cortical feature tuning and coordinate neuronal oscillations across cortex.


title: Top-down feedback enables flexible coding strategies in olfactory cortex

abstract:
In chemical sensation, multiple models have been proposed to explain how odors are represented by patterns of neuronal activity in the olfactory cortex. One hypothesis is that the identity of combinations of active neurons within specific sniff-related time windows are critical for encoding information about odors. Another model is that patterns of neural activity evolve across time and it is this temporal structure that is essential for encoding odor information. Interestingly, we found that top-down feedback to the olfactory bulb dictates what information is transmitted to the olfactory cortex by switching between these two strategies. Using a detailed model of the early olfactory system, we demonstrate that feedback control of inhibitory granule cells in the main olfactory bulb influences the balance between excitatory and inhibitory synaptic currents in mitral cells, thereby restructuring the firing patterns of piriform cortical cells across time. This resulted in performance gains in both the accuracy and reaction time of odor discrimination tasks. These findings lead us to propose a new framework for early olfactory computation, one in which top-down feedback to the bulb flexibly controls the temporal structure of neural activity in olfactory cortex, allowing the early olfactory system to dynamically switch between two distinct models of coding.

HighlightsO_LICentrifugal feedback shapes the temporal structure of neuronal firing in piriform cortical cells
C_LIO_LIFeedback controls information to piriform cortex by restructuring the ratio of excitatory and inhibitory synaptic inputs in the bulb
C_LIO_LICentrifugal feedback restructures how identity and timing of glomerular activity is represented in temporal patterns of activity in piriform cortex
C_LIO_LITemporal information improves behavioral performance in accuracy and reaction time of odor discrimination
C_LI


title: Ladostigil attenuates the oxidative and ER stress in human neuroblast-like SH-SY5Y cells

abstract: A hallmark of the aging brain is the robust inflammation mediated by microglial activation. Neuroinflammation resulting from the induction of oxidative stress in neurodegenerative diseases and following brain injury. Chronic treatment of aging rats by ladostigil, a compound with antioxidant and anti-inflammatory function, prevented microglial activation and learning deficits. In this study, we investigate the effect of ladostigil on neuronal-like SH-SY5Y cells. We show that SH-SY5Y cells exposed to acute (by H2O2) or chronic oxidative stress (by Sin1, 3-morpholinosydnonimine) induced apoptotic cell death. However, in the presence of ladostigil, the decline in cell viability and the oxidative levels were partially reversed. RNA-seq analysis showed that chronic oxidation by Sin1 resulted in coordinated suppression of endoplasmic reticulum (ER) quality control and ER stress response gene sets. Chronic oxidative stress impacted ER proteostasis and induced the expression of numerous lncRNAs. Pre-incubation with ladostigil before exposing SH-SY5Y cells to Sin1 induced Clk1 (Cdc2-like kinase 1) which was implicated in psychophysiological stress in mice and Alzheimer disease. Ladostigil also suppressed the expression of Ccpg1 (Cell cycle progression 1) and Synj1 (Synaptojanin 1) that function in ER-autophagy and endocytic pathways. We postulate that ladostigil alleviated cell damage by oxidation and ER stress. Therefore, it may attenuate neurotoxicity and cell death that accompany chronic stress conditions in the aging brain.

category: cell biology


title: Phyllotaxis Models: from the Inhibition Potential to the Real Plant

abstract: Previous phyllotaxis models allowed the initiation of new primordia when a threshold of inhibition potential is reached on the meristem front: their adequacy to botanical reality is only qualitative. We formulated the hypothesis that it is not the value of the inhibition threshold that remains constant as the meristem develops, but the difference of the inhibition thresholds during the initiation of two successive primordia. We were thus able to model with accuracy the sequence of plastochron ratios observed by Williams (1975) on the leaf meristem of flax: an outstanding result. More generally, we have shown that the evolution trajectories of the phyllotaxis modes as a function of the plastochron ratios follow the minima of the potential under decreasing plastochron ratios constraint and bifurcate when the number of these minima increases, thus giving physicochemical foundations to the famous van Iterson diagram. This historical representation of rising phyllotaxis shows the trajectories, but doesnt give the velocity of the movement: our plastochron ratio sequence adds this major dynamical information.

category: developmental biology


title: Regeneration of the larval sea star nervous system by wounding induced respecification to the sox2 lineage

abstract: The ability to restore lost body parts following traumatic injury is a fascinating area of biology that challenges current understanding of the ontogeny of differentiation. The origin of new cells needed to regenerate lost tissue, and whether they are pluripotent stem cells, tissue-specific stem cells or have de- or trans- differentiated, remains one of the most important open questions in regeneration. Additionally, it is not clearly known whether developmental gene regulatory networks (GRNs) are reused to direct specification in these cells or whether regeneration specific networks are deployed. Echinoderms, including sea stars, have extensive ability for regeneration and have therefore been the subject of many thorough studies on the ultrastructural and molecular properties of cells needed for regeneration. However, the technologies for obtaining transgenic echinoderms are limited and tracking cells involved in regeneration, and thus identifying the cellular sources and potencies has proven challenging. In this study we develop new transgenic tools for cell tracking in the regenerating bipinnaria larva of the sea star Patira minaita. We show that the larval serotonergic nervous system can regenerate following decapitation. Using a BAC-transgenesis approach with photoconvertible fluorescent proteins, we show that expression of the pan ectodermal marker, sox2, is induced in previously sox2 minus cells at the wound site, even when cell division is inhibited. sox2+ cells give rise to new sox4+ neural precursors that then proceed along an embryonic neurogenesis pathway to reform the anterior nervous systems. sox2+ cells contribute to only neural and ectoderm lineages, indicating that these progenitors maintain their normal, embryonic lineage restriction. This indicates that sea star larval regeneration uses a combination of existing lineage restricted stem cells, as well as respecification of cells into neural lineages, and at least partial reuse of developmental GRNs to regenerate their nervous system.


title: Snrpb, the gene mutated in CCMS, is haploinsufficient and required for proper morphogenesis and splicing

abstract: Heterozygous mutations in SNRPB, an essential core component of the five small ribonucleoprotein particles of the spliceosome, are responsible for Cerebrocostomandibular Syndrome (CCMS). We show that Snrpb heterozygous embryos arrest shortly after implantation. Additionally, heterozygous deletion of Snrpb in the developing brain and neural crest cells models craniofacial malformations found in CCMS, and results in death shortly after birth. RNAseq analysis of mutant heads prior to morphological defects revealed increased exon-skipping and intron-retention in association with increased 5 splice site strength. We found increased exon-skipping in negative regulators of the P53-pathway, increased levels of nuclear P53, P53-target genes. However, removing TrP53 in Snrpb heterozygous mutant neural crest cells did not completely rescue craniofacial development. We also found a small but significant increase in exon-skipping of several transcripts required for head and midface development, including Smad2 and Rere. Furthermore, mutant embryos exhibited ectopic or missing expression of Fgf8 and Shh, which are required to coordinate face and brain development. Thus, we propose that mis-splicing of transcripts that regulate P53-activity and craniofacial-specific genes both contribute to craniofacial malformations.


title: Snrpb, the gene mutated in CCMS, is haploinsufficient and required in neural crest cells for proper splicing of developmentally important genes essential for craniofacial morphogenesis


title: Snrpb, the CCMS gene, is required in neural crest cells for proper splicing of genes essential for craniofacial morphogenesis


title: A warmer growing season triggers earlier following spring phenology

abstract: Under global warming, advances in spring phenology due to the rising temperature have been widely reported. However, the physiological mechanisms underlying the warming-induced earlier spring phenology remain poorly understood. Here, using multiple long-term and large-scale phenological datasets between 1951 and 2018, we show that warmer temperatures during the previous growing season between May and September led to earlier spring phenology in the Northern Hemisphere. We also found that warming-induced increases in maximum photosynthetic rate in the previous year advanced spring phenology, with an average of 2.50 days {degrees}C-1. Furthermore, we found a significant decline in the advancing effect of warming during the previous growing season on spring phenology from cold to warm periods over the past decades. Our results suggest that the observed warming-induced earlier spring phenology may be driven by increased photosynthetic carbon assimilation in the previous season, while the slowdown in the advanced spring phenology arise likely from decreased carbon assimilation when warming exceeding the optimal temperatures for photosynthesis. Our study highlights the vital role of photosynthetic carbon assimilation during growing season in spring phenology under global warming.

category: ecology


title: Emigration Effects on Estimates of Age- and Sex-specific Survival of Small Mammals

abstract:
O_LIAge- and sex-specific survival estimates are crucial to understanding important life-history characteristics and variation in these estimates can be a key driver of population dynamics. When estimating survival using Cormack-Jolly-Seber (CJS) models and capture-recapture data, emigration is typically assumed to have a negligible effect on estimates such that apparent survival is indistinguishable from true survival. Consequently, especially for populations or age classes with high dispersal rates, apparent survival estimates are often biased low and temporal patterns in survival might be masked when site fidelity varies temporally.
C_LIO_LIWe used 9 years of annual mark-recapture data to estimate age-, sex-, and time-specific apparent survival of Humboldts flying squirrels (Glaucomys oregonensis) and Townsends chipmunks (Neotamias townsendii). For Humboldts flying squirrels, these estimates support a small body of research investigating potential variation of survival among age and sex classes, but age- and sex-specific survival has not been evaluated for Townsends chipmunks. We also quantified the effects of age- and sex-specific emigration on confounded estimates of apparent survival.
C_LIO_LIOur estimates of juvenile flying squirrel survival were high relative to other small mammal species and estimates for both species were variable among years. We found survival differed moderately among age and sex classes for Humboldts flying squirrels, but little among age and sex classes for Townsends chipmunks, and that the degree to which emigration confounded apparent survival estimates varied substantially among years. Without correcting for emigration, apparent survival estimates were lower and temporal variation was obscured, particularly for male Humboldts flying squirrels and female Townsends chipmunks.
C_LIO_LIOur results demonstrate that emigration can influence commonly used estimates of apparent survival. Unadjusted estimates confounded the interpretation of differences in survival between age and sex classes and masked potential temporal patterns in survival because the magnitude of adjustment varied among years. We conclude that apparent survival estimators are robust during some time periods; however, when emigration rates vary in time the effects of emigration should be carefully considered and accounted for, especially in comparative studies and those with policy and conservation implications.
C_LI


title: Genomic constraints to drought adaptation

abstract: Global shifts in climate and precipitation patterns are altering the diversity and structure of forests. The ability for species to adapt is especially difficult for long lived foundation species with unknown genetic and trait diversity. We harnessed genomic, physiological, and climate data to determine adaptation constraints. We used denovo assembly and 6.5 million genomic variants with drought related traits from 432 individuals sourced from across the complete range of the foundation tree species Corymbia calophylla. We found genomic variants determining traits predominantly in gene regulatory regions. The ability for populations to adapt was limited by within population genetic diversity associated with traits, and epistatic interactions within traits and pleiotropic interactions among traits. Nevertheless, we could accurately predict adaptive traits using genomic and climate data to guide reforestation. Our study indicated that some populations may contain variation sufficient for the species to adapt to the effects of drought, while other populations will need increased variability from those sources.

category: evolutionary biology


title: Ccr4-Not complex reduces transcription efficiency in heterochromatin

abstract: Heterochromatic silencing is thought to occur through a combination of transcriptional silencing and RNA degradation, but the relative contribution of each pathway is not known. In this study we analyzed RNA Polymerase II (RNA Pol II) occupancy and levels of nascent and steady-state RNA in different strains of fission yeast, in order to quantify the contribution of each pathway to heterochromatic silencing. We found that transcriptional silencing consists of two components, reduced RNA Pol II accessibility and, unexpectedly, reduced transcriptional efficiency. Heterochromatic loci showed lower transcriptional output compared to euchromatic loci, despite the presence of comparable amounts of RNA Pol II in both types of regions. We determined that the Ccr4-Not complex and H3K9 methylation are required for reduced transcriptional efficiency in heterochromatin and that a subset of heterochromatic RNA is degraded more rapidly than euchromatic RNA. Finally, we quantified the contribution of different chromatin modifiers, RNAi and RNA degradation to each silencing pathway. Our data show that several pathways contribute to heterochromatic silencing in a locus-specific manner and reveal transcriptional efficiency as a new mechanism of silencing.

category: genomics


title: Modeling bottlenecks, modularity, and context-dependency in behavioral control

abstract: In almost all animals, the transfer of information from the brain to the motor circuitry is facilitated by a relatively small number of neurons, leading to a constraint on the amount of information that can be transmitted. Our knowledge of how animals encode information through this pathway, and the consequences of this encoding, however, is limited. In this study, we use a simple feed-forward neural network to investigate the consequences of having such a bottleneck and identify aspects of the network architecture that enable robust information transfer. We are able to explain some recently observed properties of descending neurons - that they exhibit a modular pattern of connectivity and that their excitation leads to consistent alterations in behavior that are often dependent upon the prior behavioral state (context-dependency). Our model predicts that in the presence of an information bottleneck, such a modular structure is needed to increase the efficiency of the network and to make it more robust to perturbations. However, it does so at the cost of an increase in context-dependency. Despite its simplicity, our model is able to provide intuition for the trade-offs faced by the nervous system in the presence of an information processing constraint and makes predictions for future experiments.

category: animal behavior and cognition


title: A high-accuracy and high-efficiency digital volume correlation method to characterize in-vivo optic nerve head biomechanics from optical coherence tomography

abstract: In-vivo optic nerve head (ONH) biomechanics characterization is emerging as a promising way to study eye physiology and pathology. We propose a high-accuracy and high-efficiency digital volume correlation (DVC) method for the purpose of characterizing the in-vivo ONH deformation from volumes acquired by optical coherence tomography (OCT). Using a combination of synthetic tests and analysis of OCTs from monkey ONHs subjected to acute and chronically elevated intraocular pressure, we demonstrate that our proposed methodology overcomes several challenges for conventional DVC methods. First, it accounts for large ONH rigid body motion in the OCT volumes which could otherwise lead to analysis failure; second, sub-voxel-accuracy displacement can be guaranteed despite high noise and low image contrast of some OCT volumes; third, computational efficiency is greatly improved, such that the memory consumption of our method is substantially lower than with conventional methods; fourth, we introduce a parameter measuring displacements confidence. Test of image noise effects showed that the proposed DVC method had displacement errors smaller than 0.028 voxels with speckle noise and smaller than 0.037 voxels with Gaussian noise; The absolute (relative) strain errors in the three directions were lower than 0.0018 (4%) with speckle noise and than 0.0045 (8%) with Gaussian noise. Compared with conventional DVC methods, the proposed DVC method had substantially improved overall displacement and strain errors under large body motions (lower by up to 70%), with 75% lower computation times, while saving about 30% memory. The study thus demonstrates the potential of the proposed technique to investigate ONH biomechanics.

category: bioengineering


title: A divide and conquer metacell algorithm for scalable scRNA-seq analysis

abstract: Scaling scRNA-seq to profile millions of cells is increasingly feasible. Such data is crucial for the construction of high-resolution maps of transcriptional manifolds. But current analysis strategies, in particular dimensionality reduction and two-phase clustering, offers only limited scaling and sensitivity to define such manifolds. Here we introduce Metacell-2, a recursive divide and conquer algorithm allowing efficient decomposition of scRNA-seq datasets of any size into small and cohesive groups of cells denoted as metacells. We show the algorithm outperforms current solutions in time, memory and quality. Importantly, Metacell-2 also improves outlier cell detection and rare cell type identification, as we exemplify by analysis of human bone marrow cell atlas and mouse embryonic data. Metacell-2 is implemented over the scanpy framework for easy integration in any analysis pipeline.

category: bioinformatics


title: Identification and Targeting Putative G-Quadruplex Sequences in Candida glabrata: A Route to Virulence and Pathogenesis Control?

abstract: An increase in the number of Candida species that are resistant to antifungal medication and increases worldwide. Even individuals that are never exposed to antibiotics showing resistance to antifungal drugs. The increase in resistant candida species strains requires a search for novel targets for new antifungal agents. Preventing infection caused by candida species is a tremendous challenge in medicine. Although availability and use of antifungal drugs, disseminated candidiasis leads to a high mortality rate of about 40-60%, poor diagnosis, and improper disease management. Interest in G-quadruplexes as a therapeutic target has been increased in recent years, following the implication of this non-canonical G-quadruplex secondary structure in pathological diseases. However, G-quadruplex has been reported in many pathogens contributing to virulence and pathogenesis, including bacterial pathogens such as Staphylcoccus aureus and Enterococcus spp. Etc, viruses such as SARS-CoV2, HIV, HPV etc., and Fungi such as Candida species, Aspergillus fumigatus, etc. The present aim of the study is to identify and targeting G-Quadruplex forming sequences present in the Candida glabrata. PQSFinder (R Package) identified more than 5000 putative G-quadruplex forming sequences. Out of these, we have used PQS present in the SDH1 gene of Candida glabrata. It may be a key target to ameliorate the C.glabrata infection because it encodes a protein that plays a vital role in energy production in C.glabrata cells shown in the figure. The structure was built based on the already current structure that is telomeric G4 AGGG (TTAGGG)3 {PDB ID-4G0F} for identifying its stabilization by Gold carbine derivatives. Molecular docking and ADMET analysis show that Compound A has the highest binding affinity and has the best ADMET properties among the two compounds. The present study represents PQS in the SDH1 gene could be a novel antifungal target.


title: BioPhi: A platform for antibody design, humanization and humanness evaluation based on natural antibody repertoires and deep learning

abstract:
Despite recent advances in transgenic animal models and display technologies, humanization of mouse sequences remains the primary route for therapeutic antibody development. Traditionally, humanization is manual, laborious, and requires expert knowledge. Although automation efforts are advancing, existing methods are either demonstrated on a small scale or are entirely proprietary. To predict the immunogenicity risk, the human-likeness of sequences can be evaluated using existing humanness scores, but these lack diversity, granularity or interpretability. Meanwhile, immune repertoire sequencing has generated rich antibody libraries such as the Observed Antibody Space (OAS) that offer augmented diversity not yet exploited for antibody engineering. Here we present BioPhi, an open-source platform featuring novel methods for humanization (Sapiens) and humanness evaluation (OASis). Sapiens is a deep learning humanization method trained on the OAS database using language modeling. Based on an in silico humanization benchmark of 177 antibodies, Sapiens produced sequences at scale while achieving results comparable to that of human experts. OASis is a granular, interpretable and diverse humanness score based on 9-mer peptide search in the OAS. OASis separated human and non-human sequences with high accuracy, and correlated with clinical immunogenicity. Together, BioPhi offers an antibody design interface with automated methods that capture the richness of natural antibody repertoires to produce therapeutics with desired properties and accelerate antibody discovery campaigns.

BioPhi is accessible at https://biophi.dichlab.org and https://github.com/Merck/BioPhi.


title: Topological morphogenesis of neuroepithelial organoids

abstract: Animal organs exhibit complex topologies involving cavities and tubular networks, which underlie their form and function. However, how topology emerges during organ morphogenesis remains elusive. Here, we combine tissue reconstitution and quantitative microscopy to show that trans and cis epithelial fusion govern tissue topology and shape. These two modes of topological transitions can be regulated in neuroepithelial organoids, leading to divergent topologies. The morphological space can be captured by a single control parameter which is analogous to the reduced Gaussian rigidity of an epithelial surface. Finally, we identify a pharmacologically accessible pathway that regulates the frequency of trans and cis fusion, and demonstrate the control of organoid topology and shape. The physical principles uncovered here provide fundamental insights into the self-organization of complex tissues.

category: biophysics


title: Stabilizing proteins through saturation suppressor mutagenesis.

abstract: While there have been recent, transformative advances in the area of protein structure prediction, prediction of point mutations that improve protein stability remains challenging. It is possible to construct and screen large mutant libraries for improved activity or ligand binding, however reliable screens for mutants that improve protein stability do not exist, especially for proteins that are well folded and relatively stable. We demonstrate that incorporation of a single, specific destabilizing, (parent inactive) mutation into each member of a deep mutational scanning library followed by screening for suppressors, allows for robust and accurate identification of stabilizing mutations. When coupled to FACS sorting of a yeast surface display library of the bacterial toxin CcdB, followed by deep sequencing of sorted populations, multiple stabilizing mutations could be identified after a single round of sorting. Multiple libraries with different parent inactive mutations could be pooled and simultaneously screened to further enhance the accuracy of identification of stabilizing mutations. Individual stabilizing mutations could be combined to result in a multimutant with increase in thermal melting temperature of about 20 {degrees}C and enhanced tolerance to high temperature exposure. The method employs small library sizes and can be readily extended to other display and screening formats to rapidly isolate stabilized protein mutants.


title: Improving the analysis of biological ensembles through extended similarity measures

abstract: We present new algorithms to classify structural ensembles of macromolecules, based on the recently proposed extended similarity measures. Molecular Dynamics provides a wealth of structural information on systems of biologically interest. As computer power increases we capture larger ensembles and larger conformational transitions between states. Typically, structural clustering provides the statistical mechanics treatment of the system to identify relevant biological states. The key advantage of our approach is that the newly introduced extended similiarity indices reduce the computational complexity of assessing the similarity of a set of structures from O(N2) to O(N). Here we take advantage of this favorable cost to develop several highly efficient techniques, including a linear-scaling algorithm to determine the medoid of a set (which we effectively use to select the most representative structure of a cluster). Moreover, we use our extended similarity indices as a linkage criterion in a novel hierarchical agglomerative clustering algorithm. We apply these new metrics to analyze the ensembles of several systems of biological interest such as folding and binding of macromolecules (peptide,protein,DNA -protein). In particular, we design a new workflow that is capable of identifying the most important conformations contributing to the protein folding process. We show excellent performance in the resulting clusters (surpassing traditional linkage criteria), along with faster performance and an efficient cost-function to identify when to merge clusters.


title: Direct observation of Hsp90-induced compaction in a protein chain

abstract: The chaperone Hsp90 is well known to undergo important conformational changes, which depend on nucleotide, co-chaperones, substrate interactions and post-translational modifications. Conversely, how the conformations of its unstable and disordered substrates are affected by Hsp90 is difficult to address experimentally, yet central to its function. Here, using optical tweezers and luciferase and glucocorticoid receptor substrates, we find that Hsp90 promotes local contractions in unfolded chains that drive their global compaction down to dimensions of folded states. This compaction has a gradual nature while showing small steps, is stimulated by ATP, and performs mechanical work against counteracting forces that expand the chain dimensions. The Hsp90 interactions suppress the formation of larger-scale folded, misfolded and aggregated structures. The observations support a model in which Hsp90 alters client conformations directly by promoting local intra-chain interactions while suppressing distant ones. We conjecture that chain compaction may be central to how Hsp90 protects unstable kinases and receptor clients, regulates their activity, and how Hsp90 cooperates with Hsp70.


title: An Isogenic Cell Line Panel for Sequence-based Screening of Targeted Anti-cancer Drugs

abstract: We describe the creation and characterization of an isogenic cell line panel representing common cancer pathways, with multiple features optimized for high-throughput screening. More than 1,800 cell lines from three normal human cells were generated using CRISPR-technologies. Surprisingly, we discovered most of these lines did not result in complete gene inactivation, despite integration of sgRNA at the desired genomic site. However, a subset of the lines harbored true, biallelic disruptions of the targeted tumor suppressor gene, yielding a final panel of 100 well-characterize lines covering 19 pathways frequently subject to loss of function in cancers. This panel included genetic markers optimized for sequence-based ratiometric assays for drug-based screening assays. To illustrate the potential utility of this panel, we developed a multiplexed high-throughput screen that identified Wee1 inhibitor MK-1775 as a selective growth inhibitor of cells with inactivation of TP53. These cell lines and screening approach should prove useful for researchers studying a variety of cellular and biochemical phenomena.

category: cancer biology


title: Treatment with a polyherbal extract improves fat metabolism, attenuates hepatic stellate cell activation and fibrogenesis.

abstract: Non-alcoholic steatohepatitis (NASH) involves dysregulations in denovo lipogenesis, fatty acid oxidation, and fibrogenesis. Targeting these pathways holds promise for the treatment of liver disorders. Here we test the extract of a polyherbal formulation (namely Liv.52), which is approved by the Government of Indias Drug Regulatory Authority - AYUSH. The current study evaluates the effect of Liv.52 on denovo lipogenesis, fatty acid oxidation, and fibrogenesis. Both in vivo and in vitro model systems were employed to evaluate the efficacy of this polyherbal formulation. Male Wistar rats were dosed with Liv.52 for 2 weeks (250mg/k.g) and expression levels of the genes involved in de novo lipogenesis and fatty acid oxidation pathways were analysed by quantitative real time PCR. Liv.52 treatment resulted in increased hepatic fatty acid oxidation and decreased de novo lipogenesis in these rats. It also reduced hepatic stellate cell activation in CCL4 treated Wistar rats as evidenced by histological evaluation. For in vitro experiments, HepG2 cells were cultured under lipotoxic conditions (using 200M palmitic acid) and the conditioned media from these cells were used for inducing activation and fibrogenesis in human hepatic stellate cells (HHSteC). Treatment with lipotoxic conditioned media resulted in activation of hepatic stellate cells and fibrogenesis, as evidenced by increased expression of -smooth muscle actin (-SMA), and desmin (markers of stellate cell activation) and increased levels of collagen and lumican (markers of fibrogenesis). Treatment with Liv.52 reversed the up-regulation of -SMA, collagen and lumican levels in HHSteC cells. These results indicate that Liv.52 exerts its hepatoprotective effect by improving fatty acid metabolism and fibrogenesis.

category: pharmacology and toxicology


title: A method for simultaneously monitoring phloem and xylem reconnections in grafted watermelon seedlings

abstract: Grafting is an effective way to increase watermelon tolerance to biotic and abiotic stresses. However, the survival of grafted seedlings largely depends on successful graft formation. Therefore, understanding the graft formation process, particularly the vascular reconnection process is of critical importance. This study found that lignin in watermelon stem shows strong auto-fluorescence under blue-light excitation which makes blue-light excited fluorescent tracers (FTs) such as 5(6)-carboxy fluorescein diacetate (CFDA) become unsuitable for assaying vascular connectivity in watermelon. In contrast, UV-light excited esculin and red-light excited acid fuchsin were proved to be efficient FTs for monitoring the phloem and xylem connectivity, respectively, in self-grafted watermelon. Furthermore, a combined application of esculin to the scion cotyledon and acid fuchsin to the rootstock root enabled simultaneous monitoring of the phloem and xylem connectivity in individual self-grafted watermelon seedlings. In addition, this method is also applicable in investigating the phloem and xylem reconnections in self-grafted melon and cucumber, and heterograft of watermelon, melon and cucumber onto pumpkin rootstock. Based on this established method, we found that phloem and xylem reconnections are not timely separated in self-grafted watermelon. Furthermore, low temperature and removal of the rootstock cotyledons both delayed the vascular reconnection process in watermelon. In conclusion, this new method provides a convenient, accurate and rapid way to analyze the vascular connectivity not only in watermelon, but also in other cucurbit crops.

category: plant biology


title: Temporal change in chromatin accessibility predicts regulators of nodulation in Medicago truncatula

abstract:
Rhizobia can establish symbiotic associations with legumes to provide plants with nitrogen needed in agricultural systems. Symbiosis triggers extensive genome and transcriptome remodeling in the plant, yet the extent of chromatin changes and impact on gene expression is unknown. We profiled the temporal chromatin accessibility (ATAC-seq) and transcriptome (RNA-seq) dynamics of M. truncatula roots treated with rhizobia lipo-chitooligosaccharides. Using a novel approach, Dynamic Regulatory Module Networks, we predicted gene expression as a function of chromatin accessibility and accessible cis-regulatory elements. This approach identified the cis-regulatory elements and associated transcription factors that most significantly contribute to transcriptomic changes triggered by lipo-chitooligosaccharides. Regulators involved in auxin (IAA4-5,SHY2), ethylene (EIN3, ERF1) and abscisic acid (ABI5) hormone response, as well as histone and DNA methylation (IBM1), emerged among those most predictive of transcriptome dynamics. RNAi-based knockdown of EIN3 and ERF1 reduced nodule number in M. truncatula validating the role of these predicted regulators in symbiosis between legumes and rhizobia.

Significance StatementLegumes can fix nitrogen through symbiosis with rhizobia in root nodules, a critical mutualistic relationship for crop productivity and agricultural sustainability. Introducing this symbiotic relationship into non-legume crops is of great interest, but limited knowledge of host genome modifications induced by rhizobia has hampered such efforts. We applied time-course analysis of chromatin accessibility and gene expression of M. truncatula roots treated with rhizobia lipochitooligosaccharides. We show that extensive remodeling of genome accessibility drives a large component of the temporal transcriptome dynamics. By predicting gene expression as a function of accessibility of regulatory features, we identified known and novel regulators that are associated with early nodule development, which may be critical for its engineering into crops.


title: Temporal dynamics of neural responses in human visual cortex

abstract:
Neural responses to visual stimuli exhibit complex temporal dynamics, including sub-additive temporal summation, response reduction with repeated or sustained stimuli (adaptation), and slower dynamics at low contrast. These phenomena are often studied independently. Here, we demonstrate these phenomena within the same experiment and model the underlying neural computations with a single computational model. We extracted time-varying responses from electrocorticographic (ECoG) recordings from patients presented with stimuli that varied in contrast, duration, and inter-stimulus interval (ISI). Aggregating data across patients yielded 88 electrodes with robust visual responses, covering both earlier (V1-V3) and higher-order (V3a/b, LO, TO, IPS) retinotopic maps. In all regions, the temporal dynamics of neural responses exhibit several non-linear features: peak response amplitude saturates with high contrast and longer stimulus durations; the response to a second stimulus is suppressed for short ISIs and recovers for longer ISIs; response latency decreases with increasing contrast. These features are accurately captured by a computational model comprised of a small set of canonical neuronal operations: linear filtering, rectification, exponentiation, and a delayed divisive normalization. We find that an increased normalization term captures both contrast- and adaptation-related response reductions, suggesting potentially shared underlying mechanisms. We additionally demonstrate both changes and invariance in temporal response dynamics between earlier and higher-order visual areas. Together, our results reveal the presence of a wide range of temporal and contrast-dependent neuronal dynamics in the human visual cortex, and demonstrate that a simple model captures these dynamics at millisecond resolution.

Significance StatementSensory inputs and neural responses change continuously over time. It is especially challenging to understand a system that has both dynamic inputs and outputs. Here we use a computational modeling approach that specifies computations to convert a time-varying input stimulus to a neural response time course, and use this to predict neural activity measured in the human visual cortex. We show that this computational model predicts a wide variety of complex neural response shapes that we induced experimentally by manipulating the duration, repetition and contrast of visual stimuli. By comparing data and model predictions, we uncover systematic properties of temporal dynamics of neural signals, allowing us to better understand how the brain processes dynamic sensory information.


title: Recombinant NAGLU-IGF2 prevents physical and neurological disease and improves survival in Sanfilippo B syndrome

abstract: Recombinant human alpha-N-acetylglucosaminidase-insulin-like growth factor-2 (rhNAGLU-IGF2) is an investigational enzyme replacement therapy for Sanfilippo B, a lysosomal storage disease. Because recombinant human NAGLU (rhNAGLU) is poorly mannose 6-phosphorylated, we generated a fusion protein of NAGLU with IGF2 to permit its binding to the cation-independent mannose 6-phosphate receptor. We previously administered rhNAGLU-IGF2 intracerebroventricularly to Sanfilippo B mice, and demonstrated therapeutic restoration of NAGLU, normalization of lysosomal storage, and improvement in markers of neurodegeneration and inflammation. Here, we studied intracerebroventricular rhNAGLU-IGF2 delivery in both murine and canine Sanfilippo B to determine potential effects on their behavioral phenotypes and survival. Treated mice showed improvement in disease markers such as heparan sulfate glycosaminoglycans, beta-hexosaminidase, microglial activation, and lysosomal-associated membrane protein-1. Sanfilippo B mice treated with rhNAGLU-IGF2 displayed partial normalization of their stretch attend postures, a defined fear pose in mice (p<0.001). We found a more normal dark/light activity pattern in Sanfilippo B mice treated with rhNAGLU-IGF2 compared to vehicle-treated Sanfilippo B mice (p=0.025). We also found a 61% increase in survival in Sanfilippo B mice treated with rhNAGLU-IGF2 (mean 53w, median 48w) compared to vehicle-treated Sanfilippo B mice (mean 33w, median 37w; p<0.001). In canine Sanfilippo B, we found that rhNAGLU-IGF2 administered into cerebrospinal fluid normalized HS and beta-hexosaminidase activity in gray and white matter brain regions. Proteomic analysis of cerebral cortex showed restoration of protein expression levels in pathways relevant to cognitive function, synapse, and the lysosome. These data suggest that treatment with rhNAGLU-IGF2 may improve the phenotype of Sanfilippo B disease.


title: The looping lullaby: closed-loop neurostimulation decreases sleepers' sensitivity to environmental noise

abstract:
Study ObjectiveSleep is critical for physical and mental health. However, sleep disruption due to noise is a growing problem, causing long-lasting distress and fragilizing entire populations mentally and physically. Here for the first time, we tested an innovative and non-invasive potential countermeasure for sleep disruptions due to noise.

MethodsWe developed a new, modeling-based, closed-loop acoustic neurostimulation procedure (CLNS) to precisely phase-lock stimuli to slow oscillations (SO). We used CLNS to align, soft sound pulses to the start of the SO positive deflection to boost SO and sleep spindles during non-rapid eye movement (NREM) sleep. Participants underwent three overnight EEG recordings. The first night served to determine each participants individual noise arousal threshold. The remaining two nights occurred in counterbalanced order: in the "Disturbing night", loud, real-life noises were repeatedly presented; in the "Intervention night", similar loud noises were played while using the CLNS to boost SO. All experimental manipulations were performed in the first three hours of sleep; participants slept undisturbed for the rest of the night.

ResultsIn contrast to the Disturbing night, the probability of arousals caused by noise was significantly decreased in the Intervention night. Moreover, the CLNS intervention increased NREM duration and sleep spindle power across the night.

ConclusionsThese results show that our CLNS procedure can effectively protect sleep from disruptions caused by noise. Remarkably, even in the presence of loud environmental noise, CLNS soft and precisely timed sound pulses played a beneficial role in protecting sleep continuity. This represents the first successful attempt at using CLNS in a noisy environment.

Statement of SignificanceExposure to noise during sleep impairs sleep quality, leading to impaired biological functioning, decreased day-time performance and higher occurrence of accidents. Slow wave sleep is hallmarked by slow oscillations (SO) and sleep spindles, both of which are markers of sleep stability. In previous experiments, acoustic stimulation has shown potential in enhancing SO and sleep spindles. Whether these manipulations also work to stabilize sleep against the disrupting effects of environmental noise remains unclear. Here we use a new closed-loop approach to precisely align subtle acoustic stimuli with SO phase. We show that this method effectively stabilizes sleep in the presence of noise. Our results bear crucial relevance for improving sleep in the general public and in many high-risk professions.


title: Neuronal responses to focused ultrasound are gated by pre-stimulation brain rhythms

abstract:
BackgroundOwing to its high spatial resolution and penetration depth, transcranial focused ultrasound stimulation (tFUS) is one of the most promising approaches to non-invasive neuromodulation. Identifying the impact of the stimulation waveform and endogenous neural activity on neuromodulation outcome is critical to harnessing the potential of tFUS.

ObjectiveHere we tested a new form of tFUS where the amplitude of the ultrasonic waveform is modulated at a rate much slower than the operating frequency. Moreover, we sought to identify the relationship between pre-stimulation neural activity and the neuronal response to tFUS.

MethodsWe applied three minutes of amplitude modulated (AM) tFUS at 40 Hz to the rat hippocampus while recording local field potentials (LFP) and multi-unit activity (MUA)from the sonicated region. To assess the role of AM, we also tested continuous-wave (CW) stimulation.

ResultsAM tFUS reduced firing rate during and immediately after stimulation. On the other hand, CW tFUS produced an acute firing rate increase that was abolished after sonication. For both waveforms, firing rate changes were stronger in units exhibiting high baseline LFP power, particularly in the gamma band (30-250 Hz). The neuromodulatory effect was also influenced by the prevalence of sharp wave ripples (SWR) during the pre-stimulation period, with firing rates modulated by up to 33% at units showing frequent baseline SWR.

ConclusionOur findings suggest that AM and CW tFUS produce qualitatively different neuronal outcomes, and that baseline rhythms may effectively "gate" the response to tFUS.


title: Visual association cortex links cues with conjunctions of reward and locomotor contexts

abstract: Postrhinal cortex (POR) and neighboring lateral visual association areas are necessary for identifying objects and interpreting them in specific contexts, but how POR neurons encode the same object across contexts remains unclear. Here, we imaged excitatory neurons in mouse POR across tens of days throughout initial cue-reward learning and reversal learning. As such, neurons were tracked across sessions/trials where the same cue was rewarded or unrewarded, during both locomotor and stationary contexts. Surprisingly, a large class of POR neurons were minimally cue-driven prior to learning. After learning, distinct clusters within this class responded selectively to a given cue when presented in a specific conjunction of reward and locomotion contexts. In addition, another class involved clusters of neurons whose cue responses were more transient, insensitive to reward learning, and adapted over thousands of presentations. These two classes of POR neurons may support context-dependent interpretation and context-independent identification of sensory cues.


title: Columnar processing of border ownership in primate visual cortex

abstract:
To understand a visual scene, the brain segregates figures from background by assigning borders to foreground objects. Neurons in primate visual cortex encode which object owns a border (border ownership), but the underlying circuitry is not understood. Here we used multielectrode probes to record from border ownership selective units in different layers in macaque visual area V4 to study the laminar organization and timing of border ownership selectivity.

We find that border ownership selectivity occurs first in deep layer units, in contrast to spike latency for small stimuli in the classical receptive field. Units on the same penetration typically share the preferred side of border ownership, also across layers, similar to orientation preference. Units are often border ownership selective for a range of border orientations, where the preferred sides of border ownership are systematically organized in visual space.

Together our data reveal a columnar organization of border ownership in V4 where the earliest border ownership signals are not simply inherited from upstream areas, but computed by neurons in deep layers, and may thus be part of signals fed back to upstream cortical areas or the oculomotor system early after stimulus onset. The finding that preferred border ownership is clustered and can cover a wide range of spatially contiguous locations, suggests that the asymmetric context integrated by these neurons is provided in a systematically clustered manner, possibly through corticocortical feedback and horizontal connections.


title: Binding of Monomeric and Polymeric Alzheimers Aβ peptides to Exosomes

abstract: Exosomes are secreted by every cell in our body under both physiological and pathological conditions. They travel in the blood, CSF, and all studied biofluids. Their biological roles have been reported to include delivery of important physiological cargo between organs and cells, clearance of toxic proteins; maintenance of cellular stasis, and the propagation of disease pathology. In the case of Alzheimers disease (AD) exosomes have been shown to carry pathological proteins such as amyloid, yet the specificity of this association of amyloid and exosomes is unclear. To address this deficiency, we utilized Isothermal Titration Calorimetry (ITC) to measure the binding of amyloid to exosomes. Here we report that A{beta}40 and A{beta}42 bind to exosomes in a saturable and endothermic manner, a phenomenon not observed with the scrambled versions of either peptide. This points to this interaction being more specific than previously understood, and to amyloid associated with exosomes as an important pool of this peptide in the plasma.


title: Task space exploration improves adaptation after incompatible virtual surgeries

abstract:
Humans have a remarkable capacity to learn new motor skills, a process that requires novel muscle activity patterns. Muscle synergies may simplify the generation of muscle patterns through the selection of a small number of synergy combinations. Learning new motor skills may then be achieved by acquiring novel muscle synergies. In a previous study, we used myoelectric control to construct virtual surgeries that altered the mapping from muscle activity to cursor movements. After compatible virtual surgeries, which could be compensated by recombining subject-specific muscle synergies, participants adapted quickly. In contrast, after incompatible virtual surgeries, which could not be compensated by recombining existing synergies, participants explored new muscle patterns, but failed to adapt. Here, we tested whether task space exploration can promote learning of novel muscle synergies, required to overcome an incompatible surgery. Participants performed the same reaching task as in our previous study, but with more time to complete each trial, thus allowing for exploration. We found an improvement in trial success after incompatible virtual surgeries. Remarkably, improvements in movement direction accuracy after incompatible surgeries occurred faster for corrective movements than for the initial movement, suggesting that learning of new synergies is more effective when used for feedback control. Moreover, reaction time was significantly higher after incompatible than after compatible virtual surgeries, suggesting an increased use of an explicit adaptive strategy to overcome incompatible surgeries. Taken together, these results indicate that exploration is important for skill learning and suggest that human participants, with sufficient time can learn new muscle synergies.

NEW & NOTEWORTHYMotor skill learning requires the acquisition of novel muscle patterns, a slow adaptive process. Here we show that learning to control a cursor after an incompatible virtual surgery, a complex skill requiring new muscle synergies, is possible when enough time for task space exploration is provided. Our results suggest that learning new synergies is related to the exceptional human capacity to acquire a wide variety of novel motor skills with practice.


title: Gpr19 is a circadian clock-controlled orphan GPCR with a role in modulating free-running period and light resetting capacity of the circadian clock

abstract:
Background and PurposeGpr19 encodes an evolutionarily conserved orphan G-protein-coupled receptor (GPCR) with no established physiological function in vivo. The purpose of this study was to determine the role of Gpr19 in the circadian clock system.

Experimental ApproachWe examined whether and how the master circadian clock neurons in the suprachiasmatic nucleus (SCN) express Gpr19. By analysing Gpr19-deficient (Gpr19-/-) mice, we asked whether Gpr19 has a role in modulating free-running period and light resetting capacity of the circadian clock.

Key ResultsCompared with the known common core clock genes, Gpr19 was identified to show several distinct yet limited features related to the circadian clock. Gpr19 mRNA was mainly expressed in the middle-to-dorsal region of the SCN. A conserved cAMP-responsive element within the Gpr19 promoter drove the circadian expression of Gpr19. Gpr19-/- mice exhibited a prolonged circadian period and a delayed initiation of daily locomotor activity in a 12-h light/12-h dark cycle. Gpr19 deficiency caused the downregulation of several genes that normally peak during the night, including Bmal1 and Gpr176. Gpr19-/- mice had a reduced capacity for phase shift to early subjective night light. The defect was only observed for phase-delay, but not phase-advance, and accompanied by reduced response of c-Fos expression in the dorsal region of the SCN, while apparently normal in the ventral part of the SCN, in Gpr19-/- mice.

Conclusion and ImplicationsGpr19 is an SCN-enriched orphan GPCR with a distinct role in circadian regulation and thus may be a potential target for alleviating circadian clock disorders.

O_LIWhat is already known: Gpr19 is an evolutionarily conserved class-A orphan receptor with no established physiological role in vivo.
The SCN is a light-entrainable master circadian pacemaker governing daily rhythms of behaviour and physiology.

C_LIO_LIWhat this study adds: Gpr19 is an SCN-enriched orphan GPCR whose levels fluctuate in a circadian fashion.
Gpr19 is a functional clock modulator involved in period determination and phase resetting.

C_LIO_LIClinical significance: Targeting the orphan receptor Gpr19 may provide a therapeutic approach for alleviating circadian clock disorders.

C_LI


title: Augmenting Neutralization breadth against Diverse HIV-1 by increasing the Ab-Ag interface on V2

abstract:
Understanding maturation pathways of broadly neutralizing antibodies (bnAbs) against HIV-1 in non-human primates can be highly informative for HIV-1 vaccine development. We now obtained a lineage of J038 from Chinese rhesus macaques after 7-years of SHIV infection. J038 has short complementary determining loops and neutralizes 54% of global circulating HIV-1 strains. Its binding induces a unique "up" conformation for one of the V2 loops in the trimeric envelope glycoprotein (Env) and is heavily dependent on glycan, which provides nearly half of the binding surface. The unmutated common ancestor of the J038 lineage antibodies binds monomeric gp120 and neutralizes the autologous virus. Continuous maturation enhances neutralization potency and breadth of J038 lineage antibodies via expanding antibody-Env contact areas surrounding the core region contacted by germline-encoded residues. Developmental details and recognition features of J038 lineage antibodies revealed here provide a new pathway for maturation elicitation of V2-targeting bnAbs.

Highlights* Long-term infected NHPs develop antibodies neutralizing up to 54% of HIV-1 strains
* Antibody J038 binds one V2 loop on HIV-1 Env trimer in a unique "up" position
* UCA of the J038 lineage effectively neutralizes the autologous virus
* J038 lineage antibodies mature through gradually increased contact to glycans

category: immunology


title: Metataxonomic review to elucidate the role of the microbiome in Celiac disease across the gastrointestinal tract

abstract: BackgroundDysbiosis of the microbiome has been related to the Celiac disease (CeD) progress, an autoimmune disease characterised by gluten intolerance developed in genetically susceptible individuals under certain environmental factors. The microbiome contributes to CeD pathophysiology modulating the immune response by the action of short-chain fatty acids (SCFA), affecting gut barrier integrity allowing the entrance of gluten derived proteins, and degrading immunogenic peptides of gluten through endoprolyl peptidase enzymes. We reviewed state of the art in taxonomic composition for CeD and compiled the larger dataset of 16S prokaryotic ribosomal RNA (rRNA) gene high-throughput sequencing for consensus profiling. We present for the first time an integrative analysis of metataxonomic data from CeD patients, including samples from different body sites (saliva, pharynx, duodenum, and stool). We found the presence of coordinated changes through the gastrointestinal tract characterised by an increase in Actinobacteria species in the upper tract (pharynx and duodenum), and an increase in Proteobacteria in the lower tract (duodenum and stool), as well as site- specific changes evidencing a dysbiosis in CeD patients microbiota. Moreover, we described the effect of adherence to a gluten-free diet (GFD) evidenced by an increase in beneficial bacteria and a decrease in some Betaproteobacteriales but not fully restoring CeD-related dysbiosis. We illustrate that the gut microbiota acts as an enhancer of immune response in CeD through the production of lipopolysaccharides and other bacterial components that activate the immune response and by decrease SCFA producers bacteria. Furthermore, microbial changes observed through the gastrointestinal tract of CeD patients may help manage the disease and follow-up GFD treatment.

category: microbiology


title: Physiologically relevant miRNAs in mammalian oocytes are rare and highly abundant.

abstract: miRNAs, ~22nt small RNAs associated with Argonaute (AGO) proteins, are important negative regulators of gene expression in mammalian cells. However, mammalian maternal miRNAs show negligible repressive activity and the miRNA pathway is dispensable for oocytes and maternal-to-zygotic transition. The stoichiometric hypothesis proposed that this is caused by dilution of maternal miRNAs during oocyte growth. As the dilution affects miRNAs but not mRNAs, it creates unfavorable miRNA:mRNA stoichiometry for efficient repression of cognate mRNAs. Here we report that porcine ssc-miR-205 and bovine bta-miR-10b are exceptional miRNAs, which resist the diluting effect of oocyte growth and can efficiently suppress gene expression. Additional analysis of ssc-miR-205 shows that it has higher stability, reduces expression of endogenous targets, and contributes to porcine oocyte-to-embryo transition. Consistent with the stoichiometric hypothesis, our results show that the endogenous miRNA pathway in mammalian oocytes is intact and that maternal miRNAs can efficiently suppress gene expression when a favorable miRNA:mRNA stoichiometry is established.

category: molecular biology


title: ATR facilitates the degradation of Api5 through the ubiquitin-proteasome pathway via FBXW2 to regulate apoptosis upon DNA damage

abstract: Apoptosis inhibitor 5 (Api5) is an inhibitor of apoptosis, which is found to be upregulated in several cancers and promotes invasion as well as metastasis. Over-expression of Api5 is positively co-related with poor survival of cancers and inhibition of DNA damage induced apoptosis in cancerous cells. Acetylation at lysine 251 (K251) on Api5 facilitates the stability of the protein and thus functionally provides resistance to cancer cells against chemotherapeutic or anti-cancerous agents. However, the regulation of Api5 upon DNA damage is not yet known. In this study, we demonstrate that Api5 undergoes degradation following DNA damage via the ubiquitin-proteasome system. Upon DNA damage, ATR was observed to phosphorylate Api5 at serine 138 which led to the cytoplasmic localisation of Api5. The E3-ubiquitin ligase, SCF-FBXW2 ubiquitinates Api5 leading to its proteasomal degradation.


title: Temporally auto-correlated predator attacks structure ecological communities

abstract: For species primarily regulated by a common predator, the P* rule of Holt and Lawton [1993] predicts that the prey species that supports the highest mean predator density (P*) excludes the other prey species. This prediction is re-examined in the presence of temporal fluctuations in the vital rates of the interacting species including predator attack rates. When the fluctuations in predator attack rates are temporally uncorrelated, the P* rule still holds even when the other vital rates are temporally auto-correlated. However, when temporal auto-correlations in attack rates are positive but not too strong, the prey species can coexist due to the emergence of a positive covariance between predator density and prey vulnerability. This coexistence mechanism is similar to the storage effect for species regulated by a common resource. Negative or strongly positive auto-correlations in attack rates generate a negative covariance between predator density and prey vulnerability and a stochastic priority effect can emerge: with non-zero probability either prey species is excluded. These results highlight how temporally auto-correlated species interaction rates impact the structure and dynamics of ecological communities.


title: Representing the sediment retention ecosystem service within an ecosystem account of the Mitchell catchment

abstract: If soil resources and the benefits derived from water quality are to be maintained, the on- and off-site effects of soil erosion must be adequately represented and communicated to decision makers so that appropriate management responses can be identified. The System of Environmental-Economic Accounting - Ecosystem Accounting (SEEA-EA), is one approach to quantify both the contributions that ecosystems make to the economy, and the impacts of economic activity on ecosystems. However, due to the difficulty of obtaining empirical data on ecosystem service flows, in many cases such quantification is informed by ecosystem service models. Previous research on the Mitchell catchment, Queensland Australia allowed us to explore the implications of using a model of hillslope erosion and sediment delivery in isolation (as represented in one of the most frequently used ecosystem service models - InVEST) by comparing such estimates against multiple lines of local empirical data, and a more comprehensive representation of locally important erosion and deposition processes through a sediment budget model. Estimates of the magnitude of hillslope erosion modelled using an approach similar to InVEST and the calibrated sediment budget differed by an order of magnitude. If an uncalibrated InVEST-type model was used to inform the relative distribution of erosion magnitude, findings suggest the incorrect erosion process would be identified as the dominant contributor to suspended sediment loads. However, the sediment budget model could only be calibrated using data on sediment sources and sinks that had been collected in the catchment through a sustained research effort. A comparable level of research investment may not be available to inform ecosystem service assessments elsewhere. The results summarised here for the Mitchell catchment demonstrate that practitioners must exercise caution when using estimates of the sediment retention ecosystem service flow which have not been calibrated and validated against locally collected empirical data.


title: Challenges in modelling the sediment retention ecosystem service for an ecosystem account, with examples from the Mitchell catchment


title: Enhanced soil quality after forest conversion to vegetable cropland and tea plantation has contrasting effects on soil microbial structure and functions

abstract: Land-use changes could potentially exert a strong influence on soil quality and soil microbial communities. Moreover, microbial taxa are also important drivers of soil ecological functions. However, the linkage between soil quality and soil microbial communities is in need of deeper understanding. In this study, we examined the effects of soil quality on microbial community structure and functions after forest conversion to vegetable cropland and tea plantations. Soil quality index was significantly increased after natural forest conversion to vegetable cropland and tea plantations. Soil bacterial beta diversity significantly correlated to soil quality, but the sensitivity of individual microbial groups varied in response to changes in soil quality. Higher soil quality promoted bacterial diversity in vegetable cropland but decreased it in tea plantations, which implied soil quality was a structural factor in bacterial community composition but had contrasting effects for croplands versus plantations. Agricultural management played a negative role in maintaining microbial interactions, as identified by the network analysis, and furthermore the analysis revealed key functions of the microbial communities. After land-use change, the abundance (e.g., level, intensity) of microbial N-cycling function increased in tea plantations but decreased in vegetable cropland. The abundance of C-cycling function featured an opposite trend. Higher level of N-fixation in tea plantations but the higher abundance of N-oxidation in vegetable cropland was demonstrated. Higher abundance of ammonia-oxidizing bacteria and ammonia-oxidizing archaea as identified by qPCR in vegetable cropland corroborated the FAPROTAX function prediction. Therefore, the key taxa of soil microbial communities and microbial functions were largely dependent on changes in soil quality and determined responses to specific agricultural management.


title: Mitonuclear co-introgression opposes genetic differentiation between phenotypically divergent songbirds

abstract: Comparisons of genomic variation among closely related species often show more differentiation in mitochondrial DNA (mtDNA) and sex chromosomes than in autosomes, a pattern expected due to the differing effective population sizes and evolutionary dynamics of these genomic components. Yet, introgression can cause species pairs to deviate dramatically from general differentiation trends. The yellowhammer (Emberiza citrinella) and the pine bunting (E. leucocephalos) are hybridizing avian sister species that differ greatly in appearance and moderately in nuclear DNA, but that show no mtDNA differentiation. This mitonuclear discordance is best explained by adaptive mtDNA introgression--a process that can select for co-introgression at nuclear genes with mitochondrial functions (mitonuclear genes). To better understand the extent of mitonuclear discordance and characterize nuclear differentiation patterns in this system, we investigated genome-wide differentiation between allopatric yellowhammers and pine buntings and compared it to what was seen previously in mtDNA. We found significant nuclear differentiation that was highly heterogeneous across the genome, with a particularly wide differentiation peak on the sex chromosome Z. We further tested for preferential introgression of mitonuclear genes and found statistical support for this process in yellowhammers. A role for mitonuclear coevolution in this system is supported by a stronger signal of co-introgression in genes coding for subunits of the mitoribosome and electron transport chain complexes. Altogether, our study emphasizes the extreme variation seen in differentiation across genomic components and study systems as well as highlights the ramifications of this variation in species evolution.


title: Extreme Heterogeneity in Genomic Differentiation between Phenotypically Divergent Songbirds: A Test of Mitonuclear Co-introgression


title: Establishment of a Prognostic Prediction and Drug Selection Model for Patients with Clear Cell Renal Cell Carcinoma by Multi-Omics Data Analysis

abstract:
RationalePatients with clear cell renal cell cancer (ccRCC) may have completely different treatment choices and prognoses due to the wide range of heterogeneity of the disease. However, there is a lack of effective models for risk stratification, treatment decision making and prognostic prediction of renal cancer patients. The aim of the present study was to establish a model to stratify ccRCC patients in terms of prognostic prediction and drug selection based on multi-omics data analysis.

MethodsThis study was based on the multi-omics data (including mRNA, lncRNA, miRNA, methylation and WES) of 258 ccRCC patients from TCGA database. Firstly, we screened the feature values that had impact on the prognosis and obtained two subtypes. Then, we used 10 algorithms to achieve multi-omics clustering, and conducted pseudo-timing analysis to further validate the robustness of our clustering method, based on which the two subtypes of ccRCC patients were further subtyped. Meanwhile, the immune infiltration was compared between the two subtypes, and drug sensitivity and potential drugs were analyzed. Furthermore, to analyze the heterogeneity of patients at the multi-omics level, biological functions between two subtypes were compared. Finally, Boruta and PCA methods were used for dimensionality reduction and cluster analysis to construct a renal cancer risk model based on mRNA expression.

ResultsA prognosis predicting model of ccRCC was established by dividing patients into high- and low-risk groups. It was found that overall survival (OS) and progression-free interval (PFI) were significantly different between the two groups (p<0.01). The area under the OS time dependent ROC curve for 1, 3, 5 and 10 years in the training set was 0.75, 0.72, 0.71 and 0.68 respectively.

ConclusionThe model could precisely predict the prognosis of ccRCC patients and may have implications for drug selection for ccRCC patients.


title: The initial response of females towards congeneric males matches the propensity to hybridize in Ophthalmotilapia.

abstract: AO_SCPLOWBSTRACTC_SCPLOWCichlid radiations often harbour closely related species with overlapping niches and distribution ranges. Such species sometimes hybridise in nature, which raises the question how can they coexist. This also holds for the Tanganyika mouthbrooders Ophthalmotilapia ventralis and O. nasuta. Earlier studies found indications of asymmetrical hybridisation with females of O. ventralis accepting males of O. nasuta, but not the other way around. We hypothesised that this was due to differences in the capacity for species recognition. Given the higher propensity of O. ventralis females towards hybridisation, we expect a reduced ability for species recognition in O. ventralis females, compared to O. nasuta females. We staged two experiments, one focusing on 22 female O. nasuta and one on 21 female O. ventralis. These fish were placed in one half of a tank and briefly exposed to a conspecific or a heterospecific male, a conspecific female, or nothing (control). Female response was evaluated by scoring six tracking parameters and by noting the occurrence of ten discrete behaviours before and during the encounter. Females always responded to the presence of another fish by approaching it. Remarkably, for both O. nasuta and O. ventralis, we did not find a different response between encounters with conspecific males and females. However, in agreement with our hypothesis, O. nasuta females behaved differently towards conspecific or heterospecific males, whereas O. ventralis females did not. When presented with a heterospecific male, O. nasuta females performed a lower number of  ram behaviours. Additionally, they never displayed the  flee behaviour, a component of the species mating repertoire that was seen in all but one of the presentations with a conspecific male. Our findings show that differences in species recognition at first encounter predict to a large degree the outcome of the mating process, even in the absence of mating behaviour.


title: The initial response of females towards congeneric males matches the propensity to hybridise in Ophthalmotilapia.


title: Automated tracking reveals the social network of beach mice and their burrows

abstract: Evolutionary biologists have long sought to understand the selective pressures driving phenotypic evolution. While most experimental data come from the study of morphological evolution, we know much less about the ultimate drivers of behavioral variation. Among the most striking examples of behavioral evolution are the long, complex burrows constructed by oldfield mice (Peromyscus polionotus ssp.). Yet how these mice use burrows in the wild, and whether burrow length may affect fitness, remains unknown. A major barrier to studying behavior in the wild has been the lack of technologies to continuously monitor - in this case, nocturnal and underground - behavior. Here, we designed and implemented a novel radio frequency identification (RFID) system to track patterns of burrow use in a natural population of beach mice. We combine RFID monitoring with burrow measurements, genetic data, and social network analysis to uncover how these monogamous mice use burrows under fully natural ecological and social conditions. We first found that long burrows provide a more stable thermal environment and have higher juvenile activity than short burrows, underscoring the likely importance of long burrows for rearing young. We also find that adult mice consistently use multiple burrows throughout their home range and tend to use the same burrows at the same time as their genetic relatives, suggesting that inclusive fitness benefits may accrue for individuals that construct and maintain multiple burrows. Our study highlights how new automated tracking approaches can provide novel insights into animal behavior in the wild.


title: CRISPR-mediated insertion of a chimeric antigen receptor produces nonviral T cell products capable of inducing solid tumor regression

abstract:
BackgroundChimeric antigen receptor (CAR) T cells traditionally harbor viral vector-based sequences that encode the CAR transgene in the genome. These T cell products have yet to show consistent anti-tumor activity in patients with solid tumors. Further, viral vector manufacturing is resource intensive, suffers from batch-to-batch variability, and includes several animal components, adding regulatory and supply chain pressures.

MethodsAnti-GD2 CAR T cells were generated using CRISPR/Cas9 within nine days using recombinant Cas9 protein and nucleic acids, without any viral vectors or animal components. The CAR was specifically targeted to the T Cell Receptor Alpha Constant gene (TRAC). T cell products were characterized at the level of the genome, transcriptome, proteome, and secretome using CHANGE-seq, scRNA-seq, spectral cytometry, and ELISA assays. Functionality was evaluated in vivo in an NSG xenograft neuroblastoma model.

ResultsIn comparison to traditional retroviral CAR T cells, virus-free CRISPR CAR (VFC-CAR) T cells exhibit TRAC-targeted genomic integration of the CAR transgene, elevation of transcriptional and protein characteristics associated with a memory phenotype, and low tonic signaling prior to infusion arising in part from the the knockout of the TCR. Upon exposure to the GD2 target antigen, anti-GD2 VFC-CAR T cells exhibited specific cytotoxicity against GD2+ cells in vitro and induced solid tumor regression in vivo, with robust homing, persistence, and low exhaustion against a human neuroblastoma xenograft model.

ConclusionsThis proof-of-principle study leveraging virus-free genome editing technology could enable flexible manufacturing of clinically relevant, high-quality CAR T cells to treat cancers, including solid tumors.


title: Production and characterization of virus-free, CRISPR-CAR T cells capable of inducing solid tumor regression


title: ChemPLAN-Net: A deep learning framework to find novel inhibitor fragments for proteins

abstract: Fragment-based drug design plays an important role in the drug discovery process by reducing the complex small-molecule space into a more manageable fragment space. We leverage the power of deep learning to design ChemPLAN-Net; a model that incorporates the pairwise association of physicochemical features of both the protein drug targets and the inhibitor and learns from thousands of protein co-crystal structures in the PDB database to predict previously unseen inhibitor fragments. Our novel protocol handles the computationally challenging multi-label, multi-class problem, by defining a fragment database and using an iterative featurepair binary classification approach. By training ChemPLAN-Net on available co-crystal structures of the protease protein family, excluding HIV-1 protease as a target, we are able to outperform fragment docking and recover the targets inhibitor fragments found in co-crystal structures or identified by in-vitro cell assays.


title: AMES: Automated evaluation of sarcomere structures in cardiomyocytes

abstract:
BackgroundArrhythmias are severe cardiac diseases and lethal if untreated. To serve as an in vitro drug testing option for anti-arrhythmic agents, cardiomyocytes are being generated in vitro from induced pluripotent stem cells (iPSCs). Unfortunately, these generated cardiomyocytes resemble fetal cardiac tissue rather than adult cardiomyocytes. An automated tool for an unbiased evaluation of cardiomyocytes would highly facilitate the establishment of new differentiation protocols to increase cellular maturity.

ResultsIn this work, a novel deep learning-based approach for this task is presented and evaluated. Different convolutional neural networks (CNNs) including 2D and 3D models were trained on fluorescence images of human iPSC-derived cardiomyocytes, which were rated based on their sarcomere content (sarcomerisation) and the orientation of sarcomere filaments (directionality) beforehand by a domain expert. The CNNs were trained to perform classifications on sarcomerisation, directionality ratings, and cell source, including primary adult and differentiated cardiomyocytes. The best accuracies are reached by a 3D model with a classification accuracy of about 90 % for sarcomerisation classification, 63 % for directionality classification, and 80 % for cell source classification. The trained models were additionally evaluated using two explanatory algorithms, IGrad and Grad-CAM. The heatmaps computed by those explainability algorithms show that the important regions in the image occur inside the cell and at the cellular borders for the classifier, and, therefore, validate the calculated regions.

ConclusionIn summary, we showed that cellular fluorescence images can be analyzed with CNNs and subsequently used to predict different states of sarcomere maturation. Our developed prediction tool AMES (https://github.com/maxhillemanns/AMES) can be used to make trustworthy predictions on the quality of a cardiomyocyte, which ultimately facilitates the optimized generation of cardiomyocytes from iPSCs and improves the quality control in an automated, unbiased manner. The applied workflow of testing different CNN models, adjusting parameters, and using a variety of explanatory algorithms can be easily transferred to further image based quality control, stratification, or analysis setups.


title: Predictive functionality of bacteria in naturally fermented milk products of India using PICRUSt2 and Piphillin pipelines

abstract: Naturally fermented milk (NFM) products are popular food delicacies in Indian states of Sikkim and Arunachal Pradesh. Bacterial communities in these NFM products of India were previously analysed by high-throughput sequence method. However, predictive gene functionality of NFM products of India has not been studied. In this study, raw sequences of NFM products of Sikkim and Arunachal Pradesh were accessed from MG-RAST/NCBI database server. PICRUSt2 and Piphillin tools were applied to study microbial functional gene prediction. MUSiCC-normalized KOs and mapped KEGG pathways from both PICRUSt2 and Piphillin resulted in higher percentage of the former in comparison to the latter. Though, functional features were compared from both the pipelines, however, there were significant differences between the predictions. Therefore, a consolidated presentation of both the algorithms presented an overall outlook into the predictive functional profiles associated with the microbiota of the NFM products of India.


title: Exploration of novel αβ-protein folds through de novo design

abstract: Most naturally occurring protein folds have likely been discovered1-3. The question is whether natural evolution has exhaustively sampled almost all possible protein folds4, or whether a large fraction of the possible folds remains unexplored5-7. To address this question, we introduce a set of rules for {beta}-sheet topology to predict novel folds, and carry out the systematic de novo protein design for the novel folds predicted by the rules. The rules predicted eight novel {beta}-folds with a four-stranded {beta}-sheet, including a knot-forming one. We designed proteins for all the predicted {beta}-folds and found that all the designs are monomeric with high thermal stability and fold into the structures close to the design models, demonstrating the ability of the set of rules to predict novel {beta}-folds. The rules also predicted about twelve thousand novel {beta}-folds with five- to eight-stranded {beta}-sheets; the number is far exceeding the number of {beta}-folds observed so far. This result suggests that the enormous number of {beta}-folds are possible but have not emerged or become extinct due to evolutionary bias. The predicted novel folds should open up the possibility of designing functional proteins of our interests.


title: Comparative untargeted metabolomic profiling of induced mitochondrial fusion in pancreatic cancer

abstract: Mitochondria are dynamic organelles that constantly alter their shape through the recruitment of specialized proteins, like mitofusin-2 (Mfn2) and dynamin-related protein 1 (Drp1). Mfn2 induces the fusion of nearby mitochondria, while Drp1 mediates mitochondrial fission. We previously found that the genetic or pharmacological activation of mitochondrial fusion was tumor suppressive against pancreatic ductal adenocarcinoma (PDAC) in several model systems [1]. The mechanisms of how these different inducers of mitochondrial fusion reduce pancreatic cancer growth are still unknown. Here, we characterized and compared the metabolic reprogramming of these three independent methods of inducing mitochondrial fusion in KPC cells: overexpression of Mfn2, genetic editing of Drp1, or treatment with leflunomide. We identified significantly altered metabolites via robust, orthogonal statistical analyses and found that mitochondrial fusion consistently produces alterations in the metabolism of amino acids. Our unbiased methodology revealed that metabolic perturbations were similar across all these methods of inducing mitochondrial fusion, proposing a common pathway for metabolic targeting with other drugs.


title: Ploidetect enables pan-cancer analysis of the causes and impacts of chromosomal instability

abstract: Cancers routinely exhibit chromosomal instability, resulting in the accumulation of changes in the abundance of genomic material, known as copy number variants (CNVs). Unfortunately, the detection of these variants in cancer genomes is difficult. We developed Ploidetect, a software package that effectively identifies CNVs within whole-genome sequenced tumors. Ploidetect was more sensitive to CNVs in cancer related genes within advanced, pre-treated metastatic cancers than other tools, while also segmenting the most contiguously. Chromosomal instability, as measured by segment contiguity, was associated with several biological and clinical variables, including tumor mutation burden, tumor type, duration of therapy and immune microenvironment, highlighting the relevance of measuring CNV across the cancer genome. Investigation of gene mutations in samples revealed and the mutation status of several genes including ROCK2 and AC074391.1. Leveraging our heightened ability to detect CNVs, we identified 282 genes which were recurrently homozygously deleted in metastatic tumors. Further analysis of one recurrently deleted gene, MACROD2, identified a putative fragile tumor suppressor locus associated with response to chromosomal instability and chemotherapeutic agents. Our results outline the multifaceted impacts of CNVs in cancer by providing evidence of their involvement in tumorigenic behaviors and their utility as biomarkers for biological processes. We propose that increasingly accurate determination of CNVs is critical for their productive study in cancer, and our work demonstrates advances made possible by progress in this regard.


title: Mathematical Modeling Reveals Quantitative Properties of KEAP1-NRF2 Signaling

abstract: In response to oxidative and electrophilic stresses, cells launch an NRF2-mediated transcriptional antioxidant program. The activation of NRF2 depends on a redox sensor, KEAP1, which acts as an E3-ligase adaptor to promote the ubiquitination and degradation of NRF2. While a great deal has been learned about this redox duo, its quantitative signaling properties are still largely unexplored. In the present study, we examined these properties including response time, half-life, maximal activation, and response steepness (ultrasensitivity) of NRF2, through mathematical modeling. The models describe, with increasing complexity, the reversible binding of KEAP1 dimer and NRF2 via the ETGE and DLG motifs, NRF2 production, KEAP1-dependent and independent NRF2 degradation, and perturbations by different classes of NRF2 activators. Our simulations revealed that at the basal condition, NRF2 molecules are largely sequestered by KEAP1, with the KEAP1-NRF2 complex distributed comparably in either an ETGE-bound only (open) state or an ETGE and DLG dual-bound (closed) state, corresponding to the unlatched and latched configurations of the conceptual hinge-latch model. With two-step ETGE binding, the open and closed states operate in cycle mode at the basal condition and transition to equilibrium mode at stressed conditions. Class I-V, electrophilic NRF2 activators, which modify redox-sensing cysteine residues of KEAP1, shift the balance to a closed state that is unable to degrade NRF2 effectively. When total NRF2 accumulates to a level that nearly saturates existing KEAP1, ultrasensitive NRF2 activation, i.e., a steep rise in the free NRF2 level, can occur. The ultrasensitivity results from two simultaneous mechanisms, zero-order degradation mediated by DLG binding and protein sequestration (molecular titration) mediated by ETGE binding. These response characteristics of class I-V activators do not require disruption of DLG binding to unlatch the KEAP1-NRF2 complex. In comparison, class VI NRF2 activators, which directly compete with NRF2 for KEAP1 binding, cause a shift to the unlatched, open state of KEAP1-NRF2 complex and ultimately its complete dissociation (unhinged), resulting in a fast release of free NRF2 followed by stabilization. Although class VI activators may induce free NRF2 to higher levels, ultrasensitivity is lost due to lower free KEAP1 and thus its NRF2-sequestering effect. Stress-induced nuclear NRF2 accumulation is enhanced when basal nuclear NRF2 turnover constitutes a small load of NRF2 production. Our simulation further demonstrated that optimal abundances of cytosolic and nuclear KEAP1 exist to maximize ultrasensitivity. In summary, by simulating the dual role of KEAP1 in repressing NRF2, i.e., sequestration and promoting degradation, our mathematical modeling provides key novel quantitative insights into the signaling properties of the KEAP1-NRF2 system, which may help with the design of novel classes of NRF2 activators and inhibitors and understanding of the toxic actions of environmental oxidative stressors.

category: systems biology


title: Charting Brain Growth and Aging at High Spatial Precision

abstract: Defining reference models for population variation, and the ability to study individual deviations is essential for understanding inter-individual variability and its relation to the onset and progression of medical conditions. In this work, we assembled a reference cohort of neuroimaging data from 82 sites (N=58,836; ages 2-100) and use normative modeling to characterize lifespan trajectories of cortical thickness and subcortical volume. Models are validated against a manually quality checked subset (N=24,354) and we provide an interface for transferring to new data sources. We showcase the clinical value by applying the models to a transdiagnostic psychiatric sample (N=1,985), showing they can be used to quantify variability underlying multiple disorders whilst also refining case-control inferences. These models will be augmented with additional samples and imaging modalities as they become available. This provides a common reference platform to bind results from different studies and ultimately paves the way for personalized clinical decision making.

category: neuroscience


title: Integrated single-cell transcriptomic and epigenetic analyses of cell-state transition and lineage commitment in the embryonic mouse cerebellum

abstract: Recent studies using single-cell RNA-seq have revealed cellular heterogeneity in the developing mammalian cerebellum, yet the regulatory logic underlying this cellular diversity remains to be elucidated. Using integrated single-cell RNA and ATAC analyses, we resolved developmental trajectories of cerebellar progenitors and identified putative trans- and cis- elements that control cell state transition. We reverse-engineered gene regulatory networks (GRNs) of each cerebellar cell type. Through in silico simulations and in vivo experiments, we validated the efficacy of GRN analyses and uncovered the molecular control of a newly identified stem zone, the posterior transitory zone (PTZ), which contains multipotent progenitors for granule neurons, Bergmann glia, and choroid plexus epithelium. Importantly, we showed that perturbing cell fate specification of PTZ progenitors causes posterior cerebellar vermis hypoplasia, the most common cerebellar birth defect in humans. Our study provides a foundation for comprehensive studies of developmental programs of the mammalian cerebellum.