title: Oocyte spindle assembly depends on multiple interactions between HP1 and the CPC

abstract: The chromosomes in the oocytes of many animals appear to promote bipolar spindle assembly. In Drosophila oocytes, spindle assembly requires the chromosome passenger complex (CPC), which consists of INCENP, Borealin, Survivin and Aurora B. To determine what recruits the CPC to the chromosomes and its role in spindle assembly, we developed a strategy to manipulate the function and localization of INCENP, which is critical for recruiting the Aurora B kinase. We found that an interaction between Borealin and the chromatin is crucial for the recruitment of the CPC to the chromosomes and is sufficient to build kinetochores and recruit spindle microtubules. We also found that HP1 moves from the chromosomes to the spindle microtubules along with the CPC. We propose that the interaction with HP1 promotes the movement of the CPC from the chromosomes to the microtubules. In addition, within the central spindle, rather than at the centromeres, the CPC and HP1 are required for homologous chromosome bi-orientation.


title: MINTIE: identifying novel structural and splice variants in transcriptomes using RNA-seq data

abstract:
Genomic rearrangements can modify gene function by altering transcript sequences, and have been shown to be drivers in both cancer and rare diseases. Although there are now many methods to detect structural variants from Whole Genome Sequencing (WGS), RNA sequencing (RNA-seq) remains under-utilised as a technology for the detection of gene altering structural variants. Calling fusion genes from RNA-seq data is well established, but other transcriptional variants such as fusions with novel sequence, tandem duplications, large insertions and deletions, and novel splicing are difficult to detect using existing approaches.

To identify all types of variants in transcriptomes, we developed MINTIE, an integrated pipeline for RNA-seq data. We take a reference free approach, which combines de novo assembly of transcripts with differential expression analysis, to identify up-regulated novel variants in a case sample.

We validated MINTIE on simulated and real data sets and compared it with eight other approaches for finding novel transcriptional variants. We found MINTIE was able to detect >85% of variants while no other method was able to achieve this.

We applied MINTIE to RNA-seq data from a cohort of acute lymphoblastic leukemia (ALL) patient samples and identified several clinically relevant variants, including a recurrent unpartnered fusion involving the tumour suppressor gene RB1, and variants in ALL-associated genes: tandem duplications in IKZF1 and PAX5, and novel splicing in ETV6. We further demonstrate the utility of MINTIE to identify rare disease variants using RNA-seq, including the discovery of an inter-chromosomal translocation in the DMD gene in a patient with muscular dystrophy. We posit that MINTIE will be able to identify new disease variants across a range of cancers and other disease types.

category: bioinformatics


title: The Sequenced Genomes of Non-Flowering Land Plants Reveal the (R)Evolutionary History of Peptide Signaling

abstract: An understanding of land plant evolution is a prerequisite for in-depth knowledge of plant biology. Here we extract and explore information hidden in the increasing number of sequenced plant genomes, from bryophytes to angiosperms, to elucidate a specific biological question - how peptide signaling evolved. To conquer land and cope with changing environmental conditions, plants have gone through transformations that must have required a revolution in cell-to-cell communication. We discuss peptides mediating endogenous and exogenous changes by interaction with receptors activating intracellular molecular signaling. Signaling peptides were discovered in angiosperms and operate in tissues and organs like flowers, seeds, vasculature, and 3D meristems that are not universally conserved across land plants. Nevertheless, orthologues of angiosperm peptides and receptors have been identified in non-flowering plants. These discoveries provoke questions regarding the co-evolution of ligands and their receptors, and whether de novo interactions in peptide signaling pathways may have contributed to generate novel traits in land plants. The answers to such questions will have profound implications for the understanding of evolution of cell-to-cell communication and the wealth of diversified terrestrial plants. Under this perspective we have generated, analyzed and reviewed phylogenetic, genomic, structural, and functional data to elucidate the evolution of peptide signaling.

category: evolutionary biology


title: Polygenic Scores for Cognitive Abilities and their Association with Different Aspects of General Intelligence - a Deep Phenotyping Approach

abstract: Intelligence is a highly polygenic trait and genome-wide association studies (GWAS) have identified thousands of DNA variants contributing with small effects. Polygenic scores (PGS) can aggregate those effects for trait prediction in independent samples. As large-scale light-phenotyping GWAS operationalized intelligence as performance in rather superficial tests, the question arises which intelligence facets are actually captured. We used deep-phenotyping to investigate the molecular determinants of individual differences in cognitive ability. We therefore studied the association between PGS of intelligence (IQ-PGS), cognitive performance (CP-PGS) and educational attainment (EA-PGS) with a wide range of intelligence facets in a sample of 557 healthy adults. IQ-PGS, CP-PGS and EA-PGS had the highest incremental R2s for general (2.71%; 4.27%; 2.06%), verbal (3.30%; 4.64%; 1.61%) and numerical intelligence (3.06%; 3.24%; 1.26%) and the weakest for non-verbal intelligence (0.89%; 1.47%; 0.70%) and memory (0.80%; 1.06%; 0.67%). These results indicate that PGS derived from light-phenotyping GWAS do not reflect different facets of intelligence equally well, and thus should not be interpreted as genetic indicators of intelligence per se. The findings refine our understanding of how PGS are related to other traits or life outcomes.

category: genetics


title: Information theoretic approaches to deciphering the neural code with functional fluorescence imaging

abstract: Information theoretic metrics have proven useful in quantifying the relationship between behaviorally relevant parameters and neuronal activity with relatively few assumptions. However, these metrics are typically applied to action potential recordings and were not designed for the slow timescales and variable amplitudes typical of functional fluorescence recordings (e.g. calcium imaging). The lack of research guidelines on how to apply and interpret these metrics with fluorescence traces means the neuroscience community has yet to realize the power of information theoretic metrics. Here, we used computational methods to create mock action potential traces with known amounts of information. From these, we generated fluorescence traces and examined the ability of different information metrics to recover the known information values. We provide guidelines for how to use information metrics when applying them to functional fluorescence and demonstrate their appropriate application to GCaMP6f population recordings from mouse hippocampal neurons imaged during virtual navigation.

category: neuroscience


title: Genetic and species-level biodiversity patterns are linked by demography and ecological opportunity

abstract: The processes that give rise to species richness gradients are not well understood, but may be linked to resource-based limits on the number of species a region can support. Ecological limits placed on regional species richness would also limit population sizes, suggesting that these processes could also generate genetic diversity gradients. If true, we might better understand how broad-scale biodiversity patterns are formed by identifying the common causes of genetic diversity and species richness. We develop a hypothetical framework based on the consequences of regional variation in ecological limits to simultaneously explain spatial patterns of species richness and neutral genetic diversity. Repurposing raw genotypic data spanning 38 mammal species sampled across 801 sites in North America, we show that estimates of genome-wide genetic diversity and species richness share spatial structure. Notably, species richness hotspots tend to harbor lower levels of within-species genetic variation. A structural equation model encompassing eco-evolutionary processes related to resource availability, habitat heterogeneity, and human disturbance explained 78% of variation in genetic diversity and 74% of the variation in species richness. These results suggest we can infer broad-scale patterns of species and genetic diversity using two simple environmental measures of resource availability and ecological opportunity.


title: CAMAP: Artificial neural networks unveil the role of codon arrangement in modulating MHC-I peptides presentation

abstract:
MHC-I associated peptides (MAPs) play a central role in the elimination of virus-infected and neoplastic cells by CD8 T cells. However, accurately predicting the MAP repertoire remains difficult, because only a fraction of the transcriptome generates MAPs. In this study, we investigated whether codon arrangement (usage and placement) regulates MAP biogenesis. We developed an artificial neural network called Codon Arrangement MAP Predictor (CAMAP), predicting MAP presentation solely from mRNA sequences flanking the MAP-coding codons (MCCs), while excluding the MCC per se. CAMAP predictions were significantly more accurate when using original codon sequences than shuffled codon sequences which reflect amino acid usage. Furthermore, predictions were independent of mRNA expression and MAP binding affinity to MHC-I molecules and applied to several cell types and species. Combining MAP ligand scores, transcript expression level and CAMAP scores was particularly useful to increaser MAP prediction accuracy. Using an in vitro assay, we showed that varying the synonymous codons in the regions flanking the MCCs (without changing the amino acid sequence) resulted in significant modulation of MAP presentation at the cell surface. Taken together, our results demonstrate the role of codon arrangement in the regulation of MAP presentation and support integration of both translational and post-translational events in predictive algorithms to ameliorate modeling of the immunopeptidome.

Author summaryMHC-I associated peptides (MAPs) are small fragments of intracellular proteins presented at the surface of cells and used by the immune system to detect and eliminate cancerous or virus-infected cells. While it is theoretically possible to predict which portions of the intracellular proteins will be naturally processed by the cells to ultimately reach the surface, current methodologies have prohibitively high false discovery rates. Here we introduce an artificial neural network called Codon Arrangement MAP Predictor (CAMAP) which integrates information from mRNA-to-protein translation to other factors regulating MAP biogenesis (e.g. MAP ligand score and transcript expression levels) to improve MAP prediction accuracy. While most MAP predictive approaches focus on MAP sequences per se, CAMAPs novelty is to analyze the MAP-flanking mRNA sequences, thereby providing completely independent information for MAP prediction. We show on several datasets that the integration of CAMAP scores with other known factors involved in MAP presentation (i.e. MAP ligand score and mRNA expression) significantly improves MAP prediction accuracy, and further validate CAMAP learned features using an in-vitro assay. These findings may have major implications for the design of vaccines against cancers and viruses, and in times of pandemics could accelerate the identification of relevant MAPs of viral origins.


title: Brief sensory deprivation triggers plasticity of dopamine-synthesising enzyme expression in genetically labelled olfactory bulb dopaminergic neurons

