DeLUCS: Deep Learning for Unsupervised Clustering of DNA Sequences We present a novel Deep Learning method for the Unsupervised Clustering of DNA Sequences (DeLUCS) that does not require sequence alignment, sequence homology, or (taxonomic) identifiers. DeLUCS uses Frequency Chaos Game Representations (FCGR) of primary DNA sequences, and generates "mimic" sequence FCGRs to self-learn data patterns (genomic signatures) through the optimization of multiple neural networks. A majority voting scheme is then used to determine the final cluster assignment for each sequence. The clusters learned by DeLUCS match true taxonomic groups for large and diverse datasets, with accuracies ranging from 77% to 100%: 2,500 complete vertebrate mitochondrial genomes, at taxonomic levels from sub-phylum to genera; 3,200 randomly selected 400 kbp-long bacterial genome segments, into clusters corresponding to bacterial families; three viral genome and gene datasets, averaging 1,300 sequences each, into clusters corresponding to virus subtypes. DeLUCS significantly outperforms two classic clustering methods (K-means++ and Gaussian Mixture Models) for unlabelled data, by as much as 47%. DeLUCS is highly effective, it is able to cluster datasets of unlabelled primary DNA sequences totalling over 1 billion bp of data, and it bypasses common limitations to classification resulting from the lack of sequence homology, variation in sequence length, and the absence or instability of sequence annotations and taxonomic identifiers. Thus, DeLUCS offers fast and accurate DNA sequence clustering for previously intractable datasets.
Ultrafast, alignment-free detection of repeat expansions in next-generation DNA and RNA sequencing data Short tandem repeat expansions are an established cause of diseases such as Huntingtons disease. Bioinformatic methods for detecting repeat expansions in short-read sequencing have revealed new repeat expansions in humans. Current bioinformatic methods to detect repeat expansions require alignment information to identify repetitive motif enrichment at genomic locations. We present superSTR, an ultrafast method that does not require alignment. We demonstrate superSTRs ability to efficiently process both whole-genome and whole-exome sequencing data. Using superSTR we perform the first analysis of the UK Biobank to efficiently screen the exomes of 49,953 biobank participants for repeat expansions. We identify known mutations, as well as diseases not previously associated with REs. We further demonstrate the first bioinformatic screening of RNA sequencing data to detect repeat expansions in patients with spinocerebellar ataxia and Fuchs endothelial corneal dystrophy, and mouse models of myotonic dystrophy. superSTR is a highly computationally-efficient repeat expansion tool screening and detection tool for genomewide novel repeat expansion analysis, significantly outperforming existing methods. superSTR is available from https://github.com/bahlolab/superSTR.
Cell type-specific differential expression in spatial transcriptomics Spatial transcriptomics enables spatially resolved gene expression measurements at near single-cell resolution. There is a pressing need for computational tools to enable the detection of genes that are differentially expressed (DE) within specific cell types across tissue context. We show that current approaches cannot learn cell type-specific DE due to changes in cell type composition across space and the fact that measurement units often detect transcripts from more than one cell type. Here, we introduce a statistical method, Cell type-Specific Inference of Differential Expression (C-SIDE), that identifies cell type-specific patterns of differential gene expression while accounting for localization of other cell types. We model spatial transcriptomics gene expression as an additive mixture across cell types of general log-linear cell type-specific expression functions. This approach provides a unified framework for defining and identifying gene expression changes in a wide-range of relevant contexts: changes due to pathology, anatomical regions, physical proximity to specific cell types, and cellular microenvironment. Furthermore, our approach enables statistical inference across multiple samples and replicates when such data is available. We demonstrate, through simulations and validation experiments on Slide-seq and MER-FISH datasets, that our approach accurately identifies cell type-specific differential gene expression and provides valid uncertainty quantification. Lastly, we apply our method to characterize spatially-localized tissue changes in the context of disease. In an Alzheimers mouse model Slide-seq dataset, we identify plaque-dependent patterns of cellular immune activity. We also find a putative interaction between tumor cells and myeloid immune cells in a Slide-seq tumor dataset. We make our C-SIDE method publicly available as part of the open source R package https://github.com/dmcable/spacexr. Chromosome counting in the mouse zygote using low-invasive super-resolution live-cell imaging In preimplantation embryos, an abnormal chromosome number causes developmental failure and a reduction in the pregnancy rate. Conventional chromosome testing methods requiring biopsy reduce the risk of associated genetic diseases; nevertheless, the reduction in cell number also reduces the pregnancy rate. Therefore, we attempted to count the chromosomes in mouse (Slc:ICR) embryos using super-resolution live-cell imaging as a new method of chromosome counting that does not reduce the cell number or viability. We counted the forty chromosomes at the first mitosis by injecting embryos with histone H2B-mCherry mRNA under conditions by which pups could be obtained; however, the results were often an underestimation of chromosome number and varied by embryo and time point. Therefore, we developed a method to count the chromosomes via CRISPR/dCas-mediated live-cell fluorescence in situ hybridization targeting the sequence of the centromere region, enabling us to count the chromosomes more accurately in mouse embryos. The methodology presented here may provide useful information for assisted reproductive technologies, such as those used in livestock animals/humans, as a technique for assessing the chromosomal integrity of embryos prior to transfer.

