Temporal Scaling of Human Scalp-Recorded Potentials Much of human behaviour is governed by common processes that unfold over varying timescales. Standard event-related potential analysis assumes fixed-latency responses relative to experimental events. However, recent single unit recordings in animals have revealed neural activity scales to span different durations during behaviours demanding flexible timing. Here, we employed a general linear modelling approach using a novel combination of fixed-duration and variable-duration regressors to unmix fixed-time and scaled-time components in human magneto/electroencephalography (M/EEG) data. We use this to reveal consistent temporal scaling of human scalp-recorded potentials across four independent EEG datasets, including interval perception, production, prediction and value-based decision making. Between-trial variation in the temporally scaled response predicts between-trial variation in subject reaction times, demonstrating the relevance of this temporally scaled signal for temporal variation in behaviour. Our results provide a general approach for studying flexibly timed behaviour in the human brain.
Age-related differences in negative cognitive empathy but similarities in positive affective empathy Empathy, among other social-cognitive processes, changes across adulthood. More specifically, cognitive components of empathy (understanding anothers perspective) appear to decline with age, while findings for affective empathy (sharing anothers emotional state) are rather mixed. Structural and functional correlates underlying cognitive and affective empathy in aging and the extent to which valence affects empathic response in brain and behavior are not well understood yet. To fill these research gaps, younger and older adults completed a modified version of the Multifaceted Empathy Test, which measures both cognitive and affective empathy as well as empathic responding to both positive and negative stimuli (i.e., positive vs. negative empathy). Adopting a multimodal imaging approach and applying multivariate analysis, the study found that for cognitive empathy to negative emotions, regions of the salience network including the anterior insula and anterior cingulate were more involved in older than younger adults. For affective empathy to positive emotions, in contrast, younger and older adults recruited a similar brain network including main nodes of the default mode network. Additionally, increased structural integrity (fractional anisotropy values) of the posterior, but not the anterior, cingulum bundle was related to activation of default mode regions during affective empathy for positive emotions in both age groups. These findings provide novel insights into the functional networks subserving cognitive and affective empathy in younger and older adults and highlight the importance of considering valence in empathic response in aging research. Further this study, for the first time, underscores the role of the posterior cingulum bundle in higher-order social-cognitive processes such as empathy, specifically for positive emotions, in aging.
Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. A unified mechanism for innate and learned visual landmark guidance in the insect central complex Insects can navigate efficiently in both novel and familiar environments, and this requires flexiblity in how they are guided by sensory cues. A prominent landmark, for example, can ellicit strong innate behaviours (attraction or menotaxis) but can also be used, after learning, as a specific directional cue as part of a navigation memory. However, the mechanisms that allow both pathways to co-exist, interact or override each other are largely unknown. Here we propose a model for the behavioural integration of innate and learned guidance based on the neuroanatomy of the central complex (CX), adapted to control landmark guided behaviours. We consider a reward signal provided either by an innate attraction to landmarks or a long-term visual memory in the mushroom bodies (MB) that modulates the formation of a local vector memory in the CX. Using an operant strategy for a simulated agent exploring a simple world containing a single visual cue, we show how the generated short-term memory can support both innate and learned steering behaviour. In addition, we show how this architecture is consistent with the observed effects of unilateral MB lesions in ants that cause a reversion to innate behaviour. We suggest the formation of a directional memory in the CX can be interpreted as transforming rewarding (positive or negative) sensory signals into a mapping of the environment that describes the geometrical attractiveness (or repulsion). We discuss how this scheme might represent an ideal way to combine multisensory information gathered during the exploration of an environment and support optimized cue integration.

