Imaging sensory transmission and neuronal plasticity in primary sensory neurons with genetically-encoded voltage indicator, ASAP4.4-Kv Detection of somatosensory inputs requires conversion of external stimuli into electrical signals by activation of primary sensory neurons. The mechanisms by which heterogeneous primary sensory neurons encode different somatosensory inputs remains unclear. In vivo dorsal root ganglia (DRG) imaging using genetically-encoded Ca2+ indicators (GECIs) is currently the best technique for this purpose by providing an unprecedented spatial and populational resolution. It permits the simultaneous imaging of >1800 neurons/DRG in live mice. However, this approach is not ideal given that Ca2+ is a second messenger and has inherently slow response kinetics. In contrast, genetically-encoded voltage indicators (GEVIs) have the potential to track voltage changes in multiple neurons in real time but often lack the brightness and dynamic range required for in vivo use. Here, we used soma-targeted ASAP4.4-Kv, a novel GEVI, to dissect the temporal dynamics of noxious and non-noxious neuronal signals during mechanical, thermal, or chemical stimulation in DRG of live mice. ASAP4.4-Kv is sufficiently bright and fast enough to optically characterize individual neuron coding dynamics. Notably, using ASAP4.4-Kv, we uncovered cell-to-cell electrical synchronization between adjacent DRG neurons and robust dynamic transformations in sensory coding following tissue injury. Finally, we found that a combination of GEVI and GECI imaging empowered in vivo optical studies of sensory signal processing and integration mechanisms with optimal spatiotemporal analysis.

HighlightsO_LIIn vivo ultra fast and sensitive dynamic voltage imaging of peripheral primary sensory neurons by a newly generated genetically-encoded voltage indicator.
C_LIO_LIIdentification of mechanical, thermal, or chemical stimuli-evoked voltage signals with superior temporal resolution.
C_LIO_LISingle-cell detection of changes in sub- and suprathreshold voltage dynamics across different disease conditions.
C_LIO_LICombination of voltage (by ASAP4.4-Kv) and Ca2+ (by Pirt-GCaMP3) signals to facilitate the understanding of signal processing and integration of primary sensory neurons, especially for noxious versus non-noxious sensation.
C_LI
Reset Networks: Emergent Topography by Composition of Convolutional Neural Networks We introduce the Reset model, a composition of neural networks - typically several levels of convolutional neural networks - whose outputs at one level are gathered and reshaped into a spatial input for the next level. We show that units in Reset networks self-organize into clusters when trained on MNIST, Fashion MNIST, CIFAR-10 and CIFAR-100. We then show that a stronger type of self-organization, reminiscent of the topography found for numbers in parietal cortex, arises when number images are mapped onto developmentally realistic number codes. We outline the implications of this model for theories of the cortex and developmental neuroscience.
Contrastive learning on protein embeddings enlightens midnight zone Experimental structures are leveraged through multiple sequence alignments, or more generally through homology-based inference (HBI), facilitating the transfer of information from a protein with known annotation to a query without any annotation. A recent alternative expands the concept of HBI from sequence-distance lookup to embedding-based annotation transfer (EAT). These embeddings are derived from protein Language Models (pLMs). Here, we introduce using single protein representations from pLMs for contrastive learning. This learning procedure creates a new set of embeddings that optimizes constraints captured by hierarchical classifications of protein 3D structures defined by the CATH resource. The approach, dubbed ProtTucker, has an improved ability to recognize distant homologous relationships than more traditional techniques such as threading or fold recognition. Thus, these embeddings have allowed sequence comparison to step into the "midnight zone" of protein similarity, i.e., the region in which distantly related sequences have a seemingly random pairwise sequence similarity. The novelty of this work is in the particular combination of tools and sampling techniques that ascertained good performance comparable or better to existing state-of-the-art sequence comparison methods. Additionally, since this method does not need to generate alignments it is also orders of magnitudes faster. The code is available at https://github.com/Rostlab/EAT. Increasing triacylglycerol formation and lipid storage by unsaturated lipids protects renal proximal tubules in diabetes In diabetic patients, dyslipidemia frequently contributes to organ damage such as diabetic kidney disease (DKD). DKD is associated with excessive renal deposition of triacylglycerol (TAG) in lipid droplets (LD). Yet, it is unclear whether LDs play a protective or damaging role and how this might be influenced by dietary patterns. By using a diabetes mouse model, we find here that high fat diet enriched in the unsaturated oleic acid (OA) caused more lipid storage in LDs in renal proximal tubular cells (PTC) but less tubular damage than a corresponding butter diet with the saturated palmitic acid (PA). Mechanistically, we identify endoplasmic reticulum (ER) stress as the main cause of PA-induced PTC injury. ER stress is caused by elevated cellular levels of saturated TAG precursors and to higher membrane order in the ER. The resulting cell death is preceded by a transcriptional rewiring of phospholipid metabolism. Simultaneous addition of OA rescues the cytotoxic effects by normalizing membrane order and by increasing the total TAG amount. The latter also stimulates the formation of LDs that in turn can release unsaturated lipids upon demand by lipolysis. Our study thus clarifies mechanisms underlying PA-induced cell stress in PTCs and emphasizes the importance of olive oil for the prevention of DKD. Artificial seawater based long-term culture of colonial ascidians Tunicates are highly diverse marine invertebrate filter-feeders that are vertebrates closest relatives. These organisms, despite a drastically different body plan during their adulthood, have a tissue complexity related to that of vertebrates. Ascidians, which compose most of the Tunicata, are benthic sessile hermaphrodites that reproduce sexually through a motile tadpole larval stage. Over half of the known ascidians species are able to reproduce asexually by budding, typically leading to the formation of colonies where animals, called zooids, are interconnected through an external vascular system. In addition, colonial ascidians are established models for important biological processes including allorecognition, immunobiology, aging, angiogenesis and whole-body regeneration. However, the current paucity in breeding infrastructures limits the study of these animals to coastal regions.

