Cell type-specific differential expression in spatial transcriptomics Spatial transcriptomics enables spatially resolved gene expression measurements at near single-cell resolution. There is a pressing need for computational tools to enable the detection of genes that are differentially expressed (DE) within specific cell types across tissue context. We show that current approaches cannot learn cell type-specific DE due to changes in cell type composition across space and the fact that measurement units often detect transcripts from more than one cell type. Here, we introduce a statistical method, Cell type-Specific Inference of Differential Expression (C-SIDE), that identifies cell type-specific patterns of differential gene expression while accounting for localization of other cell types. We model spatial transcriptomics gene expression as an additive mixture across cell types of general log-linear cell type-specific expression functions. This approach provides a unified framework for defining and identifying gene expression changes in a wide-range of relevant contexts: changes due to pathology, anatomical regions, physical proximity to specific cell types, and cellular microenvironment. Furthermore, our approach enables statistical inference across multiple samples and replicates when such data is available. We demonstrate, through simulations and validation experiments on Slide-seq and MER-FISH datasets, that our approach accurately identifies cell type-specific differential gene expression and provides valid uncertainty quantification. Lastly, we apply our method to characterize spatially-localized tissue changes in the context of disease. In an Alzheimers mouse model Slide-seq dataset, we identify plaque-dependent patterns of cellular immune activity. We also find a putative interaction between tumor cells and myeloid immune cells in a Slide-seq tumor dataset. We make our C-SIDE method publicly available as part of the open source R package https://github.com/dmcable/spacexr.
Expansion of mutation-driven haematopoietic clones is associated with insulin resistance and low HDL-cholesterol in individuals with obesity AimsHaematopoietic clones caused by somatic mutations with [&ge;]2% variant allele frequency (VAF), known as clonal haematopoiesis of indeterminate potential (CHIP), increase with age and have been linked to risk of haematological malignancies and cardiovascular disease. Recent observations suggest that smaller clones are also associated with adverse clinical outcomes. Our aims were to determine the prevalence of clonal haematopoiesis driven by clones of variable sizes, and to examine the development of clones over time in relation to age and metabolic dysregulation over up to 20 years in individuals with obesity.

Methods and ResultsWe used an ultrasensitive single-molecule molecular inversion probe sequencing assay to identify clonal haematopoiesis driver mutations (CHDMs) in blood samples from individuals with obesity from the Swedish Obese Subjects study. In a single-timepoint dataset with samples from 1050 individuals, we identified 273 candidate CHDMs in 216 individuals, with VAF ranging from 0.01% to 31.15% and CHDM prevalence and clone sizes increasing with age. Longitudinal analysis over 20 years in CHDM-positive samples from 40 individuals showed that small clones can grow over time and become CHIP. VAF increased on average by 7% (range -4% to 27%) per year. Rate of clone growth was positively associated with insulin resistance (R=0.40, P=0.025) and low circulating levels of high-density lipoprotein-cholesterol (HDL-C) (R=-0.68, P=1.74E-05).

ConclusionOur results show that haematopoietic clones can be detected and monitored before they become CHIP and indicate that insulin resistance and low HDL-C, well-established cardiovascular risk factors, are associated with clonal expansion in individuals with obesity.

Translational perspectivesClonal haematopoiesis-driver mutations are somatic mutations in haematopoietic stem cells that lead to clones detectable in peripheral blood. Haematopoietic clones with a variant allele frequency (VAF) [&ge;]2%, known as clonal haematopoiesis of indeterminate potential (CHIP), are recognized as an independent cardiovascular risk factor. Here, we show that smaller clones are prevalent, and also correlate with age. Our longitudinal observations in individuals with obesity over 20 years showed that more than half of all clone-positive individuals show growing clones and clones with VAF <2% can grow and become CHIP. Importantly, clone growth was accelerated in individuals with insulin resistance and low high-density lipoprotein-cholesterol (HDL-C).

Translational outlook 1: Haematopoietic clones can be detected and monitored before they become CHIP.

