title: Meclizine and metabotropic glutamate receptor agonists attenuate severe pain and primary sensory neuron Ca2+ activity in chemotherapy-induced peripheral neuropathy

abstract: Chemotherapy-induced peripheral neuropathy (CIPN) affects about 68% of patients undergoing chemotherapy and causes severe neuropathic pain which is debilitating health problem and greatly reduces quality of life. Cisplatin is a commonly used platinum-based chemotherapeutic drug known to cause CIPN, possibly by causing oxidative stress damage to primary sensory neurons. Metabotropic glutamate receptors (mGluRs) are widely hypothesized to be involved in pain processing. Meclizine is an H1 histamine receptor antagonist which is known to have neuroprotective effects including anti-oxidative effect. Here, we used a mouse model of cisplatin-induced CIPN to test agonists of mGluR8 and group II mGluR as well as meclizine as interventions to reduce cisplatin-induced pain. We performed behavioral pain tests and in vivo entire DRG neurons Ca2+ imaging using genetically-encoded Ca2+ indicator, Pirt-GCaMP3 to monitor different drug interventions on a populational ensemble level. CIPN induced increased spontaneous Ca2+ activity in DRG neurons, increased Ca2+ transient amplitudes, and hyperresponses to mechanical, thermal, and chemical stimuli. We found mGluR8 agonist, DCPG, group II mGluR agonist, LY379268, and Histamine1 receptor antagonist, meclizine all significantly attenuated mechanical and thermal pain caused by CIPN. LY379268 and meclizine, but not DCPG, attenuated DRG neuronal Ca2+ activity elevated by CIPN. Furthermore, meclizine attenuated cisplatin-induced weight loss. These results suggest group II mGluR agonist, mGluR8 agonist, and meclizine are excellent candidates to study for mechanisms and new treatment option for CIPN.


title: In vivo direct imaging of neuronal activity at high temporo-spatial resolution

abstract: There has been a longstanding demand for noninvasive neuroimaging methods capable of detecting neuronal activity at both high temporal and spatial resolution. Here, we propose a novel method that enables Direct Imaging of Neuronal Activity for functional MRI (termed DIANA-fMRI) that can dynamically image spiking activity in milliseconds precision, while retaining the original benefit of high spatial resolution of MRI. DIANA-fMRI was demonstrated through in vivo mice brain imaging at 9.4 T applying electrical whisker-pad stimulation, directly imaging the spiking activity as well as capturing its sequential propagation along the thalamocortical pathway, as further confirmed through in vivo spike recording and optogenetics. DIANA-fMRI will open up new avenues in brain science by providing a deeper understanding of the brains functional organization including neural networks.

category: neuroscience


title: Intranasal vasopressin like oxytocin increases social attention by influencing top-down control, but additionally enhances bottom-up control

abstract: The respective roles of the neuropeptides arginine vasopressin (AVP) and oxytocin (OXT) in modulating social cognition and for therapeutic intervention in autism spectrum disorder have not been fully established. In particular, while numerous studies have demonstrated effects of oxytocin in promoting social attention the role of AVP has not been examined. The present study employed a randomized, double-blind, placebo (PLC)-controlled between-subject design to explore the social- and emotion-specific effects of AVP on both bottom-up and topdown attention processing with a validated emotional anti-saccade eye-tracking paradigm in 80 male subjects (PLC = 40, AVP = 40). Our findings showed that AVP increased the error rate for social (angry, fearful, happy, neutral and sad faces) but not non-social (oval shapes) stimuli during the anti-saccade condition and reduced error rates in the pro-saccade condition. Comparison of these findings with a previous study (sample size: PLC = 33, OXT = 33) using intranasal oxytocin revealed similar effects of the two peptides on anti-saccade errors but a significantly greater effect of AVP on pro-saccades. Both peptides also produced a post-task anxiolytic effect by reducing state anxiety. Together these findings suggested that both AVP and OXT decrease goal-directed top-down attention control to social salient stimuli but that AVP more potently increased bottom-up social attentional processing.


title: Met is required for oligodendrocyte progenitor cell migration in Danio rerio

abstract: During vertebrate central nervous system development, most oligodendrocyte progenitor cells (OPCs) are specified in the ventral spinal cord and must migrate throughout the neural tube until they become evenly distributed, occupying non-overlapping domains. While this process of developmental OPC migration is well characterized, the nature of the molecular mediators that govern it remain largely unknown. Here, using zebrafish as a model, we demonstrate that Met signaling is required for initial developmental migration of OPCs, and, using cell-specific knock-down of Met signaling, show that Met acts cell-autonomously in OPCs. Taken together, these findings demonstrate in vivo, the role of Met signaling in OPC migration and provide new insight into how OPC migration is regulated during development.


title: Coordinated regulation of CB1 cannabinoid receptors and anandamide metabolism stabilize network activity during homeostatic scaling down

abstract:
Neurons express overlapping homeostatic mechanisms to regulate synaptic function and network properties in response to perturbations of neuronal activity. Endocannabinoids (eCBs) are bioactive lipids synthesized in the post-synaptic compartments that regulate synaptic transmission, plasticity, and neuronal excitability throughout much of the brain, by activating pre-synaptic cannabinoid receptor CB1. The eCB system is well situated to regulate neuronal network properties and coordinate pre- and post-synaptic activity. However, the role of the eCB system in homeostatic adaptations to neuronal hyperactivity is unknown. We show that in mature cultured rat cortical neurons, chronic bicuculline treatment, known to induce homeostatic scaling-down, induces a coordinated adaptation to enhance tonic eCB signaling. Hyperexcitation triggers down regulation of fatty acid amide hydrolase (FAAH), the lipase that degrades the eCB anandamide. Subsequently, we measured an accumulation of anandamide and related metabolites, and an upregulation of total and cell surface CB1. We show that bicuculline induced downregulation of surface AMPA-type glutamate receptors and upregulation of CB1 occur through independent mechanisms. Finally, using live-cell microscopy of neurons expressing an extracellular fluorescent glutamate reporter (iGluSnFR), we confirm that cortical neurons in vitro exhibit highly synchronized network activity, reminiscent of cortical up-states in vivo. Up-state like activity in mature cortical cultures requires CB1 signaling under both control conditions and following chronic bicuculline treatment. We propose that during the adaptation to chronic neuronal hyperexcitation, tonic eCB signaling is enhanced through coordinated changes in anandamide metabolism and cell-surface CB1 expression to maintain synchronous network activity.

Significance statementNeurons are remarkably resilient to perturbations in network activities thanks to the expression of overlapping homeostatic adaptations. In response to network hyperactivity or silencing, neurons can respond through changes in excitatory and inhibitory post-synaptic neurotransmitter receptors, probability of pre-synaptic neurotransmitter release, or changes in membrane excitability. The endocannabinoid system is a prominent signaling pathway at many synapses that is known to be involved in multiple forms of short and long-term synaptic plasticity. Here we find that components of the endocannabinoid system are upregulated in response to chronic hyperexcitation of cultured cortical neurons, and that endocannabinoid signaling is required to maintain synchronized network activity. This work supports a novel tonic homeostatic function for the endocannabinoid system in neurons.


title: ChAlPred: A Web Server for Prediction of Allergenicity of Chemical Compounds

abstract: Allergy is the abrupt reaction of the immune system that may occur after the exposure with allergens like protein/peptide or chemical allergens. In past number of methods of have been developed for classifying the protein/peptide based allergen. To the best of our knowledge, there is no method to classify the allergenicity of chemical compound. Here, we have proposed a method named "ChAlPred", which can be used to fill the gap for predicting the chemical compound that might cause allergy. In this study, we have obtained the dataset of 403 allergen and 1074 non-allergen chemical compounds and used 2D, 3D and FP descriptors to train, test and validate our prediction models. The fingerprint analysis of the dataset indicates that PubChemFP129 and GraphFP1014 are more frequent in the allergenic chemical compounds, whereas KRFP890 is highly present in non-allergenic chemical compounds. Our XGB based model achieved the AUC of 0.89 on validation dataset using 2D descriptors. RF based model has outperformed other classifiers using 3D descriptors (AUC = 0.85), FP descriptors (AUC = 0.92), combined descriptors (AUC = 0.93), and hybrid model (AUC = 0.92) on validation dataset. In addition, we have also reported some FDA-approved drugs like Cefuroxime, Spironolactone, and Tioconazole which can cause the allergic symptoms. A user user-friendly web server named "ChAlPred" has been developed to predict the chemical allergens. It can be easily accessed at https://webs.iiitd.edu.in/raghava/chalpred/.

category: immunology


title: Maturation of Human Long-lived Plasma Cells Results in Resistance to Apoptosis by Transcriptional and Epigenetic Regulation

abstract: Antibody secreting cells (ASC) circulate after vaccination and migrate to the bone marrow (BM) where a subset known as long-lived plasma cells (LLPC) persist and secrete antibodies for a lifetime. The mechanisms of how circulating ASC become LLPC are not well elucidated. Here, we show that human blood ASCs have distinct morphology, transcriptomes, and epigenetics compared to BM LLPC. LLPC acquire transcriptional and epigenetic changes in the apoptosis pathway to support their survival. Upregulation of pro-survival gene expression accompanies downregulation of pro-apoptotic gene expression in LLPC. While pro-apoptotic gene loci are less accessible, pro-survival gene loci are not always accompanied by accessibility changes. Importantly, we show similar LLPC morphological and transcriptional maturation of blood ASC in response to the novel in vitro BM mimetic. In all, our study demonstrates that blood ASC in the BM microniche must undergo morphological and molecular changes to mature into apoptotic-resistant LLPC.


title: The endoplasmic reticulum associated degradation adaptor Sel1L regulates T cell homeostasis and function