1 Author summaryIn this paper, we modeled the neural pathway allowing insects to perform landmark guided behaviours using their internal compass. First, we observed the intrinsic property of the connectome, extracted from drosophila online database, between the internal compass neurons and the steering neurons to support an oriented behaviour towards a single landmark. Then, we proposed and evaluated an adaptation of the bees path integration neural circuit, to sustain flexible landmark guidance behaviours such as attraction or menotaxis. We showed the model ability to form a memory during the exploration of the local environment to support both innate or learned navigation behaviour using a single landmark in the environment. In addtion, we demonstrated the transformation of a simple goodness/badness signal, from innate or long-term memory pathways, into an oriented steering signal that could be applied to other sensory pathways. Furthermore, by reproducing lesion experiments in the mushroom bodies of wood ants we highlight the consistency of the model with biological observation. We then discuss the different emergent properties and the potential outcome that this local, and operant, memory supports. Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. An Extended Motif in the SARS-CoV-2 Spike Modulates Binding and Release of Host Coatomer in Retrograde Trafficking {beta}-Coronaviruses such as SARS-CoV-2 hijack coatomer protein-I (COPI) for spike protein retrograde trafficking to the progeny assembly site in endoplasmic reticulum-Golgi intermediate compartment (ERGIC). However, limited residue-level details are available into how the spike interacts with COPI. Here we identify an extended COPI binding motif in the spike that encompasses the canonical K-x-H dibasic sequence. This motif demonstrates selectivity for COPI subunit. Guided by an in silico analysis of dibasic motifs in the human proteome, we employ mutagenesis and binding assays to show that the spike motif terminal residues are critical modulators of complex dissociation, which is essential for spike release in ERGIC. COPI residues critical for spike motif binding are elucidated by mutagenesis and crystallography and found to be conserved in the zoonotic reservoirs, bats, pangolins, camels, and in humans. Collectively, our investigation on the spike motif identifies key COPI binding determinants with implications for retrograde trafficking. Fatal neuroinvasion and SARS-CoV-2 tropism in K18-hACE2 mice is partially independent on hACE2 expression Animal models recapitulating distinctive features of severe COVID-19 are critical to enhance our understanding of SARS-CoV-2 pathogenesis. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. The precise mechanisms of lethality in this mouse model remain unclear. Here, we evaluated the spatiotemporal dynamics of SARS-CoV-2 infection for up to 14 days post-infection. Despite infection and moderate pneumonia, rapid clinical decline or death of mice was invariably associated with viral neuroinvasion and direct neuronal injury (including brain and spinal neurons). Neuroinvasion was observed as early as 4 dpi, with virus initially restricted to the olfactory bulb supporting axonal transport via the olfactory neuroepithelium as the earliest portal of entry. No evidence of viremia was detected suggesting neuroinvasion occurs independently of entry across the blood brain barrier. SARS-CoV-2 tropism was not restricted to ACE2-expressing cells (e.g., AT1 pneumocytes), and some ACE2-positive lineages were not associated with the presence of viral antigen (e.g., bronchiolar epithelium and brain capillaries). Detectable ACE2 expression was not observed in neurons, supporting overexpression of ACE2 in the nasal passages and neuroepithelium as more likely determinants of neuroinvasion in the K18-hACE2 model. Although our work incites caution in the utility of the K18-hACE2 model to study global aspects of SARS-CoV-2 pathogenesis, it underscores this model as a unique platform for exploring the mechanisms of SARS-CoV-2 neuropathogenesis that may have clinical relevance acknowledging the growing body of evidence that suggests COVID-19 may result in long-standing neurologic consequences.