Translational outlook 2: The association between insulin resistance and low HDL-C with growth of haematopoietic clones opens the possibility that treatments improving metabolism, such as weight loss, may reduce growth of clones and thereby cardiovascular risk.

One Sentence SummaryIn obesity, the growth rate of mutation-driven haematopoietic clones increased with insulin resistance and low HDL-C, both known risk factors for cardiovascular disease.
Micro topographical instruction of bacterial attachment, biofilm formation and in vivo host response Bio-instructive materials that prevent bacterial biofilm formation and drive an appropriate host immune response have the potential to significantly reduce the burden of medical device-associated infections. As bacterial surface attachment is sensitive to surface topography, we experimentally surveyed 2,176 combinatorially generated shapes using an unbiased high-throughput micro topographical polystyrene polymer chip screen. This identified topographies that reduced colonization in vitro by up to 15-fold compared with a flat surface for both motile and non-motile bacterial pathogens. Equivalent reductions were achieved on polyurethane, a polymer commonly used in medical devices. Using machine learning methods, a set of design rules based on generalisable topo-descriptors was established for predicting bacteria-resistant micro topographies. In a murine foreign body infection model, anti-attachment topographies were shown to be refractory to P. aeruginosa colonization and serendipitously, the fibrotic response to the implant was reduced, offering control of foreign body response using simple topographical patterning of non-eluting implants. Comprehensive-immunopeptidomics analysis reveals presentation of potential neoantigens carrying cancer driver mutations Understanding the properties of human leukocyte antigen (HLA) peptides (HLAp) is essential for cancer precision medicine, while the direct identification of HLAp from a tiny piece of clinical tissues by mass spectrometry (MS) is still confronted with technical challenges. Here, to overcome that hindrances, high-field asymmetric waveform ion mobility spectrometry (FAIMS) is introduced to conduct differential ion mobility (DIM)-MS by seamless gas-phase fractionation optimal for scarce samples. The global-immunopeptidomics analysis enables in-depth identification of 44,785 HLAp, including 2 neoantigens with KRAS-G12V and CPPED1-R228Q from stage IV colorectal cancer (CRC) tissues obtained from 17 subjects. Comparison of tissue-based personal immunopeptidome disclosed cancer-specific processing of HLAp. Since the direct identification of neoantigens from CRC tissues suggests more potential neoantigens have yet to been identified, the targeted-immunopeptidomics by screening-oriented parallel reaction monitoring (PRM) has been established and identified two more neoantigens with oncogenic KRAS. The targeted-immunopeptidomics analysis clarifies the presentation of aimed neoantigens, including neoantigens that are unpredictable by computational algorithms. The comprehensive-immunopeptidomics approach combining the FAIMS-assisted DIM-MS and the screening-oriented PRM contributes the direct and effective identification of HLAp and emerging neoantigens from clinical tissue samples, as an advanced strategy to identify the authentic and recurrent neoantigens for cancer immunotherapy. Applications of the indole-alkaloid gramine modulate the assembly of individual members of the barley rhizosphere microbiota Microbial communities proliferating at the root-soil interface, collectively referred to as the rhizosphere microbiota, represent an untapped beneficial resource for plant growth, development and health. Integral to a rational manipulation of the microbiota for sustainable agriculture is the identification of the molecular determinants of these communities. In plants, biosynthesis of allelochemicals is centre stage in defining inter-organismal relationships in the environment. Intriguingly, this process has been moulded by domestication and breeding selection. The indole-alkaloid gramine, whose occurrence in barley (Hordeum vulgare L.) is widespread among wild genotypes but has been counter selected in several modern varieties, is a paradigmatic example of this phenomenon. This prompted us to investigate how exogenous applications of gramine impacted on the rhizosphere microbiota of two, gramine-free, elite barley varieties grown in a reference agricultural soil. High throughput 16S rRNA gene amplicon sequencing revealed that applications of gramine interfere with the proliferation of a subset of soil microbes with a relatively broad phylogenetic assignment. Strikingly, growth of these bacteria appeared to be rescued by barley plants in a genotype- and dosage-independent manner. In parallel, we discovered that host recruitment cues can interfere with the impact of gramine application in a host genotype-dependent manner. Interestingly, this latter effect displayed a bias for members of the phyla Proteobacteria. These initial observations indicate that gramine can act as a determinant of the prokaryotic communities inhabiting the root-soil interface. 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. Carbon Dot-based Fluorescent Antibody Nanoprobes as Brain Tumour Glioblastoma Diagnostics. The development of efficient and sensitive tools for the detection of brain cancer in patients is of the utmost importance particularly because many of these tumours go undiagnosed until the disease has advanced and when treatment is less effective. Current strategies employ antibodies (Abs) to detect Glial Fibrillary Acid Protein (GFAP) in tissue samples, since GFAP is unique to the brain and not present in normal peripheral blood, and it relies on fluorescent reporters.