abstract: In the glomerular layer of the olfactory bulb, local dopaminergic interneurons play a key role in regulating the flow of sensory information from nose to cortex. These dual dopamine- and GABA-releasing cells are capable of marked experience-dependent changes in the expression of neurotransmitter-synthesising enzymes, including tyrosine hydroxylase (TH). However, such plasticity has most commonly been studied in cell populations identified by their expression of the enzyme being studied, and after long periods of sensory deprivation. Here, instead, we used brief 1- or 3-day manipulations of olfactory experience in juvenile mice, coupled with a conditional genetic approach that labelled neurons contingent upon their expression of the dopamine transporter (DAT-tdTomato). This enabled us to evaluate the potential for rapid changes in neurotransmitter-synthesising enzyme expression in an independently identified neuronal population. Our labelling strategy showed good specificity for olfactory bulb dopaminergic neurons, whilst revealing a minority sub-population of non-dopaminergic DAT-tdTomato cells that expressed the calcium-binding protein calretinin. Crucially, the proportions of these neuronal subtypes were not affected by brief alterations in sensory experience. Short-term olfactory manipulations also produced no significant changes in immunofluorescence or whole-bulb mRNA for the GABA-synthesising enzyme GAD67/Gad1. However, in bulbar DAT-tdTomato neurons brief sensory deprivation was accompanied by a transient, small drop in immunofluorescence for the dopamine-synthesising enzyme dopa decarboxylase (DDC), and a sustained decrease for TH. Deprivation also produced a sustained decrease in whole-bulb Th mRNA. Careful characterisation of an independently identified, genetically labelled neuronal population therefore enabled us to uncover rapid experience-dependent changes in dopamine-synthesising enzyme expression.


title: Multi-Omic Analysis Reveals Cannabidiol Disruption of Cholesterol Homeostasis in Human Cell Lines

abstract: The non-psychoactive cannabinoid, cannabidiol (CBD), is FDA-approved for treatment of two drug-resistant epileptic disorders, and is seeing increased use among the general public, yet the mechanisms that underlie its therapeutic effects and side-effect profiles remain unclear. Here, we report a systems-level analysis of CBD action in human cell lines using temporal multi-omic profiling. FRET-based biosensor screening revealed that CBD treatment resulted in a sharp rise in cytosolic calcium, and activation of AMPK and ERK kinases in human keratinocyte and neuroblastoma cell lines. CBD treatment led to alterations in the abundance of metabolites, mRNA transcripts, and proteins consistent with activation of cholesterol biosynthesis, transport and storage. We found that CBD rapidly incorporated into cellular membranes and altered cholesterol chemical activity, suggesting direct perturbation of cholesterol-dependent membrane properties. CBD treatment induced apoptosis in a dose-dependent manner in multiple human cell lines, which was rescued by inhibition of cholesterol synthesis, and potentiated by compounds that disrupt cholesterol trafficking and storage. Our data point to a pharmacological interaction of CBD with cholesterol homeostasis pathways, with potential implications in its therapeutic use.

category: systems biology


title: In silico analysis of imprinted gene expression in the mouse skin.

abstract: 1.Imprinted genes help mediate embryonic cell proliferation and differentiation, but their roles after birth are far less well understood. A subset of 16 imprinted gene network (IGN) genes is expressed at higher levels in stem cell progenitor cells of adult skeletal muscle and epidermis compared to their differentiated counterparts. While these genes function in muscle regeneration, their role in the skin is poorly understood. We assessed the expression profiles of these 16 IGN genes in publicly available datasets and revealed elevated expression of IGN genes in the telogen and early anagen phases in mouse skin. We also identified IGN genes among a list of previously identified hair cycle-associated genes. Furthermore, our results suggest that IGN genes form part of a larger network and function predominantly as upstream regulators of hair cycle-regulated genes. Based on these in silico data, we propose a potential novel role of these 16 IGN genes as upstream regulators of hair cycle-associated genes. We speculate that IGN gene dysregulation participates in syndromes characterized by an impaired hair cycle. Thus, IGN gene expression might serve as a point of therapeutic intervention for patients suffering from cutaneous pathologies such as common hair-loss disorders.


title: SpikeShip: A method for fast, unsupervised discovery of high-dimensional neural spiking patterns

abstract: Neural coding and memory formation depend on temporal spiking sequences that span high-dimensional neural ensembles. The unsupervised discovery and characterization of these spiking sequences requires a suitable dissimilarity measure to compare any two spiking patterns. Here, we present a new measure of dissimilarity between multineuron spike sequences based on optimal transport theory, called SpikeShip. SpikeShip compares any two spiking patterns based on all the relative spike-timing relationships among neurons. Specifically, SpikeShip computes the optimal transport cost, i.e. the Earth Movers Distance, to make all the relative spike-timing relationships (across neurons) identical between two spiking patterns. It achieves this by computing the optimal transport of spikes across time, per neuron separately, and then decomposing this transport cost in a temporal rigid translation term and a vector of neuron-specific transport flows. We show that SpikeShip can be effectively computed for high-dimensional neuronal ensembles, has a low computational cost, and is explicitly based on the higher-order structure in the spiking patterns. We then applied SpikeShip to large-scale Neuropixel recordings during spontaneous activity and visual encoding. Using SpikeShip, we show that high-dimensional spiking sequences reliably distinguish between different natural images and different behavioral states. These spiking sequences carried complementary information to conventional firing rate codes. SpikeShip opens new avenues for studying neural coding and memory consolidation by finding patterns in high-dimensional neural ensembles.


title: Diverse viruses in deep-sea hydrothermal vent fluids have restricted dispersal across ocean basins

abstract:
In the ocean, viruses impact microbial mortality, regulate biogeochemical cycling, and alter the metabolic potential of microbial lineages. At deep-sea hydrothermal vents, abundant viruses infect a wide range of hosts among the archaea and bacteria that inhabit these dynamic habitats. However, little is known about viral diversity, host range, and biogeography across different vent ecosystems, which has important implications for how viruses manipulate microbial function and evolution. Here, we examined viral diversity, viral and host distribution, and viral-host interactions in microbial metagenomes generated from venting fluids from several vent sites within three different geochemically and geographically distinct hydrothermal systems: Piccard and Von Damm vent fields at the Mid-Cayman Rise in the Caribbean Sea, and at several vent sites within Axial Seamount in the Pacific Ocean. Analysis of viral sequences and Clustered Regularly InterSpaced Palindromic Repeats (CRISPR) spacers revealed highly diverse viral assemblages and evidence of active infection. Network analysis revealed that viral host range was relatively narrow, with very few viruses infecting multiple microbial lineages. Viruses were largely endemic to individual vent sites, indicating restricted dispersal, and in some cases viral assemblages persisted over time. Thus, we show that hydrothermal vent fluids are home to novel, diverse viral assemblages that are highly localized to specific regions and taxa.

ImportanceViruses play important roles in manipulating microbial communities and their evolution in the ocean, yet not much is known about viruses in deep-sea hydrothermal vents. However, viral ecology and evolution are of particular interest in hydrothermal vent habitats because of their unique nature: previous studies have indicated that most viruses in hydrothermal vents are temperate rather than lytic, and it has been established that rates of horizontal gene transfer (HGT) are particularly high among thermophilic vent microbes, and viruses are common vectors for HGT. If viruses have broad host range or are widespread across vent sites, they have increased potential to act as gene-sharing "highways" between vent sites. By examining viral diversity, distribution, and infection networks across disparate vent sites, this study provides the opportunity to better characterize and constrain the viral impact on hydrothermal vent microbial communities. We show that viruses in hydrothermal vents are diverse and apparently active, but most have restricted host range and are not widely distributed among vent sites. Thus, the impacts of viral infection are likely to be highly localized and constrained to specific taxa in these habitats.

category: microbiology


title: The evolution of Sox gene repertoires and regulation of segmentation in arachnids

abstract: The Sox family of transcription factors regulate many processes during metazoan development, including stem cell maintenance and nervous system specification. Characterising the repertoires and roles of these genes can therefore provide important insights into animal evolution and development. We further characterised the Sox repertoires of several arachnid species with and without an ancestral whole genome duplication (WGD), and compared their expression between the spider Parasteatoda tepidariorum and the harvestman Phalangium opilio. We also found that most Sox families have been retained as ohnologs after WGD and evidence for potential subfunctionalisation and/or neofunctionalization events. Our results also suggest that Sox21b-1 likely regulated segmentation ancestrally in arachnids, playing a similar role to the closely related SoxB gene, Dichaete, in insects. We previously showed that Sox21b-1 is required for the simultaneous formation of prosomal segments and sequential addition of opisthosomal segments in P. tepidariorum. We studied the expression and function of Sox21b-1 further in this spider and found that while this gene regulates the generation of both prosomal and opisthosomal segments, it plays different roles in the formation of these tagmata reflecting their contrasting modes of segmentation and deployment of gene regulatory networks with different architectures.

category: developmental biology


title: A RIF1/KAP1-based toggle switch stabilises the identities of the inactive and active X chromosomes during X inactivation

abstract: The onset of random X inactivation in mouse requires the switch from a symmetric to an asymmetric state, where the identities of the future inactive and active X chromosomes are assigned. Here we show that RIF1 and KAP1 are two fundamental factors for the definition of the asymmetry. Our data show that at the onset of mESC differentiation, upregulation of the long non-coding RNA Tsix weakens the symmetric RIF1 association with the Xist promoter, and opens a window of opportunity for a more stable association of KAP1. KAP1 is required to sustain high levels of Tsix, thus reinforcing and propagating the asymmetry, and, as a result, marking the future active X chromosome. Furthermore, we show that RIF1 association with the future inactive X chromosome is essential for Xist upregulation. This double-bookmarking system, based on the mutually exclusive relationships of Tsix and RIF1, and RIF1 and KAP1, thus coordinates the identification of the inactive and active X chromosomes and initiates a self-sustaining loop that transforms an initially stochastic event into a stably inherited asymmetric X chromosome state.

category: molecular biology


title: Fructose Metabolism Contributes to the Warburg effect

