A class of identifiable phylogenetic birth-death models In a striking result, Louca and Pennell (2020) recently proved that a large class of phylogenetic birth-death models are statistically unidentifiable from lineage-through-time (LTT) data: any pair of sufficiently smooth birth and death rate functions is "congruent" to an infinite collection of other rate functions, all of which have the same likelihood for any LTT vector of any dimension. As Louca and Pennell argue, this fact has distressing implications for the thousands of studies that have utilized birth-death models to study evolution.

In this paper, we qualify their finding by proving that an alternative and widely used class of birth-death models is indeed identifiable. Specifically, we show that piecewise constant birth-death models can, in principle, be consistently estimated and distinguished from one another, given a sufficiently large extant time tree and some knowledge of the present-day population. Subject to mild regularity conditions, we further show that any unidentifiable birth-death model class can be arbitrarily closely approximated by a class of identifiable models. The sampling requirements needed for our results to hold are explicit, and are expected to be satisfied in many contexts such as the phylodynamic analysis of a global pandemic.
Natural selection promotes the evolution of recombination 1: between the products of natural selection Shuffling ones genetic material with another individual seems a risky endeavor more likely to decrease than to increase offspring fitness. This intuitive argument is commonly employed to explain why the ubiquity of sex and recombination in nature is enigmatic. It is predicated on the notion that natural selection assembles selectively well-matched combinations of genes that recombination would break up resulting in low-fitness offspring - a notion often stated in the literature as a self-evident premise. We show however that, upon closer examination, this premise is flawed: we find to the contrary that natural selection in fact has an encompassing tendency to assemble selectively mismatched gene combinations; recombination breaks up these selectively mismatched combinations (on average), assembles selectively matched combinations, and should thus be favored. The new perspective our findings offer suggests that sex and recombination are not so enigmatic but are instead unavoidable byproducts of natural selection.
Positive Selection and Heat-Response Transcriptomes Reveal Adaptive Features of the Brassicaceae Desert Model, Anastatica hierochuntica O_LIPlant adaptation to a desert environment and its endemic heat stress is poorly understood at the molecular level. The naturally heat-tolerant Brassicaceae species Anastatica hierochuntica is an ideal extremophyte model to identify genetic adaptations that have evolved to allow plants to tolerate heat stress and thrive in deserts.
C_LIO_LIWe generated an A. hierochuntica reference transcriptome and pinpointed extremophyte adaptations by comparing Arabidopsis thaliana and A. hierochuntica transcriptome responses to heat and identifying positively selected genes in A. hierochuntica.
C_LIO_LIThe two species exhibit similar transcriptome adjustment in response to heat and the A. hierochuntica transcriptome does not exist in a constitutive heat "stress-ready" state. Furthermore, the A. hierochuntica global transcriptome as well as heat-responsive orthologs, display a lower basal and higher heat-induced expression than in A. thaliana. Genes positively selected in multiple extremophytes are associated with stomatal opening, nutrient acquisition, and UV-B induced DNA repair while those unique to A. hierochuntica are consistent with its photoperiod-insensitive, early-flowering phenotype.
C_LIO_LIWe suggest that evolution of a flexible transcriptome confers the ability to quickly react to extreme diurnal temperature fluctuations characteristic of a desert environment while positive selection of genes involved in stress tolerance and early flowering could facilitate an opportunistic desert lifestyle.
C_LI Catchet-MS identifies IKZF1-targeting Thalidomide analogues as novel HIV-1 latency reversal agents A major pharmacological strategy toward HIV cure aims to reverse latency in infected cells as a first step leading to their elimination. While the unbiased identification of molecular targets physically associated with the latent HIV-1 provirus would be highly valuable to unravel the molecular determinants of HIV-1 transcriptional repression and latency reversal, due to technical limitations, this has not been possible. Here we use a dCas9 targeted chromatin and histone enrichment strategy coupled to mass spectrometry (Catchet-MS) to describe the protein composition of the latent and activated HIV-1 5LTR. Catchet-MS identified known and novel latent 5LTR-associated host factors. Among these, IKZF1 is a novel HIV-1 transcriptional repressor, required for Polycomb Repressive Complex 2 recruitment to the LTR. We find the clinically advanced thalidomide analogue iberdomide, and the FDA approved analogues lenalidomide and pomalidomide, to be novel LRAs that, by targeting IKZF1 for degradation, reverse HIV-1 latency in CD4+T-cells isolated from virally suppressed people living with HIV-1.