Herein we describe a low cost, practical and general method for the labelling of proteins and antibodies with fluorescent carbon dots (CD) to generate diagnostic probes that are robust, photostable and applicable to the clinical setting. The two-step protocol relies on the conjugation of a dibenzocyclooctyne (DBCO)-functionalised CD with azide functionalised proteins by combining amide conjugation and strain promoted alkyne-azide cycloaddition (SPAAC) ligation chemistry. The new class of Ab-CD conjugates developed using this strategy was successfully used for the immunohistochemical staining of human brain tissues of patients with glioblastoma (GBM) validating the approach. Overall, these novel fluorescent probes offer a promising and versatile strategy in terms of costs, photostability and applicability which can be extended to other Abs and protein systems. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. Naming of human diseases on the wrong side of history BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. Epitopedia: identifying molecular mimicry between pathogens and known immune epitopes Upon infection, foreign antigenic proteins stimulate the hosts immune system to produce antibodies targeting the pathogen. These antibodies bind to regions on the antibody called epitopes. Structural similarity (molecular mimicry) of epitopes between an infecting pathogen and host proteins or other pathogenic proteins the host has previously encountered can impact the host immune response to the pathogen and may lead to cross-reactive antibodies. The ability to identify potential regions of molecular mimicry in a pathogen can illuminate immune effects which are especially important to pathogen treatment and vaccine design. Here we present Epitopedia, a software pipeline that facilitates the identification of regions that may exhibit potential three-dimensional molecular mimicry between an antigenic pathogen protein and known immune epitopes as catalogued by the immune epitope database (IEDB). Epitopedia is open-source software released under the MIT license and is freely available on GitHub, including a Docker container with all other software dependencies preinstalled. We performed an analysis describing how various secondary structure states, identity between pentapeptide pairs, and identity between the parent sequences of pentapeptide pairs affects RMSD. We found that pentapeptides pairs in a helical conformation had considerably lower RMSD values than those in Extended or Coil conformations. We also found that RMSD is significantly increased when pentapeptide pairs are from non-homologous sequences. GCN5L1 promotes diastolic dysfunction by inhibiting cardiac pyruvate oxidation Left ventricular diastolic dysfunction is a structural and functional condition that precedes the development of heart failure with preserved ejection fraction (HFpEF). The etiology of diastolic dysfunction includes alterations in fuel substrate metabolism that negatively impact cardiac bioenergetics, and may precipitate the eventual transition to heart failure. To date, the molecular mechanisms that regulate early changes in fuel metabolism leading to diastolic dysfunction remain unclear. In this report, we use a diet-induced obesity model and quantitative acetylproteomics in aged mice to show that inhibitory lysine acetylation of the pyruvate dehydrogenase (PDH) complex promotes energetic deficits and diastolic dysfunction in mouse hearts. Cardiomyocyte-specific deletion of the mitochondrial lysine acetylation regulatory protein GCN5L1 prevented hyperacetylation of the PDH complex subunit PDHA1, allowing aged obese mice to continue using pyruvate as a bioenergetic substrate in the heart. Our findings suggest that changes in mitochondrial protein lysine acetylation represent a key metabolic component of diastolic dysfunction that precedes the development of heart failure. Bacterial H-NS contacts DNA at the same irregularly spaced sites in both bridged and hemi-sequestered linear filaments Gene silencing in bacteria is mediated by chromatin proteins, of which Escherichia coli H-NS is a paradigmatic example. H-NS forms nucleoprotein filaments with either one or two DNA duplexes. However, the structures, arrangements of DNA-binding domains (DBDs), and positions of DBD-DNA contacts in linear and bridged filaments are uncertain. To characterize the contacts that silence transcription by RNA polymerase, we combined {middle dot}OH footprinting, molecular dynamics, statistical modeling, and DBD mapping using a chemical nuclease (Fe2+-EDTA) tethered to the DBDs (TEN-map). We find that H-NS DBDs contact DNA at indistinguishable locations in bridged or linear filaments and that the DBDs vary in orientation and position with ~10-bp average spacing. Our results support a hemi-sequestration model of linear-to-bridged H-NS switching in which linear filaments able to inhibit only transcription initiation switch to bridged filaments able to inhibit both initiation and elongation using the same irregularly spaced DNA contact sites.