abstract: Suppressor/Enhancer of Lin-12-like (Sel1L) is a critical adaptor for endoplasmic reticulum-associated degradation (ERAD), a process that maintains cellular protein quality control through degradation of misfolded proteins. Here we investigate the role of Sel1L in T cell homeostasis and function. T cell-specific deletion of Sel1L profoundly impairs peripheral T cell survival and promotes apoptotic cell death. Furthermore, Sel1L is required to maintain naive CD8+ T cell homeostasis in a cell-intrinsic manner with loss of quiescence as evidenced by increased proliferation. Sel1L-deficient T cells exhibit enhanced activation of the mammalian target of rapamycin (mTOR) pathway and altered cellular metabolism, including increased cellular reactive oxygen species, mitochondrial mass and mitochondrial membrane potential in the naive CD8+ T cell compartment. Furthermore, loss of Sel1L impaired CD8+ T cell immune responses following bacterial infection. These results demonstrate a novel role for Sel1L/ERAD in T cell homeostasis and function.


title: Ebola virus glycoprotein domains associated with protective efficacy

abstract: Ebola virus (EBOV) is the cause of sporadic outbreaks of human hemorrhagic disease in Africa, and the best-characterized virus in the filovirus family. The West Africa epidemic accelerated the clinical development of vaccines and therapeutics leading to licensure of vaccines and antibody-based therapeutics for human use in recent years. The most widely used vaccine is based on vesicular stomatitis virus (VSV) expressing the EBOV glycoprotein (GP)(VSV-EBOV). Due to its favorable immune cell targeting, this vaccine has also been used as base-vector for the development of second generation VSV-based vaccines against Influenza, Nipah, and Zika viruses. However, in these situations it may be beneficial if the immunogenicity against EBOV GP is minimized to induce a better protective immune response against the other foreign immunogen. Here, we analyzed if EBOV GP can be truncated to be less immunogenic yet still able to drive replication of the vaccine vector. We found that the EBOV GP glycan cap and the mucin-like domain are both dispensable for VSV-EBOV replication. The glycan cap domain, however, appears critical for mediating a protective immune response against lethal EBOV challenge in mice.

category: microbiology


title: Mitigating antimicrobial resistance with aspirin (acetylsalicylic acid) and paracetamol (acetaminophen): Conversion of doxycycline and minocycline resistant bacteria into sensitive in presence of aspirin and paracetamol

abstract: The emergence of antimicrobial resistance (AMR) stimulated research for alternatives antimicrobials, repurposing of other drugs, antibiotic adjuvants and alternative therapies for infections. Antimicrobial activity of NSAIDs is often reported and this study evaluated the antimicrobial potential of the two most common NSAIDs, aspirin (acetylsalicylic acid) and paracetamol (acetaminophen), against 293 clinical strains of bacteria. The ability of aspirin and paracetamol to convert minocycline and doxycyclin-resistant bacteria into sensitivity was also tested using micro-broth dilution assays used for determining minimum inhibitory concentration (MIC). Aspirin inhibited all 293 bacterial strains at [≤]10.24 mg/ mL concentration. Except for one strain each of Serratia grimaceae and S, aureus paracetamol inhibited none of the 293 strains at 10.24 mg/ mL. Of the 293 strains 116 (39.59%) were sensitive (MIC [≤]4 g/mL) to doxycycline and 127 (43.34%) to minocycline. Of the selected 57 minocycline-resistant (MIC >4 g/mL) strains aspirin converted 32 (56.14%) to minocycline-sensitive. Of the 49 doxycycline-resistant (MIC >4 g/mL) strains tested in presence of aspirin 30 (61.22%) turned sensitive. Of the 34 doxycycline-resistant strains tested in presence of paracetamol 11 (32.35%) become sensitive. The study concluded that most of the bacterial strains were not susceptible to aspirin and paracetamol at their concentrations often available in plasma at maximum therapeutic dose levels and had no significant change in their susceptibility to doxycycline and minocycline. The study indicated the potential of aspirin and its combination with antibiotics in the development of therapeutically useful topical antimicrobial formulations.


title: Silencing Doublesex expression triggers three-level pheromonal feminization in Nasonia males

abstract: The transcription factor Doublesex (Dsx) has a conserved function in controlling sexual morphological differences in insects, but our knowledge on its role in regulating sexual behavior is widely limited to Drosophila. Here, we show in the parasitoid wasp Nasonia vitripennis that males whose Dsx gene had been silenced by RNA interference (NvDsx-i) underwent a three-level pheromonal feminization: (1) NvDsx-i males were no longer able to attract females from a distance, owing to drastically reduced titers of the abdominal long-range sex pheromone. (2) NvDsx-i males were courted by wild-type males like females which correlated with a lower abundance of alkenes in their cuticular hydrocarbon (CHC) profiles. Supplementation of NvDsx-i male CHC profiles with realistic amounts of synthetic (Z)-9-hentriacontene (Z9C31), the most significantly reduced alkene in NvDsx-i males, interrupted courtship by wild-type conspecific males. Supplementation of female CHC profiles with Z9C31 reduced courtship and mating attempts by wild-type males. These results prove that Z9C31 is crucial for sex discrimination in Nasonia. (3) Nvdsx-i males were hampered in eliciting female receptivity during courtship and thus experienced severely reduced mating success, suggesting that they are unable to produce the hitherto unidentified oral aphrodisiac pheromone reported in N. vitripennis males. We conclude that Dsx is a multi-level key regulator of pheromone-mediated sexual communication in N. vitripennis. Silencing Dsx by RNA interference provides a new avenue for unraveling the molecular mechanisms underlying the pheromone-mediated sexual communication in insects.

category: molecular biology


title: Gains of 12p13.31 delay WNT-mediated initiation of hPSC differentiation and promote residual pluripotency in a cell cycle dependent manner

abstract: Though gains of chromosome 12p13.31 are highly recurrent in hPSC, their impact on differentiation is poorly understood. We identify a reduction in differentiation capacity towards all three germ layers and a subpopulation of residual pluripotent cells that appear during hepatic specification. These cells form as a result of the overexpression of NANOG and GDF3, whereby NANOG as the primary driver delays activation of WNT signaling, partly as a result of a direct physical interaction with TCF7. Entry into the residual state is determined by cell cycle position at the onset of differentiation and is maintained by a feedback loop between NANOG and GDF3. These findings highlight the ability of genetically abnormal hPSC to escape correct differentiation and to form residual pluripotent cells, an important risk in the safe clinical translation of hPSC. Our results further refine the molecular mechanisms that underpin the exit from pluripotency and onset of differentiation.

category: cell biology


title: Constraints and Opportunities for Detecting Land Surface Phenology in Drylands

abstract: Land surface phenology (LSP) enables global scale tracking of ecosystem processes, but its utility is limited in drylands due to low vegetation cover and resulting low annual amplitudes of vegetation indices (VIs). Due to the importance of drylands for biodiversity, food security, and the carbon cycle it is necessary to understand limitations in measuring dryland dynamics. Here, using simulated data and multi-temporal unmanned aerial vehicle (UAV) imagery of a desert shrubland, we explore the feasibility of detecting LSP with respect to fractional vegetation cover, plant functional types, VI uncertainty, and two different detection algorithms. Using simulated data we found that plants with distinct VI signals, such as deciduous shrubs, can require up to 60% fractional cover to consistently detect LSP. Evergreen plants, with lower seasonal VI amplitude, require considerably higher cover and can have undetectable phenology even with 100% vegetation cover. Our evaluation of two algorithms showed that neither performed the best in all cases. Even with adequate cover, biases in phenological metrics can still exceed 20 days, and can never be 100% accurate due to VI uncertainty from shadows, sensor view angle, and atmospheric interference. We showed how high-resolution UAV imagery enables LSP studies in drylands, and highlighted important scale effects driven by within canopy VI variation. With high-resolution imagery the open canopies of drylands are beneficial as they allow for straightforward identification of individual plants, enabling the tracking of phenology at the individual level. Drylands thus have the potential to become an exemplary environment for future LSP research.

category: ecology


title: Linking animal migration and ecosystem processes: data-driven simulation of propagule dispersal by migratory herbivores

abstract:
O_LIAnimal migration is a key process underlying active subsidies and species dispersal over long distances, which affects the connectivity and functioning of ecosystems. Despite much research describing patterns of where animals migrate, we still lack a framework for quantifying and predicting how animal migration affects ecosystem processes.
C_LIO_LIIn this study, we aim to integrate animal movement behavior and ecosystem functioning by developing a predictive modeling framework that can inform ecosystem management and conservation. Our framework models individual-level migration trajectories between populations seasonal ranges as well as the resulting dispersal and fate of propagules carried by the migratory animals, and it can be calibrated using empirical data at every step of the modeling process.
C_LIO_LIAs a case study, we applied our framework to model the spread of guava seeds, Psidium guajava, by a population of migratory Galapagos tortoises, Chelonoidis porteri, across Santa Cruz Island. Galapagos tortoises are large herbivores that transport seeds and nutrients across the island, while Guava is one of the most problematic invasive species in the Galapagos archipelago.
C_LIO_LIOur model is able to predict the pattern of spread of guava seeds alongside tortoises downslope migration range, and it identified areas most likely to see germination success and establishment. Our results show that Galapagos tortoises seed dispersal may particularly contribute to guava range expansion on Santa Cruz Island, due to both long gut retention time and tortoises long-distance migration across vegetation zones. In particular, we predict that tortoises are dispersing a significant amount of guava seeds into the Galapagos National Park, which has important consequences for the native flora.
C_LIO_LIThe flexibility and modularity of our framework allows for the integration of multiple data sources. It also allows for a wide range of applications to investigate how migratory animals affect ecosystem processes, including propagule dispersal but also other processes such as nutrient transport across ecosystems. Our framework is also a valuable tool for predicting how animal-mediated propagule dispersal can be affected by environmental change. These different applications can have important conservation implications for the management of ecosystems that include migratory animals.
C_LI


title: Intrapopulation adaptive variance supports selective breeding in a reef-building coral

