Computational Simulations of the 4-D Micro-Circulatory Network in Zebrafish Tail Amputation and Regeneration Wall shear stress (WSS) in the micro-vasculature contributes to biomechanical cues that regulate mechanotransduction underlying vascular development, regeneration, and homeostasis. We hereby elucidate the interplay between hemodynamic shear forces and luminal remodeling in response to vascular injury and regeneration in the zebrafish model of tail amputation. Using the transgenic Tg(fli1:eGFP; Gata1:ds-red) line, we were able to track the enhanced green-fluorescent protein (eGFP)-labeled endothelial lining of the 3-D microvasculature for post-image segmentation and reconstruction of fluid domain for computational fluid dynamics (CFD) simulation. At 1 day post amputation (dpa), dorsal aorta (DA) and posterior cardinal vein (PCV) were severed, and vasoconstriction developed in the dorsal longitudinal anastomotic vessel (DLAV) with a concomitant increase in WSS in the segmental vessels (SV) proximal to the amputation site and a decrease in WSS in SVs distal to amputation. Simultaneously, we observed angiogenesis commencing at the tips of the amputated DLAV and PCV where WSS was minimal in the absence of blood flow. At 2 dpa, vasodilation occurred in a pair of SVs proximal to amputation, resulting in increased flow rate and WSS, whereas in the SVs distal to amputation, WSS normalized to the baseline. At 3 dpa, the flow rate in the arterial SV proximal to amputation continued to rise and merged with DLAV that formed a new loop with PCV. Thus, our CFD modeling uncovered a well-coordinated micro-vascular adaptation process following tail amputation, accompanied by the rise and fall of WSS and dynamic changes in flow rate during vascular regeneration.
A flexible electronic strain sensor for the real-time monitoring of tumor regression Assessing the efficacy of cancer therapeutics in mouse models is a critical step in treatment development. However, low resolution measurement tools and small sample sizes make determining drug efficacy in vivo a difficult and time-intensive task. Here we present a commercially scalable wearable electronic sensor that automates the in vivo testing of cancer therapeutics by continuously monitoring micrometer-scale tumor volume progression or regression in real-time. Histology, caliper measurements, and bioluminescence imaging over a one-week treatment period validated the sensors recordings, but the sensors higher time and length scale resolutions provided the only measurement technique capable of continuously monitoring the immediate hour-scale pharmacodynamic response of a given drug. In mice with two subcutaneously implanted tumor models our sensors recorded significant volume reductions in tumors just 5 hours after small molecule or immunotherapy treatment initiation. We anticipate that real-time tumor regression datasets could help expedite and automate the process of screening cancer therapies in vivo.
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. Feed efficiency of lactating Holstein cows was not as repeatable across diets as within diet over subsequent lactation stages BackgroundImproving feed efficiency has become a common target for dairy farmers to meet the requirement of producing more milk with fewer resources. To improve feed efficiency, a prerequisite is to ensure that the cows identified as most or least efficient will remain as such, independently of diet composition. Therefore, the current research analysed the ability of lactating dairy cows to maintain their feed efficiency while changing the energy density of the diet by changing its concentration in starch and fibre. A total of 60 lactating Holstein cows, including 33 primiparous cows, were first fed a high starch diet-low fibre (diet S+F-), then switched over to a low starch diet-high fibre (diet S-F+). To know if diet affect feed efficiency, we compared the ability of feed efficiency to be maintained within a diet over subsequent lactation stages, known as repeatability of feed efficiency, with its ability to be maintained across diets, known as reproducibility of feed efficiency. To do so we used two indicators: the estimation of the error of repeatability/reproducibility, which is commonly used in metrology, and the coefficient of correlation of concordance (CCC), which is used in biology. The effect of diet change could also lead to a change in cows sorting behaviour which could potentially affect feed efficiency if for example the most efficient cows select more concentrate than the least efficient. We therefore analysed the relationship between the differences in individual feed refusals composition and the differences in feed efficiency. To do so, the composition of each feed refusal was described with its near infra-red (NIR) spectroscopy and was performed on each individual feed ingredient, diet and refusals and used as composition variable. The variability of the NIR spectra of the refusals was described with its principal components thanks to a principal component analysis (PCA). The Pearson correlation was estimated to check the relationship between feed efficiency and refusals composition, i.e. sorting behaviour.