HighlightsO_LITethered-nuclease mapping (TEN-map) of H-NS DNA-binding domains detects DNA contacts
C_LIO_LIBridged and linear H-NS filaments use the same DNA contact sites
C_LIO_LIH-NS-DNA contacts are unevenly spaced with ~10 bp average separation
C_LIO_LIAT-steps, minor groove width, and electrostatic potential best predict contact sites
C_LI Manipulating the microbiome alters regenerative outcomes in Xenopus laevis tadpoles via lipopolysaccharide signalling Xenopus laevis tadpoles can regenerate functional tails, containing spinal cord, notochord, muscle, fin, blood vessels and nerves, except for a brief refractory period at around one week of age. At this stage, amputation of the tadpoles tail may either result in scarless wound healing, or the activation of a regeneration programme, which replaces the lost tissues. We recently demonstrated a link between bacterial lipopolysaccharides and successful tail regeneration in refractory stage tadpoles, and proposed that this could result from lipopolysaccharides binding to Toll-like receptor 4 (TLR4). Here, we have used 16S rRNA sequencing to show that the tadpole skin microbiome is highly variable between sibships and that the community can be altered by raising embryos in the antibiotic gentamicin. Six gram-negative genera, including Delftia and Chryseobacterium, were over-represented in tadpoles that underwent tail regeneration. Lipopolysaccharides purified from a commensal Chryseobacterium spp. XDS4, an exogenous Delftia spp. or Escherichia coli could significantly increase the number of antibiotic-raised tadpoles that attempted regeneration. Conversely, the quality of regeneration was impaired in native-raised tadpoles exposed to the antagonistic lipopolysaccharide of Rhodobacter sphaeroides. Editing TLR4 using CRISPR/Cas9 also reduced regeneration quality, but not quantity, at the level of the cohort. However, we found that the editing level of individual tadpoles was a poor predictor of regenerative outcome. In conclusion, our results suggest that variable regeneration in refractory stage tadpoles depends at least in part on the skin microbiome and lipopolysaccharide signalling, but that signalling via TLR4 cannot account for all of this effect. Toxicity Profile of Artesunate in Rats and Dogs AO_SCPLOWBSTRACTC_SCPLOWO_ST_ABSOBJECTIVESC_ST_ABSThe objectives of these studies were to investigate the toxicity, safety and toxicokinetics of single and multiple doses of artesunate for injection in rats and dogs.

METHODSSprague-Dawley rats and Beagle dogs were treated intravenously or intramuscularly for 28 consecutive days with doses of up to 30 mg/kg artesunate, evaluating toxicity, kinetics, genotoxicity, and cardiovascular and central nervous safety parameters after single and 4-week repeated administrations. Furthermore, respiratory parameters were evaluated after a single intravenous administration in rats.

RESULTSArtesunate was well tolerated with no mortality and only minor effects on clinical pathology parameters. Following repeated intramuscular administration, local reactions at the injection site became evident.

CONCLUSIONSThe results obtained in these studies support the safe use of artesunate in humans.

