title: Plasmid analysis of NDM metallo-β-lactamase-producing Enterobacterales isolated in Vietnam

abstract: Carbapenem-resistant Enterobacterales (CRE) represent a serious threat to public health due to the lack of treatment and high mortality. The rate of antimicrobial resistance of Enterobacterales isolates to major antimicrobials, including carbapenems, is much higher in Vietnam than in Western countries, but the reasons remain unknown due to the lack of genomic epidemiology research. A previous study suggested that carbapenem resistance genes, such as the carbapenemase gene blaNDM, spread via plasmids among Enterobacterales in Vietnam. In this study, we characterized blaNDM-carrying plasmids in Enterobacterales isolated in Vietnam, and identified several possible cases of horizontal transfer of plasmids both within and among species of bacteria. Twenty-five carbapenem-nonsusceptible isolates from a medical institution in Hanoi were sequenced on Illumina short-read sequencers, and 13 blaNDM-positive isolates, including isolates of Klebsiella pneumoniae, Escherichia coli, Citrobacter freundii, Morganella morganii, and Proteus mirabilis, were further sequenced on an Oxford Nanopore Technologies long-read sequencer to obtain complete plasmid sequences. Almost identical 73 kb IncFII(pSE11)::IncN hybrid plasmids carrying blaNDM-1 were found in a P. mirabilis isolate and an M. morganii isolate. A 112 kb IncFII(pRSB107)::IncN hybrid plasmid carrying blaNDM-1 in an E. coli isolate had partially identical sequences with a 39 kb IncR plasmid carrying blaNDM-1 and an 88 kb IncFII(pHN7A8)::IncN hybrid plasmid in a C. freundii isolate. 148-149 kb IncFIA(Hl1)::IncA/C2 plasmids and 75-76 kb IncFII(Yp) plasmids, both carrying blaNDM-1 were shared among three sequence type 11 (ST11) isolates and three ST395 isolates of K. pneumoniae, respectively. Most of the plasmids co-carried genes conferring resistance to clinically relevant antimicrobials, including third-generation cephalosporins, aminoglycosides, and fluoroquinolones, in addition to blaNDM-1. These results provide insight into the genetic basis of CRE in Vietnam, and could help control nosocomial infections.


title: Generation of excitatory and inhibitory neurons from common progenitors via Notch signaling in the cerebellum

abstract:
Brain neurons arise from relatively few progenitors capable of giving rise to an enormous diversity of neuronal types. Nonetheless, a cardinal feature of mammalian brain neurogenesis in both the cortex and the cerebellum is that excitatory neurons and inhibitory neurons derive from separate, spatially segregated, progenitors. Whether bi-potential progenitors with an intrinsic capacity to generate both excitatory and inhibitory lineages exist and how such a fate decision may be regulated is unknown. Using cerebellar development as a model, we discover that individual embryonic cerebellar progenitors give rise to both inhibitory and excitatory lineages. We find that gradations of Notch activity levels determine the fates of the progenitors and their daughters. Daughters with the highest levels of Notch activity retain the progenitor fate. Daughters with intermediate levels of Notch activity become fate restricted to generate inhibitory neurons, while daughters with very low levels of Notch signaling adopt the excitatory fate. Therefore, Notch mediated binary cell fate choice is a mechanism for regulating the ratio of excitatory to inhibitory neurons from common progenitors.

Graphical summary

O_FIG O_LINKSMALLFIG WIDTH=144 HEIGHT=200 SRC="FIGDIR/small/997205v2_ufig1.gif" ALT="Figure 1">
View larger version (25K):
org.highwire.dtl.DTLVardef@c477f8org.highwire.dtl.DTLVardef@a079bforg.highwire.dtl.DTLVardef@1992877org.highwire.dtl.DTLVardef@2fd5de_HPS_FORMAT_FIGEXP  M_FIG C_FIG

category: developmental biology


title: Different temporal trends in vascular plant and bryophytes communities along elevational gradients over four decades of warming

abstract: Despite many studies showing biodiversity responses to warming, the generality of such responses across taxonomic groups remains unclear. Very few studies have tested for evidence of bryophyte community responses to warming, even though bryophytes are major contributors to diversity and functioning in many ecosystems. Here we report an empirical study comparing long-term change of bryophyte and vascular plant communities in two sites with contrasting long-term warming trends, using "legacy" botanical records as a baseline for comparison with contemporary resurveys. We hypothesized that ecological changes would be greater in sites with a stronger warming trend, and that vascular plant communities, with narrower climatic niches, would be more sensitive than bryophyte communities to climate warming. For each taxonomic group in each site, we quantified the magnitude of changes in species distributions along the elevation gradient, species richness, and community composition. We found contrasted temporal changes in bryophyte vs. vascular plant communities, which only partially supported the warming hypothesis. In the area with a stronger warming trend, we found a significant increase of local diversity and beta-diversity for vascular plants, but not for bryophytes. Presence absence data did not provide sufficient power to detect elevational shifts in species distributions. The patterns observed for bryophytes are in accordance with recent literature showing that local diversity can remain unchanged despite strong changes in composition. Regardless of whether one taxon is systematically more or less sensitive to environmental change than another, our results suggest that vascular plants cannot be used as a surrogate for bryophytes in terms of predicting the nature and magnitude of responses to warming. Thus, to assess overall biodiversity responses to global change, abundance data from different taxonomic groups and different community properties need to be synthesized.

category: ecology


title: The p53-p21 axis plays a central role in lymphatic homeostasis and disease

abstract:
Activation of the transcription factor p53 has been associated with several developmental syndromes. In normal tissues, p53 is kept at very low undetectable physiological levels. When triggered by cellular stressors, p53 prompts important anti-proliferative and apoptotic programs part of its tumor suppressor activity or as the guardian of tissue homeostasis.

We generated two murine models that display cutaneous hemorrhaging, severe edema, and distended blood-filled lymphatic vessels at late-gestation due to overactive p53 uniquely affecting lymphatic endothelial cells during development. Overactive p53 operated distinctively through anti-proliferative route in this tissue resulting in a decrease in initial lymphatics that normally absorb interstitial fluid. Remarkably, genetic or pharmacologic normalization of p53 restored lymphatic homeostasis and reversed lymphatic phenotypes. In parallel, several human lymphatic disease tissues exhibited high p53 levels exclusively in the lymphatic endothelium while p53 remained undetectable in surrounding arterial or venous vessels.

We report here, for the first time, an extended role that the p53 pathway plays in the genesis of lymphatic homeostasis deficiencies opening the way for new therapeutic avenues for these rare, poorly understood, and incurable lymphatic maladies.


title: Clinal and seasonal change are correlated in Drosophila melanogaster natural populations

abstract: Spatial and seasonal variation in the environment are ubiquitous. Environmental heterogeneity can affect natural populations and lead to covariation between environment and allele frequencies. Drosophila melanogaster is known to harbor polymorphisms that change both with latitude and seasons. Identifying the role of selection in driving these changes is not trivial, because non-adaptive processes can cause similar patterns. Given the environment changes in similar ways across seasons and along the latitudinal gradient, one promising approach may be to look for parallelism between clinal and seasonal change. Here, we test whether there is a genome-wide correlation between clinal and seasonal change, and whether the pattern is consistent with selection. Allele frequency estimates were obtained from pooled samples from seven different locations along the east coast of the US, and across seasons within Pennsylvania. We show that there is a genome-wide correlation between clinal and seasonal variation, which cannot be explained by linked selection alone. This pattern is stronger in genomic regions with higher functional content, consistent with natural selection. We derive a way to biologically interpret these correlations and show that around 3.7% of the common, autosomal variants could be under parallel seasonal and spatial selection. Our results highlight the contribution of natural selection in driving fluctuations in allele frequencies in natural fly populations and point to a shared genomic basis to climate adaptation which happens over space and time in D. melanogaster.

category: evolutionary biology


title: Meta-analysis suggests the microbiome responds to Evolve and Resequence experiments in Drosophila melanogaster

abstract:
BackgroundExperimental evolution has a long history of uncovering fundamental insights into evolutionary processes, but has largely neglected one underappreciated component--the microbiome. As eukaryotic hosts evolve, the microbiome may also evolve. However, the microbial contribution to host evolution remains poorly understood. Here, we re-analyzed genomic data to characterize the metagenomes from ten Evolve and Resequence (E&R) experiment in Drosophila melanogaster to determine how the microbiome changed in response to host selection.

ResultsBacterial diversity was significantly different in 5/10 studies, primarily in traits associated with metabolism or immunity. Duration of selection did not significantly influence bacterial diversity, highlighting the importance of associations with specific host traits. Additionally, we find that excluding reads from a facultative symbiont, Wolbachia, in the analysis of bacterial diversity changes the inference, raising important questions for future E&R experiments in the D. melanogaster microbiome.

ConclusionsOur genomic re-analysis suggests the microbiome often responds to host selection; thus, the microbiome may contribute to the response of Drosophila in E&R experiments. We outline important considerations for incorporating the microbiome into E&R experiments. The E&R approach may provide critical insights into host-microbiome interactions and fundamental insight into the genomic basis of adaptation.


title: The histone variant H2A.W and linker histone H1 co-regulate heterochromatin accessibility and DNA methylation

abstract: In flowering plants, heterochromatin is demarcated by the histone variant H2A.W, elevated levels of the linker histone H1, and specific epigenetic modifications, such as high levels of DNA methylation at both CG and non-CG sites. How H2A.W regulates heterochromatin organization and interacts with other heterochromatic features is unclear. To analyze the in vivo function of H2A.W, we created a h2a.w null mutant via CRISPR-Cas9, h2a.w-2. We find that H2A.W antagonizes deposition of H1 at heterochromatin and that non-CG methylation and accessibility are moderately decreased in h2a.w-2 heterochromatin. Compared to H1 loss alone, combined loss of H1 and H2A.W greatly increases accessibility and facilitates non-CG DNA methylation in heterochromatin, suggesting co-regulation of heterochromatic features by H2A.W and H1. Our results suggest that H2A.W helps maintain optimal heterochromatin accessibility and DNA methylation by promoting chromatin compaction together with H1, while also inhibiting excessive H1 incorporation.

category: genomics


