High performance single-cell gene regulatory network inference at scale: The Inferelator 3.0 MotivationGene regulatory networks define regulatory relationships between transcription factors and target genes within a biological system, and reconstructing them is essential for understanding cellular growth and function. Methods for inferring and reconstructing networks from genomics data have evolved rapidly over the last decade in response to advances in sequencing technology and machine learning. The scale of data collection has increased dramatically; the largest genome-wide gene expression datasets have grown from thousands of measurements to millions of single cells, and new technologies are on the horizon to increase to tens of millions of cells and above.

ResultsIn this work, we present the Inferelator 3.0, which has been significantly updated to integrate data from distinct cell types to learn context-specific regulatory networks and aggregate them into a shared regulatory network, while retaining the functionality of the previous versions. The Inferelator is able to integrate the largest single-cell datasets and learn cell-type specific gene regulatory networks. Compared to other network inference methods, the Inferelator learns new and informative Saccharomyces cerevisiae networks from single-cell gene expression data, measured by recovery of a known gold standard. We demonstrate its scaling capabilities by learning networks for multiple distinct neuronal and glial cell types in the developing Mus musculus brain at E18 from a large (1.3 million) single-cell gene expression dataset with paired single-cell chromatin accessibility data.

AvailabilityThe inferelator software is available on GitHub (https://github.com/flatironinstitute/inferelator) under the MIT license and has been released as python packages with associated documentation (https://inferelator.readthedocs.io/).
A theoretical model of mitochondrial ATP synthase deficiencies. The role of mitochondrial carriers. The m.8993T>G mutation of the mitochondrial MT-ATP6 gene is associated with NARP syndrome (Neuropathy, Ataxia and Retinitis Pigmentosa). The equivalent point mutation introduced in yeast Saccharomyces cerevisiae mitochondrial DNA considerably reduced the activity of ATP synthase and of cytochrome-c-oxidase preventing yeast growth on oxidative substrates. The over-expression of the mitochondrial oxodicarboxylate carrier (Odc1p) is able to rescue the growth on oxidative substrate by increasing the substrate-level phosphorylation of ADP coupled to conversion of -ketoglutarate (AKG) into succinate with an increase in Complex IV activity. Previous studies showed that equivalent point mutations in ATP synthase behave similarly and can be rescued by Odc1p overexpression and/or uncoupler of OXPHOS from ATP synthesis. In order to better understand the mechanism of ATP synthase mutation bypass, we developed a core model of mitochondrial metabolism based on AKG as respiratory substrate. We describe the different possible metabolite output and the ATP/O ratio values as a function of ATP synthase inhibition.
Enhancers regulate 3' end processing activityto control expression of alternative 3'UTR isoforms Multi-UTR genes are widely transcribed but show cell type-specific 3'UTR isoform expression. As transcriptional enhancers regulate mRNA expression, we investigated if they also regulate mRNA isoform expression. Deletion of an endogenous enhancer of a multi-UTR gene did not impair transcript production but prevented a switch in 3'UTR isoform expression. Also, the same enhancers are able to increase transcript production in the context of single-UTR gene promoters, but they increase 3' end processing activity when paired with multi-UTR gene promoters. We show that transcription factors regulate processing activity of weak polyadenylation sites to control cell type-specific alternative 3'UTR isoform expression of widely expressed genes. Inter-tissue convergence of gene expression during ageing suggests age-related loss of tissue and cellular identity Developmental trajectories of gene expression may reverse in their direction during ageing, a phenomenon previously linked to cellular identity loss. Our analysis of cerebral cortex, lung, liver and muscle transcriptomes of 16 mice, covering development and ageing intervals, revealed widespread but tissue-specific ageing-associated expression reversals. Cumulatively, these reversals create a unique phenomenon: mammalian tissue transcriptomes diverge from each other during postnatal development, but during ageing, they tend to converge towards similar expression levels, a process we term Divergence followed by Convergence, or DiCo. We found that DiCo was most prevalent among tissue-specific genes and associated with loss of tissue identity, which is confirmed using data from independent mouse and human datasets. Further, using publicly available single-cell transcriptome data, we showed that DiCo could be driven both by alterations in tissue cell type composition and also by cell-autonomous expression changes within particular cell types. Circulating multimeric immune complexes drive immunopathology in COVID-19 A dysregulated immune response with high levels of SARS-CoV-2 specific IgG antibodies characterizes patients with severe or critical COVID-19. Although a robust IgG response is traditionally considered to be protective, excessive triggering of activating Fc-gamma-receptors (Fc{gamma}Rs) could be detrimental and cause immunopathology. Here, we document that patients who develop soluble circulating IgG immune complexes (sICs) during infection are subject to enhanced immunopathology driven by Fc{gamma}R activation. Utilizing cell-based reporter systems we provide evidence that sICs are predominantly formed prior to a specific humoral response against SARS-CoV-2. sIC formation, together with increased afucosylation of SARS-CoV-2 specific IgG eventually leads to an enhanced CD16 (Fc{gamma}RIII) activation of immune cells reaching activation levels comparable active systemic lupus erythematosus (SLE) disease. Our data suggest a vicious cycle of escalating immunopathology driven by an early formation of sICs in predisposed patients. These findings reconcile the seemingly paradoxical findings of high antiviral IgG responses and systemic immune dysregulation in severe COVID-19.

Clinical implicationsThe identification of sICs as drivers of an escalating immunopathology in predisposed patients opens new avenues regarding intervention strategies to alleviate critical COVID-19 progression.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=78 SRC="FIGDIR/small/449893v4_ufig1.gif" ALT="Figure 1">
View larger version (17K):
org.highwire.dtl.DTLVardef@b4616corg.highwire.dtl.DTLVardef@682d1aorg.highwire.dtl.DTLVardef@16946cborg.highwire.dtl.DTLVardef@a6ef7d_HPS_FORMAT_FIGEXP  M_FIG C_FIG A vicious cycle of immunopathology in COVID-19 patients is driven by soluble multimeric immune complexes (sICs). SARS-CoV-2 infection triggers sIC formation in prone individuals. Activation of Fc{gamma}RIII/CD16 expressing immune cells by sICs precedes a humoral response to SARS-CoV2 infection. sICs and infection add to IgG afucosylation, further enhancing Fc{gamma}RIII/CD16 activation by opsonized targets. High inflammation induces further sIC mediated immune cell activation ultimately leading to an escalating immunopathology. 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. Genetic interactions drive heterogeneity in causal variant effect sizes for gene expression and complex traits Despite the growing number of genome-wide association studies (GWAS), it remains unclear to what extent gene-by-gene and gene-by-environment interactions influence complex traits in humans. The magnitude of genetic interactions in complex traits has been difficult to quantify because GWAS are generally underpowered to detect individual interactions of small effect. Here, we develop a method to test for genetic interactions that aggregates information across all trait-associated loci. Specifically, we test whether SNPs in regions of European ancestry shared between European American and admixed African American individuals have the same causal effect sizes. We hypothesize that in African Americans, the presence of genetic interactions will drive the causal effect sizes of SNPs in regions of European ancestry to be more similar to those of SNPs in regions of African ancestry. We apply our method to two traits: gene expression in 296 African Americans and 482 European Americans in the Multi-Ethnic Study of Atherosclerosis (MESA) and low-density lipoprotein cholesterol (LDL-C) in 74K African Americans and 296K European Americans in the Million Veteran Program (MVP). We find significant evidence for genetic interactions in our analysis of gene expression; for LDL-C, we observe a similar point estimate although this is not significant, likely due to lower statistical power. These results suggest that gene-by-gene or gene-by-environment interactions modify the effect sizes of causal variants in human complex traits. Tourmaline: a containerized workflow for rapid and iterable amplicon sequence analysis using QIIME 2 and Snakemake BackgroundAmplicon sequencing (metabarcoding) is a common method to survey diversity of environmental communities whereby a single genetic locus is amplified and sequenced from the DNA of whole or partial organisms, organismal traces (e.g., skin, mucus, feces), or microbes in an environmental sample. Several software packages exist for analyzing amplicon data, among which QIIME 2 has emerged as a popular option because of its broad functionality, plugin architecture, provenance tracking, and interactive visualizations. However, each new analysis requires the user to keep track of input and output file names, parameters, and commands; this lack of automation and standardization is inefficient and creates barriers to meta-analysis and sharing of results.

FindingsWe developed Tourmaline, a Python-based workflow that implements QIIME 2 and is built using the Snakemake workflow management system. Starting from a configuration file that defines parameters and input files--a reference database, a sample metadata file, and a manifest or archive of FASTQ sequences--it uses QIIME 2 to run either the DADA2 or Deblur denoising algorithm, assigns taxonomy to the resulting representative sequences, performs analyses of taxonomic, alpha, and beta diversity, and generates an HTML report summarizing and linking to the output files. Features include support for multiple cores, automatic determination of trimming parameters using quality scores, representative sequence filtering (taxonomy, length, abundance, prevalence, or ID), support for multiple taxonomic classification and sequence alignment methods, outlier detection, and automated initialization of a new analysis using previous settings. The workflow runs natively on Linux and macOS or via a Docker container. We ran Tourmaline on a 16S rRNA amplicon dataset from Lake Erie surface water, showing its utility for parameter optimization and the ability to easily view interactive visualizations through the HTML report, QIIME 2 viewer, and R- and Python-based Jupyter notebooks.

ConclusionsAutomated workflows like Tourmaline enable rapid analysis of environmental and biomedical amplicon data, decreasing the time from data generation to actionable results. Tourmaline is available for download at github.com/aomlomics/tourmaline. Opposing mechanisms of competition determine species invasions and functional diversity Understanding how species phenotypic differences affect competition is key to explaining community assembly and predicting biodiversity responses. Many studies overlook the variable effects that species trait differences can have on the direction of competitive exclusion, which reverses depending on the specific mechanism at play. We performed a comprehensive trait-based study of an ant invasion integrating morphological, dietary, physiological and behavioral analyses. We found that trait differences between invasive and resident species were not only associated with niche differences which promoted the coexistence of dissimilar species, but also competitive ability differences which acted in the opposite fashion. Furthermore, competition along separate trait axes led to complex and contrasting patterns in the invaded assemblages, where species were at once similar (clustered) in some traits but also dissimilar (overdispersed) in others. Our results reveal that different aspects of phenotype may distinctly modulate the effect of competition in structuring ecological communities and functional diversity. 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. Bacterial glycyl tRNA synthetase offers glimpses of ancestral protein topologies Aminoacyl tRNA synthetases (aaRSs) are among the proposed proteins present in the Last Universal Common Ancestor (LUCA). There are two types of glycyl tRNA synthetases (GlyRSs), from which the archaeal-eukaryal type is the one suggested to be present in LUCA. Here we solved the crystal structure of a complete bacterial glycyl tRNA synthetase (bacGlyRS) and show that indeed, bacGlyRS carries several structural signals that point it at the origin of all aaRSs. Furthermore, if bacGlyRS is ancestral, it should help to build a reliable Tree of Life (ToL). Given the modular nature of protein evolution, we used only two sub-domain segments with duplicated ancestral topologies, no detected orthologs and an assumed limited horizontal gene transfer (HGT). These motifs correspond to the non-specific RNA binding regions of contemporary bacGlyRS, archaeal CCA-adding enzyme (arch-CCAadd), and eukaryotic rRNA processing enzyme (euk-rRNA). The calculated, rooted bacterial ToL agrees with several phyla relationships unaccounted by the available trees. Once-daily feeding is associated with better health in companion dogs: Results from the Dog Aging Project A variety of diets have been studied for possible anti-aging effects. In particular, studies of intermittent fasting and time-restricted feeding in laboratory rodents have found evidence of beneficial health outcomes. Companion dogs represent a unique opportunity to study diet in a large mammal that shares human environments. The Dog Aging Project has been collecting data on thousands of companion dogs of all different ages, sizes, and breeds since 2019. We leveraged this diverse cross-sectional dataset to investigate associations between feeding frequency and cognitive function (n = 10,474) as well as nine broad categories of health conditions (n = 24,238). Controlling for sex, age, breed, and other potential confounders, we found that dogs fed once daily rather than more frequently had lower mean scores on a cognitive dysfunction scale, and lower odds of having gastrointestinal, dental, orthopedic, kidney/urinary, and liver/pancreas disorders. Therefore, we find that once-daily feeding is associated with better health in multiple domains. Future research with longitudinal data can provide stronger evidence for a possible causal effect of feeding frequency on health in companion dogs. An early signaling GARP transcription factor regulates the entry point for differentiation in Giardia Transcriptional regulation of differentiation is critical for parasitic pathogens to adapt to environmental changes and regulate transmission. In response to encystation stimuli, Giardia lamblia cells shift from G1+S phase of the cell cycle to G2+M. By 2-4 hours, cyst wall proteins are upregulated, indicating that key regulatory steps occur within the first four hours of encystation. However, all characterized transcription factors (TFs) in Giardia have only been investigated at later time points of encystation. How TFs initiate encystation and link it to the cell cycle remains enigmatic. Here, we systematically screened six putative early upregulated TFs for nuclear localization, established their dynamic expression profiles, and determined their functional role in regulating encystation. We found a critical repressor, GLP4 that increases rapidly after 30 min of encystation stimuli and downregulates encystation specific markers, including Cyst Wall Proteins (CWPs) and enzymes in the cyst N-acetyl galactosamine (GalNAc) pathway. Depletion of GLP4 increases cyst production. Importantly, we observe that G2+M cells exhibit higher levels of CWP1 resulting from the activation of MYB2, a TF previously linked to encystation in Giardia. GLP4 upregulation occurs in G1+S cells, suggesting a role in repressing MYB2 and encystation specific genes in the G1+S phase of the cell cycle. Furthermore, we demonstrate that depletion of GLP4 upregulates MYB2 and promotes encystation while overexpression of GLP4 downregulates MYB2 and represses encystation. Together, these results suggest that Giardia employs a dose-dependent transcriptional response that involves the cell cycle regulated repressor GLP4 to orchestrate MYB2 and entry into the encystation pathway.

ImportanceTransition between life cycle stages is a common feature among parasitic pathogens and its regulation must be optimized to balance persistence of infection with transmission. The early transcription factors (TFs) regulating commitment to differentiate are totally unknown in Giardia. In this work, we identified GLP4, a previously uncharacterized GARP-like TF, as an early-acting transcriptional repressor that inhibits G1+S cells from entering the encystation pathway. GLP4 is therefore a key regulator controlling the balance between proliferative growth and terminal differentiation into infective cysts. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Reliable  in silico  ranking of engineered therapeutic TCR binding affinities with MMPBSA and MMGBSA. Accurate and efficient in silico ranking of protein-protein binding affinities is useful for protein design with applications in biological therapeutics. One popular approach to rank binding affinities is to apply the molecular mechanics Poisson Boltzmann/generalized Born surface area (MMPB/GBSA) method to molecular dynamics trajectories. Here, we identify protocols that enable the reliable evaluation of T-cell receptor (TCR) variants binding to their target, peptide-human leukocyte antigens (pHLAs). We suggest different protocols for variant sets with few ([&le;]4) or many mutations, with entropy corrections important for the latter. We demonstrate how potential outliers could be identified in advance and that just 5-10 replicas of short (4 ns) MD simulations may be sufficient for reproducible and accurate ranking of TCR variants. The protocols developed here can be applied towards in silico screening during the optimization of therapeutic TCRs, potentially reducing both the cost and time taken for biologic development. S-acylation stabilizes ligand-induced receptor kinase complex formation during plant pattern-triggered immune signalling. Plant receptor kinases are key transducers of physical extracellular stimuli, such as the presence of beneficial or pathogenic microbes or secreted signalling molecules. Receptor kinases are regulated by numerous post-translational modifications. Here, using the bacterial flagellin perceiving receptor kinase FLS2, we show that S-acylation at an evolutionarily conserved cysteine is crucial for function. S-acylation involves the addition of long-chain fatty acids to cysteine residues within proteins, altering their biophysical properties and behaviour within the membrane environment. We observe S-acylation of FLS2 at C-terminal kinase domain cysteine residues within minutes of treatment with flg22 ligand, and in a BAK1 co-receptor dependent manner. We demonstrate that S-acylation is essential for FLS2-mediated immune signalling, including anti-bacterial immunity. Similarly, mutating the corresponding conserved cysteine residue in the immune receptor kinase EFR supressed elf18 mediated signalling. Biochemical analysis of unstimulated and activated FLS2 containing complexes using detergents and native membrane DIBMA nanodiscs indicates that S-acylation assists the stabilization of activated receptor kinase complexes withing the membrane environment to increase signalling efficiency. Combinatorial drug-microenvironment interaction mapping reveals cell-extrinsic drug resistance mechanisms and clinically relevant patient subgroups in CLL The tumour microenvironment and genetic alterations collectively influence drug efficacy in cancer, but current evidence is limited to small scale studies and systematic analyses are lacking. We chose Chronic Lymphocytic Leukaemia (CLL), the most common leukaemia in adults, as a model disease to study this complex interplay systematically. We performed a combinatorial assay using 12 drugs individually co-applied with each of 17 microenvironmental stimuli in 192 primary CLL samples, generating a comprehensive map of drug-microenvironment interactions in CLL. This data was combined with whole-exome sequencing, DNA-methylation, RNA-sequencing and copy number variant annotation. Our assay identified four distinct CLL subgroups that differed in their responses to the panel of microenvironmental stimuli. These subgroups were characterized by distinct clinical outcomes independently of known prognostic markers. We investigated the effect of CLL- specific recurrent genetic alterations on microenvironmental responses and identified trisomy 12 as an amplifier of multiple microenvironmental stimuli. We further quantified the impact of microenvironmental stimuli on drug response, confirmed known interactions such as Interleukin (IL) 4 mediated resistance to B cell receptor (BCR) inhibitors, and identified new interactions such as Interferon-{gamma} induced resistance to BCR inhibitors. Finally, we identified interactions which were limited to genetic subgroups. Resistance to chemotherapeutics, such as Fludarabine, induced by Toll-Like Receptor (TLR) agonists could be observed in IGHV unmutated patient samples and IGHV mutated samples with trisomy 12. In-vivo relevance was investigated in CLL-infiltrated lymph nodes, which showed increased IL4 and TLR signalling activity compared to healthy samples (p<0.001). High IL4 activity in lymph nodes correlated with faster disease progression (p=0.038).

We provide a publicly available resource (www.dietrichlab.de/CLL_Microenvironment/) which uncovers tumour cell extrinsic influences on drug response and disease progression in CLL, and how these interactions are modulated by cell intrinsic molecular features. Dimerization of kringle 1 domain from hepatocyte growth factor/scatter factor provides a potent minimal MET receptor agonist. Hepatocyte Growth Factor/Scatter Factor (HGF/SF) and its cognate receptor MET play several essential roles in embryogenesis and regeneration in post-natal life of epithelial organs such as liver, kidney, lung, and pancreas, prompting a strong interest in harnessing HGF/SF-MET signalling for regeneration of epithelial organs after acute or chronic damage. The limited stability and tissue diffusion of native HGF/SF, however, which reflect the tightly controlled, local mechanism of action of the morphogen, have led to a major search for HGF/SF mimics for therapy. In this work, we describe the rational design, production and characterisation of K1K1, a novel minimal MET agonist consisting of two copies of the kringle 1 domain of HGF/SF placed in tandem. K1K1 is highly stable and displays biological activities equivalent or superior to native HGF/SF in a variety of in vitro assay systems and in a mouse model of liver disease. These data suggest that this engineered ligand may find wide applications in acute and chronic diseases of the liver and other epithelial organs dependent on MET activation. Electrophysiological correlates of mood and reward dynamics in human adolescents Despite its omnipresence in everyday interactions and its importance for mental health, mood and its neuronal underpinnings are poorly understood. Computational models can help identify parameters affecting self-reported mood during mood induction tasks. Here we test if computationally modelled dynamics of self-reported mood during monetary gambling can be used to identify trial-by-trial variations in neuronal activity. To this end, we shifted mood in healthy (N=24) and depressed (N=30) adolescents by delivering individually tailored reward prediction errors whilst recording magnetoencephalography (MEG) data. Following a pre-registered analysis, we hypothesize that expectation (defined by previous reward outcomes) would be predictive of beta-gamma oscillatory power (25-40Hz), a frequency shown to modulate to reward feedback. We also hypothesize that trial variations in the evoked response to the presentation of gambling options and in source localized responses to reward feedback. Through our multilevel statistical analysis, we found confirmatory evidence that beta-gamma power is positively related to reward expectation during mood shifts, with possible localized sources in the posterior cingulate cortex. We also confirmed reward prediction error to be predictive of trial-level variations in the response of the paracentral lobule and expectation to have an effect on the cerebellum after presentation of gambling options. To our knowledge, this is the first study to relate fluctuations in mood on a minute timescale to variations in neural oscillations with noninvasive electrophysiology.

Significance StatementBrain mechanisms underlying mood and its relationship with changes in reward contingencies in the environment are still elusive but could have a strong impact on our understanding and treatment of debilitating mood disorders. Building on a previously proposed computational mood model we use multilevel statistical models to find relationship between trial-by-trial variations in model components of mood and neural responses to rewards measured with non-invasive electrophysiology (MEG). Through confirmatory analysis we show that it is possible to observe relationships between trial variations in neural responses and computational parameters describing mood dynamics. Identifying the dynamics of mood and the neural processes it affects could pave the way for more effective neuromodulation treatments. Impacts of an urban sanitation intervention on fecal indicators and the prevalence of human fecal contamination in Mozambique Fecal source tracking (FST) may be useful to assess pathways of fecal contamination in domestic environments and to estimate the impacts of water, sanitation, and hygiene (WASH) interventions in low-income settings. We measured two non-specific and two human-associated fecal indicators in water, soil, and surfaces before and after a shared latrine intervention from low-income households in Maputo, Mozambique participating in the Maputo Sanitation (MapSan) trial. Up to a quarter of households were impacted by human fecal contamination, but trends were unaffected by improvements to shared sanitation facilities. The intervention reduced E. coli gene concentrations in soil but did not impact culturable E. coli or the prevalence of human FST markers in a difference-in-differences analysis. Using a novel Bayesian hierarchical modeling approach to account for human marker diagnostic sensitivity and specificity, we revealed a high amount of uncertainty associated with human FST measurements and intervention effect estimates. The field of microbial source tracking would benefit from adding measures of diagnostic accuracy to better interpret findings, particularly when FST analyses convey insufficient information for robust inference. With improved measures, FST could help identify dominant pathways of human and animal fecal contamination in communities and guide implementation of effective interventions to safeguard health.

SYNOPSISAn urban sanitation intervention had minimal and highly uncertain effects on human fecal contamination after accounting for fecal indicator sensitivity and specificity.

TOC GRAPHIC/ABSTRACT ART

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=108 SRC="FIGDIR/small/432000v2_ufig1.gif" ALT="Figure 1">
View larger version (25K):
org.highwire.dtl.DTLVardef@312d2dorg.highwire.dtl.DTLVardef@8cf2caorg.highwire.dtl.DTLVardef@197b549org.highwire.dtl.DTLVardef@fefeb3_HPS_FORMAT_FIGEXP  M_FIG C_FIG N-VEGF, the forgotten isoform of VEGF-A, activate cellular stress survival circuits Vascular endothelial growth factor A (VEGF-A) is a secreted protein that stimulates angiogenesis in response to hypoxia. Under hypoxic conditions, a non-canonical long isoform called L-VEGF is concomitantly expressed with VEGF-A. Once translated, L-VEGF it is proteolytically cleaved to generate N-VEGF and VEGF-A. Interestingly, while VEGF-A is secreted and affects the surrounding cells, N-VEGF is mobilized to the nucleus. This suggests that N-VEGF participates in transcriptional response to hypoxia. In this study, we performed a series of complementary experiments to examine the functional role of N-VEGF. Strikingly, we found that the mere expression of N-VEGF followed by its hypoxia-independent mobilization to the nucleus was sufficient to induce key genes associated with angiogenesis, such as Hif1, VEGF-A isoforms, as well as genes associated with cell survival under hypoxia. Complementarily, when N-VEGF was genetically depleted, key hypoxia-induced genes were downregulated and cells were significantly susceptible to hypoxia-mediated apoptosis. This is the first reports of N-VEGF serving as an autoregulatory arm of VEGF-A. Further experiments will be needed to determine the role of N-VEGF in cancer and embryogenesis. 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. 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. B1 SINE-binding ZFP266 impedes reprogramming through suppression of chromatin opening mediated by pioneering factors Induced pluripotent stem cell reprogramming is inherently inefficient and understanding the molecular mechanisms underlying this inefficiency holds the key to successfully control cellular identity. Here, we report 16 novel reprogramming roadblock genes identified by CRISPR/Cas9-mediated genome-wide knockout (KO) screening. Of these, depletion of the predicted KRAB zinc finger protein (KRAB-ZFP) Zfp266 strongly and consistently enhanced iPSC generation in several iPSC reprogramming settings, emerging as the most robust roadblock. Further analyses revealed that ZFP266 binds Short Interspersed Nuclear Elements (SINEs) adjacent to binding sites of pioneering factors, OCT4 (POU5F1), SOX2 and KLF4, and impedes chromatin opening. Replacing the KRAB co-suppressor with a co-activator domain converted ZFP266 from a reprogramming inhibitor to a potent reprogramming facilitator. This work proposes SINE-KRAB-ZFP interaction to be a critical regulator of chromatin accessibility at enhancers for efficient cellular identity changes and also serves as a resource to further illuminate molecular mechanisms hindering reprogramming. Sex determination gene transformer regulates the male-female difference in Drosophila fat storage via the adipokinetic hormone pathway Sex differences in whole-body fat storage exist in many species. For example, Drosophila females store more fat than males. Yet, the mechanisms underlying this sex difference in fat storage remain incompletely understood. Here, we identify a key role for sex determination gene transformer (tra) in regulating the male-female difference in fat storage. Normally, a functional Tra protein is present only in females, where it promotes female sexual development. We show that loss of Tra in females reduced whole-body fat storage, whereas gain of Tra in males augmented fat storage. Tras role in promoting fat storage was largely due to its function in neurons, specifically the Adipokinetic hormone (Akh)-producing cells (APCs). Our analysis of Akh pathway regulation revealed a male bias in APC activity and Akh pathway function, where this sex-biased regulation influenced the sex difference in fat storage by limiting triglyceride accumulation in males. Importantly, Tra loss in females increased Akh pathway activity, and genetically manipulating the Akh pathway rescued Tra-dependent effects on fat storage. This identifies sex-specific regulation of Akh as one mechanism underlying the male-female difference in whole-body triglyceride levels, and provides important insight into the conserved mechanisms underlying sexual dimorphism in whole-body fat storage. 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.