Unifying Single-Channel Permeability from Rare-Event Sampling and Steady-State Flux Various all-atom molecular dynamics (MD) simulation methods have been developed to compute free energies and crossing rates of ions and small molecules through ion channels. However, a systemic comparison across different methods is scarce. Using a carbon nanotube as a model of small conductance ion channel, we computed the single-channel permeability for potassium ion using umbrella sampling, Markovian milestoning, and steady-state flux under applied voltage. We show that a slightly modified inhomogeneous solubility-diffusion equation yields a single-channel permeability consistent with the mean first passage time (MFPT)-based method. For milestoning, applying cylindrical and spherical bulk boundary conditions yield consistent MFPT if factoring in the effective bulk concentration. The sensitivity of the MFPT to the output frequency of collective variables is highlighted using the convergence and symmetricity of the inward and outward MFPT profiles. The consistent transport kinetic results from all three methods demonstrated the robustness of MD-based methods in computing ion channel permeation. The advantages and disadvantages of each technique are discussed, focusing on the future applications of milestoning in more complex systems.
Conformational dynamics and allosteric modulation of the SARS-CoV-2 spike Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects cells through binding to angiotensin-converting enzyme 2 (ACE2). This interaction is mediated by the receptor-binding domain (RBD) of the viral spike (S) glycoprotein. Structural and dynamic data have shown that S can adopt multiple conformations, which controls the exposure of the ACE2-binding site in the RBD. Here, using single-molecule Forster resonance energy transfer (smFRET) imaging we report the effects of ACE2 and antibody binding on the conformational dynamics of S from the Wuhan-1 strain and the B.1 variant (D614G). We find that D614G modulates the energetics of the RBD position in a manner similar to ACE2 binding. We also find that antibodies that target diverse epitopes, including those distal to the RBD, stabilize the RBD in a position competent for ACE2 binding. Parallel solution-based binding experiments using fluorescence correlation spectroscopy (FCS) indicate antibody-mediated enhancement of ACE2 binding. These findings inform on novel strategies for therapeutic antibody cocktails.
Non-invasive imaging of gene expression and protein secretion dynamics in living mice uncovers ectopic prothrombin expression as driver of thrombosis in cancer The topology of gene expression and protein localization is a crucial characteristic of life, where the spatiotemporal dynamic of secretory proteins instruct higher order organization, including the orchestration of developmental and adaptive programs. However tools to non-invasively interrogate the fate of secretory proteins in vivo are scarce. Here we introduce a genetic tagging strategy for in vivo imaging of the secretion and expression dynamics of secretory proteins in living animals. Applying this to a prototypical liver-derived secretory protein, we demonstrate that this approach, combined with optical in-vivo imaging, uncovers extrahepatic prothrombin expression in multiple novel anatomical sites (including testes, placenta, brain, kidney, heart and lymphatic system) and in emerging tumors, resulting in significant amounts of tumor-derived prothrombin in the blood with procoagulant properties. Syngeneic cell lines from this mouse model enable unravelling regulatory mechanisms in high resolution, and in a scalable format ex vivo. Beyond discovering new functions of proteins in a targeted manner, this model allows identifying rheostats in the cross-talk between gene expression and availability of a secretory protein. It is also a valuable resource for uncovering novel (tissue-specific) therapeutic vulnerabilities. Accuracy and Power Analysis of Social Networks Built From Count Data O_LIPower analysis is used to estimate the probability of correctly rejecting a null hypothesis for a given statistical model and dataset. Conventional power analyses assume complete information, but the stochastic nature of behavioural sampling can mean that true and estimated networks are poorly correlated. Power analyses do not currently take the effect of sampling into account. This could lead to inaccurate estimates of statistical power, potentially yielding misleading results.
C_LIO_LIHere we develop a method for computing network correlation: the correlation between an estimated social network and its true network, using a Gamma-Poisson model of social event rates for networks constructed from count data. We use simulations to assess how the level of network correlation affects the power of nodal regression analyses. We also develop a generic method of power analysis applicable to any statistical test, based on the concept of diminishing returns.
C_LIO_LIWe demonstrate that our network correlation estimator is both accurate and moderately robust to its assumptions being broken. We show that social differentiation, mean social event rate, and the harmonic mean of sampling times positively impacts the strength of network correlation. We also show that the required level of network correlation to achieve a given power level depends on many factors, but that 0.80 network correlation usually corresponds to around 0.80 power for nodal regression in ideal circumstances.
C_LIO_LIWe provide guidelines for using our network correlation estimator to verify the accuracy of networks built from count data, and to conduct power analysis. This can be used prior to data collection, in post hoc analyses, or even for subsetting networks in dynamic network analysis. The network correlation estimator and custom power analysis methods have been made available as an R package.
C_LI A vascularized tumoroid model for human glioblastoma angiogenesis Glioblastoma (GBM) angiogenesis is critical for tumor growth and recurrence, making it a compelling therapeutic target. Here, a disease-relevant, vascularized tumoroid in vitro model with stem-like features and stromal surrounds is reported. The model is used to recapitulate how individual components of the GBMs complex brain microenvironment such as hypoxia, vasculature-related stro-mal cells and growth factors support GBM angiogenesis. It is scalable, tractable, cost-effective and can be used with biologically-derived or biomimetic matrices. Patient-derived primary GBM cells are found to closely participate in blood vessel formation in contrast to a GBM cell line containing differentiated cells. Exogenous growth factors amplify this effect under normoxia but not at hypoxia suggesting that a significant amount of growth factors is already being produced under hypoxic conditions. Under hypoxia, primary GBM cells and umbilical vein endothelial cells are found to strongly co-localize in a mosaic pattern to form sprouting vascular networks, with GBM cells acquiring an endothelial-like behaviour, which has been reported to occur in vivo. These findings demonstrate that our 3D tumoroid in vitro model exhibits biomimetic attributes that may permit its use as a preclinical model in studying microenvironment cues of tumor angiogenesis. Virus-host interactions and genetic diversity of Antarctic sea ice bacteriophages Despite generally appreciated significant roles of microbes in sea ice and polar waters, detailed studies of virus-host systems from such environments have been so far limited by only a few isolates. Here, we investigated infectivity under various conditions, infection cycles, and genetic diversity of Antarctic sea ice bacteriophages: PANV1, PANV2, OANV1, and OANV2. The phages infect common sea ice bacteria belonging to the genera Paraglaciecola or Octadecabacter. Although the phages are marine and cold-active, replicating at 0-5 {degrees}C, they all survived temporal incubations at [&ge;]30 {degrees}C and remained infectious without any salts or supplemented only with magnesium, suggesting a robust virion assembly maintaining integrity under a wide range of conditions. Host recognition in the cold proved to be effective, and the release of progeny viruses occurred as a result of cell lysis. The analysis of viral genome sequences showed that nearly half of the gene products of each virus are unique, highlighting that sea ice harbors unexplored virus diversity. Based on predicted genes typical for tailed double-stranded DNA phages, we suggest placing the four studied viruses in the class Caudoviricetes. Searching against viral sequences from metagenomic assemblies revealed that related viruses are not restricted to Antarctica, but also found in distant marine environments.