To promote a wider scientific spreading and popularity of colonial ascidians, we have developed a flexible recirculating husbandry setup for their long-term in-lab culture. Our system is inspired both by the flow-through aquariums used by coastal ascidian labs, as well as by the recirculating in-lab systems used for zebrafish research. Our hybrid system thus combines colony breeding, water filtering and food culturing in a semi-automated system where specimens develop on hanging microscopy glass slides. Temperature, light/dark cycles, flow speed and feeding rates can be controlled independently in four different breeding environments to provide room for species-specific optimization as well as for running experiments. This setup is complemented with a quarantine for the acclimatization of wild isolates.

Herein we present our success in breeding Botrylloides diegensis, a species of colonial ascidians, for more than 3 years in recirculating artificial seawater over 600 km away from their natural habitat. We show that colonies adapt well to in-lab culturing provided that a suitable marine microbiome is present, and that a specific strain can be isolated, propagated and efficiently used for research over prolonged periods of time. The flexible and modular structure of our system can be scaled and adapted to the needs of specific species, such as Botryllus schlosseri, as well as of particular laboratory spaces. Overall, we show that Botrylloides diegensis can be proficiently bred in-land and suggest that our results can be extended to other species of colonial ascidians to promote research on these fascinating animals.

HighlightsO_LIFirst artificial seawater based recirculating aquaculture for colonial ascidians
C_LIO_LIOver 3 years of continuous breeding
C_LIO_LISemi-automated setup with minimized maintenance
C_LIO_LIGood biomass production for strain propagation
C_LIO_LI4 different culture conditions for optimized breeding for species of interest
C_LI Protein-protein interactions on membrane surfaces analysed using pull-downs with supported bilayers on silica beads Discovery-based proteomics workflows that identify novel interactors rely on immunoprecipitations or pull-downs with genetically-tagged bait proteins immobilized on appropriate matrices. But strategies to analyse protein interactions on a diffusible membrane surface combined with the practical ease of pull-downs remain unavailable. Such strategies are important to analyse protein complexes that mature in composition and stability because of diffusion-based encounter between participant proteins. Here, we describe a generic pull-down strategy to analyse such complexes using chelating lipid-containing supported bilayers formed on silica beads. These templates can display desired His-tagged bait proteins on a diffusible membrane surface. Using clathrin-mediated endocytosis as a paradigm, we find that the clathrin-binding adaptor protein epsin1 displayed as bait on these templates pulls down significantly higher amounts of clathrin from brain lysates than when immobilized on conventional matrices. Together, our results establish the potential of such templates as superior matrices for analysing protein-protein interactions and resultant complexes formed on membrane surfaces. Heterogeneous recruitment abilities to RNA polymerases generate nonlinear scaling of gene expression level with cell volume While most genes expression levels are proportional to cell volumes, some genes exhibit nonlinear scaling between their expression levels and cell volume. Therefore, their mRNA and protein concentrations change as the cell volume increases, which often have crucial biological functions such as cell-cycle regulation. However, the biophysical mechanism underlying the nonlinear scaling between gene expression and cell volume is still unclear. In this work, we show that the nonlinear scaling is a direct consequence of the heterogeneous recruitment abilities of promoters to RNA polymerases based on a gene expression model at the whole-cell level. Those genes with weaker (stronger) recruitment abilities than the average ability spontaneously exhibit superlinear (sublinear) scaling with cell volume. Analysis of the promoter sequences and the nonlinear scaling of Saccharomyces cerevisiaes mRNA levels shows that motifs associated with transcription regulation are indeed enriched in genes exhibiting nonlinear scaling, in concert with our model. Survival Genie, a web platform for survival analysis across pediatric and adult cancers. The genomics data-driven identification of gene signatures and pathways has been routinely explored for predicting cancer survival and making decisions related to targeted treatments. A large number of packages and tools have been developed to correlate gene expression/mutations to the clinical outcome but lack the ability to perform such analysis based on pathways, gene sets, and gene ratios. Furthermore, in this single-cell omics era, the cluster markers from cancer single-cell transcriptomics studies remain an underutilized prognostic option. Additionally, no bioinformatics online tool evaluates the associations between the enrichment of canonical cell types and survival across cancers. Here we have developed Survival Genie, a web tool to perform survival analysis on single-cell RNA-seq (scRNA-Seq) data and a variety of other molecular inputs such as gene sets, genes ratio, tumor infiltrating immune cells proportion, gene expression profile scores, and tumor mutation burden. For a comprehensive analysis, Survival Genie contains 53 datasets of 27 distinct malignancies from 11 different cancer programs related to adult and pediatric cancers. Users can upload scRNA-Seq data or gene sets and select a gene expression partitioning method (i.e., mean, median, quartile, cutp) to determine the effect of expression levels on survival outcomes. The tool provides comprehensive results including box plots of low and high-risk groups, Kaplan-Meier plots with univariate Cox proportional hazards model, and correlation of immune cell enrichment and molecular profile.