Impact StatementLow-toxic super-resolution observation enables chromosome counting in preimplantation embryos without cell collection. Delayed induction of type I and III interferons mediates nasal epithelial cell permissiveness to SARS-CoV-2 The nasal epithelium is a plausible entry point for SARS-CoV-2, a site of pathogenesis and transmission, and may initiate the host response to SARS-CoV-2. Antiviral interferon (IFN) responses are critical to outcome of SARS-CoV-2. Yet little is known about the interaction between SARS-CoV-2 and innate immunity in this tissue. Here we applied single-cell RNA sequencing and proteomics to a primary cell model of human nasal epithelium differentiated at air-liquid interface. SARS-CoV-2 demonstrated widespread tropism for nasal epithelial cell types. The host response was dominated by type I and III IFNs and interferon-stimulated gene products. This response was notably delayed in onset relative to viral gene expression and compared to other respiratory viruses. Nevertheless, once established, the paracrine IFN response began to impact on SARS-CoV-2 replication. When provided prior to infection, recombinant IFN{beta} or IFN{lambda}1 induced an efficient antiviral state that potently restricted SARS-CoV-2 viral replication, preserving epithelial barrier integrity. These data suggest that the IFN-I/III response to SARS-CoV-2 initiates in the nasal airway and suggest nasal delivery of recombinant IFNs to be a potential chemoprophylactic strategy. Erythropoietin supplementation ameliorates the immunological and hematological deficiencies of lysinuric protein intolerance in mice SLC7A7 encodes for y+LAT1, a cationic amino acid transporter expressed on the basolateral membrane of epithelial cells and macrophages. Mutations in human SLC7A7 cause lysinuric protein intolerance (LPI), a rare autosomal disorder with a variety of clinical symptoms. Citrulline and a low protein diet are widely used to treat the main symptoms of LPI; however, the hematological and immunological abnormalities in LPI remain poorly understood and difficult to manage by established treatments. Here, we show that erythropoietin levels are significantly diminished in mice lacking Slc7a7, which reduces the number of erythroblast precursors, resulting in mature red blood cells with a smaller mean corpuscular volume. Ex vivo analysis demonstrated that erythrocytes from Slc7a7-deficient mice were phagocytosed more by bone marrow-derived macrophages. Gene expression analysis of sorted red pulp macrophages (RPMs) from Slc7a7-deficient mice indicated that this cell population had higher erythrophagocytic activity, but preserved heme catabolic activity. The increase in erythrophagocytosis resulted in iron accumulation in macrophages and tissues, and ultimately led to RPM apoptosis. Surprisingly, ablation of Slc7a7 in myeloid cells failed to trigger tissue iron accumulation or changes in the RPM population, indicating that the expression of this gene in macrophages is not key for the development of iron accumulation in tissue. Finally, restoration of serum erythropoietin levels in Slc7a7-deficient mice re-established normal erythropoiesis and erythrophagocytosis and prevented the accumulation of iron in tissues. Thus, erythropoietin has a fundamental role in the immune-hematological alterations of LPI in mice, revealing a new crucial role for this hormone in this disease. Heterogeneous recruitment abilities to RNA polymerases generate nonlinear scaling of gene expression with cell volume While most genes expression levels are proportional to cell volumes, some genes exhibit nonlinear scaling between their expression levels and cell volume. Therefore, their mRNA and protein concentrations change as the cell volume increases, which often have crucial biological functions such as cell-cycle regulation. However, the biophysical mechanism underlying the nonlinear scaling between gene expression and cell volume is still unclear. In this work, we show that the nonlinear scaling is a direct consequence of the heterogeneous recruitment abilities of promoters to RNA polymerases based on a gene expression model at the whole-cell level. Those genes with weaker (stronger) recruitment abilities than the average ability spontaneously exhibit superlinear (sublinear) scaling with cell volume. Analysis of the promoter sequences and the nonlinear scaling of Saccharomyces cerevisiaes mRNA levels shows that motifs associated with transcription regulation are indeed enriched in genes exhibiting nonlinear scaling, in concert with our model. Re-convolving the compositional landscape of primary and recurrent glioblastoma reveals prognostic andtargetable tissue states Glioblastoma is an aggressive diffusely infiltrating neoplasm that spreads beyond surgical resection margins, where it intermingles with non-neoplastic brain cells. This complex microenvironment harbouring infiltrating glioma and non-neoplastic brain cells is the origin of tumor recurrence. Thus, understanding the cellular and molecular features of the glioma microenvironment is therapeutically and prognostically important. Here, we used single-nucleus RNA sequencing (snRNAseq) of primary and recurrent glioma to define three compositional  tissue-states rooted in cohabitation of cell-types and transcriptional states. These comprise a state featuring A) abundant normal brain microenvironment cells, B) reactive/inflamed brain microenvironment, and C) cellular/proliferative tumor. All these states exhibited variable degrees of infiltration by glioma cells. Spatial transcriptomics confirmed that the cell-types and transcriptomics states which compositionally cohabitate tissue states also colocalize in space. Tissue states are clinically significant because they correlate with radiographic, histopathologic, and prognostic features. We found that enrichment of tissue state B signature correlated with shorter survival. Importantly, we found that our compositionally defined tissue states are enriched in distinct metabolic pathways. One such pathway is fatty acid biosynthesis, which was enriched in tissue state B - a state enriched in recurrent glioblastoma and composed of astrocyte-like/mesenchymal glioma cells, reactive astrocytes resembling those seen in neurodegeneration, and monocyte-like myeloid cells. We showed that treating acute slices of GBM with a fatty acid synthesis inhibitor is sufficient to deplete the transcriptional signature of tissue state B. Our findings define a novel compositional approach to the glioma-infiltrated tissue which allows us to discover prognostic and targetable features, paving the way to new mechanistic and therapeutic discoveries. Fronto-Parietal Network Segregation Predicts Maintained Processing Speed in the Cognitively Healthy Oldest-old Evaluating brain network interactions in cognitively healthy older adults can help us understand how brain characteristics vary with age and how these variations affect cognitive functioning. Functional connections among groups of brain areas give insight into the brains organization, and the cognitive effects of aging may relate to this large-scale organization. We investigated functional network properties in 146 cognitively healthy participants aged 85+ in the McKnight Brain Aging Registry. We found that the segregation of both the cortical association system and the segregation of the fronto-parietal network (FPN) were strong predictors of cognition and processing speed. We also provide a healthy oldest-old (85+) cortical parcellation that can be used in future work in this age group. This study shows that network segregation of the oldest-old brain is closely linked with cognitive performance. In particular, segregation of the FPN plays an important role in supporting overall cognition and processing speed in this aging cohort. This work adds to the growing body of knowledge about age-related dedifferentiation by demonstrating that cognitive ability is associated with differentiated functional networks even in very old individuals experiencing successful cognitive aging. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. Differential fate of acellular vascular scaffolds in vivo involves macrophages and T-cell subsets Biological scaffold or implant is a popular choice for the preparation of tissue-engineered organs and has the potential to address donor shortage in the clinics. However, biological scaffolds prepared by physical or chemical agents cause damage to the extracellular matrix by potentially inducing immune responses after implantation. The current study explores an alternative route for the preparation of acellular scaffolds and explores the fate of the prepared scaffolds in a milieu of immune cells following implantation without using immunosuppressant. Using the syngeneic (Lewis male-Lewis female) and allogeneic (Brown Norway male-Lewis female) models and different tissue routes (subcutaneous vs omentum) for transplantation, normal blood vascular scaffolds were implanted which was converted to acellular vascular scaffolds by in vivo natural decellularization at the end of 2 months of observation. We also prepared chemically decellularized acellular scaffolds from normal untreated blood vascular scaffolds using a cocktail of chemicals which was also similarly placed in subcutaneous and omentum sites. Here, we applied in-depth quantitative proteomics along with histology and image analysis to comprehensively describe and compare the proteome of the natural and chemically decellularized scaffold. Our data confirm that site-specific advantages exist in modulating the ECM and regulating the immune responses (macrophage and T cells) following implantation, which possibly led to the production of an acellular scaffold (natural decellularization) under in vivo conditions. The current approach opens up the possibility to create tailor-made acellular scaffolds to build functional blood vessels. In addition, the identification of different tissue sites facilitates differential immune response against the scaffolds. This study provides a rich resource aimed toward an enhanced mechanistic understanding to study immune responses under similar settings in the field of transplantation and regenerative medicine.

