title: Generation, characterization and drug sensitivities of twelve patient-derived IDH1 mutant glioma cell cultures

abstract:
BackgroundMutations of the isocitrate dehydrogenase (IDH) gene occur in over 80% of low-grade gliomas and secondary glioblastomas. Despite considerable efforts, endogenous in vitro IDH-mutated glioma models remain scarce. Availability of these models is key for the development of new therapeutic interventions.

MethodsCell cultures were established from fresh tumor material and expanded in serum-free culture media. D-2-Hydroxyglutarate levels were determined by mass-spectrometry. Genomic and transcriptomic profiling were carried out on the Illumina Novaseq platform, methylation profiling was performed with the Infinium MethylationEpic BeadChip array. Mitochondrial respiration was measured with the Seahorse XF24 Analyzer. Drug screens were performed with an NIH FDA-approved anti-cancer drug set and two IDH-mutant specific inhibitors.

ResultsA set of twelve patient-derived IDHmt cell cultures was established. We confirmed high concordance in driver mutations, copy number and methylation profiles between the tumors and derived cultures. Homozygous deletion of CDKN2A/B was observed in all cultures. IDH-mutant cultures had lower mitochondrial reserve capacity. IDH-mutant specific inhibitors did not affect cell viability or global gene expression. Screening of 107 FDA-approved anti-cancer drugs identified nine compounds with potent activity against IDHmt gliomas, including three compounds with favorable pharmacokinetic characteristics for CNS penetration: teniposide, omacetaxine mepesuccinate, and marizomib.

ConclusionsOur twelve IDH-mutant cell cultures show high similarity to the parental tissues and offer a unique tool to study the biology and drug sensitivities of high-grade IDHmt gliomas in vitro. Our drug screening studies reveal lack of sensitivity to IDHmt inhibitors, but sensitivity to a set of nine available anti-cancer agents.

Key pointsO_LIIDHmt glioma cultures closely resemble their parental tumors
C_LIO_LIMicroscopic monitoring of early passages and colony isolation increases IDH1mt culture success
C_LIO_LIDrug screening identified nine candidate repurposed drugs for IDHmt glioma
C_LI

Importance of the studyIDH-mutations are highly prevalent in low grade and secondary high-grade gliomas. Despite this high frequency however, very few in vitro models have been reported for IDH-mutated gliomas. In this manuscript we describe and characterize in detail twelve primary cultures from IDH-mutant astrocytomas. We show that these cultures retain most of the genetic, epigenetic and metabolic features of their respective parental tumors. Because of these similarities, these independent model systems will not only help understand the molecular defects driven by the mutation, but are also vital to identify means to target these tumors. Screening of 107 FDA-approved anti-cancer agents on these cultures identified a set of highly effective agents that may offer candidates for either systemic or assisted delivery treatment of this tumor subtype.

category: cancer biology


title: Intergenerational Stress Transmission is Associated with Brain Metabotranscriptome Remodeling andMitochondrial Dysfunction

abstract:
Intergenerational stress increases lifetime susceptibility to depression and other psychiatric disorders. Whether intergenerational stress transmission is a consequence of in utero neurodevelopmental disruptions vs early-life mother-infant interaction is largely unknown. Here, we demonstrated that exposure to traumatic stress in mice during pregnancy, through predator scent exposure, induces in the offspring social deficits and depressive-like behavior. We found, through cross-fostering experiments, that raising of normal pups by traumatized mothers produced a similar behavioral phenotype to that induced in pups raised by their biological traumatized mothers. Good caregiving (by non-traumatized mothers), however, did not completely protect against the prenatal trauma-induced behavioral deficits. These findings support a two-hit stress mechanism of both in utero and early-life parenting (poor caregiving by the traumatized mothers) environments. Associated with the behavioral deficits, we found profound changes in brain metabolomics and transcriptomic (metabotranscriptome). Striking increases in the mitochondrial hypoxia marker and epigenetic modifier 2-hydroxyglutaric acid, in the brains of neonatal and adult pups whose mothers were exposed to stress during pregnancy, indicated mitochondrial metabolism dysfunctions and epigenetic mechanisms. Bioinformatic analyses revealed mechanisms involving stress- and hypoxia-response metabolic pathways in the brains of the neonatal mice, which appear to lead to long-lasting alterations in mitochondrial-energy metabolism, and epigenetic processes pertaining to DNA and chromatin modifications. Most strikingly, we demonstrated that an early pharmacological intervention that can correct mitochondria metabolism - lipid metabolism and epigenetic modifications with acetyl-L-carnitine (ALCAR) supplementation - produces long-lasting protection against the behavioral deficits associated with intergenerational transmission of traumatic stress.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=129 SRC="FIGDIR/small/438868v1_ufig1.gif" ALT="Figure 1">
View larger version (40K):
org.highwire.dtl.DTLVardef@1fb2500org.highwire.dtl.DTLVardef@13a5c6forg.highwire.dtl.DTLVardef@8aa512org.highwire.dtl.DTLVardef@5b6906_HPS_FORMAT_FIGEXP  M_FIG C_FIG

category: neuroscience


title: Efficacy and breadth of adjuvanted SARS-CoV-2 receptor-binding domain nanoparticle vaccine in macaques

abstract:
Emergence of novel variants of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) underscores the need for next-generation vaccines able to elicit broad and durable immunity. Here we report the evaluation of a ferritin nanoparticle vaccine displaying the receptor-binding domain of the SARS-CoV-2 spike protein (RFN) adjuvanted with Army Liposomal Formulation QS-21 (ALFQ). RFN vaccination of macaques using a two-dose regimen resulted in robust, predominantly Th1 CD4+ T cell responses and reciprocal peak mean neutralizing antibody titers of 14,000-21,000. Rapid control of viral replication was achieved in the upper and lower airways of animals after high-dose SARS-CoV-2 respiratory challenge, with undetectable replication within four days in 7 of 8 animals receiving 50 {micro}g RFN. Cross-neutralization activity against SARS-CoV-2 variant B.1.351 decreased only [~]2-fold relative to USA-WA1. In addition, neutralizing, effector antibody and cellular responses targeted the heterotypic SARS-CoV-1, highlighting the broad immunogenicity of RFN-ALFQ for SARS-like betacoronavirus vaccine development.

Significance StatementThe emergence of SARS-CoV-2 variants of concern (VOC) that reduce the efficacy of current COVID-19 vaccines is a major threat to pandemic control. We evaluate a SARS-CoV-2 Spike receptor-binding domain ferritin nanoparticle protein vaccine (RFN) in a nonhuman primate challenge model that addresses the need for a next-generation, efficacious vaccine with increased pan-SARS breadth of coverage. RFN, adjuvanted with a liposomal-QS21 formulation (ALFQ), elicits humoral and cellular immune responses exceeding those of current vaccines in terms of breadth and potency and protects against high-dose respiratory tract challenge. Neutralization activity against the B.1.351 VOC within two-fold of wild-type virus and against SARS-CoV-1 indicate exceptional breadth. Our results support consideration of RFN for SARS-like betacoronavirus vaccine development.

category: immunology


title: Atypical N-glycosylation of SARS-CoV-2 impairs the efficient binding of Spike-RBM to the human-host receptor hACE2

abstract: SARS-CoV-2 internalization by human host cells relies on the molecular binding of its spike glycoprotein (SGP) to the angiotensin-converting-enzyme-2 (hACE2) receptor. It remains unknown whether atypical N-glycosylation of SGP modulates SARS-CoV-2 tropism for infections. Here, we address this question through an extensive bioinformatics analysis of publicly available structural and genetic data. We identified two atypical sequons (sequences of N-glycosylation: NGV 481-483 and NGV 501-503), strategically located on the receptor-binding motif (RBM) of SGP and facing the hACE2 receptor. Interestingly, the cryo-electron microscopy structure of trimeric SGP in complex with potent-neutralizing antibodies from convalescent patients revealed covalently-linked N-glycans in NGV 481-483 atypical sequons. Furthermore, NGV 501-503 atypical sequon involves the asparagine-501 residue, whose highly-transmissible mutation N501Y is present in circulating variants of major concerns and affects the SGP-hACE2 binding-interface through the well-known hotspot-353. These findings suggest that atypical SGP post-translational modifications modulate the SGP-hACE2 binding-affinity affecting consequently SARS-CoV-2 transmission and pathogenesis.

category: molecular biology


title: Subjective confidence reflects representation of Bayesian probability in cortex

abstract: What gives rise to the human sense of confidence? Here, we tested the Bayesian hypothesis that confidence is based on a probability distribution represented in neural population activity. We implemented several computational models of confidence, and tested their predictions using psychophysics and fMRI. Using a generative model-based fMRI decoding approach, we extracted probability distributions from neural population activity in human visual cortex. We found that subjective confidence tracks the shape of the decoded distribution. That is, when sensory evidence was more precise, as indicated by the decoded distribution, observers reported higher levels of confidence. We furthermore found that neural activity in the insula, anterior cingulate, and prefrontal cortex was linked to both the shape of the decoded distribution and reported confidence, in ways consistent with the Bayesian model. Altogether, our findings support recent statistical theories of confidence and suggest that probabilistic information guides the computation of ones sense of confidence.


title: Escherichia coli Cryptic Prophages Sense Nutrients to Control Persister Cell Resuscitation

abstract: We determined previously that some cryptic prophages are not genomic junk but instead enable cells to combat myriad stresses as part of an active stress response. However, how these phage fossils affect the extreme stress response of dormancy; i.e., how cryptic prophages affect persister cell formation and resuscitation, has not been fully explored. Persister cells form as a result of stresses such as starvation, antibiotics, and oxidative conditions, and resuscitation of these persister cells likely causes recurring infections such as those associated with tuberculosis, cystic fibrosis, and Lyme disease. Unlike for the active stress response, here we find that deletion of each of the nine Escherichia coli cryptic prophages has no effect on persister cell formation. Strikingly, elimination of each cryptic prophage results in an increase in persister cell resuscitation with a dramatic increase in resuscitation upon deleting all nine prophages. This increased resuscitation includes eliminating the need for a carbon source and is due to activation of the phosphate import system as a result of inactivating transcriptional regulator AlpA of the CP4-57 cryptic prophage, since we found {Delta}alpA increases persister resuscitation, and AlpA represses phosphate regulator PhoR. Therefore, we report a novel cellular stress mechanism controlled by cryptic prophages: regulation of phosphate uptake which controls the exit of the cell from dormancy and prevents premature resuscitation in the absence of nutrients.