The analytical options and comprehensive collection of cancer datasets make Survival Genie a unique resource to correlate gene sets, pathways, cellular enrichment, and single-cell signatures to clinical outcomes to assist in developing next-generation prognostic and therapeutic biomarkers. Survival Genie is open-source and available online at https://bbisr.shinyapps.winship.emory.edu/SurvivalGenie/. Protective host-dependent antagonism among 	Pseudomonas in the Arabidopsis phyllosphere The plant microbiome is a rich biotic environment, comprising numerous taxa. The community structure of these colonizers is constrained by multiple factors, including host-microbe and microbe-microbe interactions, as well as the interplay between the two. While much can be learned from pairwise relationships between individual hosts and microbes, or individual microbes with themselves, the ensemble of interrelations between the host and microbial consortia may lead to different outcomes that are not easily predicted from the individual interactions. Their study can thus provide new insights into the complex relationship between plants and microbes. Of particular importance is how strain-specific such plant-microbe-microbe interactions are, and how they eventually affect plant health. Here, we test strain-level interactions in the phyllosphere between groups of co-existing commensal and pathogenic Pseudomonas among each other and with A. thaliana, by employing synthetic communities of genome-barcoded isolates. We found that commensal Pseudomonas prompted a host response leading to a selective inhibition of a specific pathogenic lineage, resulting in plant protection. The extent of plant protection, however, was dependent on plant genotype, indicating that these effects are host-mediated. There were similar genotype-specific effects on the microbe side, as we could pinpoint an individual Pseudomonas isolate as the predominant cause for this differential interaction. Collectively, our work highlights how within-species genetic differences on both the host and microbe side can have profound effects on host-microbe-microbe dynamics. The paradigm that we have established provides a platform for the study of host-dependent microbe-microbe competition and cooperation in the A. thaliana-Pseudomonas system. 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. Substrate traits shape the structure of microbial community engaged in metabolic division of labor Metabolic division of labor (MDOL) is widespread in nature, whereby a complex metabolic pathway is shared between different strains within a community for mutual benefit. However, little is known about how the mutual interactions in the microbial community engaged in MDOL are regulated. We hypothesized that when degradation of an organic compound is carried out via MDOL, the substrate traits (i.e., concentration and its toxicity) modulate the benefit allocation between the two microbial populations, thus affecting the structure of this community. We tested this hypothesis by combining mathematical modelling with experiments using engineered synthetic microbial consortia. Numerous modelling analyses suggested that the proportion of the population executing the first metabolic step can be simply estimated by Monod-like formulas governed by substrate traits. The model and the proposed formula quantitatively predicted the structure of our synthetic consortia composed of two strains degrading salicylate through MDOL. Individual-based modelling and colony pattern formation assays further indicated that our rule is also applicable to estimating community structure in spatially structured environments. Our results demonstrate that the structure of the microbial communities can be quantitatively predicted from simple environmental factors, such as substrate concentration and its toxicity, which provides novel perspectives on understanding the assembly of natural communities, as well as insights into how to manage artificial microbial systems. Neuroinvasion and neurotropism by SARS-CoV-2 variants in the K18-hACE2 mouse Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) not only affects the respiratory tract but also causes neurological symptoms such as loss of smell and taste, headache, fatigue or severe cerebrovascular complications. Using transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) we investigated the spatiotemporal distribution and pathomorphological features in the CNS following intranasal infection with SARS-CoV-2 variants, also after prior influenza A virus infection. Apart from Omicron, we found all variants to frequently spread to and within the CNS. Infection was restricted to neurons and appeared to spread from the olfactory bulb mainly in basally orientated regions in the brain and into the spinal cord, independent of ACE2 expression and without evidence of neuronal cell death, axonal damage or demyelination. However, microglial activation, microgliosis and a mild macrophage and T cell dominated inflammatory response was consistently observed, accompanied by apoptotic death of endothelial, microglial and immune cells, without their apparent infection. Microgliosis and immune cell apoptosis indicate a potential role of microglia for pathogenesis and viral effect in COVID-19 and possible impairment of neurological functions, especially in long COVID. These data may also be informative for the selection of therapeutic candidates, and broadly support investigation of agents with adequate penetration into relevant regions of the CNS. 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 perceptual relevance of the speech envelope-following response in normal hearing and age-related hearing impaired listeners The aim of this work was to investigate the perceptual relevance of the frequency following response to the syllable /da/ for speech intelligibility in noise based on age and hearing deficits. Recordings of the auditory evoked potential from young normal hearing (NH) and older individuals with both normal hearing and high-frequency (HF) hearing loss were analyzed. EFR metrics obtained in quiet and noise condition were calculated and correlated with speech reception. The envelope following responses were analyzed in terms of amplitude, latency and noise robustness. The response was first simulated to form predictions on the effect of cochlear synaptopathy and outer hair cell loss on the EFR. The experimental findings were in line with the computational predictions in the found observation that the EFR was reduced as a consequence of ageing and HF hearing loss. Both the audiogram and the speech EFR magnitude fell short in the individual prediction of SRT in stationary noise, but they accounted well for group performance. We also obtained within-group EFR latency with a cross covariance matrix. Validation of the method confirmed that speech EFR latency was predictive of click ABR Wave V peak latency. Moreover, statistical analysis not only showed that the robustness of the EFR obtained in the noise condition was dependent on the degree of high-frequency hearing loss in the older NH adults, but also dependent on the EFR magnitude in the NH younger adults. These findings provide evidence towards the important role of the EFR in speech-in-noise perception. Relative contributions of sex hormones, sex chromosomes, and gonads to sex differences in tissue gene regulation Sex differences in physiology and disease in mammals result from the effects of three classes of factors that are inherently unequal in males and females: reversible (activational) effects of gonadal hormones, permanent (organizational) effects of gonadal hormones, and cell-autonomous effects of sex chromosomes, as well as genes driven by these classes of factors. Often, these factors act together to cause sex differences in specific phenotypes, but the relative contribution of each and the interactions among them remain unclear. Here, we used the Four Core Genotypes (FCG) mouse model with or without hormone replacement to distinguish the effects of each class of sex-biasing factors on transcriptome regulation in liver and adipose tissues. We found that the activational hormone levels have the strongest influence on gene expression, followed by the organizational gonadal sex effect and, lastly, sex chromosomal effect, along with interactions among the three factors. Tissue specificity was prominent, with a major impact of estradiol on adipose tissue gene regulation, and of testosterone on the liver transcriptome. The networks affected by the three sex-biasing factors include development, immunity and metabolism, and tissue-specific regulators were identified for these networks. Furthermore, the genes affected by individual sex-biasing factors and interactions among factors are associated with human disease traits such as coronary artery disease, diabetes, and inflammatory bowel disease. Our study offers a tissue-specific account of the individual and interactive contributions of major sex-biasing factors to gene regulation that have broad impact on systemic metabolic, endocrine, and immune functions. The "Alu-ome" shapes the epigenetic environment of regulatory elements controlling cellular defense Promoters and enhancers are sites of transcription initiation (TSSs) and carry active histone modifications, including H3K4me1, H3K4me3, and H3K27ac. Yet, the principles governing the boundaries of such regulatory elements are still poorly characterized. Alu elements are good candidates for a boundary function, being highly abundant in gene-rich regions, while essentially excluded from regulatory elements. Here, we show that the interval from the TSS to the first upstream Alu accommodates all H3K4me3 and most H3K27ac marks, while excluding DNA methylation. Remarkably, the average length of this intervals greatly varies in-between tissues, being longer in stem-and shorter in immune-cells. The very shortest TSS-to-Alu intervals were observed at promoters active in T cells, particularly at immune genes, correlating with RNA polymerase II transcription through the first Alu and accumulation of H3K4me1 signal on this first Alu. Finally, DNA methylation at first-Alus was found to evolved with age, regressing from young to middle-aged, then recovering later in life. Thus, the first Alus upstream of TSSs appear as dynamic boundaries marking the transition from DNA methylation to active histone modifications at regulatory elements, while also participating in the recording of immune gene transcriptional events by positioning H3K4me1-modified nucleosomes. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Severe COVID-19 associated variants linked to chemokine receptor gene control in monocytes and macrophages Genome-wide association studies have identified 3p21.31 as the main risk locus for severe disease in COVID-19 patients, although underlying biological mechanisms remain elusive. We performed a comprehensive epigenomic dissection of the 3p21.31 locus, identifying a CTCF-dependent tissue-specific 3D regulatory chromatin hub that controls the activity of several tissue-homing chemokine receptor (CCR) genes in monocytes and macrophages. Risk SNPs colocalized with regulatory elements and were linked to increased expression of CCR1, CCR2 and CCR5 in monocytes and macrophages. As excessive organ infiltration of inflammatory monocytes and macrophages is a hallmark of severe COVID-19, our findings provide a rationale for the genetic association of 3p21.31 variants with elevated risk of hospitalization upon SARS-CoV-2 infection. The compound of baicalein, wogonin, and oroxylin-A inhibits EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. Modeling microRNA-driven post-transcriptional regulation by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Trade-offs in biosensor optimization for dynamic pathway engineering Recent progress in synthetic biology allows the construction of dynamic control circuits for metabolic engineering. This technology promises to overcome many challenges encountered in traditional pathway engineering, thanks to their ability to self-regulate gene expression in response to bioreactor perturbations. The central components in these control circuits are metabolite biosensors that read out pathway signals and actuate enzyme expression. However, the construction of metabolite biosensors is a major bottleneck for strain design, and a key challenge is to understand the relation between biosensor dose-response curves and pathway performance. Here we employ multiobjective optimization to quantify performance trade-offs that arise in the design and calibration of metabolite biosensors. Our approach reveals strategies for tuning dose-response curves along an optimal trade-off between production flux and the cost of an increased expression burden on the host. We explore properties of control architectures built in the literature, and identify their advantages and caveats in terms of performance and robustness to growth conditions and leaky promoters. We demonstrate the optimality of a control circuit for glucaric acid production in Escherichia coli, which has been shown to increase titer by 2.5-fold as compared to static designs. Our results lay the groundwork for the automated design of control circuits for pathway engineering, with applications in the food, energy and pharmaceutical sectors. Decreased Interfacial Dynamics Caused by the N501Y Mutation in the SARS-CoV-2 S1 Spike:ACE2 Complex Corona Virus Disease of 2019 (COVID-19) caused by Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2) has caused a massive health crisis across the globe, with some genetic variants gaining enhanced infectivity and competitive fitness, and thus significantly aggravating the global health concern. In this regard, the recent SARS-CoV-2 alpha variant, B.1.1.7 lineage, reported from the United Kingdom (UK), is of great significance in that it contains several mutations that increase its infection and transmission rates as evident from clinical reports. Specifically, the N501Y mutation in the SARS-CoV-2 spike S1 receptor binding domain (S1-RBD) has been shown to possess an increased affinity for ACE2, although the basis for this is not entirely clear yet. Here, we dissect the mechanism underlying the increased affinity using molecular dynamics (MD) simulations of the available ACE2-S1-RBD complex structure (6M0J) and show a prolonged and stable interfacial interaction of the N501Y mutant S1-RBD with ACE2 compared to the wild type S1-RBD. Additionally, we find that the N501Y mutant S1-RBD displays altered dynamics that likely aids in its enhanced interaction with ACE2. By elucidating a mechanistic basis for the increased affinity of the N501Y mutant S1-RBD for ACE2, we believe that the results presented here will aid in developing therapeutic strategies against SARS-CoV-2 including designing drugs targeting the ACE2-S1-RBD interaction.