ImportanceVery little is known about sea ice microbes despite the significant role of sea ice in the global oceans as well as microbial input into biogeochemical cycling. Studies on the sea ice viruses have been typically limited to -omics-based approaches and microscopic examinations of sea ice samples. Up to date, only four cultivable viruses have been isolated from Antarctic sea ice. Our study of these unique isolates advances the understanding of the genetic diversity of viruses in sea ice environments, their interactions with host microbes, and possible links to other biomes. Such information contributes to more accurate future sea ice biogeochemical models. A Eukaryotic Circuit for Secretion-coupled Cellular Autonomy Cancers represent complex autonomous robust systems, displaying interconnectivity with feedback control. Autonomous growth is fueled by a cancer cells ability to  secrete-and-sense growth factors: a poorly understood phenomenon. Using an integrated systems and experimental approach, here we dissect the impact of a feedback-coupled GTPase circuit within the secretory pathway that imparts secretion-coupled autonomy. The circuit is assembled when the Ras-superfamily monomeric GTPase Arf1, and the heterotrimeric GTPase Giabg and their corresponding GAPs and GEFs are coupled by GIV/Girdin, a protein that is known to fuel aggressive traits in diverse cancers. One forward and two key negative feedback loops within the circuit create closed-loop control (CLC), allow the two GTPases to coregulate each other, and convert the expected switch-like behavior of Arf1-dependent secretion into an unexpected dose response alignment behavior of sensing and secretion. Such behavior translates into cell survival that is self-sustained by stimuli-proportionate secretion. Findings highlight how enhanced coupling of two biological switches in cancer cells is critical for multiscale feedback control to achieve secretion-coupled autonomy.