category: microbiology


title: RNA-seq of the medusavirus suggests remodeling of the host nuclear environment at an early infection stage

abstract:
Nucleo-cytoplasmic large DNA viruses (NCLDVs) undergo a cytoplasmic or nucleo-cytoplasmic cycle, and the latter involves both nuclear and cytoplasmic compartments to proceed viral replication. Medusavirus, a recently isolated NCLDV, has a nucleo-cytoplasmic replication cycle in amoebas during which the host nuclear membrane apparently remains intact, a unique feature among amoeba-infecting giant viruses. The medusavirus genome lacks most transcription genes but encodes a full set of histone genes. To investigate the infection strategy, we performed a time-course RNA-seq experiment. All the viral genes were transcribed and classified into five temporal expression clusters. The immediate early genes (cluster 1, 42 genes) were mostly (83%) of unknown functions, frequently (95%) associated with a palindromic promoter-like motif, and enriched (45%) in putative nuclear-targeting genes. The later genes (clusters 2-5) were assigned to various functional categories. The viral linker histone H1 gene was in cluster 1, whereas the four core histone genes were in cluster 3, suggesting they had distinct roles during the course of the virus infection. The transcriptional profile of the host amoeba, Acanthamoeba castellanii, genes was greatly altered post-infection. Several encystment-related host genes showed increased representation levels at 48 hours post-infection, which is consistent with the previously reported amoeba encystment upon medusavirus infection. Overall, the transcriptional landscape during the course of medusavirus infection suggests that the virus modifies the host nuclear environment immediately after the initiation of infection. -

ImportanceMedusavirus is an amoeba-infecting giant virus that was isolated from a hot spring in Japan. It belongs to the proposed family "Medusaviridae" in the phylum Nucleocytoviricota. Unlike other amoeba-infecting giant viruses, medusavirus initiates its DNA replication in the host nucleus without disrupting the nuclear membrane. Our RNA-seq analysis of its infection course uncovered ordered viral gene expression profiles. We identified temporal expression clusters of viral genes and associated putative promoter motifs. The subcellular localization prediction showed a clear spatiotemporal correlation between gene expression timing and localization of the encoded proteins. Notably, the immediate early expression cluster was enriched in genes targeting the nucleus, suggesting the priority of remodeling the host intra-nuclear environment during infection. The transcriptional profile of the amoeba genes was greatly altered post-infection. Notably, the expression of encystment-related genes increased 48 hours post-infection, suggesting that encystment may be an antiviral strategy of amoeba.


title: Plasmid fitness costs are caused by specific genetic conflicts

abstract: Plasmids play an important role in bacterial genome evolution by transferring genes between lineages. Fitness costs associated with plasmid acquisition are expected to be a barrier to gene exchange, but the causes of plasmid fitness costs are poorly understood. Single compensatory mutations are often sufficient to completely ameliorate plasmid fitness costs, suggesting that such costs are caused by specific genetic conflicts rather than generic properties of plasmids, such as their size, metabolic burden, or expression level. Here we show -- using a combination of experimental evolution, reverse genetics, and transcriptomics -- that fitness costs of two divergent large plasmids in Pseudomonas fluorescens are caused by inducing maladaptive expression of a chromosomal tailocin toxin operon. Mutations in single genes unrelated to the toxin operon, and located on either the chromosome or the plasmid, ameliorated the disruption associated with plasmid acquisition. We identify one of these compensatory loci, the chromosomal gene PFLU4242, as the key mediator of the fitness costs of both plasmids, with the other compensatory loci either reducing expression of this gene or mitigating its deleterious effects by upregulating a putative plasmid-borne ParAB operon. The chromosomal mobile genetic element Tn6291, which uses plasmids for transmission, remained upregulated even in compensated strains, suggesting that mobile genetic elements communicate through pathways independent of general physiological disruption. Plasmid fitness costs caused by specific genetic conflicts are unlikely to act as a long-term barrier to horizontal gene transfer due to their propensity for amelioration by single compensatory mutations, explaining why plasmids are so common in bacterial genomes.


title: Discovery of bactericides as an acute mitochondrial membrane damage inducer

abstract: Mitochondria evolved from endosymbiotic bacteria to become essential organelles of eukaryotic cells. The defined lipid composition and structure of mitochondrial membranes are critical for the proper functioning of mitochondria. However, mitochondrial stress responses that help maintain the integrity of mitochondrial membranes against internal or external insults are not well understood. One reason for this lack of insight is the absence of efficient tools to specifically damage mitochondrial membranes. Here, through a compound screen originally aimed at identifying inhibitors of the inner mitochondrial membrane (IMM)-resident protease OMA1, we found that two bis-biguanide compounds, Chlorhexidine and Alexidine, modified OMA1 activity by altering the integrity of the IMM. Interestingly, these compounds are well-known bactericides whose mechanism of action has centered on their damage-inducing activity on bacterial membranes. We found Alexidine binds to the IMM likely through the electrostatic interaction driven by the membrane potential as well as an affinity for anionic phospholipids. Electron microscopic analysis revealed that Alexidine severely perturbated the IMM, especially the cristae structure. Along with this, we observed the altered localization of IMM-resident membrane-shaping proteins, including Mic60. Notably, Alexidine evoked a specific transcriptional/proteostasis signature that was not induced by other typical mitochondrial stressors, highlighting the unique property of Alexidine as a novel mitochondrial membrane stressor. Our findings provide a chemical-biological tool that can induce acute and selective perturbation of the IMM integrity, which should enable the delineation of mitochondrial stress-signaling pathways required to maintain the mitochondrial membrane homeostasis.

category: cell biology


title: Nej1(XLF) inhibits Sae2(CTIP)-Dna2(DNA2) mediated resection at DNA double strand breaks

abstract:
The two major pathways of DNA double strand break (DSB) repair, non-homologous end-joining (NHEJ) and homologous recombination (HR), are highly conserved from yeast to mammals. The regulation of 5 DNA resection controls repair pathway choice and influences repair outcomes. Nej1 was first identified as a canonical NHEJ factor involved in stimulating the ligation of broken DNA ends and more recently, it was shown to be important for DNA end-bridging and inhibiting 5 resection mediated by Dna2-Sgs1. Nej1 interacts with Sae2 and this impacts DSB repair in three ways. First, Nej1 inhibits MRX-Sae2 interactions and Sae2 localization to a DSB. Second, Nej1 inhibits Sae2-dependent recruitment of Dna2 in the absence of Sgs1. Third, NEJ1 and SAE2 showed an epistatic relationship for DNA end-bridging, an event that restrains the broken ends and reduces the frequency of genomic deletions from developing at the DSB. Deletion of NEJ1 suppressed the synthetic lethality of sae2{Delta} sgs1{Delta} and was dependent on the nuclease activity of Dna2. These Nej1 functions promote end-joining DSB repair, but could also be relevant for controlling resection initiation during HR repair.

HighlightsO_LINej1 physically interacts with Sae2 and inhibits end-resection at a DSB.
C_LIO_LINej1 inhibits Sae2 interactions with the MRX complex.
C_LIO_LINej1 inhibits Sae2-dependent recruitment of Dna2 to a DSB.
C_LIO_LINEJ1 and SAE2 are epistatic for DNA end-bridging.
C_LIO_LIDeletion of NEJ1 suppresses the synthetic lethality of sae2{Delta} sgs1{Delta}, which is dependent on Dna2 nuclease activity.
C_LI

category: genetics


title: Sgs1BLM independent role of Dna2DNA2 nuclease at DNA double strand break is inhibited by Nej1XLF


title: Nej1XLF interacts with Sae2CTIP to inhibit Dna2 mediated resection at DNA doublestrand break


title: Genetic regulation of homeostatic immune architecture in the lungs of Collaborative Cross mice

abstract: Variation in immune homeostasis, immune system stability, in organ systems such as the lungs is likely to shape the host response to infection at these exposed tissues. We evaluated immune homeostasis in immune cell populations in the lungs of the Collaborative Cross (CC) mouse genetic reference population. We found vast heritable variation in leukocyte populations with the frequency of many of these cell types showing distinct patterns relative to classic inbred strains C57BL/6J and BALB/cJ. We identified 28 quantitative trait loci (QTL) associated with variation in baseline lung immune cell populations, including several loci that broadly regulate the abundance of immune populations from distinct developmental lineages, and found that many of these loci have predictive value for influenza disease outcomes, demonstrating that genetic determinants of homeostatic immunity in the lungs regulate susceptibility to virus-induced disease. All told, we highlight the need to assess diverse mouse strains in understanding immune homeostasis and resulting immune responses.


title: Barcoded Competitive Clone-Initiating Cell (BC-CIC) Analysis Reveals Differences in Ovarian Cancer Cell Genotype and Niche Specific Clonal Fitness During Growth and Metastasis In Vivo

abstract: During oncogenesis, pathogenic clones develop which contain cells capable of spreading throughout the body, ultimately compromising vital organ functions and physiology. Understanding how metastatic clones develop and spread is critical for improving cancer treatments. However, our understanding of these processes has been hampered by a paucity of quantitative methodologies to comprehensively map, track and characterize such clones. To address this shortcoming, we have developed a DNA barcoding and next-generation sequencing based system-wide clonal tracking technology integrated with a computational data analysis pipeline called Clone-Initiating Cell (CIC) Calculator. The CIC Calculator interfaces with the CIC Morbus Mandala (CIC-MM) plot, a novel tool to visually comprehend and detect four distinct categories that explains their complex relationships with various tissues/organ sites. Further, we describe machine learning approaches to study CIC number, frequency, and estimate clone size and distribution demonstrating distinct growth patterns, and their inter-relationships and their routes of metastatic spread at clonal resolution. We demonstrate these methodologies, using our novel multifunctional lentiviral barcode libraries, and specifically barcoded tubal-ovarian metastatic OVCAR5 cell lines (engineered to express varying levels of metastasis promoting LRRC15 gene) and co-injected cells in a competitive CIC assay into tubal or ovarian sites in highly immunodeficient NSG mice. DNA was isolated from primary tumors, omental/bowel metastasis and system-wide anatomical site/organs. Amplicon sequencing libraries were constructed with spike-in-control barcodes (serving as internal calibration controls) to estimate absolute clone sizes. The computational pipeline CIC Calculator was then used to deconvolute and filter the data, set stringent thresholds, and generate high-quality information on CIC numbers and frequencies, clone sizes, linkages across sites and classify clones based on their extent of metastatic activity. Using of CIC-MM plot, statistical models and machine learning approaches, we generated high-resolution clonal maps of metastasis for each animal. The information generated included clone types and system-wide metastasis, similar and dissimilar clonal patterns of dominance at heterotopic sites and their routes of metastases. The data revealed previously unknown influences of cellular genotype and their implanted sites on selecting certain clones with specific system-wide clonal patterns, and identified rare LRRC15 expressor clones (classified as CIC.Toti) predisposed to exploit  all sites, albeit at varying degrees of dominance. The genomic technology and computational methodology described here are tissue-agnostic. They enable rapid adoption for an investigation into various stages of system-wide metastasis and growth of transplantable malignant cells at the highest clonal resolution.