abstract: Fructose metabolism is increasingly recognized as a preferred energy source for cancer cell proliferation. However, it remains unclear why cancer cells favor fructose metabolism and how they acquire a sufficient amount of fructose. Here we report that cancer cells convert glucose into fructose through intra- and extracellular polyol pathways. The fructose metabolism bypasses normal aerobic respirations self-control to supply excessive metabolites to glycolysis and promotes the Warburg effect. Inhibition of fructose production drastically suppressed glycolysis and ATP production in cancer. Furthermore, we determined that a glucose transporter, SLC2A8/GLUT8, exports intracellular fructose to other cells in the tumor microenvironment. Taken together, our study suggests that the Warburg effect is achieved by means of fructose metabolism, instead of glucose metabolism alone.

category: cancer biology


title: Linear-Nonlinear Cascades Capture Synaptic Dynamics

abstract:
Short-term synaptic dynamics differ markedly across connections and strongly regulate how action potentials communicate information. To model the range of synaptic dynamics observed in experiments, we have developed a flexible mathematical framework based on a linear-nonlinear operation. This model can capture various experimentally observed features of synaptic dynamics and different types of heteroskedasticity. Despite its conceptual simplicity, we show that it is more adaptable than previous models. Combined with a standard maximum likelihood approach, synaptic dynamics can be accurately and efficiently characterized using naturalistic stimulation patterns. These results make explicit that synaptic processing bears algorithmic similarities with information processing in convolutional neural networks.

Author summaryUnderstanding how information is transmitted relies heavily on knowledge of the underlying regulatory synaptic dynamics. Existing computational models for capturing such dynamics are often either very complex or too restrictive. As a result, effectively capturing the different types of dynamics observed experimentally remains a challenging problem. Here, we propose a mathematically flexible linear-nonlinear model that is capable of efficiently characterizing synaptic dynamics. We demonstrate the ability of this model to capture different features of experimentally observed data.


title: Two FtsZ proteins orchestrate archaeal cell division through distinct functions in ring assembly and constriction

abstract: The tubulin homolog FtsZ assembles a cytokinetic ring in bacteria and plays a key role in the machinery that constricts to divide the cells. Many archaea encode two FtsZ proteins from distinct families, FtsZ1 and FtsZ2, of previously unclear functions. Here we show that Haloferax volcanii cannot divide properly without either or both, but DNA replication continues, and cells proliferate in alternative ways via remarkable envelope plasticity. FtsZ1 and FtsZ2 co-localize to form the dynamic division ring. However, FtsZ1 can assemble rings independently of FtsZ2, and stabilizes FtsZ2 in the ring, whereas FtsZ2 functions primarily in the constriction mechanism. FtsZ1 also influenced cell shape suggesting it forms a hub-like platform at midcell for the assembly of shape-related systems too. Both FtsZ1 and FtsZ2 are widespread in archaea with a single S-layer envelope, but archaea with a pseudomurein wall and division septum only have FtsZ1. FtsZ1 is therefore likely to provide a fundamental recruitment role in diverse archaea, and FtsZ2 is required for constriction of a flexible S-layer envelope, where an internal constriction force might dominate the division mechanism, in contrast to the single-FtsZ bacteria and archaea that divide primarily by wall ingrowth.


title: Acquisition of ionic copper by a bacterial outer membrane protein

abstract:
Copper, while toxic in excess, is an essential micronutrient in all kingdoms of life due to its essential role in the structure and function of many proteins. Proteins mediating ionic copper import have been characterised in detail for eukaryotes, but much less so for prokaryotes. In particular, it is still unclear whether and how Gram-negative bacteria acquire ionic copper. Here we show that Pseudomonas aeruginosa OprC is an outer membrane, TonB-dependent transporter that is conserved in many Proteobacteria and which mediates acquisition of both reduced and oxidised ionic copper via an unprecedented CxxxM-HxM metal binding site. Crystal structures of wild type and mutant OprC variants with silver and copper suggest that acquisition of Cu(I) occurs via a surface-exposed "methionine track" leading towards the principal metal binding site. Together with whole-cell copper quantitation and quantitative proteomics in a murine lung infection model, our data identify OprC as an abundant component of bacterial copper biology that may enable copper acquisition under a wide range of conditions.

SignificanceCopper is an essential metal in biology due to its role in the structure and function of many proteins. Despite this, it is not very clear how bacteria acquire copper, especially for Gram-negative organisms. In this study we show that the outer membrane protein OprC has an unusual metal binding site that allows OprC to bind both reduced and oxidised ionic copper near-irreversibly. Given the versatility of OprC, its presence in many Proteobacteria and its abundance during lung infection in mice, our study shows that OprC is an important component of prokaryote copper biology that warrants further study to uncover its regulation and to assess its role in bacterial virulence.


title: MYBL2 regulates ATM to control replication initiation and prevent replication stress in pluripotent stem cells

abstract: Replication stress, a major cause of genome instability in cycling cells, is mainly prevented by the ATR-dependent replication stress response pathway in somatic cells. However, the replication stress response pathway in embryonic stem cells (ESCs) may be substantially different due to alterations in cell cycle phase length. The transcription factor MYBL2, which is implicated in cell cycle regulation, is expressed between hundred to thousand-fold more highly in ESCs compared to somatic cells. Here we show that MYBL2 functions to activate ATM and suppress replication stress in ESCs. Consequently, loss of MYBL2 or inhibition of ATM or Mre11 in ESCs, results in replication fork slowing, increased fork stalling and elevated origin firing. Additionally, we demonstrate that inhibition of CDC7 activity rescues replication stress induced by MYBL2 loss and ATM inhibition, suggesting that uncontrolled new origin firing may underlie the replication stress phenotype resulting from loss/inhibition of MYBL2 and ATM. Overall, this study proposes that in addition to ATR, a MYBL2-MRN-ATM replication stress response pathway functions in ESCs to control DNA replication initiation and prevent genome instability.

category: cell biology


title: Signal Sequences Encode Information for Protein Folding in the Endoplasmic Reticulum

abstract: One-third of newly synthesized proteins in mammals are translocated into the endoplasmic reticulum (ER) through the Sec61 translocon. It remains unclear how protein translocation coordinates with the chaperone availability to promote protein folding in the ER. We find that signal sequences cause a translocation pause at the Sec61 translocon until nascent chains engage with luminal chaperones for efficient translocation and folding in the ER. Using a substrate-trapping proteomic approach, we identify that nascent proteins with marginally hydrophobic signal sequences accumulate on the cytosolic side of the Sec61 translocon, which recruits BiP chaperone through Sec63 to bind onto nascent chains. Surprisingly, BiP binding not only releases translocationally paused nascent chains into the ER lumen but also ensures protein folding. Increasing hydrophobicity of signal sequences bypasses Sec63/BiP-dependent protein translocation but translocated nascent chains misfold and aggregate under conditions of limited BiP availability in the ER. Thus, signal sequence-dependent protein folding explains why signal sequences are diverse and use multiple protein translocation pathways.


title: High Resolution and Multidimensional Protein Condensate Phase Diagrams with a Combinatorial Microdroplet Platform

abstract: The assembly of intracellular proteins into biomolecular condensates through phase separation is a fundamental process underlying organisation of the intracellular space and regulation of many central cellular processes. Elucidating the physico-chemical parameters that modulate phase separation is crucial in order to understand its effects in health and disease, as well as to develop therapeutic strategies that target biomolecular condensate systems. Phase behaviour is described by the phase diagram that defines conditions at which distinct phases exist, as well as the boundaries between these phases. Phase diagrams are typically built up from numerous individual measurements sampling different parts of the phase space. For biomolecular condensates, this process is slow, low throughput and associated with significant sample consumption even when performed in microwell format. Due to the rapidly increasing number of biologically and disease-relevant condensate systems, experimental techniques that enable high-throughput analysis of protein phase behaviour are required. To address this challenge, we present here a combinatorial droplet microfluidic platform, termed PhaseScan, for the rapid and high-resolution acquisition of protein phase diagrams in multiple chemical and physical dimensions. Using this platform, we characterise the phase behaviour of a wide range of protein systems under a variety of conditions and show that this platform can be used to characterise quantitatively the effect of small molecules on biomolecular phase transitions.

category: biophysics


title: A burning question: Can savannah fire management generate enough carbon revenue to help save the lion from extinction?

abstract:
Lions (Panthera leo) in Africa have lost nearly half their population in just the last two decades, and effective management of the protected areas (PAs) where lions live will cost an estimated USD >$1 B/year in new funding. We explore the potential for launching a fire management and habitat restoration carbon-financing program to help fill this PA management funding gap. We demonstrate how introducing early dry season fire management programs could produce potential carbon revenues (PCR) from either a single carbon-financing method (avoided emissions) or from multiple sequestration methods of USD $59.6-$655.9 M/year (at USD $5/ton) or USD $155.0 M-$1.7 B/year (at USD $13/ton). We highlight variable but significant PCR for PAs with the greatest potential for restoring lion numbers between USD $1.5-$44.4 M per PA. We suggest investing in lion-centric fire management programs to jump-start the United Nations Decade of Ecological Restoration and help preserve African lions across their range.

SCIENCE FOR SOCIETYThe United Nations recently launched the Decade of Ecological Restoration in response to planet-wide land degradation. This study analyses the potential for savanna fire management programs to restore fire regimes that can generate new sources of revenue from carbon financing for chronically under-funded protected areas in Africa with lions, as lions are a key indicator of savanna ecosystem health. We estimated the amount of carbon saved by shifting fires that normally burn in the late dry season (and emit more carbon) to the early dry season (that accrue more carbon in the soil and woody vegetation). Based on current carbon market values we found substantial potential to eliminate or significantly reduce the $>1B annual funding gap needed to save the lion from extinction. Given additional benefits for nature and people from new savanna fire programs, we recommend integrated conservation and development projects direct more funding to some of the least developed countries with high biodiversity and support fire management programs in Africa.

category: ecology


title: Biases in demographic modelling affect our understanding of recent divergence

