Dynamic functional brain networks underlying the temporal inertia of negative emotions Affective inertia represents the lasting impact of transient emotions at one time point on affective state at a subsequent time point. Here we describe the neural underpinnings of inertia following negative emotions elicited by sad events in movies. Using a co-activation pattern analysis of dynamic functional connectivity, we examined the temporal expression and reciprocal interactions among brain-wide networks during movies and subsequent resting periods. Our findings revealed distinctive spatiotemporal expression of visual (VIS), default mode (DMN), central executive (CEN), and frontoparietal control (FPCN) networks both in negative movies and in rest periods following these movies. We also identified different reciprocal relationships among these networks, in transitions from movie to rest. While FPCN and DMN expression increased during and after negative movies, respectively, FPCN occurrences during the movie predicted lower DMN and higher CEN expression during subsequent rest after neutral movies, but this relationship was reversed after the elicitation of negative emotions. Changes in FPCN and DMN activity correlated with more negative subjective affect. These findings provide new insights into the transient interactions of intrinsic brain networks underpinning the inertia of negative emotions. More specifically, they describe a major role of FPCN in emotion elicitation processes, with prolonged impact on DMN activity in subsequent rest, presumably involved in emotion regulation and restoration of homeostatic balance after negative events.

HighlightsO_LIModulations of dynamic functional brain connectivity are associated to the temporal inertia of negative emotions.
C_LIO_LIFunctional co-activations patterns (CAPs) during emotional episodes predict changes in spontaneous brain dynamics during subsequent resting state.
C_LIO_LIClassical "task-rest" anticorrelations in network activity are reversed by negative emotions.
C_LI
Beyond linear regression: mapping models in cognitive neuroscience should align with research goals Many cognitive neuroscience studies use large feature sets to predict and interpret brain activity patterns. Feature sets take many forms, from human stimulus annotations to representations in deep neural networks. Of crucial importance in all these studies is the mapping model, which defines the space of possible relationships between features and neural data. Until recently, most encoding and decoding studies have used linear mapping models. Increased availability of large datasets and computing resources has allowed some researchers to employ more flexible nonlinear mapping models instead; however, the question of whether nonlinear mapping models can yield meaningful scientific insights remains debated. Here, we discuss the choice of a mapping model in the context of three overarching desiderata: predictive accuracy, interpretability, and biological plausibility. We show that, contrary to popular intuition, these desiderata do not map cleanly onto the linear/nonlinear divide; instead, each desideratum can refer to multiple research goals, each of which imposes its own constraints on the mapping model. Moreover, we argue that, instead of categorically treating the mapping models as linear or nonlinear, we should instead aim to estimate the complexity of these models. We show that, in many cases, complexity provides a more accurate reflection of restrictions imposed by various research goals. Finally, we outline several complexity metrics that can be used to effectively evaluate mapping models.
Co-infection of Candidatus Piscichlamydia trichopodus (order Chlamydiales) and Henneguya sp (Myxosporea, Myxobolidae) in snakeskin gourami Trichopodus pectoralis (Regan 1910) The present study describes a simultaneous infection of a novel Chlamydia-like organism (CLO) with a Myxozoa parasite, Henneguya sp. in snakeskin gourami Trichopodus pectoralis in Thailand. A new CLO is proposed "Candidatus Piscichlamydia trichopodus" (CPT) based on 16S rRNA phylogenetic analysis. Systemic intracellular CPT infection was confirmed by histological examination, in situ hybridization, PCR assay, and sequencing of 16S rRNA. This novel pathogen belongs to the order Chlamydiales but differs in certain aspects from other species. The histopathological changes associated with CPT infection were different from the typical pathological lesions of epitheliocystis caused by previously known CLO. Unlike other CLO, CPT localized in the connective tissue rather than in the epithelial cells and formed smaller clumps of intracellular bacteria that stained dark blue with hematoxylin. On the other hand, typical myxospores of the genus Henneguya with tails were observed in the gill sections. Infection with Henneguya sp. resulted in extensive destruction of the gill filaments, most likely leading to respiratory distress. Due to the frequency of co-infections and the unavailability of culture methods for CLO and Henneguya sp, it was difficult to determine which pathogens were directly responsible for the associated mortality. However, co-infections may increase the negative impact on the host and the severity of the disease. Given the commercial importance of the snakeskin gourami and its significant aquaculture potential, the findings of this study are important for further studies on disease prevention. Actomyosin activity-dependent apical targeting of Rab11 vesicles reinforces apical constriction During tissue morphogenesis, cell shape changes resulting from cell-generated forces often require active regulation of intracellular trafficking. How mechanical stimuli influence intracellular trafficking and how such regulation impacts tissue mechanics are not fully understood. In this study, we identify an actomyosin dependent mechanism involving Rab11- mediated trafficking in regulating apical constriction in the Drosophila embryo. During Drosophila mesoderm invagination, apical actin and Myosin II (actomyosin) contractility induces apical accumulation of Rab11-marked vesicle-like structures ("Rab11 vesicles") by promoting a directional bias in dynein mediated vesicle transport. At the apical domain, Rab11 vesicles are enriched near the adherens junctions (AJs). The apical accumulation of Rab11 vesicles is essential to prevent fragmented apical AJs, breaks in the supracellular actomyosin network and a reduction in the apical constriction rate. This Rab11 function is separate from its role in promoting apical Myosin II accumulation. These findings suggest a feedback mechanism between actomyosin activity and Rab11-mediated intracellular trafficking that regulates the force generation machinery during tissue folding. Modeling miRNA-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. Phases of cortical actomyosin dynamics coupled to the neuroblast polarity cycle The Par complex dynamically polarizes to the apical cortex of asymmetrically dividing Drosophila neuroblasts where it directs fate determinant segregation. Previously we showed that apically directed cortical movements that polarize the Par complex require F-actin (Oon and Prehoda, 2019). Here we report the discovery of cortical actomyosin dynamics that begin in interphase when the Par complex is cytoplasmic but ultimately become tightly coupled to cortical Par dynamics. Interphase cortical actomyosin dynamics are unoriented and pulsatile but rapidly become sustained and apically-directed in early mitosis when the Par protein aPKC accumulates on the cortex. Apical actomyosin flows drive the coalescence of aPKC into an apical cap that is depolarized in anaphase when the flow reverses direction. Together with the previously characterized role of anaphase flows in specifying daughter cell size asymmetry, our results indicate that multiple phases of cortical actomyosin dynamics regulate asymmetric cell division. SARS-CoV-2 Beta Variant Substitutions Alter Spike Glycoprotein Receptor Binding Domain Structure and Stability The emergence of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) and the subsequent COVID-19 pandemic has visited a terrible cost on the world in the forms of disease, death, and economic turmoil. The rapid development and deployment of extremely effective vaccines against SARS-CoV-2 have seemingly brought within reach the end of the pandemic. However, the virus has acquired mutations; and emerging variants of concern (VOC) are more infectious and reduce the efficacy of existing vaccines. While promising efforts to combat these variants are underway, the evolutionary pressures leading to these variants are poorly understood. To that end, here we have studied the effects on the structure and function of the SARS-CoV-2 spike glycoprotein receptor-binding domain of three amino-acid substitutions found in several variants of concern, including alpha (B.1.1.7), beta (B.1.351), and gamma (P.1). We found that these substitutions alter the RBD structure, stability, and ability to bind to ACE2, in such a way as to possibly have opposing and compensatory effects. These findings provide new insights into how these VOC may have been selected for infectivity while maintaining the structure and stability of the receptor binding domain. Development of an exosomal gene signature to detect residual disease in dogs with osteosarcoma using a novel xenograft platform and machine learning Osteosarcoma has a guarded prognosis. A major hurdle in developing more effective osteosarcoma therapies is the lack of disease-specific biomarkers to predict risk, prognosis, or therapeutic response. Exosomes are secreted extracellular microvesicles emerging as powerful diagnostic tools. However, their clinical application is precluded by challenges in identifying disease-associated cargo from the vastly larger background of normal exosome cargo. We developed a method using canine osteosarcoma in mouse xenografts to distinguish tumor-derived from host-response exosomal mRNAs. The model allows for the identification of canine osteosarcoma-specific gene signatures by RNA sequencing and a species-differentiating bioinformatics pipeline. An osteosarcoma-associated signature consisting of five gene transcripts (SKA2, NEU1, PAF1, PSMG2, and NOB1) was validated in dogs with spontaneous osteosarcoma by qRT-PCR, while a machine learning model assigned dogs into healthy or disease groups. Serum/plasma exosomes were isolated from 53 dogs in distinct clinical groups ("healthy", "osteosarcoma", "other bone tumor", or "non-neoplastic disease"). Pre-treatment samples from osteosarcoma cases were used as the training set and a validation set from post-treatment samples was used for testing, classifying as "osteosarcoma-detected" or "osteosarcoma-NOT detected". Dogs in a validation set whose post-treatment samples were classified as "osteosarcoma-NOT detected" had longer remissions, up to 15 months after treatment. In conclusion, we identified a gene signature predictive of molecular remissions with potential applications in the early detection and minimal residual disease settings. These results provide proof-of-concept for our discovery platform and its utilization in future studies to inform cancer risk, diagnosis, prognosis, and therapeutic response.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=147 SRC="FIGDIR/small/429432v2_ufig1.gif" ALT="Figure 1">
View larger version (41K):
org.highwire.dtl.DTLVardef@1b612c4org.highwire.dtl.DTLVardef@8a1847org.highwire.dtl.DTLVardef@3ffd64org.highwire.dtl.DTLVardef@114df8c_HPS_FORMAT_FIGEXP  M_FIG C_FIG Conformational dynamics and allosteric modulation of the SARS-CoV-2 spike Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects cells through binding to angiotensin-converting enzyme 2 (ACE2). This interaction is mediated by the receptor-binding domain (RBD) of the viral spike (S) glycoprotein. Structural and dynamic data have shown that S can adopt multiple conformations, which controls the exposure of the ACE2-binding site in the RBD. Here, using single-molecule Forster resonance energy transfer (smFRET) imaging we report the effects of ACE2 and antibody binding on the conformational dynamics of S from the Wuhan-1 strain and the B.1 variant (D614G). We find that D614G modulates the energetics of the RBD position in a manner similar to ACE2 binding. We also find that antibodies that target diverse epitopes, including those distal to the RBD, stabilize the RBD in a position competent for ACE2 binding. Parallel solution-based binding experiments using fluorescence correlation spectroscopy (FCS) indicate antibody-mediated enhancement of ACE2 binding. These findings inform on novel strategies for therapeutic antibody cocktails. HLA-SPREAD: A Natural Language Processing based resource for curating HLA association from PubMed abstracts Extreme complexity in the Human Leukocyte Antigens (HLA) system and its nomenclature makes it difficult to interpret and integrate relevant information for HLA associations with diseases, Adverse Drug Reactions (ADR) and Transplantation. PubMed search displays [~]144,000 studies on HLA reported from multiple diseases in diverse locations. Currently, IPD-IMGT/HLA database houses data on 28,320 HLA alleles. We developed an automated pipeline with a unified graphical user interface HLA-SPREAD that provides a structured information on SNPs, Populations, REsources, ADRs and Diseases information. Information on HLA was extracted from [~]24 million PubMed abstracts extracted using Natural Language Processing (NLP). Python scripts were used to mine and curate information on diseases, filter false positives and categorize to 24 tree hierarchical groups and named Entity Recognition (NER) algorithms followed by semantic analysis to infer HLA association(s). This resource from 112 countries and 32 ethnic groups provides interesting insights on: markers associated with allelic/haplotypic association in autoimmune, cancer, viral and skin diseases, transplantation outcome and ADRs for hypersensitivity. Summary information on clinically relevant biomarkers related to HLA disease associations with mapped susceptible/risk alleles are readily retrievable from HLASPREAD. The resource is available at URL http://hla-spread.igib.res.in/. This resource is first of its kind that can help uncover novel patterns in HLA gene-disease associations. Different contribution of the monkey prefrontal and premotor dorsal cortex in decision making during a transitive inference task Several studies have reported similar neural modulations between brain areas of the frontal cortex, such as the dorsolateral prefrontal (DLPFC) and the premotor dorsal (PMd) cortex, in tasks requiring encoding of the abstract rules for selecting the proper action. Here we compared the neuronal modulation of the DLPFC and PMd of monkeys trained to choose the higher rank from a pair of abstract images (target item), selected from an arbitrarily rank-ordered set (A>B>C>D>E>F) in the context of a transitive inference task. Once acquired by trial-and-error, the ordinal relationship between pairs of adjacent images (i.e., A>B; B>C; C>D; D>E; E>F), monkeys were tested in indicating the ordinal relation between items of the list not paired during learning. During these decisions, we observed that the choice accuracy increased and the reaction time decreased as the rank difference between the compared items enhanced. This result is in line with the hypothesis that after learning, the monkeys built an abstract mental representation of the ranked items, where rank comparisons correspond to the items position comparison on this representation. In both brain areas, we observed higher neuronal activity when the target item appeared in a specific location on the screen with respect to the opposite position and that this difference was particularly enhanced at lower degrees of difficulty. By comparing the time evolution of the activity of the two areas, we observed that the neural encoding of target item spatial position occurred earlier in the DLPFC than in the PMd. Evidence for spatiotemporal shift in demersal fisheries management priority areas in the western Mediterranean Marine Protected Areas (MPAs) are a promising management tool for the conservation and recovery of marine ecosystems, as well as fisheries management. MPAs are generally established as permanent closures but marine systems are dynamic, which has generated debate in favour of more dynamic designs. As a consequence, the identification of priority areas should assess their persistence in space and time. Here, we develop a step-by-step approach to assess the spatiotemporal dynamics of fisheries management priority areas using standard fishery-independent survey data. To do so, we fit Bayesian hierarchical spatiotemporal SDM models to different commercially important demersal species and use the resulting maps to fit different spatial prioritisation configurations. We use these results to assess the spatiotemporal dynamics of fisheries priority areas. The proposed method is illustrated through a western Mediterranean case study using fishery-independent trawl survey data on six commercially important species collected over 17 years. We identified two fisheries priority area patterns in the study area, each predominant during a different time-period of the study, asserting the importance of regularly re-assessing MPA designs. Liver ductal organoids reconstruct intrahepatic biliary trees in decellularized liver grafts Three-dimensional scaffolds decellularized from native organs are a promising technique to establish engineered liver grafts and overcome the current shortage of donor organs. However, limited sources of bile duct cells and inappropriate cell distribution in bioengineered liver grafts have hindered their practical application. Organoid technology is anticipated to be an excellent tool for the advancement of regenerative medicine. In the present study, we reconstructed intrahepatic bile ducts in a rat decellularized liver graft by recellularization with liver ductal organoids. Using an ex vivo perfusion culture system, we demonstrated the biliary characteristics of repopulated mouse liver organoids, which maintained bile duct markers and reconstructed biliary tree-like networks with luminal structures. We also established a method for the co-recellularization with engineered bile ducts and primary hepatocytes, revealing the appropriate cell distribution to mimic the native liver. We then utilized this model in human organoids to demonstrate the reconstructed bile ducts. Our results show that liver ductal organoids are a potential cell source for bile ducts from bioengineered liver grafts using three-dimensional scaffolds. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. A generic hidden Markov model for multi-parent populations A common step in the analysis of multi-parent populations is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of multi-parent populations. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly-available data on Diversity Outbred and Collaborative Cross mice Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Deep learning-based quantification of arbuscular mycorrhizal fungi in plant roots Soil fungi establish mutualistic interactions with the roots of most vascular land plants. Arbuscular mycorrhizal (AM) fungi are among the most extensively characterised mycobionts to date. Current approaches to quantifying the extent of root colonisation and the abundance of hyphal structures in mutant roots rely on staining and human scoring involving simple, yet repetitive tasks prone to variations between experimenters.

