Single cell genomics reveals plastid-lacking Picozoa are close relatives of red algae The endosymbiotic origin of plastids from cyanobacteria gave eukaryotes photosynthetic capabilities and launched the diversification of countless forms of algae. These primary plastids are found in members of the eukaryotic supergroup Archaeplastida. All known archaeplastids still retain some form of primary plastids, which are widely assumed to have a single origin. Here, we used single-cell genomics from natural samples combined with phylogenomics to infer the evolutionary origin of the phylum Picozoa, a globally distributed but seemingly rare group of marine microbial heterotrophic eukaryotes. Strikingly, the analysis of 43 single-cell genomes shows that Picozoa belong to Archaeplastida, specifically related to red algae and the phagotrophic rhodelphids. These picozoan genomes support the hypothesis that Picozoa lack a plastid, and further reveal no evidence of an early cryptic endosymbiosis with cyanobacteria. These findings change our understanding of plastid evolution as they either represent the first complete plastid loss in a free-living taxon, or indicate that red algae and rhodelphids obtained their plastids independently of other archaeplastids.
The important choice of reference environment in microevolutionary climate response predictions It is well documented that individuals of wild populations can adjust to climate change by means of phenotypic plasticity, but few reports on adaptation by means of genetically based microevolution caused by selection. Disentanglement of these separate effects requires that the reference environment (the environmental zero-point) is defined, and this should not be done arbitrarily. The problem is that an error in the reference environment may lead to large errors in predicted microevolution. Together with parameter values and initial mean trait values, the reference environment can be estimated from environmental, phenotypic and fitness data. A prediction error method for this purpose is described, with the feasibility shown by simulations. As shown in a toy example, an estimated reference environment may have large errors, especially for small populations. This may still be a better choice than use of an initial environmental value in a recorded time series, or the mean value, which is often used. Another alternative may be to use the mean value of a past and stationary stochastic environment, which the population is judged to have been fully adapted to, in the sense that the expected geometric mean fitness was at a global maximum. Exceptions are cases with constant phenotypic plasticity, where the microevolutionary changes per generation follow directly from phenotypic and environmental data, independent of the chosen reference environment.
Temporal Scaling of Human Scalp-Recorded Potentials Much of human behaviour is governed by common processes that unfold over varying timescales. Standard event-related potential analysis assumes fixed-latency responses relative to experimental events. However, recent single unit recordings in animals have revealed neural activity scales to span different durations during behaviours demanding flexible timing. Here, we employed a general linear modelling approach using a novel combination of fixed-duration and variable-duration regressors to unmix fixed-time and scaled-time components in human magneto/electroencephalography (M/EEG) data. We use this to reveal consistent temporal scaling of human scalp-recorded potentials across four independent EEG datasets, including interval perception, production, prediction and value-based decision making. Between-trial variation in the temporally scaled response predicts between-trial variation in subject reaction times, demonstrating the relevance of this temporally scaled signal for temporal variation in behaviour. Our results provide a general approach for studying flexibly timed behaviour in the human brain. An AI-guided signature reveals the nature of the shared proximal pathways of host immune response in MIS-C and Kawasaki disease A significant surge in cases of multisystem inflammatory syndrome in children (MIS-C, also called Pediatric Inflammatory Multisystem Syndrome - PIMS) has been observed amidst the COVID-19 pandemic. MIS-C shares many clinical features with Kawasaki disease (KD), although clinical course and outcomes are divergent. We analyzed whole blood RNA sequences, serum cytokines, and formalin fixed heart tissues from these patients using a computational toolbox of two gene signatures, i.e., the 166-gene viral pandemic (ViP) signature, and its 20-gene severe (s)ViP subset that were developed in the context of SARS-CoV-2 infection and a 13-transcript signature previously demonstrated to be diagnostic for KD. Our analyses revealed that KD and MIS-C are on the same continuum of the host immune response as COVID-19. While both the pediatric syndromes converge upon an IL15/IL15RA-centric cytokine storm, suggestive of shared proximal pathways of immunopathogenesis, they diverge in other laboratory parameters and cardiac phenotypes. The ViP signatures also revealed unique targetable cytokine pathways in MIS-C, place MIS-C farther along in the spectrum in severity compared to KD and pinpoint key clinical (reduced cardiac function) and laboratory (thrombocytopenia and eosinopenia) parameters that can be useful to monitor severity. Naming Human Diseases: Ethical Principles of Curating Exclusive Substitute for Inopportune Nosology BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. 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. The compound of baicalein, wogonin, and oroxylin-A inhibiting EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. Entorhinal cortex epigenome-wide association study highlights four novel loci showing differential methylation in Alzheimer's disease BackgroundStudies on DNA methylation (DNAm) in Alzheimers disease (AD) have recently highlighted several genomic loci showing association with disease onset and progression.

