Profiles of expressed mutations in single cells reveal subclonal expansion patterns and therapeutic impact of intratumor heterogeneity Advances in single-cell RNA sequencing (scRNAseq) technologies uncovered an unexpected complexity in tumors, underlining the relevance of intratumor heterogeneity to cancer progression and therapeutic resistance. Heterogeneity in the mutational composition of cancer cells is a result of distinct (sub)clonal expansions, each with a distinct metastatic potential and resistance to specific treatments. Unfortunately, due to their low read coverage per cell, scRNAseq datasets are too sparse and noisy to be used for detecting expressed mutations in single cells. Additionally, the large number of cells and mutations present in typical scRNAseq datasets are too large for available computational tools to, e.g., infer distinct subclones, lineages or trajectories in a tumor. Finally, there are no principled methods to assess distinct subclones inferred through single-cell sequencing data and the genomic alterations that seed and potentially cause them. Here we present Trisicell, a computational toolkit for scalable mutational intratumor heterogeneity inference and assessment from scRNAseq as well as single-cell genome or exome sequencing data. Trisicell allows reliable identification of distinct clonal lineages of a tumor, offering the ability to focus on the most important subclones and the genomic alterations that are associated with tumor proliferation. We comprehensively assessed Trisicell on a melanoma model by comparing distinct lineages and subclones it identifies on scRNAseq data, to those inferred using matching bulk whole exome (bWES) and transcriptome (bWTS) sequencing data from clonal sublines derived from single cells. Our results demonstrate that distinct lineages and subclones of a tumor can be reliably inferred and evaluated based on mutation calls from scRNAseq data through the use of Trisicell. Additionally, they reveal a strong correlation between aggressiveness and mutational composition, both across the inferred subclones, and among human melanomas. We also applied Trisicell to infer and evaluate distinct subclonal expansion patterns of the same mouse melanoma model after treatment with immune checkpoint blockade (ICB). After integratively analyzing our cell-specific mutation calls with their expression profiles, we observed that each subclone with a distinct set of novel somatic mutations is strongly associated with a specific developmental status. Moreover, each subclone had developed a unique ICB-resistance mechanism. These results demonstrate that Trisicell can robustly utilize scRNAseq data to delineate intratumor heterogeneity and help understand biological mechanisms underlying tumor progression and resistance to therapy.
Bayesian Additive Regression Trees for Genotype by Environment Interaction Models We propose a new class of models for the estimation of genotype by environment (GxE) interactions in plant-based genetics. Our approach, named AMBARTI, uses semi-parametric Bayesian additive regression trees to accurately capture marginal genotypic and environment effects along with their interaction in a fully Bayesian model. We demonstrate that our approach is competitive or superior to similar models widely used in the literature via both simulation and a real world data set. Furthermore, we introduce new types of visualisation to properly assess both the marginal and interactive predictions from the model. An R package that implements our approach is available at https://github.com/ebprado/ambarti.
Two Types of Cinnamoyl-CoA Reductase Function Divergently in Tissue Lignification, Phenylpropanoids Flux Control, and Inter-pathway Cross-talk with Glucosinolates as Revealed in Brassica napus Cinnamoyl-CoA reductase (CCR) is the entry point of lignin pathway and a crucial locus in dissection and manipulation of associated traits, but its functional dissection in Brassicaceae plants is largely lagged behind though Arabidopsis thaliana CCR1 has been characterized to certain extent. Here, 16 CCR genes are identified from Brassica napus and its parental species B. rapa and B. oleracea. Brassicaceae CCR genes are divided into CCR1 subfamily and CCR2 subfamily with divergent organ-specificity, yellow-seed trait participation and stresses responsiveness. CCR1 is preferential in G- and H-lignins biosynthesis and vascular development, while CCR2 has a deviation to S-lignin biosynthesis and interfascicular fiber development. CCR1 has stronger effects on lignification-related development, lodging resistance, phenylpropanoid flux control and seed coat pigmentation, whereas CCR2 controls sinapates levels. CCR1 upregulation could delay bolting and flowering time, while CCR2 upregulation weakens vascular system in leaf due to suppressed G lignin accumulation. Besides, CCR1 and CCR2 are deeply but almost oppositely linked with glucosinolates metabolism through inter-pathway crosstalk. Strangely, upregulation of both CCR1 and CCR2 could not enhance resistance to UV-B and S. sclerotiorum though CCR2 is sharply induced by them. These results provide systemic dissection on Brassica CCRs and CCR1-CCR2 divergence in Brassicaceae.