abstract: The long-term persistence of coral reefs under climate change requires heritable drivers of thermal tolerance which support adaptation. The genomic basis of thermal tolerance has been evaluated across strong spatial and environmental gradients, but this variation also exists within populations due to neutral evolutionary processes. Small scale heterogeneity in coral bleaching is ubiquitous, so we used corals from a single reef to examine genomic signatures of bleaching performance, their biochemical correlates and the downstream consequences of selective breeding. In the absence of directional selection due to environmental differences, adult corals from a single population exhibit strong genomic patterns related to natural bleaching tolerance and symbiosis state, including functional differentiation in signaling pathways, protein and amino acid modification and metabolism. Conversely, growth, development and innate immune responses did not distinguish bleaching tolerance in adult corals. The genomic signatures of these gene ontologies influence biochemical patterns in healthy corals, primarily via cell-signaling pathway impacts on peptides and amino acids. Thermal tolerance in this population is highly heritable, with significantly higher survivorship under temperature stress in larvae and juveniles reared from thermally tolerant parents than those from sensitive parents. Using a select and re-sequence approach, certain gene ontologies were reproducibly impacted, while antioxidant activity and cell signaling ontologies were disproportionately selected in thermally tolerant corals, demonstrating the genomic drivers of successful selective breeding. These data show that intrapopulation variance in the absence of historical selection supports the adaptive capacity of coral reefs under climate change.


title: Use of an enclosed elk population to assess two non-invasive methods for estimating population size

abstract: Non-invasive genetic sampling and spatially explicit capture-recapture (SCR) models are used increasingly to estimate abundance of wildlife populations, but have not been adequately tested on gregarious animals such as elk (Cervus canadensis), for which correlated space use and movements violate model assumptions of independence. To evaluate the robustness and accuracy of SCR, and to assess the utility of an alternative non-invasive method for estimating density of gregarious ungulates, we utilized a tule elk (Cervus canadensis nannodes) population of known size within a fenced enclosure on the San Luis National Wildlife Refuge in central California. We evaluated fecal genetic SCR to camera trap-based random encounter model (REM) approaches to density estimation based on comparison to the true abundance. We also subsampled the dataset to explore the effects of varying search effort and elk density on the precision and accuracy of results. We found that SCR outperformed REM methods in the full datasets, and reliably provided accurate (relative bias <10%) and reasonably precise (relative standard error [≤]20%) estimates of density at moderately low to high densities (6-17 elk/km2), when the subsampling scenarios yielded a minimum average of 20 recaptures. We also found that the number of samples used to construct detection histories was a reliable predictor of precision, and could be used to establish minimum sampling requirements in future population surveys of elk. Although field-testing in free-ranging populations is needed, our results suggest that non-invasive genetic SCR is a promising tool for future population studies and monitoring of elk and potentially other gregarious ungulates. In contrast, the REM estimate of density was highly inaccurate, imprecise, and highly sensitive to camera parameters.


title: Diversification through sexual selection on gustatorial courtship traits in dwarf spiders

abstract: Sexual dimorphism can evolve under sexual selection or ecological factors. Sexually dimorphic male prosomal modifications are associated with gustatorial courtship in erigonines. The modifications vary from moderate elevations to bizarre shapes. Males transfer substances from these structures to females, which affect mate acceptance and fecundity. Here, we explore lability of these traits by investigating if modified prosomata are inherently linked to secretory glands, if glands evolved prior to prosomal modifications, and the possibility of convergent evolution and cryptic differentiation, aiming at assessing the possible role of this trait complex in speciation. We reconstructed the positions of glands and the musculature in the anterior part of prosomata of 76 erigonines and three outgroups using micro-CT. We incorporated these characters into an existing morphological character matrix and reanalyzed the phylogeny. Our results support the possession of glands as the ancestral state. The manifold modifications of the prosomal shape have evolved convergently. Differences in glandular positions between species with modified/unmodified prosomata suggest high lability of these traits. Cases of gland loss suggest considerable costs of gustatorial courtship. Our findings demonstrate divergent evolutionary patterns of these traits, and a likely facilitating effect of this type of sexual selection on speciation.

category: evolutionary biology


title: Synergistic epistasis of the deleterious effects of transposable elements

abstract: The replicative nature and generally deleterious effects of transposable elements (TEs) give rise to an outstanding question about how TE copy number is stably contained in host populations. Classic theoretical analyses predict that, when the decline in fitness due to each additional TE insertion is greater than linear, or when there is synergistic epistasis, selection against TEs can result in a stable equilibrium of TE copy number. While several mechanisms are predicted to yield synergistic deleterious effects of TEs, we lack empirical investigations of the presence of such epistatic interactions. Purifying selection with synergistic epistasis generates repulsion linkage between deleterious alleles and, accordingly, an underdispersed distribution for the number of deleterious mutations among individuals. We investigated this population genetic signal in an African Drosophila melanogaster population and found evidence for synergistic epistasis among TE insertions, especially those expected to have large fitness impacts. Curiously, even though ectopic recombination has long been predicted to generate nonlinear fitness decline with increased TE copy number, TEs predicted to suffer higher rates of ectopic recombination are not more likely to be underdispersed. On the other hand, underdispersed TE families are more likely to show signatures of deleterious epigenetic effects and stronger ping-pong signals of piRNA amplification, a hypothesized source from which synergism of TE-mediated epigenetic effects arises. Our findings set the stage for investigating the importance of epistatic interactions in the evolutionary dynamics of TEs.

category: genetics


title: Leveraging a founder population to identify novel rare-population genetic determinants of lipidome

abstract: Identifying the genetic determinants of inter-individual variation in lipid species (lipidome) may provide deeper understanding and new insight into the mechanistic effect of complex lipidomic pathways in CVD risk and progression beyond simple traditional lipids. Previous studies have been largely population based and thus only powered to discover associations with common genetic variants. Founder populations represent a powerful resource to accelerate discovery of novel biology associated with rare population alleles that have risen to higher frequency due to genetic drift. We performed a GWAS of 355 lipid species in 650 individuals from the Old Order Amish founder population including 127 lipid species not previously tested. We report for the first time the lipid species associated with two rare-population but Amish-enriched lipid variants: APOB_rs5742904 and APOC3_rs76353203. We also identified novel associations for 3 rare-population Amish-enriched loci with several sphingolipids and with proposed potential functional/causal variant in each locus including GLPTD2_rs536055318, CERS5_rs771033566, and AKNA_rs531892793. We replicated 7 previously known common loci including novel associations with two sterols: androstenediol with UGT locus on chromosome 2 and estriol with SLC22A8/A24 locus on chromosome 11. Our results show the power of founder populations to discover novel biology due to genetic drift that can increase the frequency of an allele from only few copies in large sample cohorts such as the UK Biobank to dozens of copies in sample size as small as 650.


title: Histone H3K27 methylation perturbs transcriptional robustness and underpins dispensability of highly conserved genes in fungi

abstract: Epigenetic modifications are key regulators of gene expression and underpin genome integrity. Yet, how epigenetic changes affect the evolution and transcriptional robustness of genes remains largely unknown. Here, we show how the repressive histone mark H3K27me3 influences the trajectory of highly conserved genes in fungi. We first performed transcriptomic profiling on closely related species of the plant pathogen Fusarium graminearum species complex. We determined transcriptional responsiveness of genes across environmental conditions to determine expression robustness. To infer evolutionary conservation of coding sequences, we used a comparative genomics framework of 23 species across the Fusarium genus. We integrated histone methylation data from three Fusarium species across the phylogenetic breadth of the genus. Gene expression variation is negatively correlated with gene conservation confirming that highly conserved genes show higher expression robustness. Furthermore, we show that highly conserved genes marked by H3K27me3 deviate from the typical housekeeping gene archetype. Compared to the genomic background, H3K27me3 marked genes encode smaller proteins, exhibit lower GC content, weaker codon usage bias, higher levels of hydrophobicity and are enriched for functions related to regulation and membrane transport. The evolutionary age of conserved genes with H3K27me3 histone marks falls typically within the origins of the Fusarium genus. We show that highly conserved genes marked by H3K27me3 are more likely to be dispensable for survival. Lastly, we show that conserved genes exposed to repressive H3K27me3 marks across distantly related fungi predict transcriptional perturbation at the microevolutionary scale in Fusarium fungi. In conclusion, we establish how repressive histone marks determine the evolutionary fate of highly conserved genes across evolutionary timescales.

category: genomics


title: Cell type-specific signal analysis in EWAS

abstract: Hundreds of epigenome-wide association studies (EWAS) have been performed, successfully identifying replicated epigenomic signals in processes such as ageing and smoking. Despite this progress, it remains a major challenge in EWAS to detect both cell type-specific and cell type confounding effects impacting study results. One way to identify these effects is through eFORGE (experimentally derived Functional element Overlap analysis of ReGions from EWAS), a published tool that uses 815 datasets from large-scale mapping studies to detect enriched tissues, cell types and genomic regions. Here, I show that eFORGE analysis can be extended to EWAS differentially variable positions (DVPs), identifying target cell types and tissues. In addition, I also show that eFORGE tissue-specific enrichment can be detected for sites below EWAS significance threshold. I develop on these and other analysis examples, extending our knowledge of eFORGE cell type- and tissue-specific enrichment results for different EWAS.


title: Interrogation of cancer gene dependencies reveals novel paralog interactions of autosome and sexchromosome encoded genes

abstract: Genetic networks are characterized by extensive buffering. During tumour evolution, disruption of these functional redundancies can create de novo vulnerabilities that are specific to cancer cells. In this regard, paralog genes are of particular interest, as the loss of one paralog gene can render tumour cells dependent on a remaining paralog. To systematically identify cancer-relevant paralog dependencies, we searched for candidate dependencies using CRISPR screens and publicly available loss-of-function datasets. Our analysis revealed >2,000 potential candidate dependencies, several of which were subsequently experimentally validated. We provide evidence that DNAJC15-DNAJC19, FAM50A-FAM50B and RPP25-RPP25L are novel cancer relevant paralog dependencies. Importantly, our analysis also revealed unexpected redundancies between sex chromosome genes. We show that chrX- and chrY- encoded paralogs, as exemplified by ZFX-ZFY, DDX3X-DDX3Y and EIF1AX-EIF1AY, are functionally linked so that tumour cell lines from male patients with Y-chromosome loss become exquisitely dependent on the chrX-encoded gene. We therefore propose genetic redundancies between chrX- and chrY- encoded paralogs as a general therapeutic strategy for human tumours that have lost the Y-chromosome.


title: Duplex-Repair enables highly accurate sequencing, despite DNA damage

