High resolution, serial imaging of early mouse and human liver bud morphogenesis in three dimensions BackgroundLiver organogenesis has thus far served as a paradigm for solid organ formation. The developing liver bud is a well established model of organogenesis, and murine genetic studies demonstrate key molecules involved in key morphogenetic changes. However, the analysis of the liver bud is typically limited to 2D tissue sections, which precludes extensive visualization, quantitation, and analysis. Further, the lack of human liver bud data has further hindered our understanding of human liver organogenesis. Therefore, new analytical and visualization approaches are needed to elicit further morphogenetic details of liver organogenesis and to elucidate differences between mouse and human liver bud growth.

ResultsTo address this need, we focused on high resolution imaging, visualization, and analysis of early liver growth by using available online databases for both mouse (EMAP, Toronto Phenogenomics center) and human (3D Atlas of Human Embryology), noninvasive multimodality imaging studies of the murine embryo, and mouse/human liver weight data. First, we performed three-dimensional (3D reconstructions) of stacked, digital tissue sections that had been initially segmented for the liver epithelium and the septum transversum mesenchyme (STM). 3D reconstruction of both mouse and human data sets enabled visualization and analysis of the dynamics of liver bud morphogenesis, including hepatic cord formation and remodeling, mechanisms of growth, and liver-epithelial STM interactions. These studies demonstrated potentially under-appreciated mechanisms of growth, including rapid exponential growth that is matched at the earliest stages by STM growth, and unique differences between mouse and human liver bud growth. To gain further insight into the exponential liver bud growth that was observed, we plotted volumetric data from 3D reconstruction together with fetal liver growth data from multimodality (optical projection tomography, magnetic resonance imaging, micro-CT) and liver weight data to compose complete growth curves during mouse (E8.5-E18) and human (day 25-300) liver development. For further analysis, we performed curve fitting and parameter estimation, using Gompertzian models, which enables the comparison between mouse and human liver bud growth, as well as comparisons to processes like liver regeneration. To demonstrate the importance of mesenchyme in rapid liver bud growth and morphogenesis in the human liver bud, we performed functional analysis in which human pluripotent stem cell (hPSC)-derived hepatic organoids were used to model collective migration that occurs in vivo, demonstrating that migration is strongly dependent upon mesenchyme.