abstract: Testing among competing demographic models of divergence has become an important component of evolutionary research in model and non-model organisms. However, the effect of unaccounted demographic events on model choice and parameter estimation remains largely unexplored. Using extensive simulations, we demonstrate that under realistic divergence scenarios, failure to account for population size (Ne) changes in daughter and ancestral populations leads to strong biases in divergence time estimates as well as model choice. We illustrate these issues reconstructing the recent demographic history of North Sea and Baltic Sea turbots (Schopthalmus maximus) by testing 16 Isolation with Migration (IM) and 16 Secondary Contact (SC) scenarios, modelling changes in Ne as well as the effects of linked selection and barrier loci. Failure to account for changes in Ne resulted in selecting SC models with long periods of isolation and divergence times preceding the formation of the Baltic Sea. In contrast, models accounting for Ne changes suggest recent (<6 kya) divergence with constant gene flow. We further show how interpreting genomic landscapes of differentiation can help discerning among competing models. For example, in the turbots data islands of differentiation show signatures of recent selective sweeps, rather than old divergence resisting secondary introgression. The results have broad implications for the study of population divergence by high-lighting the potential effects of unmodeleld changes in Ne on demographic inference. Tested models should aim at representing realistic divergence scenarios for the target taxa, and extreme caution should always be exercised when interpreting results of demographic modelling.


title: PIVOTAL: Prioritizing variants of uncertain significance with spatial genomic patterns in the 3D proteome