title: Differential immunopeptidome analysis revealed cancer specific amino acid usage of HLA class-I antigens and novel neoantigens of colorectal cancer

abstract: Understanding the properties of human leukocyte antigen (HLA) peptides (HLAp) is essential for cancer precision medicine, while the direct identification of HLAp from a tiny piece of clinical tissues by mass spectrometry (MS) is still confronted with technical challenges. Here, to overcome that hindrances, high-field asymmetric waveform ion mobility spectrometry (FAIMS) is introduced to conduct differential ion mobility (DIM)-MS by seamless gas-phase fractionation optimal for scarce samples. The global-immunopeptidomics analysis enables in-depth identification of 44,785 HLAp, including 2 neoantigens with KRAS-G12V and CPPED1-R228Q from stage IV colorectal cancer (CRC) tissues obtained from 17 subjects. Comparison of tissue-based personal immunopeptidome disclosed cancer-specific processing of HLAp. Since the direct identification of neoantigens from CRC tissues suggests more potential neoantigens have yet to been identified, the targeted-immunopeptidomics by screening-oriented parallel reaction monitoring (PRM) has been established and identified two more neoantigens with oncogenic KRAS. The targeted-immunopeptidomics analysis clarifies the presentation of aimed neoantigens, including neoantigens that are unpredictable by computational algorithms. The comprehensive-immunopeptidomics approach combining the FAIMS-assisted DIM-MS and the screening-oriented PRM contributes the direct and effective identification of HLAp and emerging neoantigens from clinical tissue samples, as an advanced strategy to identify the authentic and recurrent neoantigens for cancer immunotherapy.


title: Comprehensive-immunopeptidomics analysis reveals presentation of potential neoantigens carrying cancer driver mutations


title: State dependence of neural networks on past history and stimulus presentation in Caenorhabditis elegans

abstract:
Neural responses are influenced by both external stimuli and internal network states. While network states have been linked to behavioral and stimulus states, little is known about how sensory inputs are filtered by whole-brain activity to affect motor and command neurons. Here, we recorded whole-brain activity of Caenorhabditis elegans experiencing bacterial food stimuli, and modeled how sensory inputs affect motor and command neurons in a network state-dependent manner. First, we classified active neurons into six functional clusters: two sensory neuron clusters (ON, OFF), and four motor/command neuron clusters (AVA, RME, SMDD, SMDV). Using encoding models, we found that ON and OFF sensory neurons that respond to onset and removal of bacteria, respectively, employ different adaptation strategies. Next, we used decoding models to show that bacterial onset and removal differentially drive AVA and RME cluster activity. To explore state-dependent effects on AVA and RME clusters, we developed a model that identified network states and fitted submodels for each state to predict how each of the six functional clusters drive AVA and RME cluster activity. We also identified network states in which AVA and RME clusters were either largely unperturbed by or receptive to bacterial sensory input. Furthermore, this model allowed us to disentangle the state-dependent contributions of stimulus timescales and bacterial content to neural activity. Collectively, we present an interpretable approach for modeling network dynamics that goes beyond implication of neurons in particular states, and moves toward explicitly dissecting how neural populations work together to produce state dependence.

Significance StatementA major function of the brain is to transform sensory information into behavior. As the first receiver of sensory input, sensory neuron activity is often most correlated with stimulus features. However, this high-fidelity representation of sensory input becomes diluted as it travels to downstream neurons, where sensory information is integrated with network activity. By the time sensory information reaches motor neurons, it is often difficult to dissociate the influence of sensory input from the influence of network activity. Here, we describe a method that is fully interpretable such that we can show how neural populations on a whole-brain scale interact to produce network states. From there, we can attribute motor neuron activity to network history and sensory input.


title: State-dependent network interactions differentially gate sensory input at the motor and command neuron level in Caenorhabditis elegans


title: Intrinsic Mechanisms in the Gating of Resurgent Na+ Currents

abstract:
The resurgent component of the voltage-gated sodium current (INaR) is a depolarizing conductance, revealed on membrane hyperpolarizations following brief depolarizing voltage steps, which has been shown to contribute to regulating the firing properties of numerous neuronal cell types throughout the central and peripheral nervous systems. Although mediated by the same voltage-gated sodium (Nav) channels that underlie the transient and persistent Nav current components, the gating mechanisms that contribute to the generation of INaR remain unclear. Here, we characterized Nav currents in mouse cerebellar Purkinje neurons, and used tailored voltage-clamp protocols to define how the voltage and the duration of the initial membrane depolarization affect the amplitudes and kinetics of INaR. Using the acquired voltage-clamp data, we developed a novel Markov kinetic state model with parallel (fast and slow) inactivation pathways and, we show that this model reproduces the properties of the resurgent, as well as the transient and persistent, Nav currents recorded in (mouse) cerebellar Purkinje neurons. Based on the acquired experimental data and the simulations, we propose that resurgent Na+ influx occurs as a result of fast inactivating Nav channels transitioning into an open/conducting state on membrane hyperpolarization, and that the decay of INaR reflects the slow accumulation of recovered/opened Nav channels into a second, alternative and more slowly populated, inactivated state. Additional simulations reveal that extrinsic factors that affect the kinetics of fast or slow Nav channel inactivation and/or impact the relative distribution of Nav channels in the fast- and slow-inactivated states, such as the accessory Nav{beta}4 channel subunit, can modulate the amplitude of INaR.

SUMMARYThe resurgent component of the voltage-gated sodium current (INaR) is revealed on membrane hyperpolarizations following brief depolarizing voltage steps that activate the rapidly activating and inactivating, transient Nav current (INaT). To probe the mechanisms contributing to the generation and properties of INaR, we combined whole-cell voltage-clamp recordings from mouse cerebellar Purkinje neurons with computational modeling to develop a novel, blocking particle-independent, model for the gating of INaR that involves two parallel inactivation pathways, and we show that this model recapitulates the detailed biophysical properties of INaR measured in mouse cerebellar Purkinje neurons.


title: Enteric pathogens induce tissue tolerance and prevent neuronal loss from subsequent infections

abstract: The enteric nervous system (ENS) controls several intestinal functions including motility and nutrient handling, which can be disrupted by infection-induced neuropathies or neuronal cell death. We investigated possible tolerance mechanisms preventing neuronal loss and disruption in gut motility after pathogen exposure. We found that following enteric infections, muscularis macrophages (MMs) acquire a tissue-protective phenotype that prevents neuronal loss and dysmotility during subsequent challenge with unrelated pathogens. Bacteria-induced neuroprotection relied on activation of gut-projecting sympathetic neurons and signaling via {beta}2-adrenergic receptors ({beta}2AR) on MMs. In contrast, helminth-mediated neuroprotection was dependent on T cells and systemic production of interleukin (IL)-4 and -13 by eosinophils, which induced arginase-expressing MMs that prevented neuronal loss from an unrelated infection located in a different intestinal region. Collectively, these data suggest that distinct enteric pathogens trigger a state of disease- or tissue tolerance that preserves ENS number and functionality.


title: Eosinophils are dispensable for development of MOG35-55-induced experimental autoimmune encephalomyelitis in mice

abstract: Experimental autoimmune encephalomyelitis (EAE) represents the mouse model of multiple sclerosis, a devastating neurological disorder. EAE development and progression involves the infiltration of different immune cells into the brain and spinal cord. However, less is known about a potential role of eosinophil granulocytes for EAE disease pathogenesis. In the present study, we found enhanced eosinophil abundance accompanied by increased concentration of the eosinophil chemoattractant eotaxin-1 in the spinal cord in the course of EAE induced in C57BL/6 mice by immunization with MOG35-55 peptide. However, the absence of eosinophils did not affect neuroinflammation, demyelination and clinical development or severity of EAE, as assessed in {Delta}dblGATA1 eosinophil-deficient mice. Taken together, despite their enhanced abundance in the inflamed spinal cord during disease progression, eosinophils were dispensable for EAE development.


title: Loss of Usp22 enhances histone H2B monoubiquitination and stimulates intracellular and systemic interferon immunity

abstract: Interferons protect from virus infections by inducing hundreds of interferon-stimulated genes (ISG) which orchestrate anti-viral adaptive and innate immunity. Upon viral infection or type I interferon (IFN) stimulation of cell lines, a histone modification, monoubiquitinated histone 2B (H2Bub1), increases at ISG loci, raising the possibility that a specific chromatin state can broadly stimulate IFN immunity in vivo. Here we show that, in the absence of virus infection or elevated interferon levels, mice lacking the relevant deubiquitinase, Usp22, in immune cells have elevated H2Bub1 levels. Hypermonoubiquitinated H2B is physically associated with dozens of ISG loci, and expression of large numbers of ISG is upregulated. This epigenetic state promotes intracellular and systemic immune phenotypes akin to adaptive and innate interferon immunity, and thereby identifies Usp22 as a negative regulator of interferon immunity.


title: Loss of TDP-43 function and rimmed vacuoles persist after T cell depletion in a xenograft model of sporadic inclusion body myositis