HighlightBrassicaceae contains two types of Cinnamoyl-CoA reductase. As revealed in Brassica napus, they are divergently involved in lignin monomer biosynthesis, tissue lignification, phenylpropanoid flux control, and inter-pathway crosstalk with glucosinolates. Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. Inhibition of microbial deconjugation of micellar bile acids protects against intestinal permeability and liver injury Altered host-microbe interactions and increased intestinal permeability have been implicated in the pathogenesis of a range of diseases. However, the mechanisms by which gut microbes affect epithelial barrier integrity remain unclear. Few host-produced metabolites that protect against epithelial damage have been identified, and whether microbial metabolism of host factors alters intestinal barrier function is unknown. Here, we investigate the effects of bacterial metabolism of host-produced bile acid (BA) metabolites on epithelial barrier integrity. We observe that rats fed a choline-deficient, high-fat diet (CDAHFD) exhibit reduced abundance of host-produced conjugated BAs in the intestine at early timepoints coinciding with increased permeability. We show that in vitro, conjugated BAs protect gut epithelial monolayers from damage caused by bacterially produced unconjugated BAs through micelle formation. We then demonstrate that inhibition of BA deconjugation using a small molecule inhibitor of gut bacterial bile salt hydrolase (BSH) enzymes prevents development of pathologic intestinal permeability and hepatic inflammation in CDAHFD-fed rats. Finally, we show that the predominant conjugated BAs in humans protect against epithelial barrier disruption in vitro. Our study identifies a protective role for conjugated BAs in intestinal epithelial barrier function and suggests that rational manipulation of microbial BA metabolism could be leveraged to regulate gut barrier integrity. Priority conservation areas and a global population estimate for the Critically Endangered Philippine Eagle derived from modelled range metrics using remote sensing habitat characteristics Many range-restricted taxa are currently experiencing population declines yet lack fundamental information regarding distribution and population size. Establishing baseline estimates for both these key biological parameters is however critical for directing conservation planning for at-risk range-restricted species. The International Union for the Conservation of Nature (IUCN) Red List uses three range metrics that define species distributions and inform extinction risk assessments: extent of occurrence (EOO), area of occupancy (AOO) and area of habitat (AOH). However, calculating all three metrics using standard IUCN approaches relies on a geographically representative sample of locations, which for rare species is often spatially biased. Here, we apply model-based interpolation using Species Distribution Models (SDMs), correlating occurrences with remote-sensing covariates, to calculate IUCN range metrics, protected area coverage and a global population estimate for the Critically Endangered Philippine Eagle (Pithecophaga jefferyi). Our final range wide continuous SDM had high predictive accuracy (Continuous Boyce Index = 0.927) and when converted to a binary model estimated an AOH = 23,185 km2, a maximum EOO = 605,759 km2, a minimum EOO = 272,272 km2, with an AOO = 53,867 km2. Based on inferred habitat from the AOH metric, we estimate a global population of 318 breeding pairs (range: 258-362 pairs), or 636 mature individuals, across the Philippine Eagle global range. Protected areas covered 34 % of AOH, 15 % less than the target representation, with the continuous model identifying key habitat as priority conservation areas. We demonstrate that even when occurrences are geographically biased, robust habitat models can be built that enable quantification of baseline IUCN range metrics, protected area coverage, and a population size estimate. In the absence of adequate location data for many rare and threatened taxa, our method is a promising spatial modelling tool with widespread applications, in particular for island endemics facing high extinction risk. Re-convolving the compositional landscape of primary and recurrent glioblastoma reveals prognostic andtargetable tissue states Glioblastoma is an aggressive diffusely infiltrating neoplasm that spreads beyond surgical resection margins, where it intermingles with non-neoplastic brain cells. This complex microenvironment harbouring infiltrating glioma and non-neoplastic brain cells is the origin of tumor recurrence. Thus, understanding the cellular and molecular features of the glioma microenvironment is therapeutically and prognostically important. Here, we used single-nucleus RNA sequencing (snRNAseq) of primary and recurrent glioma to define three compositional  tissue-states rooted in cohabitation of cell-types and transcriptional states. These comprise a state featuring A) abundant normal brain microenvironment cells, B) reactive/inflamed brain microenvironment, and C) cellular/proliferative tumor. All these states exhibited variable degrees of infiltration by glioma cells. Spatial transcriptomics confirmed that the cell-types and transcriptomics states which compositionally cohabitate tissue states also colocalize in space. Tissue states are clinically significant because they correlate with radiographic, histopathologic, and prognostic features. We found that enrichment of tissue state B signature correlated with shorter survival. Importantly, we found that our compositionally defined tissue states are enriched in distinct metabolic pathways. One such pathway is fatty acid biosynthesis, which was enriched in tissue state B - a state enriched in recurrent glioblastoma and composed of astrocyte-like/mesenchymal glioma cells, reactive astrocytes resembling those seen in neurodegeneration, and monocyte-like myeloid cells. We showed that treating acute slices of GBM with a fatty acid synthesis inhibitor is sufficient to deplete the transcriptional signature of tissue state B. Our findings define a novel compositional approach to the glioma-infiltrated tissue which allows us to discover prognostic and targetable features, paving the way to new mechanistic and therapeutic discoveries. Chemogenetically decreasing activity of the dorsal cochlear nucleus can ameliorate noise-induced tinnitus in mice The dorsal cochlear nucleus (DCN) is a region known to integrate somatosensory and auditory inputs and is identified as a potential key structure in the generation of phantom sound perception, especially noise-induced tinnitus. Yet, how altered homeostatic plasticity of the DCN induces and maintains the sensation of tinnitus is not clear. Here, we chemogenetically decrease activity of a subgroup of DCN neurons, Ca2+/Calmodulin kinase 2 (CaMKII) positive DCN neurons, using Gi-coupled human M4 Designer Receptors Exclusively Activated by Designer Drugs (hM4Di DREADDs), to investigate their role in noise-induced tinnitus. Mice were exposed to loud noise (9-11kHz, 90dBSPL, 1h, followed by 2h of silence) and auditory brainstem responses (ABRs) and gap prepulse inhibition of acoustic startle (GPIAS) were recorded two days before and two weeks after noise exposure to identify animals with a significantly decreased inhibition of startle, indicating tinnitus but without permanent hearing loss. Neuronal activity of CaMKII+ neurons expressing hM4Di in the DCN was lowered by administration of clozapine-N-oxide (CNO). We found that acutely decreasing firing rate of CaMKII+ DCN units decrease tinnitus-like responses (p = 0.038, n = 11 mice), compared to the control group that showed no improvement in GPIAS (control virus; CaMKII-YFP + CNO, p = 0.696, n = 7 mice). Extracellular recordings confirmed CNO to decrease unit firing frequency of CaMKII-hM4Di+ mice and alter best frequency and tuning width of response to sound. However, these effects were not seen if CNO had been previously administered during the noise exposure (n = 6 experimental and 6 control mice). Our results suggest that CaMKII-hM4Di positive cells in the DCN are not crucial for tinnitus induction but play a significant role in maintaining tinnitus perception in mice. Anaplerotic flux into the Calvin-Benson cycle. Hydrogen isotope evidence for in vivo occurrence in C3 metabolism - As the central carbon uptake pathway in photosynthetic cells, the Calvin-Benson cycle is among the most important biochemical cycles for life on Earth. A carbon flux of anaplerotic origin (i.e., through the chloroplast-localised oxidative branch of the pentose phosphate pathway) into the Calvin-Benson cycle was proposed recently.
- Here, we measured intramolecular deuterium abundances in leaf starch of Helianthus annuus grown at varying ambient CO2 concentrations, Ca. Additionally, we modelled deuterium fractionations expected for the anaplerotic pathway and compared modelled with measured fractionations.
- We report deuterium fractionation signals at H1 and H2 of starch glucose. Below a Ca change point, these signals increase with decreasing Ca consistent with modelled fractionations by anaplerotic flux. Under standard conditions (Ca=450 ppm corresponding to intercellular CO2 concentrations, Ci, of 328 ppm), we estimate negligible anaplerotic flux. At Ca=180 ppm (Ci=140 ppm), more than 10% of the glucose 6-phosphate entering the starch biosynthesis pathway is diverted into the anaplerotic pathway.
- In conclusion, we report evidence consistent with anaplerotic carbon flux into the Calvin-Benson cycle in vivo. We propose the flux may help to (i) maintain high levels of ribulose 1,5-bisphosphate under source-limited growth conditions to facilitate photorespiratory nitrogen assimilation required to build-up source strength and (ii) counteract oxidative stress. The genome sequence of the avian vampire fly (Philornis downsi), an invasive nest parasite of Darwin's finches in Galapagos The invasive avian vampire fly (Philornis downsi) is considered one of the greatest threats to the unique and endemic avifauna of the Galapagos Islands, Ecuador. The fly parasitizes nearly every passerine species, including Darwins finches, in the Galapagos. The fly is thought to have been introduced from mainland Ecuador, although the full pathway of invasion is not yet known. The majority of research to date has focused on the effects of the fly on the fitness of avian hosts and explorations of mitigation methods. A lag in research related to the genetics of this invasion demonstrates, in part, a need to develop full-scale genomic resources with which to address further questions within this system. In this study, an adult P. downsi collected from San Cristobal Island within the Galapagos archipelago was sequenced to generate a high-quality genome assembly. We examined various features of the genome (e.g., coding regions, non-coding transposable elements) and carried out comparative genomics analysis against other dipteran genomes. We identified lists of gene families that are significantly expanding/contracting in P. downsi that are related to insecticide resistance, detoxification, and potential feeding ecology and counter defense against host immune responses. The P. downsi genome assembly provides an important foundational resource for studying the molecular basis of its successful invasion in the Galapagos and the dynamics of its population across multiple islands. The findings of significantly changing gene families associated with insecticide resistance and immune responses highlight the need for further investigations into the role of different gene families in aiding the flys successful invasion. Furthermore, this genomic resource will also better help inform future research studies and mitigation strategies aimed at minimizing the flys impact on the birds of the Galapagos. The Consequences of Differential Origin Licensing Dynamics in Distinct Chromatin Environments MCM complexes are loaded onto chromosomes to license DNA replication origins in G1 phase of the cell cycle, but it is not yet known how mammalian MCM complexes are adequately distributed to both euchromatin and heterochromatin. To address this question, we combined time-lapse live-cell imaging with fixed cell immunofluorescence imaging of single human cells to quantify the relative rates of MCM loading in heterochromatin and euchromatin at different times within G1. We report here that MCM loading in euchromatin is faster than in heterochromatin in very early G1, but surprisingly, heterochromatin loading accelerates relative to euchromatin loading in middle and late G1. These different loading dynamics require ORCA-dependent differences in ORC distribution during G1. A consequence of heterochromatin origin licensing dynamics is that cells experiencing a truncated G1 phase from premature cyclin E expression enter S phase with under-licensed heterochromatin, and DNA damage accumulates preferentially in heterochromatin in the subsequent S/G2 phase. Thus G1 length is critical for sufficient MCM loading, particularly in heterochromatin, to ensure complete genome duplication and to maintain genome stability. Multi-omics investigation of tacrolimus nephrotoxicity BackgroundTacrolimus, an immunosuppressive drug prescribed to a majority of transplanted patients is nephrotoxic, through still unclear mechanisms. This study aims to evaluate the impact of tacrolimus on a lineage of proximal tubular cells using a multi-omics approach.

