HTSQualC is a Flexible and One-Step Quality Control Software for High-throughput Sequencing Data Analysis Use of high-throughput sequencing (HTS) has become indispensable in life science research. Raw HTS data contains several sequencing artifacts, and as a first step it is imperative to remove the artifacts for reliable downstream bioinformatics analysis. Although there are multiple stand-alone tools available that can perform the various quality control steps separately, availability of an integrated tool that can allow one-step, automated quality control analysis of HTS datasets will significantly enhance handling large number of samples parallelly. Here, we developed HTSQualC, a stand-alone, flexible, and easy-to-use software for one-step quality control analysis of raw HTS data. HTSQualC can evaluate HTS data quality and perform filtering and trimming analysis in a single run. We evaluated the performance of HTSQualC for conducting batch analysis of HTS datasets with 322 samples with an average [~]1M (paired end) sequence reads per sample. HTSQualC accomplished the QC analysis in [~]3 hours in distributed mode and [~]31 hours in shared mode, thus underscoring its utility and robust performance. In addition to command-line execution, we integrated HTSQualC into the free, open-source, CyVerse cyberinfrastructure resource as a GUI interface, for wider access to experimental biologists who have limited computational resources and/or programming abilities.
Accommodating site variation in neuroimaging data using normative and hierarchical Bayesian models AO_SCPLOWBSTRACTC_SCPLOWThe potential of normative modeling to make individualized predictions from neuroimaging data has enabled inferences that go beyond the case-control approach. However, site effects are often confounded with variables of interest in a complex manner and can bias estimates of normative models, which has impeded the application of normative models to large multi-site neuroimaging data sets. In this study, we suggest accommodating for these site effects by including them as random effects in a hierarchical Bayesian model. We compared the performance of a linear and a non-linear hierarchical Bayesian model in modeling the effect of age on cortical thickness. We used data of 570 healthy individuals from the ABIDE (autism brain imaging data exchange) data set in our experiments. In addition, we used data from individuals with autism to test whether our models are able to retain clinically useful information while removing site effects. We compared the proposed single stage hierarchical Bayesian method to several harmonization techniques commonly used to deal with additive and multiplicative site effects using a two stage regression, including regressing out site and harmonizing for site with ComBat, both with and without explicitly preserving variance related to age and sex as biological variation of interest. In addition, we made predictions from raw data, in which site has not been accommodated for. The proposed hierarchical Bayesian method showed the best predictive performance according to multiple metrics. Beyond that, the resulting z-scores showed little to no residual site effects, yet still retained clinically useful information. In contrast, performance was particularly poor for the regression model and the ComBat model in which age and sex were not explicitly modeled. In all two stage harmonization models, predictions were poorly scaled, suffering from a loss of more than 90 % of the original variance. Our results show the value of hierarchical Bayesian regression methods for accommodating site variation in neuroimaging data, which provides an alternative to harmonization techniques. While the approach we propose may have broad utility, our approach is particularly well suited to normative modelling where the primary interest is in accurate modelling of inter-subject variation and statistical quantification of deviations from a reference model.