We developed AMFinder which allows for automatic computer vision-based identification and quantification of AM fungal colonisation and intraradical hyphal structures on ink-stained root images using convolutional neural networks.

AMFinder delivered high-confidence predictions on image datasets of roots of multiple plant hosts (Nicotiana benthamiana, Medicago truncatula, Lotus japonicus, Oryza sativa) and captured the altered colonisation in ram1-1, str and smax1 mutants. A streamlined protocol for sample preparation and imaging allowed us to quantify mycobionts from the genera Rhizophagus, Claroideoglomus, Rhizoglomus and Funneliformis via flatbed scanning or digital microscopy including dynamic increases in colonisation in whole root systems over time.

AMFinder adapts to a wide array of experimental conditions. It enables accurate, reproducible analyses of plant root systems and will support better documentation of AM fungal colonisation analyses. AMFinder can be accessed here: https://github.com/SchornacklabSLCU/amfinder.git Analyses of symbiotic bacterial communities in the plant pest Bemisia tabaci reveal high prevalence of Candidatus Hemipteriphilus asiaticus on the African continent AO_SCPLOWBSTRACTC_SCPLOWMicrobial symbionts are widespread in insects and some of them have been associated to adaptive changes. Primary symbionts (P-symbionts) have a nutritional role that allows their hosts to feed on unbalanced diets (plant sap, wood, blood). Most of them have undergone genome reduction, but their genomes still retain genes involved in pathways that are necessary to synthesize the nutrients that their hosts need. However, in some P-symbionts, essential pathways are incomplete and secondary symbionts (S-symbionts) are required to complete parts of their degenerated functions. The P-symbiont of the phloem sap-feeder Bemisia tabaci, Candidatus Portiera aleyrodidarium, lacks genes involved in the synthesis of vitamins, cofactors, and also of some essential amino-acids. Seven S-symbionts have been detected in the B. tabaci species complex. Phenotypic and genomic analyses have revealed various effects, from reproductive manipulation to fitness benefits, notably some of them have complementary metabolic capabilities to Candidatus Portiera aleyrodidarium, suggesting that their presence may be obligatory. In order to get the full picture of the symbiotic community of this pest, we investigated, through metabarcoding approaches, the symbiont content of individuals from Burkina Faso, a West African country where B. tabaci induces severe crop damage. While no new putative B. tabaci S-symbiont was identified, Candidatus Hemipteriphilus asiaticus, a symbiont only described in B. tabaci populations from Asia, was detected for the first time on this continent. Phylogenetic analyses however reveal that it is a different strain than the reference found in Asia. Specific diagnostic PCRs showed a high prevalence of these S-symbionts and especially of Candidatus Hemipteriphilus asiaticus in different genetic groups. These results suggest that Candidatus Hemipteriphilus asiaticus may affect the biology of B. tabaci and provide fitness advantage in some B. tabaci populations. A Histoplasma capsulatum lipid metabolic map identifies antifungal targets Lipids play a fundamental role in fungal cell biology, being essential cell membrane components and major targets of antifungal drugs. A deeper knowledge of lipid metabolism is key for developing new drugs and a better understanding of fungal pathogenesis. Here we built a comprehensive map of the Histoplasma capsulatum lipid metabolic pathway by incorporating proteomic and lipidomic analyses. We performed genetic complementation and overexpression of H. capsulatum genes in Saccharomyces cerevisiae to validate reactions identified in the map and to determine enzymes responsible for catalyzing orphan reactions. The map led to the identification of both the fatty acid desaturation and the sphingolipid biosynthesis pathways as targets for drug development. We found that the sphingolipid biosynthesis inhibitor myriocin, the fatty acid desaturase inhibitor thiocarlide and the fatty acid analog 10-thiastearic acid inhibit H. capsulatum growth in nanomolar to low micromolar concentrations. These compounds also reduced the intracellular infection in an alveolar macrophage cell line. Overall, this lipid metabolic map revealed pathways that can be targeted for drug development. 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. 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. Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

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

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. Tuft cell-produced cysteinyl leukotrienes and IL-25 synergistically initiate lung type 2 inflammation Aeroallergen sensing by airway epithelial cells can trigger pathogenic immune responses leading to chronic type 2 inflammation, the hallmark of airway diseases such as asthma. Airway tuft cells are specialized chemosensory epithelial cells and the dominant source of the epithelial cytokine IL-25 in the trachea and of cysteinyl leukotrienes (CysLTs) in the naive murine nasal mucosa. The interaction of IL-25 and CysLTs and the contribution of tuft cell-derived CysLTs to the development of allergen-triggered inflammation in the airways has not been clarified. Here we show that inhalation of LTC4 in combination with a subthreshold dose of IL-25 leads to dramatic synergistic induction of type 2 inflammation throughout the lungs, causing rapid eosinophilia, dendritic cell (DC) and inflammatory type 2 innate lymphoid cell (ILC2) expansion, and goblet cell metaplasia. While lung eosinophilia is dominantly mediated through the classical CysLT receptor CysLT1R, type 2 cytokines and activation of innate immune cells require signaling through both CysLT1R and CysLT2R. Tuft cell-specific deletion of the terminal enzyme requisite for CysLT production, Ltc4s, was sufficient to reduce both the innate immune response in the lung - eosinophilia, KLRG1+ ILC2 activation and DC recruitment - and the systemic immune response in the draining lymph nodes after inhalation of the mold aeroallergen Alternaria. Our findings identify surprisingly potent synergy of CysLTs and IL-25 downstream of aeroallergen-trigged activation of airway tuft cells leading to a highly polarized type 2 immune response and further implicate airway tuft cells as powerful modulators of type 2 immunity in the lungs.