abstract: Accurate DNA sequencing is crucial in biomedicine. Underlying the most accurate methods is the assumption that a mutation is true if altered bases are present on both strands of the DNA duplex. We now show that this assumption can be wrong. We establish that current methods to prepare DNA for sequencing, via  End Repair/dA-Tailing, may substantially resynthesize strands, leading amplifiable lesions or alterations on one strand to become indiscernible from true mutations on both strands. Indeed, we discovered that 7-17% and 32-57% of interior  duplex base pairs from cell-free DNA and formalin-fixed tumor biopsies, respectively, could be resynthesized in vitro and potentially introduce false mutations. To address this, we present Duplex-Repair, and show that it limits interior duplex base pair resynthesis by 8- to 464-fold, rescues the impact of induced DNA damage, and affords up to 8.9-fold more accurate duplex sequencing. Our study uncovers a major Achilles heel in sequencing and offers a solution to restore high accuracy.


title: ROS-Degradable Polythioketal Urethane Foam Dressings to Promote Porcine Skin Wound Repair

abstract: Impaired skin healing and progression into chronic wounds is a prevalent and growing medical problem. Porous, resorbable biomaterials can be used as temporary substrates placed into skin defects to support cell infiltration, neo-tissue formation, and remodeling of nonhealing wounds. Naturally-derived biomaterials have promising healing benefits, but their low mechanical properties and exuberant costs limit their performance and use. Synthetic materials can be affordably manufactured and tuned across a broader range of physiochemical properties, but opportunities remain for tailoring them for ideal host immune and regenerative responses. Polyesters are the most clinically-tested class of synthetic biomaterials, but their hydrolysis releases acidic degradation products that can cause autocatalytic degradation processes that are poorly controlled and are not tied to cellular or other biologic activities. Here, we systemically explored a series of ROS-degradable polythioketal (PTK) urethane (UR) foams with varied hydrophilicity as an alternative class of synthetic biomaterials for wound healing. It was found that the most hydrophilic PTK- UR variant, which had 7 ethylene glycol (EG7) repeats flanking each side of each thioketal bond, had the highest ROS reactivity of the PTK-URs tested. In an in vivo porcine excisional skin wound healing model, hydrophilic EG7 PTK-UR foams more effectively promoted tissue integration, ECM deposition, and re- epithelialization of full-thickness skin wound compared to more hydrophobic PTK-UR variants. Resolution of type 1 inflammation and lower foreign body response to scaffold remnants was also observed for EG7 versus more hydrophobic PTK-UR scaffolds. Finally, porcine wound healing studies showed that EG7 PTK-UR foams had similar wound healing response to a collagen-based clinical gold standard product, Integra Bilayer Wound Matrix (BWM), while outperforming polyester UR foam-based NovoSorb Biodegradable Temporizing Matrix (BTM) with respect to increased ECM production, vascularization, and biomaterial-associated immune phenotype. In sum, PTK-UR foams warrant further development toward a new class of synthetic biomaterial foams for skin wound healing applications.

category: bioengineering


title: Demuxalot: scaled up genetic demultiplexing for single-cell sequencing

abstract: Demultiplexing methods have facilitated the widespread use of single-cell RNA sequencing (scRNAseq) experiments by lowering costs and reducing technical variations. Here, we present demuxalot: a method for probabilistic genotype inference from aligned reads, with no assumptions about allele ratios and efficient incorporation of prior genotype information from historical experiments in a multi-batch setting. Our method efficiently incorporates additional information across reads originating from the same transcript, enabling up to 3x more calls per read relative to naive approaches. We also propose a novel and highly performant tradeoff between methods that rely on reference genotypes and methods that learn variants from the data, by selecting a small number of highly informative variants that maximize the marginal information with respect to reference single nucleotide variants (SNVs). Our resulting improved SNV-based demultiplex method is up to 3x faster, 3x more data efficient, and achieves significantly more accurate doublet discrimination than previously published methods. This approach renders scRNAseq feasible for the kind of large multi-batch, multi-donor studies that are required to prosecute diseases with heterogeneous genetic backgrounds.

category: bioinformatics


title: Distribution-free complex hypothesis testing for single-cell RNA-seq differential expression analysis

abstract: SO_SCPLOWUMMARYC_SCPLOWState-of-the-art methods for single-cell RNA sequencing (scRNA-seq) Differential Expression Analysis (DEA) often rely on strong distributional assumptions that are difficult to verify in practice. Furthermore, while the increasing complexity of clinical and biological single-cell studies calls for greater tool versatility, the majority of existing methods only tackle the comparison between two conditions. We propose a novel, distribution-free, and flexible approach to DEA for single-cell RNA-seq data. This new method, called ccdf, tests the association of each gene expression with one or many variables of interest (that can be either continuous or discrete), while potentially adjusting for additional covariates. To test such complex hypotheses, ccdf uses a conditional independence test relying on the conditional cumulative distribution function, estimated through multiple regressions. We provide the asymptotic distribution of the ccdf test statistic as well as a permutation test (when the number of observed cells is not sufficiently large). ccdf substantially expands the possibilities for scRNA-seq DEA studies: it obtains good statistical performance in various simulation scenarios considering complex experimental designs (i.e. beyond the two condition comparison), while retaining competitive performance with state-of-the-art methods in a two-condition benchmark. We apply ccdf to a large publicly available scRNA-seq dataset of 84,140 SARS-CoV-2 reactive CD8+ T cells, in order to identify the diffentially expressed genes across 3 groups of COVID-19 severity (mild, hospitalized, and ICU) while accounting for seven different cellular subpopulations.


title: Novel feature selection via kernel tensor decomposition for improved multi-omics data analysis

abstract:
BackgroundFeature selection of multi-omics data analysis remains challenging owing to the size of omics datasets, comprising approximately 102-105 features. In particular, appropriate methods to weight individual omics datasets are unclear, and the approach adopted has substantial consequences for feature selection. In this study, we extended a recently proposed kernel tensor decomposition (KTD)-based unsupervised feature extraction (FE) method to integrate multi-omics datasets obtained from common samples in a weight-free manner.

MethodKTD-based unsupervised FE was reformatted as the collection of kernelized tensors sharing common samples, which was applied to synthetic and real datasets.

ResultsThe proposed advanced KTD-based unsupervised FE method showed comparative performance to that of the previously proposed KTD method, as well as tensor decomposition-based unsupervised FE, but required reduced memory and central processing unit time. Moreover, this advanced KTD method, specifically designed for multi-omics analysis, attributes P-values to features, which is rare for existing multi-omics-oriented methods.

ConclusionsThe sample R code is available at https://github.com/tagtag/MultiR/


title: Novel feature selection method via kernel tensor decomposition for improved multi-omics data analysis


title: PARROT: a flexible recurrent neural network framework for analysis of large protein datasets

abstract: The rise of high-throughput experiments has transformed how scientists approach biological questions. The ubiquity of large-scale assays that can test thousands of samples in a day has necessitated the development of new computational approaches to interpret this data. Among these tools, machine learning approaches are increasingly being utilized due to their ability to infer complex non-linear patterns from high-dimensional data. Despite their effectiveness, machine learning (and in particular deep learning) approaches are not always accessible or easy to implement for those with limited computational expertise. Here we present PARROT, a general framework for training and applying deep learning-based predictors on large protein datasets. Using an internal recurrent neural network architecture, PARROT is capable of tackling both classification and regression tasks while only requiring raw protein sequences as input. We showcase the potential uses of PARROT on three diverse machine learning tasks: predicting phosphorylation sites, predicting transcriptional activation function of peptides generated by high-throughput reporter assays, and predicting the fibrillization propensity of amyloid beta with data generated by deep mutational scanning. Through these examples, we demonstrate that PARROT is easy to use, performs comparably to state-of-the-art computational tools, and is applicable for a wide array of biological problems.


title: Fast processing of environmental DNA metabarcoding sequence data using convolutional neural networks

abstract: 1The intensification of anthropogenic pressures have increased consequences on biodiversity and ultimately on the functioning of ecosystems. To monitor and better understand biodiversity responses to environmental changes using standardized and reproducible methods, novel high-throughput DNA sequencing is becoming a major tool. Indeed, organisms shed DNA traces in their environment and this "environmental DNA" (eDNA) can be collected and sequenced using eDNA metabarcoding. The processing of large volumes of eDNA metabarcoding data remains challenging, especially its transformation to relevant taxonomic lists that can be interpreted by experts. Speed and accuracy are two major bottlenecks in this critical step. Here, we investigate whether convolutional neural networks (CNN) can optimize the processing of short eDNA sequences. We tested whether the speed and accuracy of a CNN are comparable to that of the frequently used OBITools bioinformatic pipeline. We applied the methodology on a massive eDNA dataset collected in Tropical South America (French Guiana), where freshwater fishes were targeted using a small region (60pb) of the 12S ribosomal RNA mitochondrial gene. We found that the taxonomic assignments from the CNN were comparable to those of OBITools, with high correlation levels and a similar match to the regional fish fauna. The CNN allowed the processing of raw fastq files at a rate of approximately 1 million sequences per minute which was 150 times faster than with OBITools. Once trained, the application of CNN to new eDNA metabarcoding data can be automated, which promises fast and easy deployment on the cloud for future eDNA analyses.


title: Comparison of human and mouse tissues with focus on genes with no 1-to-1 homology

abstract:
We assembled a panel of 28 tissue pairs of human and mouse with RNA-Seq data on gene expression. We focused on genes with no 1-to-1 homology, because they pose special challenges. In this way, we identified expression patterns that identify and explain differences between the two species and suggest target genes for therapeutic applications. Here we mention three examples.

One pattern is observed by defining the aggregate expression of immunoglobulin genes (which have no homology) as a measure of different levels of an immune response. In Lung, we used this statistic to find genes that have significantly higher expression in low/moderate response, and thus they may be therapy targets: increasing their expression or mimicking their function with medications may help in recovery from inflammation in the lungs. Some of the observed associations are common to human and mouse; other associations involve genes involved in cell-to-cell signaling or in regeneration but were not known to be important in Lung.

Second pattern is that in the Small Intestine, mouse expresses much less antimicrobial defensins, while it has much higher expression of enzymes that are found to improve adaptive immune response. Such enzymes may be tested if they improve probiotic supplements that help in gut inflammation and other diseases.