MethodsLLC-PK1 cells were exposed to 5 M of tacrolimus for 24h. Intracellular proteins and metabolites, and extracellular metabolites were extracted and analysed by LC-MS/MS. The transcriptional expression of the dysregulated proteins PCK-1, FBP1 and FBP2 was measured using RT-qPCR.

ResultsIn our cell model, tacrolimus impacted different metabolic pathways including those of arginine (e.g., citrulline, ornithine) (p < 0.0001), amino acids (e.g., valine, isoleucine, aspartic acid) (p < 0.0001) and pyrimidine (p<0.01). In addition, it induced oxidative stress (p < 0.01) as shown by a decrease in total cell glutathione quantity. It impacted cell energy through an increase in Krebs cycle intermediates (e.g., citrate, aconitate, fumarate) (p < 0.01) and down-regulation of PCK-1 (p < 0.05) and FPB1 (p < 0.01), which are key enzymes in gluconeogenesis. Apart from glucose synthesis, gluconeogenesis is an important process in kidney mediated acid-base balance control.

ConclusionThe variations found using this multi-omics approach clearly point towards a dysregulation of energy production in epithelial cells of the renal tubule, and potentially of their functions, that may be implicated in tacrolimus nephrotoxicity in the clinics. Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less-functional Misfolded protein conformations with decreased functionality can bypass the proteostasis machinery and remain soluble in vivo. This is an unexpected phenomenon as several cellular quality control mechanisms have evolved to rid cells of misfolded proteins. Three questions, then, are: how is it structurally possible for long-lived, soluble, misfolded proteins to bypass the proteostasis machinery and processes? How widespread are these soluble, misfolded states across the proteome? And how long do they persist for? Here, we address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of all proteins exhibit subpopulations of misfolded conformations that are likely to bypass molecular chaperones, avoid aggregation, and not be rapidly degraded. These misfolded states may persist for months or longer for some proteins. Structurally characterizing these misfolded states, we observe they have a large amount of native structure, but also contain localized misfolded regions from non-native changes in entanglement, in which a protein segment threads through a loop formed by another portion of the protein that is not found in the native state. The surface properties of these misfolded states are native like, suggesting they may bypass the proteostasis machinery and its regulatory processes to remain soluble, while their entanglements make these states long-lived kinetic traps, as disentanglement requires unfolding of already folded portions of the protein. In terms of function, we predict that one-third of proteins have subpopulations that misfold into less-functional states that have structurally perturbed functional sites yet remain soluble. Data from limited-proteolysis mass spectrometry experiments, which interrogate the misfolded conformations populated by proteins upon unfolding and refolding, are consistent with the structural changes seen in the entangled states of glycerol-3-phosphate dehydrogenase upon misfolding. These results provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass cellular quality controls, and indicate, unexpectedly, this may be a wide-spread phenomenon in proteomes. Such entanglements are observed in many native structures, suggesting the non-native entanglements we observe are plausible. More broadly, these near-native entangled structures suggest a hypothesis for how synonymous mutations can modulate downstream protein structure and function, with these mutations partitioning nascent proteins between these kinetically trapped states. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. Clp protease and antisense RNA jointly regulate the global regulator CarD to mediate mycobacterial starvation response Under starvation conditions, bacteria tend to slow down their translation rate by reducing rRNA synthesis, but the way they accomplish that may vary in different bacteria. In Mycobacterium species, transcription of rRNA is activated by the RNA polymerase (RNAP) accessory transcription factor CarD, which interacts directly with RNAP to stabilize the RNAP-promoter open complex formed on rRNA genes. The functions of CarD have been extensively studied, but the mechanisms that control its expression remain obscure. Here, we report that the level of CarD was tightly regulated when mycobacterial cells switched from nutrient-rich to nutrient-deprived conditions. At the translational level, an antisense RNA of carD (AscarD) was induced in a SigF-dependent manner to bind with carD mRNA and inhibit CarD translation, while at the post-translational level, the residual intracellular CarD was quickly degraded by the Clp protease. AscarD thus worked synergistically with Clp protease to decrease the CarD level to help mycobacterial cells cope with the nutritional stress. Altogether, our work elucidates the regulation mode of CarD and delineates a new mechanism for the mycobacterial starvation response, which is important for the adaptation and persistence of mycobacterial pathogens in the host environment. 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. Genomically Mined Acoustic Reporter Genes Enable In Vivo Monitoring of Tumors and Tumor-Homing Bacteria A major outstanding challenge in the fields of biological research, synthetic biology and cell-based medicine is the difficulty of visualizing the function of natural and engineered cells noninvasively inside opaque organisms. Ultrasound imaging has the potential to address this challenge as a widely available technique with a tissue penetration of several centimeters and spatial resolution below 100 {micro}m. Recently, the first genetically encoded acoustic reporters were developed based on bacterial gas vesicles to link ultrasound signals to molecular and cellular function. However, the properties of these first-generation acoustic reporter genes (ARGs) resulted in limited sensitivity and specificity for imaging in the in vivo context. Here, we describe second-generation ARGs with greatly improved acoustic properties and expression characteristics, identified through a phylogenetic screen of candidate gene clusters from diverse bacteria and archaea. The resulting constructs offer major qualitative and quantitative improvements, including much stronger ultrasound contrast, the ability to produce nonlinear signals distinguishable from background tissue, and stable long-term expression. We demonstrate the utility of these next-generation ARGs by imaging in situ gene expression in mouse models of breast cancer and tumor-homing probiotic bacteria, revealing the unique spatial distributions of tumor growth and colonization by therapeutic cells noninvasively in living subjects. Gradients in cellular metabolism regulate development of tonotopy along the cochlea. In vertebrates with elongated auditory organs, mechanosensory hair cells (HCs) are organised such that complex sounds are broken down into their component frequencies along the basal-to-apical long (tonotopic) axis. To generate the frequency-specific characteristics required at the appropriate positions, nascent HCs must determine their tonotopic properties during development. This relies on complex signalling between the developing HC and its local niche along the axis within the auditory organ. Here we apply NAD(P)H fluorescence lifetime imaging (FLIM) and live imaging of mitochondria to reveal metabolic gradients along the tonotopic axis of the developing cochlea. We further show that re-shaping these gradients during development, by inhibiting cytosolic glucose metabolism, alters normal Bmp7 and Chordin like-1 signalling leading to flattening of tonotopic patterning along the axis. Our work supports a causal link between morphogen signalling and metabolic reprogramming in specifying tonotopic morphologies in auditory HCs. Contribution of arm movements to balance recovery after tripping in older adults. Falls are common in daily life, often caused by trips and slips and, particularly in older adults, with serious consequences. Although arm movements play an important role in balance control, there is limited research into the role of arm movements during balance recovery after tripping in older adults. We investigated how older adults use their arms to recover from a trip and the difference in the effects of arm movements between fallers (n=5) and non-fallers (n=11).