SignificanceThe emergence of the new SARS-CoV-2 lineage in the UK in December 2020 has further aggravated the COVID-19 pandemic due to an increased ability of the variant to infect human hosts, likely due to mutations in the viral S1 spike protein including the N501Y S1-RBD mutation that is located at the interface of S1-RBD and ACE2, the host cell receptor for SARS-CoV-2. Given its location at the interface, N501Y S1-RBD mutation can therefore potentially alter the interfacial interaction. Multiple, all-atom, explicit solvent MD simulations of the ACE2-S1-RBD complex carried here indicated a more stable interaction between the N501Y mutant S1-RBD and ACE2 through stabilizing interfacial interactions of residues at one end of the interface that are either sequentially or physically near the mutation site. These mechanistic details will aid in better understanding the mechanism by which the alpha variant has increased infectivity as well as in designing better therapeutics including ACE2-S1 spike protein inhibitors that will, in turn, help thwarting the current and future pandemic.

HighlightsO_LIN501 in the wild type SARS-CoV-2 S1-RBD forms unsustained hydrogen bonds with residues in the ACE2, namely Y41 and K353
C_LIO_LIY501 in the N501Y mutant SARS-CoV-2 S1-RBD is not capable of forming substantial hydrogen bonds with ACE2 within the time span of the current simulation
C_LIO_LIEvidence from analyzing the simulation results suggests that Y501 of S1-RBD could form other types of non-covalent interactions with ACE2, such as van der Waals interactions
C_LIO_LIN501Y S1-RBD mutation stabilizes the position of interfacial residues neighboring to the mutation site, as well as other non-interfacial residues that are distant from the mutation site
C_LIO_LIThese altered dynamics results in more stable interaction of S1-RBD with ACE2 which could be the main reason underlying the reported enhanced affinity of S1-RBD in the SARS-CoV-2 alpha variant (UK B.1.1.7 lineage) to ACE2
C_LI Modelling interference between vectors of non-persistently transmitted plant viruses to identify effective control strategies Aphids are the primary vector of plant viruses. Transient aphids, which probe several plants per day, are considered to be the principal vectors of non-persistently transmitted (NPT) viruses. However, resident aphids, which can complete their life cycle on a single host and are affected by agronomic practices, can transmit NPT viruses as well. Moreover, they can interfere both directly and indirectly with transient aphids, eventually shaping plant disease dynamics. By mean of an epidemiological model, originally accounting for ecological principles and agronomic practices, we explore the consequences of fertilization and irrigation, pesticide deployment and roguing of infected plants on the spread of viral diseases in crops. Our results indicate that the spread of NPT viruses can be i) both reduced or increased by fertilization and irrigation, depending on whether the interference is direct or indirect; ii) counter-intuitively increased by pesticide application and iii) reduced by roguing infected plants. We show that a better understanding of vectors interactions would enhance our understanding of disease transmission, supporting the development of disease management strategies.