abstract:
Sporadic inclusion body myositis (IBM) is the most common acquired muscle disease in adults over age 50, yet it remains unclear whether the disease is primarily driven by T cell-mediated autoimmunity. IBM muscle biopsies exhibit nuclear clearance and cytoplasmic aggregation of TDP-43 in muscle cells, a pathologic finding observed initially in neurodegenerative disease, and nuclear loss of TDP-43 in neurons causes aberrant RNA splicing. Here, we show that loss of TDP-43 splicing repression, as determined by inclusion of cryptic exons, occurs in skeletal muscle of IBM patients. Out of 119 muscle biopsies tested, RT-PCR-mediated detection of cryptic exon expression is 84% sensitive and 99% specific for diagnosing IBM, indicating utility as a functional and diagnostic biomarker. To determine the role of T cells in pathogenesis, we generated a novel xenograft model by transplanting human IBM muscle into the hindlimb of immunodeficient mice. Xenografts from IBM patients display robust regeneration of human myofibers and recapitulate both inflammatory and degenerative features of the disease. Myofibers in IBM xenografts are invaded by human, oligoclonal CD8+ T cells and exhibit MHC-I upregulation, rimmed vacuoles, mitochondrial pathology, p62-positive inclusions, and nuclear clearance and cytoplasmic aggregation of TDP-43, resulting in expression of cryptic exons. Depletion of human T cells within IBM xenografts by treating mice intraperitoneally with anti-CD3 (OKT3) suppresses MHC-I upregulation, but rimmed vacuoles and loss of TDP-43 function persist. These data suggest that myofiber degeneration occurs independent of T cells, and muscle cell-intrinsic mechanisms, such as loss of TDP-43 splicing repression, drive IBM pathogenesis.

One Sentence SummaryDepletion of T cells in a xenograft model of sporadic inclusion body myositis suppresses inflammation but not TDP-43 pathology or muscle degeneration.

category: pathology


title: The Hepcidin Regulator Erythroferrone is a New Member of the Erythropoiesis-Iron-Bone Circuitry

abstract: Erythroblast erythroferrone (ERFE) secretion inhibits hepcidin expression by sequestering several bone morphogenetic protein (BMP) family members to increase iron availability for erythropoiesis. We report that ERFE expression in osteoblasts is higher compared with erythroblasts, is independent of erythropoietin, and functional in suppressing hepatocyte hepcidin expression. Erfe-/- mice display low-bone-mass arising from increased bone resorption despite a concomitant increase in bone formation. Consistently, Erfe-/- osteoblasts exhibit enhanced mineralization, Sost and Rankl expression, and BMP-mediated signaling ex vivo. The ERFE effect on osteoclasts is mediated through increased osteoblastic RANKL and sclerostin expression, increasing osteoclastogenesis in Erfe-/- mice. Importantly, Erfe loss in {beta}-thalassemic (Hbbth3/+) mice, a disease model with increased ERFE expression, triggers profound osteoclastic bone resorption and bone loss. Together, ERFE exerts an osteoprotective effect by modulating BMP signaling in osteoblasts, decreasing RANKL production to limit osteoclastogenesis, and prevents excessive bone loss during expanded erythropoiesis in {beta}-thalassemia.

category: physiology


title: Regulation of Plant Phototropic Growth by NPH3/RPT2-like Substrate Phosphorylation and 14-3-3 Binding

abstract: Polarity underlies all plant physiology and directional growth responses such as phototropism. Yet, our understanding of how plant tropic responses are established is far from complete. The plasma-membrane associated BTB-containing protein, NON-PHOTOTROPIC HYPOCOTYL 3 (NPH3) is a key determinant of phototropic growth which is regulated by AGC kinases known as the phototropins (phots). However, the mechanism by which phots initiate phototropic signalling via NPH3, and other NPH3/RPT2-like (NRL) members, has remained unresolved. Here we demonstrate that NPH3 is directly phosphorylated by phot1 both in vitro and in vivo. Light-dependent phosphorylation within a conserved consensus sequence (RxS) located at the extreme C-terminus of NPH3 is necessary to promote its functionality for phototropism and petiole positioning in Arabidopsis. Phosphorylation of this region by phot1 also triggers 14-3-3 binding combined with changes in NPH3 phosphorylation and localisation status. Seedlings expressing mutants of NPH3 that are unable to bind or constitutively bind 14-3-3s show compromised functionality that is consistent with a model where signalling outputs arising from a gradient in NPH3 RxS phosphorylation/localisation across the stem are a major contributor to phototropic responsiveness. Our current findings provide further evidence that 14-3-3 proteins are instrumental components regulating auxin-dependent growth and show for the first time that NRL proteins are direct phosphorylation targets for plant AGC kinases. Moreover, the C-terminal phosphorylation site/14-3-3-binding motif of NPH3 is conserved in several members of the NRL family, suggesting a common mechanism of regulation.

category: plant biology


title: Light-triggered and phosphorylation-dependent 14-3-3 association with NONPHOTOTROPIC HYPOCOTYL 3 is required for hypocotyl phototropism

abstract: NON-PHOTOTROPIC HYPOCOTYL 3 (NPH3) is a key component of the phototropic response, acting downstream of the primary photoreceptor phototropin and upstream of auxin redistribution. Despite the obvious physiological significance of the blue light-induced differential growth process, the molecular mode of NPH3 action is poorly understood. Light-triggered dephosphorylation of NPH3, however, is thought to constitute a major signaling event. Here, we show that NPH3 directly binds to polyacidic phospholipids via a polybasic motif in its C-terminal domain, allowing for plasma membrane association in darkness. We further demonstrate that blue light induces phosphorylation of a C-terminal 14-3-3 binding motif in NPH3. Subsequent binding of 14-3-3 to the phosphorylated NPH3 in turn is required for light-triggered release of NPH3 from the plasma membrane. In the cytosol, NPH3 undergoes a dynamic transition from a dilute to a condensed state. Intriguingly, the dephosphorylated state of the 14-3-3 binding site as well as NPH3 plasma membrane association are recoverable in darkness. Given that NPH3 variants constitutively localizing either to the plasma membrane or to cytosolic condensates are non-functional, the phototropin-triggered and 14-3-3 mediated dynamic change in the subcellular localization of NPH3 seems to be crucial for its function. Taken together, our data demonstrate a fundamental role for 14-3-3 members in regulating NPH3 localization and auxin-dependent phototropic responses.


title: Integrative systems biology approach identified crucial genes and transcription factors associated with gallbladder cancer pathogenesis

abstract: Gallbladder cancer (GBC) has a lower incidence rate among the population relative to other cancer types but majorly contributes to the total cancer cases of the biliary tract system. GBC is distinguished from other malignancies due to its high mortality, marked geographical variation and poor prognosis. To date no systemic targeted therapy is available for GBC. The main objective of this study is to determine the molecular signatures correlated with GBC development using integrative system level approaches. We performed analysis of publicly available transcriptomic data to identify differentially regulated genes and pathways. Co-expression network analysis and differential regulatory network analysis identified hub genes and hub transcription factors (TFs) associated with GBC pathogenesis and progression. We then assessed the epithelial-mesenchymal transition (EMT) status of the hub genes using a combination of three scoring methods. The hub genes such as; CDC6, MAPK15, CCNB2, BIRC7, L3MBTL1 identified are regulators of cell cycle components which suggests that cell cycle regulatory genes are significantly linked to GBC pathogenesis and progression.

category: systems biology


title: A complex of BRCA2 and PP2A-B56 is required for DNA repair by homologous recombination

abstract: Mutations in the tumour suppressor gene BRCA2 are associated with predisposition to breast and ovarian cancers. BRCA2 has a central role in maintaining genome integrity by facilitating the repair of toxic DNA double-strand breaks (DSBs) by homologous recombination (HR). BRCA2 acts by promoting RAD51 nucleoprotein filament formation on resected single-stranded DNA, but how BRCA2 activity is regulated during HR is not fully understood. Here, we delineate a pathway where ATM and ATR kinases phosphorylate a highly conserved region in BRCA2 in response to DSBs. These phosphorylations stimulate the binding of the protein phosphatase PP2A-B56 to BRCA2 through a conserved binding motif. We show that the phosphorylation-dependent formation of the BRCA2-PP2A-B56 complex is required for efficient RAD51 loading to sites of DNA damage and HR-mediated DNA repair. Moreover, we find that several cancer-associated mutations in BRCA2 deregulate the BRCA2-PP2A-B56 interaction and sensitize cells to PARP inhibition. Collectively, our work uncovers PP2A-B56 as a positive regulator of BRCA2 function in HR with clinical implications for BRCA2 and PP2A-B56 mutated cancers.


title: Direct and indirect cues can enable dual-adaptation, but through different learning processes

abstract:
Switching between motor tasks requires accurate adjustments for changes in dynamics (grasping a cup) or sensorimotor transformations (moving a computer mouse). Dual-adaptation studies have investigated how learning of context-dependent dynamics or transformations is enabled by sensory cues. However, certain cues, such as color, have shown mixed results. We propose that these mixed results may arise from two major classes of cues: "direct" cues, which are part of the dynamic state and "indirect" cues, which are not. We hypothesized that explicit strategies would primarily account for adaptation for an indirect color cue but would be limited to simple tasks while a direct visual separation cue would allow implicit adaptation regardless of task complexity. To test this idea, we investigated the relative contribution of implicit and explicit learning in relation to contextual cue type (colored or visually shifted workspace) and task complexity (one or eight targets) in a dual-adaptation task. We found that the visual workspace location cue enabled adaptation across conditions primarily through implicit adaptation. In contrast, we found that the color cue was largely ineffective for dual-adaptation, except in a small subset of participants who appeared to use explicit strategies. Our study suggests that the previously inconclusive role of color cues in dual-adaptation may be explained by differential contribution of explicit strategies across conditions.

New & NoteworthyWe present evidence that learning of context-dependent dynamics proceeds via different processes depending on the type of sensory cue used to signal the context. Visual workspace location enabled learning different dynamics implicitly, presumably because it directly enters the dynamic state estimate. In contrast, a color cue was only successful where learners were apparently able to leverage explicit strategies to account for changed dynamics. This suggests a unification for the previously inconclusive role of color cues.


title: Motor usage imprints microtubule stability in the shaft

abstract: Tubulin dimers assemble into dynamic microtubules which are used by molecular motors as tracks for intracellular transport. Organization and dynamics of the microtubule network is commonly thought to be regulated at the polymer ends, where tubulin-dimers can be added or removed. Here we show that molecular motors running on microtubules cause exchange of dimers along the shaft. These sites of dimer exchange act as rescue sites where depolymerising microtubules stop shrinking and start re-growing. Consequently, the average length of microtubules increases depending on how frequently they are used as motor tracks. An increase of motor activity densifies the cellular microtubule network and enhances cell polarity. Running motors leave marks in the shaft serving as traces of microtubule usage to organize the polarity landscape of the cell.