ResultsThe error of reproducibility of feed efficiency across diets was 2.95 MJ/d. This error was significantly larger than the errors of repeatability estimated within diet, which were 2.01 MJ/d within diet S-F+ and 2.40 MJ/d within diet S+F-. The CCC was 0.64 between feed efficiency estimated within diet S+F- and feed efficiency estimated within diet S-F+. This CCC was smaller than the one observed for feed efficiency estimated within diet between two subsequent lactation stages (CCC = 0.72 within diet S+F- and 0.85 within diet S-F+). Feed efficiency was poorly correlated to the first two principal components, which explained 90% of the total variability of the NIR spectra of the individual refusals. This suggests that feed sorting behaviour did not explain differences in feed efficiency.

ConclusionsFeed efficiency was significantly less reproducible across diets than repeatable within the same diet over subsequent lactation stages, but cows ranking for feed efficiency was not significantly affected by diet change. This loss in repeatability across diets could be due to a more pronounced feed sorting subsequent to the change in diet composition. However, the differences in sorting behaviour between cows were not associated to feed efficiency differences in this trial neither with the S+F- diet nor with the S-F+ diet. Those results have to be confirmed on diets having different forage to concentrate ratios to ensure that the least and most efficient cows will not change. No evidence of human genome integration of SARS-CoV-2 found by long-read DNA sequencing A recent study proposed severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) hijacks the LINE-1 (L1) retrotransposition machinery to integrate into the DNA of infected cells. If confirmed, this finding could have significant clinical implications. Here, we applied deep (>50x) long-read Oxford Nanopore Technologies (ONT) sequencing to HEK293T cells infected with SARS-CoV-2, and did not find the virus integrated into the genome. By examining ONT data from separate HEK293T cultivars, we completely resolved 78 L1 insertions arising in vitro in the absence of L1 overexpression systems. ONT sequencing applied to hepatitis B virus (HBV) positive liver cancer tissues located a single HBV insertion. These experiments demonstrate reliable resolution of retrotransposon and exogenous virus insertions via ONT sequencing. That we found no evidence of SARS-CoV-2 integration suggests such events are, at most, extremely rare in vivo, and therefore are unlikely to drive oncogenesis or explain post-recovery detection of the virus. Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector. Induction of exaggerated cytokine production in human peripheral blood mononuclear cells by a recombinant SARS-CoV-2 spike glycoprotein S1 and its inhibition by dexamethasone An understanding of the pathological inflammatory mechanisms involved in SARS-CoV-2 virus infection is necessary in order to discover new molecular pharmacological targets for SARS-CoV-2 cytokine storm. In this study, the effects of a recombinant SARS-CoV-2 spike glycoprotein S1 was investigated in human peripheral blood mononuclear cells (PBMCs). Stimulation of PBMCs with spike glycoprotein S1 (100 ng/mL) resulted in significant elevation in the production of TNF, IL-6, IL-1{beta} and IL-8. However, pre-treatment with dexamethasone (100 nM) caused significant reduction in the release of these cytokines. Further experiments revealed that S1 stimulation of PBMCs increased phosphorylation of NF-{kappa}B p65 and I{kappa}B, and I{kappa}B degradation. DNA binding of NF-{kappa}B p65 was also significantly increased following stimulation with spike glycoprotein S1. Treatment of PBMCs with dexamethasone (100 nM) or BAY11-7082 (1 {micro}M) resulted in inhibition of spike glycoprotein S1-induced NF-{kappa}B activation. Activation of p38 MAPK by S1 was blocked in the presence of dexamethasone and SKF 86002. CRID3, but not dexamethasone pre-treatment produced significant inhibition of S1-induced activation of NLRP3/caspase-1. Further experiments revealed that S1-induced increase in the production of TNF, IL-6, IL-1{beta} and IL-8 was reduced in the presence of BAY11-7082 and SKF 86002, while CRID3 pre-treatment resulted in the reduction of IL-1{beta} production. These results suggest that SARS-CoV-2 spike glycoprotein S1 stimulated PBMCs to release pro-inflammatory cytokines through mechanisms involving activation of NF-{kappa}B, p38 MAPK and NLRP3 inflammasome. It is proposed that the clinical benefits of dexamethasone in COVID-19 is possibly due to its anti-inflammatory activity in reducing SARS-CoV-2 cytokine storm. 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. Proteasome complexes experience profound structural and functional rearrangements throughout mammalian spermatogenesis During spermatogenesis, spermatogonia undergo a series of mitotic and meiotic divisions on their path to spermatozoa. To achieve this, a succession of processes requiring high proteolytic activity are in part orchestrated by the proteasome. The spermatoproteasome (s20S) is specific to the developing gametes, in which the gamete-specific 4s subunit replaces the 4 isoform found in the constitutive proteasome (c20S). Although the s20S is conserved across species and was shown to be crucial for germ cell development, its mechanism, function and structure remain incompletely characterized. Here, we used advanced mass spectrometry (MS) methods to map the composition of proteasome complexes and their interactomes throughout spermatogenesis. We observed that the s20S becomes highly activated as germ cells enter meiosis, mainly through a particularly extensive 19S activation and, to a lesser extent, PA200 binding. Additionally, the proteasome population shifts from c20S (98%) to s20S (>82-92%) during differentiation, presumably due to the shift from 4 to 4s expression. We demonstrated that s20S, but not c20S, interacts with components of the meiotic synaptonemal complex, where it may localize via association with the PI31 adaptor protein. In vitro, s20S preferentially binds to 19S, and displays higher trypsin- and chymotrypsin-like activities, both with and without PA200 activation. Moreover, using MS methods to monitor protein dynamics, we identified significant differences in domain flexibility between 4 and 4s. We propose that these differences induced by 4s incorporation result in significant changes in the way the s20S interacts with its partners, and dictate its role in germ cell differentiation. Testing for phylogenetic signal in single-cell RNA-seq data AO_SCPLOWBSTRACTC_SCPLOWPhylogenetic methods are emerging as a useful tool to understand cancer evolutionary dynamics, including tumor structure, heterogeneity, and progression. Most currently used approaches utilize either bulk whole genome sequencing (WGS) or single-cell DNA sequencing (scDNA-seq) and are based on calling copy number alterations and single nucleotide variants (SNVs). Here we explore the potential of single-cell RNA sequencing (scRNA-seq) to reconstruct cancer evolutionary dynamics. scRNA-seq is commonly applied to explore differential gene expression of cancer cells throughout tumor progression. The method exacerbates the single-cell sequencing problem of low yield per cell with uneven expression levels. This accounts for low and uneven sequencing coverage and makes SNV detection and phylogenetic analysis challenging. In this paper, we demonstrate for the first time that scRNA-seq data contains sufficient evolutionary signal and can be utilized in phylogenetic analyses. We explore and compare results of such analyses based on both expression levels and SNVs called from scRNA-seq data. Both techniques are shown to be useful for reconstructing phylogenetic relationships between cells, reflecting the clonal composition of a tumor. Both standardized expression values and SNVs appear to be equally capable of reconstructing a similar pattern of phylogenetic relationship. This pattern is stable even when phylogenetic uncertainty is taken in account. Our results open up a new direction of somatic phylogenetics based on scRNA-seq data. Further research is required to refine and improve these approaches to capture the full picture of somatic evolutionary dynamics in cancer. Control of bacterial immune signaling by a WYL domain transcription factor Bacteria use diverse immune systems to defend themselves from ubiquitous viruses termed bacteriophages (phages). Many anti-phage systems function by abortive infection to kill a phage-infected cell, raising the question of how they are regulated to avoid cell killing outside the context of infection. Here, we identify a transcription factor associated with the widespread CBASS bacterial immune system, that we term CapW. CapW forms a homodimer and binds a palindromic DNA sequence in the CBASS promoter region. Two crystal structures of CapW suggest that the protein switches from an unliganded, DNA binding-competent state to a ligand-bound state unable to bind DNA. We show that CapW strongly represses CBASS gene expression in uninfected cells, and that phage infection causes increased CBASS expression in a CapW-dependent manner. Unexpectedly, this CapW-dependent increase in CBASS expression is not required for robust anti-phage activity, suggesting that CapW may mediate CBASS activation and cell death in response to a signal other than phage infection. Our results parallel concurrent reports on the structure and activity of BrxR, a transcription factor associated with the BREX anti-phage system, suggesting that CapW and BrxR are members of a family of universal defense signaling proteins. Erythropoietin supplementation ameliorates the immunological and hematological deficiencies of lysinuric protein intolerance in mice SLC7A7 encodes for y+LAT1, a cationic amino acid transporter expressed on the basolateral membrane of epithelial cells and macrophages. Mutations in human SLC7A7 cause lysinuric protein intolerance (LPI), a rare autosomal disorder with a variety of clinical symptoms. Citrulline and a low protein diet are widely used to treat the main symptoms of LPI; however, the hematological and immunological abnormalities in LPI remain poorly understood and difficult to manage by established treatments. Here, we show that erythropoietin levels are significantly diminished in mice lacking Slc7a7, which reduces the number of erythroblast precursors, resulting in mature red blood cells with a smaller mean corpuscular volume. Ex vivo analysis demonstrated that erythrocytes from Slc7a7-deficient mice were phagocytosed more by bone marrow-derived macrophages. Gene expression analysis of sorted red pulp macrophages (RPMs) from Slc7a7-deficient mice indicated that this cell population had higher erythrophagocytic activity, but preserved heme catabolic activity. The increase in erythrophagocytosis resulted in iron accumulation in macrophages and tissues, and ultimately led to RPM apoptosis. Surprisingly, ablation of Slc7a7 in myeloid cells failed to trigger tissue iron accumulation or changes in the RPM population, indicating that the expression of this gene in macrophages is not key for the development of iron accumulation in tissue. Finally, restoration of serum erythropoietin levels in Slc7a7-deficient mice re-established normal erythropoiesis and erythrophagocytosis and prevented the accumulation of iron in tissues. Thus, erythropoietin has a fundamental role in the immune-hematological alterations of LPI in mice, revealing a new crucial role for this hormone in this disease. Comparative Analysis of T Cell Spatial Proteomics and the Influence of HIV Expression As systems biology approaches to virology have become more tractable, highly studied viruses such as HIV can now be analyzed in new, unbiased ways, including spatial proteomics. We employed here a differential centrifugation protocol to fractionate Jurkat T cells for proteomic analysis by mass spectrometry; these cells contain inducible HIV-1 genomes, enabling us to look for changes in the spatial proteome induced by viral gene expression. Using these proteomics data, we evaluated the merits of several reported machine learning pipelines for classification of the spatial proteome and identification of protein translocations. From these analyses we found that classifier performance in this system was organelle-dependent, with Bayesian t-augmented Gaussian mixture modeling outperforming support vector machine (SVM) learning for mitochondrial and ER proteins, but underperforming on cytosolic, nuclear, and plasma membrane proteins by QSep analysis. We also observed a generally higher performance for protein translocation identification using a Bayesian model, BANDLE, on SVM-classified data. Comparative BANDLE analysis of cells induced to express the wild-type viral genome vs. cells induced to express a genome unable to express the accessory protein Nef identified known Nef-dependent interactors such as TCR signaling components and coatomer complex. Lastly, we found that SVM classification showed higher consistency and was less sensitive to HIV-dependent noise. These findings illustrate important considerations for studies of the spatial proteome following viral infection or viral gene expression and provide a reference for future studies of HIV-gene-dropout viruses. oCEM: Automatic detection and analysis of overlapping co-expressed gene modules BackgroundWhen it comes to the co-expressed gene module detection, its typical challenges consist of overlap between identified modules and local co-expression in a subset of biological samples. The nature of module detection is the use of unsupervised clustering approaches and algorithms. Those methods are advanced undoubtedly, but the selection of a certain clustering method for sample- and gene-clustering tasks is separate, in which the latter task is often more complicated.