title: Reward learning and working memory: effects of massed versus spaced training and post-learning delay period

abstract: Neuroscience research has illuminated the mechanisms supporting learning from reward feedback, demonstrating a critical role for the striatum and midbrain dopamine system. However, in humans, short-term working memory that is dependent on frontal and parietal cortices can also play an important role, particularly in commonly-used paradigms in which learning is relatively condensed in time. Given the growing use of reward-based learning tasks in translational studies in computational psychiatry, it is important to understand the degree of working memory contributions and whether gradual learning mechanisms can be better isolated. In our experiments, we manipulated the spacing between repetitions along with a post-learning delay preceding a test phase. We found that learning was slower for stimuli repeated after a long delay (spaced-trained) compared to those repeated immediately (massed-trained), likely reflecting the remaining contribution of feedback learning mechanisms when working memory is not available. Brief interruptions of massed learning led to drops in subsequent performance. Further, individual differences in working memory capacity positively correlated with massed learning performance. Critically, after a delay period but not immediately, relative preferences decayed in the massed condition and increased in the spaced condition. Overall, our results provide additional support for a large role of working memory in reward-based learning in temporally condensed designs. We suggest that spacing training within or between sessions is a promising approach to better isolate and understand mechanisms supporting gradual reward-based learning, with particular importance for understanding potential learning dysfunctions in addiction and psychiatric disorders.

category: neuroscience


title: celldeath: a tool for detection of cell death in transmitted light microscopy images by deep learning-based visual recognition

abstract: Cell death experiments are routinely done in many labs around the world, these experiments are the backbone of many assays for drug development. Cell death detection is usually performed in many ways, and requires time and reagents. However, cell death is preceded by slight morphological changes in cell shape and texture. In this paper, we trained a neural network to classify cells undergoing cell death. We found that the network was able to highly predict cell death after one hour of exposure to camptothecin. Moreover, this prediction largely outperforms human ability. Finally, we provide a simple python tool that can broadly be used to detect cell death.

category: cell biology


title: Genetic analyses identify widespread sex-differential participation bias

abstract: Genetic association results are often interpreted with the assumption that study participation does not affect downstream analyses. Understanding the genetic basis of this participation bias is challenging as it requires the genotypes of unseen individuals. However, we demonstrate that it is possible to estimate comparative biases by performing GWAS contrasting one subgroup versus another. For example, we show that sex exhibits autosomal heritability in the presence of sex-differential participation bias. By performing a GWAS of sex in ~3.3 million males and females, we identify over 158 autosomal loci significantly associated with sex and highlight complex traits underpinning differences in study participation between sexes. For example, the body mass index (BMI) increasing allele at the FTO locus was observed at higher frequency in males compared to females (OR 1.02 [1.02-1.03], P=4.4x10-36). Finally, we demonstrate how these biases can potentially lead to incorrect inferences in downstream analyses and propose a conceptual framework for addressing such biases. Our findings highlight a new challenge that genetic studies may face as sample sizes continue to grow.

category: genetics


title: The glycolysis regulator PFKFB4 interacts with ICMT and activates RAS/AKT signaling-dependent cell migration in melanoma.

abstract:
Cell migration is a complex process, tightly regulated during embryonic development and abnormally activated during cancer metastasis. RAS-dependent signaling is a major nexus controlling essential cell parameters including proliferation, survival and migration, utilizing downstream effectors such as the PI3K/AKT signaling pathway. In melanoma, oncogenic mutations frequently enhance RAS, PI3K/AKT or MAP kinase signaling, and trigger other cancer hallmarks among which the activation of metabolism regulators. PFKFB4 is one of these critical regulators of glycolysis and of the Warburg effect. Here however, we explore a novel function of PFKFB4 in melanoma cell migration. We find that PFKFB4 interacts with ICMT, a post-translational modifier of RAS. PFKFB4 promotes ICMT/RAS interaction, controls RAS localization at the plasma membrane, activates AKT signaling and enhances cell migration. We thus provide evidence of a novel and glycolysis-independent function of PFKFB4 in human cancer cells. This unconventional activity links the metabolic regulator PFKFB4 to RAS-AKT signaling and impacts melanoma cell migration.

Highlights- PFKFB4, a known regulator of glycolysis, also displays an unconventional role in melanoma cell migration.
- PFKFB4 interacts with ICMT and promotes RAS localization at the plasma membrane.
- PFKFB4 and ICMT cooperation modulates AKT signaling and controls melanoma cell migration.


title: A piRNA-lncRNA regulatory network initiates responder and trailer piRNA formation during mosquito embryonic development

abstract: PIWI-interacting (pi)RNAs are small silencing RNAs that are crucial for the defense against transposable elements in germline tissues of animals. In Aedes aegypti mosquitoes, the piRNA pathway also contributes to gene regulation in somatic tissues, illustrating additional roles for piRNAs and PIWI proteins besides transposon repression. Here, we identify a highly abundant endogenous piRNA (propiR1) that associates with both Piwi4 and Piwi5. PropiR1-mediated target silencing requires base pairing in the seed region with supplemental base pairing at the piRNA 3 end. Yet, propiR1 represses a limited set of targets, among which the lncRNA AAEL027353 (lnc027353). Slicing of lnc027353 initiates production of responder and trailer piRNAs from the cleavage fragment. Expression of propiR1 commences early during embryonic development and mediates degradation of maternally provided lnc027353. Both propiR1 and its lncRNA target are conserved in the closely related Aedes albopictus mosquito, underscoring the importance of this regulatory network for mosquito development.

category: molecular biology


title: Kappa but not delta or mu opioid receptors form homodimers at low membrane densities

abstract: Opioid receptors (ORs) have been observed as homo- and heterodimers, but it is unclear if the dimers are stable under physiological conditions, and whether monomers or dimers comprise the predominant fraction in a cell. Here we use three live-cell imaging approaches to assess dimerization of ORs at different expression levels. At high membrane densities, a split GFP assay reveals that {kappa}OR dimerizes, while OR and {delta}OR stay monomeric. In contrast, single-molecule imaging showed no {kappa}OR dimers at low receptor densities. To reconcile our seemingly contradictory results, we used a high-density single-molecule assay to assess membrane protein interactions at densities up to 100x higher than conventional single-molecule imaging. We observe that {kappa}OR is monomeric at low densities and forms dimers at densities that are considered physiological. In contrast, OR and {delta}OR stay monomeric even at the highest densities covered by our approach. The observation of long-lasting {kappa}OR dimers but not higher-order aggregates suggests that ORs dimerize through a single, specific interface.

category: biochemistry


title: Flagellar dynamics reveal fluctuations and kinetic limit in the Escherichia coli chemotaxis network

abstract: Biochemical signaling networks allow living cells to adapt to a changing environment, but these networks must cope with unavoidable number fluctuations ("noise") in their molecular constituents. Escherichia coli chemotaxis, by which bacteria modulate their random run/tumble swimming pattern to navigate their environment, is a paradigm for the role of noise in cell signaling. The key signaling protein, CheY, when activated by (reversible) phosphorylation, causes a switch in the rotational direction of the flagellar motors propelling the cell, leading to tumbling. CheY-P concentration, [CheY-P], is thus a measure of the chemotaxis networks output, and temporal fluctuations in [CheY-P] provide a proxy for network noise. Here we quantify the fluctuations in [CheY-P] from the switching dynamics of individual flagella, observed using time-resolved fluorescence microscopy of optically trapped E. coli cells. This approach reveals [CheY-P] fluctuations at steady state, which may play a critical role in driving flagellar switching and cell tumbling. A stochastic theoretical model, inspired by work on gene expression noise, points to CheY activation occurring in bursts, driving the [CheY-P] fluctuations. When the network is stimulated chemically to higher activity, we observe a dramatic decrease in [CheY-P] fluctuations. Our stochastic model shows that an intrinsic kinetic ceiling on network activity places an upper limit on [CheY-P], which when encountered suppresses its fluctuations. This limit may also prevent cells from tumbling unproductively in steep gradients.

category: biophysics


title: T606-phosphorylation deprives the function of Kaiso as a transcription repressor

abstract: It is well known that Kaiso protein encoded by ZBTB33 gene is a transcription repressor and that Kaiso-P120ctn interaction increases the shift of Kaiso from the nucleus into the cytoplasm. However, the regulatory mechanisms of Kaiso compartmentalization are far from clear. Here, we reported that AKT1 could phosphorylate 606-threonine residue (T606) within the RSSTIP motif of Kaiso in the cytoplasm. The T606-phosphorylated Kaiso (pT606-Kaiso) could directly bind to 14-3-3 family proteins and the depletion of T606 phosphorylation by T606A mutation abolished most of the Kaiso-14-3-3 binding. In addition, the Kaiso-P120ctn interaction was essential for the pT606-Kaiso accumulation in the cytoplasm. Notably, enforced 14-3-3{sigma} (SFN) overexpression could increase the pT606-Kaiso accumulation in the cytoplasm and de-repress the transcription of Kaiso target gene CDH1. Decreased amounts of both pT606-Kaiso and CDH1 proteins were frequently observed in human gastric cancer tissues relative to paired normal controls. The mRNA levels of 14-3-3{sigma} and Kaiso target gene CDH1 were positively and significantly correlated with each other in bioinformatics analyses using publicly available RNA-seq datasets for human normal tissues (n=11688, r=0.60, p<0.001) in the GTEx project and for cancer cell lines (n=1156, r=0.41, p<0.001) in the CCLE project. Furthermore, Kaiso T606A mutant (unable to be phosphorylated) significantly increased the migration and invasion of cancer cells in vitro as well as boosted the growth of these cells in vivo. In conclusion, Kaiso could be phosphorylated by AKT1 at the T606 and the pT606-Kaiso accumulates in the cytoplasm through binding to 14-3-3/P120ctn that de-represses the expression of Kaiso target gene CDH1 in normal tissues. Decreased Kaiso phosphorylation may contribute to the development of gastrointestinal cancer. The status of Kaiso phosphorylation is a determinant factor for the role of Kaiso in the development of cancer.


title: T606-phosphorylation deprives the function of Kaiso as a transcription factor and cancer repressor


title: T606-phosphorylation deprives the function of Kaiso as a transcription and oncogenic factor