title: Single-cell transcriptomics reveals a conserved metaplasia program in pancreatic injury

abstract:
BACKGROUND & AIMSAcinar to ductal metaplasia (ADM) occurs in the pancreas in response to tissue injury and is a potential precursor for adenocarcinoma. The goal of these studies was to define the populations arising from ADM, the associated transcriptional changes, and markers of disease progression.

METHODSAcinar cells were lineage-traced with enhanced yellow fluorescent protein (EYFP) to follow their fate upon injury. Transcripts of over 13,000 EYFP+ cells were determined using single-cell RNA sequencing (scRNA-seq). Developmental trajectories were generated. Data were compared to gastric metaplasia, KrasG12D-induced neoplasia, and human pancreatitis. Results were confirmed by immunostaining and electron microscopy. KrasG12D was expressed in injury-induced ADM using several inducible Cre drivers. Surgical specimens of chronic pancreatitis from 15 patients were evaluated by immunostaining.

RESULTSscRNA-seq of ADM revealed emergence of a mucin/ductal population resembling gastric pyloric metaplasia. Lineage trajectories suggest that some pyloric metaplasia cells can generate tuft and enteroendocrine cells (EECs). Comparison to KrasG12D-induced ADM identifies populations associated with disease progression. Activation of KrasG12D expression in HNF1B+ or POU2F3+ ADM populations leads to neoplastic transformation and formation of MUC5AC+ gastric-pit-like cells. Human pancreatitis samples also harbor pyloric metaplasia with a similar transcriptional phenotype.

CONCLUSIONSUnder conditions of chronic injury, acinar cells undergo a pyloric-type metaplasia to mucinous progenitor-like populations, which seed disparate tuft cell and EEC lineages. ADM-derived EEC subtypes are diverse. KrasG12D expression is sufficient to drive neoplasia from injury-induced ADM and offers an alternative origin for tumorigenesis. This program is conserved in human pancreatitis, providing insight into early events in pancreas diseases.


title: TBX5 drives Aldh1a2 expression to regulate a RA-Hedgehog-Wnt gene regulatory network coordinating cardiopulmonary development

abstract:
The gene regulatory networks that coordinate the development of the cardiac and pulmonary systems are essential for terrestrial life but poorly understood. The T-box transcription factor Tbx5 is critical for both pulmonary specification and heart development, but how these activities are mechanistically integrated remains unclear. We show that Tbx5 regulates an evolutionarily conserved retinoic acid (RA)-Hedgehog-Wnt signaling cascade coordinating cardiopulmonary development. We demonstrate that Tbx5 directly maintains expression of the RA-synthesizing enzyme Aldh1a2 in the foregut lateral plate mesoderm via an intronic enhancer that is evolutionarily conserved among terrestrial vertebrates. Tbx5 promotes posterior second heart field identity in a positive feedback loop with RA, antagonizing a Fgf8-Cyp regulatory module and restricting FGF activity to the anterior. Tbx5/Aldh1a2-dependent RA signaling also directly activates Shh transcription in the adjacent foregut endoderm through the conserved MACS1 enhancer. Epithelial Hedgehog then signals back to the mesoderm, where together with Tbx5 it activates expression of Wnt2/2b that ultimately induce pulmonary fate in the foregut endoderm. These results provide mechanistic insight into the interrelationship between heart and lung development informing cardiopulmonary evolution and birth defects.

KEY FINDINGSO_LITbx5 regulates second heart field patterning and pulmonary development via retinoic acid (RA) and Hedgehog (Hh) signaling.
C_LIO_LITbx5 directly maintains transcription of the RA-synthesizing enzyme Aldh1a2 in the posterior second heart field mesoderm via an evolutionarily conserved intronic enhancer.
C_LIO_LIDownstream of Tbx5, RA directly promotes Shh transcription through the evolutionarily conserved MACS1 endoderm enhancer.
C_LIO_LIDownstream of Tbx5, RA suppresses FGF signaling to pattern the second heart field while promoting a Hedgehog-Wnt2/2b signaling cascade that induces pulmonary fate.
C_LI

SUMMARY STATEMENTTbx5-dependent Retinoic Acid signaling regulates an evolutionarily conserved gene regulatory network that coordinates cardiac and pulmonary development.

category: developmental biology


title: Cardiac Differentiation of Human Pluripotent Stem Cells Using Defined Extracellular Matrix Proteins Reveals Essential Role of Fibronectin

abstract: Research and therapeutic applications using human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) require robust differentiation strategies. Efforts to improve hPSC-CM differentiation have largely overlooked the role of extracellular matrix (ECM). The present study investigates the ability of defined ECM proteins to promote hPSC cardiac differentiation. Fibronectin, laminin-111, and laminin-521 enabled hPSCs to attach and expand; however, fibronectin ECM either endogenously produced or exogenously added promoted, while laminins inhibited, cardiac differentiation in response to growth factors Activin A, BMP4, and bFGF. Inducible shRNA knockdown of fibronectin prevented Brachyury+ mesoderm formation and subsequent hPSC-CM differentiation. Antibodies blocking fibronectin binding to integrin {beta}1, but not 5, inhibited cardiac differentiation. Furthermore, inhibition of integrin-linked kinase blocked cardiac differentiation. These results identify fibronectin, laminin-111 and laminin-521 as defined substrates enabling cardiac differentiation of hPSCs and uncover the essential role of fibronectin and downstream signaling pathways in the early stage of hPSC-CM differentiation.


title: Insm1-expressing neurons and secretory cells develop from a common pool of progenitors in the sea anemone Nematostella vectensis

abstract: Neurons are highly specialized cells present in nearly all animals, but their evolutionary origin and relationship to other cell types is not well understood. We use here the sea anemone Nematostella vectensis as a model system for early-branching animals to gain fresh insights into the evolutionary history of neurons. We generated a transgenic reporter line to show that the transcription factor NvInsm1 is expressed in post-mitotic cells that give rise to various types of neurons and secretory cells. Expression analyses, double transgenics and gene knockdown experiments show that the NvInsm1-expressing neurons and secretory cells derive from a common pool of NvSoxB(2)-positive progenitor cells. These findings, together with the requirement for Insm1 for the development of neurons and endocrine cells in vertebrates, support a close evolutionary relationship of neurons and secretory cells.


title: Abundance-occupancy relationships along taxonomic ranks reveal a consistency of niche differentiation in marine bacterioplankton with distinct lifestyles

abstract:
Abundance-occupancy relationships (AORs) are an important determinant of biotic community dynamics and habitat suitability. However, little is known about their role in complex bacterial communities, either within a phylogenetic framework or as a function of niche breadth. Based on data obtained in a field study in the St. Lawrence Estuary, we used 16S rRNA gene sequencing to examine the vertical patterns, strength, and character of AORs for particle-attached and free-living bacterial assemblages. Free-living communities were phylogenetically more diverse than particle-attached communities. The dominant taxa were consistent in terms of their presence/absence but population abundances differed in surface water vs. the cold intermediate layer. Significant, positive AORs characterized all of the surveyed communities across all taxonomic ranks of bacteria, thus demonstrating an ecologically conserved trend for both free-living and particle-attached bacteria. The strength of the AORs was low at the species level but higher at the genus and phylum levels. These results demonstrate that an assessment of the distributions and population densities of finely resolved taxa does not necessarily improve determinations of apparent niche differences in marine bacterioplankton communities at regional scales compared with the information inferred from a broad taxonomic classification.

Subject Categorymicrobial population and community ecology

category: ecology


title: Double-observer approach with camera traps: Towards an unbiased density estimation of unmarked animal populations

abstract: Camera traps are a powerful tool for wildlife surveys. However, camera traps may not always detect animals passing in front. This constraint may create a substantial bias in estimating critical parameters such as the density of unmarked populations. We proposed the  double-observer approach with camera traps to counter the constraint, which involves setting up a paired camera trap at a station and correcting imperfect detection with a hierarchal capture-recapture model for stratified populations. We performed simulations to evaluate this approachs reliability and determine how to obtain desirable data for this approach. We then applied it to 12 mammals in Japan and Cameroon. The results showed that the approach could correct imperfect detection as long as paired camera traps detect animals nearly independently (Correlation coefficient < 0.2). Camera traps should be installed to monitor a predefined small focal area from different directions to satisfy this requirement. The field surveys showed that camera trap might miss animals by 3 %-40%, suggesting that current density estimation models relying on perfect detection may underestimate animal density by the same order of magnitude. We hope that our approach will be incorporated into existing density estimation models to improve their accuracy.


title: Improving estimation of species distribution from citizen-science records using data-integration models

abstract: Mapping species distributions is a crucial but challenging requirement of wildlife management. The frequent need to sample vast expanses of potential habitat increases the cost of planned surveys and rewards accumulation of opportunistic observations. In this paper, we integrate planned survey data from roost counts with opportunistic samples from eBird, WikiAves and Xeno-canto citizen-science platforms to map the geographic range of the endangered Vinaceous-breasted Parrot. We demonstrate the estimation and mapping of species occurrence based on data integration while accounting for specifics of each data set, including observation technique and uncertainty about the observations. Our analysis illustrates 1) the incorporation of sampling effort, spatial autocorrelation, and site covariates in a joint-likelihood, hierarchical, data-integration model; 2) the evaluation of the contribution of each data set, as well as the contribution of effort covariates, spatial autocorrelation, and site covariates to the predictive ability of fitted models using a cross-validation approach; and 3) how spatial representation of the latent occupancy state (i.e. realized occupancy) helps identify areas with high uncertainty that should be prioritized in future field work. Our results reveal a Vinaceous-breasted Parrot geographic range of 434,670 km2, which is three times larger than the  Extant area previously reported in the IUCN Red List. The exclusion of one data set at a time from the analyses always resulted in worse predictions by the models of truncated data than by the full model, which included all data sets. Likewise, exclusion of spatial autocorrelation, site covariates, or sampling effort resulted in worse predictions. The integration of different data sets into one joint-likelihood model produced a more reliable representation of the species range than any individual data set taken on its own improving the use of citizen science data in combination with planned survey results.