abstract: Variants of uncertain significance (VUS) have posed an increasingly prominent challenge to clinicians due to their growing numbers and difficulties in making clinical responses to them. Currently there are no existing methods that leverage the spatial relationship of known disease mutations and genomic properties for prioritizing variants of uncertain significance. More importantly, disease genes often associate with multiple clinically distinct diseases, but none of the existing variant prioritization methods provide clues as to the specific type of disease potentially associated with a given variant. We present PIVOTAL, a spatial neighborhood-based method using three-dimensional structural models of proteins, that significantly improves current variant prioritization tools and identifies potential disease etiology of candidate variants on a proteome scale. Using PIVOTAL, we made pathogenicity predictions for over 140,000 VUS and deployed a web application (http://pivotal.yulab.org) that enables users both to explore these data and to perform custom calculations.

category: genomics


title: Chaperone-mediated stress-sensing in Mycobacterium tuberculosis enables fast activation and sustained response

abstract:
Dynamical properties of gene-regulatory networks are tuned to ensure bacterial survival. In mycobacteria, MprAB-{sigma}E network responds to the presence of stressors, such as surfactants causing surface stress. Positive feedback loops in this network were previously predicted to cause hysteresis, i.e. different responses to identical stressor levels for pre-stressed and unstressed cells. Here we show that hysteresis does not occur in non-pathogenic Mycobacterium smegmatis but occurs in Mycobacterium tuberculosis. However, the observed rapid temporal response in M. tuberculosis is inconsistent with the model predictions. To reconcile these observations, we implement a recently proposed mechanism for stress-sensing: the release of MprB from the inhibitory complex with chaperone DnaK upon the stress exposure. Using modeling and parameter fitting, we demonstrate that this mechanism can accurately describe the experimental observations. Furthermore, we predict perturbations in DnaK expression that can strongly affect dynamical properties. Experiments with these perturbations agree with model predictions, confirming the role of DnaK in fast and sustained response.

IMPORTANCEGene-regulatory networks controlling stress response in mycobacterial species have been linked to persistence switches enabling the bacterial dormancy within a host. However, the mechanistic basis of switching and stress sensing is not fully understood. In this paper, combining quantitative experiments and mathematical modeling, we uncover how interactions between two master regulators of stress response -- MprAB two-component system and alternative sigma factor{sigma} E - shape the dynamical properties of surface-stress network. The result show hysteresis (history dependence) in the response of pathogenic bacteria M. tuberculosis to surface stress and lack of the hysteresis in non-pathogenic M. smegmatis. Furthermore, to resolve the apparent contradiction between the existence of hysteresis and fast activation of the response, we utilize a recently proposed role of chaperone DnaK in stress-sensing. The results leads to a novel system-level understanding of bacterial stress-response dynamics.


title: Influence of the large-Z effect during contact between butterfly sister species

abstract: Comparisons of genomes from recently diverged butterfly populations along a suture zone in central Texas have revealed high levels of divergence on the Z chromosome relative to autosomes, as measured by fixation index, Fst. The pattern of divergence appears to result from accumulation of incompatible alleles, obstructing introgression on the Z chromosome in hybrids. However, it is unknown whether this mechanism is sufficient to explain the data. Here, we simulate the effects of hybrid incompatibility on interbreeding butterfly populations using a model in which populations accumulate cross-incompatible alleles in allopatry prior to contact. We compute statistics for introgression and population divergence during contact between model butterfly populations and compare them to statistics obtained for 15 pairs of butterfly species interbreeding along the Texas suture zone. For populations that have evolved sufficiently in allopatry, the model exhibits high levels of divergence on the Z chromosome relative to autosomes in populations inter-breeding on time scales comparable to periods of interglacial contact between butterfly populations in central Texas. Levels of divergence on the Z chromosome increase when interacting groups of genes are closely linked, consistent with interacting clusters of functionally related genes in butterfly genomes. Results for various periods in allopatry are in qualitative agreement with the pattern of data for butterflies, supporting a picture of speciation in which populations are subjected to cycles of divergence in glacial isolation, and partial fusion during interglacial contact.


title: Abiotic treatment to common bean plants results in an altered seed microbiome.

abstract:
There has been a growing interest in the seed microbiome due to its important role as an end and starting point of plant microbiome assembly that can have consequences for plant health. However, the effect of abiotic conditions on the seed microbial community remains unknown. We performed a pilot study in a controlled growth chamber to investigate how the endophytic seed microbiome of the common bean (Phaseolus vulgaris L. (var. Red Hawk)) was altered under abiotic treatments relevant for crop management with changing climate. Bean plants were subjected to one of three treatments: 66% water withholding to simulate mild drought, 50% Hoagland nutrient solution to simulate fertilization, or control with sufficient water and baseline nutrition. We performed 16S rRNA gene amplicon sequencing and ITS1 amplicon sequencing of the endophytic DNA to assess seed bacterial/archaeal and fungal community structure, respectively. We found that variability in the seed microbiome structure was high while alpha diversity was low, with tens of taxa present. Water withholding and nutrient addition altered the seed microbiome structure for bacterial/archaeal communities as compared to the control, and each treatment resulted in a distinct microbiome structure. There were no statistically supported differences in the fungal microbiome across treatments. While we discuss several limitations of this study, the promising results suggest that further investigation is needed to better understand abiotic or stress-induced changes in the seed microbiome, the mechanisms that drive those changes, and their implications for the health and stress responses of the next plant generation.

ImportanceSeed microbiome members initiate the assembly of plant-associated microbial communities, but the environmental drivers of endophytic seed microbiome composition are unclear. Here, we exposed plants to short-term drought and fertilizer treatments during early vegetative growth and quantified the microbiome composition of the seeds that were ultimately produced. We found that seeds produced by plants stressed by water limitation or receiving nutrient addition had statistically different endophytic bacterial/archaeal microbiome compositions from each other and from seeds produced by control plants. This work suggests that the abiotic experience of a parental plant can influence the composition of its seed microbiome, with unknown consequences for the next plant generation.


title: Post-Translational Modifications Soften Vimentin Intermediate Filaments

abstract: The mechanical properties of biological cells are determined by the cytoskeleton, a composite biopolymer network consisting of microtubules, actin filaments and intermediate filaments (IFs). By differential expression of cytoskeletal proteins, modulation of the network architecture and interactions between the filaments, cell mechanics may be adapted to varying requirements on the cell. Here, we focus on the intermediate filament protein vimentin and introduce post-translational modifications as an additional, much faster mechanism for mechanical modulation. We study the impact of phosphorylation on filament mechanics by recording force-strain curves using optical traps. Partial phosphorylation softens the filaments. We show that binding of the protein 14-3-3 to phosphorylated vimentin IFs further enhances this effect and speculate that in the cell 14-3-3 may serve to preserve the softening and thereby the altered cell mechanics. We explain our observation by the additional charges introduced during phosphorylation.


title: Ancient viral genomes reveal introduction of HBV and B19V into Mexico during the transatlantic slave trade.

abstract: After the European colonization of the Americas there was a dramatic population collapse of the Indigenous inhabitants caused in part by the introduction of new pathogens. Although there is much speculation on the etiology of the Colonial epidemics, direct evidence for the presence of specific viruses during the Colonial era is lacking. To uncover the diversity of viral pathogens during this period, we designed an enrichment assay targeting ancient DNA (aDNA) from viruses of clinical importance and applied it on DNA extracts from individuals found in a Colonial hospital and a Colonial chapel (16th c. - 18th c.) where records suggest victims of epidemics were buried during important outbreaks in Mexico City. This allowed us to reconstruct three ancient human parvovirus B19 genomes, and one ancient human hepatitis B virus genome from distinct individuals. The viral genomes are similar to African strains, consistent with the inferred morphological and genetic African ancestry of the hosts as well as with the isotopic analysis of the human remains, suggesting an origin on the African continent. This study provides direct molecular evidence of ancient viruses being transported to the Americas during the transatlantic slave trade and their subsequent introduction to New Spain. Altogether, our observations enrich the discussion about the etiology of infectious diseases during the Colonial period in Mexico.


title: ANANSE: An enhancer network-based computational approach for predicting key transcription factors in cell fate determination

abstract: Proper cell fate determination is largely orchestrated by complex gene regulatory networks centered around transcription factors. However, experimental elucidation of key transcription factors that drive cellular identity is currently often intractable. Here, we present ANANSE (ANalysis Algorithm for Networks Specified by Enhancers), a network-based method that exploits enhancer-encoded regulatory information to identify the key transcription factors in cell fate determination. As cell type-specific transcription factors predominantly bind to enhancers, we use regulatory networks based on enhancer properties to prioritize transcription factors. First, we predict genome-wide binding profiles of transcription factors in various cell types using enhancer activity and transcription factor binding motifs. Subsequently, applying these inferred binding profiles, we construct cell type-specific gene regulatory networks, and then predict key transcription factors controlling cell fate transitions using differential networks between cell types. This method outperforms existing approaches in correctly predicting major transcription factors previously identified to be sufficient for trans-differentiation. Finally, we apply ANANSE to define an atlas of key transcription factors in 18 normal human tissues. In conclusion, we present a ready-to-implement computational tool for efficient prediction of transcription factors in cell fate determination and to study transcription factor-mediated regulatory mechanisms. ANANSE is freely available at https://github.com/vanheeringen-lab/ANANSE.


title: Dopaminergic Modulation of Working Memory and Cognitive Flexibility in a Zebrafish Model of Aging-Related Cognitive Decline

abstract: Healthy aging is associated with a decline in memory and executive function, which have both been linked with aberrant dopaminergic signalling. We examined the relationship between cognitive performance and dopamine function of young and aging zebrafish (Danio rerio). We revealed age-related decreases in working memory and cognitive flexibility in the Free-Movement Pattern (FMP) Y-maze. An increase in drd5 gene expression in aging adults coincided with a decrease in cognitive performance. Treatment with a D1/D5 receptor agonist (SKF-38393, 35 M) 30 minutes prior to behavioural assessment resulted in improved working memory in aging zebrafish, but no effect in younger adults. However, an  overdosing effect caused by agonist treatment resulted in downregulation of dat expression in 6-month old, treated zebrafish. The translational relevance of these findings was tested in humans by analysing exploratory behaviour in young-adult, 18-35-year olds, and aged adults, 70+ year olds, in a virtual FMP Y-maze. Our findings revealed similar age-related decline in working memory. Thus, strongly supporting zebrafish as a translational model of aging and cognitive decline.

category: animal behavior and cognition


title: Anesthetics fragment hippocampal network activity, alter spine dynamics and affect memory consolidation

abstract: General anesthesia is characterized by reversible loss of consciousness accompanied by transient amnesia. Yet, long-term memory impairment is an undesirable side-effect. How different types of general anesthetics (GAs) affect the hippocampus, a brain region central to memory formation and consolidation, is poorly understood. Using extracellular recordings, chronic 2-photon imaging and behavioral analysis, we monitor the effects of isoflurane (Iso), medetomidine/midazolam/fentanyl (MMF), and ketamine/xylazine (Keta/Xyl) on network activity and structural spine dynamics in the hippocampal CA1 area of adult mice. GAs robustly reduced spiking activity, decorrelated cellular ensembles, albeit with distinct activity signatures, and altered spine dynamics. CA1 network activity under all three anesthetics was different to natural sleep. Iso anesthesia most closely resembled unperturbed activity during wakefulness and sleep, and network alterations recovered more readily than with Keta/Xyl and MMF. Correspondingly, memory consolidation was impaired after exposure to Keta/Xyl and MMF, but not Iso. Thus, different anesthetics distinctly alter hippocampal network dynamics, synaptic connectivity, and memory consolidation, with implications for GA strategy appraisal in animal research and clinical settings.


title: ScRNA-seq discover cell cluster change under OAB: ACE2 expression reveal possible alternation of 2019-nCoV infectious pathway

abstract:
ObjectivePrevious study indicated that bladder cells which express ACE2 were a potential infection route of 2019-nCov. This study observed some differences of bladder cell cluster and their ACE2 expression between OAB mice and healthy mice, indicating the change of infectious possibility and pathway under overactive bladder (OAB) circumstance.

Material and methodPubic dataset acquisition was used to get ACE2 expression in normal human bladder and mice bladder (GSE129845). We built up over OAB model and studied the impact on cell typing and ACE2 expression. By way of using single-cell RNA sequencing (scRNA-seq) technique, bladder cell clustering and ACE2 expression in various cell types were measured respectively.

ResultIn pubic database (healthy human and mice bladder), ACE2 expression in humans and mice is concentrated in bladder epithelial cells. The disappearance of umbrella cells, a component of bladder epithelial, was found in our OAB model. In the two mouse bladder samples, ACE2 expression of epithelial cells is 34.1%, also the highest of all cell types.

ConclusionThe disappearance of umbrella cell may alternate the infection pathway of 2019-nCov and relate to the onset and progression of OAB.


title: Light sets the brain's daily clock by regional quickening and slowing of the molecular clockworks at dawn and dusk

abstract: How daily clocks in the brain are set by light to local environmental time and encode the seasons is not fully understood. The suprachiasmatic nucleus (SCN) is a central circadian clock that orchestrates physiology and behavior in tune with daily and seasonal light cycles. Here we have found that optogenetically simulated light input to explanted SCN changes the waveform of the molecular clockworks from sinusoidal in free-running conditions to highly asymmetrical shapes with accelerated synthetic (rising) phases and extended degradative (falling) phases marking clock advances and delays at simulated dawn and dusk. Daily waveform changes are sufficient to entrain to simulated winter and summer photoperiods, and to non-24h periods. SCN imaging further reveals that acute waveform shifts are greatest in the ventrolateral SCN, while period effects are greatest in the dorsomedial SCN. Thus, circadian entrainment is encoded by SCN clock gene waveform changes that arise from spatiotemporally distinct intrinsic responses within the SCN neural network.


title: Speech perception under the tent: A domain-general predictive role for the cerebellum

abstract: The role of the cerebellum in speech perception remains a mystery. Given its uniform architecture, we tested the hypothesis that it implements a domain-general mechanism whose role in speech is determined by connectivity. We collated all neuroimaging studies reporting cerebellar activity in the Neurosynth database (n = 8,206). From this set, we found all studies involving passive speech and sound perception (n = 72, 64% speech, 12.5% sounds, 12.5% music, and 11% tones) and speech production and articulation (n = 175). Standard and coactivation neuroimaging meta-analyses were used to compare cerebellar and associated cortical activations between passive perception and production. We found distinct regions of perception-and production-related activity in the cerebellum and regions of perception-production overlap. Each of these regions had distinct patterns of cortico-cerebellar connectivity. To test for domain generality versus specificity, we identified all psychological and task-related terms in the Neurosynth database that predicted activity in cerebellar regions associated with passive perception and production. Regions in the cerebellum activated by speech perception were associated with domain-general terms related to prediction. One hallmark of predictive processing is metabolic savings (i.e., decreases in neural activity when events are predicted). To test the hypothesis that the cerebellum plays a predictive role in speech perception, we examined cortical activation between studies reporting cerebellar activation and those without cerebellar activation during speech perception. When the cerebellum was active during speech perception there was far less cortical activation than when it was inactive. The results suggest that the cerebellum implements a domain-general mechanism related to prediction during speech perception.


title: Intracranial self-stimulation and concomitant behaviors following systemic methamphetamine administration in Hnrnph1 mutant mice

abstract:
RationaleAddiction to methamphetamine (MA) is a major public health issue in the United States. While psychostimulant use disorders are heritable, their genetic basis remains poorly understood. We previously identified heterogeneous nuclear ribonucleoprotein H1 (Hnrnph1; H1) as a quantitative trait gene underlying sensitivity to MA-induced locomotor activity. Mice heterozygous for a frameshift deletion in the first coding exon of H1 (H1+/-) showed reduced MA phenotypes including oral self-administration, locomotor activity, dopamine release, and dose-dependent differences in MA conditioned place preference. However, the effects of H1+/- on innate and MA-modulated reward sensitivity are not known.

ObjectivesWe examined innate reward sensitivity and modulation by MA in H1+/- mice via intracranial self-stimulation (ICSS).

MethodsWe used intracranial self-stimulation (ICSS) of the medial forebrain bundle to assess shifts in reward sensitivity following acute, ascending doses of MA (0.5-4.0 mg/kg, i.p.) using a within-subjects design. We also assessed video-recorded behaviors during ICSS testing sessions.

ResultsH1+/- mice displayed reduced normalized maximum response rates, H1+/- females showed lower normalized M50 values compared to wild-type females following MA, and H1+/- influenced ICSS responding relative to maximum baseline rates. There was a dose-dependent reduction in distance to the response wheel following MA administration, providing an additional measure of reward-related behavior.

ConclusionsH1+/- mice displayed reduced reward facilitation following MA in a sex- and dose-dependent manner. This result expands upon the set of MA-induced phenotypes observed in H1+/- mice.


title: A new pipeline to automatically segment and semi-automatically measure bone length on 3D models obtained by Computed Tomography

abstract:
The characterization of developmental phenotypes often relies on the accurate linear measurement of structures that are small and require laborious preparation. This is tedious and prone to errors, especially when repeated for the multiple replicates that are required for statistical analysis, or when multiple distinct structures have to be analysed. To address this issue, we have developed a pipeline for characterization of long-bone length and inter-vertebral distance using X-ray microtomography (XMT) scans. The pipeline involves semi-automated algorithms for automatic thresholding and fast interactive isolation and 3D-model generation of the main limb bones, using either the open-source ImageJ plugin BoneJ or the commercial Mimics Innovation Suite package. The tests showed the appropriate combination of scanning conditions and analysis parameters yields fast and comparable length results, highly correlated with the measurements obtained via ex vivo skeletal preparations. Moreover, since XMT is not destructive, the samples can be used afterwards for histology or other applications. Our new pipelines will help developmental biologists and evolutionary researchers to achieve fast, reproducible and non-destructive length measurement of bone samples from multiple animal species.

Summary statementBeltran Diaz et al. present a semi-automated pipeline for fast and versatile characterization of bone length from micro-CT images of mouse developmental samples.


title: Single-cell transcriptomic analysis of antiviral responses and viral antagonism in Chikungunya virus-infected synovial fibroblasts

abstract: In recent years, (re-)emerging arboviruses including Chikungunya virus (CHIKV) and Mayaro virus (MAYV) have caused growing concern due to expansion of insect vector ranges. No protective vaccine or specific antiviral strategies are currently available. Long-term morbidity after CHIKV infection includes debilitating chronic joint pain, which has associated health and economic impact. Here, we analyzed the early cell-intrinsic response to CHIKV and MAYV infection in primary human synovial fibroblasts. This interferon-competent cell type represents a potential source of polyarthralgia induced by CHIKV infection. Synovial fibroblasts from healthy and osteoarthritic donors were similarly permissive to CHIKV and MAYV infection ex vivo. Using RNA-seq, we defined a CHIKV infection-induced transcriptional profile with several hundred interferon-stimulated and arthralgia-mediating genes upregulated. Type I interferon was both secreted by infected fibroblasts and protective when administered exogenously. IL-6 secretion, which mediates chronic synovitis, however, was not boosted by infection. Single-cell RNA-seq and flow cytometric analyses uncovered an inverse correlation of activation of innate immunity and productive infection at the level of individual cells. In summary, primary human synovial fibroblasts serve as bona-fide ex vivo primary cell model of CHIKV infection and provide a valuable platform for studies of joint tissue-associated aspects of CHIKV immunopathogenesis.


title: A molecular phylogeny of forktail damselflies (genus Ischnura) reveals a dynamic macroevolutionary history of female colour polymorphisms

abstract: Colour polymorphisms are popular study systems among biologists interested in evolutionary dynamics, genomics, sexual selection and sexual conflict. In many damselfly groups, such as in the globally distributed genus Ischnura (forktails), sex-limited female colour polymorphisms occur in multiple species. Female-polymorphic species contain two or three female morphs, one of which phenotypically matches the male (androchrome or male mimic) and the other(s) which are phenotypically distinct from the male (heterochrome). These female colour polymorphisms are thought to be maintained by frequency-dependent sexual conflict, but their macroevolutionary histories are unknown, due to the lack of a robust molecular phylogeny. Here, we present the first time-calibrated phylogeny of Ischnura, using a multispecies coalescent approach (StarBEAST2) and incorporating both molecular and fossil data for 41 extant species (55% of the genus). We estimate the age of Ischnura to be between 13.8 and 23.4 millions of years, i.e. Miocene. We infer the ancestral state of this genus as female monomorphism with heterochrome females, with multiple gains and losses of female polymorphisms, evidence of trans-species female polymorphisms and a significant positive relationship between female polymorphism incidence and current geographic range size. Our study provides a robust phylogenetic framework for future research on the dynamic macroevolutionary history of this clade with its extraordinary diversity of sex-limited female polymorphisms.


title: The First complete Zoroastrian-Parsi Mitochondria Reference Genome: Implications of mitochondrial signatures in an endogamous, non-smoking population

abstract:
The present-day Zoroastrian-Parsis have roots in ancient pastoralist migrations from circumpolar regions leading to their settlement on the Eurasian Steppes and later, as Indo-Iranians in the Fertile Crescent. After migrating from the Persian province of Pars to India, the Zoroastrians from Pars ("Parsis") practiced endogamy, thereby preserving their genetic identity and social practices. The study was undertaken to gain an insight into the genetic consequences of migration on the community, the practice of endogamy, to decipher the phylogenetic relationships with other groups, and elucidate the disease linkages to their individual haplotypes

We generated the de novo the Zoroastrian-Parsi Mitochondrial Reference Genome (AGENOME-ZPMS-HV2a-1), which is the first complete mitochondrial reference genome assembled for this group. Phylogenetic analysis of an additional 99 Parsi mitochondrial genome sequences showed the presence of HV, U, T, A and F (belonging to the macrohaplogroup N) and Z and other M descendents of the macrohaplogroup M (M5, M39, M33, M4452, M24, M3, M30, M2, M430, M2, M35 and M27) and a largely Persian origin for the Parsi community. We assembled individual reference genomes for each major haplogroup and the Zoroastrian-Parsi Mitochondrial Consensus Genome (AGENOME-ZPMCG V1.0), which is the first consensus genome assembled for this group. We report the existence of 420 mitochondrial genetic variants, including 12 unique variants, in the 100 Zoroastrian-Parsi mitochondrial genome sequences. Disease association mapping showed 217 unique variants linked to longevity and 41 longevity-associated disease phenotypes across the majority of haplogroups.

Analysis of the coding genes, tRNA genes, and the D-loop region revealed haplogroup-specific disease associations for Parkinsons disease, Alzheimers disease, cancers, and rare diseases. No known mutations linked to lung cancer were found in our study. Mutational signatures linked to tobacco carcinogens, specifically, the C>A and G>T transitions, were observed at extremely low frequencies in the Parsi cohort, suggestive of an association between the cultural norm prohibiting smoking and its reflection in the genetic signatures. In sum, the Parsi mitochondrial genome provides an exceptional resource for determining details of their migration and uncovering novel genetic signatures for wellness and disease.


title: Explaining face representation in the primate brain using different computational models

abstract: Understanding how the brain represents the identity of complex objects is a central challenge of visual neuroscience. The principles governing object processing have been extensively studied in the macaque face patch system, a sub-network of inferotemporal (IT) cortex specialized for face processing. A previous study reported that single face patch neurons encode axes of a generative model called the "active appearance" model, which transforms 50-d feature vectors separately representing facial shape and facial texture into facial images. However, a systematic investigation comparing this model to other computational models, especially convolutional neural network models that have shown success in explaining neural responses in the ventral visual stream, has been lacking. Here, we recorded responses of cells in the most anterior face patch AM to a large set of real face images and compared a large number of models for explaining neural responses. We found that the active appearance model better explained responses than any other model except CORnet-Z, a feedforward deep neural network trained on general object classification to classify non-face images, whose performance it tied on some face image sets and exceeded on others. Surprisingly, deep neural networks trained specifically on facial identification did not explain neural responses well. A major reason is that units in the network, unlike neurons, are less modulated by face-related factors unrelated to facial identification such as illumination.


title: Connectivity dynamics of the Default Mode Network encode the stream of affect at multiple timescales

abstract:
In everyday life the stream of affect results from the interaction between past experiences, expectations, and the unfolding of events. How the brain represents the relationship between time and affect has been hardly explored, as it requires modeling the complexity of everyday life in the laboratory setting. Movies condense into hours a multitude of emotional responses, synchronized across subjects and characterized by temporal dynamics alike real-world experiences.

Here, we use time-varying intersubject brain synchronization and real-time behavioral reports to test whether connectivity dynamics track changes in affect during movie watching. Results show that polarity and intensity of experiences relate to connectivity of the default mode and control networks and converge in the right temporo-parietal cortex. We validate these results in two experiments including four independent samples, two movies, and alternative analysis workflows. Lastly, we reveal chronotopic connectivity maps within temporo-parietal and prefrontal cortex, where adjacent areas preferentially encode affect at specific timescales.


title: Default and Control networks connectivity dynamics track the stream of affect at multiple timescales


title: BioDynaMo: a general platform for scalable agent-based simulation

abstract:
MotivationAgent-based modeling is an indispensable tool for studying complex biological systems. However, existing simulators do not always take full advantage of modern hardware and often have a field-specific software design.

ResultsWe present a novel simulation platform called BioDynaMo that alleviates both of these problems. BioDynaMo features a general-purpose and high-performance simulation engine. We demonstrate that BioDynaMo can be used to simulate use cases in: neuroscience, oncology, and epidemiology. For each use case we validate our findings with experimental data or an analytical solution. Our performance results show that BioDynaMo performs up to three orders of magnitude faster than the state-of-the-art baseline. This improvement makes it feasible to simulate each use case with one billion agents on a single server, showcasing the potential BioDynaMo has for computational biology research.

AvailabilityBioDynaMo is an open-source project under the Apache 2.0 license and is available at www.biodynamo.org. Instructions to reproduce the results are available in supplementary information.

Contactlukas.breitwieser@inf.ethz.ch, a.s.hesam@tudelft.nl, omutlu@ethz.ch, r.bauer@surrey.ac.uk

Supplementary informationAvailable at https://doi.org/10.5281/zenodo.4501515.


title: Norepinephrine as a Memory Reset Signal: Switching the System from Retrieval to Encoding During a Spatial Memory Task can be Both Adaptive and Maladaptive

abstract:
Contextual information is represented in the hippocampus (HPC) partially through the recruitment of distinct neuronal ensembles. It is believed that reactivation of these ensembles underlies memory retrieval processes. Recently, we showed that norepinephrine (NE) input from phasic locus coeruleus (LC) activation induces hippocampal plasticity resulting in the recruitment of new neurons and a disengagement from previously established representations. We hypothesize that NE may provide a neuromodulatory, mnemonic switch signaling the HPC to move from a state of retrieval to encoding in the presence of novelty, and therefore, plays a role in memory updating. Here, we tested whether bilateral dorsal dentate gyrus (DG) infusions of the {beta}-adrenergic receptor (BAR) agonist isoproterenol (ISO), administered prior to encoding or retrieval, would impair spatial working and reference memory by reverting the system to encoding (thereby recruiting new neurons) potentially interfering with retrieval of the previously established spatial ensemble. We also investigated whether dDG infusions of ISO could promote cognitive flexibility by switching the system to encoding when it is adaptive (i.e. when new information is presented e.g. reversal learning). We found that intra-dDG infusions of ISO given prior to retrieval caused deficits in working and reference memory which was blocked by pre-treatment with the BAR-antagonist, propranolol (PRO). In contrast, ISO administered prior to reversal learning led to improved performance. These data support our hypothesis that NE serves as a novelty signal to update HPC contextual representations via BAR activation-facilitated recruitment of new neurons. This can be both maladaptive and adaptive depending on the situation.

SIGNIFICANCE STATEMENTThe current work highlights the involvement of hippocampal BARs in determining the flexibility of contextual representations to promote new learning in a way that supports adaptive behavior. This work builds upon previous work showing that noradrenergic input to the hippocampus is involved in recruiting new neurons resulting in new contextual representations and may be involved in the underlying neural mechanisms that support memory updating. These data suggest targets for anxiety disorders such as PTSD, which are characterized by noradrenergic dysregulation, and may also involve impairments in memory updating mechanisms where the incorporation of new information is not effectively encoded. The further understanding of the neurobiological mechanisms involved in updating memories may provide insight into novel treatment strategies.


title: Reconciling Conformational Heterogeneity and Substrate Recognition in Cytochrome P450

abstract:
Cytochrome P450, the ubiquitous metalloenzyme involved in detoxification of foreign components, has remained one of the most popular systems for substrate-recognition process. However, despite being known for its high substrate specificity, the mechanistic basis of substrate-binding by archetypal system cytochrome P450cam has remained at odds with the contrasting reports of multiple diverse crystallographic structures of its substrate-free form. Here we address this issue by elucidating the probability of mutual dynamical transition to the other crystallographic pose of cytochrome P450cam and vice versa via unbiased all-atom computer simulation. A robust Markov state model (MSM), constructed using adaptively sampled 84 microsecond-long Molecular dynamics simulation trajectories, maps the broad and heterogenous P450cam conformational landscape into five key sub-states. In particular, the MSM identifies an intermediate-assisted dynamic equilibrium between a pair of conformations of P450cam, in which the substrate-recognition sites remain  closed and  open respectively. However, the estimate of a significantly high stationary population of closed conformation, coupled with faster rate of open [->] closed transition than its reverse process, dictates that the net conformational equilibrium would be swayed in favour of  closed conformation. Together, the investigation quantitatively infers that while a potential substrate of cytochrome P450cam would in principle explore a diverse array of conformations of substrate-free protein, it would mostly encounter a  closed or solvent-occluded conformation and hence would follow an induced-fit based recognition process. Overall, the work reconciles multiple precedent crystallographic, spectroscopic investigations and establishes how a statistical elucidation of conformational heterogeneity in protein would provide crucial insights in the mechanism of potential substrate-recognition process.

STATEMENT OF SIGNIFICANCEConformational heterogeneity plays an important role in defining the structural and functional dynamics of the enzymes. While the static three-dimensional crystallographic structures of enzymes solved in different conditions and/or environments are crucial to provide the conformational sub-states of enzymes, these are not sufficient to understand the kinetics and thermodynamics of these sub-states and their role in substrate recognition process. Cytochrome P450cam, the archtypal metalloenzyme, presents such a complex scenario due to prevalent reports of contrasting crystallographic structures of its substrate-free form. This work quantifies the conformational heterogeneity of substrate-free P450cam by exploring the possibility of mutual transition among the crystallographic poses at an atomic resolution and in the process elucidates its possible substrate-recognition mechanism.

category: biochemistry


title: Impact of trimethoprim on the river microbiome and antimicrobial resistance.

abstract: Recent evidence suggests that anthropogenic activity can increase the levels of antimicrobial resistance (AMR) in the environment. Rivers and waterways are significant examples of environmental settings that have become repositories of antibiotics and antibiotic resistance genes (ARGs). Our recent study quantified drug concentrations in freshwater samples taken at a range of sites located on the Thames catchment; the highest levels of antibiotics and other drugs were recorded downstream of waste water treatment plants (WWTPs). One specific antibiotic: Trimethoprim (TMP) was shown at elevated concentrations reaching 2000ng/L at particular sites. We have also shown a correlative relationship between the residue of TMP and the prevalence of sulfonamide antibiotic resistance genes such as sul1. Despite this, there is still no evidence of a causative relationship between TMP concentrations and the prevalence of the ARGs at river sites. The aim of the current study was to conduct in-depth analysis using a combination of large metagenomic, geospatial and chemical datasets, in order to conduct a comparison between those sites with the highest TMP and lowest TMP levels across the Thames catchment. We aimed to establish the proximity of these sites to WWTPs, their population equivalence (PE) and land coverage. A secondary aim was to investigate seasonal variation in TMP and ARGs. Exploring these factors will help to decipher the clinical relevance of ARG accumulation at river sites. A significant correlation was shown between TMP levels at river sites and their distance downstream from a WWTP. Three sites located on the Rivers Cut and Ray showed significantly higher TMP concentrations in winter compared to summer. The population equivalence (PE) for sites with the highest TMP levels was significantly higher than those with the lowest levels. The land coverage of sites with the highest TMP levels was significantly more urban/suburban than sites with the lowest TMP concentrations, which were found to be significantly more arable. Five ARGs relevant to TMP and sulfonamides were identified across the Thames catchment. The most prevalent ARG was sul1, which was significantly more prevalent in winter compared to summer. By contrast sul2 was found to be significantly more prevalent in summer compared to winter at a site on the River Coln. The prevalence of the class 1 integron marker gene (inti1) did not differ significantly by season or between sites with the highest/lowest TMP levels.


title: The neural basis of intelligence in fine-grained cortical topographies

abstract: Intelligent thought is the product of efficient neural information processing, which is embedded in fine-grained, topographically-organized population responses and supported by fine-grained patterns of connectivity among cortical fields. Previous work on the neural basis of intelligence, however, has focused on coarse-grained features of brain anatomy and function, because cortical topographies are highly idiosyncratic at a finer scale, obscuring individual differences in fine-grained connectivity patterns. We used a computational algorithm, hyperalignment, to resolve these topographic idiosyncrasies, and found that predictions of general intelligence based on fine-grained (vertex-by-vertex) connectivity patterns were markedly stronger than predictions based on coarse-grained (region-by-region) patterns. Intelligence was best predicted by fine-grained connectivity in the default and frontoparietal cortical systems, both of which are associated with self-generated thought. Previous work overlooked fine-grained architecture because existing methods couldnt resolve idiosyncratic topographies, preventing investigation where the keys to the neural basis of intelligence are more likely to be found.


title: Neural systems underlying the learning of cognitive effort costs

abstract: People balance the benefits of cognitive work against the costs of cognitive effort. Models that incorporate prospective estimates of the costs of cognitive effort into decision making require a mechanism by which these costs are learned. However, it remains open what brain systems are important for this learning, particularly when learning is not tied explicitly to a decision about what task to perform. In this fMRI experiment, we parametrically manipulated the level of effort a task requires by increasing task switching frequency across six task contexts. In a scanned learning phase, participants implicitly learned about the task switching frequency in each context. In a subsequent test phase, participants made selections between pairs of these task contexts. We modeled learning within a reinforcement learning framework, and found that effort expectations that derived from task-switching probability and response time (RT) during learning were the best predictors of later choice behavior. Prediction errors (PE) from these two models were associated with FPN during distinct learning epochs. Specifically, PE derived from expected RT was most correlated with the fronto-parietal network early in learning, whereas PE derived from expected task switching frequency was correlated with the fronto-parietal network late in learning. These results suggest that multiple task-related factors are tracked by the brain while performing a task that can drive subsequent estimates of effort costs.


title: A consensus-based ensemble approach to improve transcriptome assembly

abstract:
BackgroundSystems-level analyses, such as differential gene expression analysis, co-expression analysis, and metabolic pathway reconstruction, depend on the accuracy of the transcriptome. Multiple tools exist to perform transcriptome assembly from RNAseq data. However, assembling high quality transcriptomes is still not a trivial problem. This is especially the case for non-model organisms where adequate reference genomes are often not available. Different methods produce different transcriptome models and there is no easy way to determine which are more accurate. Furthermore, having alternative-splicing events exacerbates such difficult assembly problems. While benchmarking transcriptome assemblies is critical, this is also not trivial due to the general lack of true reference transcriptomes.

ResultsIn this study, we first provide a pipeline to generate a set of the benchmark transcriptome and corresponding RNAseq data. Using the simulated benchmarking datasets, we compared the performance of various transcriptome assembly approaches including both de novo and genome-guided methods. The results showed that the assembly performance deteriorates significantly when alternative transcripts (isoforms) exist or for genome-guided methods when the reference is not available from the same genome. To improve the transcriptome assembly performance, leveraging the overlapping predictions between different assemblies, we present a new consensus-based ensemble transcriptome assembly approach, ConSemble.

ConclusionsWithout using a reference genome, ConSemble using four de novo assemblers achieved an accuracy up to twice as high as any de novo assemblers we compared. When a reference genome is available, ConSemble using four genome-guided assemblies removed many incorrectly assembled contigs with minimal impact on correctly assembled contigs, achieving higher precision and accuracy than individual genome-guided methods. Furthermore, ConSemble using de novo assemblers matched or exceeded the best performing genome-guided assemblers even when the transcriptomes included isoforms. We thus demonstrated that the ConSemble consensus strategy both for de novo and genome-guided assemblers can improve transcriptome assembly. The RNAseq simulation pipeline, the benchmark transcriptome datasets, and the script to perform the ConSemble assembly are all freely available from: http://bioinfolab.unl.edu/emlab/consemble/.


title: Model-based prediction of muscarinic receptor function from auditory mismatch negativity responses

abstract:
AO_SCPLOWBSTRACTC_SCPLOWDrugs affecting neuromodulation, for example by dopamine or acetylcholine, take centre stage among therapeutic strategies in psychiatry. These neuromodulators can change both neuronal gain and synaptic plasticity and therefore affect electrophysiological measures. An important goal for clinical diagnostics is to exploit this effect in the reverse direction, i.e., to infer the status of specific neuromodulatory systems from electrophysiological measures.

In this study, we provide proof-of-concept that the functional status of cholinergic (specifically muscarinic) receptors can be inferred from electrophysiological data using generative (dynamic causal) models. To this end, we used epidural EEG recordings over two auditory cortical regions during a mismatch negativity (MMN) paradigm in rats. All animals were treated, across sessions, with muscarinic receptor agonists and antagonists at different doses. Together with a placebo condition, this resulted in five levels of muscarinic receptor status. Using a dynamic causal model - embodying a small network of coupled cortical microcircuits - we estimated synaptic parameters and their change across pharmacological conditions. The ensuing parameter estimates associated with (the neuromodulation of) synaptic efficacy showed both graded muscarinic effects and predictive validity between agonistic and antagonistic pharmacological conditions.

This finding illustrates the potential utility of generative models of electrophysiological data as computational assays of muscarinic function. In application to EEG data of patients from heterogeneous spectrum diseases, e.g. schizophrenia, such models might help identify subgroups of patients that respond differentially to cholinergic treatments.

SO_SCPLOWIGNIFICANCEC_SCPLOW SO_SCPLOWTATEMENTC_SCPLOWIn psychiatry, the vast majority of pharmacological treatments affect actions of neuromodulatory transmitters, e.g. dopamine or acetylcholine. As treatment is largely trial-and-error based, one of the goals for computational psychiatry is to construct mathematical models that can serve as "computational assays" and infer the status of specific neuromodulatory systems in individual patients. This neuro-modelling strategy has great promise for electrophysiological data in particular but requires careful validation. The present study demonstrates that the functional status of cholinergic (muscarinic) receptors can be inferred from electrophysiological data using dynamic causal models of neural circuits. While accuracy needs to be enhanced and our results must be replicated in larger samples, our current results provide proof-of-concept for computational assays of muscarinic function using EEG.


title: Defective cyclophilin A induces TDP-43 proteinopathy: implications for amyotrophic lateral sclerosis and frontotemporal dementia

abstract: Aggregation and cytoplasmic mislocalization of TDP-43 are pathological hallmarks of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) spectrum. However, the molecular mechanism by which TDP-43 aggregates form and cause neurodegeneration remains poorly understood. Cyclophilin A, also known as peptidyl-prolyl cis-trans isomerase A (PPIA), is a foldase and molecular chaperone. We previously found that PPIA interacts with TDP-43 and governs some of its functions, and its deficiency accelerates disease in a mouse model of ALS. Here we characterized PPIA knock-out mice throughout their lifespan and found that they develop a neurodegenerative disease with key behavioural features of FTD, marked TDP-43 pathology and late-onset motor dysfunction. In the mouse brain, deficient PPIA induces aggregation of the GTP-binding nuclear protein Ran, a PPIA substrate required for TDP-43 nucleocytoplasmic trafficking. Moreover, in absence of PPIA, TDP-43 autoregulation is perturbed and TDP-43 and proteins involved in synaptic function are downregulated, leading to impairment of synaptic plasticity. Finally, we found that PPIA was downregulated in several ALS and ALS-FTD patients and identified a PPIA loss-of-function mutation in a sporadic ALS patient. The mutant PPIA has low stability, altered structure and impaired interaction with TDP-43. These findings strongly implicate that defective PPIA function causes TDP-43 mislocalization and dysfunction and should be considered in future therapeutic approaches.