IMPORTANCECOVID-19 is predominantly a respiratory disease caused by SARS-CoV-2 that has infected more than 191 million people with over 4 million fatalities (2021-07-20). The development of animal models recapitulating distinctive features of severe COVID-19 is critical to enhancing our understanding of SARS-CoV-2 pathogenesis and in the evaluation of vaccine and therapeutic efficacy. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. Here, we show lethality of this model is invariably associated with viral neuroinvasion linked with viral replication and assembly. Importantly, pneumonia albeit invariably present was generally moderate with the absence of culturable infectious virus at peak neuroinvasion. The dynamics of viral neuroinvasion and pneumonia were only partially dependent on hACE2. Overall, this study provides an in-depth sequential characterization of the K18-hACE2 model following SARS-CoV-2 infection, highlighting its significance to further study the mechanisms of SARS-CoV-2 neuropathogenesis. Mobile element integration reveals a chromosome dimer resolution system in Legionellales. In bacteria, the mechanisms used to repair DNA lesions during genome replication include homologous recombination between sister chromosomes. This can lead to the formation of chromosome dimers if an odd number of crossover events occurs. The concatenated DNA must be resolved before cell separation to ensure genomic stability and cell viability. The broadly conserved dif/Xer system counteracts the formation of dimers by catalyzing one additional crossover event immediately prior to cell separation. While dif/Xer systems have been characterized or predicted in the vast majority of proteobacteria, no homologs to dif or xer have been identified in the order Legionellales. Here we report the discovery of a distinct single-recombinase dif/Xer system in the intracellular pathogen Legionella pneumophila. The dif site was uncovered by our analysis of Legionella mobile element-1 (LME-1), which harbors a dif site-mimic and integrates into the L. pneumophila genome via site-specific recombination. We demonstrate that lpg1867 (herein named xerL) encodes a tyrosine recombinase that is necessary and sufficient for catalyzing recombination at the dif site, and that deletion of dif or xerL causes filamentation along with extracellular and intracellular growth defects. We show that the dif/XerL system is present throughout Legionellales and that Coxiella burnetii XerL and its cognate dif site can functionally substitute for the native system in L. pneumophila. Lastly, we describe an unexpected link between C. burnetii dif/Xer and the maintenance of its virulence plasmids.

Significance StatementThe maintenance of circular chromosomes depends on the ability to resolve aberrant chromosome dimers as they form. In most proteobacteria, broadly conserved Xer recombinases catalyze single crossovers at short, species-specific dif sites located near the replication terminus. Chromosomal dimerization leads to the formation of two copies of dif, leading to rapid site-specific recombination and elimination of the duplicated intervening sequence. The apparent absence of chromosome-dimer resolution mechanisms in Legionellales has been a mystery to date. By studying a phage-like mobile genetic element, LME-1, we have identified a previously unknown single-recombinase dif/Xer system that is not only widespread across Legionellales but whose activity is linked to virulence in two important human pathogens. Standardized genome-wide function prediction enables comparative functional genomics: a new application area for Gene Ontologies in plants BackgroundGenome-wide gene function annotations are useful for hypothesis generation and for prioritizing candidate genes potentially responsible for phenotypes of interest. We functionally annotated the genes of 18 crop plant genomes across 14 species using the GOMAP pipeline.

ResultsBy comparison to existing GO annotation datasets, GOMAP-generated datasets cover more genes, contain more GO terms, and produce datasets similar in quality (based on precision and recall metrics using existing gold standards as the basis for comparison). From there, we sought to determine whether the datasets across multiple species could be used together to carry out comparative functional genomics analyses in plants. To test the idea and as a proof of concept, we created dendrograms of functional relatedness based on terms assigned for all 18 genomes. These dendrograms were compared to well-established species-level evolutionary phylogenies to determine whether trees derived were in agreement with known evolutionary relationships, which they largely are. Where discrepancies were observed, we determined branch support based on jack-knifing then removed individual annotation sets by genome to identify the annotation sets causing unexpected relationships.