Another pattern involves a many-to-many homology group of defensins that did not have a described function. In human tissues, expression of its genes was found only in a study of a disease of hair covered skin, but several of its genes are highly expressed in two tissues of our panel: mouse Skin and to a lesser degree mouse Vagina. This suggests that those genes or their homologs in other species may provide non-antibiotic medications for hair covered skin and other tissues with microbiome that includes fungi.


title: A proteomics sample metadata representation for multiomics integration, and big data analysis.

abstract:
The amount of public proteomics data is increasing at an extraordinary rate. Hundreds of datasets are submitted each month to ProteomeXchange repositories, representing many types of proteomics studies, focusing on different aspects such as quantitative experiments, post-translational modifications, protein-protein interactions, or subcellular localization, among many others. For every proteomics dataset, two levels of data are captured: the dataset description, and the data files (encoded in different file formats). Whereas the dataset description and data file formats are supported by all ProteomeXchange partner repositories, there is no standardized format to properly describe the sample metadata and their relationship with the dataset files in a way that fully allows their understanding or re-analysis. It is left to the users choice whether to provide or not an ad hoc document containing this information. Therefore, in many cases, understanding the study design and data requires going back to the associated publication. This can be tedious and may be restricted in the case of non-open access publications. In many cases, this problem limits the generalization and reuse of public proteomics data.

Here we present a standard representation for sample metadata tailored to proteomics datasets produced by the HUPO Proteomics Standards Initiative and supported by ProteomeXchange resources. We repurposed the existing data format MAGE-TAB used routinely in the transcriptomics field to represent and annotate proteomics datasets. MAGETAB-Proteomics defines a set of annotation rules that the datasets submitted to ProteomeXchange should follow, ranging from sample properties to data analysis protocols. We also introduce a crowdsourcing project that enabled the manual curation of over 200 public datasets using MAGE-TAB-Proteomics. In addition, we describe an ecosystem of tools and libraries that were developed to validate and submit sample metadata-related information to ProteomeXchange. We expect that these tools will improve the reproducibility of published results and facilitate the reanalysis and integration of public proteomics datasets.


title: Stromal HIF2 Regulates Immune Suppression in the Pancreatic Cancer Microenvironment

abstract:
Background & AimsPancreatic ductal adenocarcinoma (PDAC) has a hypoxic, immunosuppressive stroma, which contributes to its resistance to immune checkpoint blockade therapies. The hypoxia-inducible factors (HIFs) mediate the cellular response to hypoxia, but their role within the PDAC tumor microenvironment remains unknown.

MethodsWe used a dual recombinase mouse model to delete Hif1 or Hif2 in -smooth muscle actin (SMA)-expressing cancer-associated fibroblasts (CAFs) arising within spontaneous pancreatic tumors. The effects of CAF-Hif2 expression on tumor progression and composition of the tumor microenvironment were evaluated by Kaplan-Meier analysis, quantitative real-time polymerase chain reaction, histology, immunostaining, and by both bulk and single-cell RNA sequencing. CAF-macrophage crosstalk was modeled ex vivo using conditioned media from CAFs after treatment with hypoxia and PT2399, a HIF2 inhibitor currently in clinical trials. Syngeneic flank and orthotopic PDAC models were used to assess whether HIF2 inhibition improves response to immune checkpoint blockade.

ResultsCAF-specific deletion of HIF2, but not HIF1, suppressed PDAC tumor progression and growth, and improved survival of mice by 50% (n = 21-23 mice/group, Log-rank P = 0.0009). Deletion of CAF-HIF2 modestly reduced tumor fibrosis and significantly decreased the intratumoral recruitment of immunosuppressive M2 macrophages and regulatory T cells. Treatment with the clinical HIF2 inhibitor PT2399 significantly reduced in vitro macrophage chemotaxis and M2 polarization, and improved tumor responses to immunotherapy in both syngeneic PDAC mouse models.

ConclusionsTogether, these data suggest that stromal HIF2 is an essential component of PDAC pathobiology and is a druggable therapeutic target that could relieve tumor microenvironment immunosuppression and enhance immune responses in this disease.

category: cancer biology


title: Transcriptional and Functional Activity of Canine Hemangiosarcoma to Support Hematopoiesis Demonstrate Bone Marrow Nurse Cell Ontogeny

abstract: Hemangiosarcoma and angiosarcoma are soft tissue sarcomas of malignant blood vessel-forming cells in dogs and humans, respectively. These vasoformative sarcomas are aggressive and highly metastatic, with disorganized, irregular blood-filled vascular spaces. Our objective was to define molecular programs that support the niche, enabling progression of canine hemangiosarcoma and human angiosarcoma. Here, we show that the transcriptional landscape of canine hemangiosarcoma and human angiosarcoma included comparable angiogenic and inflammatory programs. Dog-in-mouse hemangiosarcoma xenografts recapitulated the vasoformative and highly angiogenic morphology and molecular characteristics of primary tumors. Blood vessels in the tumors were complex and disorganized, and they were lined by both donor and host cells, a trait that was not observed in xenografts from canine osteosarcoma and lymphoma. In some cases, the xenografted hemangiosarcoma cells created exuberant myeloid hyperplasia and gave rise to lymphoproliferative tumors of mouse origin. We did not uncover a definitive transmissible etiology, but our data indicate that transcriptional programs of hemangiosarcoma cells resemble those of hematopoietic nurse cells, and these malignant cells support expansion and differentiation of human hematopoietic progenitors. We conclude that canine hemangiosarcoma, and possibly human angiosarcoma, originate from nurse cells that make up the stromal bone marrow niche, and that these cells may also support the growth of hematopoietic tumors.


title: Transcriptional and Functional Activity of Hemangiosarcoma Support Bone Marrow Nurse Cell Ontogeny


title: Rapamycin induced hyperglycemia is associated with exacerbated age-related osteoarthritis

abstract:
BackgroundThe objective of this study was to determine if mechanistic target of rapamycin (mTOR) inhibition with or without AMP-activated protein kinase (AMPK) activation can protect against primary, age-related OA.

DesignDunkin-Hartley guinea pigs develop mild primary OA pathology by 5-months of age that progresses to moderate OA by 8-months of age. At 5-months, guinea pigs sacrificed as young control (n=3) or were fed either a control diet (n=8), a diet enriched with the mTOR-inhibitor rapamycin (Rap, 14ppm, n=8), or Rap with the AMPK-activator metformin (Rap+Met, 1000ppm, n=8) for 12 weeks. Knee joints were evaluated by OARSI scoring, micro-computed tomography, and immunohistochemistry. Glenohumeral articular cartilage was collected for western blotting.

ResultsRap and Rap+Met treated guinea pigs displayed lower body weight than control. Rap and Rap+Met inhibited articular cartilage mTORC1 but not mTORC2 signaling. Rap+Met, but not Rap alone, stimulated AMPK. Despite lower body weight and articular cartilage mTORC1 inhibition, Rap and Rap+Met treated guinea pigs had greater OA severity in the medial tibial plateau due to articular cartilage structural damage and/or proteoglycan loss. Rap and Rap+Met increased plasma glucose compared to control. Plasma glucose concentration was positively correlated with proteoglycan loss, suggesting hyperglycemic stress may have contributed to worsened OA.

ConclusionsThis is the first study to show that Rap induced increase in plasma glucose was associated with greater OA severity. Further, articular cartilage mTORC1 inhibition and bodyweight reduction by dietary Rap and Rap+Met did not protect against primary OA during the prevailing hyperglycemia.

category: physiology


title: Machine learning-guided acyl-ACP reductase engineering for improved in vivo fatty alcohol production

abstract: Fatty acyl reductases (FARs) catalyze the reduction of thioesters to alcohols and are key enzymes for the microbial production of fatty alcohols. Many existing metabolic engineering strategies utilize these reductases to produce fatty alcohols from intracellular acyl-CoA pools; however, acting on acyl-ACPs from fatty acid biosynthesis has a lower energetic cost and could enable more efficient production of fatty alcohols. Here we engineer FARs to preferentially act on acyl-ACP substrates and produce fatty alcohols directly from the fatty acid biosynthesis pathway. We implemented a machine learning-driven approach to iteratively search the protein fitness landscape for enzymes that produce high titers of fatty alcohols in vivo. After ten design-test-learn rounds, our approach converged on engineered enzymes that produce over twofold more fatty alcohols than the starting natural sequences. We further characterized the top identified sequence and found its improved alcohol production was a result of an enhanced catalytic rate on acyl-ACP substrates, rather than enzyme expression or KM effects. Finally, we analyzed the sequence-function data generated during the enzyme engineering to identify sequence and structure features that influence fatty alcohol production. We found an enzymes net charge near the substrate-binding site was strongly correlated with in vivo activity on acyl-ACP substrates. These findings suggest future rational design strategies to engineer highly active enzymes for fatty alcohol production.

category: synthetic biology


title: Schmidtea sp., from the S-W Romania (Platyhelminthes, Tricladida, Dugesiidae)

abstract: The morphology and the anatomy of the copulatory apparatus in a Dugesiidae population from the SW Romania are presented. The copulatory apparatus is characterized by intermingled bursal canal musculature and two distinct penis bulbs with two large seminal vesicles. Based on these morphological characters, the population is assigned to the "lugubris-polychroa" group of species, now belonging to the genus Schmidtea (de Vries & Sluys 1991). The copulatory apparatus is also characterised by the presence of an atrial fold, characteristic of S. mediterranea. The assign of the morphotype here presented to the species level is delayed until integrative molecular analysis.

category: zoology


title: Automated Pipeline for Infants Continuous EEG (APICE): a flexible pipeline for developmental studies

abstract: Infant electroencephalography (EEG) presents several challenges compared with adult data. Recordings are typically short. Motion artifacts heavily contaminate the data. The EEG neural signal and the artifacts change throughout development. Traditional data preprocessing pipelines have been developed mainly for event-related potentials analyses, and they required manual steps, or use fixed thresholds for rejecting epochs. However, larger datasets make the use of manual steps infeasible, and new analytical approaches may have different preprocessing requirements. Here we propose an Automated Pipeline for Infants Continuous EEG (APICE). APICE is fully automated, flexible, and modular. Artifacts are detected using multiple algorithms and adaptive thresholds, making it suitable to different age groups and testing procedures without redefining parameters. Artifacts detection and correction of transient artifacts is performed on continuous data, allowing for better data recovery and providing flexibility (i.e., the same preprocessing is usable for different analyses). Here we describe APICE and validate it using two infant datasets of different ages tested in different experimental paradigms. We also tested the combination of APICE with common data cleaning methods such as Independent Component Analysis and Denoising Source Separation. APICE uses EEGLAB and compatible custom functions. It is freely available at https://github.com/neurokidslab/eeg_preprocessing, together with example scripts.