GRAPHIC ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=163 HEIGHT=200 SRC="FIGDIR/small/436048v3_ufig1.gif" ALT="Figure 1">
View larger version (51K):
org.highwire.dtl.DTLVardef@183ff5corg.highwire.dtl.DTLVardef@f21de7org.highwire.dtl.DTLVardef@1893683org.highwire.dtl.DTLVardef@159d966_HPS_FORMAT_FIGEXP  M_FIG C_FIG In BriefThis work delivers an experimentally validated dynamical systems model of the cooperativity between two distinct classes of biological switches in eukaryotic cells and reveals the basis of secretion-coupled cellular autonomy in cancer.

HIGHLIGHTSO_LIModeling and experimental approaches were used to dissect a coupled GTPase circuit
C_LIO_LICoupling enables closed loop feedback and mutual control of GTPases
C_LIO_LICoupling generates dose response alignment behavior of sensing and secretion.
C_LIO_LICoupling is critical for multiscale feedback control to achieve secretion-coupled autonomy.
C_LI GABA facilitates spike propagation at branch points of sensory axons in the spinal cord Movement and posture depend on sensory feedback that is regulated by specialized GABAergic neurons (GAD2+) that form axo-axonic contacts onto myelinated proprioceptive sensory axons and are thought to be inhibitory. However, we report here that activating GAD2+ neurons, directly with optogenetics or indirectly by cutaneous stimulation, facilitates sensory feedback to motoneurons in awake rodents and humans. GABAA receptors and GAD2+ innervation at or near nodes of Ranvier of sensory axons cause this facilitation, preventing spike propagation failure at the many axon branch points, which is otherwise common without GABA. In contrast, GABAA receptors are generally lacking from axon terminals (unlike GABAB) and do not presynaptically inhibit transmitter release onto motoneurons. GABAergic innervation near nodes and branch points allows individual branches to function autonomously, with GAD2+ neurons regulating which branches conduct, adding a computational layer to the neuronal networks generating movement and likely generalizing to other CNS axons. IRGM/Irgm1 Aggravates Progression of Atherosclerosis by Inducing Macrophage Apoptosis through the MAPK Signaling Pathway AimsAtherosclerosis underlies most cardiovascular diseases, among which acute coronary syndrome (ACS) caused by plaque rupture (PR) often leads to death. Immune-related GTPases (IRGM/Irgm1) have been extensively studied in inflammatory diseases, but their role in atherosclerosis is unclear. Determining how IRGM/Irgm1 promotes atherosclerotic plaque vulnerability will provide information for new biomarkers and/or therapeutic targets.

Methods and resultsWe identified ruptured and unruptured plaques by optical coherence tomography, and found that serum IRGM was highly expressed in patients with ST-segment elevation myocardial infarction in PR. We used ApoE-/-Irgm1+/+, ApoE-/-Irgm1+/- mice and chimeric mice to establish a model of advanced atherosclerosis. The results of pathological experiments showed that Irgm1 caused plaque necrosis. The ratio of neutral lipids and cholesterol crystals increases, while the content of collagen fibers decreases, aggravating the destabilization of atherosclerotic plaques. In vitro, we used multiple approaches to confirm that Irgm1 promotes macrophage apoptosis by promoting the production of reactive oxygen species and activating the MAPK signaling pathway.

ConclusionsIRGM may be a potential risk factor for PR. Mechanistic studies have shown that IRGM/Irgm1 contributes to the formation and rupture of fragile plaques. This is partly mediated by the induction of macrophage apoptosis via the MAPK signaling pathway. IRGM may offer new strategies for early treatment of ACS.