DiscussionOur data demonstrates improved visualization with 3D images, under-appreciated and potentially new mechanisms of growth, complete liver growth curves with quantitative analysis through embryonic and fetal stages, and a new functional human stem cell-derived liver organoid assay demonstrating mesenchyme-driven collective migration. These data enhance our understanding of liver organogenesis.
Loss of Resf1 reduces the efficiency of embryonic stem cell self-renewal and germline entry Retroelement silencing factor 1 (RESF1) interacts with the key regulators of mouse embryonic stem cells (ESCs) OCT4 and NANOG, and its absence results in sterility of mice. However, the function of RESF1 in ESCs and germ line specification is poorly understood. In this study, we used Resf1 knockout cell lines to determine the requirements of RESF1 for ESCs self-renewal and for in vitro specification of ESCs into primordial germ cell-like cells (PGCLCs). We found that deletion of Resf1 in ESCs cultured in serum and LIF reduces self-renewal potential whereas episomal expression of RESF1 has a modest positive effect on ESC self-renewal. In addition, RESF1 is not required for the capacity of NANOG and its downstream target ESRRB to drive self-renewal in the absence of LIF. However, Resf1 deletion reduces efficiency of PGCLC differentiation in vitro. These results identify Resf1 as a novel player in the regulation of pluripotent stem cells and germ cell specification.
ZIKV disrupts placental ultrastructure and drug transporter expression in mice Congenital Zika virus (ZIKV) infection can induce fetal brain abnormalities. Here, we investigated whether maternal ZIKV infection affects placental physiology and metabolic transport potential and impacts the fetal outcome, regardless of viral presence in the fetus at term. Low (103 PFU-ZIKVPE243; low ZIKV) and high (5x107 PFU-ZIKVPE243; high ZIKV) virus titers were injected into immunocompetent (ICompetent C57BL/6) and immunocompromised (ICompromised A129) mice at gestational day (GD) 12.5 for tissue collection at GD18.5 (term). High ZIKV elicited fetal death rates of 66% and 100%, whereas low ZIKV induced fetal death rates of 0% and 60% in C57BL/6 and A129 dams, respectively. All surviving fetuses exhibited intrauterine growth restriction (IUGR) and decreased placental efficiency. High-ZIKV infection in C57BL/6 and A129 mice resulted in virus detection in maternal spleens and placenta, but only A129 fetuses presented virus RNA in the brain. Nevertheless, pregnancies in both strains produced fetuses with decreased head sizes (p<0.05). Low-ZIKV-A129 dams had higher IL-6 and CXCL1 levels (p<0.05), and their placentas showed increased CCL-2 and CXCL-1 contents (p<0.05). In contrast, low-ZIKV-C57BL/6 dams had an elevated CCL2 serum level and increased type I and II IFN expression in the placenta. Notably, less abundant microvilli and mitochondrial degeneration were evidenced in the placental labyrinth zone (Lz) of ICompromised and high-ZIKV-ICompetent mice but not in low-ZIKV-C57BL/6 mice. In addition, decreased placental expression of the drug transporters P-glycoprotein (P-gp) and breast cancer resistance protein (Bcrp) and the lipid transporter Abca1 was detected in all ZIKV-infected groups, but Bcrp and Abca1 were only reduced in ICompromised and high-ZIKV ICompetent mice. Our data indicate that gestational ZIKV infection triggers specific proinflammatory responses and affects placental turnover and transporter expression in a manner dependent on virus concentration and maternal immune status. Placental damage may impair proper fetal-maternal exchange function and fetal growth/survival, likely contributing to congenital Zika syndrome. Staphylococcus aureus-serine protease-like protein B (SplB) activates PAR2 and induces endothelial barrier dysfunction Staphylococcus aureus (S. aureus) is a major cause of life-threatening systemic infection in humans. To cause blood stream infections such as sepsis and endocarditis, the bacteria must overcome the hosts endothelial barrier. The serine protease-like proteins (Spls) of S. aureus are known to contribute to pneumonia and allergic airway inflammation in animal models, but their role in endothelial damage is unknown. Here we demonstrate that SplB induces proinflammatory cytokine release in primary human vascular endothelial cells (HUVECs) in vitro. Mechanistically, we show that SplB selectively cleaves and activates human proteinase-activated receptor-2 (PAR2), and induces biased signaling via G13, Gi1/oB (NF-k{beta}), and {beta}-arrestin-1/-2. This activation did not trigger Gq/11-mediated calcium release nor ERK phosphorylation. Inhibition of PAR2 in HUVECs reduced the SplB-mediated cytokine release. Intravital microscopy of cremaster muscles in mice demonstrated that administration of SplB causes microvascular leakage. Genetic deletion of PAR2 in mice or neutralization of SplB with a monoclonal antibody preserved the endothelial barrier. This study identifies PAR2 as a receptor and substrate for SplB and highlights its role in mediating endothelial damage.



O_FIG O_LINKSMALLFIG WIDTH=199 HEIGHT=200 SRC="FIGDIR/small/424670v2_ufig1.gif" ALT="Figure 1">
View larger version (32K):
org.highwire.dtl.DTLVardef@16a60a2org.highwire.dtl.DTLVardef@ee4e41org.highwire.dtl.DTLVardef@153cd7borg.highwire.dtl.DTLVardef@d746eb_HPS_FORMAT_FIGEXP  M_FIG Graphical abstract C_FIG Mate sharing in male stump-tailed macaques as a possible case of coalition-like behavior to modify the group-wise fitness distribution In multi-male multi-female groups of animals, male reproductive success is often skewed toward top-ranking males. Coalition formation by middle- to low-rankers can be seen as a collaborative effort to modify the distribution of reproductive success within the group, so that they gain more than they would do otherwise. It has been pointed out, on the other hand, that a coalition of top-ranking males could hardly be profitable in the sense that they would gain little additional benefit from making collaborative effort. Here we report our novel observation of facultative sharing of mating opportunities among males in a wild population of stump-tailed macaques (Macaca arctoides) as a possible case of coalition-like behavior in which dominant males jointly guard females from mating with subordinate males and actively share mating opportunities within the allies. First, we report our novel observation of facultative sharing of mating opportunities in male stump-tailed macaques, where two or more males remain in close proximity to and copulate with a female in turn without contesting or sneaking. Second, considering the kind of coalition formation in which dominant males collaboratively exclude subordinates from mating competition and thereby strengthen the reproductive skew that already exists, we specify, by means of mathematical modeling, the condition for this kind of coalition formation to be selectively favored. Finally, we derive predictions about the occurrence of the coalition-like behavior depending on ecological and demographic factors, and test them empirically using data from the five groups of stump-tailed macaques in our study population. Re-convolving the compositional landscape of primary and recurrent glioblastoma reveals prognostic andtargetable tissue states Glioblastoma is an aggressive diffusely infiltrating neoplasm that spreads beyond surgical resection margins, where it intermingles with non-neoplastic brain cells. This complex microenvironment harbouring infiltrating glioma and non-neoplastic brain cells is the origin of tumor recurrence. Thus, understanding the cellular and molecular features of the glioma microenvironment is therapeutically and prognostically important. Here, we used single-nucleus RNA sequencing (snRNAseq) of primary and recurrent glioma to define three compositional  tissue-states rooted in cohabitation of cell-types and transcriptional states. These comprise a state featuring A) abundant normal brain microenvironment cells, B) reactive/inflamed brain microenvironment, and C) cellular/proliferative tumor. All these states exhibited variable degrees of infiltration by glioma cells. Spatial transcriptomics confirmed that the cell-types and transcriptomics states which compositionally cohabitate tissue states also colocalize in space. Tissue states are clinically significant because they correlate with radiographic, histopathologic, and prognostic features. We found that enrichment of tissue state B signature correlated with shorter survival. Importantly, we found that our compositionally defined tissue states are enriched in distinct metabolic pathways. One such pathway is fatty acid biosynthesis, which was enriched in tissue state B - a state enriched in recurrent glioblastoma and composed of astrocyte-like/mesenchymal glioma cells, reactive astrocytes resembling those seen in neurodegeneration, and monocyte-like myeloid cells. We showed that treating acute slices of GBM with a fatty acid synthesis inhibitor is sufficient to deplete the transcriptional signature of tissue state B. Our findings define a novel compositional approach to the glioma-infiltrated tissue which allows us to discover prognostic and targetable features, paving the way to new mechanistic and therapeutic discoveries. 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. Joint representations of color and form in mouse visual cortex described by random pooling from rods and cones Spatial transitions in color can aid any visual perception task, and its neural representation - the "integration of color and form" - is thought to begin at primary visual cortex (V1). An integration of color and form is untested in mouse V1, yet studies show that the ventral retina provides the necessary substrate from green-sensitive rods and UV-sensitive cones. Here, we used two-photon imaging in V1 to measure spatial frequency (SF) tuning along four axes of rod and cone contrast space, including luminance and color. We first reveal that V1s sensitivity to color is similar to luminance, yet average SF tuning is significantly shifted lowpass for color. Next, guided by linear models, we used SF tuning along all four color axes to estimate the proportion of neurons that fall into classic models of color opponency - "single-", "double-", and "non-opponent". Few neurons (~6%) fit criteria for double-opponency, which are uniquely tuned for chromatic borders. Most of the population can be described as a unimodal distribution ranging from strongly single-opponent to non-opponent. Consistent with recent studies of the rodent and primate retina, our V1 data is well-described by a simple model in which ON and OFF channels to V1 sample the photoreceptor mosaic randomly. A nontoxic tetrodotoxin analog, 5,6,11-trideoxytetrodotoxin, as an olfactory chemoattractant for the grass puffer Toxic puffers contain the potent neurotoxin, tetrodotoxin (TTX). Although TTX is considered to serve as a defense substance, previous behavioral studies have demonstrated that TTX acts as an attractive pheromone for some toxic puffers. To elucidate the physiological mechanism of putative pheromonal action of TTX, we examined whether grass puffers Takifugu alboplumbeus can detect TTX. Electroolfactogram (EOG) results suggest that the olfactory epithelium (OE) of grass puffers responded to a type of TTX analog (5,6,11-trideoxyTTX), although it did not respond to TTX. We also examined the attractive action of 5,6,11-trideoxyTTX on grass puffers by recording their swimming behavior under dark conditions. Grass puffers preferred to stay on the side of the aquarium where 5,6,11-trideoxyTTX was administered, and their swimming speed decreased. Additionally, odorant-induced labeling of olfactory sensory neurons by immunohistochemistry against neural activity marker (phosphorylated extracellular signal regulated kinase; pERK) revealed that labeled olfactory sensory neurons were localized in the region surrounding "islets" where there was considered as nonsensory epithelium. 5,6,11-trideoxyTTX has been known to accumulate in grass puffers, but its toxicity is much lower (almost nontoxic) than TTX. Our results suggest that toxic puffers may positively use this unknown function TTX analog as an olfactory chemoattractant. Estimation of differential cell cycle kinetics in higher plant root meristem with cellular fate and positional resolution. Plant root development is a complex spatial-temporal process that originates in the root apical meristem (RAM). To shape the organs structure signaling between the different cells and cell files must be highly coordinated. Thereby, diverging kinetics of chromatin remodeling and cell growth in these files need to be integrated and balanced by differential cell growth and local differences in cell proliferation frequency. Understanding the local differences in cell cycle duration in the RAM and its correlation with chromatin organization is crucial to build a holistic view on the different regulatory processes and requires a quantitative estimation of the chromatin geometry and underlying mitotic cell cycle phases timing at every cell file and every position. Unfortunately, so far precise methods for such analysis are missing.

This study presents a robust and straightforward pipeline to determine in parallel the duration of cell cycles key stages in all cell layers of a plants root and their nuclei organization. The methods combine marker-free techniques based on the detection of the nucleus, deep analysis of the chromatin phase transition, incorporation of 5-ethynyl-2'-deoxyuridine (EdU), and mitosis with a deep-resolution plant phenotyping platform to analyze all key cell cycle events kinetics.

In the Arabidopsis thaliana L. RAM S-phase duration was found to be as short as 20-30 minutes in all cell files. The subsequent G2-phase duration however depends on the cell type/position and varies from 3.5 hours in the pericycle to more than 4.5 hours in the epidermis. Overall, S+G2+M duration in Arabidopsis under our condition is 4 hours in the pericycle and up to 5.5 hours in the epidermis.

Endocycle duration was determined as the time required to achieve 100% EdU index in the transition zone and estimated to be in the range of 3-4 hours.

Besides Arabidopsis, we show that the presented technique is applicable also to root tips of other dicot and monocot plants (tobacco (Nicotiana tabacum L.), tomato (Lycopersicon esculentum L.) and wheat (Triticum aestivum L.). Decreased cell stiffness facilitates detachment and migration of breast cancer cells in 3D collagen matrices: An exploratory study BackgroundTumour-cell detachment is a critical early event in the metastatic cascade. However, the role of the cells mechanical properties in detachment and migration is not well understood. This exploratory study aimed to assess how intracellular stiffness changes these processes.

MethodsMDA-MB-231 cells were embedded as 10,000-cell spheroids in 2 and 4 mg/ml collagen matrices. Intracellular stiffness was assessed with mitochondria tracking microrheology of cells that migrated distances equivalent to four and six times the cell diameter (dC) from the spheroid and compared to cells at the spheroid surface (0dC), representing medium, high and no migration, respectively.

Results and discussionThe mitochondrial mean square displacement and intracellular stiffness decreased during detachment and migration for both collagen concentrations (i.e. rigidities). The mean square displacement of 4dC and 6dC cells was similar, whereas cell stiffness was lower for 4dC than for 6dC cells. With increasing matrix rigidity, the intracellular stiffness decreased for 0dC cells and did not change for 4dC and 6dC cells. It is proposed that decreased cell stiffness drives detachment and migration and increased matrix rigidity physically hinders migration, and cells need to become softer or remodel the environment to migrate. The independence of the stiffness of migrated cells from matrix rigidity suggests that cells remodel their environment through matrix proteins cleavage to migrate.

ConclusionsThe study revealed the collective effects of enhanced migratory conditions and increased matrix rigidity on the mechanical properties of the cells. The expression of matrix metalloproteinases and transforming growth factor {beta} and the role of cell volume on detachment and migration in matrices with varying pore sizes are proposed targets for further studies on metastatic cancer cells. Characterization of regulatory transcriptional mechanisms in hepatocyte lipotoxicity Non-alcoholic fatty liver disease is a continuum of disorders among which non-alcoholic steatohepatitis (NASH) is particularly associated with a negative prognosis. Hepatocyte lipotoxicity is one of the main pathogenic factors of liver fibrosis and NASH. However, the molecular mechanisms regulating this process are poorly understood. The main aim of this study was to dissect transcriptional mechanisms regulated by lipotoxicity in hepatocytes. We achieved this aim by combining transcriptomic, proteomic and chromatin accessibility analyses from human liver and mouse hepatocytes. This integrative approach revealed several transcription factor networks deregulated by NASH and lipotoxicity. To validate these predictions, genetic deletion of the transcription factors MAFK and TCF4 was performed, resulting in hepatocytes that were better protected against saturated fatty acid oversupply. MAFK- and TCF4-regulated gene expression profiles suggest a mitigating effect against cell stress, while promoting cell survival and growth. Moreover, in the context of lipotoxicity, some MAFK and TCF4 target genes were to the corresponding differentially regulated transcripts in human liver fibrosis. Collectively, our findings comprehensively profile the transcriptional response to lipotoxicity in hepatocytes, revealing new molecular insights and providing a valuable resource for future endeavours to tackle the molecular mechanisms of NASH. Interplay between receptor binding, immune escape and protein stability determines the natural selection of SARS-CoV-2 variants Emergence of new SARS-CoV-2 variants has raised concerns at the effectiveness of vaccines and antibody therapeutics developed against the unmutated wild-type virus. We examined the effect of 12 most commonly occurring mutations in the receptor binding domain on its expression, stability, activity, and antibody escape potential-some of the factors that may influence the natural selection of mutants. Recombinant proteins were expressed in human cells. Stability was measured using thermal denaturation melts. Activity and antibody escape potential were measured using isothermal titration calorimetry in terms of binding to ACE2 and to a neutralizing human antibody CC12.1, respectively. Our results show that variants differ in their expression levels with the two least stable variants showing lesser expression. Out of the 8 well-expressed mutants, only 2 (N501Y and K417T/E484K/N501Y) showed stronger affinity to ACE2, 4 (Y453F, S477N, T478I and S494P) have similar affinity, whereas the other 2 (K417N and E484K) have weaker affinity when compared to the wild-type. In terms of CC12.1 binding, when compared to the wild-type, 4 variants (K417N, Y453F, N501Y and K417T/E484K/N501Y) have weaker affinity, 2 (S477N and S494P) have similar affinity, and 2 (T478I and E484K) have stronger affinity. Taken together, these results indicate that multiple factors contribute towards the natural selection of variants, and all these factors need be considered to understand the evolution of the virus. In addition, since not all variants can escape a given neutralizing antibody, antibodies to treat new variants can be chosen based on the specific mutations in that variant. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Cytokinin oxidase/dehydrogenase family genes exhibit functional divergence and overlap in rice growth and development Cytokinins play key roles in plant growth and development; hence, cytokinin biosynthesis and degradation have been extensively studied. Cytokinin oxidase/dehydrogenases (CKXs) are a group of enzymes that regulate oxidative cleavage to maintain cytokinin homeostasis. In rice, 11 OsCKX genes have been identified to date; however, most of their functions remain unknown. Here, we comprehensively analyzed the expression patterns and functions of OsCKX genes. Using CRISPR/Cas9 technology, we constructed mutants of all OsCKX genes to determine the functions of OsCKXs in rice development. The results revealed that the single osckx and higher-order osckx4 osckx9 mutant lines showed functional overlap and subfunctionalization. Notably, the osckx1 osckx2 and osckx4 osckx9 double mutants displayed contrasting phenotypic changes in tiller number and panicle size compared to the wild type. Moreover, we identified several genes with significantly altered expression in osckx4 and osckx9 single and double mutant plants. Many differentially expressed genes were found to be associated with auxin and cytokinin pathways. Additionally, the cytokinins in osckx4 osckx9 mutants were increased compared to the wild type. Overall, our findings provide new insights into the functions of OsCKX genes in rice growth and may be used as a foundation for future studies aimed at improving rice yield. Key words: Cytokinin, expression pattern, OsCKX, panicle, phenotype, rice, tiller

HighlightThe osckx4 osckx9 double mutant had a significantly greater number of tillers, whereas the osckx1 osckx2 double mutant showed the opposite phenotypic change, compared to the wild type 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. Toxicity Profile of Artesunate in Rats and Dogs AO_SCPLOWBSTRACTC_SCPLOWO_ST_ABSOBJECTIVESC_ST_ABSThe objectives of these studies were to investigate the toxicity, safety and toxicokinetics of single and multiple doses of artesunate for injection in rats and dogs.

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

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

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

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

With the studies performed in rats and dogs for a treatment period of 28-days including measures of safety pharmacology, toxicokinetic evaluation and in vivo genotoxicity data, a comparison of toxicity and exposure of intravenous versus intramuscular artesunate for a prolonged treatment period is possible, delivering important information, when a decision between both routes has to be made for use in patients. Results of the peripheral blood micronucleus assay in rats were positive, suggesting genotoxic potential for the intravenous route of administration of artesunate, what is, beside patients, likewise relevant for production and health care personnel. Revisiting regulatory decoherence: Accounting for temporal bias in a co-expression analysis reveals novel candidates controlling environmental response Environment can alter the degree of phenotypic variation and covariation, potentially influencing evolutionary trajectories. However, environment-driven changes in phenotypic variation remain understudied. In an effort to exploit the abundance of RNASequencing data now available, an increasing number of ecological studies rely on population-level correlation to characterize the plastic response of the entire transcriptome and to identify environmentally responsive molecular pathways. These studies are fundamentally interested in identifying groups of genes that respond in concert to environmental shifts. We show that population-level differential co-expression exhibits biases when capturing changes of regulatory activity and strength in rice plants responding to elevated temperature. One possible cause of this bias is regulatory saturation, the observation that detectable co-variance between a regulator and its target may be low as their transcript abundances are induced. This phenomenon appears to be particularly acute for rapid-onset environmental stressors. However, our results suggest that temporal correlations may be a reliable means to detect transient regulatory activity following rapid onset environmental perturbations such as temperature stress. Such temporal bias is likely to confound the studies of phenotypic integration, where high-order organismal traits are hypothesized to be more integrated with strong correlation under stressful conditions, while recent transcriptome studies exhibited weaker coexpression between genes under stressful conditions. Collectively, our results point to the need to account for the nuances of molecular interactions and the possibly confounding effects that these can introduce into conventional approaches to study transcriptome datasets. Prenatal environmental conditions underlie alternative reproductive tactics that drive the formation of a mixed-kin cooperative society Although animal societies often evolve due to limited natal dispersal that results in kin clustering and facilitates cooperation among relatives, many species form cooperative groups with low and variable kin structure. Groups in such mixed-kin societies often comprise residents and immigrants of the same sex that compete for breeding opportunities. To understand how such mixed-kin societies form despite their potential for social conflict, we investigated the environmental causes and subsequent fitness consequences of dispersal decisions in male cooperatively breeding superb starlings (Lamprotornis superbus) that live in a climatically unpredictable environment. We show that the two alternative reproductive tactics--natal dispersal or philopatry--exhibit reproductive tradeoffs resulting in equal lifetime inclusive fitness. The tactic an individual adopts is governed by the environment its parents experience prior to laying rather than the environment it experiences during its juvenile stage. When individuals adopt the tactic not predicted by early life environmental conditions, their fitness is reduced. Ultimately, climate-driven oscillating selection may help stabilize mixed-kin animal societies, despite their reduced kin structure and potential for social conflict. Number and proportion of Plasmodium falciparum gametocytes vary from acute infection to chronic carriage despite unaltered sexual commitment rate In areas of seasonal malaria, resuming transmission every wet season relies on persistent asexual parasites during dry season that maintain the ability to produce gametocytes until the mosquito population resurges, following the start of annual rains. Although human asymptomatic P. falciparum reservoirs in the dry season are widely recognized, the longitudinal dynamics of parasite sexual commitment and gametocytogenesis are unclear. Here, we compared the density and proportion of P. falciparum gametocytes in blood of asymptomatic subjects during the dry season versus subjects with symptomatic malaria in the wet season. Blood concentrations of phospholipids potentially involved in gametocytogenesis and parasite transcriptomes were also compared. While blood densities of gametocytes and asexual parasites were lower during the dry season, we found that gametocytes were proportionally higher during the dry versus wet season. Levels of parasite transcripts involved in sexual commitment were unaltered throughout the year, and phospholipid content in the plasma throughout the year did not affect sexual commitment in vitro. Our experimental data and mathematical modelling support that gametocyte density and proportion diverge as infections progress from recently transmitted to chronic carriage, without significant alterations in the rate of sexual commitment over time. Global structural stability and the role of cooperation in mutualistic system Dynamical systems on graphs allow to describe multiple phenomena from different areas of Science. In particular, many complex systems in Ecology are studied by this approach. In this paper we analize the mathematical framework for the study of the structural stability of each stationary point, feasible or not, introducing a generalization for this concept, defined as Global Structural Stability. This approach would fit with the proper mathematical concept of structural stability, in which we find a full description of the complex dynamics on the phase space due to nonlinear dynamics. This fact can be analyzed as an informational field grounded in a global attractor whose structure can be completely characterized. These attractors are stable under perturbation and suppose the minimal structurally stable sets. We also study in detail, mathematically and computationally, the zones characterizing different levels of biodiversity in bipartite graphs describing mutualistic antagonistic systems of population dynamics. In particular, we investigate the dependence of the region of maximal biodiversity of a system on its connectivity matrix. On the other hand, as the network topology does not completely determine the robustness of the dynamics of a complex network, we study the correlation between structural stability and several graph measures. A systematic study on synthetic and biological graphs is presented, including 10 mutualistic networks of plants and seed-dispersal and 1000 random synthetic networks. We compare the role of centrality measures and modularity, concluding the importance of just cooperation strength among nodes when describing areas of maximal biodiversity. Indeed, we show that cooperation parameters are the central role for biodiversity while other measures act as secondary supporting functions.

Author summaryWe introduce the concept of Global Structural Stability as a proper mathematical concept to fully understand biodiversity in some ecological systems. Our concept retakes the definitions in the classical works of R. Thom [1] and Andronov-Pontryagin [2]. Moreover, there exists a close relation between the structure of a complex network, described as a graph, and its associated dynamics. Mutualistic networks introduce cooperation links between two groups of species, as plant and pollinators or seed-dispersal. The understanding of organizational aspects leading to maximizing biodiversity is one of the more important research areas in Theoretical and Applied Ecology. In this work we introduce a systematic study on different graph measures in order to identify optimal organization for maximal biodiversity (defined as structural stability). Our results conclude that, for mutualistic systems, the strength in cooperation parameters are the core fact, i.e., cooperation is the real fact optimizing biodiversity among other possible structural configurations. A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Stitching and registering highly multiplexed whole slide images of tissues and tumors using ASHLAR MotivationStitching microscope images into a mosaic is an essential step in the analysis and visualization of large biological specimens, particularly human and animal tissues. Recent approaches to highly-multiplexed imaging generate high-plex data from sequential rounds of lower-plex imaging. These multiplexed imaging methods promise to yield precise molecular single-cell data and information on cellular neighborhoods and tissue architecture. However, attaining mosaic images with single-cell accuracy requires robust image stitching and image registration capabilities that are not met by existing methods.

ResultsWe describe the development and testing of ASHLAR, a Python tool for coordinated stitching and registration of 103 or more individual multiplexed images to generate accurate whole-slide mosaics. ASHLAR reads image formats from most commercial microscopes and slide scanners, and we show that it performs better than existing open source and commercial software. ASHLAR outputs standard OME-TIFF images that are ready for analysis by other open-source tools and recently developed image analysis pipelines.

Availability and implementationASHLAR is written in Python and available under an MIT license at https://github.com/labsyspharm/ashlar. An informational website with user guides and test data is available at https://labsyspharm.github.io/ashlar/. Using total RNA quality metrics for time since deposition estimates in degrading bloodstains Determining the time since deposition (TSD) of bloodstains would provide forensic scientists with critical information regarding the timeline of the events involving bloodshed. The physicochemical changes occurring to biomolecules as a bloodstain degrades can be used to approximate the TSD of bloodstains. Our study aims to quantify the timewise degradation trends and temperature dependence found in total RNA from bloodstains without the use of amplification, expanding the scope of the RNA TSD research which has previously targeted mRNA molecules. Whole bovine blood was stored in plastic microcentrifuge tubes at 21{degrees}C or 4{degrees}C and tested over different timepoints spanning one week. Nine RNA metrics were visually assessed and quantified using linear and mixed models; the RNA Integrity Number equivalent (RINe) and the DV200 demonstrated strong negative trends over time and statistical independence. The RINe model fit was high (R2 = 0.60), and while including the biological replicate as a random effect increased the fit for all RNA metrics, no significant differences were found between biological replicates stored at the same temperature for the RINe and DV200 metrics. Importantly, this suggests that these standardized metrics can likely be directly compared between scenarios and individuals, with DV200 having an inflection point at [~]28 hrs. This study provides a novel approach for blood TSD estimates, producing metrics that are not affected by inter-individual variation and improving our understanding of the rapid degradation occurring in bloodstains.

HIGHLIGHTSO_LIAmplification-free analysis of total RNA in degrading bloodstains.
C_LIO_LIShort-term RNA degradation assessment using high-resolution size measurements.
C_LIO_LITotal RNA quality and quantity metrics were assessed across a one-week.
C_LIO_LITotal RNA quality metrics demonstrated the strongest timewise trends.
C_LIO_LIBiological replicates produced similar results for RNA quality metrics.
C_LI Robustness in regulatory networks of Epithelial Mesenchymal Plasticity as a function of positive and negative feedback loops Epithelial-Mesenchymal plasticity (EMP) is a key arm of cancer metastasis and is observed across many contexts. Cells undergoing EMP can reversibly switch between three classes of phenotypes: Epithelial (E), Mesenchymal (M), and Hybrid E/M. While a large number of multistable regulatory networks have been identified to be driving EMP in various contexts, the exact mechanisms and design principles that enable robustness in driving EMP across contexts are not yet fully understood. Here we investigated dynamic and structural robustness in EMP networks with regards to phenotypic distribution and plasticity. We use two different approaches to simulate these networks: a computationally inexpensive, parameter-independent continuous state space boolean model, and an ODE-based parameter-agnostic framework (RACIPE), both of which yield similar phenotypic distributions. Using perturbations to network topology and by varying network parameters, we show that multistable EMP networks are structurally and dynamically more robust as compared to their randomized counterparts, thereby highlighting their topological hallmarks. These features of robustness are governed by a balance of positive and negative feedback loops embedded in these networks. Using a combination of the number of negative and positive feedback loops weighted by their lengths and sign, we identified a metric that can explain the structural and dynamical robustness of these networks. This metric enabled us to compare networks across multiple sizes, and the network principles thus obtained can be used to identify fragilities in large networks without simulating their dynamics. Our analysis highlights a network topology-based approach to quantify robustness in multistable EMP networks.

Significance StatementEpithelial-Mesenchymal plasticity (EMP) is a key arm of cancer metastasis. Despite extensive intra- and inter-tumor heterogeneity, the characteristics of EMP have been observed to be robust across multiple contexts. We hypothesize that topology of EMP regulatory networks contributes towards this robustness. Here, we measure the robustness of EMP in the form of its phenotypic heterogeneity and multistability and show that EMP networks are more robust to dynamical (change in kinetic parameters) and structural (change in network topology) perturbations as compared to their random network counterparts. Furthermore, we propose a network topology-based metric using the nature and length of feedback loops that explains the observed robustness. Our metric hence serves to quantify robustness in multistable EMP networks without simulating their dynamics. DeGlyPHER: an ultrasensitive method for analysis of viral spike N-glycoforms Viruses can evade the host immune system by displaying numerous glycans on their surface "spike-proteins" that cover immune epitopes. We have developed an ultrasensitive "single pot" method to assess glycan occupancy and the extent of glycan processing from high-mannose to complex forms at each N-glycosylation site. Though aimed at characterizing glycosylation of viral spike-proteins as potential vaccines, this method is applicable for analysis of site-specific glycosylation of any glycoprotein.