Author summaryA range of both experimental and theoretical studies show that the behaviour and population dynamics of insects depend strongly upon interactions with other insect species. These interactions have the potential to greatly affect the dynamics of insect-vectored plant disease, as transmission of viruses is intimately dependent on the local density of vectors, as well as how they select and move between potential host plants. Surprisingly, the effects of interaction between vector species on epidemics remains little studied and even worse understood, probably because experimentation is costly and difficult. Here, we present a model which permits to investigate the effect of interaction between a virus, two vector species, and the host plant on the spread of viral disease in crops. In this study, our model is used to explore the consequences of common agronomic practices on epidemics. Our study highlights the importance of exploring vectors interactions to enhance the understanding of disease transmission, supporting the development of disease management strategies. Phylogenetic distribution and expression pattern analyses identified a divergent basal body assembly protein involved in land plant spermatogenesis Oogamy is a form of sexual reproduction and evolved independently in animals, fungi, and plants. In streptophyte plants, Charophyceae, Coleochaetophyceae, bryophytes, lycophytes, ferns (monilophytes), and some gymnosperms (Cycads and Ginkgo) utilize spermatozoids as the male gamete. Plant spermatozoids commonly possess characteristic structures such as the spline, which consists of a microtubule array, the multilayered structure (MLS) in which the uppermost layer is continuum of the spline, and multiple flagella. However, the molecular mechanisms underpinning plant spermatogenesis remain to be elucidated. To identify the genes involved in plant spermatogenesis, we performed computational analyses and successfully found deeply divergent BLD10s by combining multiple methods and omics-data. We then validated the functions of candidate genes in the liverwort Marchantia polymorpha and the moss Physcomitrium patens and found that MpBLD10 and PpBLD10 are required for normal basal body and flagella formation. Mpbld10 mutants exhibited defects in remodeling of the cytoplasm and nucleus during spermatozoid formation, thus MpBLD10 should be involved in chromatin reorganization and elimination of the cytoplasm during spermiogenesis. Streptophyte BLD10s are orthologous to BLD10/CEP135 family proteins, which function in basal body assembly, but we found that BLD10s evolved especially fast in land plants and MpBLD10 might obtain additional functions in spermatozoid formation through the fast molecular evolution. This study provides a successful example of combinatorial study from evolutionary and molecular genetic perspectives that elucidated a function of the key protein of the basal body formation that fast evolved in land plants. Prediction of survival after partial hepatectomy using a physiologically based pharmacokinetic model of indocyanine green liver function tests The evaluation of hepatic function and functional capacity of the liver are essential tasks in hepatology as well as in hepatobiliary surgery. Indocyanine green (ICG) is a widely applied test compound that is used in clinical routine to evaluate hepatic function. Important questions for the functional evaluation with ICG in the context of hepatectomy are how liver disease such as cirrhosis alters ICG elimination, and if postoperative survival can be predicted from preoperative ICG measurements. Within this work a physiologically based pharmacokinetic (PBPK) model of ICG was developed and applied to the prediction of the effects of a liver resection under various degrees of cirrhosis. For the parametrization of the computational model and validation of model predictions a database of ICG pharmacokinetic data was established. The model was applied (i) to study the effect of liver cirrhosis and liver resection on ICG pharmacokinetics; and (ii) to evaluate the model-based prediction of postoperative ICG-R15 as a measure for postoperative outcome. Key results are the accurate prediction of changes in ICG pharmacokinetics caused by liver cirrhosis and postoperative changes of ICG-elimination after liver resection, as validated with a wide range of data sets. Based on the PBPK model, individual survival after liver resection could be classified, demonstrating its potential value as a clinical tool. Sperm Histone H3 Lysine 4 tri-methylation serves as a metabolic sensor of paternal obesity and is associated with the inheritance of metabolic dysfunction ObjectiveParental environmental exposures can strongly influence descendant risks for adult disease. How paternal obesity changes the sperm chromatin leading to the acquisition of metabolic disease in offspring remains controversial and ill-defined. The objective of this study was to assess: (1) whether obesity induced by a high-fat diet alters sperm histone methylation; (2) whether paternal obesity can induce metabolic disturbances across generations; (3) whether there could be cumulative damage to the sperm epigenome leading to enhanced metabolic dysfunction in descendants; and (4) whether obesity-sensitive regions associate with embryonic epigenetic and transcriptomic profiles. Using a genetic mouse model of epigenetic inheritance, we investigated the role of histone H3 lysine 4 methylation (H3K4me3) in the paternal transmission of metabolic dysfunction. This transgenic mouse overexpresses the histone demethylase enzyme KDM1A in the developing germline and has an altered sperm epigenome at the level of histone H3K4 methylation. We hypothesized that challenging transgenic sires with a high-fat diet would further erode the sperm epigenome and lead to enhanced metabolic disturbances in the next generations.