ConclusionsGOMAP-derived functional annotations used together across multiple species generally retain sufficient biological signal to recover known phylogenetic relationships based on genome-wide functional similarities, indicating that comparative functional genomics across species based on GO data hold promise for generating novel hypotheses about comparative gene function and traits. Cis-Regulatory Element Hijacking by Structural Variants Overshadows Higher-Order Topological Changes in Primary Prostate Cancer Prostate cancer is a heterogeneous disease whose progression is linked to genome instability. However the impact of this instability on the three-dimensional chromatin organization and how this drives progression is unclear. Using primary benign and tumour tissue, we find a high concordance in the higher-order three-dimensional genome organization across normal and prostate cancer cells. This concordance argues for constraints to the topology of prostate tumour genomes. Nonetheless, we identify changes to focal chromatin interactions and show how structural variants can induce these changes to guide cis-regulatory element hijacking. Such events result in opposing differential expression on genes found at antipodes of rearrangements. Collectively, our results argue that cis-regulatory element hijacking from structural variant-induced altered focal chromatin interactions overshadows higher-order topological changes in the development of primary prostate cancer. Bioactive constituents of Verbena officinalis regulate killing efficiency of primary human natural killer cells by accelerating killing processes BackgroundNatural killer (NK) cells play a key role in eliminating tumorigenic and pathogen-infected cells. Verbena officinalis (V. officinalis) has been used as a medical plant in traditional and modern medicine, exhibiting anti-tumor and anti-inflammation activity, but its role in NK cells remains unclear.

PurposeThe impact of bioactive constituents of V. officinalis on immune responses remains largely elusive. In this work we investigated the potential targets of V. officinalis and focused on killing efficiency and related functions of NK cells regulated by bioactive constituents of V. officinalis.

Study design/methodsWe used primary human NK cells from peripheral blood mononuclear cells. Potential regulatory roles of selected compounds in NK killing function-related genes were first analyzed by network pharmacology approaches. Potential engagement between compounds and NK receptors was examined by molecular docking assay. Killing efficiency of NK cells was determined with real-time killing assay and live-cell imaging. Proliferation was examined by CFSE staining. Expression of cytotoxic proteins was analyzed using flow cytometry. Lytic granule release was quantified by CD107a degranulation assay. Contact time required for killing and determination of serial killers were analyzed using live cell imaging results.

ResultsResults from network pharmacology approaches and molecular docking analyses predicted that the five bioactive compounds (acteoside, apigenin, kaempferol, verbenalin and hastatoside) from V. officinalis can regulated various key genes essential for NK effector functions and can bind to NK inhibitory receptors as well as activating receptors, indicating their potential role in regulating NK killing capacity. Analyses of NK killing efficiency revealed that out of these five compounds verbenalin exhibited an encouraging potency to enhance NK killing efficiency. Further investigation showed that verbenalin did not affect proliferation, expression of cytotoxic proteins, or lytic granule degranulation, but rather reduced contact time required for killing, enhanced total killing events per NK cell, and elevated the fraction of killing competent NK cells.

ConclusionsIn this work, we have identified verbenalin as the most potent bioactive constituent from V. officinalis to enhance NK cell killing efficiency via accelerating killing processes and converting non-killers to killing competent NK cells without affecting the lytic granule pathway. The SLIM1 transcription factor is required for arsenic resistance in Arabidopsis thaliana The transcriptional regulators of arsenic-induced gene expression remain largely unknown. Sulfur assimilation is tightly linked with arsenic detoxification. Here we report that mutant alleles in the SLIM1 transcription factor are substantially more sensitive to arsenic than cadmium. Arsenic treatment caused high levels of oxidative stress in the slim1 mutants, and slim1 alleles were impaired in both thiol and sulfate accumulation. We further found enhanced arsenic accumulation in roots of slim1 mutants. Transcriptome analyses indicate an important role for SLIM1 in arsenic-induced tolerance mechanisms. The present study identifies the SLIM1 transcription factor as an essential component in arsenic tolerance and arsenic-induced gene expression. Our results suggest that the severe arsenic sensitivity of the slim1 mutants is caused by altered redox status. 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. Bone marrow haematopoietic stem cells influence liver homeostatic networks and cancer development In fatty liver disease, systemic homeostasis is perturbed. While pre-clinical models are used to understand its pathogenesis, translating this knowledge to patients is difficult. However, by focusing on the most preserved homeostasis systems between species and models, novel disease dimensions can be unearthed. We interrogated core liver gene co-expression networks in a mouse model of liver cancer following dietary challenge. Immunometabolic modules showed temporal changes under the influence of diet duration and aging. The behaviour of immune network in tumours mirrored their counterparts in non-tumour liver. A high immune response network was associated with a lower tumour burden in mice and humans. This module in mice was enriched for genes related to haematopoietic cell differentiation. Consistently, the bone marrow haematopoietic stem and progenitor cells response was reflective of the liver immune response. Linking haematopoiesis to hepatic homeostasis uncovers a hitherto unexplored dimension of tissue crosstalk that can inform pathogenesis. Critical re-evaluation of the slope factor of the operational model of agonism Background and PurposeAlthough being a relative term, agonist efficacy is a cornerstone in the proper assessment of agonist selectivity and signalling bias. The operational model of agonism (OMA) has become successful in the determination of agonist efficacies and ranking them. In 1985, Black et al. introduced the slope factor to the OMA to make it more flexible and allow for fitting steep as well as flat concentration-response curves.