title: Split-HaloTag Imaging Assay for Sophisticated Microscopy of Protein-Protein Interactions in planta

abstract:
An ever-increasing number of intracellular multi-protein networks have been identified in plant cells. Split-GFP based protein-protein interaction assays combine the advantages of in vivo interaction studies in a native environment with additional visualisation of protein complex localisation. Due to its simple protocols, it has become one of the most frequently used methods. However, standard fluorescent proteins entail several drawbacks for sophisticated microscopy.

With the HaloTag(R) system, these drawbacks can be overcome as this reporter forms covalent irreversible bonds with synthetic photostable fluorescent ligands. Dyes can be used in adjustable concentrations and are suitable for advanced microscopy methods. Therefore, we established the Split-HaloTag(R) imaging assay in plants which is based on the reconstitution of a functional HaloTag(R) protein upon protein-protein interaction and subsequent covalent binding of an added fluorescent ligand. Its suitability and robustness were demonstrated using well-characterised interactions as an example for protein-protein interaction at cellular structures: the molybdenum cofactor biosynthesis complex anchoring to filamentous actin. Additionally, a specific interaction was visualised with subdiffractional polarisation microscopy in a more distinctive manner as example for sophisticated imaging.

Split-GFP and Split-HaloTag(R) can complement one another as Split-HaloTag(R) represents an alternative option and an addition to the large toolbox of in vivo methods. Therefore, this promising new Split-HaloTag(R) imaging assay provides a unique and sensitive approach for more detailed characterization of protein-protein interaction with specific microscopic techniques such as 3D-imaging, single molecule tracking and super-resolution microscopy.

