Host Chitinase 3-like-1 is a Universal Therapeutic Target for SARS-CoV-2 Viral Variants in COVID 19 COVID 19 is the disease caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2; SC2) which has caused a world-wide pandemic with striking morbidity and mortality. Evaluation of SC2 strains demonstrated impressive genetic variability and many of these viral variants are now defined as variants of concern (VOC) that cause enhanced transmissibility, decreased susceptibility to antibody neutralization or therapeutics and or the ability to induce severe disease. Currently, the delta ({delta}) and omicron (o) variants are particularly problematic based on their impressive and unprecedented transmissibility and ability to cause break through infections. The delta variant also accumulates at high concentrations in host tissues and has caused waves of lethal disease. Because studies from our laboratory have demonstrated that chitinase 3-like-1 (CHI3L1) stimulates ACE2 and Spike (S) priming proteases that mediate SC2 infection, studies were undertaken to determine if interventions that target CHI3L1 are effective inhibitors of SC2 viral variant infection. Here we demonstrate that CHI3L1 augments epithelial cell infection by pseudoviruses that express the alpha, beta, gamma, delta or omicron S proteins and that the CHI3L1 inhibitors anti-CHI3L1 and kasugamycin inhibit epithelial cell infection by these VOC pseudovirus moieties. Thus, CHI3L1 is a universal, VOC-independent therapeutic target in COVID 19.
SARS-CoV-2 variants resist antibody neutralization and broaden host ACE2 usage New SARS-CoV-2 variants continue to emerge from the current global pandemic, some of which can replicate faster and with greater transmissibility and pathogenicity. In particular, UK501Y.V1 identified in UK, SA501Y.V2 in South Africa, and BR501Y.V3 in Brazil are raising serious concerns as they spread quickly and contain spike protein mutations that may facilitate escape from current antibody therapies and vaccine protection. Here, we constructed a panel of 28 SARS-CoV-2 pseudoviruses bearing single or combined mutations found in the spike protein of these three variants, as well as additional nine mutations that within or close by the major antigenic sites in the spike protein identified in the GISAID database. These pseudoviruses were tested against a panel of monoclonal antibodies (mAbs), including some approved for emergency use to treat SARS-CoV-2 infection, and convalescent patient plasma collected early in the pandemic. SA501Y.V2 pseudovirus was the most resistant, in magnitude and breadth, against mAbs and convalescent plasma, followed by BR501Y.V3, and then UK501Y.V1. This resistance hierarchy corresponds with Y144del and 242-244del mutations in the N-terminal domain as well as K417N/T, E484K and N501Y mutations in the receptor binding domain (RBD). Crystal structural analysis of RBD carrying triple K417N-E484K-N501Y mutations found in SA501Y.V2 bound with mAb P2C-1F11 revealed a molecular basis for antibody neutralization and escape. SA501Y.V2 and BR501Y.V3 also acquired substantial ability to use mouse and mink ACE2 for entry. Taken together, our results clearly demonstrate major antigenic shifts and potentially broadening the host range of SA501Y.V2 and BR501Y.V3, which pose serious challenges to our current antibody therapies and vaccine protection.
Modeling microRNA-driven post-transcriptional regulation by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Source-sink connectivity: A novel interictal EEG marker for seizure localization Over 15 million epilepsy patients worldwide have drug-resistant epilepsy (DRE). Successful surgery is a standard of care treatment for DRE but can only be achieved through complete resection or disconnection of the epileptogenic zone (EZ), the brain region(s) where seizures originate. Surgical success rates vary between 20-80% because no clinically validated biological markers of the EZ exist. Localizing the EZ is a costly and time-consuming process beginning with non-invasive neuroimaging and often followed by days to weeks of intracranial EEG (iEEG) monitoring. Clinicians visually inspect iEEG data to identify abnormal activity (e.g., low-voltage high frequency activity) on individual channels occurring immediately before seizures or spikes that occur on interictal iEEG (i.e., between seizures). In the end, the clinical standard mainly relies on a small proportion of the iEEG data captured to assist in EZ localization (minutes of seizure data versus days of recordings), missing opportunities to leverage these largely ignored interictal data to better diagnose and treat patients.