Experimental ApproachFunctional analysis of OMA demonstrates that the slope factor implemented by Black et al. affects relationships among parameters of the OMA, leading to wrong estimates of operational efficacy and affinity of an agonist. In contrast, functional analysis of OMA with Hill coefficient as a slope factor has not this flaw.

Key ResultsFitting of the OMA with Black et al. slope factor to concentration-response curves of several models (homotropic allosteric modulation, substrate inhibition and non-competitive auto-inhibition) gives wrong estimates of operational efficacy and affinity. In contrast, fitting of the OMA modified by the Hill coefficient to concentration-response curves of these several models gives correct estimates of operational efficacy and affinity.

Conclusion and ImplicationsThe OMA modified by the Hill coefficient can be used for correct ranking of agonism and subsequent analysis, like quantification of signalling bias. Phenotypic Heterogeneity Facilitates Survival While Hindering the Evolution of Drug Resistance Due to Intraspecific Competition Rising rates of resistance to antimicrobial drugs threatens the effective treatment of infections across the globe. Drug resistance has been established to emerge from non-genetic mechanisms, such as "persistence" in quiescent microbes and fluctuations in gene expression in actively replicating cells, as well as from genetic mutations. However, it is still unclear how non-genetic drug resistance affects the evolution of genetic drug resistance. We develop deterministic and stochastic population models that incorporate resource competition to quantitatively investigate the transition from non-genetic to genetic resistance during the exposure to static and cidal drugs. We find that non-genetic resistance facilitates the survival of cell populations during drug treatment, but that it hinders the development of genetic resistance due to the competition between the non-genetically and genetically resistant subpopulations. Non-genetic drug resistance in the presence of subpopulation competition is found to increase the first-appearance and fixation times of drug resistance mutations, while increasing the probability of mutation before population extinction during cidal drug treatment. Intense intraspecific competition during drug treatment leads to extinction of the susceptible and non-genetically resistant subpopulations. These findings advance our fundamental understanding of the evolution of drug resistance and may guide novel treatment strategies for patients with drug-resistant infections.

SIGNIFICANCEDrug resistance is predicted to kill as many as 10 million people per year and cost over 100 trillion USD in cumulative lost production globally by 2050. To mitigate these socio-economic costs, we need to fundamentally understand the drug resistance process. We investigate the effect that different forms of resistance have on the evolution of drug resistance using mathematical modeling and computer simulations. We find that the presence of non-genetically drug-resistant cells (whose resistance is temporary and not encoded in a genetic mutation) allows the population to survive drug treatment, while competition between these subpopoulations simultaneously slows down the evolution of permanent genetic drug resistance and in some cases drives them extinct. These findings have important implications for advancing evolutionary theory and for developing effective "resistance-proof" treatments. Production and characterization of virus-free, CRISPR-CAR T cells capable of inducing solid tumor regression BackgroundChimeric antigen receptor (CAR) T cells traditionally harbor viral vector-based sequences that encode the CAR transgene in the genome. These T cell products have yet to show consistent anti-tumor activity in patients with solid tumors. Further, viral vector manufacturing is resource intensive, suffers from batch-to-batch variability, and includes several animal components, adding regulatory and supply chain pressures.