category: plant biology


title: Morpho-physiological and gene expression responses of wheat by Aegilops cylindrica amphidiploids to salt stress

abstract: Aegilops cylindrica Host is one of the most salt-tolerant species in the Triticeae tribe. Amphidiploid plants derived from hybridization of  Roshan x Aegilops cylindrica and  Chinese Spring x Ae. cylindrica genotypes contrasting in salt tolerance were assessed for their morpho-physiological responses and the expression patterns of three genes related to ion homeostasis under 250 mM NaCl. Results showed that salt stress caused significant declines in both their morphological and phenological traits. Moreover, salt stress reduced not only their chlorophyll content but also their root and shoot K contents and K/Na ratios, while it led to significant enhancements in the remaining traits. Similar to Ae. cylindrica, the amphidiploids subjected to salt stress exhibited significantly higher H2O2 levels, root and shoot K contents, and root and shoot K/Na ratios accompanied by lower root and shoot Na contents and MDA concentrations when compared with the same traits in the wheat parents. Quantitative Real-Time PCR showed significant differential expression patterns of the HKT1;5, NHX1, and SOS1 genes between the amphidiploids and their parents. The transcript level of HKT1;5 was found to be higher in the roots than in the shoots of both the amphidiploids and Ae. cylindrica while NHX1 exhibited a higher expression in the shoot tissues. The consistency of these data provides compelling support for the hypothesis that active exclusion of Na from the roots and elevated vacuolar sequestration of Na in the leaves might explain the declining Na levels in the shoots and roots of both the amphidiploids and Ae. cylindrica relative to those measured in wheat parents. It is concluded that the hybridized amphiploids are potentially valuable resources for salt improvement in bread wheat through the backcrossing approach.