1 HighlightsO_LIDevelopment and presentation of normative modeling approach based on hierarchical Bayesian modeling that can be applied to large multi-site neuroimaging data sets.
C_LIO_LIComparison of performance of Hierarchical Bayesian model including site as predictor to several common ways to harmonize for multi-site effects.
C_LIO_LIPresentation of normative modeling as site correction tool.
C_LI
A cell-based screen identifies HDAC inhibitors as activators of RIG-I signaling Enhancing the immune microenvironment in cancer by targeting the nucleic acid sensors is becoming a potent therapeutic strategy. Among the nucleic acid sensors, activation of the RNA sensor Retinoic Acid-inducible Gene (RIG-I) using small hairpin RNAs has been shown to elicit powerful innate and adaptive immune responses. Given the challenges inherent in pharmacokinetics and delivery of RNA based agonists, we set out to discover small molecule agonists of RIG-I using a cell-based assay. To this end, we established and validated a robust high throughput screening assay based on a commercially available HEK293 reporter cell line with a luciferase reporter downstream of tandem interferon stimulated gene 54 (ISG54) promoter elements. We first confirmed that the luminescence in this cell line is dependent on RIG-I and the interferon receptor using a hairpin RNA RIG-I agonist. We established a 96-well and a 384-well format HTS based on this cell line and performed a proof-of-concept screen using an FDA approved drug library of 1200 compounds. Surprisingly, we found two HDAC inhibitors Entinostat, Mocetinostat and the PLK1 inhibitor Volasertib significantly enhanced ISG-luciferase activity. This luminescence was substantially diminished in the null reporter cell line indicating the increase in signaling was dependent on RIG-I expression. Combination treatment of tumor cell lines with Entinostat increased RIG-I induced cell death in a mammary carcinoma cell line that is resistant to either Entinostat or RIG-I agonist alone. Taken together, our data indicates an unexpected role for HDAC1,-3 inhibitors in enhancing RIG-I signaling and highlight potential opportunities for therapeutic combinations. SERCA function is impaired in skeletal and cardiac muscles from young DBA/2J mdx mice The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. Force profile of the two-handed hardstyle kettlebell swing performed by an RKC-certified Instructor BackgroundThe effects of hardstyle kettlebell training are increasingly cited in strength and conditioning research, yet reference data from a proficient swing is scarce. The aim of this exploratory study was to investigate the force profile of a two-handed hardstyle swing performed by a Russian Kettlebell Challenge (RKC) instructor.

MethodsThe subject is a 44-year-old male, body mass 75.6 kg, height 173.5 cm, with six years of regular hardstyle training experience. Two-handed hardstyle swings were performed with a series of incremental mass kettlebells (8-68 kg). Ground reaction force (GRFs) was obtained from a floor-mounted force platform. Force-time curves (FTCs), peak force, forward force, rate of force development (RFD) and swing cadence were investigated.

ResultsData revealed the FTC of a proficient swing is highly consistent and dominated by a single force peak (mean SD = 47 N), with a profile that remained largely unchanged to 24 kg. Pearson correlation analysis revealed a very strong positive correlation in peak force with kettlebell mass (r = 0.95), which increased disproportionately from the lightest to heaviest kettlebells; net peak force increased from 8.36 {+/-} 0.75 N.kg-1 (0.85 x BW) to 12.82 {+/-} 0.39 N.kg-1 (1.3x BW). There was a strong negative correlation between RFD and kettlebell mass (r = 0.82) that decreased from 39.2 N.s-1.kg-1 to 21.5 N.s-1.kg-1. There was a very strong positive correlation in forward ground reaction force with kettlebell mass (r = 0.99), expressed as a ratio of vertical ground reaction, that increased from 0.092 (9.2%) to 0.205 (20.5%). Swing cadence exceeded 40 swings per minute (SPM) with all kettlebells.