title: Revisiting giraffe photo-identification using deep learning and network analysis

abstract:
An increasing number of ecological monitoring programs rely on photographic capture-recapture of individuals to study distribution, demography and abundance of species. Photo-identification of individuals can sometimes be done using idiosyncratic coat or skin patterns, instead of using tags or loggers. However, when performed manually, the task of going through photographs is tedious and rapidly becomes too time consuming as the number of pictures grows.

Computer vision techniques are an appealing and unavoidable help to tackle this apparently simple task in the big-data era. In this context, we propose to revisit animal re-identification using image similarity networks and metric learning with convolutional neural networks (CNNs), taking the giraffe as a working example.

We first developed an end-to-end pipeline to retrieve a comprehensive set of re-identified giraffes from about 4, 000 raw photographs. To do so, we combined CNN-based object detection, SIFT pattern matching, and image similarity networks. We then quantified the performance of deep metric learning to retrieve the identity of known individuals and detect unknown individuals never seen in the previous years of monitoring.

After a data augmentation procedure, the re-identification performance of the CNN reached a Top-1 accuracy of about 90%, despite the very small number of images per individual in the training data set. While the complete pipeline succeeded in re-identifying known individuals, it slightly under-performed with unknown individuals.

Fully based on open-source software packages, our work paves the way for further attempts to build automatic pipelines for re-identification of individual animals, not only in giraffes but also in other species.


title: HumGut: A comprehensive Human Gut prokaryotic genomes collection filtered by metagenome data

abstract:
BackgroundA major bottleneck in the use of metagenome sequencing for human gut microbiome studies has been the lack of a comprehensive genome collection to be used as a reference database. Several recent efforts have been made to re-construct genomes from human gut metagenome data, resulting in a huge increase in the number of relevant genomes. In this work, we aimed to create a collection of the most prevalent healthy human gut prokaryotic genomes, to be used as a reference database, including both MAGs from the human gut and ordinary RefSeq genomes.

ResultsWe screened > 5,700 healthy human gut metagenomes for the containment of > 490,000 publicly available prokaryotic genomes sourced from RefSeq and the recently announced UHGG collection. This resulted in a pool of > 379,000 genomes that were subsequently scored and ranked based on their prevalence in the healthy human metagenomes. The genomes were then clustered at subspecies resolution, and cluster representatives were retained to comprise the HumGut collection. Using the Kraken2 software for classification, we find superior performance in the assignment of metagenomic reads, classifying on average 94.5% of the reads in a metagenome, as opposed to 86% with UHGG and 44% when using standard Kraken2 database. HumGut, half the size of standard Kraken2 database and directly comparable to the UHGG size, outperforms them both.

ConclusionsThe HumGut collection contains > 30,000 genomes clustered at subspecies resolution and ranked by human gut prevalence. We demonstrate how metagenomes from IBD-patients map equally well to this collection, indicating this reference is relevant also for studies well outside the metagenome reference set used to obtain HumGut. We believe this is a valuable resource in a field in dire need of method standardization. All data and metadata, as well as helpful code, are available at http://arken.nmbu.no/~larssn/humgut/.

category: microbiology


title: Sponging of glutamate at the outer plasma membrane surface reveals roles for glutamate in development

abstract:
Plants use electrical and chemical signals for systemic communication. Herbivory, for instance, appears to trigger local apoplasmic glutamate accumulation, systemic electrical signals and calcium waves that travel to report insect damage to neighboring leaves and initiate defense. To monitor extra- and intracellular glutamate concentrations in plants, we generated Arabidopsis lines expressing genetically encoded fluorescent glutamate sensors. In contrast to cytosolically localized sensors, extracellularly displayed variants inhibited plant growth and proper development. Phenotypic analyses of high-affinity display sensor lines revealed that root meristem development, particularly the quiescent center (QC), number of lateral roots, vegetative growth and floral architecture were impacted. Notably, the severity of the phenotypes was positively correlated with the affinity of the display sensors, intimating that their ability to sequester glutamate at the surface of the plasma membrane was responsible for the defects. Root growth defects were suppressed by supplementing culture media with low levels of glutamate. Together, the data indicate that sequestration of glutamate at the cell surface either disrupts the supply of glutamate to meristematic cells and/or impairs localized glutamatergic signaling important for developmental processes.

SIGNIFICANCE STATEMENTGlutamate is an important signaling molecule in animals and plants. The affinity-dependent effects of surface-displayed glutamate sensors on plant development intimate that glutamate levels in the vicinity of the plasma membrane play important roles in signaling in plants.

category: plant biology


title: MeSSII interplays the amylose and amylopectin biosynthesis by protein interactions in cassava storage root

abstract:
Starch is a glucose polymer synthesized by green plants for energy storage, and is crucial for plant growth and reproduction. The biosynthesis of starch polysaccharides is mediated by members of the large starch synthase (SS) protein superfamily. Here, we showed that in cassava storage roots, soluble starch synthase II (MeSSII) plays an important role in starch biosynthesis via forming protein complexes with other starch biosynthetic enzymes by directly interacting with MeSSI, MeSBEII and MeISAII. The MeSSII-RNAi cassava lines showed increased amylose content and reduced intermediate chain of amylopectin (B1 type) biosynthesis in their storage roots, leading to altered starch physico-chemical properties. Further gel permeation chromatography analysis of starch biosynthetic enzymes between the wild type and MeSSII-RNAi lines confirmed the key rule of MeSSII in the organization of heteromeric starch synthetic protein complexes, including the MeSSII-MeSSI-MeGBSSI and MeSBEII-MeSSI-MeSSII-MeISAII-MeISAI complexes. The lack of MeSSII in cassava also reduced the binding capacity of the MeSSI, MeSBEII, MeISAI, and MeISAII to starch granules. Our results show a close coordination between granule-bound starch synthase and amylopectin biosynthetic enzymes, which implies that the processes of amylose synthesis and amylopectin synthesis are directly interrelated. These findings shed light on the key components of the starch biosynthesis machinery in root crops.

One Sentence SummaryHere, by focusing on cassava SSII, we elucidated its function and the molecular mechanism by which it chaperones the starch synthase complex to regulate starch biosynthesis in cassava storage roots.


title: The cytosolic Arabidopsis thaliana cysteine desulfurase ABA3 delivers sulfur to the sulfurtransferase STR18

abstract: The biosynthesis of many sulfur-containing molecules depends on cysteine as a sulfur source. Cysteine desulfurase (CD) and rhodanese (Rhd) domain-containing protein families participate in the trafficking of sulfur for various metabolic pathways in bacteria and human, but their connection is not yet described in plants. The existence of natural chimeric proteins, however, containing both CD and Rhd domains in specific bacterial genera suggests a general interaction between both proteins. We report here the biochemical relationships between two cytosolic proteins from Arabidopsis thaliana, a Rhd domain containing protein, the sulfurtransferase 18 (STR18), and a CD isoform referred to as ABA3, and compare these biochemical features to those of a natural CD-Rhd fusion protein from the bacterium Pseudorhodoferax sp.. We observed that the bacterial enzyme is bifunctional exhibiting both CD and STR activities using L-cysteine and thiosulfate as sulfur donors but preferentially uses L-cysteine to catalyze trans-persulfidation reactions. In vitro activity assays and mass spectrometry analyses revealed that STR18 stimulates the CD activity of ABA3 by reducing the intermediate persulfide on its catalytic cysteine thereby accelerating the overall transfer reaction. Both proteins interact in planta and form an efficient sulfur relay system whereby STR18 catalyzes trans-persulfidation reactions from ABA3 to the model acceptor protein roGFP2. In conclusion, the ABA3-STR18 couple likely represents an uncharacterized pathway of sulfur trafficking in the cytosol of plant cells, independent of ABA3 function in molybdenum cofactor maturation.


title: Age-dependent selection on MHC class 2 variation in a free-living ruminant

abstract: Genes within the major histocompatibility complex (MHC) are the most variable identified in vertebrates. Pathogen-mediated selection (PMS) is believed to be the main force maintaining diversity at MHC class I and II genes, but it has proven hard to demonstrate the exact PMS regime that is acting in natural populations. Demonstrating contemporary selection on MHC alleles is not trivial, and previous work has been constrained by limited genetic tools, low sample sizes and short time scales and has sometimes involved anticonservative statistical approaches. Here, we use appropriate statistical approaches to examine associations between MHC variation and several fitness measurements including total fitness (lifetime breeding success) and five fitness components, in 3400 wild Soay sheep (Ovis aries) monitored over their lifetimes between 1989 and 2012. We found haplotypes C and D were associated with decreased and increased male total fitness respectively. In terms of fitness components, juvenile survival was positively associated with haplotype divergence. Of the eight MHC haplotypes (A-H), haplotypes C and F were associated with decreased adult male breeding success and decreased adult female life span respectively. Consistent with the increased male total fitness, haplotype D, which is the rarest, has increased in frequency throughout the study period. Our results suggest that contemporary balancing selection is acting on MHC class II genes in Soay sheep and that different selection mechanisms are acting between juveniles and adults.


title: Non-Parametric Analysis of Inter-Individual Relations Using an Attention-Based Neural Network

abstract:
O_LISocial network analysis, which has been widely adopted in animal studies over the past decade, enables the revelation of global characteristic patterns of animal social systems from pairwise inter-individual relations. Animal social networks are typically drawn based on geometric proximity and/or frequency of social behaviors (e.g., grooming), but the appropriate metric for inter-individual relationship is not clear, especially when prior knowledge on the species/data is limited.
C_LIO_LIIn this study, researchers explored a non-parametric analysis of inter-individual relations using a neural network with the attention mechanism, which plays a central role in natural language processing. The high interpretability of the attention mechanism and flexibility of the entire neural network allow for automatic detection of inter-individual relations included in the raw data, without requiring prior knowledge/assumptions about what modes/types of relations are included in the data. For these case studies, three-dimensional location data collected from simulated agents and real Japanese macaques were analyzed.
C_LIO_LIThe proposed method successfully recovered the latent relations behind the simulated data and discovered female-oriented relations in the real data, which are in accordance with previous generalizations about the macaque social structure.
C_LIO_LIThe proposed method does not exploit any behavioral patterns that are particular to Japanese macaques, and researchers can use it for location data of other animals. The exibility of the neural network would also allow for its application to a wide variety of data with interacting components, such as vocal communication.
C_LI