MethodsHere, we conducted an epigenome-wide association study (EWAS) using DNAm profiles in entorhinal cortex (EC) from 149 AD patients and control brains and combined these with two previously published EC datasets by meta-analysis (total n=337).

ResultsWe identified 12 cytosine-phosphate-guanine (CpG) sites showing epigenome-wide significant association with either case-control status or Braaks tau-staging. Four of these CpGs, located in proximity to CNFN/LIPE, TENT5A, PALD1/PRF1, and DIRAS1, represent novel findings. Integrating DNAm levels with RNA sequencing-based mRNA expression data generated in the same individuals showed significant DNAm-mRNA correlations for 6 of the 12 significant CpGs. Lastly, by calculating rates of epigenetic age acceleration using two recently proposed "epigenetic clock" estimators we found a significant association with accelerated epigenetic aging in AD patients vs. controls.

ConclusionIn summary, our study represents the hitherto most comprehensive EWAS in AD using EC and highlights several novel differentially methylated loci with potential effects on gene expression. Low photorespiratory capacity is sufficient to block the induction of photosynthetic capacity in high light The induction of high photosynthetic capacity in high light (HL) is a common response among many herbaceous dicot plants, however, the signals that control this response remain largely unknown. Here, multiple independent lines of evidence are presented in support of the conclusion that low photorespiratory capacity acts a negative signal to limit photosynthetic capacity acclimation in HL in Arabidopsis thaliana. Using a panel of natural accessions, primary nitrogen (N) assimilation and photorespiration rates early after a shift to growth in HL, as well as activities for key enzymes in these pathways, were shown to positively correlate with the magnitude of the subsequent induction of photosynthetic capacity, which occurred several days later. Time-resolved metabolomic data during acclimation to HL were collected using a strongly acclimating ecotype and a weakly acclimating ecotype, revealing in greater detail the differences in N assimilation, photorespiration, and triose-phosphate utilization pathways underlying efficient photosynthetic capacity acclimation. When shifted into HL growth conditions under non-photorespiratory conditions, weakly acclimating ecotypes and even photorespiratory mutants gained the ability to strongly induce high photosynthetic capacity in HL. Thus, a negative, photorespiration-dependent signal early in the HL shift appears to block photosynthetic capacity acclimation in accessions with low photorespiratory capacity, whereas accessions with high photorespiratory capacity are licensed to increase photosynthetic capacity. Climate warming changes synchrony of plants and pollinators Climate warming changes the phenology of many species. When interacting organisms respond differently, climate change may disrupt their interactions and affect the stability of ecosystems. Here, we used GBIF occurrence records to examine phenology trends in plants and their associated insect pollinators in Germany since the 1980s. We found strong phenological advances in plants, but differences in the extent of shifts among pollinator groups. The temporal trends in plant and insect phenologies were generally associated with interannual temperature variation, and thus likely driven by climate change. When examining the synchrony of species-level plant-pollinator interactions, their temporal trends differed among pollinator groups. Overall, plant-pollinator interactions become more synchronized, mainly since the phenology of plants responded more strongly to climate change than that of the pollinators. However, if the observed trends continue, many interactions may become more asynchronous again in the future. Our study suggests that climate change affects the phenologies of both plants and insects, and that it also influences the synchrony of plant-pollinator interactions. Clp protease and antisense RNA jointly regulate the global regulator CarD to mediate mycobacterial starvation response Under starvation conditions, bacteria tend to slow down their translation rate by reducing rRNA synthesis, but the way they accomplish that may vary in different bacteria. In Mycobacterium species, transcription of rRNA is activated by the RNA polymerase (RNAP) accessory transcription factor CarD, which interacts directly with RNAP to stabilize the RNAP-promoter open complex formed on rRNA genes. The functions of CarD have been extensively studied, but the mechanisms that control its expression remain obscure. Here, we report that the level of CarD was tightly regulated when mycobacterial cells switched from nutrient-rich to nutrient-deprived conditions. At the translational level, an antisense RNA of carD (AscarD) was induced in a SigF-dependent manner to bind with carD mRNA and inhibit CarD translation, while at the post-translational level, the residual intracellular CarD was quickly degraded by the Clp protease. AscarD thus worked synergistically with Clp protease to decrease the CarD level to help mycobacterial cells cope with the nutritional stress. Altogether, our work elucidates the regulation mode of CarD and delineates a new mechanism for the mycobacterial starvation response, which is important for the adaptation and persistence of mycobacterial pathogens in the host environment. LabGym: quantification of user-defined animal behaviors using learning-based holistic assessment Quantifying animal behavior is important for many branches of biological research. Current computational tools for quantifying behaviors typically rely on a few pre-defined, simplified features to identify a behavior. However, such an approach underuses the behavioral information and restricts the tools applicability to limited behavior types or animal species. Here we report a new tool, LabGym, for quantifying animal behaviors without such limitations. Combining a novel approach for effective evaluation of animal motion with customizable convolutional recurrent networks for capturing spatiotemporal details, LabGym provides holistic behavioral assessment and accurately identify user-defined animal behaviors without restrictions on behavior types or animal species. It then calculates various parameters to provide quantitative measurements of the intensity and dynamics of each behavior. LabGym requires neither any intermediate step that causes information loss nor programming knowledge from users for post-hoc analysis. It tracks multiple animals simultaneously in various experimental settings for high-throughput and versatile analysis, and provides users a way to generate visualizable behavioral datasets that are valuable resources for research in neuroscience, behavioral science, and computer science. We demonstrate its efficacy in capturing subtle behavioral changes in animals ranging from soft-bodied invertebrates to mammals. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Tracheal tube fusion in Drosophila involves release of extracellular vesicles from multivesicular bodies Extracellular vesicles (EVs) comprise diverse types of cell-released membranous structures that are thought to play important roles in intercellular communication. Despite extensive work on the formation and functions of EVs in cultured cells, studies of EVs in vivo have remained scarce. We report here that EVs are present in the developing lumen of tracheal tubes in Drosophila embryos. We defined two distinct EV subpopulations, one of which contains the Munc13-4 homologue Staccato (Stac) and is spatially and temporally associated with tracheal tube fusion events. The formation of Stac-positive luminal EVs depends on the tip-cell-specific GTPase Arl3, which is also required for the formation of Stac-positive multivesicular bodies, suggesting that Stac-EVs derive from fusion of Stac-MVBs with the luminal membrane in tip cells during anastomosis formation. The GTPases Rab27 and Rab35 cooperate downstream of Arl3 to promote Stac-MVB formation and tube fusion. We propose that Stac-MVBs act as membrane reservoirs that facilitate tracheal lumen fusion in a process regulated by Arl3, Rab27, Rab35, and Stac/Munc13-4. Modulations in motor unit discharge are related to changes in fascicle length during isometric contractions The integration of electromyography (EMG) and ultrasound imaging has provided important information about the mechanisms of muscle activation and contraction. Unfortunately, EMG does not allow an accurate assessment of the interplay between the neural drive received by muscles, changes in fascicle length (FL) and the force/torque produced. We aimed to assess the relationship between modulations in tibialis anterior (TA) motor unit (MU) firing rate, FL and dorsiflexion torque (DT) using ultrasound-transparent high-density EMG electrodes. EMG and ultrasound images were recorded simultaneously from TA, using a 32-electrode silicon matrix, while performing isometric dorsiflexion, at diverse ankle joint positions (0{degrees} and 30{degrees} plantar flexion) and torques (20% and 40% of maximum). EMG signals were decomposed into individual MUs and changes in FL were assessed with a fascicle-tracking algorithm. MU firings were converted into a cumulative spike train (CST) that was cross-correlated with DT (CST-DT) and FL (CST-FL). High cross-correlations were found for CST-FL, 0.60 (range: 0.31-0.85) and CST-DT 0.71 (range: 0.31-0.88). Cross-correlation lags revealed that the delay between CST-FL (~75ms) was significantly smaller than CST-DT (~150ms, p<0.001). These delays affected the interpretation of MU recruitment/de-recruitment thresholds, with FL showing similar lengths for both recruitment and de-recruitment. This study is the first to demonstrate that changes in TA FL are closely related to both modulations in MU firing frequency and DT. These relationships allow assessment of the interplay between neural drive, muscle contraction and resultant torque, thereby providing a better understanding of the mechanisms responsible for the generation of muscle force.