Impact statementThe development of a scaffold helps in the preparation of a functional organ in the clinics. In the current study, we prepared an acellular vascular scaffold by utilizing site specific tissue changes and vis-a-vis compared with a conventionally chemically prepared biological scaffold at genomic and protein level, which helped us to identify immunological trigger following implantation. The current study which was carried out without any immunosuppressive agents could help to establish (a) alternative strategies for preparing biological scaffolds as well as (b) implantable sites as potential bioreactors to circumvent any adverse immune reactions for acceptance of the scaffold/implant post implantation. Reconstruction of nitrogenase predecessors suggests origin from maturase-like proteins The evolution of biological nitrogen fixation, uniquely catalyzed by nitrogenase enzymes, has been one of the most consequential biogeochemical innovations over lifes history. Though understanding the early evolution of nitrogen fixation has been a longstanding goal from molecular, biogeochemical, and planetary perspectives, its origins remain enigmatic. In this study, we reconstructed the evolutionary histories of nitrogenases, as well as homologous maturase proteins that participate in the assembly of the nitrogenase active-site cofactor but are not able to fix nitrogen. We combined phylogenetic and ancestral sequence inference with an analysis of predicted functionally divergent sites between nitrogenases and maturases to infer the nitrogen-fixing capabilities of their shared ancestors. Our results provide phylogenetic constraints to the emergence of nitrogen fixation and are consistent with a model wherein nitrogenases emerged from maturase-like predecessors. Though the precise functional role of such a predecessor protein remains speculative, our results highlight evolutionary contingency as a significant factor shaping the evolution of a biogeochemically essential enzyme.