MethodsAnti-GD2 CAR T cells were generated using CRISPR/Cas9 within nine days using recombinant Cas9 protein and nucleic acids, without any viral vectors or animal components. The CAR was specifically targeted to the T Cell Receptor Alpha Constant gene (TRAC). T cell products were characterized at the level of the genome, transcriptome, proteome, and secretome using CHANGE-seq, scRNA-seq, spectral cytometry, and ELISA assays. Functionality was evaluated in vivo in an NSG xenograft neuroblastoma model.

ResultsIn comparison to traditional retroviral CAR T cells, virus-free CRISPR CAR (VFC-CAR) T cells exhibit TRAC-targeted genomic integration of the CAR transgene, elevation of transcriptional and protein characteristics associated with a memory phenotype, and low tonic signaling prior to infusion arising in part from the the knockout of the TCR. Upon exposure to the GD2 target antigen, anti-GD2 VFC-CAR T cells exhibited specific cytotoxicity against GD2+ cells in vitro and induced solid tumor regression in vivo, with robust homing, persistence, and low exhaustion against a human neuroblastoma xenograft model.

ConclusionsThis proof-of-principle study leveraging virus-free genome editing technology could enable flexible manufacturing of clinically relevant, high-quality CAR T cells to treat cancers, including solid tumors. The cell growth cycle: resource allocation and subcellular pattern formation in budding yeast. Metabolic oscillation percolate throughout cellular physiology. The transcriptome oscillates between gene cohorts encoding for biosynthesis and growth, and for catabolism and stress response. However, long protein half-lives dampen effects on protein abundance, and the function of periodic transcription remains elusive. To test three prevailing hypotheses, we performed RNA-seq analysis during a dynamic state of the system, just preceding a bifurcation of system dynamics to a longer period. A few enzymes with short protein half-live and high transcript abundance amplitude could play a direct role in metabolic dynamics (hypothesis H1), e.g., carbonic anhydrase, sulfate uptake and the glyoxylate cycle. Substantial changes in the relative duration of expression of the antagonistic cohorts precede the bifurcation, supporting the idea of a function in cellular resource allocation (H2). The pulse-width modulation model, a mathematical formulation of this idea, explains the dependence of the cycle period on the culture growth rate. This pulse-like model of cell growth provides a first theoretical framework, where the phenomenon is understood as a mechanism of cellular resource allocation and protein homeostasis. The universal temporal program of transcription encodes the spatial structure of the cell, and we further suggest a function in pattern formation of the growing cell (H3). Beta-cell specific Insr deletion promotes insulin hypersecretion and improves glucose tolerance prior to global insulin resistance Abstract