MethodsTo assess whether paternal obesity can have inter- or transgenerational impacts, and if so, to identify potential mechanisms of this non-genetic inheritance, we used wildtype C57BL/6NCrl and transgenic males with a pre-existing altered sperm epigenome. To induce obesity, sires were fed either a control or high-fat diet (10% or 60% kcal fat, respectively) for 10-12 weeks, then bred to wildtype C57BL/6NCrl female fed a regular diet. F1 and F2 descendants were characterized for metabolic phenotypes by examining the effects of paternal obesity by sex, on body weight, fat mass distribution, the liver transcriptome, intraperitoneal glucose and insulin tolerance tests. To determine whether obesity altered the F0 sperm chromatin, native chromatin immunoprecipitation-sequencing targeting H3K4me3 was performed. To gain insight into mechanisms of paternal transmission, we compared our sperm H3K4me3 profiles with embryonic and placental chromatin states, histone modification and gene expression profiles.

ResultsObesity-induced alterations in H3K4me3 occurred at genes implicated in metabolic, inflammatory, and developmental processes. These processes were associated with offspring metabolic dysfunction and corresponded to genes enriched for H3K4me3 in embryos, and overlapped embryonic and placenta gene expression profiles. Transgenerational susceptibility to metabolic disease was only observed when obese F0 had a pre-existing modified sperm epigenome. This coincided with increased H3K4me3 alterations in sperm and more severe phenotypes affecting their offspring.

ConclusionsOur data suggest sperm H3K4me3 might serve as a metabolic sensor that connects paternal diet with offspring phenotypes via the placenta. This non-DNA based knowledge of inheritance has the potential to improve our understanding of how environment shapes heritability and may lead to novel routes for the prevention of disease. This study highlights the need to further study the connection between the sperm epigenome, placental development and childrens health.