NEW AND NOTEWORTHYBy employing ultrasound-transparent high-density surface EMG electrodes, we show that modulations in tibialis anterior motor unit discharge rate were closely related to both changes in its fascicle length and resultant torque. These relationships allowed quantifying delays between neural drive and muscle shortening as well as muscle shortening and torque during submaximal isometric contractions, providing an accurate estimation of the time required to generate muscle force and subsequent production of torque via the tendon. Autophagosomes fuse to phagosomes and facilitate the degradation of apoptotic cells in Caenorhabditis elegans Autophagosomes are double-membrane intracellular vesicles that degrade protein aggregates, intracellular organelles, and other cellular components. During the development of the nematode Caenorhabditis elegans, many somatic and germ cells undergo apoptosis. These cells are engulfed and degraded by their neighboring cells. We discovered a novel role of autophagosomes in facilitating the degradation of apoptotic cells using a real-time imaging technique. Specifically, the double-membrane autophagosomes in engulfing cells are recruited to the surfaces of phagosomes containing apoptotic cells and subsequently fuse to phagosomes, allowing the inner vesicle to enter the phagosomal lumen. Mutants defective in the production of autophagosomes display significant defects in the degradation of apoptotic cells, demonstrating the importance of autophagosomes to this process. The signaling pathway led by the phagocytic receptor CED-1, the adaptor protein CED-6, and the large GTPase dynamin (DYN-1) promotes the recruitment of autophagosomes to phagosomes. Moreover, the subsequent fusion of autophagosomes with phagosomes requires the functions of the small GTPase RAB-7 and the HOPS complex components. Further observations suggest that autophagosomes provide apoptotic cell-degradation activities in addition to and in parallel of lysosomes. Our findings reveal that, unlike the single-membrane, LC3-associated phagocytosis (LAP) vesicles reported for mammalian phagocytes, the canonical double-membrane autophagosomes facilitate the clearance of C. elegans apoptotic cells. These findings add autophagosomes to the collection of intracellular organelles that contribute to phagosome maturation, identify novel crosstalk between the autophagy and phagosome maturation pathways, and discover the upstream signaling molecules that initiate this crosstalk. Enhancing CDK4/6 inhibitor therapy for medulloblastoma using nanoparticle delivery and scRNA-seq-guided combination with sapanisertib CDK4/6 inhibitors hold promise for brain tumor treatment, but efficacy has been limited by recurrence in both preclinical models and clinical trials. To address recurrence, we tested a nanoparticle formulation of the CDK4/6 inhibitor palbociclib (POx-palbo) in mice genetically-engineered to develop SHH-driven medulloblastoma. We then analyzed medulloblastomas in mice receiving palbociclib treatment, and compared the efficacy of combining palbociclib with specific inhibitors suggested by our analysis. POx-Palbo showed reduced toxicity compared to conventional palbociclib, was tolerable in parenteral administration, improved CNS pharmacokinetics, and extended survival of mice with medulloblastoma. Recurrence, however, remained problematic as fractions of tumor cells proliferated during therapy. ScRNA-seq identified a gene expression pattern unique to proliferating medulloblastoma cells in POx-Palbo-treated mice, marked by up-regulation of the glutamate transporter Slc1a2 and down-regulation of diverse ribosomal genes. Reduced mTORC1 signaling, suggested by ribosomal suppression in POx-Palbo-treated tumors was confirmed by decreased 4EBP1 phosphorylation (p4EBP1). Further reducing mTORC1 activity by combining POx-Palbo with the mTORC1 inhibitor sapanisertib produced mutually enhancing effects, with increased suppression of both pRB and p4EBP1, and prolonged mouse survival compared to either agent alone. In contrast, targeting cell cycle progression by combining POx-Palbo with the SHH-pathway inhibitor vismodegib, or with the replication-targeting agents gemcitabine or etoposide, failed to enhance efficacy. Our data show the potential of nanoparticle formulation and scRNA-seq analysis of resistance to improve brain tumor treatment, and identify POx-palbo plus sapanisertib as effective combinatorial therapy for SHH medulloblastoma. This combination may be appropriate for testing in patients with recurrence, who need new options. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics. SARS-CoV-2 Beta Variant Substitutions Alter Spike Glycoprotein Receptor Binding Domain Structure and Stability The emergence of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) and the subsequent COVID-19 pandemic has visited a terrible cost on the world in the forms of disease, death, and economic turmoil. The rapid development and deployment of extremely effective vaccines against SARS-CoV-2 have seemingly brought within reach the end of the pandemic. However, the virus has acquired mutations; and emerging variants of concern (VOC) are more infectious and reduce the efficacy of existing vaccines. While promising efforts to combat these variants are underway, the evolutionary pressures leading to these variants are poorly understood. To that end, here we have studied the effects on the structure and function of the SARS-CoV-2 spike glycoprotein receptor-binding domain of three amino-acid substitutions found in several variants of concern, including alpha (B.1.1.7), beta (B.1.351), and gamma (P.1). We found that these substitutions alter the RBD structure, stability, and ability to bind to ACE2, in such a way as to possibly have opposing and compensatory effects. These findings provide new insights into how these VOC may have been selected for infectivity while maintaining the structure and stability of the receptor binding domain. Bias correction for inverse variance weighting Mendelian randomization BackgroundInverse-variance weighted two-sample Mendelian randomization (IVW-MR) is the most widely used approach that utilizes genome-wide association studies (GWAS) summary statistics to infer the existence and the strength of the causal effect between an exposure and an outcome. Estimates from this approach can be subject to different biases due to the use of weak instruments and winners curse, which can change as a function of the overlap between the exposure and outcome samples.