Insulin receptor (Insr) protein can be found at higher levels in pancreatic {beta}-cells than in most other tissues, but the consequences of {beta}-cell insulin resistance remain enigmatic. Ins1cre allele was used to delete Insr specifically in {beta}-cells of both female and male mice. Experimental mice were compared to Ins1cre-containing littermate controls at multiple ages and on multiple diets. RNA-seq of purified recombined {beta}-cells revealed transcriptomic consequences of Insr loss, which differed between female and male mice. Action potential and calcium oscillation frequencies were increased in Insr knockout {beta}- cells from female, but not male mice, whereas only male {beta}InsrKO mice had reduced ATP-coupled oxygen consumption rate and reduced expression of genes involved in ATP synthesis. Female {beta}InsrKO and {beta}InsrHET mice exhibited elevated insulin release in perifusion experiments, during hyperglycemic clamps, and following i.p. glucose challenge. Deletion of Insr did not alter {beta}-cell area up to 9 months of age, nor did it impair hyperglycemia-induced proliferation. Based on our data, we adapted a mathematical model to include {beta}-cell insulin resistance, which predicted that {beta}-cell Insr knockout would improve glucose tolerance depending on the degree of whole-body insulin resistance. Indeed, glucose tolerance was significantly improved in female {beta}InsrKO and {beta}InsrHET mice when compared to controls at 9, 21 and 39 weeks, and also in insulin-sensitive 4-week old males. We did not observe improved glucose tolerance in older male mice or in high fat diet-fed mice, corroborating the prediction that global insulin resistance obscures the effects of {beta}-cell specific insulin resistance. The propensity for hyperinsulinemia was associated with mildly reduced fasting glucose and increased body weight. We further validated our main in vivo findings using the Ins1-CreERT transgenic line and found that male mice had improved glucose tolerance 4 weeks after tamoxifen-mediated Insr deletion. Collectively, our data show that loss of {beta}-cell Insr contributes to glucose-induced hyperinsulinemia, thereby improving glucose homeostasis in otherwise insulin sensitive dietary and age contexts. An orthologous gene coevolution network provides insight into eukaryotic cellular and genomic structure and function Orthologous gene coevolution--which refers to gene pairs whose evolutionary rates covary across speciation events--is often observed among functionally related genes. We present a comprehensive gene coevolution network inferred from the examination of nearly three million orthologous gene pairs from 332 budding yeast species spanning [~]400 million years of eukaryotic evolution. Modules within the network provide insight into cellular and genomic structure and function, such as genes functioning in distinct cellular compartments and DNA replication. Examination of the phenotypic impact of network perturbation across 14 environmental conditions using deletion mutant data from the bakers yeast Saccharomyces cerevisiae suggests that fitness in diverse environments is impacted by orthologous gene neighborhood and connectivity. By mapping the network onto the chromosomes of S. cerevisiae and the opportunistic human pathogen Candida albicans, which diverged [~]235 million years ago, we discovered that coevolving orthologous genes are not clustered in either species; rather, they are most often located on different chromosomes or far apart on the same chromosome. The budding yeast coevolution network captures the hierarchy of eukaryotic cellular structure and function, provides a roadmap for genotype-to-phenotype discovery, and portrays the genome as an extensively linked ensemble of genes. 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. 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. Hepatic MIR20B promotes nonalcoholic fatty liver disease by suppressing PPARA BackgroundNon-alcoholic fatty liver disease (NAFLD) is characterized by excessive lipid accumulation and imbalances in lipid metabolism in the liver. Although nuclear receptors (NRs) play a crucial role in hepatic lipid metabolism, the underlying mechanisms of NR regulation in NAFLD remain largely unclear.

MethodsUsing network analysis and RNA-seq to determine the correlation between NRs and microRNA in human NAFLD patients, we revealed that MIR20B specifically targets PPARA. MIR20B mimic and anti-MIR20B were administered to human HepG2 and Huh-7 cells and mouse primary hepatocytes as well as high fat diet (HFD)- or methionine-deficient diet (MCD)-fed mice to verify the specific function of MIR20B in NAFLD. We tested the inhibition of the therapeutic effect of a PPAR agonist, fenofibrate, by MIR20B and the synergic effect of combination of fenofibrate with anti-MIR20B.

ResultsWe revealed that MIR20B specifically targets PPARA through miRNA regulatory network analysis of nuclear receptor genes in NAFLD. The expression of MIR20B was upregulated in free fatty acid (FA)- treated hepatocytes and the livers of both obesity-induced mice and NAFLD patients. Overexpression of MIR20B significantly increased hepatic lipid accumulation and triglyceride levels. Furthermore, MIR20B significantly reduced FA oxidation and mitochondrial biogenesis by targeting PPARA. In MIR20B-introduced mice, the effect of fenofibrate to ameliorate hepatic steatosis was significantly suppressed. Finally, inhibition of MIR20B significantly increased FA oxidation and uptake, resulting in improved insulin sensitivity and a decrease in NAFLD progression. Moreover, combination of fenofibrate and anti-MIR20B exhibited the synergic effect on improvement of NAFLD in MCD-fed mice.

ConclusionsTaken together, our results demonstrate that the novel MIR20B targets PPARA, plays a significant role in hepatic lipid metabolism, and present an opportunity for the development of novel therapeutics for NAFLD.