Subchronic data for repeated intravenous and intramuscular artesunate in animals are available in the literature for a maximum duration of 14-days. Published results of genotoxicity studies in vitro in bacteria or cells of animal or human origin and in vivo in animals are equivocal, with a couple of negative results following oral administration and positive results in human cells.

With the studies performed in rats and dogs for a treatment period of 28-days including measures of safety pharmacology, toxicokinetic evaluation and in vivo genotoxicity data, a comparison of toxicity and exposure of intravenous versus intramuscular artesunate for a prolonged treatment period is possible, delivering important information, when a decision between both routes has to be made for use in patients. Results of the peripheral blood micronucleus assay in rats were positive, suggesting genotoxic potential for the intravenous route of administration of artesunate, what is, beside patients, likewise relevant for production and health care personnel. Changes to the mtDNA copy number during yeast culture growth We show that the mtDNA copy number in growing cultures of the yeast Saccharomyces cerevisiae increases by a factor of up to 4, being lowest ([~]10 per haploid genome) and stable during rapid fermentative growth, and highest at the end of the respiratory phase. When yeast are grown on glucose, the onset of the mtDNA copy number increase coincides with the early stages of the diauxic shift, and the increase continues through respiration. A lesser yet still substantial copy number increase occurs when yeast are grown on a nonfermentable carbon source, i.e. when there is no diauxic shift. The mtDNA copy number increase during and for some time after the diauxic shift is not driven by an increase in cell size. The copy number increase occurs in both haploid and diploid strains, but is markedly attenuated in a diploid wild isolate that is a ready sporulator. Strain-to-strain differences in mtDNA copy number are least apparent in fermentation and most apparent in late respiration or stationary phase. While changes in mitochondrial morphology and function were previously known to accompany changes in physiological state, it had not been previously shown that the mtDNA copy number changes substantially over time in a clonal growing culture. The mtDNA copy number in yeast is therefore a highly dynamic phenotype. Environment dependent benefits of inter-individual variation in honey bee recruitment Inter-individual differences in behaviour within the members of a social group can affect the groups productivity. In eusocial insects, individual differences amongst workers in a colony play a central role division of labour and task allocation. Extensive empirical and theoretical work has highlighted variation in response thresholds as a proximate mechanism underlying individual behavioural differences and hence division of labour. However, other response parameters, like response probability and intensity, can affect these differences. In this study, we first extended a previously published agent-based model on honey bee foraging to understand the relative importance of response (dance) probability and response (dance) intensity in the task of recruitment. Comparing variation obtained from the simulations with previously published empirical data, we found that response intensity plays a more important role than probability in producing consistent inter-individual differences in recruitment behaviour. We then explored the benefits provided by this individual variation in recruitment behaviour to the colonys collective foraging effort under different environmental conditions. Our results revealed that individual variation leads to a greater energetic yield per forager, but only when food is abundant. Our study highlights the need to consider all response parameters while studying division of labour and adds to the growing body of evidence linking individual variation in behavioural responses to the success of social groups. Advanced paternal age diversifies individual trajectories of vocalization patterns in neonatal mice Infant crying is an innate communicative behavior that is frequently impaired in certain neurodevelopmental disorders (NDDs). Since advanced paternal age is a reported risk factor for NDDs in offspring, we evaluated the impact of a fathers age on early vocal development in C57BL/6J mice. We recorded and applied a unique combination of computational analyses to ultrasonic vocalizations (USVs) emitted by mouse pups sired by young and aged fathers. Our data showed that advanced paternal age reduced the number and duration of USVs, and altered the syllable composition in pups. Moreover, pups born to young fathers showed convergent vocal characteristics with a rich repertoire during postnatal development, while those born to aged fathers exhibited more divergent vocal patterns with limited repertoire. Principal component analysis in conjunction with clustering analysis demonstrated that pups from aged fathers deviated from typical trajectories of vocal development, which were considered as atypical individuals. Thus, our study indicates that advanced paternal age has a significant effect on offsprings early vocal development. It is suggested that the trajectories of vocal development could be a useful marker of the NDD-like phenotype associated with the advanced paternal age. In addition, our comprehensive computational analysis described here is an effective approach to characterize the altered individual diversity relevant to neurodevelopmental disorders.