Intracranial EEG offers a unique opportunity to observe epileptic cortical network dynamics but waiting for seizures increases patient risks associated with invasive monitoring. In this study, we aim to leverage interictal iEEG data by developing a new network-based interictal iEEG marker of the EZ. We hypothesize that when a patient is not clinically seizing, it is because the EZ is inhibited by other regions. We developed an algorithm that identifies two groups of nodes from the interictal iEEG network: those that are continuously inhibiting a set of neighboring nodes ("sources") and the inhibited nodes themselves ("sinks"). Specifically, patient-specific dynamical network models (DNMs) were estimated from minutes of iEEG and their connectivity properties revealed top sources and sinks in the network, with each node being quantified by source-sink metrics (SSMs). We validated the SSMs in a retrospective analysis of 65 patients by using the SSMs of the annotated EZ to predict surgical outcomes. The SSMs predicted outcomes with an accuracy of 79% compared to an accuracy of 43% for clinicians predictions (surgical success rate of this dataset). In failed outcomes, we identified regions of the brain with high SSMs that were untreated. When compared to high frequency oscillations, the most commonly proposed interictal iEEG feature for EZ localization, SSMs outperformed in predictive power (by a factor of 1.2) suggesting SSMs may be an interictal iEEG fingerprint of the EZ. Behavioral and postural analyses establish sleep-like states for mosquitoes that can impact host landing and blood feeding Sleep is an evolutionarily conserved process that has been described in different animal systems. For insects, sleep characterization has been primarily achieved using behavioral and electrophysiological correlates in a few systems. Sleep in mosquitoes, which are important vectors of disease-causing pathogens, has not been directly examined. This is surprising as circadian rhythms, which have been well studied in mosquitoes, influence sleep in other systems. In this study, we characterized sleep in mosquitoes using body posture analysis and behavioral correlates, and quantified the effect of sleep deprivation on sleep rebound, host landing and blood-feeding propensity. Body and appendage position metrics revealed a clear distinction between the posture of mosquitoes in their putative sleep and awake states for multiple species, which correlate with a reduction in responsiveness to host cues. Sleep assessment informed by these posture analyses indicated significantly more sleep during periods of low activity. Nighttime and daytime sleep deprivation resulting from the delivery of vibration stimuli induced sleep rebound in the subsequent phase in day and night active mosquitoes, respectively. Lastly, sleep deprivation suppressed host landing in both laboratory and field settings and also impaired blood feeding of a human host when mosquitoes would normally be active. These results suggest that quantifiable sleep states occur in mosquitoes, and highlight the potential epidemiological importance of mosquito sleep. Assessment and Optimization of the Interpretability of Machine Learning Models Applied to Transcriptomic Data Explainable artificial intelligence aims to interpret how the machine learning models make decisions, and many model explainers have been developed in the computer vision field. However, the understandings of the applicability of these model explainers to biological data are still lacking. In this study, we comprehensively evaluated multiple explainers by interpreting pretrained models of predicting tissue types from transcriptomic data, and by identifying top contributing genes from each sample with the greatest impacts on model prediction. To improve the reproducibility and interpretability of results generated by model explainers, we proposed a series of optimization strategies for each explainer on two different model architectures of Multilayer Perceptron (MLP) and Convolutional Neural Network (CNN). We observed three groups of explainer and model architecture combinations with high reproducibility. Group II, which contains three model explainers on aggregated MLP models, identified top contributing genes in different tissues that exhibited tissue-specific manifestation and were potential cancer biomarkers. In summary, our work provides novel insights and guidance for exploring biological mechanisms using explainable machine learning models. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Chemotherapy coupled to macrophage inhibition leads to enhanced T and B cell infiltration and durable triple negative breast cancer regression Immunosuppressive elements within the tumor microenvironment such as Tumor Associated Macrophages (TAMs) can present a barrier to successful anti-tumor responses by cytolytic T cells. We employed preclinical syngeneic p53 null mouse models of triple negative breast cancer (TNBC) to develop a treatment regimen that harnessed the immunostimulatory effects of low-dose cyclophosphamide coupled with the pharmacologic inhibition of TAMs using either a small molecule CSF1R inhibitor or an anti-CSF1R antibody. This therapeutic combination was used to successfully treat several highly aggressive TNBC murine mammary tumors and lung metastasis. Using this regimen and single cell RNA sequencing we characterized tumor infiltrating lymphocytes (TILs) including helper T cells and antigen-presenting B cells that were highly enriched in good responders to combination therapy. Using high dimensional imaging techniques, we identified the close spatial localization of B220+ CD86+ activated B cells and CD4+ T cells in tertiary lymphoid structures that were present up to 6 weeks post-treatment in one model that also exhibited long-term tumor regression post-treatment. We also characterized the transcriptional and metabolic heterogeneity of TAMS in these two closely related claudin-low/mesenchymal subtype tumor models with differential treatment responses. A murine TAM signature derived from the T12 model is highly expressed and conserved in human claudin-low breast cancers, and high expression of the T12 signature correlated with reduced overall survival. This T12 tumor TAM signature may help identify human claudin-low breast cancer patients that will benefit from the combination of cyclophosphamide and anti-CSF1R therapy. These studies illustrate the complexity of the tumor immune microenvironment and highlight different immune responses that result from rationale combinations of immunotherapy.

SignificanceA treatment regimen that harnessed the immunostimulatory effects of cyclophosphamide coupled with the inhibition of CSF1R was used to successfully treat several highly aggressive claudin-low TNBC murine mammary tumors and lung metastasis. The relative importance of individual attributes varies across preference elicitation methods When choosing between different options, we tend to consider specific attribute qualities rather than deliberating over some general sense of the options overall values. The importance of each attribute together with its quality will determine our preference rankings over the available alternatives. Here, we test the hypothesis that the most prominent class of model for simple decisions - sequential sampling or evidence accumulation to bound - can be bolstered by explicitly incorporating variables related to individual attributes in addition to the standard usage of overall value estimates. We examine six datasets in which participants evaluated snack foods both in terms of overall value and individual attributes, then chose between pairs of the same snacks, and show that only models that explicitly incorporate information about the individual attributes are able to reproduce fundamental patterns in the choice data, such as the influence of attribute disparity on decisions, and such models provide quantitatively better fits to the choice outcomes, response times, and confidence ratings compared to models based on overall value alone. Our results provide important evidence that incorporating attribute-level information into computational models helps us to better understand the cognitive processes involved in value-based decision-making. Trehalose increases tomato drought tolerance, induces defenses, and increases resistance to bacterial wilt disease Ralstonia solanacearum causes bacterial wilt disease, leading to severe crop losses. Xylem sap from R. solanacearum-infected tomato is enriched in the disaccharide trehalose. Water-stressed plants also accumulate trehalose, which increases drought tolerance via abscisic acid (ABA) signaling. Because R. solanacearum-infected plants suffer reduced water flow, we hypothesized that bacterial wilt physiologically mimics drought stress, which trehalose could mitigate. We found that R. solanacearum-infected plants differentially expressed drought-associated genes, including those involved in ABA and trehalose metabolism, and had more ABA in xylem sap. Consistent with this, treating tomato roots with ABA reduced both stomatal conductance and stem colonization by R. solanacearum. Treating roots with trehalose increased xylem sap ABA and reduced plant water use by lowering stomatal conductance and temporarily improving water use efficiency. Trehalose treatment also upregulated expression of salicylic acid (SA)-dependent tomato defense genes; increased xylem sap levels of SA and other antimicrobial compounds; and increased bacterial wilt resistance of SA-insensitive NahG tomato plants. Additionally, trehalose treatment increased xylem concentrations of jasmonic acid and related oxylipins. Finally, trehalose-treated plants were substantially more resistant to bacterial wilt disease. Together, these data show that exogenous trehalose reduced both water stress and bacterial wilt disease and triggered systemic disease resistance, possibly through a Damage Associated Molecular Pattern (DAMP) response pathway. This suite of responses revealed unexpected linkages between plant responses to biotic and abiotic stress and suggested that R. solanacearum-infected plants increase trehalose to improve water use efficiency and increase wilt disease resistance. The pathogen may degrade trehalose to counter these efforts. Together, these results suggest that treating tomatoes with exogenous trehalose could be a practical strategy for bacterial wilt management. Bulk-assembly of monodisperse coacervates and giant unilamellar vesicles with programmable hierarchical complexity Compartmentalization is an essential step for the emergence of protocellular life on the early earth and the endeavor to bottom-up construct synthetic cells in the lab. Among manifold natural strategies for compartmentalization, assembly of lipid vesicles or coacervate droplets with resemblance to cell envelope or cytoplasm have emerged as two dominant paradigms. However, the spontaneous assembly of lipid vesicles or coacervates generally results in polydisperse microcompartments and lacks the ability to efficiently integrate functional building blocks towards defined complexity, which give rise to the difficulty to achieve the order and complexity towards the scope of that of life via spontaneous assembly. Herein, we show that the interplay of coacervates and colloidal particles provides a bulk-assembly approach to form monodisperse coacervates and giant unilamellar vesicles (GUVs) with programmable hierarchical complexity. The hierarchical complexity of the synthetic cell microcompartments can be delicately engineered from the control of the spatial or temporal organization of components in entity unit to manifold forms of interconnected synthetic cell consortia, achieving a compartment hierarchy analogous to the organelle-cell-tissue structure. The increase in complexity gives rise to emergent properties, for instance, collective morphology evolution of coacervate assemblages, gated permeability of GUVs without sophisticated protein machinery, and remarkable structural and functional stability of GUVs at extreme conditions. This work paves an unprecedented step to form monodisperse, hierarchical coacervates and GUVs via bulk-assembly, which provides a prebiotically plausible approach for the proper integration of inanimate matter towards the order and complexity of life and monodisperse, hierarchical microdroplets for follow-up scientific and technological applications. 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. Effects of phenotypic variation on consumer coexistence and prey community structure A popular idea in ecology is that trait variation among individuals from the same species may promote the coexistence of competing species. However, theoretical and empirical tests of this idea have yielded inconsistent findings. We manipulated intraspecific trait diversity in a ciliate competing with a nematode for bacterial prey in experimental microcosms. We found that intraspecific trait variation inverted the original competitive hierarchy to favor the consumer with variable traits, ultimately resulting in competitive exclusion. This competitive outcome was driven by foraging traits (size, speed, and directionality) that increased the ciliates fitness ratio and niche overlap with the nematode. The interplay between consumer trait variation and competition resulted in non-additive cascading effects - mediated through prey defense traits - on prey community assembly. Our results suggest that predicting consumer competitive population dynamics and the assembly of prey communities will require understanding the complexities of trait variation within consumer species. Bacterial community dynamics during embryonic development of the little skate (Leucoraja erinacea) BackgroundMicrobial transmission from parent to offspring is hypothesized to be widespread in vertebrates. However, evidence for this is limited as many evolutionarily important clades remain unexamined. There is currently no data on the microbiota associated with any Chondrichthyan species during embryonic development, despite the global distribution, ecological importance, and phylogenetic position of this clade. In this study, we take the first steps towards filling this gap by investigating the microbiota associated with embryonic development in the little skate, Leucoraja erinacea, a common North Atlantic species and popular system for chondrichthyan biology.

MethodsTo assess the potential for bacterial transmission in an oviparous chondrichthyan, we used 16S rRNA amplicon sequencing to characterize the microbial communities associated with the skin, gill, and egg capsule of the little skate, at six points during ontogeny. Community composition was analyzed using the QIIME2 pipeline and microbial continuity between stages was tracked using FEAST.

ResultsWe identify site-specific and stage-specific microbiota dominated by the bacterial phyla Proteobacteria and Bacteroidetes. This composition is similar to, but distinct from, that of previously published data on the adult microbiota of other chondrichthyan species. Our data reveal that the skate egg capsule harbors a highly diverse bacterial community-particularly on the internal surface of the capsule-and facilitates intergenerational microbial transfer to the offspring. Embryonic skin and external gill tissues host similar bacterial communities; the skin and gill communities later diverge as the internal gills and skin denticles develop.

ConclusionsOur study is the first exploration of the chondrichthyan microbiota throughout ontogeny and provides the first evidence of vertical transmission in this group. Uncovering non-random binary patterns within sequences of intrinsically disordered proteins Sequence-ensemble relationships of intrinsically disordered proteins (IDPs) are governed by binary patterns such as the linear clustering or mixing of specific residues or residue types with respect to one another. To enable the discovery of potentially important, shared patterns across sequence families, we describe a computational method referred to as NARDINI for Non-random Arrangement of Residues in Disordered Regions Inferred using Numerical Intermixing. This work was partially motivated by the observation that parameters that are currently in use for describing different binary patterns are not interoperable across IDPs of different amino acid compositions and lengths. In NARDINI, we generate an ensemble of scrambled sequences to set up a composition-specific null model for the patterning parameters of interest. We then compute a series of pattern-specific z-scores to quantify how each pattern deviates from a null model for the IDP of interest. The z-scores help in identifying putative non-random linear sequence patterns within an IDP. We demonstrate the use of NARDINI derived z-scores by identifying sequence patterns in three well-studied IDP systems. We also demonstrate how NARDINI can be deployed to study archetypal IDPs across homologs and orthologs. Overall, NARDINI is likely to aid in designing novel IDPs with a view toward engineering new sequence-function relationships or uncovering cryptic ones. We further propose that the z-scores introduced here are likely to be useful for theoretical and computational descriptions of sequence-ensemble relationships across IDPs of different compositions and lengths. Social integration predicts survival in female white-faced capuchin monkeys Across multiple species of social mammals, a growing number of studies have found that that individual sociality is associated with survival. In long-lived species, like primates, lifespan is one of the main components of fitness. We used 18 years of data from the Lomas Barbudal Monkey project to quantify social integration in 11 capuchin (Cebus capucinus) groups and tested whether female survivorship was associated with females tendencies to interact with three types of partners: (1) all group members, (2) adult females, (3) adult males. We found strong evidence that females who engaged more with other females in affiliative interactions and foraged in close proximity experienced increased survivorship. We found some weak evidence that females might also benefit from engaging in more support in agonistic contexts with other females. These benefits were evident in models that account for the females rank and group size. Female interactions with all group members also increased survival, but the estimates of the effects were more uncertain. In interactions with adult males, only females who provided more grooming to males survived longer. The results presented here suggest that social integration may result in survival-related benefits. Females might enjoy these benefits through exchanging grooming for other currencies, such as coalitionary support or tolerance. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation. Non-targeted multimodal metabolomics and lipidomics data from ovine rumen fluid fractions BackgroundMetabolomics is a powerful and sensitive approach for investigating low molecular weight metabolite profiles present in rumen biofluids to identify potential roles of metabolites in rumen microbiome and understanding host-level regulatory mechanisms associated with animal production.

FindingsRumen samples from sheep grazed on a mixed ryegrass and clover pasture diet, were fractionated based on molecular weight and analysed using metabolomics and lipidomics to detect the small molecules present in ovine rumen fluid fractions.

ConclusionsUntargeted metabolomics provides a detailed snapshot of the ovine ruminal fluid metabolome that can be used as reference for future studies on ovine rumen fluid or as a comparator for other ruminant species. All data and metadata are available for download in the MetaboLights database. Determining significant correlation between pairs of extant characters in a small parsimony framework When studying the evolutionary relationships among a set of species, the principle of parsimony states that a relationship involving the fewest number of evolutionary events is likely the correct one. Due to its simplicity, this principle was formalized in the context of computational evolutionary biology decades ago by, e.g., Fitch and Sankoff. Because the parsimony framework does not require a model of evolution, unlike maximum likelihood or Bayesian approaches, it is often a good starting point when no reasonable estimate of such a model is available.

In this work, we devise a method for determining if pairs of discrete characters are significantly correlated across all most parsimonious reconstructions, given a set of species on these characters, and an evolutionary tree. The first step of this method is to use Sankoffs algorithm to compute all most parsimonious assignments of ancestral states (of each character) to the internal nodes of the phylogeny. Correlation between a pair of evolutionary events (e.g., absent to present) for a pair of characters is then determined by the (co-) occurrence patterns between the sets of their respective ancestral assignments. The probability of obtaining a correlation this extreme (or more) under a null hypothesis where the events happen randomly on the evolutionary tree is then used to assess the significance of this correlation. We implement this method: parcours (PARsimonious CO-occURrenceS) and use it to identify significantly correlated evolution among vocalizations and morphological characters in the Felidae family.

The parcours tool is freely available at https://github.com/murraypatterson/parcours Mitofusin 1 and 2 regulation of mitochondrial DNA content is a critical determinant of glucose homeostasis The dynamin-like GTPases Mitofusin 1 and 2 (Mfn1 and Mfn2) are essential for mitochondrial function, which has been principally attributed to their regulation of fission/fusion dynamics. Here, we report that Mfn1 and 2 are critical for glucose-stimulated insulin secretion (GSIS) primarily through control of mitochondrial DNA (mtDNA) content. Whereas Mfn1 and Mfn2 individually were dispensable for glucose homeostasis, combined Mfn1/2 deletion in {beta}-cells reduced mtDNA content, impaired mitochondrial morphology and networking, and decreased respiratory function, ultimately resulting in severe glucose intolerance. Importantly, gene dosage studies unexpectedly revealed that Mfn1/2 control of glucose homeostasis was dependent on maintenance of mtDNA content, rather than mitochondrial structure. Mfn1/2 maintain mtDNA content by regulating the expression of the crucial mitochondrial transcription factor Tfam, as Tfam overexpression ameliorated the reduction in mtDNA content and GSIS in Mfn1/2-deficient {beta}-cells. Thus, the primary physiologic role of Mfn1 and 2 in {beta}-cells is coupled to the preservation of mtDNA content rather than mitochondrial architecture, and Mfn1 and 2 may be promising targets to overcome mitochondrial dysfunction and restore glucose control in diabetes. Cell Rearrangement Generates Pattern Emergence as a Function of Temporal Morphogen Exposure. As tissues elongate, cell rearrangement alters positional information in manner that must be accounted for to generate gene expression pattern. How this is achieved during paraxial mesoderm elongation is unknown. By reverse-engineering gene regulatory networks that predict single cell expression trajectories across the tissue, we find a network capable of recapitulating the full range of dynamic differentiation profiles observed both in vivo and in vitro. Simulating gene expression profiles on in toto cell tracking data sets reveal that temporal exposure to Wnt and FGF is generated by cell movement. The absence of reversal in gene expression towards a more premature gene expression state predicts the generation of aberrant tbx6 expression in the posterior progenitor zone that we then confirm by quantitative single cell imaging. Taken together, these results demonstrate cell rearrangement tunes the dynamics of mesoderm progenitor differentiation to generate pattern emergence as a function of temporal Wnt and FGF exposure. Inhibitors of ERp44, PDIA1, and AGR2 induce disulfide-mediated oligomerization of Death Receptors 4 and 5 and cancer cell death Breast cancer mortality remains unacceptably high, indicating a need for safer and more effective therapeutic agents. Disulfide bond Disrupting Agents (DDAs) were previously identified as a novel class of anticancer compounds that selectively kill cancers that overexpress the Epidermal Growth Factor Receptor (EGFR) or its family member HER2. DDAs kill EGFR+ and HER2+ cancer cells via the parallel downregulation of EGFR, HER2, and HER3 and activation/oligomerization of Death Receptors 4 and 5 (DR4/5). However, the mechanisms by which DDAs mediate these effects are unknown. Affinity purification analyses employing biotinylated-DDAs reveal that the Protein Disulfide Isomerase (PDI) family members AGR2, PDIA1, and ERp44 are DDA target proteins. Further analyses demonstrate that shRNA-mediated knockdown of AGR2 and ERp44, or expression of ERp44 mutants, enhance basal and DDA-induced DR5 oligomerization. DDA treatment of breast cancer cells disrupts PDIA1 and ERp44 mixed disulfide bonds with their client proteins. Together, the results herein reveal DDAs as the first small molecule, active site inhibitors of AGR2 and ERp44, and demonstrate roles for AGR2 and ERp44 in regulating the activity, stability, and localization of DR4 and DR5, and activation of Caspase 8. ATP-binding cassette protein ABCF1 couples gene transcription with maintenance of genome integrity in embryonic stem cells OCT4 and SOX2 confer pluripotency by recruiting coactivators to activate stem cell-specific gene expression programs. However, the composition of coactivator complexes and their roles in maintaining stem cell fidelity remain unclear. Here we report the identification of ATP-binding cassette subfamily F member 1 (ABCF1) as a critical coactivator for OCT4/SOX2. ABCF1 is required for pluripotency gene expression and stem cell self-renewal. ABCF1 binds co-dependent coactivators XPC and DKC1 via its intrinsically disordered region and stimulates transcription by linking SOX2 to the transcription machinery. Furthermore, in response to pathogen infection and DNA damage, ABCF1 binds intracellular DNAs accumulated in cells, concomitant with loss of SOX2 interaction and pluripotency gene transcription. This results in spontaneous differentiation of compromised stem cells and elimination from the self-renewing population. Thus, ABCF1 directly couples pluripotency gene transcription with sensing aberrant DNAs and acts as a checkpoint for self-renewal to safeguard stem cell fidelity and genome integrity. Multi-omics-based label-free metabolic flux inference reveals obesity-associated dysregulatory mechanisms in liver glucose metabolism Glucose homeostasis is maintained by modulation of metabolic flux. Enzymes and metabolites regulate the involved metabolic pathways. Dysregulation of glucose homeostasis is a pathological event in obesity. Analyzing metabolic pathways and the mechanisms contributing to obesity-associated dysregulation in vivo is challenging. Here, we introduce OMELET: Omics-Based Metabolic Flux Estimation without Labeling for Extended Trans-omic Analysis. OMELET uses metabolomic, proteomic, and transcriptomic data to identify changes in metabolic flux, and to quantify contributions of metabolites, enzymes, and transcripts to the changes in metabolic flux. By evaluating the livers of fasting ob/ob mice, we found that increased metabolic flux through gluconeogenesis resulted primarily from increased transcripts, whereas that through the pyruvate cycle resulted from both increased transcripts and changes in substrates of metabolic enzymes. With OMELET, we identified mechanisms underlying the obesity-associated dysregulation of metabolic flux in liver.

HighlightsO_LIWe created OMELET to infer metabolic flux and its regulation from multi-omic data.
C_LIO_LIGluconeogenic and pyruvate cycle fluxes increased in fasting ob/ob mice.
C_LIO_LITranscripts increases mediated the increase in gluconeogenic fluxes in ob/ob mice.
C_LIO_LIIncreases in transcripts and substrates enhanced pyruvate cycle flux in ob/ob mice.
C_LI Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. wnt16 exerts pleiotropic effects on bone and lean mass in zebrafish Osteoporosis, a disease of bone fragility, and sarcopenia, a condition of reduced muscle mass and strength, commonly occur in the same individual--a condition termed osteosarcopenia. Genetic variants at the CPED1-WNT16 locus have pleiotropic effects on human bone mineral density (BMD) and lean tissue mass. However, our current understanding of Wnt16 and its influence on bone fails to explain effects on lean mass at this locus. In this study, we sought to better understand the biological mechanism underlying pleiotropy at the CPED1-WNT16 locus on BMD and lean mass. Using single cell analysis, microCT imaging, and genetic approaches in zebrafish, we show that wnt16 exerts dual influence on muscle and bone, and provide evidence that this occurs through a function in somite development, a period in which precursors for both of these tissues are specified. We show that mutations in wnt16 alter vertebral size and lean mass during axial skeletogenesis. We show that wnt16+ cells transcriptionally resemble cells within the amniote dermomyotome, and provide evidence they function as myogenic precursors during embryonic development. We also provide evidence that WNT16 is a gene of major effect on lean mass at the CPED-WNT16 locus. Our studies support the potential for wnt16 to exert pleiotropic effects on bone and lean mass through a function in myogenic precursors.