ConclusionOur findings challenge some of the popular beliefs of the hardstyle kettlebell swing. Consistent with hardstyle practice, and previous kinematic analysis of expert and novice, force-time curves show a characteristic single large force peak, differentiating passive from active shoulder flexion. Ground reaction force did not increase proportionate to kettlebell mass, with a magnitude of forward force smaller than described in practice. These results could be useful for coaches and trainers wanting to improve athletic performance, and healthcare providers using the kettlebell swing for therapeutic purposes. Findings from this study were used to inform the BELL Trial, a pragmatic controlled trial of kettlebell training with older adults. www.anzctr.org.au ACTRN12619001177145. 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. Conformational decoupling in acid-sensing ion channels uncovers mechanism and stoichiometry of PcTx1-mediated inhibition Acid-sensing ion channels (ASICs) are trimeric proton-gated cation channels involved in fast synaptic transmission. Pharmacological inhibition of ASIC1a reduces neurotoxicity and stroke infarct volumes, with the cysteine knot toxin Psalmotoxin-1 (PcTx1) being one of the most potent and selective inhibitors. PcTx1 binds at the subunit interface in the extracellular domain (ECD), but the mechanism and conformational consequences of the interaction, as well as the number of toxin molecules required for inhibition remain unknown. Here we use voltage-clamp fluorometry and subunit concatenation to decipher the mechanism and stoichiometry of PcTx1 inhibition of ASIC1a. Besides the known inhibitory binding mode, we propose PcTx1 to have at least two additional binding modes that are decoupled from the pore. One of these modes induces a long-lived ECD conformation that reduces the activity of an endogenous neuropeptide. This long-lived conformational state is proton-dependent and can be destabilized by a mutation that decreases PcTx1 sensitivity. Lastly, the use of concatemeric channel constructs reveal that disruption of a single PcTx1 binding site is sufficient to destabilize the toxin-induced conformation, while functional inhibition is not impaired until two or more binding sites are mutated. Together, our work provides insight into the mechanism of PcTx1 inhibition of ASICs and uncovers a prolonged conformational change with possible pharmacological implications. Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Human 14-3-3 proteins site-selectively bind the mutational hotspot region of SARS-CoV-2 nucleoprotein modulating its phosphoregulation Phosphorylated within its Ser/Arg-rich region, the SARS-CoV-2 nucleoprotein (N) recruits the phosphopeptide-binding human 14-3-3 proteins that play a well-recognized role in replication of many viruses. Here we use genetic code expansion to demonstrate that phosphorylation of SARS-CoV-2 N at either of two pseudo-repeats centered at Ser197 and Thr205 is sufficient for 14-3-3 binding. According to fluorescence anisotropy, the pT205-motif, present in SARS-CoV-2 but not in SARS-CoV, is preferred over the pS197-motif by all seven human 14-3-3 isoforms, which display unforeseen pT205/pS197 binding selectivity hierarchy. Crystal structures demonstrate that pS197 and pT205 are mutually exclusive 14-3-3-binding sites, whereas SAXS and biochemical data indicate 14-3-3 binding occludes the Ser/Arg-rich region, inhibiting its dephosphorylation. This Ser/Arg-rich region of N is highly prone to mutations, as exemplified by the Omicron and Delta variants, with our data suggesting how the strength of its 14-3-3 binding can be linked with the replicative fitness of the virus. The Hippocampus uses cryptographic principles rather than plasticity to enable animals to secretly cache and retrieve their food For animals, the ability to hide and retrieve valuable information, such as the location of food, can mean the difference between life and death. Here, we propose that to achieve this, their brain uses spatial cells similarly to how we utilize encryption for data security. Some animals are able to cache hundreds of thousands of food items annually by each individual and later retrieve most of what they themselves stashed. Rather than memorizing their cache locations as previously suggested, we propose that they use a single cryptographic-like mechanism during both caching and retrieval. The model we developed is based on hippocampal spatial cells, which respond to an animals positional attention, such as when the animal enters a specific region (place-cells) or gazes at a particular location (spatial-view-cells). We know that the region that activates each spatial cell remains consistent across subsequent visits to the same area but not between areas. This remapping, combined with the uniqueness of cognitive maps, produces a persistent crypto-hash function for both food caching and retrieval. We also show that the model stores temporal information that helps animals in food caching order preference as previously observed. This behavior, which we refer to as crypto-taxis, might also explain consistent differences in decision-making when animals are faced with a large number of alternatives such as in foraging. Ultrahigh-Throughput Ambient MS: Direct Analysis at 22 Samples per Second by Infrared Matrix-Assisted Desorption Electrospray Ionization Mass Spectrometry Infrared Matrix-Assisted Laser Desorption Electrospray Ionization (IR-MALDESI) mass spectrometry is an ambient-direct sampling method being developed for high-throughput, label-free, biochemical screening of large-scale compound libraries. Here, we report the development of an ultrahigh-throughput continuous motion IR-MALDESI sampling approach capable of acquiring data at rates up to 22.7 samples per second in a 384-well microtiter plate. At top speed, less than 1% analyte carryover is observed from well-to-well and signal intensity relative standard deviations (RSD) of 11.5% and 20.9% for 3 M 1-hydroxymidazolam and 12 M dextrorphan, respectively, are achieved. The ability to perform parallel kinetics studies on 384 samples with ~30s time resolution using an isocitrate dehydrogenase 1 (IDH1) enzyme assay is shown. Finally, we demonstrate the repeatability and throughput of our approach by measuring 115,200 samples from 300 microtiter plate reads consecutively over 5.54 hours with RSDs under 8.14% for each freshly introduced plate. Taken together, these results demonstrate the use of IR-MALDESI at sample acquisition rates that surpass other currently reported direct sampling mass spectrometry approaches used for high throughput compound screening.