One Sentence SummaryAdvanced paternal age affects vocal development in early postnatal mice, with more pups showing atypical developmental trajectories. Macroscopic control of synchronous electrical signaling with chemically-excited gene expression Excitable cells can convert electrical signals into chemical outputs to facilitate the active transport of information across larger distances. This electrical-to-chemical conversion requires a tightly regulated expression of ion channels. Alterations of ion channel expression provide landmarks of numerous pathological diseases, such as cardiac arrhythmia, epilepsy, or cancer. Although the activity of ion channels can be locally regulated by external light or chemical stimulus, it remains challenging to coordinate the expression of ion channels on extended spatial-temporal scales in a non-invasive manner. Here, we have engineered yeast S. cerevisiae to read and convert local chemical concentrations into a dynamic electrical field distributed across cell populations. The core mechanism encodes a chemically-excitable dual-feedback gene circuit that precisely tunes the expression domain of potassium channels, globally coordinating cyclic firing of the plasma membrane potential (PMP). We demonstrate that this mechanism leverages an engineered constitutively open bacterial potassium channel KcsA to directly couple chemical stimuli with ion flux through gene expression and it can interface with the host ion channels through the pulsatile production of toxins. Our study provides a robust synthetic transcriptional toolbox underlying the conversion of local chemical environments into spatiotemporally organized electrical impulses for various cellular engineering, synthetic biology, and potential therapeutic applications. Kinetics and Mechanism of Fentanyl Dissociation from the mu-Opioid Receptor Driven by illicit fentanyl, opioid related deaths have reached the highest level in 2020. Currently, an opioid over-dose is resuscitated by the use of naloxone, which competitively binds and antagonizes the {micro}-opioid receptor (mOR). Thus, knowledge of the residence times of opioids at mOR and the unbinding mechanisms is valuable for assessing the effectiveness of naloxone. In the present study we calculate the fentanyl-mOR dissociation time and elucidate the mechanism by applying an enhanced sampling molecular dynamics (MD) technique. Two sets of metadynamics simulations with different initial structures were performed while accounting for the protonation state of the conserved H2976.52, which has been suggested to modulate the ligand-mOR affinity and binding mode. Surprisingly, with the N{delta}-protonated H2976.52, fentanyl can descend as much as 10 [A] below the level of the conserved D1473.32 before escaping the receptor, and has a calculated residence time{tau} of 38 s. In contrast, with the N - and doubly protonated H2976.52, the calculated{tau} are 2.6 s and 0.9 s, respectively. Analysis suggests that formation of the piperidine-Hid297 hydrogen bond strengthens the hydrophobic contacts with the transmembrane helix (TM) 6, allowing fentanyl to explore a deep pocket. Considering the experimental{tau} of [~]4 min for fentanyl and the role of TM6 in mOR activation, we suggest that the deep insertion mechanism may be biologically relevant. The work paves the way for large-scale computational predictions of opioid dissociation rates to inform the regulatory decision. The profound role of the histidine protonation state found here may shift the paradigm in computational studies of ligand-receptor kinetics. To knot or not to knot: Multiple conformations of the SARS-CoV-2 frameshifting RNA element The SARS-CoV-2 frameshifting RNA element (FSE) is an excellent target for therapeutic intervention against Covid-19. This small gene element employs a shifting mechanism to pause and backtrack the ribosome during translation between Open Reading Frames 1a and 1b, which code for viral polyproteins. Any interference with this process has profound effect on viral replication and propagation. Pinpointing the structures adapted by the FSE and associated structural transformations involved in frameshifting has been a challenge. Using our graph-theory-based modeling tools for representing RNA secondary structures, "RAG" (RNA-As-Graphs), and chemical structure probing experiments, we show that the 3-stem H-type pseudoknot (3_6 dual graph), long assumed to be the dominant structure has a viable alternative, an HL-type 3-stem pseudoknot (3_3) for longer constructs. In addition, an unknotted 3-way junction RNA (3_5) emerges as a minor conformation. These three conformations share Stems 1 and 3, while the different Stem 2 may be involved in a conformational switch and possibly associations with the ribo-some during translation. For full-length genomes, a stem-loop motif (2_2) may compete with these forms. These structural and mechanistic insights advance our understanding of the SARS-CoV-2 frameshifting process and concomitant virus life cycle, and point to three avenues of therapeutic intervention. 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. The effect of resource depletion on the thermal response of mosquito population fitness The temperature-dependencies of life history traits are increasingly being used to predict how climatic warming will affect vector-borne disease dynamics, partially by affecting the abundance dynamics of the vector population. Such predictions generally arise from mathematical models that incorporate the temperature dependence of traits measured under laboratory conditions. These temperature-trait relationships are typically estimated from juvenile populations reared under optimal resource conditions, even though natural populations experience intermittent resource depletion. Using laboratory experiments on the mosquito Aedes aegypti, combined with a stage-structured population model, we show that resource depletion in the juvenile habitat can significantly depress the vectors maximal population growth rate (rm) across the entire temperature range, cause it to peak at a lower temperature, and narrow its thermal niche width. Our results provide compelling evidence for future studies to consider resource depletion when predicting the effects of global change on vector-borne disease transmission, disease vectors and other arthropods. The activity and functions of soil microbial communities in the Finnish sub-Arctic vary across vegetation types Increased microbial activity in high-latitude soils due to climate change might lead to higher greenhouse gas (GHG) emissions. However, mechanisms of microbial GHG production and consumption in tundra soils are not thoroughly understood. We analyzed 116 soil metatranscriptomes from 73 sites in the Finnish sub-Arctic to investigate how the diversity and functional potential of bacterial and archaeal communities vary across vegetation types and soil layers. Soils differed in physicochemical conditions, with meadow soils being characterized by higher pH and low soil organic matter (SOM) and carbon/nitrogen ratio whereas dwarf shrub-dominated ecosystems with high SOM and low pH. Actinobacteria, Acidobacteria, Alphaproteobacteria, and Planctomycetes predominated all communities but there were significant differences on genus level between vegetation types, as plant polymer degrading groups were more active in shrub-dominated soils compared to meadows. Given that climate change scenarios predict expansion in dwarf shrubs at high latitudes, our results indicate that the rate of carbon turnover in tundra soils may increase in the future. Additionally, transcripts of methanotrophs were detected in the mineral layer of all soils, potentially moderating methane fluxes from deeper layers. In all, this study provides new insights into possible shifts in tundra microbial diversity and activity with climate change. The important choice of reference environment in microevolutionary climate response predictions It is well documented that individuals of wild populations can adjust to climate change by means of phenotypic plasticity, but few reports on adaptation by means of genetically based microevolution caused by selection. Disentanglement of these separate effects requires that the reference environment (the environmental zero-point) is defined, and this should not be done arbitrarily. The problem is that an error in the reference environment may lead to large errors in predicted microevolution. Together with parameter values and initial mean trait values, the reference environment can be estimated from environmental, phenotypic and fitness data. A prediction error method for this purpose is described, with the feasibility shown by simulations. As shown in a toy example, an estimated reference environment may have large errors, especially for small populations. This may still be a better choice than use of an initial environmental value in a recorded time series, or the mean value, which is often used. Another alternative may be to use the mean value of a past and stationary stochastic environment, which the population is judged to have been fully adapted to, in the sense that the expected geometric mean fitness was at a global maximum. Exceptions are cases with constant phenotypic plasticity, where the microevolutionary changes per generation follow directly from phenotypic and environmental data, independent of the chosen reference environment. DNA damage models phenotypes of beta cell senescence in Type 1 Diabetes ObjectiveType 1 Diabetes (T1D) is characterized by progressive loss of insulin-producing pancreatic {beta} cells as a result of autoimmune destruction. In addition to {beta} cell death, recent work has shown that subpopulations of {beta} cells acquire dysfunction during T1D. We previously reported that some {beta} cells adopt a senescent fate involving a DNA damage response (DDR) during the pathogenesis of T1D, however, the question of how senescence develops in {beta} cells has not been investigated.