Sixteen older males and females (69.7{+/-}2.3 years) walked along a walkway and were occasionally tripped over suddenly appearing obstacles. We analysed the first trip using a biomechanical model based on full-body kinematics and force-plate data to calculate whole body orientation during the trip and recovery phase. With this model, we simulated the effects of arm movements at foot-obstacle impact and during trip recovery on body orientation.

Apart from an increase in sagittal plane forward body rotation at touchdown in fallers, we found no significant differences between fallers and non-fallers in the effects of arm movements on trip recovery. Like earlier studies in young adults, we found that arm movements during the recovery phase had most favourable effects in the transverse plane: by delaying the transfer of angular momentum of the arms to the body, older adults rotated the tripped side more forward thereby allowing for a larger recovery step. Older adults that are prone to falling might improve their balance recovery after tripping by learning to prolong ongoing arm movements. High performance single-cell gene regulatory network inference atscale: The Inferelator 3.0 MotivationGene regulatory networks define regulatory relationships between transcription factors and target genes within a biological system, and reconstructing them is essential for understanding cellular growth and function. Methods for inferring and reconstructing networks from genomics data have evolved rapidly over the last decade in response to advances in sequencing technology and machine learning. The scale of data collection has increased dramatically; the largest genome-wide gene expression datasets have grown from thousands of measurements to millions of single cells, and new technologies are on the horizon to increase to tens of millions of cells and above.