title: The dual nature of metacommunity variability

abstract: There is increasing interest in measuring ecological stability to understand how communities and ecosystems respond to broad-scale global changes. One of the most common approaches is to quantify the variation through time in community or ecosystem aggregate attributes (e.g., total biomass), referred to as aggregate variability. It is now widely recognized that aggregate variability represents only one aspect of communities and ecosystems, and compositional variability, the changes in the relative frequency of species in an assemblage, is equally important. Recent contributions have also begun to explore ecological stability at regional spatial scales, where interconnected local communities form metacommunities, a key concept in managing complex landscapes. However, the conceptual frameworks and measures of ecological stability in space have only focused on aggregate variability, leaving a conceptual gap. Here, we address this gap with a novel framework for quantifying the aggregate and compositional variability of communities and ecosystems through space and time. We demonstrate that the compositional variability of a metacommunity depends on the degree of spatial synchrony in compositional trajectories among local communities. We then provide a conceptual framework in which compositional variability of (i) the metacommunity through time and (ii) among local communities combine into four archetype scenarios: spatial stasis (low/low); spatial synchrony (high/low); spatial asynchrony (high/high) and spatial compensation (low/high). We illustrate this framework based on numerical examples and a case study of a macroalgal metacommunity in which low spatial synchrony reduced variability in aggregate biomass at the metacommunity scale, while masking high spatial synchrony in compositional trajectories among local communities. Finally, we discuss the role of dispersal, environmental heterogeneity, species interactions and suggest future avenues. We believe this framework will be helpful for considering both aspects of variability simultaneously which is important to better understand ecological stability in natural and complex landscapes in response to environmental changes.


title: Predation by protists influences the temperature response of microbial communities

abstract: Temperature strongly influences microbial community structure and function, which in turn contributes to the global carbon cycle that can fuel further warming. Recent studies suggest that biotic interactions amongst microbes may play an important role in determining the temperature responses of these communities. However, how microbial predation regulates these communities under future climates is still poorly understood. Here we assess whether predation by one of the most important bacterial consumers globally - protists - influences the temperature response of a freshwater microbial community structure and function. To do so, we exposed these microbial communities to two cosmopolitan species of protists at two different temperatures, in a month-long microcosm experiment. While microbial biomass and respiration increased with temperature due to shifts in microbial community structure, these responses changed over time and in the presence of protist predators. Protists influenced microbial biomass and function through effects on community structure, and predation actually reduced microbial respiration rate at elevated temperature. Indicator species and threshold indicator taxa analyses showed that these predation effects were mostly determined by phylum-specific bacterial responses to protist density and cell size. Our study supports previous findings that temperature is an important driver of microbial communities, but also demonstrates that predation can mediate these responses to warming, with important consequences for the global carbon cycle and future warming.


title: Correlated evolution of parental care with dichromatism, colour, and patterns in anurans

abstract: Parental care is widespread among vertebrates, with clear fitness benefits. Caring parents however incur costs that include higher predation risk. Anurans have among the most diverse forms of parental care, and we test whether the occurrence of care is associated with morphology that minimizes predation risk. We first examine whether parental care co-occurs with sexual dichromatism, testing the hypothesis that when one sex is conspicuous, the other is cryptically patterned and cares for the young. From our phylogenetic comparative analyses of 988 anurans distributed globally, we find that parental care is less likely to co-occur with dichromatism, irrespective of the caring sex. We then examine whether colour gradients and patterns that enhance crypticity are associated with the occurrence of parental care. We found that species with male-only care were more likely to have Bars-Bands, but contrary to our expectation, other colours (Green-Brown, Red-Blue-Black, Yellow) and patterns (Plain, Spots, Mottled-Patches) were not associated with caregiving behaviours. The lack of strong correlations between dorsal morphology and parental care suggests that crypticity is not the dominant strategy to minimise predation risk for care-giving anurans, and that the evolution of body colour and parental care are driven by independent selection pressures.

category: evolutionary biology


title: Persistence and expansion of cryptic endangered red wolf genomic ancestry along the American Gulf coast

abstract: Admixture and introgression play a critical role in adaptation and genetic rescue that has only recently gained a deeper appreciation. Here, we explored the geographic and genomic landscape of cryptic ancestry of the endangered red wolf that persists within the genome of a ubiquitous sister taxon, the coyote, all the while the red wolf has been extinct in the wild since the early 1980s. We assessed admixture across 102,621 SNP loci genotyped in 293 canid genomes. We found support for increased red wolf ancestry along an east-to-west gradient across the southern United States that was associated with historical admixture in the past 100 years. Southwestern Louisiana and southeastern Texas, the geographic zone where the last red wolves were known prior to their extinction in the wild, contained the highest and oldest levels of red wolf ancestry. X-linked regions of low recombination rates were depleted of introgression, relative to the autosomes, suggestive of the large X effect and enrichment with loci involved in maintaining reproductive isolation. Recombination rate was positively correlated with red wolf ancestry across coyote genomes, consistent with theoretical predictions. The geographic and genomic extent of cryptic red wolf ancestry can provide novel and variable genomic resources for the survival of the endangered red wolf.


title: The unresolved phylogenomic tree of butterflies and moths (Lepidoptera): assessing the potential causes and consequences

abstract: The field of molecular phylogenetics is being revolutionised with next-generation sequencing technologies making it possible to sequence large numbers of genomes for non-model organisms ushering us into the era of phylogenomics. The current challenge is no longer how to get enough data, but rather how to analyse the data and how to assess the support for the inferred phylogeny. We focus on one of the largest animal groups on the planet - butterflies and moths (order Lepidoptera). We clearly demonstrate that there are unresolved issues in the inferred phylogenetic relationships of the major lineages, despite several recent phylogenomic studies of the group. We assess the potential causes and consequences of the conflicting phylogenetic hypotheses. With a dataset consisting of 331 protein-coding genes and the alignment length over 290 000 base pairs, including 200 taxa representing 81% of lepidopteran superfamilies, we compare phylogenetic hypotheses inferred from amino acid and nucleotide alignments. The resulting two phylogenies are discordant, especially with respect to the placement of the superfamily Gelechioidea, which is likely due to compositional bias of both the nucleotide and amino acid sequences. With a series of analyses, we dissect our dataset and demonstrate that there is sufficient phylogenetic signal to resolve much of the lepidopteran tree of life. Overall, the results from the nucleotide alignment are more robust to the various perturbations of the data that we carried out. However, the lack of support for much of the backbone within Ditrysia makes the current butterfly and moth tree of life still unresolved. We conclude that taxon sampling remains an issue even in phylogenomic analyses, and recommend that poorly sampled highly diverse groups, such as Gelechioidea in Lepidoptera, should receive extra attention in the future.


title: Extremely rapid development and early reproduction reduces sexual selection in D. melanogaster populations.

abstract:
The role of sexual selection in mediating levels of sexual conflict has been demonstrated in many experimental evolution studies on Drosophila spp. where competition among males for mating was the target of selection. Sexual selection has also been shown to affect the evolution of life-histories. However, the influence of divergent life-histories on reproductive strategies and, therefore, sexual selection and possibly sexual conflict, has been less well studied. We examined D. melanogaster populations selected for a short development time and early age at reproduction for changes in reproductive behaviour and traits that are proxies of sexual selection. We report a large reduction in reproductive competition experienced by the males of these populations, compared to ancestral populations that are not consciously selected for rapid development or early reproduction, potentially leading to reduced sexual selection. We show that rapidly developing and early reproducing populations have very low levels of mating in their lifetime (females are more or less monandrous), low courtship levels, shorter copulation duration, and longer time from eclosion to first mating, compared to the controls. These results are discussed in the context of the previously demonstrated reduction of inter-locus sexual conflict in these populations. We show that life-history strategies might have a large and significant impact on sexual selection, with each influencing the other and contributing to the complexities of adaptation.

Significance statementSexual conflict, often manifested as an arms-race between males and females trying to enhance their own reproductive success at some cost to the other, is of great evolutionary interest because it can maintain genetic variation in populations, prevent the independent optimization of male and female traits, and also promote speciation. Sexual selection, or variation in mating success, is well known to affect levels of sexual conflict. However, it is not so clear whether, and how, the regular evolution of life-histories also affects sexual selection. Here, we show that life-history evolution in fruit fly populations selected for traits not directly related to sexual conflict might, nevertheless, mediate the possible evolution of altered sexual conflict levels through effects on sexual selection. Populations that evolved to develop to adulthood fast, and reproduce relatively early in life, are shown to potentially experience less sexual selection, which can explain the low sexual conflict levels earlier observed in them.


title: Shorter effective lifespan in laboratory populations of D. melanogaster can reduce sexual selection.


title: Shorter effective lifespan in laboratory populations of D. melanogaster might reduce sexual selection.


title: Flashfm: A Flexible and Shared Information Fine-mapping Approach for Multiple Quantitative Traits

abstract: Joint fine-mapping that leverages information between quantitative traits could improve accuracy and resolution over single-trait fine-mapping. Using summary statistics, flashfm (FLexible And SHared information Fine-Mapping) fine-maps signals for multiple traits, allowing for missing trait measurements and use of related individuals. In a Bayesian framework, prior model probabilities are formulated to favour model combinations that share causal variants to capitalise on information between traits. Simulation studies demonstrate that both approaches produce broadly equivalent results when traits have no shared causal variants. When traits share at least one causal variant, flashfm reduces the number of potential causal variants by 30% compared with single-trait fine-mapping. In a Ugandan cohort with 33 cardiometabolic traits, flashfm gave a 20% reduction in the total number of potential causal variants from single-trait fine-mapping. Flashfm is computationally efficient and can easily be deployed across publicly available summary statistics for signals in up to six traits.


title: Stop Meta-Analyzing, Start Instrumenting: Maximizing the Predictive Power of Polygenic Scores

abstract: The conventional way of boosting the predictive power of polygenic scores is to increase the GWAS sample size by meta-analyzing GWAS results of multiple cohorts. In this paper, we challenge this convention. Through simulations, we show that Instrumental Variable (IV) regression using two polygenic scores constructed from independent GWAS summary statistics outperforms the typical Ordinary Least Squares (OLS) model employing a single meta-analysis based polygenic score in terms of bias, root mean squared error, and statistical power. We verify the empirical validity of the simulations by predicting educational attainment (EA) and height in a sample of siblings from the UK Biobank. We show that IV regression between-families approaches the SNP-based heritability, and improves the predictive power of polygenic scores by 12% (height) to 22% (EA). Furthermore, IV regression within-families provides the tightest lower bound for the direct genetic effect, increasing the lower bound for EA from 2.0% to 3.4%, and for height from 28.9% to 37.7%.