title: In vitro endoderm emergence and self-organisation in the absence of extraembryonic tissues and embryonic architecture

abstract:
The endoderm is the cell lineage which gives rise in the embryo to the organs of the respiratory and gastrointestinal system. Uniquely, endodermal tissue does not just derive from descendants of the embryo proper (the epiblast) but instead arises from their gradual incorporation into an extraembryonic substrate (the visceral endoderm). Given the configuration of the early embryo, such a paradigm requires epiblast endodermal progenitors to negotiate embryonic compartments with very diverse epithelial character, a developmental contingency reflected by the fact that key early endodermal markers such as Foxa2 and Sox17 have been consistently found to be embedded within gene programmes involved in epithelialisation.

To explore the underlying cell biology of embryonic endoderm precursors, and to explore the relationship between endoderm development, epithelial identity, and extraembryonic mixing, we leveraged Gastruloids, in vitro models of early development. These self-organising three-dimensional aggregates of mouse embryonic stem cells do not possess an extraembryonic component, nor do they appear to display typical tissue architecture. Yet, they generate cells expressing endodermal markers. By tracking these cells throughout in vitro development, we highlight a persistent and uninterrupted pairing between epithelial and endodermal identity, with FoxA2+/Sox17+ endoderm progenitors never transitioning through mesenchymal intermediates and never leaving the epithelial compartment in which they arise. We also document the dramatic morphogenesis of these progenitors into a macroscopic epithelial primordium extending along the entire anterior-posterior axis of the Gastruloid. Finally, we find that this primordium correctly patterns into broad domains of gene expression, and matures cells with anterior foregut, midgut, and hindgut identities within 7 days of culture. We thus postulate that Gastruloids may serve as a potential source of endodermal types difficult to obtain through classical 2D differentiation protocols.


title: EEG in Motion: Using an Oddball Task to Explore Motor Interference in Active Skateboarding

abstract: Recent advancements in portable computer devices have opened new avenues in the study of human cognition outside research laboratories. This flexibility in methodology has led to the publication of several Electroencephalography (EEG) studies recording brain responses in real-world scenarios such as cycling and walking outside. In the present study, we tested the classic auditory oddball task while participants moved around an indoor running track using an electric skateboard. This novel approach allows for the study of attention in motion while virtually removing body movement. Using the skateboard auditory oddball paradigm, we found reliable and expected standard-target differences in the P3 and MMN/N2b event-related potentials (ERPs). We also recorded baseline EEG activity and found that, compared to this baseline, alpha power is attenuated in frontal and parietal regions during skateboarding. In order to explore the influence of motor interference in cognitive resources during skateboarding, we compared participants preferred riding stance (baseline level of riding difficulty) vs their non-preferred stance (increased level of riding difficulty). We found that an increase in riding difficulty did not modulate the P3 and tonic alpha amplitude during skateboard motion. These results suggest that increases in motor demands might not lead to reductions in cognitive resources as shown in previous literature.


title: A ChIP-exo screen of 887 PCRP transcription factor antibodies in human cells

abstract: Antibodies offer a powerful means to interrogate specific proteins in a complex milieu. However, antibody availability and reliability are problematic and epitope tagging can be impractical in many cases. In an effort to improve this situation, the Protein Capture Reagents Program (PCRP) generated over a thousand renewable monoclonal antibodies (mAbs) against human-presumptive chromatin proteins. However, these reagents have not been widely field-tested. We therefore performed a screen to test their ability to enrich genomic regions via chromatin immunoprecipitation (ChIP) and a variety of orthogonal assays. 887 unique antibodies against 681 unique human transcription factors (TFs), were assayed by ultra-high resolution ChIP-exo/seq, primarily in a single pass in one cell type (K562). Deep systematic analyses of the resulting [~]1,200 ChIP-exo datasets can be found at www.PCRPvalidation.org. Subsets of PCRP mAbs were further tested in ChIP-seq, CUT&RUN, STORM super-resolution microscopy, immunoblots, and protein binding microarray (PBM) experiments. About 5% of the tested antibodies displayed target (i.e., cognate antigen) enrichment across at least one assay and are strong candidates for additional validation. An additional 34% produced ChIP-exo data that was distinct from background and thus warrant further testing. The remaining 61% were not substantially different from background, and likely require consideration of a much broader survey of cell types and/or assay optimizations. We demonstrate and discuss the metrics and challenges to antibody validation in chromatin-based assays.