Genome-scale metabolic modeling reveals SARS-CoV-2-induced metabolic changes and antiviral targets Tremendous progress has been made to control the COVID-19 pandemic caused by the SARS-CoV-2 virus. However, effective therapeutic options are still rare. Drug repurposing and combination represent practical strategies to address this urgent unmet medical need. Viruses, including coronaviruses, are known to hijack host metabolism to facilitate viral proliferation, making targeting host metabolism a promising antiviral approach. Here, we describe an integrated analysis of 12 published in vitro and human patient gene expression datasets on SARS-CoV-2 infection using genome-scale metabolic modeling (GEM), revealing complicated host metabolism reprogramming during SARS-CoV-2 infection. We next applied the GEM-based metabolic transformation algorithm to predict anti-SARS-CoV-2 targets that counteract the virus-induced metabolic changes. We successfully validated these targets using published drug and genetic screen data and by performing an siRNA assay in Caco-2 cells. Further generating and analyzing RNA-sequencing data of remdesivir-treated Vero E6 cell samples, we predicted metabolic targets acting in combination with remdesivir, an approved anti-SARS-CoV-2 drug. Our study provides clinical data-supported candidate anti-SARS-CoV-2 targets for future evaluation, demonstrating host metabolism-targeting as a promising antiviral strategy.
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.
LY-CoV1404 (bebtelovimab) potently neutralizes SARS-CoV-2 variants SARS-CoV-2 neutralizing monoclonal antibodies (mAbs) can reduce the risk of hospitalization when administered early during COVID-19 disease. However, the emergence of variants of concern has negatively impacted the therapeutic use of some authorized mAbs. Using a high throughput B-cell screening pipeline, we isolated a highly potent SARS-CoV-2 spike glycoprotein receptor binding domain (RBD)-specific antibody called LY-CoV1404 (also known as bebtelovimab). LY-CoV1404 potently neutralizes authentic SARS-CoV-2 virus, including the prototype, B.1.1.7, B.1.351 and B.1.617.2). In pseudovirus neutralization studies, LY-CoV1404 retains potent neutralizing activity against numerous variants including B.1.1.7, B.1.351, B.1.617.2, B.1.427/B.1.429, P.1, B.1.526, B.1.1.529, and the BA.2 subvariant and retains binding to spike proteins with a variety of underlying RBD mutations including K417N, L452R, E484K, and N501Y. Structural analysis reveals that the contact residues of the LY-CoV1404 epitope are highly conserved with the exception of N439 and N501. Notably, the binding and neutralizing activity of LY-CoV1404 is unaffected by the most common mutations at these positions (N439K and N501Y). The breadth of reactivity to amino acid substitutions present among current VOC together with broad and potent neutralizing activity and the relatively conserved epitope suggest that LY-CoV1404 has the potential to be an effective therapeutic agent to treat all known variants causing COVID-19.

In BriefLY-CoV1404 is a potent SARS-CoV-2-binding antibody that neutralizes all known variants of concern and whose epitope is rarely mutated.