category: animal behavior and cognition


title: CD19-targeting CAR T cells potently redirected to kill solid tumor cells

abstract: Successful CAR T cell therapy for the treatment of solid tumors requires exemplary CAR T cell expansion, persistence and fitness, and the ability to target tumor antigens safely. Here we address this constellation of critical attributes for successful cellular therapy by using integrated technologies that simplify development and derisk clinical translation. We have developed a CAR-CD19 T cell that secretes a CD19-anti-Her2 bridging protein. This cell therapy strategy exploits the ability of CD19-targeting CAR T cells to interact with CD19 on normal B cells to drive expansion, persistence and fitness. The secreted bridging protein potently binds to Her2-positive tumor cells, mediating CAR-CD19 T cell cytotoxicity in vitro and in vivo. Because of its short half-life, the secreted bridging protein will selectively accumulate at the site of highest antigen expression, ie. at the tumor. Bridging proteins that bind to multiple different tumor antigens have been created. Therefore, antigen-bridging CAR-CD19 T cells incorporate critical attributes for successful solid tumor cell therapy. This platform can be exploited to attack tumor antigens on any cancer.

category: cancer biology


title: A natural timeless polymorphism allowing circadian clock synchronization in 'white nights'

abstract: Daily temporal organisation of behavioural and physiological functions offers a fitness advantage for most animals. Optimized temporal niches are determined by an interplay between external environmental rhythms and internal circadian clocks. While daily light:dark cycles serve as a robust time cue (Zeitgeber) to synchronise circadian clocks, it is not clear how animals experiencing only weak environmental cues deal with this problem. Like humans, flies of the genus Drosophila originate in sub-Saharan Africa and spread North in Europe up to the polar circle where they experience extremely long days in the summer or even constant light (LL). LL is known to disrupt clock function, due to constant activation of the deep brain photoreceptor CRYPTOCHROME (CRY), which induces constant degradation of the clock protein TIMELESS (TIM). Temperature cycles are able to overcome these arrhythmia inducing effects of LL, reinstating clock protein oscillations and rhythmic behaviour. We show here that for this to occur a recently evolved natural allele (ls-tim) of the timeless gene is required, whereby the presence of this allele within the central clock neurons is sufficient. The ls-tim allele encodes a longer, less-light sensitive form of TIM (L-TIM) in addition to the shorter (S-TIM) form, the only form encoded by the ancient s-tim allele. Only after blocking light-input by removing functional CRY, s-tim flies are able to synchronise molecular and behavioural rhythms to temperature cycles in LL. Additional removal of light input from the visual system results in a phase advance of molecular and behavioural rhythms, showing that the visual system contributes to temperature synchronization in LL. We show that ls-tim, but not s-tim flies can synchronise their behavioural activity to semi-natural LL and temperature cycle conditions reflecting long Northern Europe summer days, the season when Drosophila populations massively expand. Our observations suggest that this functional gain associated with ls-tim is driving the Northern spread of this allele by directional selection.


title: Hierarchical progressive learning of cell identities in single-cell data

abstract: Supervised methods are increasingly used to identify cell populations in single-cell data. Yet, current methods are limited in their ability to learn from multiple datasets simultaneously, are hampered by the annotation of datasets at different resolutions, and do not preserve annotations when retrained on new datasets. The latter point is especially important as researchers cannot rely on downstream analysis performed using earlier versions of the dataset. Here, we present scHPL, a hierarchical progressive learning method which allows continuous learning from single-cell data by leveraging the different resolutions of annotations across multiple datasets to learn and continuously update a classification tree. We evaluate the classification and tree learning performance using simulated as well as real datasets and show that scHPL can successfully learn known cellular hierarchies from multiple datasets while preserving the original annotations. scHPL is available at https://github.com/lcmmichielsen/hierarchicalprogressivelearning.

category: bioinformatics


title: DAISM-DNNXMBD: Highly accurate cell type proportion estimation with in silico data augmentation and deep neural networks

abstract: Understanding the immune cell abundance of cancer and other disease-related tissues has an important role in guiding disease treatments. Computational cell type proportion estimation methods have been previously developed to derive such information from bulk RNA sequencing (RNA-seq) data. Unfortunately, our results show that the performance of these methods can be seriously plagued by the mismatch between training data and real-world data. To tackle this issue, we propose the DAISM-DNNXMBD1 pipeline that trains a deep neural network (DNN) with dataset-specific training data populated from a small number of calibrated samples using DAISM, a novel Data Augmentation method with an In Silico Mixing strategy. The evaluation results demonstrate that the DAISM-DNN pipeline outperforms other existing methods consistently and substantially for all the cell types under evaluation on real-world datasets.


title: Tracking bacterial lineages in complex and dynamic environments with applications to growth control and persistence

abstract: As bacteria transition from exponential to stationary phase they change greatly in size, morphology, growth and expression-profiles. These responses also vary between individual cells, but it has proven difficult to track cell lineages along the growth curve to determine the progression of events or correlations between how individual cells enter and exit dormancy. We developed a platform for tracking >105 parallel cell lineages in dense and changing cultures, independently validating that the imaged cells closely track batch populations. Initial applications show that for both Escherichia coli and Bacillus subtilis, growth changes from an  adder mode in exponential phase to mixed  adders-timers entering stationary phase, and then a near-perfect  sizer upon exit - creating broadly distributed cell sizes in stationary phase but rapidly returning to narrowly distributed sizes upon exit. Furthermore, cells that undergo more divisions entering stationary phase suffer reduced survival after long periods of dormancy but are the only cells observed that persist against antibiotics.


title: A novel platform of RNA 2-O-methylation high-throughput and site-specific quantification tools revealed its broad distribution

abstract: Ribose 2'-O-methylation is involved in critical biological processes, but its biological functions and significance remain largely unknown due to the lack of accurate and efficient identification tools. To overcome this gap, we established an accurate high-throughput single-base Nm site screening method, NJU-seq (Nm-Judged-Universally sequencing), and a site-specific 2'-O-methylation quantification tool, Nm-VAQ (RNase H-based Nm Validation and Absolute Quantification). In this study, we identified several new Nm sites on HeLa rRNA and thousands of new Nm sites on mRNA of human and mouse cell lines, and 13 of them were selected and further validated. The 2'-O-methylation ratios of all validated Nm sites on mRNA were lower than 20%, which was out of the detection abilities of previous methods. Our study suggested that the function and regulation of Nm sites are more complicated than previous studies revealed and worth further mechanistic investigation.


title: Evaluating the transcriptional fidelity of cancer models

abstract:
BackgroundCancer researchers use cell lines, patient derived xenografts, engineered mice, and tumoroids as models to investigate tumor biology and to identify therapies. The generalizability and power of a model derives from the fidelity with which it represents the tumor type under investigation, however, the extent to which this is true is often unclear. The preponderance of models and the ability to readily generate new ones has created a demand for tools that can measure the extent and ways in which cancer models resemble or diverge from native tumors.

MethodsWe developed a machine learning based computational tool, CancerCellNet, that measures the similarity of cancer models to 22 naturally occurring tumor types and 36 subtypes, in a platform and species agnostic manner. We applied this tool to 657 cancer cell lines, 415 patient derived xenografts, 26 distinct genetically engineered mouse models, and 131 tumoroids. We validated CancerCellNet by application to independent data, and we tested several predictions with immunofluorescence.

ResultsWe have documented the cancer models with the greatest transcriptional fidelity to natural tumors, we have identified cancers underserved by adequate models, and we have found models with annotations that do not match their classification. By comparing models across modalities, we report that, on average, genetically engineered mice and tumoroids have higher transcriptional fidelity than patient derived xenografts and cell lines in four out of five tumor types. However, several patient derived xenografts and tumoroids have classification scores that are on par with native tumors, highlighting both their potential as faithful model classes and their heterogeneity.

ConclusionsCancerCellNet enables the rapid assessment of transcriptional fidelity of tumor models. We have made CancerCellNet available as freely downloadable software and as a web application that can be applied to new cancer models that allows for direct comparison to the cancer models evaluated here.


title: Flexible multivariate linear mixed models for structured multiple traits

abstract: Many genetic studies collect structured multivariate traits containing rich information across traits. We present a flexible multivariate linear mixed model for quantitative trait loci mapping (FlxQTL) for multiple correlated traits that adjusts for genetic relatedness and that models information on multiple environments or multiple timepoints using trait covariates. FlxQTL handles genetic mapping of multivariate traits faster with greater flexibility compared to previous implementations.


title: Frisky CALF sometimes outruns LASSO

abstract:
Regularized regression analysis is a mature analytic approach to identify weighted sums of variables that predict outcomes. Typically, the number of subjects (N) is smaller than the number of predictors (p). Here, we present a novel coarse approximation linear function (CALF) to frugally select important predictors, build linear models, and discover causality.

CALF is a linear regression strategy applied to normalized data that employs only a few (such as 2 to 20) nonzero weights, each +1 or -1. Metrics can be Welch t-test p-value, area under curve (AUC) of receiver operating characteristic, or Pearson correlation, depending upon data type and user preferences. For quantitative approximations, a linear fit (adding an intercept value and rescaling {+/-}1 weights by a common multiplier) can be added to optimize mean squared error (MSE).

Real medical data of five types were used to generate examples with goals of binary classification or real variable approximation. Predictors considered were real, sets of integers, or ternary values of single nucleotide polymorphisms. When applied to real data, CALF approximations outperformed in p-value, AUC, correlation, or MSE a popular regularized linear regression algorithm, namely, basic LASSO. It appears that using LASSO without considering CALF might risk wasting resources.

AvailabilityR version: Comprehensive R Archive (CRAN): https://cran.r-project.org/web/packages/CALF/index.html

Python 3.x version: GitHub: https://github.com/jorufo/CALF_Python

Contactclark_jeffries@med.unc.edu


title: The Hoja loca1 Mutant and AUX/IAA Function in Phytomer Development