For Table of Contents Only

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=125 SRC="FIGDIR/small/465730v3_ufig1.gif" ALT="Figure 1">
View larger version (49K):
org.highwire.dtl.DTLVardef@5f916corg.highwire.dtl.DTLVardef@4473b4org.highwire.dtl.DTLVardef@82e8forg.highwire.dtl.DTLVardef@128b57b_HPS_FORMAT_FIGEXP  M_FIG C_FIG First evidence of in vitro cytotoxic effects of marine microlitter on Merluccius merluccius and Mullus barbatus, two Mediterranean commercial fish species Marine litter is composed mainly of plastics and is recognized as a serious threats to marine ecosystems. Ecotoxicological approaches have started elucidating the potential severity of microplastics (MPs) in controlled laboratory studies with pristine materials but no information exist on marine environmental microlitter as a whole. Here, we characterized the litter in the coastal Northern Tyrrhenian sea and in the stomach of two fish species of socio-economic importance, and exposed primary cell cultures of mucosal and lymphoid organs to marine microlitter for evaluating possible cytotoxic effects. An average of 0.30 {+/-} 0.02 microlitter items m-3 was found in water samples. FT-IR analysis revealed that plastic particles, namely HDPE, polyamide and polypropylene were present in 100% and 83.3% of Merluccius merluccius and Mullus barbatus analyzed, which overall ingested 14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual, respectively. Moreover, microlitter was confirmed as a vector of microorganisms. Lastly, the apical end-point of viability was found to be significantly reduced in splenic cells exposed in vitro to two microlitter conditions. Considering the role of the spleen in the mounting of adaptive immune responses, our results warrant more in-depth investigations for clarifying the actual susceptibility of these two species to anthropogenic microlitter.