title: An improved method for evaluating inverted encoding models

abstract: Here we present a more interpretable and flexible approach for reconstructing the contents of perception, attention, and memory from neuroimaging data. Our enhanced inverted encoding model (eIEM) incorporates methodological improvements including proper accounting of population-level tuning functions and a trial-by-trial prediction error-based metric where reconstruction quality is measured in meaningful units. Improved flexibility is further gained via eIEMs novel goodness-of-fit feature: for trial-by-trial reconstructions, goodness-of-fits are obtained independently (non-circularly) to prediction error and can be applied to any IEM procedure or decoding metric, resulting in improved reconstruction quality and brain-behavior correlations. We validate eIEM from methodological principles, simulated neuroimaging datasets, and three pre-existing fMRI datasets spanning perception, attention, and working memory. Notably, eIEM is easy to apply and broadly accessible - our Python package (https://pypi.org/project/inverted-encoding) implements eIEM in one line of code - and is easily modifiable to compare performance metrics and/or scale up to more complex models.


title: An enhanced inverted encoding model for neural reconstructions


title: Tools for accurate post hoc determination of marker location within whole-brain microscopy images

abstract: To interpret in vivo experiments designed to understand brain function, high-resolution whole-brain microscopy provides a means for post hoc determination of the location of implanted devices and recorded cells in three dimensional brain space that is a critical step for data interrogation. Here we have developed Python-based tools (brainreg and brainreg-segment) to accurately map, in a common coordinate space, the position of dye-labelled probe tracks and two-photon imaged cell populations expressing fluorescent protein. The precise location of probes and cells were validated using physiological recordings and human raters that indicate accuracy levels to less than 70{micro}m. These flexible, open-source methodologies are expected to further evolve with need and to deliver the anatomical precision that is necessary for understanding the functional architecture of the brain.


title: Afferent connections of cytoarchitectural area 6M and surrounding cortex in the marmoset: putative homologues of the supplementary and pre-supplementary motor areas

abstract: Cortical projections to the caudomedial frontal cortex were studied using retrograde tracers in marmosets. We tested the hypothesis that cytoarchitectural area 6M includes homologues of the supplementary and pre-supplementary motor areas (SMA and preSMA) of other primates. We found that, irrespective of the injection sites location within 6M, over half of the labeled neurons were located in motor and premotor areas. Other connections originated in prefrontal area 8b, ventral anterior and posterior cingulate areas, somatosensory areas (3a and 1-2), and areas on the rostral aspect of the dorsal posterior parietal cortex. Although the origin of afferents was similar, injections in rostral 6M received higher percentages of prefrontal afferents, and fewer somatosensory afferents, compared to caudal injections, compatible with differentiation into SMA and preSMA. Injections rostral to 6M (area 8b) revealed a very different set of connections, with increased emphasis in prefrontal and posterior cingulate afferents, and fewer parietal afferents. The connections of 6M were also quantitatively different from those of M1, dorsal premotor areas, and cingulate motor area 24d. These results show that the cortical motor control circuit is conserved in simian primates, indicating that marmosets can be valuable models for studying movement planning and control.


title: Mast cell-specific receptor/corticotropin-releasing factor axis regulates alcohol withdrawal-associated headache

abstract: Rehabilitation from alcohol addiction or abuse is challenging due to alcohol withdrawal symptoms. Headache is a severe alcohol withdrawal symptom that frequently contributes to rehabilitation failure. Despite the need for treating alcohol withdrawal-induced headache, there is no appropriate therapeutic option available. Development of improved therapeutics will depend on obtaining a clearer understanding of alcohol withdrawal-induced headache pain mechanisms. Here, we show that the mast cell-specific receptor MrgprB2 controls development of alcohol withdrawal-induced headache. Withdrawing alcohol from alcohol-acclimated mice induces strong headache behaviors, including facial allodynia, facial pain expressions, and reduced walking movement, symptoms often observed in humans suffering from headache. Observed pain behaviors were abolished in MrgprB2-deficient mice. We observed in vivo spontaneous activation and hypersensitization of trigeminal ganglia neurons in alcohol withdrawal mice, but not in MrgprB2-deficient mice. Corticotropin-releasing factor (CRF) was increased in dura mater after alcohol withdrawal. Injection of CRF into dura mater resulted in activation of trigeminal ganglia neurons and vasodilation, which was accompanied by headache behavior. In cells, CRF evoked Ca2+ transients via MrgprB2 or human MrgprX2. The results indicate that alcohol withdrawal causes headache via mast cell degranulation in dura mater. The process is under control of MrgprB2/MrgprX2, which would appear to represent a potential target for treating alcohol withdrawal-related headache.


title: Molecular architecture determines brain delivery of a transferrin-receptor targeted lysosomal enzyme

abstract:
Delivery of biotherapeutics across the blood-brain barrier (BBB) is a challenge. Many approaches fuse biotherapeutics to platforms that bind the transferrin receptor (TfR), a brain endothelial cell target, to facilitate receptor-mediated transcytosis across the BBB. Here, we characterized the pharmacological behavior of two distinct TfR-targeted platforms fused to iduronate 2-sulfatase (IDS), a lysosomal enzyme deficient in mucopolysaccharidosis type II (MPS II), and compared the relative brain exposures and functional activities of both approaches in mouse models. IDS fused to a moderate-affinity, monovalent TfR binding enzyme transport vehicle (ETV:IDS) resulted in widespread brain exposure, internalization by parenchymal cells, and significant substrate reduction in the CNS of an MPS II mouse model. In contrast, IDS fused to a standard high-affinity bivalent antibody (IgG:IDS) resulted in lower brain uptake, limited biodistribution beyond brain endothelial cells, and reduced brain substrate reduction. These results highlight important features likely to impact the clinical development of TfR-targeting platforms in MPS II and potentially other CNS diseases.

SummaryBrain delivery, biodistribution and pharmacodynamics of a lysosomal enzyme fused to a moderate-affinity transferrin receptor-directed blood-brain barrier enzyme transport vehicle are superior to a traditional high-affinity anti-TfR monoclonal antibody fusion.


title: Predicting behavior through dynamic modes in resting-state fMRI data

abstract: Dynamic properties of resting-state functional connectivity (FC) provide rich information on brainbehavior relationships. Dynamic mode decomposition (DMD) has been used as a method to characterize FC dynamics. However, it remains unclear whether dynamic modes (DMs), spatial-temporal coherent patterns computed by DMD, provide information about individual behavioral differences. This study established a methodological approach to predict individual differences in behavior using DMs. Furthermore, we investigated the contribution of DMs within each of seven specific frequency bands (0-0.1,...,0.6-0.7 Hz) for prediction. To validate our approach, we confirmed whether each of 59 behavioral measures could be predicted by performing multivariate pattern analysis on a gram matrix, which was created using subject-specific DMs computed from resting-state functional magnetic resonance imaging (rs-fMRI) data of individuals. The prediction was successful, and DMD outperformed temporal independent component analysis, a conventional data decomposition method for extracting spatial activity patterns. Most of the behavioral measures that showed significant prediction accuracies in a permutation test were cognitive-behavioral measures. Our results suggested that DMs within frequency bands <0.2 Hz primarily contributed to prediction. In addition, we found that DMs <0.2 Hz had spatial structures similar to several common resting-state networks. We demonstrated the effectiveness of DMs, indicating that DMD is a key approach for extracting spatiotemporal features from rs-fMRI data.


title: The Neural Correlates of Working Memory Training in Typically Developing Children

abstract: Working memory training improves childrens cognitive performance on untrained tasks; however, little is known about the underlying neural mechanisms. We investigated this in 32 typically developing children aged 10-14 years (19 girls and 13 boys; Devon, UK) using a randomized controlled design and multi-modal MRI. Training improved working memory performance and increased intrinsic functional connectivity between the bilateral intraparietal sulci. Furthermore, improvements in working memory were associated with greater recruitment of the left middle frontal gyrus on a complex span task. The repeated engagement of fronto-parietal regions during training may increase their activity and functional connectivity over time, affording greater working memory performance. We discuss the plausibility of generalizable cognitive benefits from a neurobiological perspective and implications for neurodevelopmental theory. This is not the version of record.


title: The Neural Correlates of Working Memory Training in Typically Developing Children - Working Paper


title: Repertoire of naturally acquired maternal antibodies transferred to infants for protection against shigellosis

abstract: Shigella is the second leading cause of diarrheal diseases, accounting for >200,000 infections and >50,000 deaths in children under 5 years of age annually worldwide. The incidence of Shigella-induced diarrhea is relatively low during the first year of life and increases substantially, reaching its peak between 11 to 24 months of age. This epidemiological trend hints at an early protective immunity of maternal origin and an increase in disease incidence when maternally acquired immunity wanes. The magnitude, type, antigenic diversity, and antimicrobial activity of maternal antibodies transferred via placenta that can prevent shigellosis during early infancy are not known. To address this knowledge gap, Shigella-specific antibodies directed against the lipopolysaccharide (LPS) and virulence factors (IpaB, IpaC, IpaD, IpaH, and VirG), and antibody-mediated serum bactericidal (SBA) and opsonophagocytic killing antibody (OPKA) activity were measured in maternal and cord blood sera from a longitudinal cohort of mother-infant pairs living in rural Malawi. Protein-specific (very high levels) and Shigella LPS IgG were detected in maternal and cord blood sera; efficiency of placental transfer was 100% and 60%, respectively, and had preferential IgG subclass distribution (protein-specific IgG1 > LPS-specific IgG2). In contrast, SBA and OPKA activity in cord blood was substantially lower as compared to maternal serum and varied among Shigella serotypes. LPS was identified as the primary target of SBA and OPKA activity. Maternal sera had remarkably elevated Shigella flexneri 2a LPS IgM, indicative of recent exposure. Our study revealed a broad repertoire of maternally acquired antibodies in infants living in a Shigella-endemic region and highlights the abundance of protein-specific antibodies and their likely contribution to disease prevention during the first months of life. These results contribute new knowledge on maternal infant immunity and target antigens that can inform the development of vaccines or therapeutics that can extend protection after maternally transferred immunity wanes.


title: Istradefylline, an adenosine A2a receptor antagonist, ameliorates neutrophilic airway inflammation and psoriasis in mice

abstract:
ObjectiveExtracellular adenosine is produced from secreted ATP by cluster of differentiation (CD)39 and CD73. Both are critical nucleotide metabolizing enzymes of the adenosine generating pathway and are secreted by neuronal or immune cells. Adenosine plays a role in energy processes, neurotransmission, and endogenous regulation of inflammatory responses. Istradefylline is a selective adenosine A2a receptor (A2aR) antagonist used for the treatment of Parkinsons disease. We have reported that adenosine primes hypersecretion of interleukin (IL)-17A via A2aR. Istradefylline, as well as an inhibitor of CD39 (ARL67156) and an inhibitor of CD73 (AMP-CP), suppressed IL-17A production, and the administration of istradefylline to mice with experimental autoimmune encephalomyelitis (EAE) led to the marked amelioration of the disease. These previous results suggest that adenosine is an endogenous modulator of neutrophilic inflammation. We investigated the effect of istradefylline, ARL67156 and AMP-CP on other mouse models of neutrophilic inflammation.

MethodsWe tested the effect of istradefylline, ARL67156 and AMP-CP on OVA-induced neutrophilic airway inflammation or imiquimod (IMQ)-induced psoriasis in mice. These two model mice received these drugs orally or percutaneously, respectively. The production of IL-17A in the lung and ear thickness were used as an index of the effects.

ResultsWe show that istradefylline, ARL67156 and AMP-CP suppressed the OVA-induced IL-17A production in the lung and IMQ-induced psoriasis.

ConclusionThese results indicate that adenosine-mediated IL-17A production plays a role in neutrophilic inflammation models, and moreover, istradefylline, ARL67156, and AMP-CP are effective in animal models of neutrophilic inflammation. Some clinical relevancies in COVID-19 are discussed.


title: Virus-Dependent Immune Conditioning of Tissue Microenvironments

abstract: A thorough understanding of complex spatial host-disease interactions in situ is necessary in order to develop effective preventative measures and therapeutic strategies. Here, we developed Protein And Nucleic acid IN situ Imaging (PANINI) and coupled it with Multiplexed Ion Beam Imaging (MIBI) to sensitively and simultaneously quantify DNA, RNA, and protein levels within the microenvironments of tissue compartments. The PANINI-MIBI approach was used to measure over 30 parameters simultaneously across large sections of archival lymphoid tissues from non-human primates that were healthy or infected with simian immunodeficiency virus (SIV), a model that accurately recapitulates human immunodeficiency virus infection (HIV). This enabled multiplexed dissection of cellular phenotypes, functional markers, viral DNA integration events, and viral RNA transcripts as resulting from viral infection. The results demonstrated immune coordination from an unexpected upregulation of IL10 in B cells in response to SIV infection that correlated with macrophage M2 polarization, thus conditioning a potential immunosuppressive environment that allows for viral production. This multiplexed imaging strategy also allowed characterization of the coordinated microenvironment around latently or actively infected cells to provide mechanistic insights into the process of viral latency. The spatial multi-modal framework presented here is applicable to deciphering tissue responses in other infectious diseases and tumor biology.


title: Co-option of local and systemic immune responses by the hormonal signalling system triggering metamorphosis

abstract: Insect metamorphosis is regulated by the production, secretion and degradation of two peripheral hormones: 20-hydroxyecdysone (ecdysone) and juvenile hormone (JH). In addition to their roles in developmental regulation, increasing evidence suggests that these hormones are involved in innate immunity processes, such as phagocytosis and the induction of antimicrobial peptide (AMP) production. AMP regulation includes systemic responses as well as local responses at surface epithelia that contact with the external environment. At pupariation, Drosophila melanogaster increases dramatically the expression of three AMP genes, drosomycin (drs), drosomycin-like 2 (drsl2) and drosomycin-like 5 (drsl5). We show that the systemic action of drs at pupariation is dependent on ecdysone signalling in the fat body and operates via the ecdysone downstream target, Broad-Z4. In parallel, ecdysone also regulates local responses, specifically through the activation of drsl2 expression in the gut. Finally, we confirm the relevance of this ecdysone dependent AMP expression for the control of bacterial persistence by showing that flies lacking drs expression in the fat body have higher bacterial persistence over metamorphosis. Together, our data emphasize the importance of the association between ecdysone signalling and immunity using in vivo studies and establish a new role for ecdysone at pupariation, which impacts developmental success by regulating the immune system in a stage-dependent manner. We speculate that this co-option of immune effectors by the hormonal system may constitute a pre-emptive mechanism to control bacterial numbers in the pupa at the core of metamorphosis evolution.

category: developmental biology


title: Co-option of immune effectors by the hormonal signalling system triggering metamorphosis in Drosophila melanogaster


title: Noninvasive genetic surveys before and after a megafire detect displacement of migratory mule deer

abstract: Due to climate change and past logging and fire suppression, the western US are experiencing increasingly large and frequent wildfires. Understanding how wildlife respond to these mega-fires is becoming increasingly relevant to protect and manage these populations. However, the lack of predictability inherent in such events makes studies difficult to plan. We took advantage of a large high-severity wildfire that burned adjacent to an ongoing study of mule deer (Odocoileus hemionus) on their summer range upslope of the fire to investigate their displacement onto our study area both immediately and upon their return to summer range the following year. We used spatial capture-recapture models in conjunction with noninvasive fecal DNA sampling to estimate density and non-spatial Pradel robust-design models to estimate apparent survival and recruitment rates. Compared to density before the fire, we observed an increase in deer density and an increase in per-capita recruitment rates one month after the fire. These findings suggest that the immediate response of at least some deer was to flee the fire upslope onto the study area rather than to downslope toward their winter range. These changes did not carry over into the following year, however, suggesting that deer formerly using the burned area as summer range may have returned there despite the high severity of the fire, or may have chosen new areas for their summer range. This suggests that, at least in the short term, the fire did not negatively affect the deer population.


title: Genetic isolation by distance underlies color pattern divergence in red-eyed treefrogs (Agalychnis callidryas)

abstract: Investigating the spatial distribution of genetic and phenotypic variation can provide insights into the evolutionary processes that shape diversity in natural systems. We characterized patterns of genetic and phenotypic diversity to learn about drivers of color-pattern diversification in red-eyed treefrogs (Agalychnis callidryas) in Costa Rica. Along the Pacific coast, red-eyed treefrogs have conspicuous leg color patterning that transitions from orange in the north to purple in the south. We measured phenotypic variation of frogs, with increased sampling at sites where the orange-to-purple transition occurs. At the transition zone, we discovered the co-occurrence of multiple color-pattern morphs. To explore possible causes of this variation, we generated a SNP dataset to analyze population genetic structure, measure genetic diversity, and infer the processes that mediate genotype-phenotype dynamics. We investigated how patterns of genetic relatedness correspond with individual measures of color pattern along the coast, including testing for the role of hybridization in geographic regions where orange and purple phenotypic groups co-occur. We found no evidence that color-pattern polymorphism in the transition zone arose through recent hybridization. Instead, a strong pattern of genetic isolation by distance (IBD) indicates that color-pattern variation was either retained through other processes such as ancestral color polymorphisms or ancient secondary contact, or else it was generated by novel mutations. We found that phenotype changes along the Pacific coast more than would be expected based on genetic divergence and geographic distance alone. Combined, our results suggest the possibility of selective pressures acting on color pattern at a small geographic scale.


title: A novel SERPINE1-FOSB fusion gene in a case of pseudomyogenic hemangioendothelioma results in activation of intact FOSB and the PI3K-AKT-mTOR signaling pathway and responsiveness to sirolimus

abstract: Pseudomyogenic hemangioendothelioma (PHE) is an extremely rare disease that affects mainly the young and more men than women. PHE are multicentric, locally aggressive, have low metastatic potential, and affect multiple tissue planes. Genetic aberrations are frequently detected in PHE and may play important roles in the occurrence, development, and treatment of this disease. In this study, we report a case of PHE with a novel SERPINE1-FOSB fusion gene. The fusion introduced a strong promoter near the coding region of FOSB, resulting in overexpression of intact FOSB. Immunohistochemical analysis showed overexpression of pAKT and mTOR in tumor cells, suggesting activation of the PI3K-AKT-mTOR signaling pathway. The patient responded well to targeted therapy with sirolimus, an mTOR inhibitor. Our study correlated dysregulation of a specific signaling pathway and the effectiveness of a targeted therapy to a specific genetic aberration. This information may be useful for future investigations of targeted therapeutics and provide a potential predictive biomarker for therapeutic effectiveness in cases of PHE with this genetic aberration.


title: Redundant enhancers in the iab-5 domain cooperatively activate Abd-B in the A5 and A6 abdominal segments of Drosophila

abstract:
The homeotic Abdominal-B (Abd-B) gene belongs to Bithorax complex and is regulated by four regulatory domains named iab-5, iab-6, iab-7 and iab-8, each of which is thought to be responsible for directing the expression of Abd-B in one of the abdominal segments from A5 to A8. It is assumed that male specific features of the adult cuticle in A5 is solely dependent on regulatory elements located in iab-5, while the regulatory elements in the iab-6 are both necessary and sufficient for the proper differentiation of the A6 cuticle. Unexpectedly, we found that this long held assumption is not correct. Instead, redundant tissue-specific enhancers located in the iab-5 domain are required for the proper activation of Abd-B not only in A5 but also in A6. Our study of deletions shows that the iab-5 initiator is essential for the functioning of the iab-5 enhancers in A5, as well as for the correct differentiation of A6. This requirement is circumvented by deletions that remove the initiator and most of the iab-5 regulatory domain sequences. While the remaining iab-5 enhancers are inactive in A5, they are activated in A6 and contribute to the differentiation of this segment. In this case, Abd-B stimulation by the iab-5 enhancers in A6 depends on the initiators in the iab-4 and iab-6 domains.

Summary StatementIn Drosophila, the segmental-specific expression of the homeotic gene Abdominal-B in the abdominal segments is regulated by autonomous regulatory domains. We demonstrated cooperation between these domains in activation of Abdominal-B.


title: recount3: summaries and queries for large-scale RNA-seq expression and splicing

abstract: We present recount3, a resource consisting of over 750,000 publicly available human and mouse RNA sequencing (RNA-seq) samples uniformly processed by our new Monorail analysis pipeline. To facilitate access to the data, we provide the recount3 and snapcount R/Bioconductor packages as well as complementary web resources. Using these tools, data can be downloaded as study-level summaries or queried for specific exon-exon junctions, genes, samples, or other features. Monorail can be used to process local and/or private data, allowing results to be directly compared to any study in recount3. Taken together, our tools help biologists maximize the utility of publicly available RNA-seq data, especially to improve their understanding of newly collected data. recount3 is available from http://rna.recount.bio.


title: Evolutionary puzzle: discussing the evolution of sentience in Metazoa through a phylogenetic perspective

abstract: Sentience is the capacity of organisms to feel and experience through subjective states. During the last years, several investigations have indicated that response mechanisms to harmful stimuli can be highly conserved among the Metazoa. However, there is a research bias towards vertebrates in the available studies. Here we discuss the evolution of the nervous and sensory system, pain and nociception in animals through a phylogenetic perspective testing the hypothesis of common ancestry of sentience. Our results indicate that characteristics related to sentience - morphological and molecular and behavioural -, were already present in the common ancestors of Metazoa, Eumetazoa and Bilateria. Our phylogenetic hypotheses positioned Porifera as the sister-group to all the other Metazoa, corroborating the hypothesis of a single origin of the nervous system. Our results also depict Urbilateria as the ancestor of the metazoan toolkit related to the sentience. These scenarios suggest that some attributes of the sensory system may have appeared even before the emergence of the nervous system, through possible cooptations of sensory modules of the first Metazoa.

category: animal behavior and cognition


title: No one tool to rule them all: Prokaryotic gene prediction tool performance is highly dependent on the organism of study

abstract:
MotivationThe biases in Open Reading Frame (ORF) prediction tools, which have been based on historic genomic annotations from model organisms, impact our understanding of novel genomes and metagenomes. This hinders the discovery of new genomic information as it results in predictions being biased towards existing knowledge. To date users have lacked a systematic and replicable approach to identify the strengths and weaknesses of any ORF prediction tool and allow them to choose the right tool for their analysis.

ResultsWe present an evaluation framework (ORForise) based on a comprehensive set of 12 primary and 60 secondary metrics that facilitate the assessment of the performance of ORF prediction tools. This makes it possible to identify which performs better for specific use-cases. We use this to assess 15 ab initio and model-based tools representing those most widely used (historically and currently) to generate the knowledge in genomic databases. We find that the performance of any tool is dependent on the genome being analysed, and no individual tool ranked as the most accurate across all genomes or metrics analysed. Even the top-ranked tools produced conflicting gene collections which could not be resolved by aggregation. The ORForise evaluation framework provides users with a replicable, data-led approach to make informed tool choices for novel genome annotations and for refining historical annotations.

Availabilityhttps://github.com/NickJD/ORForise

Contactnicholas@dimonaco.co.uk

Supplementary informationSupplementary data are available at bioRxiv online.


title: DecoPath: A web application for decoding pathway enrichment analysis

abstract: The past two decades have brought a steady growth of pathway databases and pathway enrichment methods. However, the advent of pathway data has not been accompanied by an improvement with regards to interoperability across databases, thus, hampering the use of pathway knowledge from multiple databases for pathway enrichment analyses. While integrative databases have attempted to address this issue by collating pathway knowledge from multiple resources, these approaches do not account for redundant information across them. On the other hand, the majority of studies that employ pathway enrichment analyses still rely upon a single database, though the use of another resource could yield differing results, which is similarly the case when different pathway enrichment methods are employed. These shortcomings call for approaches that investigate the differences and agreements across databases and enrichment methods as their selection in the experimental design of a pathway analysis can be a crucial first step in ensuring the results of such an analysis are meaningful. Here we present DecoPath, a web application to assist in the interpretation of the results of pathway enrichment analysis. DecoPath provides an ecosystem to run pathway enrichment analysis or directly upload results and facilitate the interpretation of these results with custom visualizations that highlight the consensus and/or discrepancies at the pathway- and gene-levels. DecoPath is available at https://decopath.scai.fraunhofer.de and its source code and documentation can be found on GitHub at https://github.com/DecoPath/DecoPath.


title: Comparison of Resources and Methods to infer Cell-Cell Communication from Single-cell RNA Data

abstract:
The growing availability of single-cell data has sparked an increased interest in the inference of cell-cell communication from this data. Many tools have been developed for this purpose. Each of them consists of a resource of intercellular interactions prior knowledge and a method to predict potential cell-cell communication events. Yet the impact of the choice of resource and method on the resulting predictions is largely unknown. To shed light on this, we created a framework, available at https://github.com/saezlab/ligrec_decoupler, to facilitate a comparative assessment of methods for inferring cell-cell communication from single cell transcriptomics data and then compared 15 resources and 6 methods. We found few unique interactions and a varying degree of overlap among the resources, and observed uneven coverage in terms of pathways and biological categories. We analysed a colorectal cancer single cell RNA-Seq dataset using all possible combinations of methods and resources. We found major differences among the highest ranked intercellular interactions inferred by each method even when using the same resources. The varying predictions lead to fundamentally different biological interpretations, highlighting the need to benchmark resources and methods.

FindingsO_LIBuilt a framework to systematically combine 15 resources and 6 methods to estimate cell-cell communication from single-cell RNA data
C_LIO_LICell-cell communication resources are often built from the same original databases and very few interactions are unique to a single resource. Yet overlap varies among resources and certain biological terms are unevenly represented
C_LIO_LIDifferent methods and resources provided notably different results
C_LIO_LIThe observed disagreement among the methods could have a considerable impact on the interpretation of results
C_LI


title: Vcflib and tools for processing the VCF variant call format

abstract:
Since its introduction in 2011 the variant call format (VCF) has been widely adopted for processing DNA and RNA variants in practically all population studies -- as well as in somatic and germline mutation studies. VCF can present single nucleotide variants, multi-nucleotide variants, insertions and deletions, and simple structural variants called against a reference genome. Here we present over 125 useful and much used free and open source software tools and libraries, part of vcflib tools and bio-vcf. We also highlight cyvcf2, hts-nim and slivar tools. Application is typically in the comparison, filtering, normalisation, smoothing, annotation, statistics, visualisation and exporting of variants. Our tools run daily and invisibly in pipelines and countless shell scripts. Our tools are part of a wider bioinformatics ecosystem and we consider it very important to make these tools available as free and open source software to all bioinformaticians so they can be deployed through software distributions, such as Debian, GNU Guix and Bioconda. vcflib, for example, was installed over 40,000 times and bio-vcf was installed over 15,000 times through Bioconda by December 2020. We shortly discuss the design of VCF, lessons learnt, and how we can address more complex variation that can not easily be represented by the VCF format. All source code is published under free and open source software licenses and can be downloaded and installed from https://github.com/vcflib.

Author summaryMost bioinformatics workflows deal with DNA/RNA variations that are typically represented in the variant call format (VCF) -- a file format that describes mutations (SNP and MNP), insertions and deletions (INDEL) against a reference genome. Here we present a wide range of free and open source software tools that are used in biomedical sequencing workflows around the world today.


title: Non-coding NFKBIZ 3' UTR mutations promote cell growth and resistance to targeted therapies in diffuse large B-cell lymphoma

abstract:
Amplifications and non-coding 3' UTR mutations affecting NFKBIZ have been identified as recurrent genetic events in diffuse large B-cell lymphoma (DLBCL). We confirm the prevalence and pattern of NFKBIZ 3' UTR mutations in independent cohorts and determine they are enriched in the ABC subtype as well as the recently described novel BN2/C1/NOTCH2 classes of DLBCL. Presently, the effects of and mechanism by which non-coding mutations can act as cancer drivers has been relatively unexplored. Here, we provide a functional characterization of these non-coding NFKBIZ 3' UTR mutations. We demonstrate that the resulting elevated expression of I{kappa}B-{zeta} confers growth advantage in DLBCL cell lines and primary germinal center B-cells as well as nominate novel I{kappa}B-{zeta} target genes with potential therapeutic implications. The limited responses to targeted treatments in DLBCL, particularly those targeting the NF-{kappa}B axis, led us to investigate and confirm that NFKBIZ 3' UTR mutations affect response to therapeutics and suggest it may be a useful predictive biomarker.

Statement of SignificanceThrough functional characterization we reveal that non-coding NFKBIZ 3' UTR mutations are a common driver in DLBCL, and mutation status may be a relevant biomarker to predict poor response to therapeutics targeting the NF-{kappa}B pathway.


title: Partial reprogramming restores youthful gene expression through transient suppression of cell identity

abstract: Transient induction of pluripotent reprogramming factors has been reported to reverse some features of aging in mammalian cells and tissues. However, the impact of transient reprogramming on somatic cell identity programs and the necessity of individual pluripotency factors remain unknown. Here, we mapped trajectories of transient reprogramming in young and aged cells from multiple murine cell types using single cell transcriptomics to address these questions. We found that transient reprogramming restored youthful gene expression in adipocytes and mesenchymal stem cells but also temporarily suppressed somatic cell identity programs. We further screened Yamanaka Factor subsets and found that many combinations had an impact on aging gene expression and suppressed somatic identity, but that these effects were not tightly entangled. We also found that a transient reprogramming approach inspired by amphibian regeneration restored youthful gene expression in aged myogenic cells. Our results suggest that transient pluripotent reprogramming poses a neoplastic risk, but that restoration of youthful gene expression can be achieved with alternative strategies.

category: systems biology