Translation viewOur findings indicate that IRGM/Irgm1 contributes to formation and rupture of vulnerable plaques. It suggests that IRGM may provide a new target for the early treatment of ACS. Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: the Multi-Targeting Drug DREAM Challenge A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ( polypharmacology). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Author SummaryMany modern drugs are developed with the goal of modulating a single cellular pathway or target. However, many drugs are, in fact,  dirty; they target multiple cellular pathways or targets. This phenomenon is known as multi-targeting or polypharmacology. While some strive to develop  cleaner therapeutics that eliminate secondary targets, recent work has shown that multi-targeting therapeutics have key advantages for a variety of diseases. However, while multi-targeting drugs that affect a precisely-defined profile of targets may be more effective, it is difficult to computationally predict which molecules have desirable target profiles. Here, we report the results of a competitive crowdsourcing project (the Multi-Targeting Drug DREAM Challenge), where we challenged participants to predict chemicals that have desired target profiles for cancer and neurodegenerative disease. A chemical screen based on an interruption of zebrafish gastrulation identifies the HTR2C inhibitor Pizotifen as a suppressor of EMT-mediated metastasis Metastasis is responsible for approximately 90% of cancer-associated mortality but few models exist that allow for rapid and effective screening of anti-metastasis drugs. Current mouse models of metastasis are too expensive and time consuming to use for rapid and high-throughput screening. Therefore, we created a unique screening concept utilizing conserved mechanisms between zebrafish gastrulation and cancer metastasis for identification of potential anti-metastatic drugs. We hypothesized that small chemicals that interrupt zebrafish gastrulation might also suppress metastatic progression of cancer cells and developed a phenotype-based chemical screen to test the hypothesis. The screen used epiboly, the first morphogenetic movement in gastrulation, as a marker and enabled 100 chemicals to be tested in five hours. The screen tested 1280 FDA-approved drugs and identified Pizotifen, an antagonist for serotonin receptor 2C (HTR2C) as an epiboly-interrupting drug. Pharmacologic and genetic inhibition of HTR2C suppressed metastatic progression in a mouse model. Blocking HTR2C with Pizotifen restored epithelial properties to metastatic cells through inhibition of Wnt-signaling. In contrast, HTR2C induced epithelial to mesenchymal transition (EMT) through activation of Wnt-signaling and promoted metastatic dissemination of human cancer cells in a zebrafish xenotransplantation model. Taken together, our concept offers a novel platform for discovery of anti-metastasis drugs. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Ancient Darwinian replicators nested within eubacterial genomes Mobile genetic elements (MGEs), such as transposons and insertion sequences, propagate within bacterial genomes, but persistence times in individual lineages are short. For long-term survival, MGEs must continuously invade new hosts by horizontal transfer. Theoretically, MGEs that persist for millions of years in single lineages, and are thus subject to vertical inheritance, should not exist. Here we make the case for an exception -- a class of MGE termed REPIN. REPINs are non-autonomous MGEs whose duplication depends on non-jumping RAYT transposases. Comparisons of REPINs and typical MGEs show that replication rates of REPINs are orders of magnitude lower, REPIN population size fluctuations correlate with changes in available genome space, REPIN conservation depends on RAYT function, and REPIN diversity accumulates within host lineages. From these data it follows that REPINs persist for millions of years within single host lineages. Such long-term persistence is expected to generate conflicts arising from the diverging effects of selection acting simultaneously on REPINs and host genomes. Evidence of conflict comes from analyses of REPIN abundance and diversity in two distantly related bacterial species. Taken together, our data lead to the conclusion that REPINs are ancient Darwinian replicators that have evolved enduring, beneficial relationships, with eubacterial genomes. C. elegans PEZO-1 is a Mechanosensitive Ion Channel Involved in Food Sensation PIEZO channels are force sensors essential for physiological processes, including baroreception and proprioception. The Caenorhabditis elegans genome encodes an ortholog gene of the Piezo family, pezo-1, which is expressed in several tissues, including the pharynx. This myogenic pump is an essential component of the C. elegans alimentary canal, whose contraction and relaxation are modulated by mechanical stimulation elicited by food content. Whether pezo-1 encodes a mechanosensitive ion channel and contributes to pharyngeal function remain unknown. Here, we leverage genome editing, genetics, microfluidics, and electropharyngeogram recording to establish that pezo-1 is expressed in the pharynx, including in a proprioceptive-like neuron, and regulates pharyngeal function. Knockout (KO) and gain-of-function (GOF) mutants reveal that pezo-1 is involved in fine-tuning pharyngeal pumping frequency, sensing osmolarity, and food mechanical properties. Using pressure-clamp experiments in primary C. elegans embryo cultures, we determine that pezo-1 KO cells do not display mechanosensitive currents, whereas cells expressing wild-type or GOF PEZO-1 exhibit mechanosensitivity. Moreover, infecting the Spodoptera frugiperda cell line with a baculovirus containing the G-isoform of pezo-1 (among the longest isoforms) demonstrates that pezo-1 encodes a mechanosensitive channel. Our findings reveal that pezo-1 is a mechanosensitive ion channel that regulates food sensation in worms. Vectofusin-1 based T-cell transduction approach for generating murine chimeric antigen receptor T-cells for cancer. Gene transfer into human and murine T-cells using viral-based approaches has several promising therapeutic applications including the production of chimeric antigen receptor T-cell (CAR-T) therapy. The generation of murine CAR-T is paramount to test and validate immunocompetent mouse models for CAR-T therapy. Several viral transduction enhancers already exist for gene therapy with few limitations. In this study, we tested vectofusin-1, a short cationic peptide, as a soluble transduction enhancer for gammaretroviral transduction for the generation of anti-CD19 murine CAR-T. We found that in comparison to Retronectin, Vectofusin-1 is an equally optimal transduction enhancer for the generation of murine CAR-T cells. Development of a Novel PROTAC using the Nucleic Acid Aptamer as a Targeting Ligand for Tumor Selective Degradation of Nucleolin PROteolysis TArgeting Chimeras (PROTACs) induce targeted protein degradation by hijacking the intracellular ubiquitin proteasome system, thus emerging as a new strategy for drug development. However, most PROTACs generated lack cell-type selectivity and have poorly soluble in water. To address this drawback, we developed a novel PROTAC ZL216 using aptamer AS1411 as a targeting ligand of nucleolin to conjugate with a small molecule ligand of E3 ligase VHL, which shows high aqueous solubility and serum stability. Based on the differential expression of nucleolin on the cell surface, ZL216 could bind to and internalize into breast cancer cells, but not normal breast cells. Furthermore, we revealed that ZL216 promoted the formation of a nucleolin-ZL216-VHL ternary complex in breast cancer cells and potently induced nucleolin degradation in vitro and in vivo. As a result, ZL216 inhibited the proliferation and migration of breast cancer cells. These studies demonstrate that in addition to peptides and small molecule compounds, nuclei acid aptamers can also be used to generate PROTACs, which broadens the toolbox constructing PROTACs and provides a promising strategy for development of tumor-selective PROTACs. Expression and subcellular localization of USH1C/harmonin in the human retina provides insights into pathophysiology and therapy of Usher syndrome Usher syndrome (USH) is the most common form of hereditary deafness-blindness in humans. USH is a complex genetic disorder, assigned to three clinical subtypes differing in onset, course and severity, with USH1 being the most severe. Rodent USH1 models do not reflect the ocular phenotype observed in human patients to date; hence, little is known about the pathophysiology of USH1 in the human eye. One of the USH1 genes, USH1C, exhibits extensive alternative splicing and encodes numerous harmonin protein isoforms that function as scaffolds for organizing the USH interactome. RNA-seq analysis of human retina uncovered harmonin a1 as the most abundant transcript of USH1C. Bulk RNA-seq analysis and immunoblotting showed abundant expression of harmonin in Muller glia cells (MGCs) and retinal neurons. Furthermore, harmonin was localized in the terminal endfeet and apical microvilli of MGCs, presynaptic region (pedicle) of cones, and outer segments of rods as well as at adhesive junctions of MGCs and photoreceptors in the outer limiting membrane (OLM). Our data provide evidence for the interaction of harmonin with OLM molecules in MGCs and photoreceptors and rhodopsin in photoreceptors. Subcellular expression and colocalization of harmonin correlates with the clinical phenotype observed in USH1C patients. In addition, primary cilia defects in USH1C patient-derived fibroblasts could be reverted by delivery of harmonin a1 transcript isoform. Our data provide novel insights into photoreceptor cell biology, USH1C pathophysiology and for developing gene therapy treatment. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. An IPI based immune prognostic model for diffuse large B-cell lymphoma BackgroundInternational Prognostic Index (IPI) was widely used to better discriminate prognosis of patients with diffuse large B-cell lymphoma (DLBCL). However, there is a significant need to identify novel valuable biomarkers in the context of targeted therapies, such as immune checkpoint blockade (ICB) therapy.