One Sentence SummarydCas9 targeted chromatin and histone enrichment for mass spectrometry (Catchet-MS) led to the identification of IKZF1-targeting thalidomide analogues as novel HIV-1 latency reversal agents A behavioral paradigm for cortical control of a robotic arm by freely moving rats in a one-dimensional two-target reaching task Controlling the trajectory of a neuroprosthetic device to reach multiple targets is a commonly used brain-machine interface (BMI) task in primates and has not been available for rodents yet. Here, we describe a novel experimental paradigm which enables this task for rats in one-dimensional space for reaching two distant targets depending on their limited cognitive and visual capabilities compared to primates. An online transform was used to convert the activity of a pair of primary motor cortex (M1) units into two robotic actions. The rats were shaped to adapt to the transform and direct the robotic actuator toward the selected target by modulating the activity of the M1 neurons. All three rats involved in the study were capable of achieving randomly selected targets with at least 78% accuracy. A total of 9 out of 16 pairs of units examined were eligible for exceeding this success criterion. Two out of three rats were capable of reversal learning, where the mapping between the activity of the unit pairs and the robotic actions were reversed. The present work is the first demonstration of trajectory-based control of a neuroprosthetic device by rodents to reach two distant targets using visual feedback. The paradigm introduced here may be used as a cost-effective platform for elucidating the information processing principles in the neural circuits related to neuroprosthetic control and for studying the performance of novel BMI technologies using freely moving rats. Game of theories: an agent-based model of how variation in contest cost-benefits determines the evolution of contest resolution rules Game-theory based models are used to understand rules that animals use to settle contests over indivisible resources. However, the empirical literature of contests indicates controversial support to models, with some species supporting different models and other species showing no support to any model. Since strategies used to resolve contests may have different associated costs, it is possible that different conditions have determined the evolution of distinct assessment strategies used by animals. We used an agent-based model to explore the importance of the following conditions: resource availability, probability of reproduction with resource, and damage costs on evolution of assessment strategies. We used self- and mutual-assessment models as a heurist framework to build agents with different assessment strategies. In our model, agents competed for resources in scenarios with different combinations of resource availability, probability of reproduction with resource, and damage costs. We found that agents following the self-assessment with damage strategy were prevalent in scenarios with no probability of reproduction without the resource, independently of other variables. We also found that agents following the non-aggressive strategy occurred in all scenarios. However, agents using the non-aggressive strategy were prevalent only in scenarios with probability of reproduction with the resource. Finally, we observed that agents using mutual-assessment occurred only in a scenario with high risk of damage, low availability of resources, and with probability of reproduction without the resource. These results indicate that agents following the self-assessment with damage and non-aggressive strategies may be able to stay at most scenarios. Mitochondria preserve an autarkic one-carbon cycle to confer growth-independent cancer cell migration and metastasis Progression of primary cancer to metastatic disease is the most common cause of death in cancer patients with minimal treatment options available. Canonical drugs target mainly the proliferative capacity of cancer cells, which often leaves slow-proliferating, persistent cancer cells unaffected. Metabolic determinants that contribute to growth-independent functions supporting resistance and metastatic dissemination are still poorly understood.

In the present study, we revealed that antifolate treatment results in an uncoupled and autarkic mitochondrial one-carbon (1C) metabolism allowing sustained serine catabolism and formate overflow when cytosolic 1C metabolism is impaired. Interestingly, antifolate dependent growth-arrest did not correlate with decreased migration capacity. Therefore, using the antifolate Methotrexate as a tool compound allowed us to disentangle proliferation and migration to profile the metabolic phenotype of migrating (growth-arrested) cells. Supported by an increased NAD/NADH ratio, we observed increased serine de novo synthesis and increased serine catabolism to formate.

Consequently, inhibition of serine de novo synthesis using the competitive PHGDH-inhibitor BI-4916 or direct inhibition of mitochondrial 1C metabolism reduced cancer cell migration. Using an orthotopic breast cancer model, we show that sole inhibition of mitochondrial serine catabolism does not affect primary tumor growth but strongly inhibits pulmonary metastasis.

We conclude that mitochondrial 1C metabolism, despite being dispensable for proliferative capacities, confers an advantage to cancer cells by supporting their motility potential.

Our results improve our understanding of 1C metabolism and of metabolic determinants that support the process of cancer cell migration and metastasis. Electrical Stimulated GLUT4 Signaling Attenuates Critical Illness-Associated Muscle Wasting BackgroundCritical illness myopathy (CIM) is a debilitating condition characterized by the preferential loss of the motor protein myosin. CIM is a byproduct of critical care, attributed to impaired recovery, long-term complications, and mortality. CIM pathophysiology is complex, heterogeneous and remains incompletely understood, however loss of mechanical stimuli contributes to critical illness associated muscle atrophy and weakness. Passive mechanical loading (ML) and electrical stimulation (ES) therapies augment muscle mass and function. While having beneficial outcomes, the mechanistic underpinning of these therapies is less known. Therefore, here we aimed to assess the mechanism by which chronic supramaximal ES ameliorates CIM in a unique experimental rat model of critical care.

MethodsRats were subjected to 8 days critical care conditions entailing deep sedation, controlled mechanical ventilation, and immobilization with and without direct soleus ES. Muscle size and function were assessed at the single cell level. RNAseq and Western blotting were employed to understand the mechanisms driving ES muscle outcomes in CIM.

ResultsFollowing 8 days of controlled mechanical ventilation and immobilization, soleus muscle mass, Myosin:Actin ratio and single muscle fiber maximum force normalized to cross-sectional area (specific force) were reduced by 40-50% (p< 0.0001). ES significantly reduced the loss of soleus muscle fiber cross-sectional area (CSA) and Myosin:Actin ratio by approximately 30% (p< 0.05) yet failed to effect specific force. RNAseq pathway analysis revealed downregulation of insulin signaling in the soleus muscle following critical care and GLUT4 trafficking was reduced by 55% leading to an 85% reduction of muscle glycogen content (p< 0.01). ES promoted phosphofructokinase and insulin signaling pathways to control levels (p< 0.05), consistent with the maintenance of GLUT4 translocation and glycogen levels. AMPK, but not AKT, signaling pathway was stimulated following ES, where the downstream target TBC1D4 increased 3 logFC (p= 0.029) and AMPK-specific P-TBC1D4 levels were increased approximately 2-fold (p= 0.06). Reduction of muscle protein degradation rather than protein synthesis promoted soleus CSA, as ES reduced E3 ubiquitin proteins, Atrogin-1 (p= 0.006) and MuRF1 (p= 0.08) by approximately 50%, downstream of AMPK-FoxO3.

ConclusionsES maintained GLUT4 translocation through increased AMPK-TBC1D4 signaling leading to improved muscle glucose homeostasis. Soleus CSA and myosin content was promoted through reduced protein degradation via AMPK-FoxO3 E3 ligases, Atrogin-1 and MuRF1. These results demonstrate chronic supramaximal ES reduces critical care associated muscle wasting, preserved glucose signaling and reduced muscle protein degradation in CIM. First detection of herpesvirus and frequency of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. 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. Co-infection of Candidatus Piscichlamydia trichopodus (order Chlamydiales) and Henneguya sp (Myxosporea, Myxobolidae) in snakeskin gourami Trichopodus pectoralis (Regan 1910) The present study describes a simultaneous infection of a novel Chlamydia-like organism (CLO) with a Myxozoa parasite, Henneguya sp. in snakeskin gourami Trichopodus pectoralis in Thailand. A new CLO is proposed "Candidatus Piscichlamydia trichopodus" (CPT) based on 16S rRNA phylogenetic analysis. Systemic intracellular CPT infection was confirmed by histological examination, in situ hybridization, PCR assay, and sequencing of 16S rRNA. This novel pathogen belongs to the order Chlamydiales but differs in certain aspects from other species. The histopathological changes associated with CPT infection were different from the typical pathological lesions of epitheliocystis caused by previously known CLO. Unlike other CLO, CPT localized in the connective tissue rather than in the epithelial cells and formed smaller clumps of intracellular bacteria that stained dark blue with hematoxylin. On the other hand, typical myxospores of the genus Henneguya with tails were observed in the gill sections. Infection with Henneguya sp. resulted in extensive destruction of the gill filaments, most likely leading to respiratory distress. Due to the frequency of co-infections and the unavailability of culture methods for CLO and Henneguya sp, it was difficult to determine which pathogens were directly responsible for the associated mortality. However, co-infections may increase the negative impact on the host and the severity of the disease. Given the commercial importance of the snakeskin gourami and its significant aquaculture potential, the findings of this study are important for further studies on disease prevention. Influence of delayed density and ultraviolet radiation on caterpillar baculovirus infection and mortality O_LIInfectious disease is an important potential driver of population cycles, but this must occur through delayed density-dependent infection and resulting fitness effects. Delayed density-dependent infection by baculoviruses can be caused by environmental persistence of viral occlusion bodies, which can be influenced by environmental factors. In particular, ultraviolet radiation is potentially important in reducing the environmental persistence of viruses by inactivating viral occlusion bodies.
C_LIO_LIDelayed density-dependent viral infection has rarely been observed empirically at the population level although theory predicts that it is necessary for these pathogens to drive population cycles. Similarly, field studies have not examined the potential effects of ultraviolet radiation on viral infection rates in natural animal populations. We tested if viral infection is delayed density-dependent with the potential to drive cyclic dynamics and if ultraviolet radiation influences viral infection levels.
C_LIO_LIWe censused 18 Ranchmans tiger moth (Arctia virginalis) populations across nearly 9{degrees} of latitude over two years and quantified the effects of direct and delayed density and ultraviolet radiation on baculovirus infection rates, infection severity, and survival to adulthood. Caterpillars were collected from each population in the field and reared in the laboratory. Baculovirus has not previously been described infecting Arctia virginalis, and we used genetic methods to confirm the identity of the virus.
C_LIO_LIWe found that infection rate, infection severity, and survival to adulthood exhibited delayed density-dependence. Ultraviolet radiation in the previous summer decreased infection severity, and increased survival probability of the virus. Structural equation modelling indicated that the effect of lagged density on moth survival was mediated through infection rate and infection severity, and was 2.5 fold stronger than the effect of ultraviolet radiation on survival through infection severity. We successfully amplified polh, lef-8, and lef-9 viral genes from caterpillar samples, and BLAST search results confirmed that the virus was a nucleopolyhedrovirus.
C_LIO_LIOur findings provide clear evidence that delayed density dependence can arise through viral infection rate and severity in insects, which supports the role of viral disease as a potential mechanism, among others, that may drive insect population cycles. Furthermore, our findings support predictions that ultraviolet radiation can modify viral disease dynamics in insect populations, most likely through attenuating viral persistence in the environment.
C_LI Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Toxoplasma bradyzoites exhibit physiological plasticity of calcium and energy stores controlling motility and egress. Toxoplasma gondii has evolved different developmental stages for disseminating during acute infection (i.e. tachyzoites) and for establishing chronic infection (i.e. bradyzoites). Calcium ion (Ca2+) signaling tightly regulates the lytic cycle of tachyzoites by controlling microneme secretion and motility to drive egress and cell invasion. However, the roles of Ca2+ signaling pathways in bradyzoites remain largely unexplored. Here we show that Ca2+ responses are highly restricted in bradyzoites and that they fail to egress in response to agonists. Development of dual-reporter parasites revealed dampened calcium responses and minimal microneme secretion by bradyzoites induced in vitro or harvested from infected mice and tested ex vivo. Ratiometric Ca2+ imaging demonstrated lower Ca2+ basal levels, reduced magnitude, and slower Ca2+ kinetics in bradyzoites compared with tachyzoites stimulated with agonists. Diminished responses in bradyzoites were associated with down-regulation of calcium ATPases involved in intracellular Ca2+ storage in the endoplasmic reticulum (ER) and acidocalcisomes. Once liberated from cysts by trypsin digestion, bradyzoites incubated in glucose plus calcium rapidly restored their intracellular Ca2+ and ATP stores leading to enhanced gliding. Collectively, our findings indicate that intracellular bradyzoites exhibit dampened Ca2+ signaling and lower energy levels that restrict egress, and yet upon release they rapidly respond to changes in the environment to regain motility. Comparative toxicogenomics of glyphosate and Roundup herbicides by mammalian stem cell-based genotoxicity assays and molecular profiling in Sprague-Dawley rats Whether glyphosate-based herbicides (GBHs) are more potent than glyphosate alone at activating cellular mechanisms, which drive carcinogenesis remains controversial. As GBHs are more cytotoxic that glyphosate, we reasoned they may also be more capable of activating carcinogenic pathways. We tested this hypothesis by comparing the effects of glyphosate with Roundup GBHs both in vitro and in vivo. First, glyphosate was compared with representative GBHs namely MON 52276 (EU), MON 76473 (UK) and MON 76207 (USA) using the mammalian stem cell-based ToxTracker system. Here, MON 52276 and MON 76473, but not glyphosate and MON 76207, activated oxidative stress and unfolded protein responses. Second, molecular profiling of liver was performed in female Sprague-Dawley rats exposed to glyphosate or MON 52276 (both at 0.5, 50, 175 mg/kg bw/day glyphosate) for 90 days. MON 52276 but not glyphosate increased hepatic steatosis and necrosis. MON 52276 and glyphosate altered the expression of genes in liver reflecting TP53 activation by DNA damage and circadian rhythm regulation. Genes most affected in liver were similarly altered in kidneys. Small RNA profiling in liver showed decreased amounts of miR-22 and miR-17 from MON 52276 ingestion. Glyphosate decreased mir-30 while miR-10 levels were increased. DNA methylation profiling of liver revealed 5,727 and 4,496 differentially methylated CpG sites between the control and glyphosate and MON 52276 exposed animals respectively. Apurinic/apyrimidinic DNA damage formation in liver was increased with glyphosate exposure. Altogether, our results show that Roundup formulations cause more biological changes linked with carcinogenesis than glyphosate. Computation of the mitochondrial age distribution along the axon length We describe a compartmental model of mitochondrial transport in axons, which we apply to compute mitochondrial age at different distances from the soma. The model predicts that the average mitochondrial age at the tip of an axon with a length of 1 cm is approximately 22 hours. The mitochondria are youngest closest to the soma and their age scales approximately linearly with distance from the soma. To the best of the authors knowledge, this is the first attempt to predict the spatial distribution of mitochondrial age within an axon. A sensitivity study of the mean age of mitochondria to various model parameters is also presented. Changes to the mtDNA copy number during yeast culture growth We show that the mtDNA copy number in growing cultures of the yeast Saccharomyces cerevisiae increases by a factor of up to 4, being lowest ([~]10 per haploid genome) and stable during rapid fermentative growth, and highest at the end of the respiratory phase. When yeast are grown on glucose, the onset of the mtDNA copy number increase coincides with the early stages of the diauxic shift, and the increase continues through respiration. A lesser yet still substantial copy number increase occurs when yeast are grown on a nonfermentable carbon source, i.e. when there is no diauxic shift. The mtDNA copy number increase during and for some time after the diauxic shift is not driven by an increase in cell size. The copy number increase occurs in both haploid and diploid strains, but is markedly attenuated in a diploid wild isolate that is a ready sporulator. Strain-to-strain differences in mtDNA copy number are least apparent in fermentation and most apparent in late respiration or stationary phase. While changes in mitochondrial morphology and function were previously known to accompany changes in physiological state, it had not been previously shown that the mtDNA copy number changes substantially over time in a clonal growing culture. The mtDNA copy number in yeast is therefore a highly dynamic phenotype. Structural variant-based pangenome construction has low sensitivity to variability of haplotype-resolved bovine assemblies Advantages of pangenomes over linear reference assemblies for genome research have recently been established. However, potential effects of sequence platform and assembly approach, or of combining assemblies created by different approaches, on pangenome construction have not been investigated. We generated haplotype-resolved assemblies from the offspring of three bovine trios representing increasing levels of heterozygosity that each demonstrate a substantial improvement in contiguity, completeness, and accuracy over the current Bos taurus reference genome. Diploid coverage as low as 20x for HiFi or 60x for ONT was sufficient to produce two haplotype-resolved assemblies meeting the standards set by the Vertebrate Genome Project. Structural variant-based pangenomes created from the haplotype-resolved assemblies demonstrated significant consensus regardless of sequence platform, assembler algorithm, or coverage. Inspecting pangenome topologies identified 90 thousand structural variants including 931 overlapping with coding sequences; this approach revealed variants affecting QRICH2, PRDM9, HSPA1A, TAS2R46, and GC that have potential to affect phenotype. Unbiased integration of single cell transcriptome replicates Single cell transcriptomic approaches are becoming mainstream, with replicate experiments commonly performed with the same single cell technology. Methods that enable integration of these datasets by removing batch effects while preserving biological information are required for unbiased data interpretation. Here we introduce Canek for this purpose. Canek leverages information from mutual nearest neighbor to combine local linear corrections with cell-specific non-linear corrections within a fuzzy logic framework. Using a combination of real and synthetic datasets, we show that Canek corrects batch effects while introducing the least amount of bias compared with competing methods. Canek is computationally efficient and can easily integrate thousands of single-cell transcriptomes from replicated experiments. An integrative analysis to predict the Active Ingredients and explore Polypharmacological Mechanisms of Orthosiphon stamineus Benth. for human health BackgroundOrthosiphon stamineus Benth. was a dietary supplement or a traditional folk herb with widespread clinical application, but it was lack of comprehensive understanding of its active ingredients and polypharmacological mechanisms. The aim of this work was to systematically study its natural compounds and the molecular mechanisms of the O. stamineus via network pharmacology.

MethodsCompounds from O. stamineus were collected by literature retrieval and evaluated by SwissADME with the physicochemical property (ADMET model) and the likelihood for a natural medicine. The connection of active ingredients and target genes was built and confirmed by Cytoscape and AutoDock vina. Then, Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were separately performed to obtain a more in-depth understanding of O. stamineus. Finally, the relationship among active ingredients, targets, and diseases was built to clarify the polypharmacological mechanisms and found relevant active substances for further drug discovery.

ResultsA total of 159 compounds of O. stamineus were collected, and 22 potentially active ingredients were screened out. The Ingredient-Target Interaction Network was built with 12 flavonoids, 3 diterpenes, 3 phenols, and 4 volatile oils, and 65 targets. The Docking analysis indicated that the ingredient-target interaction network was reliable; most ligand-receptor had a strong binding affinity (lowest binding energy: -6.9 kcal/mol). After pathway analysis, 185 significant biological processes and 36 signal pathways were found, and the ingredient-target-disease network of O. stamineus was constructed for polypharmacological mechanisms.

ConclusionOur study clarified the polypharmacological mechanisms via the relationship among active ingredients, targets, and diseases and provided better guidance for subsequent experiments and potential active ingredients for drug discovery or health promotion. 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. The impact of housing conditions on mesenchymal stromal/stem cell populations differ between adipose tissue and skeletal muscle BackgroundIn pigs, the ratio between lean mass and fat mass in the carcass determines production efficiency and is strongly influenced by the number and size of cells in tissues. During growth, the increase in the number of cells results from the recruitment of different populations of multipotent mesenchymal stromal/stem cells (MSCs) residing in the tissues. We hypothesized that the impact of hygiene of housing conditions during growth on the proportions of MSCs in adipose tissue and skeletal muscle may differ between pigs with different residual feed intake (RFI), a measure of feed efficiency.

MethodsAt the age of 12 weeks, Large White pigs from two lines divergently selected for low and high RFI were housed in two contrasting hygiene conditions (good vs poor). After six weeks, pigs were slaughtered (n = 30; 5-9/group). Samples of subcutaneous adipose tissue and longissimus skeletal muscle were collected, and cells from the stromal vascular fraction (SVF), which includes mesenchymal stromal/stem cell populations, were isolated from each tissue. Adipose and muscle cell populations from the SVF were phenotyped by flow cytometry using antibodies that targeted different cell surface markers (CD45 to separate hematopoietic cells from MSCs; CD34, CD38, CD56 and CD140a to identify MSC populations with adipogenic and/or myogenic potential).

ResultsAdipose tissue and muscle shared some common MSC populations although MSC diversity was higher in muscle than in adipose tissue. In muscle, the CD45-CD56+CD34-CD140a+ and CD45-CD56+CD34+CD140a+ cell populations were abundant. Of these two cell populations, only the proportions of CD45-CD56+CD34+CD140a+ cells increased (P < 0.05) in pigs housed in poor hygiene as compared with pigs in good hygiene conditions. For the CD45-CD56-CD34- cell population, present in low proportion, there was an interaction between hygiene condition and genetic line (P < 0.05) with a decrease in low RFI pigs housed in poor hygiene conditions. In adipose tissue, the two abundant MSC populations were CD45-CD56-CD34- and CD45-CD56+CD34-. The proportion of CD45-CD56-CD34- cells increased (P < 0.05) whereas the proportion of CD45-CD56+CD34- tended to decrease (P < 0.1) in pigs housed in poor conditions. This study shows that the proportions of some MSC populations were affected by hygiene of housing conditions in a tissue-dependent manner in pigs of both RFI lines. It suggests that MSCs may play a significant role in adipose tissue and skeletal muscle homeostasis and may influence later growth and body composition in growing animals. A humanized nanobody phage display library yields potent binders of SARS CoV-2 spike Neutralizing antibodies targeting the SARS-CoV-2 spike protein have shown a great preventative/therapeutic potential. Here, we report a rapid and efficient strategy for the development and design of SARS-CoV-2 neutralizing humanized nanobody constructs with sub-nanomolar affinities and nanomolar potencies. CryoEM-based structural analysis of the nanobodies in complex with spike revealed two distinct binding modes. The most potent nanobody, RBD-1-2G(NCATS-BL8125), tolerates the N501Y RBD mutation and remains capable of neutralizing the B.1.1.7 (Alpha) variant. Molecular dynamics simulations provide a structural basis for understanding the neutralization process of nanobodies exclusively focused on the spike-ACE2 interface with and without the N501Y mutation on RBD. A primary human airway air-lung interface (ALI) ex vivo model showed that RBD-1-2G-Fc antibody treatment was effective at reducing viral burden following WA1 and B.1.1.7 SARS-CoV-2 infections. Therefore, this presented strategy will serve as a tool to mitigate the threat of emerging SARS-CoV-2 variants. The cyclic dinucleotide 2'3'-cGAMP induces a broad anti-bacterial and anti-viral response in the sea anemone Nematostella vectensis In mammals, cyclic dinucleotides (CDNs) bind and activate STING to initiate an anti-viral type I interferon response. CDNs and STING originated in bacteria and are present in most animals. By contrast, interferons are believed to have emerged in vertebrates; thus, the function of CDN signaling in invertebrates is unclear. Here, we use a CDN, 23-cGAMP, to activate immune responses in a model cnidarian invertebrate, the starlet sea anemone Nematostella vectensis. Using RNA-Seq, we found that 23-cGAMP induces robust transcription of both anti-viral and anti-bacterial genes in N. vectensis. Many of the anti-viral genes induced by 23-cGAMP are homologs of vertebrate interferon stimulated genes, implying that the interferon response predates the evolution of interferons. Knockdown experiments identified a role for NF-{kappa}B in specifically inducing anti-bacterial genes downstream of 23-cGAMP. Some of these putative anti-bacterial genes were also found to be induced during Pseudomonas aeruginosa infection. We characterized the protein product of one of the putative anti-bacterial genes, the N. vectensis homolog of Dae4, and found that it has conserved anti-bacterial activity. This work suggests that a broad anti-bacterial and anti-viral transcriptional response is an evolutionarily ancestral output of 23-cGAMP signaling in animals.

Significance statementCyclic dinucleotides are signaling molecules that originated in bacteria and were subsequently acquired and co-opted by animals for immune signaling. The major cyclic dinucleotide signaling pathway in mammals results in the production of anti-viral molecules called interferons. Invertebrates such as sea anemones lack interferons, and thus it was unclear whether cyclic dinucleotide signaling would play a role in immunity in these animals. Here we report that in the anemone Nematostella vectensis, cyclic dinucleotides activate both anti-viral and anti-bacterial immune responses, and do so through a conserved pathway. These results provide insights into the evolutionary origins of innate immunity, and suggest a broader ancestral role for cyclic dinucleotide signaling that evolved toward more specialized anti-viral functions in mammals. FGF8-mediated signaling regulates tooth developmental pace during odontogenesis The developing human and mouse teeth constitute an ideal model system to study the regulatory mechanism underlying organ growth control due to the fact that their teeth share highly conserved and well-characterized developmental processes and their developmental tempo vary notably. In the current study, we manipulated heterogenous recombination between human and mouse dental tissues and demonstrate that the dental mesenchyme dominates the tooth developmental tempo and FGF8 could be a critical player during this developmental process. Forced activation of FGF8 signaling in the dental mesenchyme of mice promoted cell proliferation, prevented cell apoptosis via p38 and perhaps PI3K-Akt intracellular signaling, and impelled the transition of the cell cycle from G1- to S-phase in the tooth germ, resulting in the slowdown of the tooth developmental pace. Our results provide compelling evidence that extrinsic signals can profoundly affect tooth developmental tempo and the dental mesenchymal FGF8 could be a pivotal factor in controlling developmental pace in a non-cell-autonomous manner during mammalian odontogenesis. Comprehensive O-GlcNAc glycoproteomics on NOTCH1 EGF repeats implicated unique Lewis X epitopes in mammals The O-GlcNAc modification of Notch receptors regulates Notch ligand interactions in a manner distinct from other forms of O-glycans on epidermal growth factor-like (EGF) repeats of Notch receptors. Although many proteins, besides Notch receptors, are expected to be O-GlcNAcylated by EGF domain-specific O-GlcNAc transferase (EOGT), only a small number of proteins have been reported to be modified in vivo, and elongated O-GlcNAc glycans have not been extensively explored. To extend our view of the specificity and variety of the glycan modification, we conducted a comprehensive analysis of O-GlcNAc glycans on NOTCH1 in mammals. Mass spectrometric analysis of NOTCH1 fragments expressed in HEK293T cells revealed that several EGF domains with putative O-GlcNAcylation sites were hardly modified with O-GlcNAc. Although amino acid residues before the modification site are preferentially occupied with aromatic residues, Phe and Tyr are preferable to Trp for the apparent modification with O-GlcNAc. Furthermore, a minor form of fucosylated O-GlcNAc glycans was detected in a subset of EGF domains. Fucosylation of O-GlcNAc glycans was enhanced by FUT1, FUT2, or FUT9 expression. The FUT9-dependent Lewis X epitope was confirmed by immunoblotting using an anti-Lewis X antibody. As expected from the similarity in the glycan structures, the Lexis X antigen was detected on O-fucose glycans. Our results refined the putative consensus sequence for the EOGT-dependent extracellular O-GlcNAc modification in mammals and revealed the structural diversity of functional Notch O-glycans.