SIGNIFICANCE STATEMENTThe origin of nitrogenase-catalyzed nitrogen fixation was a transformative event in lifes history, garnering long-term study from molecular, biogeochemical, and planetary perspectives. Reconstruction of ancestral nitrogenases suggests that the protein sequence space capable of yielding a nitrogen-fixing enzyme in the past was likely more constrained than previously thought. Specifically, here we show that nitrogenases likely evolved from ancestors that resemble maturases, homologs that today participate in nitrogenase cofactor assembly, contrary to the commonly accepted view that maturases evolved from a nitrogenase ancestor. We further submit that the molecular architecture that may have been required for nitrogenase origins was unlikely to have been shaped by the same environmental drivers often implicated in the evolution of nitrogen fixation. If this decoupling is found to be a recurring pattern in metabolic origins, then the presented results would undercut the common, systems-focused rationale of using ancient environmental conditions to explain the timing of critical and singular biogeochemical innovations in lifes past. Metacommunity dynamics and the detection of species associations in co-occurrence analyses: why patch disturbance matters O_LI1: Many statistical methods attempt to detect species associations - and so infer inter-specific interactions - from species co-occurrence patterns. Habitat heterogeneity and out-of-equilibrium colonization histories are well recognized as potentially causing species associations, even when interactions are absent. The potential for patch disturbance, a classical component of metacommunity dynamics, to also drive spurious species associations has however been over-looked.
C_LIO_LI2: Using a new general metacommunity model, we derive mathematical predictions regarding how patch disturbance would affect the patterns of species associations detected in "null" co-occurrence matrices. We also conduct numerical simulations to test our predictions and to compare the performance of several widespread statistical methods, including direct tests of pairwise independence, matrix permutation approaches and joint species distribution modelling.
C_LIO_LI3: We show how classical metacommunity dynamics can produce statistical associations, both positive and negative, even when species do not interact, when there is no habitat heterogeneity, and at equilibrium. This occurs as soon as there is some rate of patch disturbance (i.e. simultaneous extinction of several species in a patch) and/or a finite life-span of patches, a common feature of a broad range of plant, animal or microbial systems.
C_LIO_LI4: Patch disturbance can compromise species co-occurrence analyses and cause the artefactual detection of species associations if not taken into account. Including patch age (i.e. the time since the last patch disturbance event) as a covariate in a joint species distribution model can resolve the artefact. However, this requires additional data that often are not available in practice. We argue that the consequences of patch disturbance should not be underestimated when analysing species distribution patterns in metacommunity-like systems.
C_LI The potential of dynamic physiological traits in young tomato plants to predict field-yield performance To address the challenge of predicting tomato yields in the field, we used whole-plant functional phenotyping to evaluate water relations under well-irrigated and drought conditions. The genotypes tested are known to exhibit variability in their yields in wet and dry fields. The examined lines included two lines with recessive mutations that affect carotenoid biosynthesis, zeta z2083 and tangerine t3406, both isogenic to the processing tomato variety M82. The two mutant lines were reciprocally grafted onto M82, and multiple physiological characteristics were measured continuously, before, during and after drought treatment in the greenhouse. A comparative analysis of greenhouse and field yields showed that the whole-canopy stomatal conductance (gsc) in the morning and cumulative transpiration (CT) were strongly correlated with field measurements of total yield (TY: r2 = 0.9 and 0.77, respectively) and plant vegetative weight (PW: r2 = 0.6 and 0.94, respectively). Furthermore, the minimum CT during drought and the rate of recovery when irrigation was resumed were both found to predict resilience. Cluster-specific gene markers enhance Shigella and Enteroinvasive Escherichia coli in silico serotyping Shigella and enteroinvasive Escherichia coli (EIEC) cause human bacillary dysentery with similar invasion mechanisms and share similar physiological, biochemical and genetic characteristics. The ability to differentiate Shigella and EIEC from each other is important for clinical diagnostic and epidemiologic investigations. The existing genetic signatures may not discriminate between Shigella and EIEC. However, phylogenetically, Shigella and EIEC strains are composed of multiple clusters and are different forms of E. coli. In this study, we identified 10 Shigella clusters, 7 EIEC clusters and 53 sporadic types of EIEC by examining over 17,000 publicly available Shigella/EIEC genomes. We compared Shigella and EIEC accessory genomes to identify the cluster-specific gene markers or marker sets for the 17 clusters and 53 sporadic types. The gene markers showed 99.63% accuracy and more than 97.02% specificity. In addition, we developed a freely available in silico serotyping pipeline named Shigella EIEC Cluster Enhanced Serotype Finder (ShigEiFinder) by incorporating the cluster-specific gene markers and established Shigella/EIEC serotype specific O antigen genes and modification genes into typing. ShigEiFinder can process either paired end Illumina sequencing reads or assembled genomes and almost perfectly differentiated Shigella from EIEC with 99.70% and 99.81% cluster assignment accuracy for the assembled genomes and mapped reads respectively. ShigEiFinder was able to serotype over 59 Shigella serotypes and 22 EIEC serotypes and provided a high specificity with 99.40% for assembled genomes and 99.38% for mapped reads for serotyping. The cluster markers and our new serotyping tool, ShigEiFinder (https://github.com/LanLab/ShigEiFinder), will be useful for epidemiologic and diagnostic investigations.

Impact statementThe differentiation of Shigella strains from enteroinvasive E. coli (EIEC) is important for clinical diagnosis and public health epidemiologic investigations. The similarities between Shigella and EIEC strains make this differentiation very difficult as both share common ancestries within E. coli. However, Shigella and EIEC are phylogenetically separated into multiple clusters, making high resolution separation using cluster specific genomic markers possible. In this study, we identified 17 Shigella or EIEC clusters including five that were newly identified through examination of over 17,000 publicly available Shigella and EIEC genomes. We further identified an individual or a set of cluster-specific gene markers for each cluster using comparative genomic analysis. These markers can then be used to classify isolates into clusters and were used to develop an in silico pipeline, ShigEiFinder (https://github.com/LanLab/ShigEiFinder) for accurate differentiation, cluster typing and serotyping of Shigella and EIEC from Illumina sequencing reads or assembled genomes. This study will have broad application from understanding the evolution of Shigella/EIEC to diagnosis and epidemiology.

Data summarySequencing data have been deposited at the National Center for Biotechnology Information under BioProject number PRJNA692536.

RepositoriesRaw sequence data are available from NCBI under the BioProject number PRJNA692536. Library-free BoxCarDIA solves the missing value problem in label-free quantitative proteomics The last decade has seen significant advances in the application of quantitative mass spectrometry-based proteomics technologies to tackle important questions in plant biology. The current standard for quantitative proteomics in plants is the use of data-dependent acquisition (DDA) analysis with or without the use of chemical labels. However, the DDA approach preferentially measures higher abundant proteins, and often requires data imputation due to quantification inconsistency between samples. In this study we systematically benchmarked a recently developed library-free data-independent acquisition (directDIA) method against a state-of-the-art DDA label-free quantitative proteomics workflow for plants. We next developed a novel acquisition approach combining MS1-level BoxCar acquisition with MS2-level directDIA analysis that we call BoxCarDIA. DirectDIA achieves a 33% increase in protein quantification over traditional DDA, and BoxCarDIA a further 8%, without any changes in instrumentation, offline fractionation, or increases in mass-spectrometer acquisition time. BoxCarDIA, especially, offers wholly reproducible quantification of proteins between replicate injections, thereby addressing the long-standing missing-value problem in label-free quantitative proteomics. Further, we find that the gains in dynamic range sampling by directDIA and BoxCarDIA translate to deeper quantification of key, low abundant, functional protein classes (e.g., protein kinases and transcription factors) that are underrepresented in data acquired using DDA. We applied these methods to perform a quantitative proteomic comparison of dark and light grown Arabidopsis cell cultures, providing a critical resource for future plant interactome studies. Our results establish BoxCarDIA as the new method of choice in quantitative proteomics using Orbitrap-type mass-spectrometers, particularly for proteomes with large dynamic range such as that of plants. 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. 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. Limitations of composability of cis-regulatory elements in messenger RNA Genes are commonly abstracted into a coding sequence and cis-regulatory elements (CREs), such as promoters and terminators, and short sequence motifs within these regions. Modern cloning techniques allow easy assembly of synthetic genetic constructs from discrete cis-regulatory modules. However, it is unclear how much the contributions of CREs to gene expression depend on other CREs in the host gene. Using budding yeast, we probe the extent of composability, or independent effects, of distinct CREs. We confirm that the quantitative effect of a terminator on gene expression depends on both promoter and coding sequence. We then explore whether individual cisregulatory motifs within terminator regions display similar context dependence, using transcriptomewide datasets of mRNA decay. To test the extent of composability, we construct reporter genes consisting of combinations of motifs within various terminator contexts, paired with different promoters. Our results show that the effect of a motif on RNA abundance depends both on its host terminator, and also on the associated promoter sequence. This emphasises the need for improved motif inference that includes both local and global context effects, which in turn could aid in the accurate use of CREs for the engineering of synthetic genetic constructs. Elucidating Microscopic Events Driven by GTP Hydrolysis Reaction in Ras-GAP System with Semi-reactive Molecular Dynamics Simulation: Alternative Role of Phosphate Binding Loop as Mechanical Energy Storage ATPase and GTPase have been widely found as chemical energy-mechanical work transducers, whereas the physicochemical mechanisms are not satisfactorily understood. We addressed the problem by examining John Ross conjecture that repulsive Coulomb interaction between ADP/GDP and inorganic phosphate (Pi) does the mechanical work upon the system. We effectively simulated the consequence of GTP hydrolysis reaction in a complex system of Rat sarcoma (Ras) and GTPase activation protein (GAP) in the framework of classical molecular dynamics by switching force field parameters between the reactant and product systems. We then observed ca. 5 kcal/mol raise of potential energy about the phosphate-binding loop (P-loop) in Ras protein, indicating that the mechanical work generated via the GTP hydrolysis is converted into the local interaction energy and stored in the P-loop. Interestingly, this local energy storage in the P-loop depends on neither impulsive nor consecutive collisions of GDP and Pi with P-loop. Instead, GTP-GDP conversion itself does work on the Ras system, elevating the potential energy. These observations encourage us to challenge a conjecture previously given by Ross. We assert that triphosphate nucleotide hydrolyses do mechanical work by producing emergent steric interaction accompanied with relaxation, namely, a shift of biomolecular system to non-equilibrium state on the reshaped potential energy landscape. Recalling the universality of the P-loop motif among GTPases and ATPases, the observations that we obtained through this study would progress physicochemical understanding for the operating principles of GTP/ATP hydrolysis-driven biological nano-machines. Serological fingerprints link antiviral activity of therapeutic antibodies to affinity and concentration We assessed the affinities of the therapeutic monoclonal antibodies (mAbs) cilgavimab, tixagevimab, sotrovimab, casirivimab, and imdevimab to the receptor binding domain (RBD) of wild type, Delta, and Omicron spike. The Omicron RBD affinities of cilgavimab, tixagevimab, casirivimab, and imdevimab decreased by at least two orders of magnitude relative to their wild type equivalents, whereas sotrovimab binding was less severely impacted. These affinity reductions correlate with reduced antiviral activities of these antibodies, suggesting that simple affinity measurements can serve as an indicator for activity before challenging and time-consuming virus neutralization assays are performed. We also compared the properties of these antibodies to serological fingerprints (affinities and concentrations) of wild type RBD specific antibodies in 74 convalescent sera. The affinities of the therapeutic mAbs to wild type and Delta RBD were in the same range as the polyclonal response in the convalescent sera indicative of their high antiviral activities against these variants. However, for Omicron RBD, only sotrovimab retained affinities that were within the range of the polyclonal response, in agreement with its high activity against Omicron. Serological fingerprints thus provide important context to affinities and antiviral activity of mAb drugs and could guide the development of new therapeutics. 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. 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. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Development and Characterization of PAR-Trackers: New Tools for Detecting Poly(ADP-ribose) In Vitro and In Vivo AO_SCPLOWBSTRACTC_SCPLOWADP-ribosylation (ADPRylation) is a reversible post-translation modification resulting in the covalent attachment of ADP-ribose (ADPR) moieties on substrate proteins. Naturally-occurring protein motifs and domains, including WWEs, PBZs, and macrodomains, act as "readers" for protein-linked ADPR. Although recombinant, antibody-like ADPR detection reagents containing these readers have facilitated the detection of ADPR, they are limited in their ability to capture the dynamic nature of ADPRylation. Herein, we describe and characterize a set of poly(ADP-ribose) (PAR) Trackers (PAR-Ts) - optimized dimerization-dependent or split-protein reassembly PAR sensors in which a naturally occurring PAR binding domain, WWE, was fused to both halves of dimerization-dependent GFP (ddGFP) or split Nano Luciferase (NanoLuc), respectively. We demonstrate that these new tools allow the detection and quantification of PAR levels in extracts, living cells, and living tissues with greater sensitivity, as well as temporal and spatial precision. Importantly, these sensors detect changes in cellular ADPR levels in response to physiological cues (e.g., hormone-dependent induction of adipogenesis without DNA damage), as well as xenograft tumor tissues in living mice. Our results indicate that PAR Trackers have broad utility for detecting ADPR in many different experimental and biological systems. Centriole distal-end proteins CP110 and Cep97 influence centriole cartwheel growth at the proximal-end Centrioles are composed of a central cartwheel tethered to nine-fold symmetric microtubule (MT) blades. The centriole cartwheel and MTs are thought to grow from opposite ends of these organelles, so it is unclear how they coordinate their assembly. We previously showed that an oscillation of Polo-like kinase 4 (Plk4) helps to initiate and time the growth of the cartwheel at the proximal end. Here, we show that CP110 and Cep97 form a complex close to the distal-end of the centriole MTs whose levels rise and fall as the new centriole MTs grow, entrained by the core Cdk/Cyclin oscillator that drives the nuclear divisions in these embryos. These CP110/Cep97 dynamics, however, do not appear to time the period of centriole MT growth directly. Instead, we find that changing the levels of CP110/Cep97 alters the Plk4 oscillation and the growth of the cartwheel at the proximal end. These findings reveal an unexpected crosstalk between factors normally concentrated at opposite ends of the growing centrioles, which may help to coordinate centriole growth.