abstract: Plant architecture is determined by the iterative building of phytomer units, each containing a portion of stem, an organ and an axillary meristem. Each phytomer can follow different developmental paths underpinning the complexity, and plasticity, of plant form. Auxin plays a central role in the coordination of phytomer development, regulating organ initiation and patterning across all axes. This diversity in auxin function results from changes in the activities and expression of auxin signaling components, including the AUX/IAA repressors. Higher land plants have multigene AUX/IAA families, which leads to functional redundancy. Dominant mutations, which prevent AUX/IAA degradation in response to auxin, have highlighted the importance of these proteins in fine-tuning organ development. Here we report a new dominant AUX/IAA mutant in maize, Hoja loca1 (Oja). Oja has a mutation in the degron motif of ZmIAA28 and specifically affects aerial organ initiation and leaf medio-lateral patterning, but phytomer initiation remains unchanged. Mutant phenotypes are variable and transcriptional profiling of individual meristems identified clusters of genes that may underpin the phenotypic differences. The unique phenotype of Oja provides evidence of species-specific sub-functionalization of the AUX/IAAs, and illustrates the crucial role of auxin signaling in the tight coordination of phytomer unit development.


title: DeepLC can predict retention times for peptides that carry as-yet unseen modifications

abstract: The inclusion of peptide retention time prediction promises to remove peptide identification ambiguity in complex LC-MS identification workflows. However, due to the way peptides are encoded in current prediction models, accurate retention times cannot be predicted for modified peptides. This is especially problematic for fledgling open modification searches, which will benefit from accurate retention time prediction for modified peptides to reduce identification ambiguity. We here therefore present DeepLC, a novel deep learning peptide retention time predictor utilizing a new peptide encoding based on atomic composition that allows the retention time of (previously unseen) modified peptides to be predicted accurately. We show that DeepLC performs similarly to current state-of-the-art approaches for unmodified peptides, and, more importantly, accurately predicts retention times for modifications not seen during training. Moreover, we show that DeepLCs ability to predict retention times for any modification enables potentially incorrect identifications to be flagged in an open modification search of CD8-positive T-cell proteome data. DeepLC is available under the permissive Apache 2.0 open source license and comes with a user-friendly graphical user interface, as well as a Python package on PyPI, Bioconda, and BioContainers for effortless workflow integration.


title: Evaluating and integrating spatial capture-recapture models with data of variable individual identifiability

abstract: Spatial capture-recapture (SCR) models have become the preferred tool for estimating densities of carnivores. Within this family of models are variants requiring identification of all individuals in each encounter (SCR), a subset of individuals only (generalized spatial mark-resight, gSMR), or no individual identification (spatial count or spatial presence-absence). Although each technique has been shown through simulation to yield unbiased results, the consistency and relative precision of estimates across methods in real-world settings are seldom considered. We tested a suite of models ranging from those only requiring detections of unmarked individuals to others that integrate remote camera, physical capture, genetic, and global positioning system (GPS) data into a  hybrid model, to estimate population densities of black bears, bobcats, cougars, and coyotes. For each species we genotyped fecal DNA collected with detection dogs during a 20-day period. A subset of individuals from each species was affixed with GPS collars bearing unique markings and resighted by remote cameras over 140 days contemporaneous with scat collection. Camera-based gSMR models produced density estimates that differed by less than 10% from genetic SCR for bears, cougars, and coyotes once important sources of variation (sex or behavioral status) were controlled for. For bobcats, SCR estimates were 33% higher than gSMR. The cause of the discrepancies in estimates was likely attributable to challenges designing a study compatible for species with disparate home range sizes and the difficulty of collecting sufficient data in a timeframe in which demographic closure could be assumed. Unmarked models estimated densities that varied greatly from SCR, but estimates became more consistent in models wherein more individuals were identifiable. Hybrid models containing all data sources exhibited the most precise estimates for all species. For studies in which only sparse data can be obtained and the strictest model assumptions are unlikely to be met, we suggest researchers use caution making inference from models lacking individual identity. For best results, we further recommend the use of methods requiring at least a subset of the population is marked and that multiple datasets are incorporated when possible.


title: Homeostatic control of synaptic rewiring in recurrent networks induces the formation of stable memory engrams

abstract: Brain networks store new memories using functional and structural synaptic plasticity. Memory formation is generally attributed to Hebbian plasticity, while homeostatic plasticity is thought to have an ancillary role in stabilizing network dynamics. Here we report that homeostatic plasticity alone can also lead to the formation of stable memories. We analyze this phenomenon using a new theory of network remodeling, combined with numerical simulations of recurrent spiking neural networks that exhibit structural plasticity based on firing rate homeostasis. These networks are able to store repeatedly presented patterns and recall them upon the presentation of incomplete cues. Storage is fast, governed by the homeostatic drift. In contrast, forgetting is slow, driven by a diffusion process. Joint stimulation of neurons induces the growth of associative connections between them, leading to the formation of memory engrams. These memories are stored in a distributed fashion throughout connectivity matrix, and individual synaptic connections have only a small influence. Although memory-specific connections are increased in number, the total number of inputs and outputs of neurons undergo only small changes during stimulation. We find that homeostatic structural plasticity induces a specific type of "silent memories", different from conventional attractor states.


title: Reconstruction of microbial haplotypes by integration of statistical and physical linkage in scaffolding

abstract: DNA sequencing technologies provide unprecedented opportunities to analyze within-host evolution of microorganism populations. Often, within-host populations are analyzed via pooled sequencing of the population, which contains multiple individuals or  haplotypes. However, current next-generation sequencing instruments, in conjunction with single-molecule barcoded linked-reads, cannot distinguish long haplotypes directly. Computational reconstruction of haplotypes from pooled sequencing has been attempted in virology, bacterial genomics, metagenomics and human genetics, using algorithms based on either cross-host genetic sharing or within-host genomic reads. Here we describe PoolHapX, a flexible computational approach that integrates information from both genetic sharing and genomic sequencing. We demonstrated that PoolHapX outperforms state-of-the-art tools tailored to specific organismal systems, and is robust to within-host evolution. Importantly, together with barcoded linked-reads, PoolHapX can infer whole-chromosome-scale haplotypes from 50 pools each containing 12 different haplotypes. By analyzing real data, we uncovered dynamic variations in the evolutionary processes of within-patient HIV populations previously unobserved in single position-based analysis.


title: Multimodal brain imaging study of 36,678 participants reveals adverse effects of moderate drinking

abstract:
Heavy alcohol consumption can have significant deleterious neural consequences, including brain atrophy, neuronal loss, poorer white matter fiber integrity, and cognitive decline. However, the effects of light-to-moderate alcohol consumption on brain structure remain unclear. Here, we examine the associations between alcohol intake and brain structure using multimodal imaging data from 36,678 generally healthy middle-aged and older adults from the UK Biobank, controlling for numerous potential confounds. We find negative associations between alcohol intake and global gray matter volume (GMV) and white matter volume (WMV), which become stronger as intake increases. An examination of the associations between alcohol intake and 139 regional GMV imaging-derived phenotypes (IDPs) and 375 WM microstructure IDPs yielded 304 (59.1%) significant findings, including 125 GMV IDPs that are spread across the brain and 179 WM microstructure IDPs across multiple tract regions. In general, findings comport with the existing literature. However, a daily alcohol intake of as little as one to two units - 250 to 500 ml of a 4% beer or 76 to 146 ml of a 13% wine - is already associated with GMV deficits and altered WMV microstructure, placing moderate drinkers at risk.

One Sentence SummaryModerate alcohol intake, consuming one or more daily alcohol units, has adverse effects on brain health.


title: Distinct ensembles in the noradrenergic locus coeruleus evoke diverse cortical states

abstract: The noradrenergic locus coeruleus (LC) is a crucial controller of brain and behavioral states. Activating LC neurons synchronously en masse by electrical or optogenetic stimulation promotes a stereotypical "activated" high-frequency cortical state. However, it has been recently reported that spontaneous LC cell-pairs have sparse yet structured time-averaged cross-correlations, which is unlike the high synchrony of en masse neuronal stimulation. This suggests the untested possibility that LC population activity may be made of distinct multi-cell ensembles each with unique temporal evolution of activity. We used non-negative matrix factorization (NMF) to analyze large populations of LC single units simultaneously recorded in the rat LC. Synthetic spike train simulations showed that NMF, unlike the traditional time-averaged pairwise correlations, detects both the precise neuronal composition and the activation time courses of each ensemble. NMF identified the existence of robust ensembles of spontaneously co-active LC neurons. Since LC neurons selectively project to specific forebrain regions, we hypothesized that individual LC ensembles produce different cortical states. To test this hypothesis, we triggered local field potentials (LFP) in cortical area 24a on the activation of distinct LC ensembles. We found four cortical states, each with different spectro-temporal LFP characteristics, that were robust across sessions and animals. While some LC ensembles triggered the activated state, others were associated with a beta oscillation-specific state or a reduced high frequency oscillation state. Thus - in contrast to the stereotypical "activated" brain state evoked by en masse LC stimulation - spontaneous activation of distinct LC ensembles can control a multitude of cortical states.


title: Highly Versatile, Non-Invasive Method for Collecting Buccal DNA from Free-Ranging Non-Human Primates

abstract: Non-invasive techniques for collection of DNA samples of suitable quality and quantity are important for improving the efficiency of genetic wildlife research. The development of a non-invasive method for collection of DNA samples from wild stump-tailed macaques (Macaca arctoides) is described herein. Sterilized polyester rope was cut into 10 cm pieces, which were then soaked in a 20% sugar solution to bait individuals. Rope swabs were immediately collected and transferred to a lysis buffer solution after subjects had picked up, chewed, and discarded them. DNA was later extracted from the buffer. Quantitative real-time PCR and both allelic dropout and genotype failure rates were used to compare the quantity and quality of the buccal DNA samples to those of intestinal slough cell DNA samples collected from freshly dropped feces. The buccal samples yielded significantly more DNA (27.1 {+/-} 33.8 ng/{micro}L) than did the fecal samples (11.4 {+/-} 15.4 ng/{micro}L) and exhibited lower allelic dropout and genotyping failure rates for the 10 autosomal microsatellites investigated. Buccal cell collection was also simple, inexpensive, reliable, and less time-consuming compared to fecal sampling. Thus, this method should facilitate genome-wide studies of non-human primates and other wildlife species.