MethodsGene expression data and clinical information of DLBCL were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets. 371 immune-related hub genes in DLBCL patients with different IPI levels were identified by weighted gene co-expression network analysis (WGCNA), and 8 genes were selected to construct an IPI-based immune prognostic model (IPI-IPM). Afterward, the genetic, somatic mutational and molecular profiles of IPI-IPM subgroups were analyzed, as well as the potential clinical response of ICB in different IPI-IPM subgroups.

ResultsThe IPI-IPM was constructed based on the expression of CMBL, TLCD3B, SYNDIG1, ESM1, EPHA3, HUNK, PTX3 and IL12A, where high-risk patients had shorter overall survival (OS) than low-risk patients, consistent with the results in the GEO cohorts. The comprehensive results showed that high IPI-IPM risk scores were correlated with immune-related signaling pathways, high KMT2D and CD79B mutation rates, as well as up-regulation of inhibitory immune checkpoints including PD-L1, BTLA and SIGLEC7, indicating more potential response to ICB therapy.

ConclusionThe IPI-IPM has independent prognostic significance for DLBCL patients, which provides an immunological perspective to elucidate the mechanisms on tumor progression, also sheds a light on developing immunotherapy for DLBCL. Analysis of single-cell RNA sequencing data basedon autoencoders BackgroundSingle-cell RNA sequencing (scRNA-Seq) experiments are gaining ground to study the molecular processes that drive normal development as well as the onset of different pathologies. Finding an effective and efficient low-dimensional representation of the data is one of the most important steps in the downstream analysis of scRNA-Seq data, as it could provide a better identification of known or putatively novel cell-types. Another step that still poses a challenge is the integration of different scRNA-Seq datasets. Though standard computational pipelines to gain knowledge from scRNA-Seq data exist, a further improvement could be achieved by means of machine learning approaches.