ResultsIn this work, we present the Inferelator 3.0, which has been significantly updated to integrate data from distinct cell types to learn context-specific regulatory networks and aggregate them into a shared regulatory network, while retaining the functionality of the previous versions. The Inferelator is able to integrate the largest single-cell datasets and learn cell-type specific gene regulatory networks. Compared to other network inference methods, the Inferelator learns new and informative Saccharomyces cerevisiae networks from single-cell gene expression data, measured by recovery of a known gold standard. We demonstrate its scaling capabilities by learning networks for multiple distinct neuronal and glial cell types in the developing Mus musculus brain at E18 from a large (1.3 million) single-cell gene expression dataset with paired single-cell chromatin accessibility data.

AvailabilityThe inferelator software is available on GitHub (https://github.com/flatironinstitute/inferelator) under the MIT license and has been released as python packages with associated documentation (https://inferelator.readthedocs.io/). Targeting CISH enhances natural cytotoxicity receptor signaling and reduces NK cell exhaustion to improve solid tumor immunity Cytokine inducible SH2-containing protein (CISH) is a natural killer (NK) cell negative regulator of cytokine signaling pathway. To further understand CISH functions in NK cells, we developed a conditional Cish-deficient mouse model in NK cells (Cishfl/flNcr1Ki/+). We detected no developmental or homeostatic difference in NK cells.

However, global gene expression of Cishfl/flNcr1Ki/+ NK cells compared to Cish+/+Ncr1Ki/+ NK cells revealed upregulation of pathways and genes associated with NK cell cycling and activation. We show that CISH does not only regulate interleukin-15 (IL-15) signaling pathways but also natural cytotoxicity receptors (NCR) pathways. Indeed, CISH protein expression level increases upon NCR triggering. Primed Cishfl/flNcr1Ki/+ NK cells display increased activation upon NCR stimulation. Cishfl/flNcr1Ki/+ NK cells display lower activation thresholds and Cishfl/flNcr1Ki/+ mice are more resistant to tumor metastasis. Remarkably, we found that Cishfl/flNcr1Ki/+ mice were also more resistant to primary breast cancer growth in addition to superior control of spontaneous tumor metastasis. CISH deletion favors NK cell accumulation to the primary tumor, optimizes NK cell killing properties and decreases TIGIT immune checkpoint receptor expression, limiting NK cell exhaustion. Finally, we argue that specifically enhancing NK cell function is sufficient to boost anti-tumor response to both primary and secondary tumor models. Using CRISPRi, we then targeted CISH in human NK-92 or primary NK cells. According to the results in our mouse model, CISH deletion favors NCR signaling and anti-tumor functions in human NK cells. Our results validate CISH as an emerging therapeutic target to enhance NK cell immunotherapy. First Report of Lactoria cornuta from Karnataka, West coast of India, Role of tropical cyclones in Range Extension India is known for its rich biodiversity and is fortunate to have several endemic species from different classes of vertebrates. India is home to 7.5% of the global fish diversity, with 91 endemic species of ray finned fishes, the actinopterygians. Some fish species have never been reported until recently, one such example is the long-horned Cowfish (Ostracidae), best known for inhabiting only certain areas of the Indo-Pacific oceans. However, this fish has been reported from Bay of Bengal, Costs of Goa and Lakshadweep. This is the first time this fish has been found on the Coasts of Karnataka, west coast of India. Given its morphology, the migration seems highly improbable since it can only endure lethargic swimming. In this article, we have explored the role of other external forces that could have contributed to its range extension in the east and west coast of India. In this paper we are exploring the role of Super Cyclones as a vital force in determining the expansion of its range. The fact of collection longhorn cowfish for the first time from the coast of Karwar coincides with the passage of Cyclone Ockhi on the west coast of India this suggests the possible role of super cyclone in its range extension. A year later, another specimen was recovered from a location very close to the first, which indicates the successful establishment of Lacturia cornuta in its new environment. Related morphometric and meristic studies of our specimens are consistent with previous studies reported in the Bay of Bengal. 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. High-throughput variant detection using a color-mixing strategy Many diseases are related to multiple genetic alterations along a single gene. Probing for highly multiple (>10) variants in a single qPCR tube is not possible due to a limited number of fluorescence channels and one variant per channel, so many more tubes are needed. Here, we experimentally validate our novel color-mixing strategy that uses fluorescence combinations as digital color codes to probe multiple variants simultaneously. The color-mixing strategy relies on a simple intra-tube assay that can probe for 15 variants as part of an inter-tube assay that can probe for an exponentially increased number of variants. The color-mixing strategy is achieved using multiplex double-stranded toehold probes modified with fluorophores and quenchers; the probes are designed to be quenched or luminous after binding to wildtype or variant templates. We used the color-mixing strategy to probe for 21 pathogenic mutations in thalassemia and to distinguish between heterozygous and homozygous variants in 6 tubes, with a specificity of 99% and a sensitivity of 94%. To support tuberculosis diagnosis, we used the same strategy to simultaneously probe in Mycobacterium tuberculosis for rifampicin-resistance mutations occurring within one 81-bp region and one 48-bp region in rpoB gene, plus five isoniazid-resistance mutations in inhA and katG genes. The Evolution of Protein Folds by Creative Destruction Mechanisms by which new protein folds emerge and diverge pose central questions in biological sciences. Incremental mutation and step-wise adaptation explain relationships between topologically similar protein folds. However, the universe of folds is diverse and riotous, suggesting roles of more potent and creative forces. Sequence and structure similarity are observed between topologically distinct folds, indicating that proteins with distinct folds may share common ancestry.

We found evidence of common ancestry between three distinct {beta}-barrel folds: OB, SH3 and cradle loop barrel (CLB). The data suggest a mechanism of fold evolution that interconverts SH3, OB and CLB. This mechanism, which we call creative destruction, can be generalized to explain other examples of fold evolution including circular permutation. In creative destruction, an open reading frame duplicates or otherwise merges with another to produce a fused polypeptide. A merger forces two ancestral domains into a new sequence and spatial context. The fused polypeptide can explore folding landscapes that are inaccessible to either of the independent ancestral domains. However, the folding landscapes of the fused polypeptide are not fully independent of those of the ancestral domains. Creative destruction is thus partially conservative in that a daughter fold would inherit some motifs from the ancestral folds. After a merger and refolding, adaptive processes such as mutation and loss of extraneous segments optimize the new daughter fold.

SignificanceMechanisms of emergence and early diversification of structured proteins present deep and difficult problems in evolutionary biology. Here we excavate the deepest evolutionary history, found within the translation machinery, which is an ancient molecular fossil and the birthplace of all proteins. We demonstrate common origins of some of the simplest, oldest and most common protein folds. Furthermore, the data suggest a mechanism, that we call creative destruction, that explains at molecular level how simple folds spawn more complex folds. In this mechanism, new folds emerge from old folds via gene duplication, expression, exploration of new folding landscapes and adaptation. Creative destruction explains the facile emergence of complex from simple architectures in a funneled exploration. Innovative GenExpA software for selecting suitable reference genes for reliable normalization of gene expression in melanoma The algorithms commonly used to select the best stable reference gene in RT-qPCR data analysis have their limitations. We showed that simple selection of the reference gene or pair of genes with the lowest stability value from the pool of potential reference genes - a commonly used approach - is not sufficient to accurately and reliably normalize the target gene transcript and can lead to biologically incorrect conclusions. For reliable assessment of changes in a target gene expression level, we propose our innovative GenExpA software, which works in a manner independent of the experimental model and the normalizer used. GenExpA software selects the best reference by combining the NormFinder algorithm with progressive removal of the least stable gene from the candidate genes in a given experimental model and in the set of daughter models assigned to it. The reliability of references is validated based on the consistency of the statistical analyses of normalized target gene expression levels through all models, described by the coherence score (CS). The use of the CS value imparts a new quality to qPCR analysis because it clarifies how low the stability value of reference must be in order for biologically correct conclusions to be drawn. We tested our method on qPCR data for the B4GALT genes family in melanoma, which is characterized by a high mutation rate, and in melanocytes. GenExpA is available at https://github.com/DorotaHojaLukowicz/GenExpA or https://www.sciencemarket.pl/baza-programow-open-source#oferty.

HighlightsGenExpA - next-generation software for normalizer selection and validation

The GenExpA tool defines how low the stability value of the reference should be

The GenExpA tool determines the level of gene expression analysis robustness

The GenExpA tool increases the speed of analysis and reduces cost