title: Dopamine Signaling in the Dorsomedial Striatum Promotes Compulsive Behavior

abstract: Compulsive behavior is a defining feature of disorders such as substance use disorder and obsessive-compulsive disorder. Current evidence suggests that corticostriatal circuits control the expression of established compulsions, but little is known about the mechanisms regulating the development of compulsions. We hypothesized that dopamine, a critical modulator of striatal synaptic plasticity, could control alterations in corticostriatal circuits leading to the development of compulsions (defined as continued reward-seeking in the face of punishment). We used dual-site fiber photometry to measure dopamine axon activity in the dorsomedial striatum (DMS) and the dorsolateral striatum (DLS) as compulsions emerged. Individual variability in the speed with which compulsions emerged was predicted by DMS dopamine axon activity. Amplifying this dopamine signal accelerated animals transitions to compulsion, whereas inhibition led to learning delays. In contrast, amplifying DLS dopamine signaling had no effect on the emergence of compulsions. These results establish DMS dopamine signaling as a key controller of the development of compulsive reward-seeking.


title: Fragmentation of Small-cell Lung Cancer Regulatory States in Heterotypic Microenvironments

abstract:
Small-cell lung cancers derive from pulmonary neuroendocrine cells, which have stemlike properties to reprogram into other cell types upon lung injury. It is difficult to uncouple transcriptional plasticity of these transformed cells from genetic changes that evolve in primary tumors or secondary metastases. Profiling of single cells also problematic if the required sample dissociation activates injury-like signaling and reprogramming. Here, we defined cell-state heterogeneities in situ through laser capture microdissection-based 10-cell transcriptomics coupled with stochastic-profiling fluctuation analysis. Using labeled cells from a small-cell lung cancer mouse model initiated by neuroendocrine deletion of Rb1-Trp53, we profiled variations in transcript abundance to identify cell-to-cell differences in regulatory state in vitro and in vivo. Fluctuating transcripts in spheroid culture were partly shared among Rb1-Trp53-null models, and heterogeneities increased considerably when cells were delivered intravenously to colonize the liver. Colonization of immunocompromised animals drove fractional appearance of alveolar type II-like markers and poised cells for paracrine stimulation from immune cells and hepatocytes. Immunocompetency further exaggerated the fragmentation of tumor states in the liver, yielding mixed stromal signatures evident in bulk sequencing from autochthonous tumors and metastases. Dozens of transcript heterogeneities recur irrespective of biological context; their mapped orthologs brought together observations of murine and human small-cell lung cancer. Candidate heterogeneities recurrent in the liver also stratified primary human tumors into discrete groups not readily explained by molecular subtype but with prognostic relevance. We conclude that heterotypic interactions in the liver and lung are an accelerant for intratumor heterogeneity in small-cell lung cancer.

Statement of significanceThe single-cell regulatory heterogeneity of small-cell lung cancer becomes increasingly elaborate in the liver, a common metastatic site for the disease.


title: Premalignant Nf1,Trp53-null Oligodendrocyte Precursor Cells BecomeStalled in a Heterogeneous State of Replication Stress BeforeGliomagenesis

abstract:
Cancer evolves from premalignant clones that accumulate mutations and adopt unusual cell states to achieve transformation. Previously, we pinpointed the oligodendrocyte precursor cell (OPC) as a cell-of-origin for glioma, but the early changes of mutant OPCs during premalignancy remained unknown. Using mice engineered for inducible Nf1-Trp53 loss in OPCs, we acutely isolated labeled mutant OPCs by laser-capture microdissection and determined gene-expression changes by bulk RNA sequencing and a fluctuation analysis, called stochastic profiling, which uses RNA-sequencing measurements from random pools of 10 mutant cells. At 12 days after Nf1-Trp53 deletion, while bulk differences were mostly limited to mitotic hallmarks and genes for ribosome biosynthesis, stochastic profiling revealed a spectrum of stem-progenitor (Axl, Aldh1a1), proneural, and mesenchymal states as potential starting points for gliomagenesis. At 90 days, bulk sequencing detected very few differentially expressed transcripts, whereas stochastic profiling revealed cell states for neurons and mural cells that do not give rise to glial tumors, suggesting cellular dead-ends for gliomagenesis. Importantly, we identified mutant OPCs that strongly expressed key effectors of nonsense-mediated decay (Upf3b) and homology-dependent DNA repair (Rad51c, Slx1b, Ercc4) along with DNA-damage markers suggesting transcription-associated replication stress. Analysis of 10-cell transcriptomes at 90 days identified a locus of elevated gene expression containing an additional repair endonuclease (Mus81) and Rin1, a Ras-Raf antagonist and possible counterbalance to Nf1 loss. At 150 days, Rin1 was microdeleted in some gliomas and downregulated in all others. Replication stress may pose a considerable bottleneck that must be resolved for gliomas to initiate.

Statement of significanceIn situ stochastic profiling of heterogeneous cell states in a mouse model of glioma uncovers regulatory confusion in a glioma cell-of-origin and defines a state of replication stress that precedes tumor initiation.


title: Phylogenetic, geographic and ecological distribution of a green-brown polymorphisms in European Orthopterans

abstract: The green-brown polymorphism among polyneopteran insects represents one of the most penetrant color polymorphisms in any group of organisms. Yet systematic overviews are lacking. I here present analyses of the phylogenetic, geographic and habitat distribution of the green-brown polymorphism across the complete European orthopteran fauna. Overall, 30% of European orthopterans are green-brown polymorphic. Polymorphic species are scattered across the entire phylogenetic tree, including roughly equal proportions of Ensifera and Caelifera. A few taxonomic groups, however, include only brown species. Polymorphic species occur more frequently in clades that contain monomorphic green species than in those without green species. The relative abundance of color morphs in polymorphic species is skewed towards green, and in particular rare/exceptional brown morphs are more common in predominantly green species than rare/exceptional green morphs in predominantly brown species. The patterns of abundances support the hypothesis that loss-of-function mutations play a role in creating polymorphic populations from green species. Polymorphic species are particularly common in moist to mesic grasslands, alpine and arboreal habitats. Dry, open, rocky and cave habitats as well as nocturnal lifestyles are dominated by monomorphic brown species. The proportion of polymorphic species increases from southern to northern latitudes. These marked habitat-dependencies also show that coloration is affected by natural selection and/or environmental filtering. Overall, the results illustrate that the occurrence of the polymorphism is phylogenetically, geographically and ecologically widespread and they suggest that polymorphism is thus potentially in mutation-selection balance across a large number of species.


title: A Mendelian randomization study of glycemic and anthropometric traits and Parkinson's disease

abstract:
BackgroundImpaired glucose and obesity are frequently observed in patients with Parkinsons disease (PD), although it is unclear whether the impairment precedes or results from the neurodegeneration.

ObjectiveWe aimed to assess whether glycemic and anthropometric traits can influence the risk of PD in 33,674 cases and 449,056 healthy controls using the Mendelian randomization (MR) framework.

MethodsWe investigated causality with a two-sample MR approach in the European population to compute effect estimates with summary statistics from available discovery meta-analyses of genome-wide association studies (GWAS) on glycemic and anthropometric traits.

ResultsWe considered a threshold of p-value=0.0038 as significant after accounting for multiple testing, and p-value<0.05 was considered to be a suggestive evidence for a potential association. We observed a protective effect of waist-hip ratio (WHR) on PD (Inverse variance-weighted (IVW): OR IVW=0.735; 95%CI= 0.622-0.868; p-value=0.0003; I2 index=22.0%; MR-Egger intercept p-value=0.1508; Cochran Q test p-value=0.0003). The association was further retained after the exclusion of overlapping UK biobank (UKB) samples between the WHR and PD datasets (ORIVW=0.791; 95%CI=0.659-0.950; p-value=0.012; I2 index=13.0%; MR-Egger intercept p-value=0.733; Cochran Q test p-value=0.035). The sensitivity analysis provided suggestive evidence of an increased risk of PD on fasting glucose (FG) ({beta} IVW=0.0188; 95%CI=0.0062-0.0313, p-value=0.0055; I2 index=0.0%; MR-Egger intercept p-value=0.0957; Cochran Q test p-value=0.4555) and protective effect of PD on T2D (Weighted median effect: ORWME=0.946; 95%CI=0.9290.983; p-value=0.0051; Weighted mode effect: ORMBE=0.943; 95%CI=0.904-0.983; p-value=0.0116).

ConclusionsOur results showed that central or abdominal obesity may be protective against PD development, independent of glucose levels.


title: Highly efficient multiplex editing: One-shot generation of 8x Nicotiana benthamiana and 12x Arabidopsis mutants

abstract: Genome editing by RNA-guided nucleases, such as SpCas9, has been used in numerous different plant species. However, to what extent multiple independent loci can be targeted simultaneously by multiplexing has not been well-documented. Here, we developed a toolkit, based on a highly intron-optimized zCas9i gene, which allows assembly of nuclease constructs expressing up to 32 sgRNAs. We used this toolkit to explore the limits of multiplexing in two major model species, and report on isolation of transgene-free octuple Nicotiana benthamiana and duodecuple (12x) Arabidopsis thaliana mutant lines in a single generation (T1 and T2, respectively). We developed novel counter-selection markers for N. benthamiana, most importantly Sl-FAST2, comparable to the well-established Arabidopsis seed fluorescence marker, and FCY-UPP, based on production of toxic 5-fluorouracyl in presence of a precursor. Targeting eight genes with an array of nine different sgRNAs and relying on FCY-UPP for selection of non-transgenic T1, we identified N. benthamiana mutant lines with astonishingly high efficiencies: All analyzed plants carried mutations in all genes (~112/116 target sites edited). Furthermore, we targeted 12 genes by an array of 24 sgRNAs in A. thaliana. Efficiency was significantly lower in A. thaliana, and our results indicate Cas9 availability is the limiting factor in such higher order multiplexing applications. We identify a duodecuple mutant line by a combination of phenotypic screening and amplicon sequencing. The resources and results presented provide new perspectives for how multiplexing can be used to generate complex genotypes or to functionally interrogate groups of candidate genes.