One Sentence SummaryTuft cells produce two highly synergistic mediators: LTC4 and IL-25 to cooperatively induce allergen-driven airway inflammation. A phenomics approach for in vitro antiviral drug discovery BackgroundThe current COVID-19 pandemic has highlighted the need for new and fast methods to identify novel or repurposed therapeutic drugs. Here we present a method for untargeted phenotypic drug screening of virus-infected cells, combining Cell Painting with antibody-based detection of viral infection in a single assay. We designed an image analysis pipeline for segmentation and classification of virus-infected and non-infected cells, followed by extraction of morphological properties.

ResultsWe show that the methodology can successfully capture virus-induced phenotypic signatures of MRC-5 human lung fibroblasts infected with Human coronavirus 229E (CoV-229E). Moreover, we demonstrate that our method can be used in phenotypic drug screening using a panel of nine host- and virus-targeting antivirals. Treatment with effective antiviral compounds reversed the morphological profile of the host cells towards a non-infected state.

ConclusionsThe method can be used in drug discovery for morphological profiling of novel antiviral compounds on both infected and non-infected cells. N-VEGF, the forgotten isoform of VEGF-A, activate cellular stress survival circuits Vascular endothelial growth factor A (VEGF-A) is a secreted protein that stimulates angiogenesis in response to hypoxia. Under hypoxic conditions, a non-canonical long isoform called L-VEGF is concomitantly expressed with VEGF-A. Once translated, L-VEGF it is proteolytically cleaved to generate N-VEGF and VEGF-A. Interestingly, while VEGF-A is secreted and affects the surrounding cells, N-VEGF is mobilized to the nucleus. This suggests that N-VEGF participates in transcriptional response to hypoxia. In this study, we performed a series of complementary experiments to examine the functional role of N-VEGF. Strikingly, we found that the mere expression of N-VEGF followed by its hypoxia-independent mobilization to the nucleus was sufficient to induce key genes associated with angiogenesis, such as Hif1, VEGF-A isoforms, as well as genes associated with cell survival under hypoxia. Complementarily, when N-VEGF was genetically depleted, key hypoxia-induced genes were downregulated and cells were significantly susceptible to hypoxia-mediated apoptosis. This is the first reports of N-VEGF serving as an autoregulatory arm of VEGF-A. Further experiments will be needed to determine the role of N-VEGF in cancer and embryogenesis. Electrical Stimulated GLUT4 Signaling Attenuates Critical Illness-Associated Muscle Wasting BackgroundCritical illness myopathy (CIM) is a debilitating condition characterized by the preferential loss of the motor protein myosin. CIM is a byproduct of critical care, attributed to impaired recovery, long-term complications, and mortality. CIM pathophysiology is complex, heterogeneous and remains incompletely understood, however loss of mechanical stimuli contributes to critical illness associated muscle atrophy and weakness. Passive mechanical loading (ML) and electrical stimulation (ES) therapies augment muscle mass and function. While having beneficial outcomes, the mechanistic underpinning of these therapies is less known. Therefore, here we aimed to assess the mechanism by which chronic supramaximal ES ameliorates CIM in a unique experimental rat model of critical care.