MethodsHere, we tested the hypothesis that unrepaired DNA damage triggers {beta} cell senescence using culture models including the mouse NIT1 {beta} cell line derived from the T1D-susceptible nonobese diabetic (NOD) strain, human donor islets and EndoC {beta} cells. DNA damage was chemically induced using etoposide or bleomycin and cells or islets were analyzed by a combination of molecular assays for senescence phenotypes including Western blotting, qRT-PCR, Luminex assays, flow cytometry and histochemical staining. RNA-seq was carried out to profile global transcriptomic changes in human islets undergoing DDR and senescence. Insulin ELISAs were used to quantify glucose stimulated insulin secretion from chemically-induced senescent islets and cells in culture.

ResultsSub-lethal DNA damage in NIT1 cells led to several classical hallmarks of senescence including sustained DDR activation, growth arrest, enlarged flattened morphology and a senescence-associated secretory phenotype (SASP) resembling what occurs in primary {beta} cells during T1D in NOD mice. Some of these phenotypes differed between NIT1 cells and the MIN6 {beta} cell line derived from a non-T1D susceptible mouse strain. RNA-seq analysis of human islets undergoing DDR and early senescence revealed a coordinated p53-p21 transcriptional program and upregulation of prosurvival signaling and SASP genes, which were confirmed in independent islet preparations at the protein level. Importantly, chemically induced DNA damage also led to DDR activation and senescence phenotypes in the EndoC-{beta}H5 human {beta} cell line, confirming that this response can occur directly in human {beta} cells. Finally, DNA damage and senescence in both mouse {beta} cell lines and human islets led to decreased insulin content.