title: The quorum-sensing systems of Vibrio campbellii DS40M4 and BB120 are genetically and functionally distinct

abstract:
Vibrio campbellii BB120 (previously classified as Vibrio harveyi) is a fundamental model strain for studying quorum sensing in vibrios. A phylogenetic evaluation of sequenced Vibrio strains in Genbank revealed that BB120 is closely related to the environmental isolate V. campbellii DS40M4. We exploited DS40M4s competence for exogenous DNA uptake to rapidly generate >30 isogenic strains with deletions of genes encoding BB120 quorum-sensing system homologs. Our results show that the quorum-sensing circuit of DS40M4 is distinct from BB120 in three ways: 1) DS40M4 does not produce an acyl homoserine lactone autoinducer but encodes an active orphan LuxN receptor, 2) the quorum regulatory small RNAs (Qrrs) are not solely regulated by autoinducer signaling through the response regulator LuxO, and 3) the DS40M4 quorum-sensing regulon is much smaller than BB120 (~100 genes vs ~400 genes, respectively). Using comparative genomics to expand our understanding of quorum-sensing circuit diversity, we observe that conservation of LuxM/LuxN proteins differs widely both between and within Vibrio species. These strains are also phenotypically distinct: DS40M4 exhibits stronger interbacterial cell killing, whereas BB120 forms more robust biofilms and is bioluminescent. These results underscore the need to examine wild isolates for a broader view of bacterial diversity in the marine ecosystem.

Originality-Significance StatementWild bacterial isolates yield important information about traits that vary within species. Here, we compare environmental isolate Vibrio campbellii DS40M4 to its close relative, the model strain BB120 that has been a fundamental strain for studying quorum sensing for >30 years. We examine several phenotypes that define this species, including quorum sensing, bioluminescence, and biofilm formation. Importantly, DS40M4 is naturally transformable with exogenous DNA, which allows for the rapid generation of mutants in a laboratory setting. By exploiting natural transformation, we genetically dissected the functions of BB120 quorum-sensing system homologs in the DS40M4 strain, including two-component signaling systems, transcriptional regulators, and small RNAs.


title: Persister state-directed transitioning and vulnerability in melanoma

abstract: Melanoma is a highly plastic tumor characterized by dynamic interconversion of different cell identities depending on the biological context. For example, melanoma cells with high expression of the H3K4 demethylase KDM5B (JARID1B) rest in a slow-cycling, yet reversible persister state. Over time, KDM5Bhigh cells can promote rapid tumor repopulation with equilibrated KDM5B expression heterogeneity. The cellular identity of KDM5Bhigh persister cells has not been studied so far, missing an important cell state-directed treatment opportunity in melanoma. Here, we have established a doxycycline-titratable system for genetic induction of permanent intratumor expression of KDM5B and screened for chemical agents that phenocopy this effect. Transcriptional profiling and cell functional assays confirmed that the dihydropyridine phenoxyethyl 4-(2-fluorophenyl)-2,7,7-trimethyl-5-oxo-1,4,5,6,7,8-hexa-hydro-quinoline-3-carboxylate (termed Cpd1) supports high KDM5B expression and directs melanoma cells towards differentiation along the melanocytic lineage and to cell cycle-arrest. The high KDM5B state additionally prevents cell proliferation through negative regulation of cytokinetic abscission. Moreover, treatment with Cpd1 promoted the expression of the melanocyte-specific tyrosinase gene specifically sensitizing melanoma cells for the tyrosinase-processed antifolate prodrug 3-O-(3,4,5-trimethoxybenzoyl)-(-)-epicatechin (TMECG). In summary, our study provides proof-of-concept for a new dual hit strategy in melanoma, in which persister state-directed transitioning limits tumor growth and plasticity and primes melanoma cells towards lineage-specific elimination.


title: Mapping shifts in nanopore signal to changes in protein and protein-DNA conformation

abstract:
Solid-state nanopores have been used extensively in biomolecular studies involving DNA and proteins. However, the interpretation of signals generated by the translocation of proteins or protein-DNA complexes remains challenging. Here, we investigate the behavior of monovalent streptavidin and the complex it forms with short biotinylated DNA over a range of nanopore sizes, salts and voltages. We describe a simple geometric model that is broadly applicable and employ it to explain observed variations in conductance blockage and dwell time with experimental conditions. The general approach developed here underscores the value of nanopore-based protein analysis and represents progress toward the interpretation of complex translocation signals.

STATEMENT OF SIGNIFICANCENanopore sensing allows investigation of biomolecular structure in aqueous solution, including electricfield-induced changes in protein conformation. This nanopore-based study probes: (1) the tetramerdimer transition of streptavidin, observing the effects of increasing voltage with varying salt type and concentration; (2) the possible conformational states of DNA-streptavidin complexes when confined inside a pore. We describe a broadly applicable geometric approach that maps stepwise changes in the nanopore signal to real-time conformational transitions. These results represent progress toward accurate interpretation of nanopore signals generated by molecular complexes.


title: Pan-cancer Drivers are Recurrent Transcriptional Regulatory Heterogeneities in Early-stage Luminal Breast Cancer

abstract:
The heterogeneous composition of solid tumors is known to impact disease progression and response to therapy. Malignant cells coexist in different regulatory states that can be accessed transcriptomically by single-cell RNA sequencing, but these methods have many caveats related to sensitivity, noise, and sample handling. We revised a statistical fluctuation analysis called stochastic profiling to combine with 10-cell RNA sequencing, which was designed for laser-capture microdissection (LCM) and extended here for immuno-LCM. When applied to a cohort of late-onset, early-stage luminal breast cancers, the integrated approach identified thousands of candidate regulatory heterogeneities. Intersecting the candidates from different tumors yielded a relatively stable set of 710 recurrent heterogeneously expressed genes (RHEGs) that were significantly variable in >50% of patients. RHEGs were not strongly confounded by dissociation artifacts, cell cycle oscillations, or driving mutations for breast cancer. Rather, we detected RHEG enrichments for epithelial-to-mesenchymal transition genes and, unexpectedly, the latest pan-cancer assembly of driver genes across cancer types other than breast. Heterogeneous transcriptional regulation conceivably provides a faster, reversible mechanism for malignant cells to evaluate the effects of potential oncogenes or tumor suppressors on cancer hallmarks.

Statement of significanceProfiling intratumor heterogeneity of luminal breast carcinoma cells identifies a recurrent set of genes suggesting sporadic activation of pathways known to drive other types of cancer.


title: Cell-type specialization in the brain is encoded by specific long-range chromatin topologies

abstract: Neurons and oligodendrocytes are terminally differentiated cells that sustain cascades of gene activation and repression to execute highly specialized functions, while retaining homeostatic control. To study long-range chromatin folding without disturbing the native tissue environment, we developed Genome Architecture Mapping in combination with immunoselection (immunoGAM), and applied it to three cell types from the adult murine brain: dopaminergic neurons (DNs) from the midbrain, pyramidal glutamatergic neurons (PGNs) from the hippocampus, and oligodendroglia (OLGs) from the cortex. We find cell-type specific 3D chromatin structures that relate with patterns of gene expression at multiple genomic scales, including extensive reorganization of topological domains (TADs) and chromatin compartments. We discover the loss of TAD insulation, or  TAD melting, at long genes (>400 kb) when they are highly transcribed. We find many neuron-specific contacts which contain accessible chromatin regions enriched for putative binding sites for multiple neuronal transcription factors, and which connect cell-type specific genes that are associated with neurodegenerative disorders such as Parkinsons disease, or specialized functions such as synaptic plasticity and memory. Lastly, sensory receptor genes exhibit increased membership in heterochromatic compartments that establish strong contacts in brain cells. However, their silencing is compromised in a subpopulation of PGNs with molecular signatures of long-term potentiation. Overall, our work shows that the 3D organization of the genome is highly cell-type specific, and essential to better understand mechanisms of gene regulation in highly specialized tissues such as the brain.


title: Perceptual decisions and oculomotor responses rely on temporally distinct streams of evidence

abstract: Perceptual decisions often require the integration of noisy sensory evidence over time. This process is formalized with sequential sampling models, where evidence is accumulated up to a decision threshold before a choice is made. Although intuition suggests that decision formation must precede the preparation of a motor response (i.e., the action used to communicate the choice), neurophysiological findings have suggested that these two processes might be one and the same. To test this idea, we developed a reverse-correlation protocol in which the visual stimuli that influence decisions can be distinguished from those guiding motor responses. In three experiments, we found that the temporal weighting function of oculomotor responses did not overlap with the relatively early weighting function of stimulus properties having an impact on decision formation. These results support a timeline in which perceptual decisions are formed, at least in part, prior to the preparation of a motor response.


title: Homeostatic scaling is driven by a translation-dependent degradation axis that recruits miRISC remodelling

abstract: Homeostatic scaling in neurons has been attributed to the individual contribution of either translation or degradation; however there remains limited insight towards understanding how the interplay between the two processes effectuates synaptic homeostasis. Here, we report that a co-dependence between protein synthesis and degradation mechanisms drives synaptic homeostasis whereas abrogation of either prevents it. Coordination between the two processes is achieved through the formation of a tripartite complex between translation regulators, the 26S proteasome and the miRNA-induced-silencing-complex (miRISC) components such as Argonaute, MOV10 and Trim32 on actively translating transcripts or polysomes. The components of this ternary complex directly interact with each other in an RNA-dependent manner. Disruption of polysomes abolishes this ternary interaction, suggesting that translating RNAs facilitate the combinatorial action of the proteasome and the translational apparatus. We identify that synaptic downscaling involves miRISC remodeling which entails the mTORC1-dependent translation of Trim32, an E3 ligase and the subsequent degradation of its target, MOV10 via the phosphorylation of p70 S6 kinase. We find that the E3 ligase Trim32 specifically polyubiquitinates MOV10 for its degradation during synaptic downscaling. MOV10 degradation alone is sufficient to invoke downscaling by enhancing Arc translation through its 3 UTR and causing the subsequent removal of post-synaptic AMPA receptors. Synaptic scaling was occluded when we depleted Trim32 and overexpressed MOV10 in neurons, suggesting that the Trim32-MOV10 axis is necessary for synaptic downscaling. We propose a mechanism that exploits a translation-driven protein degradation paradigm to invoke miRISC remodeling and induce homeostatic scaling during chronic network activity.