ResultsAutoencoders (AEs) have been effectively used to capture the non-linearities among gene interactions of scRNA-Seq data, so that the deployment of AE-based tools might represent the way forward in this context. We introduce here scAEspy, a unifying tool that embodies: (1) four of the most advanced AEs, (2) two novel AEs that we developed on purpose, (3) different loss functions. We show that scAEspy can be coupled with various batch-effect removal tools to integrate data by different scRNA-Seq platforms, in order to better identify the cell-types. We benchmarked scAEspy against the most used batch-effect removal tools, showing that our AE-based strategies outperform the existing solutions.

ConclusionsscAEspy is a user-friendly tool that enables using the most recent and promising AEs to analyse scRNA-Seq data by only setting up two user-defined parameters. Thanks to its modularity, scAEspy can be easily extended to accommodate new AEs to further improve the downstream analysis of scRNA-Seq data. Considering the relevant results we achieved, scAEspy can be considered as a starting point to build a more comprehensive toolkit designed to integrate multi single-cell omics. First detection of herpesvirus and prevalence of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Environmental control of Pub1 (NEDD4 family E3 ligase) in S. pombe is regulated byTORC2 and Gsk3. Cells respond to changing nutrient environments by adjusting the abundance of surface nutrient transporters and receptors. This can be achieved by modulating ubiquitin-dependent endocytosis, which in part is regulated by the NEDD4 family of E3 ligases. Here we report novel regulation of Pub1, a fission yeast Schizosaccharomyces pombe member of the NEDD4-family of E3 ligases. We show that nitrogen stress inhibits Pub1 function, thereby increasing the abundance of the amino acid transporter Aat1 at the plasma membrane and enhancing sensitivity to the toxic arginine analogue canavanine. We show that TOR complex 2 (TORC2) signalling negatively regulates Pub1, thus TORC2 mutants under nutrient stress have decreased Aat1 at the plasma membrane and are resistant to canavanine. Inhibition of TORC2 signalling increases Pub1 phosphorylation, and this is dependent on Gsk3 activity. Addition of the Tor inhibitor Torin1 increases phosphorylation of Pub1 at serine 199 (S199) by 2.5-fold, and Pub1 protein levels in S199A phospho-ablated mutants are reduced. S199 is conserved in NEDD4 and is located immediately upstream of a WW domain required for protein interaction. Together, we describe how the major TORC2 nutrient-sensing signalling network regulates environmental control of Pub1 to modulate the abundance of nutrient transporters. Intervene or wait? A model evaluating the timing of intervention in conservation conflicts adaptive management under uncertainty. The timing of biodiversity managers interventions can be critical to the success of conservation, especially in situations of conflict between conservation objectives and human livelihood, i.e., conservation conflicts. Given the uncertainty associated with complex social-ecological systems and the potentially irreversible consequences of delayed action for biodiversity and livelihoods, managers tend to simply intervene as soon as possible by precaution. However, refraining from intervening when the situation allows it can be beneficial, notably by saving critical management resources. Here, we introduce a strategy for managers to decide, based on monitoring, whether intervention is required or if waiting is possible. This study evaluates the performance of this waiting strategy compared to a strategy of unconditional intervention at every opportunity. We built an individual-based model of conservation conflict between a manager aiming to conserve an animal population and farmers aiming to maximize yield by protecting their crop from wildlife damage. We then simulated a budget-constrained adaptive management over time applying each strategy, while accounting for uncertainty around population dynamics and around decision-making of managers and farmers. Our results showed that when the decision for the manager to intervene was based on a prediction of population size trajectory, the waiting strategy performed at least as well as unconditional intervention while also allowing managers to save resources by avoiding unnecessary interventions. Under difficult budgetary constraints on managers, this waiting strategy ensured as high yields as unconditional intervention while significantly improving conservation outcomes by compensating managers lack of resources with the benefits accrued over waiting periods. This suggests that waiting strategies are worth considering in conservation conflicts, as they can facilitate equitable management with a more efficient use of management resources, which are often limiting in biodiversity conservation. Rapid evolution of seed dormancy during sunflower de-domestication Hybridization between crops and their wild relatives may promote the evolution of de-domesticated (feral) weeds. Wild sunflower (Helianthus annuus L.) is typically found in ruderal environments, but crop-wild hybridization may facilitate the evolution of weedy populations. Using one crop-specific mitochondrial marker (CMS-PET1) and 14 nuclear SSR markers, we studied the origin and genetic diversity of a recently discovered weedy population of sunflower (named BRW). Then, using a resurrection approach, we tested for rapid evolution of weedy traits (seed dormancy, herbicide resistance, and competitive ability) by sampling weedy and wild populations 10 years apart (2007 and 2017). All the weedy plants present the CMS-PET1 cytotype, confirming their feral origin. At the nuclear markers, BRW showed higher genetic diversity than the cultivated lines and low differentiation with one wild population, suggesting that wild hybridization increased their genetic diversity. We found support for rapid evolution towards higher seed dormancy, but not for higher competitive ability or herbicide resistance. Our results highlight the importance of seed dormancy during the earliest stages of adaptation and show that crop-wild hybrids can evolve quickly in agricultural environments. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. A Physiometric Investigation of Inflammatory Composites: Comparison of A Priori Aggregates, Empirically-identified Factors, and Individual Proteins Most research testing the association between inflammation and health outcomes (e.g., heart disease, diabetes, depression) has focused on individual proteins; however, some studies have used summed composites of inflammatory markers without first investigating dimensionality. Using two different samples (MIDUS-2: N = 1,255 adults, MIDUS-R: N = 863 adults), this study investigates the dimensionality of eight inflammatory proteins (C-reactive protein (CRP), interleukin (IL)-6, IL-8, IL-10, tumor necrosis factor- (TNF-), fibrinogen, E-selectin, and intercellular adhesion molecule (ICAM)-1) and compared the resulting factor structure to a) an "a priori"/tau-equivalent factor structure in which all inflammatory proteins equally load onto a single dimension (comparable to the summed composites) and b) proteins modeled individually (i.e., no latent variable) in terms of model fit, replicability, reliability, and their associations with health outcomes. An exploratory factor analysis indicated a two-factor structure (Factor 1: CRP and fibrinogen; Factor 2: IL-8 and IL-10) in MIDUS-2 and was replicated in MIDUS-R. Results did not clearly indicate whether the empirically-identified factor structure or the individual proteins modeled without a latent variable had superior model fit, but both strongly outperformed the "a priori"/tau-equivalent structure (which did not achieve acceptable model fit in any models). Modeling the empirically-identified factors and individual proteins (without a latent factor) as outcomes of medical diagnoses resulted in comparable conclusions. However, modeling individual proteins resulted in findings more robust to correction for multiple comparisons despite more conservative adjustments. Further, reliability for all latent variables was poor. These results indicate that modeling inflammation as a unidimensional construct equally associated with all available proteins does not fit the data well. Instead, individual inflammatory proteins or, potentially (if empirically supported and biologically-plausible) empirically-identified inflammatory factors should be used in accordance with theory.