ResultsThis study presented an R-package, Overlapping CoExpressed gene Module (oCEM), armed with the decomposition methods to solve the challenges above. We also developed a novel auxiliary statistical approach to select the optimal number of principal components using a permutation procedure. We showed that oCEM outperformed state-of-the-art techniques in the ability to detect biologically relevant modules additionally.

ConclusionsoCEM helped non-technical users easily perform complicated statistical analyses and then gain robust results. oCEM and its applications, along with example data, were freely provided at https://github.com/huynguyen250896/oCEM. Daily torpor reduces the energetic consequences of habitat selection for a widespread bat Homeothermy requires increased metabolic rates as temperatures decline below the thermoneutral zone, so homeotherms typically select microhabitats within or near their thermoneutral zones during periods of inactivity. However, many mammals and birds are heterotherms that relax internal controls on body temperature when maintaining a high, stable body temperature is energetically costly. Such heterotherms should be less tied to microhabitats near their thermoneutral zones, and because heterotherms spend more time in torpor and expend less energy at colder temperatures, heterotherms may even select microhabitats in which temperatures are well below their thermoneutral zones. We studied how temperature and daily torpor influence selection of diurnal roosts by a heterothermic bat (Myotis thysanodes). We (1) quantified the relationship between ambient temperature and daily duration of torpor, (2) simulated daily energy expenditure over a range of microhabitat (roost) temperatures, and (3) quantified the influence of roost temperature on roost selection. While warm roosts substantially reduced energy expenditure of simulated homeothermic bats, heterothermic bats modulated their use of torpor to maintain a constant level of energy expenditure over the course of a day. Daily torpor expanded the range of energetically economical microhabitats, such that roost selection was independent of roost temperature. Our work adds to a growing literature documenting functions of torpor beyond its historical conceptualization as a last-resort measure to save energy during extended or acute energetic stress. Abscisic Acid Drives Stomatal Closure through Increases in Hydrogen Peroxide in Distinct Subcellular Compartments Including Mitochondria Abscisic acid (ABA) closes stomata to minimize transpiration and gas exchange in response to drought. ABA signaling increases accumulation of reactive oxygen species (ROS) in guard cells to drive stomatal closure. We used chemical dyes and genetically-encoded biosensors to identify the subcellular localizations where ABA increases hydrogen peroxide (H2O2) in guard cells. The general redox sensor, dichlorofluorescein (DCF), revealed ABA-driven ROS increases in distinct subcellular locations including chloroplasts, cytosol, and nuclei. ABA also increased fluorescence of both DCF and an H2O2-selective probe, Peroxy Orange 1, in mitochondria, which colocalize with Mitotracker Red. Whole cell and mitochondrial localized ABA-increased ROS signals were lost in an ABA-insensitive quintuple ABA receptor mutant and accentuated in an ABA-hypersensitive mutant. The ABA-induced ROS accumulation in mitochondria was lost in guard cells of mutants with defects in genes encoding H2O2 elevating respiratory burst oxidase homolog (RBOH) enzymes and reduced by treatment with the RBOH inhibitor VAS2870. ABA treatment of guard cells transformed with the genetically-encoded H2O2 reporter targeted to the cytoplasm, roGFP2-Orp1, and the mitochondrial targeted, mt-roGFP2-Orp1, revealed significant ABA-induced H2O2 increases. Furthermore, H2O2 treatment is sufficient to close stomata. These results reveal that ABA elevates H2O2 accumulation in guard cell cytosol and mitochondria to drive stomatal closure.