HighlightsO_LI0.30 {+/-} 0.02 microlitter items m-3 were found at the surface of coastal Northern Tyrrhenian sea
C_LIO_LI14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual were retrieved from the stomach of hakes and mullets
C_LIO_LIThe ingested microlitter contained plastic items
C_LIO_LIMicrolitter was validated as a carrier of bacteria, fungi and flagellates
C_LIO_LISplenic cells exposed to two microlitter conditions for 72 hours suffered cytotoxicity
C_LI Dynamic modulation of enhancer responsiveness by core promoter elements in living Drosophila embryo Regulatory interactions between enhancers and core promoters are fundamental for the temporal and spatial specificity of gene expression in development. The central role of core promoters is to initiate productive transcription in response to enhancers activation cues. However, it has not been systematically assessed how individual core promoter elements affect the induction of transcriptional bursting by enhancers. Here, we provide evidence that each core promoter element dynamically alters "enhancer responsiveness" by changing functional parameters of transcriptional bursting in developing Drosophila embryos. Quantitative live imaging analysis revealed that the timing and the continuity of burst induction are common regulatory steps on which core promoter elements impact. We also show that core promoters use the upstream TATA box as an amplifier of transcriptional bursting upon interaction with enhancers. Genome editing analysis of the pair-rule gene fushi tarazu further revealed that the endogenous TATA box and the downstream core promoter element are both essential for the correct expression and the function of key developmental genes in early embryos. We suggest that core promoter elements serve as a key regulatory module in converting enhancer activity into transcription dynamics during animal development. Key features of the genetic architecture and evolution of host-microbe interactions revealed by high-resolution genetic mapping of the mucosa-associated gut microbiome in hybrid mice Determining the forces that shape diversity in host-associated bacterial communities is critical to understanding the evolution and maintenance of metaorganisms. To gain deeper understanding of the role of host genetics in shaping gut microbial traits, we employed a powerful genetic mapping approach using inbred lines derived from the hybrid zone of two incipient house mouse species. Further, we uniquely performed our analysis on microbial traits measured at the gut mucosal interface, which is in more direct contact with host cells and the immune system. A high number of mucosa-associated bacterial taxa have significant heritability estimates; heritabilities are greater for 16S rRNA transcript-compared to gene copy-based traits, and interestingly, are positively correlated with cospeciation rate estimates. Genomewide association mapping identifies 443 loci influencing 123 taxa, with narrow genomic intervals pinpointing promising candidate genes and pathways. Importantly, we identified an enrichment of candidate genes associated with several human diseases, including inflammatory bowel disease, and functional categories including innate immunity and G-protein-coupled receptors. These results highlight key features of the genetic architecture of mammalian host-microbe interactions and how they diverge as new species form. A novel Fiji/ImageJ plugin for the rapid analysis of blebbing cells When confined, cells have recently been shown to undergo a phenotypic switch to what has been termed, fast amoeboid (leader bleb-based) migration. However, as this is a nascent area of research, few tools are available for the rapid analysis of cell behavior. Here, we demonstrate that a novel Fiji/ImageJ-based plugin, Analyze_Blebs, can be used to quickly obtain cell migration parameters and morphometrics from time lapse images. As validation, we show that Analyze_Blebs can detect significant differences in cell migration and morphometrics, such as the largest bleb size, upon introducing different live markers of F-actin, including F-tractin and LifeAct tagged with green and red fluorescent proteins. We also demonstrate, using flow cytometry, that live markers increase total levels of F-actin. Furthermore, that F-tractin increases cell stiffness, which was found to correlate with a decrease in migration, thus reaffirming the importance of cell mechanics as a determinant of Leader Bleb-Based Migration (LBBM).