title: Using Obviously-Related Instrumental Variables to Increase the Predictive Power of Polygenic Scores


title: What is a supercoiling-sensitive gene? Insights from topoisomerase I inhibition in the Gram-negative bacterium Dickeya dadantii

abstract: DNA supercoiling is an essential mechanism of bacterial chromosome compaction, whose level is mainly regulated by topoisomerase I and DNA gyrase. Inhibiting either of these enzymes with antibiotics leads to global supercoiling modifications and subsequent changes in global gene expression. In previous studies, genes responding to DNA relaxation induced by gyrase inhibition were categorised as "supercoiling-sensitive". Here, we studied the opposite variation of DNA supercoiling in the phytopathogen Dickeya dadantii using the non-marketed antibiotic seconeolitsine. We showed that the drug is active against topoisomerase I from this species, and analysed the first transcriptomic response of a Gram-negative bacterium to topoisomerase I inhibition. We find that the responding genes essentially differ from those observed after DNA relaxation, and further depend on the growth phase. We characterised these genes at the functional level, and also detected distinct patterns in terms of expression level, spatial and orientational organisation along the chromosome. Altogether, these results highlight that the supercoiling-sensitivity is a complex feature, which depends on the action of specific topoisomerases, on the physiological conditions, and on their genomic context. Based on previous in vitro expression data of several promoters, we propose a qualitative model of SC-dependent regulation that accounts for many of the contrasting transcriptomic features observed after gyrase or topoisomerase I inhibition.

category: genomics


title: In-depth phylogenomic analysis of arbuscular mycorrhizal fungi based on a comprehensive set of de novo genome assemblies

abstract:
O_LIMorphological characters and nuclear ribosomal DNA (rDNA) phylogenies have so far been the basis of the current classifications of arbuscular mycorrhizal (AM) fungi. Improved understanding of the phylogeny and evolutionary history of AM fungi requires extensive ortholog sampling and analyses of genome and transcriptome data from a wide range of taxa.
C_LIO_LITo circumvent the need for axenic culturing of AM fungi we gathered and combined genomic data from single nuclei to generate de novo genome assemblies covering seven families of AM fungi. Comparative analysis of the previously published Rhizophagus irregularis DAOM197198 assembly confirm that our novel workflow generates high-quality genome assemblies suitable for phylogenomic analysis. Predicted genes of our assemblies, together with published protein sequences of AM fungi and their sister clades, were used for phylogenomic analyses.
C_LIO_LIBased on analyses of sets of orthologous genes, we highlight three alternative topologies among families of AM fungi. In the main topology, Glomerales is polyphyletic and Claroideoglomeraceae, is the basal sister group to Glomeraceae and Diversisporales.
C_LIO_LIOur results support family level classification from previous phylogenetic studies. New evolutionary relationships among families where highlighted with phylogenomic analysis using the hitherto most extensive taxon sampling for AM fungi.
C_LI


title: High-speed widefield handheld swept-source OCT angiography with a VCSEL light source

abstract: Optical coherence tomography (OCT) and OCT angiography (OCTA) enable noninvasive structural and angiographic imaging of the eye. Portable handheld OCT/OCTA systems are required for imaging patients in the supine position. Examples include infants in the neonatal intensive care unit (NICU) and operating room (OR). The speed of image acquisition plays a pivotal role in acquiring high quality OCT/OCTA images, particularly with the handheld system, since both the operator hand tremor and subject motion can cause significant motion artifacts. In addition, having a large field of view and the ability of real-time data visualization are critical elements in rapid disease screening, reducing imaging time, and detecting peripheral retinal pathologies. The arrangement of optical components is less flexible in the handheld system due to the limitation of size and weight. In this paper, we introduce a 400-kHz, 55-degree field of view handheld OCT/OCTA system that has overcome many technical challenges as a portable OCT system as well as a high-speed OCTA system. We demonstrate imaging premature infants with retinopathy of prematurity (ROP) in the NICU, and patients with incontinentia pigmenti (IP) in the OR using our handheld OCT system. Our design may have potential for improving the diagnosis of retinal diseases and help provide a practical guideline for designing a flexible and portable OCT system.

category: bioengineering


title: Head model personalization: A framework for morphing lifespan brain images and brains with substantial anatomical changes

abstract: Finite element (FE) head models have emerged as a powerful tool in many fields within neuroscience, especially for studying the biomechanics of traumatic brain injury (TBI). Personalized head models are needed to account for geometric variations among subjects for more reliable predictions. However, the generation of subject-specific head models with conforming hexahedral elements suitable for studying the biomechanics of TBIs remains a significant challenge, which has been a bottleneck hindering personalized simulations. This study presents a framework capable of generating lifespan brain models and pathological brains with substantial anatomical changes, morphed from a previously developed baseline model. The framework combines hierarchical multiple feature and multimodality imaging registrations with mesh grouping, which is shown to be efficient with a heterogeneous dataset of seven brains, including a newborn, 1-year-old (1Y), 2Y, 6Y, adult, 92Y, and a hydrocephalus brain. The personalized models of the seven subjects show competitive registration accuracy, demonstrating the potential of the framework for generating personalized models for almost any brains with substantial anatomical changes. The family of head injury models generated in this study opens vast opportunities for studying age-dependent and groupwise brain injury mechanisms. The framework is equally applicable for personalizing head models in other fields, e.g., in tDCS, TMS, TUS, as an efficient approach for generating subject-specific head models than from scratch.


title: Extracellular vesicles Transfer Polarized Mitochondria and Increase Cellular Energetics in Ischemic Endothelial Cells

abstract:
We have demonstrated, for the first time that microvesicles, a sub-type of extracellular vesicles (EVs) derived from hCMEC/D3: a human brain endothelial cell (BEC) line transfer polarized mitochondria to recipient BECs in culture and to neurons in mice acute brain cortical and hippocampal slices. This mitochondrial transfer increased ATP levels by 100 to 200-fold (relative to untreated cells) in the recipient BECs exposed to oxygen-glucose deprivation, an in vitro model of cerebral ischemia. We have also demonstrated that transfer of microvesicles, the larger EV fraction, but not exosomes resulted in increased mitochondrial function in hypoxic endothelial cultures. Gene ontology and pathway enrichment analysis of EVs revealed a very high association to glycolysis-related processes. In comparison to heterotypic macrophage- derived EVs, BEC-derived EVs demonstrated a greater selectivity to transfer mitochondria and increase endothelial cell survival under ischemic conditions.

HighlightsO_LIMicrovesicles transfer mitochondria to endothelial cells and brain slice neurons
C_LIO_LIMitochondrial transfer increased ATP in ischemic brain endothelial cells (BECs)
C_LIO_LITransfer of microvesicles increased mitochondrial function in hypoxic BECs
C_LIO_LITransfer of exosomes did not affect mitochondrial function in hypoxic BECs
C_LIO_LIHomotypic BEC-derived EVs result in greater ATP levels in the recipient BECs
C_LI


title: Microvesicles Transfer Mitochondria and Increase Mitochondrial Function in Brain Endothelial Cells


title: An Engineered Contact Lens for Passive and Sustained Release of Lifitegrast, an Anti-Dry Eye Syndrome Drug

abstract:
Lifitegrast is an FDA-approved drug that inhibits T-cell mediated inflammation associated with dry eye syndrome (DES). Lifitegrast is a potent inhibitor of the interaction between LFA-1 on T-cells and ICAM-1 on endothelial cells at the ocular surface. While effective in treating DES, 5% (81.2 mM) lifitegrast has low drug utilization and elicits off-target effects. Here we engineer contact lenses to release therapeutically-relevant doses of lifitegrast to every tear film for up to 10-hours. Lifitegrast is coupled to the polymer of the soft hydrogel lens via a photolabile (caged) crosslinker. Exposures of the lens to the 400-430 nm wavelengths of indoor daylight excite the caged crosslinker molecules and trigger a bond-cleavage reaction that releases authentic lifitegrast passively to the tear film. The photoproduct of the reaction remains chemically-linked to the polymer of the single-use lens. Our studies show that passive exposures of the lens to indoor light would generate an average of 990 nM lifitegrast to every tear film in a zero-order reaction for up to 10-hours. This concentration exceeds the Kd for the interaction between ICAM-1 and LFA-1 by [~]330-fold and would sustain inhibition of inflammatory responses at the ocular surface. The amount of lifitegrast released from the lens increases during exposures to outdoor sunlight. Over a 10-hour exposure to indoor light, a single lens would release 0.44% of the lifitegrast present in two drops of commercial 5% lifitegrast. Compared to tear-drop approaches, our engineered lenses would sustain the passive delivery of therapeutically-relevant doses of lifitegrast over a longer period, and exhibit improved drug utilization at a lower cost. Our technology could easily be integrated into daily-use contact lenses in order to prevent inflammation at the ocular surface, dry-eye and contact lens-mediated discomfort.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=198 HEIGHT=200 SRC="FIGDIR/small/439289v1_ufig1.gif" ALT="Figure 1">
View larger version (75K):
org.highwire.dtl.DTLVardef@1049076org.highwire.dtl.DTLVardef@a83fceorg.highwire.dtl.DTLVardef@1f52c4borg.highwire.dtl.DTLVardef@12f0ead_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Improved data sets and evaluation methods for the automatic prediction of DNA-binding proteins

abstract:
MotivationAccurate automatic annotation of protein function relies on both innovative models and robust datasets. Due to their importance in biological processes, the identification of DNA-binding proteins directly from protein sequence has been the focus of many studies. However, the data sets used to train and evaluate these methods have suffered from substantial flaws. We describe some of the weaknesses of the data sets used in previous DNA-binding protein literature and provide several new data sets addressing these problems. We suggest new evaluative benchmark tasks that more realistically assess real-world performance for protein annotation models. We propose a simple new model for the prediction of DNA-binding proteins and compare its performance on the improved data sets to two previously published models. Additionally, we provide extensive tests showing how the best models predict across taxonomies.