MethodsRats were subjected to 8 days critical care conditions entailing deep sedation, controlled mechanical ventilation, and immobilization with and without direct soleus ES. Muscle size and function were assessed at the single cell level. RNAseq and Western blotting were employed to understand the mechanisms driving ES muscle outcomes in CIM.

ResultsFollowing 8 days of controlled mechanical ventilation and immobilization, soleus muscle mass, Myosin:Actin ratio and single muscle fiber maximum force normalized to cross-sectional area (specific force) were reduced by 40-50% (p< 0.0001). ES significantly reduced the loss of soleus muscle fiber cross-sectional area (CSA) and Myosin:Actin ratio by approximately 30% (p< 0.05) yet failed to effect specific force. RNAseq pathway analysis revealed downregulation of insulin signaling in the soleus muscle following critical care and GLUT4 trafficking was reduced by 55% leading to an 85% reduction of muscle glycogen content (p< 0.01). ES promoted phosphofructokinase and insulin signaling pathways to control levels (p< 0.05), consistent with the maintenance of GLUT4 translocation and glycogen levels. AMPK, but not AKT, signaling pathway was stimulated following ES, where the downstream target TBC1D4 increased 3 logFC (p= 0.029) and AMPK-specific P-TBC1D4 levels were increased approximately 2-fold (p= 0.06). Reduction of muscle protein degradation rather than protein synthesis promoted soleus CSA, as ES reduced E3 ubiquitin proteins, Atrogin-1 (p= 0.006) and MuRF1 (p= 0.08) by approximately 50%, downstream of AMPK-FoxO3.

ConclusionsES maintained GLUT4 translocation through increased AMPK-TBC1D4 signaling leading to improved muscle glucose homeostasis. Soleus CSA and myosin content was promoted through reduced protein degradation via AMPK-FoxO3 E3 ligases, Atrogin-1 and MuRF1. These results demonstrate chronic supramaximal ES reduces critical care associated muscle wasting, preserved glucose signaling and reduced muscle protein degradation in CIM.