ConclusionsTaken together, these findings suggest that some of the phenotypes of {beta} cell senescence during T1D can be modeled by chemically induced DNA damage in mouse {beta} cell lines and human islets and {beta} cells in culture. These culture models will be useful tools to understand some of the mechanisms of {beta} cell senescence in T1D.

HighlightsO_LIDNA damage induces senescent phenotypes in mouse {beta} cell lines
C_LIO_LIDNA damage induces a p53-p21 transcriptional program and senescent phenotypes in human islets and EndoC cells
C_LIO_LIDNA damage and senescence leads to decreased insulin content
C_LIO_LIDNA damage models some aspects of {beta} cell senescence in Type 1 Diabetes
C_LI Substrate specialization in microbes is driven by biochemical constraints of dynamic flux sensing Central carbon metabolism is highly conserved across microbial species, but can catalyze very different pathways depending on the organism and their ecological niche. Here, we study the dynamic re-organization of central metabolism after switches between the two major opposing pathway configurations of central carbon metabolism, glycolysis and gluconeogenesis in Escherichia coli, Pseudomonas aeruginosa and Pseudomonas putida. We combined growth dynamics and dynamic changes of intracellular metabolite levels with a coarse-grained model that integrates fluxes, regulation, protein synthesis and growth and uncovered fundamental limitations of the regulatory network: after nutrient shifts, metabolite concentrations collapse to their equilibrium, rendering the cell unable to sense which direction the flux is supposed to flow through the metabolic network. The cell can partially alleviate this by picking a preferred direction of regulation at the expense of increasing lag times in the opposite direction. Moreover, decreasing both lag times simultaneously comes at the cost of reduced growth rate or higher futile cycling between metabolic enzymes. These three trade-offs can explain why microorganisms specialize for either glycolytic or gluconeogenic substrates and can help elucidate the complex growth patterns exhibited by different microbial species.

Graphical synopsis O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

Standfirst textMicrobes face a series of fundamental trade-offs that limit their ability to optimize simultaneously for both glycolytic and gluconeogenic growth.

Bullet pointsO_LILag times between glycolysis and gluconeogenesis show asymmetry in many microbes: A long lag in one direction, but a short lag in the other.
C_LIO_LILong lag times are caused by an inability to sense fluxes after nutrient shifts.
C_LIO_LIWith existing regulation, lag time asymmetry can only be overcome by reducing either growth rate or increasing futile cycling in metabolism.
C_LI