ResultsOur new gradient boosting model, which uses features derived from a published protein language model, outperforms the earlier models. Perhaps surprisingly, so does a baseline nearest neighbor model using BLAST percent identity. We evaluate the sensitivity of these models to perturbations of DNA-binding regions and control regions of protein sequences. The successful data-driven models learn to focus on DNA-binding regions. When predicting across taxonomies, the best models are highly accurate across species in the same kingdom and can provide some information when predicting across kingdoms.

Code and Data AvailabilityAll the code and data for this paper can be found at https://github.com/AZaitzeff/tools_for_dna_binding_proteins.

Contactalexander.zaitzeff@twosixtech.com

category: bioinformatics


title: Biotext: Exploiting Biological-like format for text mining

abstract: The large amount of existing textual data justifies the development of new text mining tools. Bioinformatics tools can be brought to Text Mining, increasing the arsenal of resources. Here, we present Biotext, a package of strategies for converting natural language text into biological-like information data, providing a general protocol with standardized functions, allowing to share, encode and decode textual data for amino acid data. The package was used to encode the arbitrary information present in the headings of the biological sequences found in a BLAST survey. The protocol implemented in this study consists of 12 steps, which can be easily executed and/ or changed by the user, depending on the study area. Biotext empowers user to perform text mining using bioinformatics tools. Biotext is Freely available at https://pypi.org/project/biotext/ (Python package) and https://sourceforge.net/projects/biotext-tools/files/AMINOcode_GUI/ (Standalone tool).


title: DDN2.0: R and Python packages for differential dependency network analysis of biological systems

abstract:
Data-driven differential dependency network analysis identifies in a complex and often unknown overall molecular circuitry a network of differentially connected molecular entities (pairwise selective coupling or uncoupling depending on the specific phenotypes or experimental conditions) (Herrington, et al. 2018; Zhang, et al., 2009; Zhang and Wang, 2010; Zhang, et al., 2016). Such differential dependency networks are typically used to assist in the inference of potential key pathways. Based on our previously developed Differential Dependency Network (DDN) method, we report here the fully implemented R and Python software tool packages for public use. The DDN2.0 algorithm uses a fused Lasso model and block-wise coordinate descent to estimate both the common and differential edges of dependency networks. The identified DDN can help to provide plausible interpretation of data, gain new insight of disease biology, and generate novel hypotheses for further validation and investigations.

To address the imbalanced sample group problem, we propose a sample-size normalized formulation to correct systematic bias. To address high computational complexity, we propose four strategies to accelerate DDN2.0 learning. The experimental results show that new DDN2.0+ learning speed with combined four accelerating strategies is hundreds of times faster than that of DDN2.0 algorithm on medium-sized data (Fu, 2019). To detect intra-omics and inter-omics network rewiring, we propose multiDDN using a multi-layer signaling model to integrate multi-omics data. The simulation study shows that the multiDDN method can achieve higher accuracy of detecting network rewiring (Fu, 2019).


title: Nonnegative matrix factorization integrates single-cell multi-omic datasets with partially overlapping features

abstract: Single-cell genomic technologies provide an unprecedented opportunity to define molecular cell types in a data-driven fashion, but present unique data integration challenges. Integration analyses often involve datasets with partially overlapping features, including both shared features that occur in all datasets and features exclusive to a single experiment. Previous computational integration approaches require that the input matrices share the same number of either genes or cells, and thus can use only shared features. To address this limitation, we derive a novel nonnegative matrix factorization algorithm for integrating single-cell datasets containing both shared and unshared features. The key advance is incorporating an additional metagene matrix that allows unshared features to inform the factorization. We demonstrate that incorporating unshared features significantly improves integration of single-cell RNA-seq, spatial transcriptomic, SHARE-seq, and cross-species datasets. We have incorporated the UINMF algorithm into the open-source LIGER R package (https://github.com/welch-lab/liger).


title: PanClassif: Improving pan cancer classification of single cell RNA-seq using machine learning

abstract: Cancer is one of the major causes of human death per year. In recent years, cancer identification and classification using machine learning have gained momentum due to the availability of high throughput sequencing data. Using RNA-seq, cancer research is blooming day by day and new insights of cancer and related treatments are coming into light. In this paper, we propose PanClassif, a method that requires a very few and effective genes to detect cancer from RNA-seq data and is able to provide performance gain in several wide range machine learning classifiers. We have taken 22 types of cancer samples from The Cancer Genome Atlas (TCGA) having 8287 cancer samples and 680 normal samples. Firstly, PanClassif uses k-Nearest Neighbor (k-NN) smoothing to smooth the samples to handle noise in the data. Then effective genes are selected by Anova based test. For balancing the train data, PanClassif applies an oversampling method, SMOTE. We have performed comprehensive experiments on the datasets using several classification algorithms. Experimental results shows that PanClassif out-perform existing state-of-the-art methods available and shows consistent performance for two single cell RNA-seq datasets taken from Gene Expression Omnibus (GEO). PanClassif improves performances of a wide variety of classifiers for both binary cancer prediction and multi-class cancer classification. PanClassif is available as a python package (https://pypi.org/project/panclassif/). All the source code and materials of PanClassif are available at https://github.com/Zwei-inc/panclassif.


title: Integrating Rigidity Analysis into the Exploration of Protein Conformational Pathways using RRT* and MC

abstract: To understand how proteins function on a cellular level, it is of paramount importance to understand their structures and dynamics, including the conformational changes they undergo to carry out their function. For the aforementioned reasons, the study of large conformational changes in proteins has been an interest to researchers for years. However, since some proteins experience rapid and transient conformational changes, it is hard to experimentally capture the intermediate structures. Additionally, computational brute force methods are computationally intractable, which makes it impossible to find these pathways which require a search in a high-dimensional, complex space. In our previous work, we implemented a hybrid algorithm that combines Monte-Carlo (MC) sampling and RRT*, a version of the Rapidly Exploring Random Trees (RRT) robotics-based method, to make the conformational exploration more accurate and efficient, and produce smooth conformational pathways. In this work, we integrated the rigidity analysis of proteins into our algorithm to guide the search to explore flexible regions. We demonstrate that rigidity analysis dramatically reduces the run time and accelerates convergence.


title: MetaRNN: Differentiating Rare Pathogenic and Rare Benign Missense SNVs and InDels Using Deep Learning

abstract: With advances in high-throughput DNA sequencing, numerous genetic variants have been discovered in the human genome. One challenge we face is interpreting these variants to help in disease screening, diagnosis, and treatment. While multiple computational approaches have been proposed to improve our understanding of genetic variants, their ability to identify rare pathogenic variants from rare benign ones is still lacking. Using context annotations and deep learning methods, we present pathogenicity prediction models, MetaRNN and MetaRNN-indel, to help identify and prioritize rare non-synonymous single nucleotide variants (nsSNVs) and non-frameshift insertion/deletions (nfINDELs). A recurrent neural network incorporating a +/- 1 codon window around the affected codon was combined with 28 high-level annotation scores and allele frequency features to develop the two proposed models. We use independent test datasets to demonstrate that these new models outperform state-of-the-art competitors and achieve a more interpretable score distribution. Importantly, prediction scores from the nsSNV-based and the nfINDEL-based models are comparable, enabling easy adoption of integrated genotype-phenotype association analysis methods. In addition, we provide pre-computed MetaRNN scores for all possible human nsSNVs and a Linux executable file for a fast one-stop annotation of nsSNVs and nfINDELs. All the resources are available at http://www.liulab.science/MetaRNN.


title: TRIMER: Transcription Regulation Integrated with MEtabolic Regulation

abstract: Advances in bioengineering have enabled numerous bio-based commodities. Yet most traditional approaches do not extend beyond a single metabolic pathway and do not attempt to modify gene regulatory networks in order to buffer metabolic perturbations. This is despite access to near universal technologies allowing genome-scale engineering. To help overcome this limitation, we have developed a pipeline enabling analysis of Transcription Regulation Integrated with MEtabolic Regulation (TRIMER). TRIMER utilizes a Bayesian network (BN) inferred from transcriptomic data to model the transcription factor regulatory network. TRIMER then infers the probabilities of gene states that are of relevance to the metabolism of interest, and predicts metabolic fluxes resulting from deletion of transcription factors at the genome scale. Additionally, we have developed a simulation framework to mimic the TF-regulated metabolic network, capable of generating both gene expression states and metabolic fluxes, thereby providing a fair evaluation platform for benchmarking models and predictions. Here, we present this computational pipeline. We demonstrate TRIMERs applicability to both simulated and experimental data and show that it outperforms current approaches on both data types.


title: Sensitive and specific spectral library searching with COSS and Percolator

abstract: Maintaining high sensitivity while limiting false positives is a key challenge in peptide identification from mass spectrometry data. Here, we therefore investigate the effects of integrating the machine learning-based post-processor Percolator into our spectral library searching tool COSS. To evaluate the effects of this post-processing, we have used forty data sets from two different projects and have searched these against the NIST and MassIVE spectral libraries. The searching is carried out using two spectral library search tools, COSS and MSPepSearch with and without Percolator post-processing, and using sequence database search engine MS-GF+ as a baseline comparator. The addition of the Percolator rescoring step to COSS is effective and results in a substantial improvement in sensitivity and specificity of the identifications. COSS is freely available as open source under the permissive Apache2 license, and binaries and source code are found at https://github.com/compomics/COSS


title: On the mechanical regulation of epithelial tissue homeostasis

abstract: Despite recent efforts to understand homeostasis in epithelial tissues, there are many unknowns surrounding this steady state. It is considered to be regulated by mechanoresponse, but unlike for single cells, this remains heavily debated for tissues. Here, we show that changes in matrix stiffness induce a non-equilibrium transition from tubular to squamous Madin-Darby Canine Kidney II tissues. Nonetheless, despite different cell morphologies and densities, all homeostatic tissues display equivalent topologies, which, hence, must be actively targeted and regulated. On the contrary, the mechanoresponse induces dramatic changes in the large-scale organization of the colonies. On stiff gels, this yields an unreported cooperative state of motile cells displaying higher densities than in the arrested homeostatic state. This suggests a more complex relation between cell density and motility than previously anticipated. Our results unequivocally relate the mechanosensitive properties of individual cells to the evolving macroscopic structures, an effect that could be important for understanding the emergent pathologies of living tissues.

category: biophysics