One sentence summaryGenetically encoded biosensors and chemical probes revealed ABA-dependent increases in hydrogen peroxide, a reactive oxygen species with signaling activity, in guard cell cytoplasm and mitochondria. Multispecific DARPin(R) therapeutics demonstrate very high potency against SARS-CoV-2 variants in vitro SARS-CoV-2 has infected millions of people globally and continues to undergo evolution. Emerging variants can be partially resistant to vaccine induced immunity and therapeutic antibodies, emphasizing the urgent need for accessible, broad-spectrum therapeutics. Here, we report a comprehensive study of ensovibep, the first trispecific clinical DARPin candidate, that can simultaneously engage all three units of the spike protein trimer to potently inhibit ACE2 interaction, as revealed by structural analyses. The cooperative binding of the individual modules enables ensovibep to retain inhibitory potency against all frequent SARS-CoV-2 variants, including Omicron BA.1 and BA.2, as of February 2022. Moreover, viral passaging experiments show that ensovibep, when used as a single agent, can prevent development of escape mutations comparably to a cocktail of monoclonal antibodies (mAb). Finally, we demonstrate that the very high in vitro antiviral potency also translates into significant therapeutic protection and reduction of pathogenesis in Roborovski dwarf hamsters infected with either the SARS-CoV-2 wild-type or the Alpha variant. In this model, ensovibep prevents fatality and provides substantial protection equivalent to the standard of care mAb cocktail. These results support further clinical evaluation and indicate that ensovibep could be a valuable alternative to mAb cocktails and other treatments for COVID-19. Determinants of Murine Skin Microbiota Composition in Homeostasis and Wound Healing Animal microbiota have complex interactions with hosts and environment that determines its composition. Yet the ability of hosts to determine their microbiota composition is less well studied. In this study, to investigate the role host genetics in determining skin microbiota, we used 30 different mouse strains from the recombinant inbred panel, the Collaborative Cross. Murine skin microbiota composition was strongly dependent on murine strain with > 50% of the variation explained by murine strain. In particular, a quantitative trait locus on chromosome 4 associates both with Staphylococcus abundance and principal-component multi-trait analyses. Additionally, excisional wound associated changes in microbiota composition were not uniform across mouse strains and were host-specific, the genetic background accounting for about 40% of the variation in microbiota. Genetic background also had the highest effect on the healing speed of wounds accounting for over 50% of the variation while mouse age and microbiota composition change accounted only for 20% and 5% of the healing speed despite reaching statistical significance. In conclusion, host genetics has a significant impact on the skin microbiota composition during both homeostasis and wound healing. These findings have long reaching implications in our understanding of associations between microbiota dysbiosis and disease. SkewDB: A comprehensive database of GC and 10 other skews for over 28,000 chromosomes and plasmids GC skew denotes the relative excess of G nucleotides over C nucleotides on the leading versus the lagging replication strand of eubacteria. While the effect is small, typically around 2.5%, it is robust and pervasive. GC skew and the analogous TA skew are a localized deviation from Chargaffs second parity rule, which states that G and C, and T and A occur with (mostly) equal frequency even within a strand.