MethodsWe developed a method (MRlap) that simultaneously considers weak instrument bias and winners curse, while accounting for potential sample overlap. Assuming spike-and-slab genomic architecture and leveraging LD-score regression and other techniques, we could analytically derive, reliably estimate, and hence correct for the bias of IVW-MR using association summary statistics only.

ResultsWe tested our approach using simulated data for a wide range of realistic settings. In all the explored scenarios, our correction reduced the bias, in some situations by as much as 30 folds. When applied to real data on obesity-related exposures, we observed significant differences between IVW-based and corrected effects, both for non-overlapping and fully overlapping samples. Additionally, our results are consistent with the fact that the strength of the biases will decrease as the sample size increases. Finally, we also showed that the overall bias is also dependent on the genetic architecture of the exposure, and traits with low heritability and/or high polygenicity are more strongly affected.

ConclusionsOur method not only reduces bias in causal effect estimation but also enables the use of much larger GWAS sample sizes, by allowing for potentially overlapping samples. Characterization of DNA extrusion scenarios underlying aggregate Hi-C patterns by a new computational framework 1Aggregate Hi-C maps centered on several genetic elements show significant DNA extrusion features. However, these features are missing from the component Hi-C images. To resolve this discrepancy, we present a new computational framework that reconstructs DNA extrusion scenarios (DESs) and explains the enhancement of DNA extrusion features by aggregation. We showed that the reconstructed DES is essential to understand the site interactions of the relevant genetic elements. For example, CCCTC-binding factor (CTCF) binding elements (CBE) with a specific motif direction are only affected by DNA extrusion activities initiated from one side, manifested as a directional bias. About half of promoters search for enhancers in a specific direction, as indicated by the corresponding DES. Moreover, since DNA extrusion does not pause at promoter or enhancer elements, direct promoter-enhancer contacts are not long-lasting; rather, promoters and enhancers are brought in proximity by continuous, dynamic DNA extrusion. 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. SARS-CoV-2 RdRp uses NDPs as a substrate and is able to incorporate NHC into RNA from diphosphate form molnupiravir The coronavirus disease 2019 (COVID-19) has been ravaging throughout the world for more than two years and has severely impaired both human health and the economy. The causative agent, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) employs the viral RNA-dependent RNA polymerase (RdRp) complex for genome replication and transcription, making RdRp an appealing target for antiviral drug development. Here, we reveal that RdRp can recognize and utilize nucleoside diphosphates (NDPs) as a substrate to synthesize RNAs with an efficiency of about two thirds of using nucleoside triphosphates (NTPs) as a substrate. NDPs incorporation is also template-specific and has high fidelity. Moreover, RdRp can incorporate {beta}-d-N4-hydroxycytidine (NHC) into RNA while using diphosphate form molnupiravir (MDP) as a substrate. We also observed that MDP is a better substrate for RdRp than the triphosphate form molnupiravir (MTP). Collective polymerase dynamics emerge from DNA supercoiling during transcription All biological processes ultimately come from physical interactions. The mechanical properties of DNA play a critical role in transcription. RNA polymerase can over or under twist DNA (referred to as DNA supercoiling) when it moves along a gene resulting in mechanical stresses in DNA that impact its own motion and that of other polymerases. For example, when enough supercoiling accumulates, an isolated polymerase halts and transcription stops. DNA supercoiling can also mediate non-local interactions between polymerases that shape gene expression fluctuations. Here, we construct a comprehensive model of transcription that captures how RNA polymerase motion changes the degree of DNA supercoiling which in turn feeds back into the rate at which polymerases are recruited and move along the DNA. Surprisingly, our model predicts that a group of three or more polymerases move together at a constant velocity and sustain their motion (forming what we call a polymeton) whereas one or two polymerases would have halted. We further show that accounting for the impact of DNA supercoiling on both RNA polymerase recruitment and velocity recapitulates empirical observations of gene expression fluctuations. Finally, we propose a mechanical toggle switch whereby interactions between genes are mediated by DNA twisting as opposed to proteins. Understanding the mechanical regulation of gene expression provides new insights into how endogenous genes can interact and informs the design of new forms of engineered interactions.

PACS numbers: Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, LK-DFBA, that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this frameworks success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that the optimal constraint approach was generally robust to local perturbations of the system, indicating that just a single wild-type dataset could allow identification of the ideal constraint for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems.