title: Viral susceptibility across host species is largely independent of dietary protein to carbohydrate ratios

abstract:
The likelihood of a successful host shift of a parasite to a novel host species can be influenced by environmental factors that can act on both the host and parasite. Changes in nutritional resource availability have been shown to alter pathogen susceptibility and the outcome of infection in a range of systems. Here we examined how dietary protein to carbohydrate altered susceptibility in a large cross infection experiment. We infected 27 species of Drosophilidae with an RNA virus on three food types of differing protein to carbohydrate ratios. We then measured how viral load and mortality across species was affected by changes in diet. We found that changes in the protein:carbohydrate in the diet did not alter the outcomes of infection, with strong positive inter-species correlations in both viral load and mortality across diets, suggesting no species by diet interaction. Mortality and viral load were strongly positively correlated, and this association was consistent across diets. This suggests changes in diet may give consistent outcomes across host species, and may not be universally important in determining host susceptibility to pathogens.

Twitter summaryNo role of host diet in susceptibility to a novel viral pathogen across host species

Impact StatementA successful host shift of a parasite from one susceptible species to a novel host can be influenced by many ecological factors. Changes in host diet can alter the immune response and outcomes of host-parasite interactions and could potentially alter the outcome of a virus host shift. To investigate, we infected 27 species of Drosophilidae with an RNA virus (DCV) across three food types with differing protein to carbohydrate ratios. We then looked at pathogen loads and survival of infected hosts compared to uninfected controls. Changes in the ratio of protein to carbohydrate did not alter susceptibility to DCV across host species.


title: The wheat 55K SNP-based exploration of loci for spikelet number per spike from a tetraploid wheat (Triticum turgidum L.) recombinant inbred line population derived from a Chinese landrace Ailanmai and a wild emmer accession

abstract: The spikelet number per spike (SNS) is a primary factor determining wheat yield. Common wheat breeding reduces the genetic diversity among elite germplasm resources, leading to a detrimental effect on future wheat production. It is, therefore, necessary to explore new genetic resources for SNS to increase wheat yield. A tetraploid landrace  Ailanmai x wild emmer wheat recombinant inbred line (RIL) population was used to construct a genetic map using the wheat 55K single nucleotide polymorphism (SNP) array. The linkage map containing 1150 bin markers with a total genetic distance of 2411.8 cM was obtained. Based on phenotypic data from eight environments and best linear unbiased prediction (BLUP) values, five quantitative trait loci (QTL) for SNS were identified, explaining 6.71-29.40% of the phenotypic variation. Two of them, QSns.sau-AM-2B.3 and QSns.sau-AM-3B.2, were detected as the major and novel QTL. Their effects were further validated in two additional F2 populations using tightly linked Kompetitive Allele-Specific PCR (KASP) markers. Potential candidate genes within the physical intervals of corresponding QTL were predicted to participate in inflorescence development and spikelet formation. Genetic associations between SNS and other agronomic traits were also detected and analyzed. This study demonstrates the feasibility of wheat 55K SNP array developed for common wheat in genetic mapping of tetraploid population and shows the potential application of wheat related species in wheat improvement programs.

category: genetics


title: The 55K SNP-Based Exploration of QTL for Spikelet Number Per Spike in a Tetraploid Wheat (Triticum turgidum L.) Population: Chinese Landrace Ailanmai X Wild Emmer


title: Inhibitory interneurons show early dysfunction in a SOD1 mouse model of ALS

abstract:
Few studies in amyotrophic lateral sclerosis (ALS) focus on the inhibitory interneurons synapsing onto motoneurons (MNs). However, inhibitory interneurons could contribute to dysfunction, particularly if altered before MN neuropathology and establish a long-term imbalance of inhibition / excitation. We directly assessed excitability and morphology of glycinergic (GlyT2) interneurons located throughout the ventral horn in the lumbar enlargement from SOD1G93AGlyT2eGFP (SOD1) and wildtype GlyT2eGFP (WT) mice on postnatal day 6 to 10. Patch clamp revealed dampened excitability in SOD1 interneurons, including depolarized persistent inward currents (PICs), depolarized threshold for firing action potentials, and a shortened afterhyperpolarization (AHP). SOD1 inhibitory interneurons also had smaller somata but primary dendrites showed larger volume and surface area than WT. GlyT2 interneurons were then divided into 3 subgroups based on location: (1) interneurons within 100 m of the ventral white matter, where Renshaw cells (RCs) are located, (2) interneurons interspersed with MNs in lamina IX, and (3) interneurons in the intermediate ventral area including laminae VII and VIII. Ventral interneurons were the most profoundly affected, including more depolarized PICs, smaller somata and larger primary dendrites. Interneurons in lamina IX had depolarized PIC onset, smaller somata and longer primary dendrites. In lamina VII-VIII, interneurons were largely unaffected, mainly showing smaller somata. In summary, inhibitory interneurons show very early region-specific perturbations poised to impact excitatory / inhibitory balance of MNs, modify motor output, and provide early biomarkers of ALS. Therapeutics like riluzole which universally reduce CNS excitability could exacerbate the inhibitory dysfunction described here.

Abstract Figure

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=144 SRC="FIGDIR/small/348359v3_ufig1.gif" ALT="Figure 1">
View larger version (24K):
org.highwire.dtl.DTLVardef@17e639forg.highwire.dtl.DTLVardef@be91d4org.highwire.dtl.DTLVardef@1582b05org.highwire.dtl.DTLVardef@17071a5_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOAbstract FigureC_FLOATNO SOD1 glycinergic interneurons in the ventral horn show altered morphology and excitability, including depolarization of PICs, depolarized threshold, shorter AHPs, smaller somata and larger primary dendrites. Ventrally located interneurons are the most prominently affected.

C_FIG Key Points SummaryO_LISpinal inhibitory interneurons could contribute to amyotrophic lateral sclerosis (ALS) pathology, but their excitability has never been directly measured.
C_LIO_LIWe studied the excitability and morphology of glycinergic interneurons in postnatal transgenic mice (SOD1G93AGlyT2eGFP).
C_LIO_LIInterneurons were less excitable and had smaller somas but larger primary dendrites in SOD1 mice.
C_LIO_LIGlyT2 interneurons were analyzed according to their localization within the ventral spinal cord. Interestingly, the greatest differences were observed in the most ventrally-located interneurons.
C_LIO_LIWe conclude that inhibitory interneurons show presymptomatic changes that may contribute to excitatory / inhibitory imbalance in ALS.
C_LI

category: neuroscience


title: Inhibitory interneurons show early dysfunction in a SOD1 mouse model of Amyotrophic Lateral Sclerosis


title: Explainable Machine Learning Approach to Predict and Explain the Relationship between Task-based fMRI and Individual Differences in Cognition

abstract: Despite decades of costly research, we still cannot accurately predict individual differences in cognition from task-based fMRI. Moreover, aiming for methods with higher prediction is not sufficient. To understand brain-cognition relationships, we need to explain how these methods draw brain information to make the prediction. Here we applied an explainable machine-learning (ML) framework to predict cognition from task-based fMRI during the n-back working-memory task, using data from the Adolescent Brain Cognitive Development (n=3,989). We compared nine predictive algorithms in their ability to predict 12 cognitive abilities. We found better out-of-sample prediction from ML algorithms over the mass-univariate and OLS multiple regression. Among ML algorithms, Elastic Net, a linear and additive algorithm, performed either similar to or better than non-linear and interactive algorithms. We explained how these algorithms drew information, using SHapley Additive explanation, eNetXplorer, Accumulated Local Effects and Friedmans H-statistic. These explainers demonstrated benefits of ML over the OLS multiple regression. For example, ML provided some consistency in variable importance with a previous study (Sripada et al. 2020) and consistency with the mass-univariate approach in the directionality of brain-cognition relationships at different regions. Accordingly, our explainable-ML framework predicted cognition from task-based fMRI with boosted prediction and explainability over standard methodologies.


title: Compensatory hepatic adaptation accompanies permanent absence of intrahepatic biliary network due to YAP1 loss in liver progenitors

abstract: YAP1 regulates cell plasticity during liver injury, regeneration and cancer, but its role in liver development is unknown. YAP1 activity was detected in biliary cells and in cells at the hepato-biliary bifurcation in single-cell RNA-sequencing analysis of developing livers. Hepatoblast deletion of Yap1 led to no impairment in Notch-driven SOX9+ ductal plate formation, but prevented the formation of the abutting second layer of SOX9+ ductal cells, blocking the formation of a patent intrahepatic biliary tree. Intriguingly, the mice survived for 8 months with severe cholestatic injury and without any hepatocyte-to-biliary transdifferentiation. Ductular reaction in the perihilar region suggested extrahepatic biliary proliferation likely seeking the missing intrahepatic biliary network. Long-term survival of these mice occurred through hepatocyte adaptation via reduced metabolic and synthetic function including altered bile acid metabolism and transport. Overall, we show YAP1 as a key regulator of bile duct development while highlighting a profound adaptive capability of hepatocytes.

category: pathology


title: Context dependent activation and repression of enhancers by Hunchback binding sites in Drosophila embryo

abstract:
Hunchback (Hb) is considered a context-dependent transcription factor, able to activate or repress different enhancers during Drosophila embryo segmentation. The mechanism driving the contextdependent activity of Hb is however not well understood. Here we measure the activity of a large set of 20 synthetic enhancers that we design to elucidate the effect of Hb binding sites in Drosophila segmentation. We obtain quantitative data on the spatiotemporal dynamics of activity of all synthetic enhancers in-vivo, by using a quantitative and sensitive reporter system we recently developed. Our data reveal the dual role of Hb binding sites in segmentation enhancers: on the one hand, Hb act as a typical short range repressor by binding to its cognate sequences; on the other hand, we report a novel effect of a sequence containing multiple Hb binding sites, which is able to increase enhancer activity independently from Hb binding. This sequence, which contains multiple Poly-dA stretches, increases the activity of enhancers driven by different activators, possibly by disfavoring nucleosome occupancy.

AUTHOR SUMMARYThe control of gene expression is a fundamental process that allows cells to respond to external stimuli and take on various identities in complex organisms. Enhancers are DNA sequences that play a key role in this process. In the simplest model of an enhancer, small parts of its sequence can be specifically bound by proteins called transcription factors and the occupancy pattern of these proteins on the enhancer determines the expression level of a specific gene. In this research work we have studied enhancers in the context of the development of a fruit fly embryo. We have built synthetic enhancer sequences containing binding sites for a few specific factors and measured their activity in living embryos using fluorescence microscopy. Our results revealed that binding sites for a particular protein, Hunchback, are able to influence the activity of the enhancer even independently from Hunchback binding to them. This discovery might help to explain the complex effects that have been observed when studying Hunchback binding sites in natural enhancers.

category: developmental biology


title: PharmaNet: Pharmaceutical discovery with deep recurrent neural networks.

abstract: The discovery and development of novel pharmaceuticals is an area of active research mainly due to the large investments required and long payback times. As of 2016, the development of a novel drug candidate required up to $ USD 2.6 billion in investment for only 10% rate of approval by the FDA. To help decreasing the costs associated with the process, a number of in silico approaches have been developed with relatively low success due to limited predicting performance. Here, we introduced a machine learning-based algorithm as an alternative for a more accurate search of new pharmacological candidates, which takes advantage of Recurrent Neural Networks (RNN) for active molecule prediction within large databases. Our approach, termed PharmaNet was implemented here to search for ligands against specific cell receptors within 102 targets of the DUD-E database, which contains 22886 active molecules. PharmaNet comprises three main phases. First, a SMILES representation of the molecule is converted into a raw molecular image. Second, a convolutional encoder processes the data to obtain a fingerprint molecular image that is finally analyzed by a Recurrent Neural Network (RNN). This approach enables precise predictions of the molecules target on the basis of the feature extraction, the sequence analysis and the relevant information filtered out throughout the process. Molecule Target prediction is a highly unbalanced detection problem and therefore, we propose that an adequate evaluation metric of performance is the area under the Normalized Average Precision (NAP) curve. PharmaNet largely surpasses the previous state-of-the-art method with 97.7% in the Receiver Operating Characteristic curve (ROC-AUC) and 65.5% in the NAP curve. We obtained a perfect performance for human farnesyl pyrophosphate synthase (FPPS), which is a potential target for antimicrobial and anticancer treatments. We decided to test PharmaNet for activity prediction against FPPS by searching in the CHEMBL data set. We obtained three (3) potential inhibitors that were further validated through both molecular docking and in silico toxicity prediction. Most importantly, one of this candidates, CHEMBL2007613, was predicted as a potential antiviral due to its involvement on the PCDH17 pathway, which has been reported to be related to viral infections.

category: bioinformatics


title: A deep learning algorithm for 3D cell detection in whole mouse brain image datasets

abstract: Understanding the function of the nervous system necessitates mapping the spatial distributions of its constituent cells defined by function, anatomy or gene expression. Recently, developments in tissue preparation and microscopy allow cellular populations to be imaged throughout the entire rodent brain. How-ever, mapping these neurons manually is prone to bias and is often impractically time consuming. Here we present an opensource algorithm for fully automated 3D detection of neuronal somata in mouse whole-brain microscopy images using standard desktop computer hardware. We demonstrate the applicability and power of our approach by mapping the brain-wide locations of large populations of cells labeled with cytoplasmic fluorescent proteins expressed via retrograde trans-synaptic viral infection.


title: Genetic effects on brain traits impact cell-type specific gene regulation during neurogenesis

abstract: Interpretation of the function of non-coding risk loci for neuropsychiatric disorders and brain-relevant traits via gene expression and alternative splicing is generally performed in bulk post-mortem adult tissue. However, genetic risk loci are enriched in regulatory elements active during neocortical differentiation, and regulatory effects of risk variants may be masked by heterogeneity in bulk tissue. Here, we map expression quantitative trait loci (eQTLs), splicing QTLs (sQTLs), and allele specific expression in primary human neural progenitors (n=85) and their sorted neuronal progeny (n=74), identifying numerous loci not detected in either bulk developing cortical wall or adult cortex. Using colocalization and genetic imputation via transcriptome wide association, we uncover cell-type specific regulatory mechanisms underlying risk for brain-relevant traits that are active during neocortical differentiation. Specifically, we identified a progenitor-specific eQTL for CENPW co-localized with common variant associations of cortical area and educational attainment.


title: Executive function supports single-shot endowment of value to arbitrary transient goals

abstract: People often learn from the outcomes of their actions, even when these outcomes do not involve material rewards or punishments. How does our brain provide this flexibility? We combined behavior, computational modeling, and functional neuroimaging to probe whether learning from transient goals harnesses the same circuitry that supports learning from secondary reinforcers. Behavior and neuroimaging revealed that "one-shot" transient goals (abstract fractal images seen once) can act as a substitute for rewards during instrumental learning, and produce reliable reward-like signals in dopaminergic reward circuits. Moreover, we found evidence that prefrontal correlates of executive control may play a role in shaping these responses in reward circuits. These results suggest that learning from abstract goal outcomes is supported by an interplay between high-level representations in prefrontal cortex and low-level responses in subcortical reward circuits. This interaction may allow humans to perform reinforcement learning over flexible, arbitrarily abstract reward functions.


title: Forty years of change in a monomictic tropical lake

abstract: Many lakes across the world are entering novel states and experiencing altered biogeochemical cycling due to local anthropogenic stressors. In the tropics, understanding the drivers of these changes can be difficult due to a lack of documented historic conditions or an absence of continuous monitoring that can distinguish between intra- and inter-annual variation. Over the last forty years (1980-2020), Lake Yojoa (Honduras) has experienced increased watershed development as well as the introduction of a large net-pen Tilapia farm, resulting in a dramatic reduction in seasonal water clarity, increased trophic state and altered nutrient dynamics, shifting Lake Yojoa from an oligotrophic (low productivity) to mesotrophic (moderate productivity) ecosystem. To assess the changes that have occurred in Lake Yojoa as well as putative drivers for those changes, we compared Secchi depth (water clarity), dissolved inorganic nitrogen (DIN), and total phosphorus (TP) concentrations at continuous semi-monthly intervals for the three years between 1979-1983 and again at continuous 16-day intervals for 2018-2020. Between those two periods we observed the loss of a clear water phase that previously occurred in the months when the water column was fully mixed. Seasonal peaks in DIN coincident with mixing suggest that an enhanced accumulation of ammonium in the hypolimnion (the bottom layer of a stratified lake) during stratification, and release to the epilimnion (the top layer of a stratified lake) with mixing maintains high algal abundance and subsequently low secchi depth during what was previously the clear water phase. This interaction of nutrient loading and Lake Yojoas monomictic stratification regime illustrates a key phenomenon in how physical water column structure and nutrients interact in tropical monomictic lakes. This work highlights the need to consider nutrient dynamics of warm anoxic hypolimnions, not just surface water nutrient concentrations, to understand environmental change in these societally important but understudied ecosystems.

category: ecology


title: The interaction of physical structure and nutrient loading drives ecosystem change in a large tropical lake over 40 years


title: Plant Scaffolds Support Motor Recovery and Regeneration in Rats after Traumatic Spinal Cord Injury

abstract: As of yet, no standard of care incorporates biomaterials to treat traumatic spinal cord injury (SCI). However, intense development of biomaterials for treating SCI has focused on fabricating microscale channels to support the regrowth of axons while minimizing scar formation. We previously demonstrated that plant tissues could be decellularized and processed to form sterile, biocompatible and implantable biomaterials that support cell infiltration and vascularization in vivo. Vascularized plant tissues contain continuous microscale channels with geometries relevant for supporting neural regeneration. We hypothesized that decellularized vascular bundles would support neural regeneration and motor recovery in SCI. Sprague Dawley rats received a complete T8-T9 spinal cord transection and were implanted with acellular plant-derived scaffolds and allowed to recover over 28 weeks. Animals that received the scaffolds alone, with no other therapeutic compounds, demonstrated a significant and stable partial improvement in motor function compared to control animals as early as week 4 post-injury. Hind-limb motor function did not deteriorate over the remaining 28 weeks. Histological analysis revealed minimal astrocyte scarring at the spinal cord - scaffold interface, aligned axonal projection through the scaffolds, populations of serotonergic neurons and Schwann cells, laminin and collagen deposition and the presence of blood vessels. Axonal reconnection via the scaffold was also confirmed by Fluro-gold retrograde tracing. Taken together, our work defines a novel route for building upon naturally occurring plant microarchitectures to support the repair of the spinal cord post-injury. Notably, these results were achieved without the use of growth factors, stem/progenitor cells, or any other interventions.

category: bioengineering


title: Distinct Phenotypes of SARS-CoV-2 Isolates Reveal Viral Traits Critical for Replication in Primary Human Respiratory Cells

abstract: Since entering the human population, SARS-CoV-2 (the causative agent of COVID-19) has spread worldwide, causing >100 million infections and >2 million deaths. While large-scale sequencing efforts have identified numerous genetic variants in SARS-CoV-2 during its circulation, it remains largely unclear whether many of these changes impact adaptation, replication or transmission of the virus. Here, we characterized 14 different low-passage replication-competent human SARS-CoV-2 isolates representing all major European clades observed during the first pandemic wave in early 2020. By integrating viral sequencing data from patient material, virus stocks, and passaging experiments, together with kinetic virus replication data from non-human Vero-CCL81 cells and primary differentiated human bronchial epithelial cells (BEpCs), we observed several SARS-CoV-2 features that associate with distinct phenotypes. Notably, naturally-occurring variants in Orf3a (Q57H) and nsp2 (T85I) were associated with poor replication in Vero-CCL81 cells but not in BEpCs, while SARS-CoV-2 isolates expressing the Spike D614G variant generally exhibited enhanced replication abilities in BEpCs. Strikingly, low-passage Vero-derived stock preparation of 3 SARS-CoV-2 isolates selected for substitutions at positions 5/6 of E, and were highly attenuated in BEpCs, revealing a key cell-specific function to this region. Rare isolate-specific deletions were also observed in the Spike furin-cleavage site during Vero-CCL81 passage, but these were rapidly selected against in BEpCs, underscoring the importance of this site for SARS-CoV-2 replication in primary human cells. Overall, our study uncovers sequence features in SARS-CoV-2 variants that determine cell-specific virus replication, and highlights the need to monitor SARS-CoV-2 stocks carefully when phenotyping newly emerging variants or potential variants-of-concern.

category: microbiology


title: SARS-CoV-2 Variants Reveal Features Critical for Replication in Primary Human Cells


title: Automated Feature Extraction from Large Cardiac Electrophysiological Data Sets

abstract:
RationaleA new multi-electrode array-based application for the long-term recording of action potentials from electrogenic cells makes possible exciting cardiac electrophysiology studies in health and disease. With hundreds of simultaneous electrode recordings being acquired over a period of days, the main challenge becomes achieving reliable signal identification and quantification.

ObjectiveWe set out to develop an algorithm capable of automatically extracting regions of high-quality action potentials from terabyte size experimental results and to map the trains of action potentials into a low-dimensional feature space for analysis.

Methods and ResultsOur automatic segmentation algorithm finds regions of acceptable action potentials in large data sets of electrophysiological readings. We use spectral methods and support vector machines to classify our readings and to extract relevant features. We are able to show that action potentials from the same cell site can be recorded over days without detrimental effects to the cell membrane. The variability between measurements 24 h apart is comparable to the natural variability of the features at a single time point.

ConclusionsOur work contributes towards a non-invasive approach for cardiomyocyte functional maturation, as well as developmental, pathological and pharmacological studies. As the human-derived cardiac model tissue has the genetic makeup of its donor, a powerful tool for individual drug toxicity screening emerges.

category: physiology


title: Absence of postsynaptic activity on developing neurons alters gene expression profiles without preventing the refinement of inputs

abstract: Synaptic activity plays several roles as developing neurons make connections with their targets. It acts locally at synapses to influence the expression of genes needed to establish and maintain synaptic contacts. And, downstream it provides the necessary activity to strengthen and refine connections. Many studies have demonstrated how synaptic activity alters synaptic strength and increases synapse numbers. Much less is known, however, about the long-term consequences when a circuit develops without synaptic activity. To address this, we developed a mosaic model of sympathetic ganglia where synaptically-active and synaptically-inactive sympathetic neurons develop side-by-side in vivo. This model allowed us to address two issues. One is the relationship between activity and the refinement of converging inputs; the second is how synaptic activity contributes to a neurons gene expression profile. Our results indicate that converging presynaptic inputs to synaptically-silent neurons do not require postsynaptic activity to refine, provided these neurons share targets with synaptically-active neurons. Second, we demonstrate with single-cell RNA sequencing experiments that the expression of many genes by sympathetic neurons is independent of endogenous activity or local signals immediately downstream of excitatory postsynaptic potentials. An exception are genes required for neurotransmitter metabolism: We found that for a large sub-population of sympathetic neurons, synaptic activity increases the expression of adrenergic genes and supresses the expression of cholinergic genes. We conclude that signals generated locally at synapses do not initiate refinement of converging inputs, and that synaptic activitys influence on a neurons gene expression profiles is complex and depends on context.


title: Histone H2A ubiquitination resulting from Brap loss of function connects multiple aging hallmarks and accelerates neurodegeneration

abstract: Aging is an intricate process that is characterized by multiple hallmarks including stem cell exhaustion, genome instability, epigenome alteration, impaired proteostasis, and cellular senescence. While each of these traits is detrimental at the cellular level, it remains unclear how they are interconnected to cause systemic organ deterioration. Here we show that abrogating Brap, a BRCA1 associated protein important for neurogenesis, results in cellular senescence with persistent DNA double-strand breaks and elevation of histone H2A mono- and poly-ubiquitination (H2Aub). The high H2Aub initiates histone proteolysis, leading to both epigenetic alteration and proteasome overflow. These defects induce neuroinflammation, impair proteostasis, accelerate neurodegeneration, and substantially shorten lifespan in mice carrying Brap deletions in the brain. We further show H2Aub is also increased in human brain tissues of Alzheimers disease. These data together suggest that chromatin aberrations mediated by H2Aub act as a nexus of multiple aging hallmarks and promote tissue-wide degeneration.


title: Circular RNAs exhibit limited evidence for translation, or translation regulation of the mRNA-counterpart in terminal hematopoiesis

abstract: Each day, about 1012 erythrocytes and platelets are released into the blood stream. This substantial output from hematopoietic stem cells is tightly regulated by transcriptional and epigenetic factors. Whether and how circular RNAs (circRNAs) contribute to the differentiation and/or identity of hematopoietic cells is to date not known. We recently reported that erythrocytes and platelets contain the highest levels and numbers of circRNAs amongst hematopoietic cells. Here, we provide the first detailed analysis of circRNA expression during erythroid and megakaryoid differentiation. CircRNA expression not only significantly increased upon enucleation, but also had limited overlap between progenitor cells and mature cells, suggesting that circRNA expression stems from regulated processes rather than resulting from mere accumulation. To study circRNA function in hematopoiesis, we first compared the expression levels of circRNAs with the translation efficiency of their mRNA-counterpart. We found that only 1 out of 2531 (0.04%) circRNAs associated with mRNA-translation regulation. Furthermore, irrespective of 1000s of identified putative open reading frames, deep ribosome-footprinting sequencing and mass spectrometry analysis provided little evidence for translation of endogenously expressed circRNAs. In conclusion, circRNAs alter their expression profile during terminal hematopoietic differentiation, yet their contribution to regulate cellular processes remains enigmatic.

category: molecular biology


title: Evaluation of pigmentation gene knock-outs by CRISPR as potential targets for the genetic engineering of Culex quinquefasciatus

abstract: Culex quinquefasciatus mosquitoes are a globally widespread vector of several human and animal pathogens. Their biology and behavior allow them to thrive in proximity to urban areas, rendering them a constant public health threat. Their mixed bird/mammal feeding behavior further offers a vehicle for zoonotic pathogens transmission to people, and separately, poses a threat to the conservation of insular birds. The advent of CRISPR has led to the development of novel technologies for the genetic engineering of wild mosquito populations, yet research in Culex quinquefasciatus has been lagging compared to other disease vectors. Here we use this tool to disrupt a set of five pigmentation genes in Culex quinquefasciatus that, when altered, lead to visible, homozygous-viable phenotypes. We further validate this approach in separate laboratories and in two distinct strains of Culex quinquefasciatus that are relevant to potential future public health and bird conservation applications. We generate a double-mutant line, demonstrating the possibility of sequentially combining multiple such mutations in a single individual. Lastly, we target two loci, doublesex in the sex-determination pathway and proboscipedia a hox gene, demonstrating the flexibility of these methods applied to novel targets. Our work provides a platform of seven validated loci that could be used for targeted mutagenesis in Culex quinquefasciatus and the future development of genetic suppression strategies for this species. Furthermore, the mutant lines generated here could have widespread utility to the research community using this model organism, as they could be used as targets for transgene delivery, where a copy of the disrupted gene could be included as an easily-scored transgenesis marker.


title: A computationally tractable birth-death model that combines phylogenetic and epidemiological data

abstract:
Inferring the dynamics of pathogen transmission during an outbreak is an important problem in both infectious disease epidemiology. In mathematical epidemiology, estimates are often informed by time series of confirmed cases, while in phylodynamics genetic sequences of the pathogen, sampled through time, are the primary data source. Each data type provides different, and potentially complementary, insight; recent studies have recognised that combining data sources can improve estimates of the transmission rate and number of infected individuals. However, inference methods are typically highly specialised and field-specific and are either computationally prohibitive or require intensive simulation, limiting their real-time utility.

We present a novel birth-death phylogenetic model and derive a tractable analytic approximation of its likelihood, the computational complexity of which is linear in the size of the dataset. This approach combines epidemiological and phylodynamic data to produce estimates of key parameters of transmission dynamics and the number of unreported infections. Using simulated data we show (a) that the approximation agrees well with existing methods, (b) validate the claim of linear complexity and (c) explore robustness to model misspecification. This approximation facilitates inference on large datasets, which is increasingly important as large genomic sequence datasets become commonplace.

Author summaryMathematical epidemiologists typically studies time series of cases, ie the epidemic curve, to understand the spread of pathogens. Genetic epidemiologists study similar problems but do so using measurements of the genetic sequence of the pathogen which also contain information about the transmission process. There have been many attempts to unite these approaches so that both data sources can be utilised. However, striking a suitable balance between model flexibility and fidelity, in a way that is computationally tractable, has proven challenging; there are several competing methods but for large datasets they are intractable. As sequencing of pathogen genomes becomes more common, and an increasing amount of epidemiological data is collected, this situation will only be exacerbated. To bridge the gap between the time series and genomic methods we developed an approximation scheme, called TimTam, which can accurately and efficiently estimate key features of an epidemic such as the prevalence of the infection and the effective reproduction number, ie how many people are currently infected and the degree to which the infection is spreading.

category: genomics


title: Polyploidy-associated paramutation in Arabidopsis is determined by small RNAs, temperature, and allele structure

abstract:
Paramutation is a form of non-Mendelian inheritance in which the expression of a paramutable allele changes when it encounters a paramutagenic allele. This change in expression of the paramutable alleles is stably inherited even after segregation of both alleles. While the discovery of paramutation and studies of its underlying mechanism were made with alleles that change plant pigmentation, paramutation-like phenomena are known to modulate the expression of other traits and in other eukaryotes, and many cases have probably gone undetected. It is likely that epigenetic mechanisms are responsible for the phenomenon, as paramutation forms epialleles, genes with identical sequences but different expression states. This could account for the intergenerational inheritance of the paramutated allele, providing profound evidence that triggered epigenetic changes can be maintained over generations. Here, we use a case of paramutation that affects a transgenic selection reporter gene in tetraploid Arabidopsis thaliana. Our data suggest that different types of small RNA are derived from paramutable and paramutagenic epialleles. In addition, deletion of a repeat within the epiallele changes its paramutability. Further, the temperature during the growth of the epiallelic hybrids determines the degree and timing of the allelic interaction. The data further make it plausible why paramutation in this system becomes evident only in the segregating F2 population of tetraploid plants containing both epialleles. In summary, the results support a model for polyploidy-associated paramutation, with similarities as well as distinctions from other cases of paramutation.

AUTHOR SUMMARYIn 1866, Gregor Mendel formulated the general principles of inheritance based on crossing experiments with pea plants. Curiously, in 1915, the progeny from crossing pea plants with a regular and a "rogue" leaf phenotype was lacking the expected segregation and recovery of the regular phenotype. This discovery was one of the first observations of non-Mendelian genetics and later demonstrated for more traits in other plants and termed paramutation. Paramutation is due to the epigenetic switch of an active gene to a silenced version which is then maintained in the inactive state in later generations. This demonstrates that acquired epigenetic changes can become permanent. Despite its early observation and numerous studies, mainly in maize and tomato, it is barely understood how paramutation is established and which parameters influence the process. We investigated a case of paramutation in Arabidopsis thaliana, crossing plants with genetically identical but epigenetically different alleles that result in resistance or sensitivity to an antibiotic in the growth medium. Paramutation did not become manifest immediately but only in the progeny of the hybrids, and only in plants with a doubled chromosome set. These features make this paramutation distinct from other cases. Our studies revealed several parameters that influence paramutation: an important role for sRNAs to initiate silencing, the sequence of the allele itself, the environmental conditions during growth of the hybrids, the developmental stage, and the copy number ratio between the alleles.


title: Mapping the organization and dynamics of the posterior medial network during movie watching

abstract: Brain regions within a posterior medial network (PMN) are characterized by sensitivity to episodic tasks, and they also demonstrate strong functional connectivity as part of the default network. Despite its cohesive structure, delineating the intranetwork organization and functional diversity of the PMN is crucial for understanding its contributions to multidimensional event cognition. Here, we probed functional connectivity of the PMN during movie watching to identify its pattern of connections and subnetwork functions in a split-sample replication of 136 participants. Consistent with prior findings of default network fractionation, we identified distinct PMN subsystems: a Ventral PM subsystem (retrosplenial cortex, parahippocampal cortex, posterior angular gyrus) and a Dorsal PM subsystem (medial prefrontal cortex, hippocampus, precuneus, posterior cingulate cortex, anterior angular gyrus). Ventral and Dorsal PM subsystems were differentiated by functional connectivity with parahippocampal cortex and precuneus and integrated by retrosplenial cortex and posterior cingulate cortex, respectively. Finally, the distinction between PMN subsystems is functionally relevant: whereas both Dorsal and Ventral PM connectivity tracked the movie content, only Ventral PM connections increased in strength at event transitions and appeared sensitive to episodic memory. Overall, these findings reveal PMN functional pathways and the distinct functional roles of intranetwork subsystems during event cognition.


title: PEPATAC: An optimized pipeline for ATAC-seq data analysis with serial alignments

abstract:
MotivationAs chromatin accessibility data from ATAC-seq experiments continues to expand, there is continuing need for standardized analysis pipelines. Here, we present PEPATAC, an ATAC-seq pipeline that is easily applied to ATAC-seq projects of any size, from one-off experiments to large-scale sequencing projects.

ResultsPEPATAC leverages unique features of ATAC-seq data to optimize for speed and accuracy, and it provides several unique analytical approaches. Output includes convenient quality control plots, summary statistics, and a variety of generally useful data formats to set the groundwork for subsequent project-specific data analysis. Downstream analysis is simplified by a standard definition format, modularity of components, and metadata APIs in R and Python. It is restartable, fault-tolerant, and can be run on local hardware, using any cluster resource manager, or in provided Linux containers. We also demonstrate the advantage of aligning to the mitochondrial genome serially, which improves the accuracy of alignment statistics and quality control metrics. PEPATAC is a robust and portable first step for any ATAC-seq project.

AvailabilityBSD2-licensed code and documentation at https://pepatac.databio.org.


title: Cuttlefish: Fast, parallel, and low-memory compaction of de Bruijn graphs from large-scale genome collections

abstract:
MotivationThe construction of the compacted de Bruijn graph from collections of reference genomes is a task of increasing interest in genomic analyses. These graphs are increasingly used as sequence indices for short and long read alignment. Also, as we sequence and assemble a greater diversity of genomes, the colored compacted de Bruijn graph is being used as the basis for efficient methods to perform comparative genomic analyses on these genomes. Therefore, designing time and memory efficient algorithms for the construction of this graph from reference sequences is an important problem.

ResultsWe introduce a new algorithm, implemented in the toolCuttlefish, to construct the (colored) compacted de Bruijn graph from a collection of one or more genome references. Cuttlefish introduces a novel approach of modeling de Bruijn graph vertices as finite-state automata; it constrains these automatas state-space to enable tracking their transitioning states with very low memory usage. Cuttlefish is fast and highly parallelizable. Experimental results demonstrate that it scales much better than existing approaches, especially as the number and the scale of the input references grow. On our test hardware, Cuttlefish constructed the graph for 100 human genomes in under 9 hours, using ~29 GB of memory while no other tested tool completed this task. On 11 diverse conifer genomes, the compacted graph was constructed by Cuttlefish in under 9 hours, using ~84 GB of memory, while the only other tested tool that completed this construction on our hardware took over 16 hours and ~289 GB of memory.

AvailabilityCuttlefish is written in C++14, and is available under an open source license at https://github.com/COMBINE-lab/cuttlefish.

Contactrob@cs.umd.edu

Supplementary informationSupplementary text are available at Bioinformatics online.


title: A 9-kDa matricellular SPARC fragment released by cathepsin D exhibits pro-tumor activity in the triple-negative breast cancer microenvironment

abstract: Extracellular matrix (ECM) remodeling by proteases results in the release of protein fragments that promote tumor progression and metastasis. The protease cathepsin D (cath-D), a marker of poor prognosis in triple-negative breast cancer (TNBC), is aberrantly secreted in the tumor microenvironment. Using degradomic analyses by TAILS, we discovered that the matricellular protein SPARC is a substrate of extracellular cath-D. In vitro, cath-D induced limited proteolysis of SPARC C-terminal extracellular Ca2+ binding domain at acidic pH, leading to the production of SPARC fragments (34-, 27-, 16-, 9-, and 6-kDa). Similarly, cath-D secreted by TNBC cells cleaved fibroblast- and cancer cell-derived SPARC at the tumor pericellular acidic pH. SPARC cleavage also occurred in TNBC tumors. Among these fragments, only the 9-kDa SPARC fragment inhibited TNBC cell adhesion and spreading on fibronectin, and stimulated their migration, endothelial transmigration, and invasion. Our study establishes a novel crosstalk between proteases and matricellular proteins in the tumor microenvironment through limited proteolysis of SPARC, revealing a novel targetable 9-kDa bioactive SPARC fragment for new TNBC treatments.

category: cancer biology


title: Lymphatic PD-L1 expression restricts tumor-specific CD8+ T cell responses

abstract: Lymph node (LN)-resident lymphatic endothelial cells (LECs) mediate peripheral tolerance by self-antigen presentation on MHC-I and constitutive expression of T cell inhibitory molecules, including PD-L1. Tumor-associated LECs also upregulate PD-L1 but the specific role of lymphatic PD-L1 in tumor immunity is not well understood. We generated a mouse model lacking lymphatic PD-L1 expression and challenged these mice with two orthotopic tumor models, B16F10 melanoma and MC38 colorectal carcinoma. Lymphatic PD-L1 deficiency resulted in a consistent expansion of tumor-specific CD8+ T cells in tumor-draining LNs in both tumor models, reduced primary tumor growth in the MC38 model, and increased the efficacy of adoptive T cell therapy in the B16F10 model. Strikingly, lymphatic PD-L1 primarily acted via apoptosis induction in tumor-specific CD8+ central memory T cells. Our findings demonstrate that LECs restrain tumor-specific immunity via PD-L1 and may explain why some cancer patients without PD-L1 expression in the tumor microenvironment still respond to PD-L1 / PD-1 targeting immunotherapy.

category: immunology


title: A Gut Microbial Peptide and Molecular Mimicry in the Pathogenesis of Type 1 Diabetes

abstract:
Type 1 Diabetes (T1D) is an autoimmune disease characterized by the destruction of pancreatic {beta}-cells. One of the earliest aspects of this process is development of autoantibodies and T-cells directed at an epitope in the B-chain of insulin (insB:9-23). Analysis of microbial protein sequences with homology to insB:9-23 sequence revealed 17 peptides showing >50% identity to insB:9-23. Of these, one peptide, found in the normal human gut commensal Parabacteroides distasonis, activated both human T cell clones from T1D patients and T-cell hybridomas from non-obese diabetic (NOD) mice specific to insB:9-23. Immunization of NOD mice with P. distasonis insB:9-23 peptide mimic or insB:9-23 peptide verified immune cross-reactivity. Colonization of female NOD mice with P. distasonis accelerated the development of T1D, increasing macrophages, dendritic cells and destructive CD8+ T-cells, while decreasing FoxP3+ regulatory T-cells. Western blot analysis identified P. distasonis reacting antibodies in sera of NOD mice colonized with P. distasonis and human T1D patients. Furthermore, adoptive transfer of splenocytes from P. distasonis treated mice to NOD/SCID mice enhanced disease phenotype in the recipients. Finally, analysis of human infant gut microbiome data revealed that exposure of infants to P. distasonis may modulate disease pathogenesis. Taken together, these data demonstrate the potential role for an insB:9-23-mimimetic peptide from gut microbiota as a molecular trigger or modifier of T1D pathogenesis.

SIGNIFICANCE STATEMENTIn Type 1 diabetes (T1D), immune cells destroy pancreatic {beta}-cells. The trigger of this response, however, is unknown. Some sequences (epitopes) in the insulin molecule form a major target for this autoimmune response. We have identified a sequence in a human gut bacterium that can mimetic this insulin epitope. Immune cells specific to insulin cross-react with this bacterial mimetic. Further, this bacterium can accelerate diabetes onset in a mouse model of T1D, inducing destructive and decreasing protective immune cells. We found this mimetic in the gut of children developing T1D. Furthermore, T1D patients have a stronger immune response to this bacterium compared to healthy individuals. Taken together, this bacterial mimetic in human gut has the potential to trigger/modify T1D onset.


title: Thermally tolerant symbionts may explain Caribbean octocoral resilience to heat stress

abstract: Coral reef ecosystems are under threat from the frequent and severe impacts of anthropogenic climate change, particularly rising sea surface temperatures. The effects of thermal stress may be ameliorated by adaptation and/or acclimation of the host, symbiont, or holobiont (host + symbiont) to increased temperatures. We examined the role of the symbiont in promoting thermal tolerance of the holobiont, using Antillogorgia bipinnata (octocoral host) and Breviolum antillogorgium (symbiont) as a model system. We identified five distinct genotypes of B. antillogorgium from symbiont populations isolated from Antillogorgia colonies in the Florida Keys. Three symbiont genotypes were cultured and maintained at 26{degrees}C (ambient historical temperature) and two were cultured and maintained at 30{degrees}C (elevated historical temperature) for two years. Following culturing, we analyzed the growth rate and carrying capacity of each symbiont genotype at both ambient and elevated temperatures in culture (in vitro). All genotypes grew well at both temperatures, indicating thermal tolerance among these B. antillogorgium cultures. Prior culturing at the elevated temperature, however, did not result in increased thermal tolerance. We then inoculated juvenile A. bipinnata polyps with each of the five symbiont genotypes, and reared these polyps at both ambient and elevated temperatures (in hospite experiment). All genotypes were able to establish symbioses with polyps in both temperature treatments. Survivorship of polyps at 30{degrees}C was significantly lower than survivorship at 26{degrees}C, but all treatments had surviving polyps at 56 days post-infection, suggestive of broad thermal tolerance in B. antillogorgium, which may play a part in the increased resilience of Caribbean octocorals during heat stress events.

category: evolutionary biology


title: Asymmetric directed functional connectivity within the frontoparietal motor network during motor imagery and execution

abstract: Both imagery and execution of motor controls consist of interactions within a neuronal network, including frontal motor-related regions and posterior parietal regions. To reveal neural representation in the frontoparietal motor network, two approaches have been proposed thus far: one is decoding of actions/modes related to motor control from the spatial pattern of brain activity; another is to estimate directed functional connectivity, which means a directed association between two brain regions within motor regions. However, directed connectivity among multiple regions of the motor network during motor imagery (MI) or motor execution (ME) has not been investigated. Here, we attempted to characterize the directed functional connectivity within the frontoparietal motor-related networks between the MI and ME conditions. We developed a delayed sequential movement and imagery task to evoke brain activity associated with data under ME and MI via functional magnetic resonance imaging scanning. We applied a causal discovery approach, linear non-Gaussian acyclic causal model, to identify directed functional connectivity among the frontoparietal motor-related brain regions for each condition. We demonstrated higher directed functional connectivity from the contralateral dorsal premotor cortex (dPMC) to the primary motor cortex (M1) in ME than in MI. We mainly identified significant direct effects of the dPMC and ventral premotor cortex (vPMC) to the parietal regions. In particular, connectivity from the dPMC to the superior parietal lobule in the same hemisphere showed significant positive effects across all conditions. Contrastingly, interlateral connectivities from the vPMC to the superior parietal lobule showed significantly negative effects across all conditions. Finally, we found positive effects from A1 to M1 in the same hemisphere, such as the audio-motor pathway. These results indicated that the sources of motor command originated in d/vPMC influenced M1 and parietal regions as achieving ME and MI. Additionally, sequential sounds may functionally facilitate temporal motor processes.


title: PSK signaling controls ABA homeostasis and signaling genes and maintains shoot growth under osmotic stress

abstract:
Water deficit impairs growth and survival of plants. Many water stress responses are under control of abscisic acid (ABA) but little is known about growth control under osmotic stress. Based on the previously described growth-promoting activity of the peptide hormone phytosulfokine (PSK), we hypothesized that it may contribute to growth regulation under water stress conditions. To test this hypothesis, we analyzed the Arabidopsis thaliana PSK receptor (PSKR) null mutant pskr1-3 pskr2-1 under mannitol and drought stress. In particular under mild water stress, fresht weight and photosynthetic efficiency were more reduced in pskr1-3 pskr2-1 than in wild type. Hydroponic and grafting experiments showed that PSKR signaling was not required for long-distance signaling from mannitol-stressed roots to shoot but rather for cell growth promotion in the shoot. Unlike wild type, pskr1-3 pskr2-1 shoots did not accumulate ABA in response to mannitol, showed misregulation of ABA synthesis genes and elevated expression of ABI1 and ABI2, repressors of ABA signaling whereas application of ABA partially reversed shoot growth inhibition by mannitol in pskr1-3 pskr2-1. In turn, mannitol and ABA induced expression of PSK3 and PSKR1, and ABA promoted expression of PSK2 and PSK4 revealing feedback regulatory loops between PSKR and osmotic stress signaling.

HighlightPhytosulfokine receptor signaling regulates ABA synthesis and signaling genes and promotes ABA accumulation in the shoot of water-stressed plants and maintains leaf growth and photosynthetic efficiency which ensures plant health.

category: plant biology


title: In vitro gut modeling as a tool for adaptive evolutionary engineering of Lactiplantibacillus plantarum

abstract:
Research and marketing of probiotics demand for holistic strain improvement considering both, the biotic and abiotic gut environment. Here we aim to establish the continuous In vitro colonic fermentation model PolyFermS as a tool for adaptive evolutionary engineering. Immobilized fecal microbiota from adult donors were steadily cultivated up to 72 days in PolyFermS reactors, providing a long-term compositional and functional stable ecosystem akin to the donors gut. Inoculation of the gut microbiota with immobilized or planktonic Lactiplantibacillus plantarum NZ3400, a derivative of the probiotic model strain WCFS1, led to successful colonization. Whole genome sequencing of 45 recovered strains revealed mutations in 16 genes involved in signaling, metabolism, transport, and cell surface. Remarkably, mutations in LP_RS14990, LP_RS15205, and intergenic region LP_RS05100<LP_RS05095 were found in recovered strains from different adaptation experiments. Combined addition of the reference strain NZ3400 and each of those mutants to the gut microbiota resulted in increased abundance of the corresponding mutant in PolyFermS microbiota after 10 days, showing the beneficial nature of these mutations. Our data show that the PolyFermS system is a suitable technology to generate adapted mutants for colonization in colonic conditions. Analysis thereof will provide knowledge about factors involved in gut microbiota colonization and persistence.

ImportanceImprovement of bacterial strains to specific abiotic environmental factors is broadly used to enhance strain characteristics for processing and product quality. However, there is currently no multidimensional probiotic strain improvement approach for both, abiotic and biotic factors of a colon microbiota. The continuous PolyFermS fermentation model allows stable and reproducible continuous cultivation of colonic microbiota and provides conditions akin to the host gut with high control and easy sampling. This study investigated the suitability of PolyFermS for adaptive evolutionary engineering of a probiotic model organism for lactobacilli, Lactiplantibacillus plantarum, to an adult human colonic microbiota. The application of PolyFermS controlled gut microbiota environment led to adaptive evolution of L. plantarum strains for enhanced gut colonization characteristics. This novel tool for strain improvement can be used to reveal relevant factors involved in gut microbiota colonization and develop adapted probiotic strains with enhanced functionality in the gut.


title: Atlas of tissue-specific and tissue-preferential gene expression in ecologically and economically significant conifer Pinus sylvestris

abstract:
Despite their ecological and economical importance, conifers genomic resources are limited, mainly due to the large size and complexity of their genomes. Additionally, the available genomic resources lack complete structural and functional annotation. Transcriptomic resources have been commonly used to compensate for these deficiencies, though for most conifer species they are limited to a small number of tissues, or capture only a fraction of the genes present in the genome.

Here we provide an atlas of gene expression patterns for conifer Pinus sylvestris across five tissues: embryo, megagametophyte, needle, phloem, and vegetative bud. We used a wide range of tissues and focused our analyses on the expression profiles of genes at tissue level. We provide comprehensive information of the per-tissue normalized expression level, indication of tissue preferential upregulation and tissue-specificity of expression. We identified a total of 48,001 tissue preferentially upregulated and tissue specifically expressed genes, of which 28% have annotation in the Swiss-Prot database. Even though most of the putative genes identified do not have functional information in current biological databases, the tissue-specific patterns discovered provide valuable information about their potential functions for further studies, as for example in the areas of plant physiology, population genetics, and genomics in general. As we provide information on tissue specificity at both diploid and haploid life stages, our data will also contribute to the understanding of evolutionary rates of different tissue types and ploidy levels.


title: Quantifying growth perturbations over the fattening period in swine via mathematical modelling

abstract:
AO_SCPLOWBSTRACTC_SCPLOWO_ST_ABSBackgroundC_ST_ABSResilience can be defined as the capacity of animals to cope with short-term perturbations in their environment and return rapidly to their pre-challenge status. In a perspective of precision livestock farming, it is key to create informative indicators for general resilience and therefore incorporate this concept in breeding goals. In the modern swine breeding industry, new technologies such as automatic feeding system are increasingly common and can be used to capture useful data to monitor animal phenotypes such as feed efficiency. This automatic and longitudinal data collection integrated with mathematical modelling has a great potential to determine accurate resilience indicators, for example by measuring the deviation from expected production levels over a period of time.

ResultsThis work aimed at developing a modelling approach for facilitating the quantification of pig resilience during the fattening period, from approximately 34 kg to 105 kg of body weight. A total of 13 093 pigs, belonging to three different genetic lines were monitored (Pietrain, Pietrain NN and Duroc) since 2015, and body weight measures registered (approximately 11.1 million of weightings) with automatic feeding systems. We used the Gompertz model and linear interpolation on body weight data to quantify individual deviations from expected production, thereby creating a resilience index (ABC). The estimated heritabilities of ABC are low but not zero from 0.03 to 0.04 ({+/-} 0.01) depending on the breed.

ConclusionsOur model-based approach can be useful to quantify pig responses to perturbations using exclusively the growth curves and should contribute to the genetic improvement of resilience of fattening pigs by providing a resilience index.

category: systems biology


title: A new method to accurately identify single nucleotide variants using small FFPE breast samples.

abstract:
Most tissue collections of neoplasms are composed of formalin-fixed and paraffin-embedded (FFPE) excised tumor samples used for routine diagnostics. DNA sequencing is becoming increasingly important in cancer research and clinical management; however, it is difficult to accurately sequence DNA from FFPE samples. We developed and validated a new bioinformatic algorithm to robustly identify somatic single nucleotide variants (SNVs) from whole exome sequencing using small amounts of DNA extracted from archival FFPE samples of breast cancers. We optimized this strategy using 28 pairs of technical replicates. After optimization, the mean similarity between replicates increased 5-fold, reaching 88% (range 0-100%), with a mean of 21.4 SNVs (range 1-68) per sample, representing a markedly superior performance to existing algorithms. We found that the SNV-identification accuracy declined when there was less than 40ng of DNA available and that insertion-deletion variant calls are less reliable than single base substitutions. As the first application of the new algorithm, we compared samples of ductal carcinoma in situ (DCIS) of the breast to their adjacent invasive ductal carcinoma (IDC) samples. We observed an increased number of mutations (paired-samples sign test, p<0.05), and a higher genetic divergence in the invasive samples (paired-samples sign test, p<0.01). Our algorithm provides a significant improvement in detecting SNVs in FFPE samples over previous approaches.

Key PointsO_LIThe sequencing of reduced quantities of DNA extracted from FFPE samples leads to substantial sequencing errors that require correction in order to obtain accurate detection of somatic mutations.
C_LIO_LIWe developed and validated a new bioinformatic algorithm to robustly identify somatic single nucleotide variants using small amounts of DNA extracted from archival FFPE samples of breast cancers.
C_LIO_LIVariant calling software packages need to be optimized to reduce the impact of sequencing errors. Our bioinformatics pipeline represents a significant methodological advance compared to the currently available bioinformatic tools used for the analysis of small FFPE samples.
C_LI


title: Multivariate genome-wide association study identifies 780 unique genetic loci associated with cortical morphology.

abstract: Brain morphology has been shown to be highly heritable, yet only a small portion of the heritability is explained by the genetic variants discovered so far. Here we exploit the distributed nature of genetic effects across the brain and apply the Multivariate Omnibus Statistical Test (MOSTest) to genome-wide association studies (GWAS) of vertex-wise structural magnetic resonance imaging (MRI) cortical measures from N=35,657 participants in the UK Biobank. We identified 695 loci for cortical surface area and 539 for cortical thickness, in total 780 unique genetic loci associated with cortical morphology. This reflects an approximate 10-fold increase compared to the commonly applied univariate GWAS methods. Power analysis indicates that applying MOSTest to vertex-wise structural MRI data triples the effective sample size compared to conventional univariate GWAS approaches. Functional follow up including gene-based analyses implicate 10% of all protein-coding genes and point towards pathways involved in neurogenesis and cell differentiation.


title: Cell group analysis reveals changes in upper-layer neurons associated with schizophrenia

abstract: Genome-wide association studies (GWAS) of schizophrenia (SCZ) have revealed over 100 risk loci. We investigated whether these SCZ-associated variants regulate gene expression by cell type. Using a fully unsupervised deconvolution method, we calculated gene expression by clusters of estimated cell types (cell-groups, CGs). Five CGs emerged in the dorsolateral prefrontal cortices (DLPFC) of 341 donors with and without SCZ. By mapping expression quantitative trait loci (eQTL) per CG, we partitioned the heritability of SCZ risk in GWAS by CGs. CG-specific expressions and eQTLs were replicated in both a deconvoluted bulk tissue data set with a different method and also in sorted-cell expression data. Further, we characterized CG-specific gene differential expression and cell proportion changes in SCZ brains. We found upper-layer neurons in the DLPFC to be associated with SCZ based on enrichment of SCZ heritability in eQTLs, disease-related transcriptional signatures, and decreased cell proportion. Our study suggests that neurons and related anomalous circuits in the upper layers of the DLPFC may have a major contribution to SCZ risk.


title: mRNA based SARS-CoV-2 vaccine candidate CVnCoV induces high levels of virus neutralizing antibodies and mediates protection in rodents

abstract: mRNA technologies have recently proven clinical efficacy against coronavirus disease 2019 (COVID-19) and are among the most promising technologies to address the current pandemic. Here, we show preclinical data for our clinical candidate CVnCoV, a lipid nanoparticle encapsulated mRNA vaccine that encodes full length, pre-fusion stabilised severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) Spike protein. In contrast to previously published approaches, CVnCoV is exclusively composed of naturally occurring nucleotides. Immunisation with CVnCoV induced strong humoral responses with high titres of virus neutralizing antibodies and robust T cell responses. CVnCoV vaccination protected hamsters from challenge with wild type SARS-CoV-2, demonstrated by the absence of viral replication in the lungs. Hamsters vaccinated with a suboptimal dose of CVnCoV leading to breakthrough viral replication exhibited no evidence of vaccine enhanced disease. Overall, data presented here provide evidence that CVnCoV represents a potent and safe vaccine candidate against SARS-CoV-2.


title: Potential Achilles heels of SARS-CoV-2 are best displayed by the base order-dependent component of RNA folding energy

abstract: The base order-dependent component of folding energy has revealed a highly conserved region in HIV-1 genomes that associates with RNA structure. This corresponds to a packaging signal that is recognized by the nucleocapsid domain of the Gag polyprotein. Long viewed as a potential HIV-1 "Achilles heel," the signal can be targeted by a new antiviral compound. Although SARS-CoV-2 differs in many respects from HIV-1, the same technology displays regions with a high base order-dependent folding energy component, which are also highly conserved. This indicates structural invariance (SI) sustained by natural selection. While the regions are often also protein-encoding (e.g. NSP3, ORF3a), we suggest that their nucleic acid level functions can be considered potential "Achilles heels" for SARS-CoV-2, perhaps susceptible to therapies like those envisaged for AIDS. The ribosomal frameshifting element scored well, but higher SI scores were obtained in other regions, including those encoding NSP13 and the nucleocapsid (N) protein.


title: ATRX alteration contributes to tumor growth and immune escape in pleomorphic sarcomas

abstract: Whole genome and transcriptome sequencing of a cohort of 67 leiomyosarcomas revealed ATRX to be one of the most frequently mutated genes in leiomyosarcomas after TP53 and RB1. While its function is well described in the alternative lengthening of telomeres mechanism, we wondered whether its alteration could have complementary effects on sarcoma oncogenesis. ATRX alteration is associated with the down-expression of genes linked to differentiation in leiomyosarcomas, and to immunity in an additional cohort of 60 poorly differentiated sarcomas. In vitro and in vivo models showed that ATRX loss increases tumor growth rate and immune escape by decreasing the immunity load of active mast cells in sarcoma tumors. These data indicate that an alternative to unsuccessful targeting of the adaptive immune system in sarcoma could be to target the innate system. This might lead to a better outcome for sarcoma patients in terms of ATRX status.


title: A damaged informed lung ventilator model for ventilator waveforms

abstract: Motivated by desire to understand pulmonary physiology and pathophysiology, scientists have developed models of pulmonary physiology. However, pathophysiology and interactions between human lungs and ventilators, e.g., ventilator-induced lung injury (VILI), present problems for modeling efforts. Real-world injury is too complex for simple models to capture, and while complex models tend not to be estimable with clinical data, limiting both the clinical utility with existing approaches. To address this gap, we present a damaged-informed lung ventilator (DILV) model to model and quantify patient-ventilator interactions and lung health. This approach relies on systematically mathematizing the pathophysiologic knowledge clinicians use to interpret lung condition from ventilator waveform data. This is achieved by defining clinically relevant features in the ventilator waveform data that contain hypothesis-driven information about pulmonary physiology, patient-ventilator interaction, and ventilator settings. To capture these features, we develop a modelling framework where the model has enough flexibility to reproduce commonly observed variability in waveform data. We infer the model parameters with clinical (human) and laboratory (mouse) data. The DILV model can reproduce essential dynamics of differently damaged lungs for tightly controlled measurements in mice and uncontrolled human intensive care unit data in the absence and presence of respiratory effort. Estimated parameters correlate with known measures of lung physiology, including lung compliance. This method has the potential to translate laboratory physiology experiments to clinical applications, including pathways for high fidelity estimates of lung state and sources of VILI with an end goal of reducing the impact of VILI and acute respiratory distress syndrome.


title: The Oscillatory Effects of Rhythmic Median Nerve Stimulation

abstract: Entrainment of brain oscillations can be achieved using rhythmic non-invasive brain stimulation, and stimulation of the motor cortex at a frequency associated with sensorimotor inhibition can impair motor responses. Despite the potential for therapeutic application, these techniques do not lend themselves to use outside of a clinical setting. Here, the aim was to investigate whether rhythmic median nerve stimulation (MNS) could be used to entrain oscillations related to sensorimotor inhibition. MEG data were recorded from 20 participants during 400 trials, where for each trial 10 pulses of MNS were delivered either rhythmically or arrhythmically at 12 or 20Hz. Our results demonstrate a frequency specific increase in relative amplitude in the contralateral somatosensory cortex during rhythmic but not arrhythmic stimulation. This was coupled with an increase in inter-trial phase coherence at the same frequency, suggesting that the oscillations synchronised with the pulses of MNS. While the results show that 20Hz rhythmic peripheral nerve stimulation can produce entrainment, the response to 12Hz stimulation was largely due to the presence of rhythmic sensory evoked potentials. Regardless, rhythmic MNS resulted in synchronous firing of neuronal populations within the contralateral somatosensory cortex meaning these neurons were  occupied in processing of the afferent input. Therefore, MNS could prove therapeutically useful in disorders associated with hyperexcitability within the sensorimotor cortices.


title: A unified physiological framework of transitions between seizures, sustained ictal activity, and depolarization block at the single neuron level

abstract:
The majority of seizures recorded in humans and experimental animal models can be described by a generic phenomenological mathematical model, The Epileptor. In this model, seizure-like events (SLEs) are driven by a slow variable and occur via saddle node (SN) and homoclinic bifurcations at seizure onset and offset, respectively. Here we investigated SLEs at the single cell level using a biophysically relevant neuron model including a slow/fast system of four equations. The two equations for the slow subsystem describe ion concentration variations and the two equations of the fast subsystem delineate the electrophysiological activities of the neuron. Using extracellular K+ as a slow variable, we report that SLEs with SN/homoclinic bifurcations can readily occur at the single cell level when extracellular K+ reaches a critical value. In patients and experimental models, seizures can also evolve into sustained ictal activity (SIA) and, depolarization block (DB), activities which are also parts of the dynamic repertoire of the Epileptor. Increasing extracellular concentration of K+ in the model to values found during experimental status epilepticus and DB, we show that SIA and DB can also occur at the single cell level. Thus, seizures, SIA and DB, which have been first identified as network events, can exist in a unified framework of a biophysical model at the single neuron level and exhibit similar dynamics as observed in the Epileptor.

Author SummaryEpilepsy is a neurological disorder characterized by the occurrence of seizures. Seizures have been characterized in patients in experimental models at both macroscopic and microscopic scales using electrophysiological recordings. Experimental works allowed the establishment of a detailed taxonomy of seizures, which can be described by mathematical models. We can distinguish two main types of models. Phenomenological (generic) models have few parameters and variables and permit detailed dynamical studies often capturing a majority of activities observed in experimental conditions. But they also have abstract parameters, making biological interpretation difficult. Biophysical models, on the other hand, use a large number of variables and parameters due to the complexity of the biological systems they represent. Because of the multiplicity of solutions, it is difficult to extract general dynamical rules. In the present work, we integrate both approaches and reduce a detailed biophysical model to sufficiently low-dimensional equations, and thus maintaining the advantages of a generic model. We propose, at the single cell level, a unified framework of different pathological activities that are seizures, depolarization block, and sustained ictal activity.


title: Environmental enrichment enhances patterning and remodeling of synaptic nanoarchitecture revealed by STED nanoscopy

abstract: Synaptic plasticity underlies long-lasting structural and functional changes to brain circuitry and its experience-dependent remodeling can be fundamentally enhanced by environmental enrichment. It is unknown, however, whether and how environmental enrichment alters the morphology and dynamics of individual synapses. Here, we present a virtually crosstalk free, two-color in vivo STED microscope to simultaneously superresolve the dynamics of endogenous PSD95 of the post-synaptic density and spine geometry. With environmental enrichment, the size distributions of PSD95 and spine head sizes were sharper than in controls, indicating that synaptic strength is set more precisely with environmental enrichment. Spine head geometry and PSD95 assemblies were highly dynamic but their changes correlated only mildly. With environmental enrichment, the topography of the PSD95 nanoorganization was more dynamic; changes in size were smaller than in mice housed in standard cages and depended linearly on their original size. Thus, two-color in vivo time-lapse imaging of synaptic nanoorganization uncovers a unique synaptic nanoplasticity associated with the enhanced learning capabilities under environmental enrichment.


title: Limited and strain-specific transcriptional and growth responses to acquisition of a multidrug resistance plasmid in genetically diverse Escherichia coli lineages

abstract:
Multi-drug resistant (MDR) Escherichia coli are a major global threat to human health, wherein multi-drug resistance is primarily spread by MDR plasmid acquisition. MDR plasmids are not widely distributed across the entire E. coli species, but instead are concentrated in a small number of clones. Here, we test if diverse E. coli strains vary in their ability to acquire and maintain MDR plasmids, and if this relates to their transcriptional response following plasmid acquisition. We used strains from across the diversity of E. coli, including the common MDR lineage ST131, and the IncF plasmid, pLL35, encoding multiple antibiotic resistance genes. Strains varied in their ability to acquire pLL35 by conjugation, but all were able to stably maintain the plasmid. The effects of pLL35 acquisition on cefotaxime resistance and growth also varied among strains, with growth responses ranging from a small decrease to a small increase in growth of the plasmid-carrier relative to the parental strain. Transcriptional responses to pLL35 acquisition were limited in scale and highly strain specific. We observed significant transcriptional responses at the operon or regulon level, possibly due to stress responses or interactions with resident MGEs. Subtle transcriptional responses consistent across all strains were observed affecting functions, such as anaerobic metabolism, previously shown to be under negative frequency dependent selection in MDR E. coli. Overall, there was no correlation between the magnitude of the transcriptional and growth responses across strains. Together these data suggest that fitness costs arising from transcriptional disruption are unlikely to act as a barrier to dissemination of this MDR plasmid in E. coli.

ImportancePlasmids play a key role in bacterial evolution by transferring niche adaptive functions between lineages, including driving the spread of antibiotic resistance genes. Fitness costs of plasmid acquisition arising from the disruption of cellular processes could limit the spread of multidrug resistance plasmids. However, the impacts of plasmid acquisition are typically measured in lab-adapted strains rather than in more ecologically relevant natural isolates. Using a clinical multidrug resistance plasmid and a diverse collection of E. coli strains isolated from clinical infections and natural environments, we show that plasmid acquisition had only limited and highly strain-specific effects on bacterial growth and transcription. These findings suggest that fitness costs arising from transcriptional disruption are unlikely to act as a barrier to transmission of this plasmid in natural populations of E. coli.


title: Structural classification of proteins based on the computationally efficient recurrence quantification analysis and horizontal visibility graphs

abstract:
MotivationProtein structure prediction is one of the most significant problems in bioinformatics, as it has a prominent role in understanding the function and evolution of proteins. Designing a computationally efficient but at the same time accurate prediction method remains a pressing issue, especially for sequences that we cannot obtain a sufficient amount of homologous information from existing protein sequence databases. Several studies demonstrate the potential of utilizing chaos game representation (CGR) along with time series analysis tools such as recurrence quantification analysis (RQA), complex networks, horizontal visibility graphs (HVG) and others. However, the majority of existing works involve a large amount of features and they require an exhaustive, time consuming search of the optimal parameters. To address the aforementioned problems, this work adopts the generalized multidimensional recurrence quantification analysis (GmdRQA) as an efficient tool that enables to process concurrently a multidimensional time series and reduce the number of features. In addition, two data-driven algorithms, namely average mutual information (AMI) and false nearest neighbors (FNN), are utilized to define in a fast yet precise manner the optimal GmdRQA parameters.

ResultsThe classification accuracy is improved by the combination of GmdRQA with the HVG. Experimental evaluation on a real benchmark dataset demonstrates that our methods achieve similar performance with the state-of-the-art but with a smaller computational cost.

AvailabilityThe code to reproduce all the results is available at https://github.com/aretiz/protein_structure_classification/tree/main.

Contactedoutsi@ics.forth.gr

Supplementary informationSupplementary data are available at Bioinformatics online.


title: In vivo visualization of propagating α-synuclein pathologies in mouse and marmoset models by a bimodal imaging probe, C05-05

abstract: Deposition of intracellular -synuclein fibrils is implicated in neurodegenerative parkinsonian disorders, while high-contrast in vivo detection of -synuclein depositions has been unsuccessful in animal models and humans. Here, we have developed a bimodal imaging probe, C05-05, for visualizing -synuclein inclusions in the brains of living animals modeling -synuclein propagation. In vivo optical and PET imaging of a mouse model demonstrated visualization of -synuclein aggregates by C05-05, revealing a dynamic propagation of fibrillogenesis along neural pathways followed by disruptions of these structures. Moreover, longitudinal 18F-C05-05-PET of a marmoset model captured widespread dissemination of fibrillary pathologies accompanied by neurodegeneration detected by dopamine transporter PET. In addition, in vitro assays demonstrated the high-affinity binding of 18F-C05-05 to -synuclein versus other protein pathologies in human brain tissues. Collectively, we propose a new imaging technology enabling etiological and therapeutic assessments of -synuclein pathogenesis at non-clinical levels, highlighting the applicability of C05-05 to clinical PET.


title: Dynamics of sex chromosome evolution in a rapid radiation of cichlid fishes

abstract: Sex is a fundamental trait that is determined, depending on the species, by different environmental and/or genetic factors, including various types of sex chromosomes. While the functioning and emergence of sex chromosomes have been explored in species scattered across the eukaryotic tree of life, little is known about tempo and mode of sex chromosome evolution in closely related species. Here, we examine the dynamics of sex chromosome evolution in an archetypical example of adaptive radiation, the cichlid fishes of African Lake Tanganyika. Through inspection of male and female genomes from 244 cichlid taxa and the analysis of transcriptomes from 66 taxa, we identify signatures of sex chromosomes in 79 taxa, involving 12 different linkage groups. We estimate that Tanganyikan cichlids have the highest rates of sex chromosome turnover and heterogamety transitions known to date. We further show that the recruitment of chromosomes as sex chromosomes is not at random and that some chromosomes have convergently emerged as sex chromosomes in cichlids, which provides empirical support to the "limited options" hypothesis of sex chromosome evolution.


title: Functional Molecular Switches of Mammalian G Protein-Coupled Bitter-Taste Receptors

abstract:
Bitter taste receptors (TAS2Rs) are a poorly understood subgroup of G protein-coupled receptors (GPCRs). The experimental structure of these receptors has yet to be determined, and key-residues controlling their function remain mostly unknown. We designed an integrative approach to improve comparative modeling of TAS2Rs. Using current knowledge on class A GPCRs and existing experimental data in the literature as constraints, we pinpointed conserved motifs to entirely re-align the amino-acid sequences of TAS2Rs. We constructed accurate homology models of human TAS2Rs. As a test case, we examined the accuracy of the TAS2R16 model with site-directed mutagenesis and in vitro functional assays. This combination of in silico and in vitro results clarify sequence-function relationships and identify the functional molecular switches that encode agonist sensing and downstream signaling mechanisms within mammalian TAS2Rs sequences.

ClassificationBiological sciences, Computational biology, and bioinformatics


title: homologizer: Phylogenetic phasing of gene copies into polyploid subgenomes

abstract:
O_LIOrganisms such as allopolyploids and F1 hybrids contain multiple distinct subgenomes, each potentially with its own evolutionary history. These organisms present a challenge for multilocus phylogenetic inference and other analyses since it is not apparent which gene copies from different loci are from the same subgenome and thus share an evolutionary history.
C_LIO_LIHere we introduce homologizer, a flexible Bayesian approach that uses a phylogenetic framework to infer the phasing of gene copies across loci into their respective subgenomes.
C_LIO_LIThrough the use of simulation tests we demonstrate that homologizer is robust to a wide range of factors, such as incomplete lineage sorting and the phylogenetic informativeness of loci. Furthermore, we establish the utility of homologizer on real data, by analyzing a multilocus dataset consisting of nine diploids and 19 tetraploids from the fern family Cystopteridaceae.
C_LIO_LIFinally, we describe how homologizer may potentially be used beyond its core phasing functionality to identify non-homologous sequences, such as hidden paralogs or contaminants.
C_LI


title: Differences in the gut microbiomes of distinct ethnicities within the same geographic area are linked to host metabolic health

abstract:
BackgroundThe human gut microbiota exhibits marked variation around the world, which has been attributed to dietary intake and other environmental factors. However, the degree to which ethnicity-associated differences in gut microbial community structure and function are maintained following immigration or in the context of metabolic disease is poorly understood.

ResultsWe conducted a multi-omic study of 46 lean and obese East Asian and White participants living in the San Francisco Bay Area. 16S rRNA gene sequencing revealed significant differences between ethnic groups in bacterial richness and community structure. White individuals were enriched for the mucin-degrading Akkermansia muciniphila. East Asian participants had increased levels of multiple bacterial phyla, fermentative pathways detected by metagenomics, and the short-chain fatty acid end products acetate, propionate, and isobutyrate. Differences in the gut microbiota between the East Asian and White groups could not be explained by reported dietary intake, were more pronounced in lean individuals, and were associated with current geographical location. Microbiome transplantations into germ-free mice confirmed that the differences in the gut microbiota of the East Asian and White individuals we analyzed are independent of diet and that they differentially impact host body weight and adiposity in genetically identical mouse recipients.

ConclusionsThe reported findings emphasize the utility of studying diverse ethnic groups within a defined geographical location and provide a starting point for dissecting the mechanisms contributing to the complex interactions between the gut microbiome and ethnicity-associated lifestyle, demographic, metabolic, and genetic factors.


title: Viral mimicry of p65/RelA transactivation domain to inhibit NF-κB activation

abstract: Sensing of virus infection activates NF-{kappa}B to induce the expression of interferons, cytokines and chemokines to initiate the antiviral response. Viruses antagonise these antiviral defences by interfering with immune sensing and blocking the actions of antiviral and inflammatory molecules. Here, we show that a viral protein mimics the transactivation domain of the p65 subunit of NF-{kappa}B. The C terminus of vaccinia virus (VACV) protein F14 (residues 51-73) activates transcription when fused to a DNA-binding domain-containing protein and F14 associates with NF-{kappa}B co-activator CBP, disrupting p65-CBP interaction. Consequently, F14 diminishes CBP-mediated acetylation of p65 and the downstream recruitment of the transcriptional regulator BRD4 to the promoter of the NF-{kappa}B-responsive genes CXCL10 and CCL2, hence inhibiting their expression. Conversely, the recruitment of BRD4 to the promoters of NFKBIA, which encodes the inhibitor of NF-{kappa}B (I{kappa}B), and CXCL8 remains unaffected in the presence of either F14 or JQ1, a competitive inhibitor of BRD4 bromodomains, indicating its recruitment is acetylation-independent. Therefore, unlike other viral NF-{kappa}B antagonists, F14 is a selective inhibitor of NF-{kappa}B-dependent gene expression. A VACV strain lacking F14 showed that it contributes to virulence in an intradermal model of infection. Our results uncover a mechanism by which viruses disarm the antiviral defences through molecular mimicry of a conserved host protein and provide insight into the regulation of NF-{kappa}B-dependent gene expression by BRD4.


title: SRF-MYOCD axis is the targetable driver of a well differentiated aggressive subtype of leiomyosarcomas

abstract: Differentiated phenotypes in solid tumors mostly lead to a better outcome than in their more immature counterparts. In leiomyosarcoma (LMS), a very aggressive tumor for which neither chemotherapy nor targeted or immune therapies have demonstrated efficacy, a relatively transcriptionally uniform subgroup of well differentiated tumors has been described and is associated with poor survival. In this disease, differentiation does not systematically hinder tumor progression, so it is a good model for studying the relation between these two mechanisms. By applying ad hoc bioinformatic methods, we conducted a comprehensive analysis in order to detect subgroups of leiomyosarcoma and finely investigate their transcriptomic and genomic profiles. We first used unsupervised hierarchical clustering to identify the most transcriptionally homogeneous LMS subgroup, which we named  hLMS. We derived a transcriptional signature that could consistently stratify patients with similar clinical and genomic characteristics in TCGA and ICGC cohorts. Functional analysis of the signature showed over-expression of SRF and E2Fs targets harboring histone modifications that were specifically identified in smooth muscles in hLMS. By integrating multi-omics data, we identified punctual and structural variations as well as the potential cellular origins of each group of LMS using GTEX normal tissue gene expression. Altogether, our results suggest that hLMS, but not oLMS (other LMS), may originate from vascular smooth muscle cells in which oncogenesis is driven by the MYOCD/SRF axis, which is essential for hLMS survival. We propose that the phenotypic plasticity of vascular smooth muscle cells coupled with MYOCD/SRF pathway overactivation could explain how hLMS balance this uncommon interplay between differentiation and aggressiveness.


title: Distinct cellular origins and differentiation process account for distinct oncogenic and clinical behaviors of leiomyosarcomas


title: Neural Mechanism of Hunger Gated Food-seeking and evaluating

abstract: The physiological need for energy evokes motivated feeding behaviours that help to ensure survival. However, the neural mechanisms underlying the generation of food motivation remain poorly understood. We investigated these mechanisms by subdividing feeding-related motivated behaviours into food-seeking, evaluating, and swallowing. Micro-endoscopic results indicated that neurons containing leptin receptors (LepRs) in the lateral hypothalamus (LH) are the major food-specific subpopulation of LHGABA neurons. Optogenetic manipulation of LHLepR neurons bidirectionally regulated both food-seeking and evaluating. Furthermore, micro-endoscope data revealed that distinct LHLepR neurons encode seeking and evaluating. Computational modelling analysis demonstrated that LHLepR neurons encode motivation, whereas neurons containing agouti-related peptide and neuropeptide Y (AgRP/NPY) encode the need for food. Additionally, slice studies revealed that NPY decreases inhibitory input to LHLepR neurons via LHGABA interneurons. This mechanism explains the permissive gate role of hunger (food need) in seeking/evaluating motivation. Together, the present study provides a comprehensive neural mechanism of how physiological needs drive distinct motivated behaviours.


title: Neural Mechanism of Hunger Gated Food-seeking and Evaluating


title: CSF-1 maintains pathogenic but not homeostatic myeloid cells in the central nervous system during autoimmune neuroinflammation.

abstract:
The receptor for colony stimulating factor 1 (CSF-1R) is important for the survival and function of myeloid cells that mediate pathology during experimental autoimmune encephalomyelitis (EAE), an animal model of multiple sclerosis (MS). CSF-1 and IL-34, the ligands of CSF-1R, have similar bioactivities but distinct tissue and context-dependent expression patterns, suggesting that they have different roles. This could be the case in EAE, given that CSF-1 expression is upregulated in the CNS, while IL-34 remains constitutively expressed. We found that targeting CSF-1 with neutralizing mAb halted ongoing EAE, with efficacy superior to CSF-1R inhibitor BLZ945, whereas IL-34 neutralization had no effect, suggesting that pathogenic myeloid cells were maintained by CSF-1, not IL-34. Both anti-CSF-1- and BLZ945-treatment greatly reduced numbers of monocyte-derived cells and microglia in the CNS. However, anti-CSF-1 selectively depleted inflammatory microglia and monocytes in inflamed CNS areas, whereas BLZ945 depleted virtually all myeloid cells, including quiescent microglia, throughout the CNS. Anti-CSF-1 treatments reduced the size of demyelinated lesions, and microglial activation in the grey matter. Lastly, we found that bone marrow-derived immune cells were the major mediators of CSF-1R-dependent pathology, while microglia played a lesser role. Our findings suggest that targeting CSF-1 could be effective in ameliorating MS pathology, while preserving the homeostatic functions of myeloid cells, thereby minimizing risks associated with total ablation of CSF-1R-dependent cells.

Significance StatementMultiple sclerosis (MS) and its animal model, experimental autoimmune encephalomyelitis (EAE), are autoimmune diseases characterized by accumulation of myeloid immune cells into the central nervous system (CNS). Both harmful and beneficial myeloid cells are present in EAE/MS, and a goal of MS therapy is to preferentially remove harmful myeloid cells. The receptor for CSF-1 (CSF-1R) is found on myeloid cells and it is important for their survival. CSF-1R can bind two ligands, CSF-1 and IL-34, but is unknown whether their functions in EAE/MS differ. We found that blocking CSF-1 depleted only harmful myeloid cells in the CNS and suppressed EAE, whereas blocking IL-34 had no effect. Thus, we propose that blocking CSF-1 could be a novel therapy for MS.


title: Speciation in the face of gene flow within the toothed whale superfamily Delphinoidea

abstract: The importance of post-divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide-ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae, Monodontidae). To investigate the role of post-divergence gene flow in the speciation process, we used a multifaceted approach, including: (i) phylogenomics, (ii) the distribution of shared derived alleles, and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding-window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D-statistics, D-foil, and f-branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggests that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species gene pools.


title: A genomic assessment of the marine-speciation paradox within the toothed whale superfamily Delphinoidea