FundingThis research was funded by Korea Mouse Phenotyping Project (2016M3A9D5A01952411), the National Research Foundation of Korea (NRF) grant funded by the Korea government (2020R1F1A1061267, 2018R1A5A1024340, NRF-2021R1I1A2041463, 2020R1I1A1A01074940), and the Future-leading Project Research Fund (1.210034.01) of UNIST. Major transcriptomic, epigenetic and metabolic changes underly the pluripotency continuum in rabbit preimplantation embryos Despite the growing interest in the rabbit model for developmental and stem cell biology, the characterization of embryos at the molecular level is still poorly documented. We conducted a transcriptome analysis of rabbit pre-implantation embryos from E2.7 (morula stage) to E6.6 (early primitive streak stage) using bulk and single-cell RNA-sequencing. In parallel, we studied oxidative phosphorylation and glycolysis and analysed active and repressive epigenetic modifications during blastocyst formation and expansion. We generated a transcriptomic, epigenetic, and metabolic map of the pluripotency continuum in rabbit preimplantation embryos and identified novel markers of naive pluripotency that might be instrumental for deriving naive pluripotent stem cell lines. Although the rabbit is evolutionarily closer to mice than to primates, we found that the transcriptome of rabbit epiblast cells shares common features with that of humans and non-human primates.

Summary StatementRabbit preimplantation embryos share characteristics with human and monkey embryos with respect to timing of early lineage segregation and expression of marker genes for naive and primed pluripotency. A network approach to genetic circuit designs As genetic circuits become more sophisticated, the size and complexity of data about their designs increases. This data captured goes beyond monolithic genetic sequences and towards circuit modularity and functional details, which are beneficial for analyzing circuit performance and establishing design automation techniques. However, the accessibility, visualisation and usability of design data (and metadata) have received relatively little attention to date. Here, we present a method to turn circuit designs into networks and showcase its potential to enhance the utility of design data. Since networks are dynamic structures, initial graphs can be interactively shaped into sub-networks of relevant information based on requirements such as abstraction, hierarchy and protein interactions. Additionally, several visual changes can be applied, such as colouring or clustering nodes based on types (e.g., genes or promoters), resulting in easier comprehension from a user perspective. This approach allows circuit designs to be coupled to other networks, such as metabolic pathways or implementation protocols captured in graph-like formats. Therefore, we advocate using networks to structure, access and improve synthetic biology information. No impact of biocontrol agent's predation cues on development time or size of surviving Aedes albopictus under optimal nutritional availability Cyclopoid copepods have been applied successfully to limit populations of highly invasive Aedes albopictus mosquitoes that can transmit diseases of public health importance. However, there is concern that changes in certain mosquito traits, induced by exposure to copepod predation, might increase the risk of disease transmission. In this study, third instar Ae. albopictus larvae, hereafter referred to as "focal individuals", were exposed to Megacyclops viridis predator cues associated with both consumption of newly-hatched mosquito larvae and attacks on focal individuals. The number of newly-hatched larvae surrounding each focal larva was held constant to control for density effects on size, and the focal individuals day of pupation and wing length were recorded for each replicate.

Exposing late instar Ae. albopictus to predation decreased their chances of surviving to adulthood, and three focal larvae that died in the predator treatment showed signs of melanization, indicative of wounding. Among surviving focal Ae. albopictus, no significant difference in either pupation day or wing length was observed due to copepod predation. The absence of a significant sublethal impact from M. viridis copepod predation on surviving later-stage larvae in this analysis supports the use of M. viridis as a biocontrol agent against Ae. albopictus.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=106 SRC="FIGDIR/small/460697v2_ufig1.gif" ALT="Figure 1">
View larger version (9K):
org.highwire.dtl.DTLVardef@1abaabborg.highwire.dtl.DTLVardef@1561651org.highwire.dtl.DTLVardef@d325a7org.highwire.dtl.DTLVardef@4d1830_HPS_FORMAT_FIGEXP  M_FIG C_FIG