title: Exploring the understudied human kinome for research and therapeutic opportunities

abstract: The functions of protein kinases have been widely studied and over 60 kinase inhibitors are FDA-approved drugs. Membership in the human kinome is nonetheless subject to multiple overlapping and inconsistent definitions and is unevenly studied, complicating functional genomics and chemical genetics. We describe objective criteria for refining the definition of the human kinome to comprise an extended set of 710 kinase domains and a more narrowly curated set of 557 protein kinase like (PKL) domains. An online tool (www.kinome.org) makes it possible to sort these sets on multiple structural and functional criteria. Focusing on the least studied one-third of the kinome we find that many proteins are differentially expressed, essential in multiple cell lines, and mutated in the Cancer Genome Atlas. We show that some understudied kinases are high affinity off-targets of clinical-grade compounds and approved drugs and we describe an optimized small molecule library making use of this information for selective kinome perturbation. We conclude that the understudied kinome contains physiologically important proteins, including possible targets for future drug discovery campaigns.

category: systems biology


title: Pro-metastatic gene expression, immune evasion and an altered HPV spectrum characterize an aggressive subtype of cervical cancer

abstract: Human papillomavirus (HPV)-associated cervical cancer represents one of the leading causes of cancer death worldwide. Although low-middle income countries are disproportionately affected, our knowledge of the disease predominantly originates from populations in high-income countries. Using the largest multi-omic analysis of cervical squamous cell carcinoma (CSCC) to date, totalling 643 tumours and representing patient populations from the USA, Europe and Sub-Saharan Africa, we identify two CSCC subtypes (C1 and C2) with differing prognosis. C1 tumours are largely HPV16-driven, display increased cytotoxic T-lymphocyte infiltration and frequently harbour PIK3CA and EP300 mutations. C2 tumours are associated with shorter overall survival, are frequently driven by HPVs from the HPV18-containing alpha-7 clade, harbour alterations in the Hippo signalling pathway and increased expression of immune checkpoint genes, B7-H3 (also known as CD276) and NT5E (also known as CD73) and PD-L2 (also known as PDCD1LG2). In conclusion, we identify two novel, therapy-relevant CSCC subtypes that share the same defining characteristics across three geographically diverse cohorts.


title: Integrated analysis of cervical squamous cell carcinoma cohorts from three continents reveals conserved subtypes of prognostic significance.


title: Admixture has obscured signals of historical hard sweeps in humans

abstract: The role of natural selection in shaping biological diversity is an area of intense interest in modern biology. To date, studies of positive selection have primarily relied upon genomic datasets from contemporary populations, which are susceptible to confounding factors associated with complex and often unknown aspects of population history. In particular, admixture between diverged populations can distort or hide prior selection events in modern genomes, though this process is not explicitly accounted for in most selection studies despite its apparent ubiquity in humans and other species. Through analyses of ancient and modern human genomes, we show that previously reported Holocene-era admixture has masked more than 50 historic hard sweeps in modern European genomes. Our results imply that this canonical mode of selection has likely been underappreciated in the evolutionary history of humans and suggests that our current understanding of the tempo and mode of selection in natural populations may be quite inaccurate.


title: Systematic analysis of naturally occurring insertions and deletions that alter transcription factor spacing identifies tolerant and sensitive transcription factor pairs

abstract: Regulation of gene expression requires the combinatorial binding of sequence-specific transcription factors (TFs) at promoters and enhancers. Prior studies showed that alterations in the spacing between TF binding sites can influence promoter and enhancer activity. However, the relative importance of TF spacing alterations resulting from naturally occurring insertions and deletions (InDels) has not been systematically analyzed. To address this question, we first characterized the genome-wide spacing relationships of 75 TFs in K562 cells as determined by ChIP-sequencing. We found a dominant pattern of a relaxed range of spacing between collaborative factors, including 46 TFs exclusively exhibiting relaxed spacing with their binding partners. Next, we exploited millions of InDels provided by genetically diverse mouse strains and human individuals to investigate the effects of altered spacing on TF binding and local histone acetylation. Spacing alterations resulting from naturally occurring InDels are generally tolerated in comparison to genetic variants directly affecting TF binding sites. A remarkable range of tolerance was further established for PU.1 and C/EBP{beta}, which exhibit relaxed spacing, by introducing synthetic spacing alterations ranging from 5-bp increase to >30-bp decrease using CRISPR/Cas9 mutagenesis. These findings provide implications for understanding mechanisms underlying enhancer selection and for the interpretation of non-coding genetic variation.


title: STEP inhibition prevents Aβ mediated damage in dendritic complexity and spine density in Alzheimer's disease.

abstract: Loss of dendritic spines and decline of cognitive function are hallmarks of patients with Alzheimers disease (AD). Previous studies have shown that AD pathophysiology involves increased expression of a central nervous system-enriched protein tyrosine phosphatase called STEP (STriatal-Enriched protein tyrosine Phosphatase). STEP opposes the development of synaptic strengthening by dephosphorylating substrates, including GluN2B, Pyk2 and ERK1/2. Genetic reduction of STEP as well as pharmacological inhibition of STEP improves cognitive function and hippocampal memory in the 3xTg AD mouse model. Here, we show that the improved cognitive function is accompanied by an increase in synaptic connectivity in cell cultures as well as in the triple transgenic AD mouse model, further highlighting the potential of STEP inhibitors as a therapeutic agent.


title: Chromatin dynamics during hematopoiesis reveal discrete regulatory modules instructing differentiation

abstract: Lineage commitment and differentiation is driven by the concerted action of master transcriptional regulators at their target chromatin sites. Multiple efforts have characterized the key transcription factors (TFs) that determine the various hematopoietic lineages. However, the temporal interactions between individual TFs and their chromatin targets during differentiation and how these interactions dictate lineage commitment remains poorly understood. We performed dense, daily, temporal profiling of chromatin accessibility (DNase I-seq) and gene expression changes (total RNA-seq) along ex vivo human erythropoiesis to comprehensively define developmentally regulated DNase I hypersensitive sites (DHSs) and transcripts. We link both distal DHSs to their target gene promoters and individual TFs to their target DHSs, revealing that the regulatory landscape is organized in distinct sequential regulatory modules that regulate lineage restriction and maturation. Finally, direct comparison of transcriptional dynamics (bulk and single-cell) and lineage potential between erythropoiesis and megakaryopoiesis uncovers differential fate commitment dynamics between the two lineages as they exit pluripotency. Collectively, these data provide novel insights into the global regulatory landscape during hematopoiesis.


title: Size-dependent eco-evolutionary feedbacks in fisheries

abstract: AO_SCPLOWBSTRACTC_SCPLOWHarvesting may drive body downsizing along with population declines and decreased harvesting yields. These changes are commonly construed as consequences of direct harvest selection, where small-bodied, early-reproducing individuals are immediately favoured. However, together with directly selecting against a large body size, harvesting and body downsizing alter many ecological features, such as competitive and trophic interactions, and thus also indirectly reshape natural selection acting back on body sizes through eco-evolutionary feedback loops (EEFLs). We sketch plausible scenarios of simple EEFLs in which one-dimensional, density-dependent natural selection acts either antagonistically or synergistically with direct harvest selection on body size. Antagonistic feedbacks favour body-size stasis but erode genetic variability and associated body-size evolvability, and may ultimately impair population persistence and recovery. In contrast, synergistic feedbacks drive fast evolution towards smaller body sizes and favour population resilience, but may have far-reaching bottom-up or top-down effects. We illustrate the further complexities resulting from multiple environmental feedbacks using a co-evolving predator-prey pair, in which case outcomes from EEFLs depend not only on population densities, but also on whether prey sit above or below the optimal predator/prey body-size ratio, and whether prey are more or less evolvable than their predators. EEFLs improve our ability to understand and predict natures response to harvesting, but their integration into the research agenda will require a full consideration of the effects and dynamics of natural selection.


title: Distributed coding of stimulus magnitude in rodent prefrontal cortex

abstract: As we interact with the external world, we judge magnitudes from sensory information. The estimation of magnitudes has been characterized in primates, yet it is largely unexplored in non-primate species. Here, we show that gerbils that solve a time-interval reproduction task display primate-like magnitude estimation characteristics, most prominently a systematic overestimation of small stimuli and an underestimation of large stimuli, often referred to as regression effect. We investigated the underlying neural mechanisms by recording from medial prefrontal cortex and show that the majority of neurons respond either during the measurement or the reproduction of a time-interval. Cells that are active during both phases display distinct response patterns. We categorize the neural responses into multiple types and demonstrate that only populations with mixed responses can encode the bias of the regression effect. These results reveal the organizing neural principles of magnitude estimation a process important to higher cognitive functions.


title: Low-frequency neural tracking of speech envelope reflects the convolution of evoked responses to acoustic edges, not oscillatory entrainment

abstract: The amplitude envelope of speech is crucial for accurate comprehension, and several studies have shown that the phase of neural activity in the theta-delta bands (1 - 10 Hz) tracks the phase of the speech amplitude envelope during listening, a process referred to as envelope tracking. However, the mechanisms underlying envelope tracking have been heavily debated. A dominant model posits that envelope tracking reflects continuous entrainment of endogenous low-frequency oscillations to the speech envelope. However, it has proven challenging to distinguish this from the alternative model that envelope tracking reflects a convolution of evoked responses to acoustic landmarks within the envelope. To address this, we recorded magnetoencephalography while participants listened to natural and slowed speech. First, we found that peaks in the rate of envelope change, acoustic edge landmarks, induced a strong evoked response. We then contrasted central diverging predictions of evoked response and oscillatory models regarding the spectral and temporal extent of neural phase locking in the theta-delta range at different speech rates. Our analyses revealed transient theta phase locking at both speech rates and delta phase locking in slow speech. This pattern was predicted by the evoked response model but not by oscillatory entrainment. Furthermore, the encoding of acoustic edge magnitudes was invariant to contextual speech rate, demonstrating speech rate normalization of acoustic edge representations. Taken together, our results show that an evoked response model provides a better account of neural phase locking to the speech envelope than oscillatory entrainment.