Highlight summaryA new plugin, Analyze_Blebs, enables the rapid analysis of cell migration and morphometrics of fast amoeboid cells. Morphometrics combined with cell stiffness are found to be predictive of confined, fast amoeboid (leader bleb-based), migration. Optogenetic actuator/ERK biosensor circuits identify MAPK network nodes that shape ERK dynamics Combining single-cell measurements of ERK activity dynamics with perturbations provides insights into the MAPK network topology. We built circuits consisting of an optogenetic actuator to activate MAPK signaling and an ERK biosensor to measure single-cell ERK dynamics. This allowed us to conduct RNAi screens to investigate the role of 50 MAPK proteins in ERK dynamics. We found that the MAPK network is robust against most node perturbations. We observed that the ERK-RAF and the ERK-RSK2-SOS negative feedbacks operate simultaneously to regulate ERK dynamics. Bypassing the RSK2-mediated feedback, either by direct optogenetic activation of RAS, or by RSK2 perturbation, sensitized ERK dynamics to further perturbations. Similarly, targeting this feedback in a human ErbB2-dependent oncogenic signaling model increased the efficiency of a MEK inhibitor. The RSK2-mediated feedback is thus important for the ability of the MAPK network to produce consistent ERK outputs and its perturbation can enhance the efficiency of MAPK inhibitors. Metacommunity dynamics and the detection of species associations in co-occurrence analyses: why patch disturbance matters O_LI1: Many statistical methods attempt to detect species associations - and so infer inter-specific interactions - from species co-occurrence patterns. Habitat heterogeneity and out-of-equilibrium colonization histories are well recognized as potentially causing species associations, even when interactions are absent. The potential for patch disturbance, a classical component of metacommunity dynamics, to also drive spurious species associations has however been over-looked.
C_LIO_LI2: Using a new general metacommunity model, we derive mathematical predictions regarding how patch disturbance would affect the patterns of species associations detected in "null" co-occurrence matrices. We also conduct numerical simulations to test our predictions and to compare the performance of several widespread statistical methods, including direct tests of pairwise independence, matrix permutation approaches and joint species distribution modelling.
C_LIO_LI3: We show how classical metacommunity dynamics can produce statistical associations, both positive and negative, even when species do not interact, when there is no habitat heterogeneity, and at equilibrium. This occurs as soon as there is some rate of patch disturbance (i.e. simultaneous extinction of several species in a patch) and/or a finite life-span of patches, a common feature of a broad range of plant, animal or microbial systems.
C_LIO_LI4: Patch disturbance can compromise species co-occurrence analyses and cause the artefactual detection of species associations if not taken into account. Including patch age (i.e. the time since the last patch disturbance event) as a covariate in a joint species distribution model can resolve the artefact. However, this requires additional data that often are not available in practice. We argue that the consequences of patch disturbance should not be underestimated when analysing species distribution patterns in metacommunity-like systems.
C_LI Temporal responses of bumblebee gustatoryneurons to sugars The sense of taste permits the recognition of valuable nutrients and the avoidance of potential toxins. Models of gustatory coding propose that within modalities (e.g. sweet, bitter, etc.), taste ligands are not distinct stimuli. However, these models are based on data from mice or flies that have omnivorous, non-specialist diets. A specialist feeder might, however, be expected to have acuity within modality if stimulus resolution was critical to survival. Previously, we found that bumblebees have a specialized mechanism for sensing sugars whereby two gustatory receptor neurons (GRNs) within the galeal sensilla of the bees mouthparts exhibit burst of spikes. Here, we show that the temporal firing patterns of these GRNs separate sugars into four distinct groups that correlate with sugar nutritional value and palatability. We also identified a third GRN that responded to stimulation with relatively high concentrations of fructose, sucrose, and maltose. Sugars that were non-metabolizable or toxic suppressed the responses of bursting GRNs to sucrose. These abilities to encode information about sugar value are a refinement to the bumblebees sense of sweet taste that could be an adaptation that enables precise calculations of the nature and nutritional value of floral nectar. The Ensemble of Gene Regulatory Networks at Mutation-Selection Balance The evolution of diverse phenotypes both involves and is constrained by molecular interaction networks. When these networks influence patterns of expression, we refer to them as gene regulatory networks (GRNs). Here, we develop a quasi-species model of GRN evolution. With this model, we prove that-across a broad spectrum of viability and mutation functions-the dynamics converge to a stationary distribution over GRNs. Next, we show from first principles how the frequency of GRNs at equilibrium will be proportional to each GRNs eigenvector centrality in the genotype network. Finally, we determine the structural characteristics of GRNs that are favored in response to a range of selective environments and mutational constraints. Our work connects GRN evolution to quasi-species models, and thus can provide a mechanistic explanation for the topology of GRNs experiencing various evolutionary forces. Experimental evolution reveals post-transcriptional regulation as a novel driver of Leishmania fitness gain The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network driving parasite fitness, with genome instability causing highly reproducible, gene dosage-dependent changes in protein abundance linked to post-transcriptional regulation. These in turn were associated with a gene dosage-independent reduction in abundance of flagellar transcripts and a coordinated increase in abundance of coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of fitness-adapted ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings. Inflammatory response in hematopoietic stem and progenitor cells triggered by activating SHP2 mutations evokes blood defects Gain-of-function mutations in the protein-tyrosine phosphatase SHP2 are the most frequently occurring mutations in sporadic juvenile myelomonocytic leukemia (JMML) and JMML-like myeloproliferative neoplasm (MPN) associated with Noonan syndrome (NS). Hematopoietic stem and progenitor cells (HSPCs) are the disease propagating cells of JMML. Here, we explored transcriptomes of HSPCs with Shp2 mutations derived from JMML patients and a novel NS zebrafish model. In addition to major NS traits, CRISPR/Cas9 knock-in Shp2D61G mutant zebrafish recapitulated a JMML-like MPN phenotype, including myeloid lineage hyperproliferation, ex vivo growth of myeloid colonies and in vivo transplantability of HSPCs. Single cell mRNA sequencing of HSPCs from Shp2D61G zebrafish embryos and bulk sequencing of HSPCs from JMML patients revealed an overlapping inflammatory gene expression pattern. Strikingly, an anti-inflammatory agent rescued JMML-like MPN in Shp2D61G zebrafish embryos. Our results indicate that a common inflammatory response was triggered in the HSPCs from sporadic JMML patients and syndromic NS zebrafish, which potentiated MPN and may represent a future target for JMML therapies. Genetic variation associated with plastic and homeostatic growth responses to drought in Arabidopsis O_LINatural genetic variation influences plant responses to environmental stressors. However, the extent to which such variation underlies plastic versus homeostatic response phenotypes has been little studied.
C_LIO_LIWe quantified the extent of drought-induced changes in leaf area in a set of Iberian Arabidopsis accessions and then performed association studies correlating variation in this phenotype with genomic and transcriptomic variation.
C_LIO_LIDrought-induced plastic reductions in relative leaf area typified accessions originating from productive environments. In contrast, homeostasis in relative leaf area typified accessions originating from unproductive environments. Genome-Wide Association Studies (GWAS), Transcriptome Wide Association Studies (TWAS), and expression GWAS (eGWAS) highlighted the importance of auxin-related processes and, in particular, the potential role of the SMALL AUXIN UP RNA 26 (SAUR26) gene in conferring leaf area plasticity. Homeostatic responses in relative leaf area were associated with a diverse gene set and positively associated with a higher intrinsic water use efficiency (WUEi) as confirmed in a TWAS metanalysis of previously published {delta}13C measurements.
C_LIO_LIWe have identified not only candidate "plasticity genes" but also candidate "homeostasis genes" controlling leaf area. Our results demonstrate the value of a combined GWAS, TWAS, and eGWAS approach to identify mechanisms underlying phenotypic responses to stress.
C_LI

One-sentence summaryInformation on phenotype, genotype, and transcript abundance is integrated to identify both plasticity and homeostasis genes and processes associated with local adaptation to drought stress in Arabidopsis accessions of the Iberian Peninsula. Immunogenicity of SARS-CoV-2 trimetric spike protein associated to Poly(I:C) plus Alum The SARS-CoV-2 pandemic has had a social and economic impact worldwide, and vaccination is an efficient strategy for diminishing those damages. New adjuvant formulations are required for the high vaccine demands, especially adjuvant formulations that induce a Th1 phenotype. Herein we assess a vaccination strategy using a combination of Alum and polyinosinic:polycytidylic acid (Poly(I:C)) adjuvants plus the SARS-CoV-2 spike protein in a prefusion trimeric conformation by an intradermal (ID) route. We found high levels of IgG anti-spike antibodies in the serum by enzyme linked immunosorbent assay (ELISA) and high neutralizing titers against SARS-CoV-2 in vitro by neutralization assay, after one or two boosts. By evaluating the production of IgG subtypes, as expected, we found that formulations containing Poly(I:C) induced IgG2a whereas Alum did not. The combination of these two adjuvants induced high levels of both IgG1 and IgG2a. In addition, cellular immune responses of CD4+ and CD8+ T cells producing interferon-gamma were equivalent, demonstrating that the Alum + Poly(I:C) combination supported a Th1 profile. Based on the high neutralizing titers, we evaluated B cells in the germinal centers, which are specific for receptor-binding domain (RBD) and spike, and observed that more positive B cells were induced upon the Alum + Poly(I:C) combination. Moreover, these B cells produced antibodies against both RBD and non-RBD sites. We also studied the impact of this vaccination preparation (spike protein with Alum + Poly(I:C)) in the lungs of mice challenged with inactivated SARS-CoV-2 virus. We found a production of IgG, but not IgA, and a reduction in neutrophil recruitment in the bronchoalveolar lavage fluid (BALF) of mice, suggesting that our immunization scheme reduced lung inflammation. Altogether, our data suggest that Alum and Poly(I:C) together is a possible adjuvant combination for vaccines against SARS-CoV-2 by the intradermal route. Alkaline phosphatase activity of serum affects osteogenic differentiation cultures Fetal bovine serum (FBS) is a widely used supplement in cell culture medium, despite its known variability in composition which greatly affects cellular function and consequently the outcome of studies. In bone tissue engineering, the deposited mineralized matrix is one of the main outcome parameters, but using different brands of FBS can result in large variations. Alkaline phosphatase (ALP) is present in FBS. Not only is ALP used to judge the osteogenic differentiation of bone cells, it may affect deposition of mineralized matrix. The present study focused on the enzymatic activity of ALP in FBS of different suppliers and its contribution to mineralization in osteogenic differentiation cultures. It was hypothesized that culturing cells in a medium with high intrinsic ALP activity of FBS will lead to higher mineral deposition compared to media with lower ALP activity. The used FBS types were shown to have significant differences in enzymatic ALP activity. Our results indicate that the ALP activity of the medium not only affected the deposited mineralized matrix but also the osteogenic differentiation of cells as measured by a changed cellular ALP activity of human bone marrow derived mesenchymal stromal cells (hBMSC). In media with low inherent ALP activity, the cellular ALP activity was increased and played the major role in the mineralization process; while, in media with high intrinsic ALP activity contribution from the serum, less cellular ALP activity was measured and the ALP activity of the medium also contributed to mineral formation substantially. Our results highlight the diverse effects of ALP activity intrinsic to FBS on osteogenic differentiation and matrix mineralization and how FBS can determine the experimental outcomes, in particular for studies investigating matrix mineralization. Once again, the need to replace FBS with more controlled and known additives is highlighted.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=55 SRC="FIGDIR/small/462693v2_ufig1.gif" ALT="Figure 1">
View larger version (17K):
org.highwire.dtl.DTLVardef@902078org.highwire.dtl.DTLVardef@43b6fborg.highwire.dtl.DTLVardef@1a244fborg.highwire.dtl.DTLVardef@5ef819_HPS_FORMAT_FIGEXP  M_FIG C_FIG Early epigenetic influence on EGR1, KLF2 and KLF4 transcriptional network associates with alteration of HSCs quiescence Excessive fetal growth is associated with DNA methylation alterations in human hematopoietic stem and progenitor cells (HSPC), but their functional impact remains elusive. We implemented an integrative analysis combining single-cell epigenomics, single-cell transcriptomics, and in vitro analyses to functionally link DNA methylation changes to putative alterations of HSPC functions. We showed in hematopoietic stem cells (HSC) from large for gestational age neonates that both DNA hypermethylation and chromatin rearrangement target a specific network of transcription factors known to sustain stem cell quiescence. In parallel, we found a decrease expression of key genes regulating HSC differentiation including EGR1, KLF2, SOCS3, and JUNB. Our functional analyses showed that this epigenetic programming was associated with a decreased ability for HSCs to stay quiescent. Taken together, our multimodal approach using single-cell (epi)genomics showed that human fetal overgrowth affects hematopoietic stem cells quiescence maintenance via epigenetic programming. Reconstitution of minimal motility system based on Spiroplasma swimming by expressing two bacterial actins in synthetic minimal bacterium Motility is one of the most important features of life, but its evolutionary origin remains unknown. In this study, we focused on Spiroplasma, commensal, or parasitic bacteria. They swim by switching the helicity of a ribbon-like cytoskeleton that comprises six proteins, each of which evolved from a nucleosidase and bacterial actin called MreB. We expressed these proteins in a synthetic, non-motile minimal bacterium, JCVI-syn3.0B, whose reduced genome was computer-designed and chemically synthesized. The synthetic bacterium exhibited swimming motility with features characteristic of Spiroplasma swimming. Moreover, some combinations of the two proteins produced a helical cell shape and swimming, suggesting that the swimming originated from the differentiation and coupling of bacterial actins, and we also obtained a minimal system for motility of the synthetic bacterium.

One-Sentence SummaryThe minimal system comprised two bacterial actins that provided cell helicity and swimming to the synthetic minimal bacterium.