Most bacteria also show the analogous TA skew. Different phyla show different kinds of skew and differing relations between TA and GC skew.

This article introduces an open access database (https://skewdb.org) of GC and 10 other skews for over 28,000 chromosomes and plasmids. Further details like codon bias, strand bias, strand lengths and taxonomic data are also included. The SkewDB database can be used to generate or verify hypotheses. Since the origins of both the second parity rule, as well as GC skew itself, are not yet satisfactorily explained, such a database may enhance our understanding of microbial DNA. Darwin's small and medium ground finches might have taste preferences, but not for human foods Urbanization is rapidly changing ecological niches. On the inhabited Galapagos Islands, Darwins finches consume human-introduced foods preferentially; however, it remains unclear why. Here we presented pastry with flavour profiles typical of human foods (oily, salty, sweet) to small ground finches (Geospiza fuliginosa) and medium ground finches (Geospiza fortis) to test if latent taste preferences might drive selection of human foods. If human-food flavours were consumed more than a neutral or bitter control only at sites with human foods, then we predicted tastes were acquired after urbanization; however, if no site-differences were found then this would indicate latent taste preferences. Contrary to both predictions, we found little evidence that human-food flavours were preferred compared to control flavours at any site. Instead, finches showed a weak aversion to oily foods, but only at remote (no human foods present) sites. This was further supported by behavioural responses, with beak wiping occurring more often at remote sites after finches tasted flavours associated with human foods. Our results suggest, therefore, that while Darwins finches may have acquired a tolerance to the flavours of human food, latent taste preferences are unlikely to have played a major role in their dietary response to increased urbanization. Ecologically pre-trained RNNs explain suboptimal animal decisions The strategies found by animals facing a new task are determined both by individual experience and by structural priors evolved to leverage the statistics of natural environments. Rats can quickly learn to capitalize on the trial sequence correlations of two-alternative forced choice (2AFC) tasks after correct trials, but consistently deviate from optimal behavior after error trials, when they waive the accumulated evidence. To understand this outcome-dependent gating, we first show that Recurrent Neural Networks (RNNs) trained in the same 2AFC task outperform rats as they can readily learn to use across-trial information both after correct and error trials. We hypothesize that, while RNNs can optimize their behavior in the 2AFC task without any a priori restrictions, rats strategy is constrained by a structural prior adapted to a natural environment in which rewarded and non-rewarded actions provide largely asymmetric information. When pre-training RNNs in a more ecological task with more than two possible choices, networks develop a strategy by which they gate off the across-trial evidence after errors, mimicking rats behavior. Population analyses show that the pre-trained networks form an accurate representation of the sequence statistics independently of the outcome in the previous trial. After error trials, gating is implemented by a change in the network dynamics which temporarily decouples the categorization of the stimulus from the across-trial accumulated evidence. Our results suggest that the suboptimal behavior observed in rats reflects the influence of a structural prior that reacts to errors by isolating the network decision dynamics from the context, ultimately constraining the performance in a 2AFC laboratory task. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. Nuclear Progestin Receptor Mediated Linkage of Blood Coagulation and Ovulation Ovulation is a dramatic remodeling process that includes rupture of blood capillaries and clotting, but coagulation is not thought to directly regulate this process. Herein, we report remarkable increases of coagulation factors V (f5, [~]3145-fold) and tissue factor (f3a, [~]120-fold) in zebrafish ovarian follicle cells during ovulation. This increase was mediated through the nuclear progestin receptor (Pgr), which is essential for ovulation in zebrafish, and was totally abolished in ovarian follicular cells from pgr-/- mutants. In addition, promoter activities of f5 and f3a were significantly enhanced by progestin (DHP) via Pgr. Similar regulation of human F5 promoter activity was induced via human PGRB, suggesting a conserved mechanism. Site-directed mutagenesis of the zebrafish f5 promoter further demonstrated a direct regulation of coagulation factors via progestin response elements. Moreover, a stark increase of erythrocytes occurred in capillaries meshed in wildtype preovulatory follicles but was absent in pgr-/- mutants. Interestingly, anticoagulants significantly inhibited ovulation both in vitro and in vivo, respectively. Furthermore, reduced fecundity was observed in f5+/- female zebrafish. Taken together, our study provides plausible evidence for steroid regulation of coagulation factors, and a new hypothesis for blood clotting triggered ovulation in vertebrates. Modelling homing suppression gene drive in haplodiploid organisms Gene drives have shown great promise for suppression of pest populations. These engineered alleles can function by a variety of mechanisms, but the most common is the CRISPR homing drive, which converts wild-type alleles to drive alleles in the germline of heterozygotes. Some potential target species are haplodiploid, in which males develop from unfertilized eggs and thus have only one copy of each chromosome. This prevents drive conversion, a substantial disadvantage compared to diploids where drive conversion can take place in both sexes. Here, we study homing suppression gene drives in haplodiploids and find that a drive targeting a female fertility gene could still be successful. However, such drives are less powerful than in diploids and suffer more from functional resistance alleles. They are substantially more vulnerable to high resistance allele formation in the embryo due to maternally deposited Cas9 and gRNA and also to somatic cleavage activity. Examining spatial models where organisms move over a continuous landscape, we find that haplodiploid suppression drives surprisingly perform nearly as well as in diploids, possibly due to their ability to spread further before inducing strong suppression. Together, these results indicate that gene drive can potentially be used to effectively suppress haplodiploid populations. A Neurogenic Signature Involving Monoamine Oxidase-A controls Human Thermogenic Adipose Tissue Development Mechanisms that control "beige/brite" thermogenic adipose tissue development may be harnessed to improve human metabolic health. To define these mechanisms, we developed a species-hybrid model in which human mesenchymal progenitor cells were used to develop white or thermogenic/beige adipose tissue in mice. The hybrid adipose tissue developed distinctive features of human adipose tissue, such as larger adipocyte size, despite its neurovascular architecture being entirely of murine origin. Thermogenic adipose tissue recruited a denser, qualitatively distinct vascular network, differing in genes mapping to circadian rhythm pathways, and denser sympathetic innervation. The enhanced thermogenic neurovascular network was associated with human adipocyte expression of THBS4, TNC, NTRK3 and SPARCL1, which enhance neurogenesis, and decreased expression of MAOA and ACHE, which control neurotransmitter tone. Systemic inhibition of MAOA, which is present in human but absent in mouse adipocytes, induced browning of human but not mouse adipose tissue, revealing the physiological relevance of this pathway. Our results reveal species-specific cell type dependencies controlling the development of thermogenic adipose tissue and point to human adipocyte MAOA as a potential target for metabolic disease therapy. 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