HighlightsO_LILY-CoV1404 potently neutralizes SARS-CoV-2 authentic virus and known variants of concern including the B.1.1.529 (Omicron), the BA.2 Omicron subvariant, and B.1.617.2 (Delta) variants
C_LIO_LINo loss of potency against currently circulating variants
C_LIO_LIBinding epitope on RBD of SARS-CoV-2 is rarely mutated in GISAID database
C_LIO_LIBreadth of neutralizing activity and potency supports clinical development
C_LI Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable and highly reproducible in terms of microbial composition, function and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation. C4-dicarboxylate metabolons: Interaction of C4-dicarboxylate transporters of Escherichia coli with cytosolic enzymes Metabolons represent the structural organization of proteins for metabolic or regulatory pathways. Here the interaction of enzymes fumarase FumB and aspartase AspA with the C4-DC transporters DcuA and DcuB of Escherichia coli was tested by a bacterial two-hybrid (BACTH) assay in situ, or by co-chromatography (mSPINE). DcuB interacted strongly with FumB and AspA, and DcuA with AspA. The fumB-dcuB and the dcuA-aspA genes encoding the respective proteins are known for their colocalization on the genome and the production of co-transcripts. The data consistently suggest the formation of DcuB/FumB, DcuB/AspA and DcuA/AspA metabolons in fumarate respiration for the uptake of L-malate, or L-aspartate, conversion to fumarate and excretion of succinate after reduction. The DcuA/AspA metabolon catalyzes L-Asp uptake and fumarate excretion in concerted action also to provide ammonia for nitrogen assimilation. Proteasome complexes experience profound structural and functional rearrangements throughout mammalian spermatogenesis During spermatogenesis, spermatogonia undergo a series of mitotic and meiotic divisions on their path to spermatozoa. To achieve this, a succession of processes requiring high proteolytic activity are in part orchestrated by the proteasome. The spermatoproteasome (s20S) is specific to the developing gametes, in which the gamete-specific 4s subunit replaces the 4 isoform found in the constitutive proteasome (c20S). Although the s20S is conserved across species and was shown to be crucial for germ cell development, its mechanism, function and structure remain incompletely characterized. Here, we used advanced mass spectrometry (MS) methods to map the composition of proteasome complexes and their interactomes throughout spermatogenesis. We observed that the s20S becomes highly activated as germ cells enter meiosis, mainly through a particularly extensive 19S activation and, to a lesser extent, PA200 binding. Additionally, the proteasome population shifts from c20S (98%) to s20S (>82-92%) during differentiation, presumably due to the shift from 4 to 4s expression. We demonstrated that s20S, but not c20S, interacts with components of the meiotic synaptonemal complex, where it may localize via association with the PI31 adaptor protein. In vitro, s20S preferentially binds to 19S, and displays higher trypsin- and chymotrypsin-like activities, both with and without PA200 activation. Moreover, using MS methods to monitor protein dynamics, we identified significant differences in domain flexibility between 4 and 4s. We propose that these differences induced by 4s incorporation result in significant changes in the way the s20S interacts with its partners, and dictate its role in germ cell differentiation. LPS-TLR4 pathway exaggerates alcohol-induced acute-on-chronic liver failure via provoking NETs formation Backgrounds: Intrahepatic infiltration of neutrophils is a character of alcoholic acute-on-chronic liver failure (AACLF) and neutrophil extracellular traps (NETs) are an important strategy for neutrophils to fix and kill invading microorganisms. Intestinal bacteria and the gut-liver axis have been thought to play a key role in many liver diseases also including AACLF. However, whether NETs appear in AACLF and play a role in AACLF is still unsure. Methods: WT, NE KO, and TLR4 KO mice were used to build the AACLF model, and the intestinal bacteria were eliminated at the same time and LPS was given. Then the formation of NETs and AACLF related markers were detected. Results: The serum MPO-DNA and LPS concentration was increased in AACLF patients and a correlation was revealed between these two indexes. More intrahepatic NETs formed in AACLF mice by testing MPO-DNA, Cit H3, and NE. These markers decreased with gut detergent and restored markers with gut detergent plus LPS supplement. While NETs formation failed to change with gut microbiome or combine LPS supplement in TLR4 KO mice. As we tested AACLF related characters, liver injury, intrahepatic fat deposition, inflammation, and fibrosis alleviated with depletion of NE. These related marks were also attenuated with gut sterilization by antibiotics and recovered with combined treatment with antibiotics plus LPS. But the liver injury, intrahepatic fat, fibro deposition, and liver inflammation-related markers did show a significant difference in TLR4 KO mice when they received the same treatment. Conclusion: Intestinal-derived LPS promotes NETs formation in AACLF through the TLR4 pathway and further accelerates the AACLF process by NETs. How grower decision-making is modelled impacts the effect of plant disease control. While the spread of plant disease depends strongly on biological factors driving transmission, it also has a human dimension. Disease control depends on decisions made by individual growers, who are in turn influenced by a broad range of factors. Despite this, human behaviour has rarely been included in plant epidemic models. Considering Cassava Brown Streak Disease, we model how the perceived increase in profit due to disease management influences participation in clean seed systems (CSS). Our models are rooted in game theory, with growers making strategic decisions based on the expected profitability of different control strategies. We find that both the information used by growers to assess profitability and the perception of economic and epidemiological parameters influence long-term participation in the CSS. Over-estimation of infection risk leads to lower participation in the CSS, as growers perceive that paying for the CSS will be futile. Additionally, even though good disease management can be achieved through the implementation of CSS, and a scenario where all controllers use the CSS is achievable when growers base their decision on the average of their entire strategy, CBSD is rarely eliminated from the system. These results are robust to stochastic and spatial effects. Our work highlights the importance of including human behaviour in plant disease models, but also the significance of how that behaviour is included. Hierachical Resampling for Bagging in Multi-Study Prediction with Applications to Human Neurochemical Sensing We propose the "study strap ensemble," which combines advantages of two common approaches to fitting prediction models when multiple training datasets ("studies") are available: pooling studies and fitting one model versus averaging predictions from multiple models each fit to individual studies. The study strap ensemble fits models to bootstrapped datasets, or "pseudo-studies." These are generated by resampling from multiple studies with a hierarchical resampling scheme that generalizes the randomized cluster bootstrap. The study strap is controlled by a tuning parameter that determines the proportion of observations to draw from each study. When the parameter is set to its lowest value, each pseudo-study is resampled from only a single study. When it is high, the study strap ignores the multi-study structure and generates pseudo-studies by merging the datasets and drawing observations like a standard bootstrap. We empirically show the optimal tuning value often lies in between, and prove that special cases of the study strap draw the merged dataset and the set of original studies as pseudo-studies. We extend the study strap approach with an ensemble weighting scheme that utilizes information in the distribution of the covariates of the test dataset.

Our work is motivated by neuroscience experiments using real-time neurochemical sensing during awake behavior in humans. Current techniques to perform this kind of research require measurements from an electrode placed in the brain during awake neurosurgery and rely on prediction models to estimate neurotransmitter concentrations from the electrical measurements recorded by the electrode. These models are trained by combining multiple datasets that are collected in vitro under heterogeneous conditions in order to promote accuracy of the models when applied to data collected in the brain. A prevailing challenge is deciding how to combine studies or ensemble models trained on different studies to enhance model generalizability.

Our methods produce marked improvements in simulations and in this application. All methods are available in the studyStrap CRAN package. Lifespan prolonging mechanisms and insulin upregulation without fat accumulation in long-lived reproductives of a higher termite Kings and queens of eusocial termites can live for decades, while queens sustain a nearly maximal fertility. To investigate the molecular mechanisms underlying their long lifespan, we carried out transcriptomics, lipidomics and metabolomics in Macrotermes natalensis on sterile short-lived workers, long-lived kings and five stages spanning twenty years of adult queen maturation. Reproductives share gene expression differences from workers in agreement with a reduction of several aging-related processes, involving upregulation of DNA damage repair and mitochondrial functions. Anti-oxidant gene expression is downregulated, while peroxidability of membranes in queens decreases. Against expectations, we observed an upregulated gene expression in fat bodies of reproductives of several components of the IIS pathway, including an insulin-like peptide, Ilp9. This pattern does not lead to deleterious fat storage in physogastric queens, while simple sugars dominate in their hemolymph and large amounts of resources are allocated towards oogenesis. Our findings support the notion that all processes causing aging need to be addressed simultaneously in order to prevent it. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. Oct1 recruits the histone lysine demethylase Utx to promote lineage specification The pathways used by cells to transition between undifferentiated, pluripotent state and tissue-specific states are incompletely understood. Here we show that the widely-expressed transcription factor Oct1/Pou2f1 activates silent, developmental lineage-appropriate genes to "canalize" developmental progression. Using Oct1 inducible knockout embryonic stem cells, we show that that Oct1 deficiency impairs mesodermal and terminal muscle differentiation in a manner that can be rescued by Oct1 retroviral expression. Using bulk RNA-seq, we show that mesoderm-specific genes are not correctly induced early in the differentiation timecourse. Single-cell gene expression profiling reveals that Oct1-deficient cells lose coherence in temporal induction of lineage programs, and show inappropriate developmental lineage branching resulting in poorly differentiated cells state with epithelial characteristics and hallmarks of oxidative stress. In embryonic stem cells, Oct1 co-binds with Oct4 to genes critical for mesoderm induction. The Utx/Kdm6a histone lysine demethylase also binds to many of these genes, and using a prototypic Pax3 gene we show that Oct1 recruits Utx to remove inhibitory H3K27me3 marks and activate expression. The specificity of the ubiquitous Oct1 protein for mesodermal genes can be explained by cooperative interactions with lineage-driving Smad transcription factors, as we show that Smad and Oct binding sites frequently coexist mesoderm-specific genes, that Oct1 and Smad3 interact, and that the sites and factors act cooperatively at the Myog enhancer. Overall, these results identify Oct1 as a key mediator of the induction of mesoderm lineage-specific genes. Stimulus manipulations permit activation of fiber subpopulations in the mouse and rat vagus Afferent and efferent vagal fibers mediate bidirectional communication between the brain and visceral organs. Small, unmyelinated C-afferents constitute the majority of vagal fibers, play critical roles in numerous interoceptive circuits and autonomic reflexes in health and disease and may contribute to the efficacy and safety of vagus nerve stimulation (VNS). Selective engagement of C-afferents with electrical stimuli has not been feasible, due to the default fiber recruitment order: larger fibers first, smaller fibers last. Here, we determine and optimize an electrical stimulus that selectively engages vagal C-afferents. Intermittent KHz-frequency electrical stimulation (KES) activates motor and, preferentially, sensory vagal neurons in the brainstem. During KES, asynchronous activity of C-afferents increases, while that of larger fibers remains largely unchanged. In parallel, KES effectively blocks excitability of larger fibers while moderately suppressing excitability of C-afferents. By compiling selectivity indices in individual animals, we find that optimal KES parameters for C-afferents are >5KHz frequency and 7-10 times engagement threshold (xT) intensity in rats, 15-25xT in mice. These effects can be explained in computational models by how sodium channel responses to KES are shaped by axonal size and myelin. Our results indicate that selective engagement of vagal C-afferents is attainable by intermittent KES. Spatial and Texture Analysis of Root System Distribution with Earth Mover's Distance (STARSEED) PurposeRoot system architectures are complex, multidimensional, and challenging to characterize effectively for agronomic and ecological discovery.

MethodsWe propose a new method, Spatial and Texture Analysis of Root SystEm distribution with Earth movers Distance (STARSEED), for comparing root system distributions that incorporates spatial information through a novel application of the Earth Movers Distance (EMD).

ResultsWe illustrate that the approach captures the response of sesame root systems for different genotypes and soil moisture levels. STARSEED provides quantitative and visual insights into changes that occur in root architectures across experimental treatments.

ConclusionSTARSEED can be generalized to other plants and provides insight into root system architecture development and response to varying growth conditions not captured by existing root architecture metrics and models. The code and data for our experiments are publicly available: https://github.com/GatorSense/STARSEED. Adipose triglyceride lipase promotes prostaglandin-dependent actin remodeling by regulating substrate release from lipid droplets A key factor controlling oocyte quality and fertility is lipids. Even though lipid droplets (LDs) are crucial regulators of lipid metabolism, their roles in fertility are poorly understood. During Drosophila oogenesis, LD accumulation in nurse cells coincides with dynamic actin remodeling necessary for late-stage follicle morphogenesis and fertility. Loss of the LD-associated Adipose Triglyceride Lipase (ATGL) disrupts both actin bundle formation and cortical actin integrity, an unusual phenotype also seen when Pxt, the enzyme responsible for prostaglandin (PG) synthesis, is missing. Dominant genetic interactions and PG treatment of follicles in vitro reveal that ATGL and Pxt act in the same pathway to regulate actin remodeling, with ATGL upstream of Pxt. Further, lipidomic analysis detects arachidonic acid (AA) containing triglycerides in ovaries. Because AA is the substrate for Pxt, we propose that ATGL releases AA from LDs to drive PG synthesis necessary for follicle development. We also find that exogenous AA is toxic to follicles in vitro, and LDs modulate this toxicity. This leads to the model that LDs both sequester AA to limit toxicity, and release AA via ATGL to drive PG production. We speculate that the same pathways are conserved across organisms to regulate oocyte development and promote fertility. Human pluripotent stem cell-derived cardiomyocytes align under cyclic strain when guided by cardiac fibroblasts The human myocardium is a mechanically active tissue typified by anisotropic organization of the resident cells (cardiomyocytes (CMs) and cardiac fibroblasts (cFBs)) and the extracellular matrix (ECM). Upon ischemic injury, the anisotropic tissue is replaced by disorganized scar tissue, eventually resulting in loss of coordinated contraction. Efforts to re-establish tissue anisotropy in the injured myocardium are hampered by a lack of understanding on how CM and/or cFB structural organization is affected by the two major physical cues inherent in the myocardium: ECM organization and cyclic mechanical strain. Herein, we investigate the singular and combined effect of ECM (dis)organization and cyclic strain in a 2D human in vitro co-culture model of the myocardial microenvironment. We show that (an)isotropic ECM protein patterning can guide the orientation of CMs and cFBs, both in mono- and co-culture. Subsequent application of uniaxial cyclic strain - mimicking the local anisotropic deformation of beating myocardium - causes no effect when applied parallel to the anisotropic ECM. However, when cultured on isotropic substrates, cFBs, but not CMs, orient away from the direction of cyclic uniaxial strain (strain avoidance). In contrast, CMs show strain avoidance via active remodeling of their sarcomeres only when co-cultured with at least 30% cFBs. Paracrine signaling or N-cadherin-mediated communication between CMs and cFBs were no contributing factors, but our findings suggest that the mechanoresponsive cFBs provide structural guidance for CM orientation and elongation. Our study therefore highlights a synergistic mechanobiological interplay between CMs and cFBs in shaping tissue organization, which is of relevance for regenerating functionally organized myocardium. Sablefish (Anoplopoma fimbria Pallas, 1814) plasma biochemistry and hematology reference intervals including blood cell morphology Plasma biochemistry and hematology reference intervals are integral health assessment tools in all medical fields, including aquatic animal health. As sablefish (Anoplopoma fimbria) are becoming aquaculturally and economically more important, this manuscript provides essential reference intervals (RI) for their plasma biochemistry and hematology along with reference photomicrographs of blood cells in healthy, fasted sablefish. Blood cell morphology can differ between fish species. In addition, blood cell counts and blood chemistry can vary between fish species, demographics, water conditions, seasons, diets, and culture systems, which precludes the use of RIs from other fish species. For this study, blood was collected for plasma biochemistry and hematology analysis between June 20 and July 18, 2019, from healthy, yearling sablefish, hatched and reared in captivity on a commercial diet. Overnight fast of 16-18 hours did not sufficiently reduce lipids in the blood, which led to visible lipemia and frequent rupture of blood cells during analysis. Therefore, sablefish should be fasted for 24 to 36 hours before blood is collected to reduce hematology artifacts or possible reagent interference in plasma biochemistry analysis. Lymphocytes were the most dominant leukocytes (98%), while eosinophils were rare, and basophils were not detected in sablefish. Neutrophils were very large cells with Dohle bodies. In mammals and avian species, Dohle bodies are usually signs of toxic change from inflammation, but no such association was found in these fish. In conclusion, lipemia can interfere with sablefish blood analysis, and available removal methods should be evaluated as fasting for up to 36 h might not always be feasible. Also, more studies are required to establish RI for different developmental stages and rearing conditions. Differential impact of mutations in cell cycle checkpoint kinases on breast cancer presentation and prognosis. Cell cycle checkpoint kinases ATM/Chk2 and ATR/Chk1 are important tumor suppressors that are activated by DNA damage and other stressors to arrest replication. Germline variants in ATM and CHEK2 associate with increased incidence of estrogen receptor positive (ER+)/HER2- breast cancer and some evidence suggests that somatic inactivation of ATM/Chk2 induces therapy resistance in ER+/HER2- breast cancer patients. However, the relative contributions of germline and somatic inactivation of each of the four cell cycle checkpoint kinases to ER+/HER2- and ER-/PR-/HER2- (or triple negative breast cancer; TNBC) breast cancer incidence and outcomes has not yet been systematically assessed. Here, we comprehensively analyze associations between mutations in each cell cycle checkpoint kinase gene with age at diagnosis, tumor characteristics, and patient survival outcomes using six breast cancer datasets. Overall, we find that of all the cell cycle checkpoint kinases, CHEK2 (p=0.001) is the only gene that is preferentially mutated in ER+/HER2- breast cancer samples relative to TNBCs. CHEK2 mutations also associate with younger age at diagnosis of metastatic ER+/HER2- breast cancer (median age 41.5 years compared to 50 years for CHEK2 wildtype patients, p=0.03) that is highly PR positive (10% increase in PR positivity in metastatic vs primary CHEK2 mutant tumors, as opposed to a 25% decrease in PR positivity in CHEK2 wildtype tumors, p=0.03). These associations appear to be driven by germline variants, rather than somatic mutations. Whether germline or somatic, deleterious (truncating, frameshift or splice site) mutations in CHEK2 associate with rapid disease progression in response to standard first line endocrine therapy when compared to tumors with wildtype CHEK2 in two independent datasets (HR=6.15, p=0.01). Conversely, neither ATM nor ATR mutations are enriched in ER+/HER2- breast cancers relative to TNBC. However, co-incident mutation of ATR and TP53 is 2-fold enriched in this subset (p=0.002), and further enriched in metastatic (3.3%) vs primary (1%) ER+/HER2- breast cancers (p=0.002). This co-incident mutation of ATR and TP53 associates with decreased PR positivity in primary ER+/HER2- breast cancer diagnoses (p<0.0001) and significantly worse metastasis-free survival than mutations in either ATR or TP53 alone in two independent datasets (HR=2.01, p=0.007). This association is not specific to ER+/HER2- breast cancer but is observable in ER- breast cancer patients as well. These results provide the first systematic analysis of the contribution of germline and somatic mutations in cell cycle checkpoint kinases to tumor characteristics affecting patient prognosis and treatment response. They demonstrate a unique role for CHK2 in ER signaling and response to endocrine therapy, and a more universal role for ATR in combination with TP53 in metastatic progression across breast cancer subtypes. These findings provide rationale for streamlined testing of checkpoint kinase mutations as prognostic markers of metastatic potential and/or predictive markers of responsiveness to endocrine treatment in breast cancer patients. Cellular feedback to organotelluranes displays modulation of antioxidant proteins gene expression Organotelluranes RT3 and RT4 are thiol reagents that induce mitochondrial transition pore (MTP) opening in a sensitive and insensitive manner to cyclosporin A. Although RT3 and RT4 promote glutathione depletion, paradoxically, they are also an efficient antioxidant for membrane lipids. These compounds antagonistic effects elicited the challenging question of how the gene expression of antioxidant enzymes would respond to treatment with these compounds. The influence of RT3 and RT4 on antioxidant enzyme expression was investigated in cultured aortic smooth muscle cells (ASMC). RT3 and RT4 promoted disruption of ionic calcium homeostasis, mitochondrial transmembrane potential ({Delta}{Psi}), and cell death in a dose-dependent manner. The cell death mechanisms responded qualitatively to the increase of the organotellurane concentration and changed from apoptosis to necrosis. RT3 and RT4 increased the expression of thioredoxin significantly. RT3 also increased the expression of glutaredoxin and glutathione peroxidase, slightly the catalase expression without significant effects on SOD expression. The results are consistent with GSH and protein thiol depletion and discussed based on the cell toxicity mechanism exhibited by these compounds.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=150 SRC="FIGDIR/small/425411v4_ufig1.gif" ALT="Figure 1">
View larger version (73K):
org.highwire.dtl.DTLVardef@52bb82org.highwire.dtl.DTLVardef@19476d9org.highwire.dtl.DTLVardef@6aa99org.highwire.dtl.DTLVardef@1aa605e_HPS_FORMAT_FIGEXP  M_FIG C_FIG Nuclear Progestin Receptor Mediated Linkage of Blood Coagulation and Ovulation Ovulation is a dramatic remodeling process that includes rupture of blood capillaries and clotting, but coagulation is not thought to directly regulate this process. Herein, we report remarkable increases of coagulation factors V (f5, [~]3145-fold) and tissue factor (f3a, [~]120-fold) in zebrafish ovarian follicle cells during ovulation. This increase was mediated through the nuclear progestin receptor (Pgr), which is essential for ovulation in zebrafish, and was totally abolished in ovarian follicular cells from pgr-/- mutants. In addition, promoter activities of f5 and f3a were significantly enhanced by progestin (DHP) via Pgr. Similar regulation of human F5 promoter activity was induced via human PGRB, suggesting a conserved mechanism. Site-directed mutagenesis of the zebrafish f5 promoter further demonstrated a direct regulation of coagulation factors via progestin response elements. Moreover, a stark increase of erythrocytes occurred in capillaries meshed in wildtype preovulatory follicles but was absent in pgr-/- mutants. Interestingly, anticoagulants significantly inhibited ovulation both in vitro and in vivo, respectively. Furthermore, reduced fecundity was observed in f5+/- female zebrafish. Taken together, our study provides plausible evidence for steroid regulation of coagulation factors, and a new hypothesis for blood clotting triggered ovulation in vertebrates. Generalization of Cortical MOSTest Genome-Wide Associations Within and Across Samples Genome-Wide Association studies have typically been limited to single phenotypes, given that high dimensional phenotypes incur a large multiple comparisons burden: ~1 million tests across the genome times the number of phenotypes. Recent work demonstrates that a Multivariate Omnibus Statistic Test (MOSTest) is well powered to discover genomic effects distributed across multiple phenotypes. Applied to cortical brain MRI morphology measures, MOSTest has resulted in a drastic improvement in power to discover loci when compared to established approaches (min-P). One question that arises is how well these discovered loci replicate in independent data. Here we perform 10 times cross validation within 35,644 individuals from UK Biobank for imaging measures of cortical area, thickness and sulcal depth (>1,000 dimensionality for each). By deploying a replication method that aggregates discovered effects distributed across multiple phenotypes, termed PolyVertex Score (PVS), we demonstrate a higher replication yield and comparable replication rate of discovered loci for MOSTest (# replicated loci: 348-845, replication rate: 94-95%) in independent data when compared with the established min-P approach (# replicated loci: 31-68, replication rate: 65-80%). An out-of-sample replication of discovered loci was conducted with a sample of 8,336 individuals from the Adolescent Brain Cognitive Development(R) (ABCD) study, who are on average 50 years younger than UK Biobank individuals. We observe a higher replication yield and comparable replication rate of MOSTest compared to min-P. This finding underscores the importance of using well-powered multivariate techniques for both discovery and replication of high dimensional phenotypes in Genome-Wide Association studies. Sgs1BLM independent role of Dna2DNA2 nuclease at DNA double strand break is inhibited by Nej1XLF The two major pathways of DNA double strand break (DSB) repair, non-homologous end-joining (NHEJ) and homologous recombination (HR), are highly conserved from yeast to mammals. The regulation of 5 DNA resection controls repair pathway choice and influences repair outcomes. Nej1 was first identified as a canonical NHEJ factor involved in stimulating the ligation of broken DNA ends and more recently, it was shown to be important for DNA end-bridging and inhibiting 5 resection mediated by Dna2-Sgs1. Nej1 interacts with Sae2 and this impacts DSB repair in three ways. First, Nej1 inhibits MRX-Sae2 interactions and Sae2 localization to a DSB. Second, Nej1 inhibits Sae2-dependent recruitment of Dna2 in the absence of Sgs1. Third, NEJ1 and SAE2 showed an epistatic relationship for DNA end-bridging, an event that restrains the broken ends and reduces the frequency of genomic deletions from developing at the DSB. Deletion of NEJ1 suppressed the synthetic lethality of sae2{Delta} sgs1{Delta} and was dependent on the nuclease activity of Dna2. These Nej1 functions promote end-joining DSB repair, but could also be relevant for controlling resection initiation during HR repair.

HighlightsO_LINej1 physically interacts with Sae2 and inhibits end-resection at a DSB.
C_LIO_LINej1 inhibits Sae2 interactions with the MRX complex.
C_LIO_LINej1 inhibits Sae2-dependent recruitment of Dna2 to a DSB.
C_LIO_LINEJ1 and SAE2 are epistatic for DNA end-bridging.
C_LIO_LIDeletion of NEJ1 suppresses the synthetic lethality of sae2{Delta} sgs1{Delta}, which is dependent on Dna2 nuclease activity.
C_LI Failure to account for behavioral variability significantly compromises accuracy in indirect population monitoring Indirect wildlife population surveying largely depends upon counts of artefacts of behavior (e.g., nests or dung). Likelihood to encounter these artefacts is derived from both artefact production and decay, and variability in artefact production behavior is considered to contribute minimally to inaccuracy in wildlife estimation. Here we demonstrate how ignoring behavioral variability contributes to significant population misestimation, using an example of an endangered ape, the bonobo (Pan paniscus). Until now, a single estimate of (i.e., nest construction) has been used to extrapolate bonobo densities, assumed to be broadly representative of bonobo sign production behavior. We estimated nest construction rate across seasons and social groups at the Kokolopori Bonobo Reserve, DRC and find that nest construction rates in bonobos to be highly variable across populations as well as seasonal. Failure to account for this variability led to potentially severe degradation in accuracy of bonobo population estimates of abundance, accounting for a likely overestimation of bonobo numbers by 34%, and at worst as high as 80%. With this example, we demonstrate that failure to account for inter- and intra-population behavioral variation compromises the ability to estimate both relative and absolute wildlife abundances. We argue that variation in sign production is but one of several potential ways that behavioral variability can affect conservation monitoring, should be measured across contexts whenever possible, and must be considered in population estimation confidence intervals. With increasing attention to behavioral variability as a potential tool for conservation, conservationists must also account for the impact that behavioral variability can play upon wildlife population estimation. Our results underline the importance of observational research to wildlife monitoring schemes as a critical but overlooked component of conservation management. We extensively discuss the various avenues through which behavioral variability is likely to impact wildlife monitoring accuracy and precision and propose potential avenues for accounting for behavioral variability in wildlife monitoring. Evolution of bacterial persistence to antibiotics during a 50,000-generation experiment in an antibiotic-free environment Failure of antibiotic therapies causes > 700,000 deaths yearly and involves both bacterial resistance and persistence. Persistence results in the relapse of infections by producing a tiny fraction of pathogen survivors that stay dormant during antibiotic exposure. From an evolutionary perspective, persistence is either a  bet-hedging strategy that helps to cope with stochastically changing environments or an unavoidable minimal rate of  cellular errors that lock the cells in a low activity state. Here, we analyzed the evolution of persistence over 50,000 bacterial generations in a stable environment by improving a published method that estimates the number of persister cells based on the growth of the reviving population. Our results challenged our understanding of the factors underlying persistence evolution. In one case, we observed a substantial decrease in persistence proportion, suggesting that the naturally observed persistence level is not an unavoidable minimal rate of  cellular errors. However, although there was no obvious environmental stochasticity, in 11 of the 12 investigated populations, the persistence level was maintained during 50,000 bacterial generations. 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. 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/.