Vesalius: high-resolution in silico anatomization of Spatial Transcriptomic data using Image Analysis Characterization of tissue architecture promises to deliver insights into development, cell communication and disease. In silico spatial domain retrieval methods have been developed for spatial transcriptomics (ST) data assuming transcriptional similarity of neighboring barcodes. However, domain retrieval approaches with this assumption cannot work in complex tissues composed of multiple cell types. This task becomes especially challenging in cellular resolution ST methods. We developed Vesalius to decipher tissue anatomy from ST data by applying image processing technology. Vesalius uniquely detected territories composed of multiple cell types and successfully recovered tissue structures in high-resolution ST data including in mouse brain, embryo, liver, and colon. Utilizing this tissue architecture, Vesalius identified tissue morphology specific gene expression and regional specific gene expression changes for astrocytes, interneuron, oligodendrocytes, and entorhinal cells in the mouse brain.
Revisiting regulatory decoherence: Accounting for temporal bias in a co-expression analysis reveals novel candidates controlling environmental response Environment can alter the degree of phenotypic variation and covariation, potentially influencing evolutionary trajectories. However, environment-driven changes in phenotypic variation remain understudied. In an effort to exploit the abundance of RNASequencing data now available, an increasing number of ecological studies rely on population-level correlation to characterize the plastic response of the entire transcriptome and to identify environmentally responsive molecular pathways. These studies are fundamentally interested in identifying groups of genes that respond in concert to environmental shifts. We show that population-level differential co-expression exhibits biases when capturing changes of regulatory activity and strength in rice plants responding to elevated temperature. One possible cause of this bias is regulatory saturation, the observation that detectable co-variance between a regulator and its target may be low as their transcript abundances are induced. This phenomenon appears to be particularly acute for rapid-onset environmental stressors. However, our results suggest that temporal correlations may be a reliable means to detect transient regulatory activity following rapid onset environmental perturbations such as temperature stress. Such temporal bias is likely to confound the studies of phenotypic integration, where high-order organismal traits are hypothesized to be more integrated with strong correlation under stressful conditions, while recent transcriptome studies exhibited weaker coexpression between genes under stressful conditions. Collectively, our results point to the need to account for the nuances of molecular interactions and the possibly confounding effects that these can introduce into conventional approaches to study transcriptome datasets.
A chemical screen based on an interruption of zebrafish gastrulation identifies the HTR2C inhibitor Pizotifen as a suppressor of EMT-mediated metastasis Metastasis is responsible for approximately 90% of cancer-associated mortality but few models exist that allow for rapid and effective screening of anti-metastasis drugs. Current mouse models of metastasis are too expensive and time consuming to use for rapid and high-throughput screening. Therefore, we created a unique screening concept utilizing conserved mechanisms between zebrafish gastrulation and cancer metastasis for identification of potential anti-metastatic drugs. We hypothesized that small chemicals that interrupt zebrafish gastrulation might also suppress metastatic progression of cancer cells and developed a phenotype-based chemical screen to test the hypothesis. The screen used epiboly, the first morphogenetic movement in gastrulation, as a marker and enabled 100 chemicals to be tested in five hours. The screen tested 1280 FDA-approved drugs and identified Pizotifen, an antagonist for serotonin receptor 2C (HTR2C) as an epiboly-interrupting drug. Pharmacologic and genetic inhibition of HTR2C suppressed metastatic progression in a mouse model. Blocking HTR2C with Pizotifen restored epithelial properties to metastatic cells through inhibition of Wnt-signaling. In contrast, HTR2C induced epithelial to mesenchymal transition (EMT) through activation of Wnt-signaling and promoted metastatic dissemination of human cancer cells in a zebrafish xenotransplantation model. Taken together, our concept offers a novel platform for discovery of anti-metastasis drugs. The compound of baicalein, wogonin, and oroxylin-A inhibits 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. Pathways of savannization in a mesic African savanna-forest mosaic O_LIFires in savannas limit tree cover, thereby promoting flammable grass accumulation and fuelling further frequent fires. Meanwhile, forests and thickets form dense canopies that reduce C4-grass fuel loads and creating a humid microclimate, thereby excluding fires under typical climatic conditions.
C_LIO_LIHowever, extreme fires occasionally burn into these closed-canopy systems. Although these rare fires cause substantial tree mortality and can make repeat fires more likely, the long-term consequences of an extreme fire for closed canopy vegetation structure and potential to convert to savanna (hereafter "savannization") remain largely unknown.
C_LIO_LIHere, we analysed whether an extreme fire could, alone, alter species composition, vegetation structure, and fire regimes of closed-canopy ecosystems in an intact savanna-forest-thicket mosaic, or whether successive fires after an initial extreme fire were necessary to trigger a biome transition between from forest to savanna.
C_LIO_LIWe found that forests that only burned once recovered, whereas those that burned again following an initial extreme fire transitioned from closed-canopy forests towards open, grassy savannas.
C_LIO_LIWhile thickets had less tree mortality in fires than forests, repeat fires nonetheless precipitated a transition towards savannas.
C_LIO_LIColonization of the savanna tree community lagged behind the grass community, but also began to transition.
C_LI

SynthesisOur results suggest that rare extreme fires, followed by repeated burning can indeed result in savannization in places where savanna and forest represent alternative stable states. 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. Bulk-assembly of monodisperse coacervates and giant unilamellar vesicles with programmable hierarchical complexity Compartmentalization is an essential step for the emergence of protocellular life on the early earth and the endeavor to bottom-up construct synthetic cells in the lab. Among manifold natural strategies for compartmentalization, assembly of lipid vesicles or coacervate droplets with resemblance to cell envelope or cytoplasm have emerged as two dominant paradigms. However, the spontaneous assembly of lipid vesicles or coacervates generally results in polydisperse microcompartments and lacks the ability to efficiently integrate functional building blocks towards defined complexity, which give rise to the difficulty to achieve the order and complexity towards the scope of that of life via spontaneous assembly. Herein, we show that the interplay of coacervates and colloidal particles provides a bulk-assembly approach to form monodisperse coacervates and giant unilamellar vesicles (GUVs) with programmable hierarchical complexity. The hierarchical complexity of the synthetic cell microcompartments can be delicately engineered from the control of the spatial or temporal organization of components in entity unit to manifold forms of interconnected synthetic cell consortia, achieving a compartment hierarchy analogous to the organelle-cell-tissue structure. The increase in complexity gives rise to emergent properties, for instance, collective morphology evolution of coacervate assemblages, gated permeability of GUVs without sophisticated protein machinery, and remarkable structural and functional stability of GUVs at extreme conditions. This work paves an unprecedented step to form monodisperse, hierarchical coacervates and GUVs via bulk-assembly, which provides a prebiotically plausible approach for the proper integration of inanimate matter towards the order and complexity of life and monodisperse, hierarchical microdroplets for follow-up scientific and technological applications. mRNA vaccine-induced T cells respond identically to SARS-CoV-2 variants of concern but differ in longevity and homing properties depending on prior infection status While mRNA vaccines are proving highly efficacious against SARS-CoV-2, it is important to determine how booster doses and prior infection influence the immune defense they elicit, and whether they protect against variants. Focusing on the T cell response, we conducted a longitudinal study of infection-naive and COVID-19 convalescent donors before vaccination and after their first and second vaccine doses, using a high-parameter CyTOF analysis to phenotype their SARS-CoV-2-specific T cells. Vaccine-elicited spike-specific T cells responded similarly to stimulation by spike epitopes from the ancestral, B.1.1.7 and B.1.351 variant strains, both in terms of cell numbers and phenotypes. In infection-naive individuals, the second dose boosted the quantity and altered the phenotypic properties of SARS-CoV-2-specific T cells, while in convalescents the second dose changed neither. Spike-specific T cells from convalescent vaccinees differed strikingly from those of infection-naive vaccinees, with phenotypic features suggesting superior long-term persistence and ability to home to the respiratory tract including the nasopharynx. These results provide reassurance that vaccine-elicited T cells respond robustly to emerging viral variants, confirm that convalescents may not need a second vaccine dose, and suggest that vaccinated convalescents may have more persistent nasopharynx-homing SARS-CoV-2-specific T cells compared to their infection-naive counterparts. Characterization of methylation profiles in spontaneous preterm birth placental villous tissue Preterm birth is a global public health crisis which results in significant neonatal and maternal mortality. Yet little is known regarding the molecular mechanisms of idiopathic spontaneous preterm birth, and we have few diagnostic markers for adequate assessment of placental development and function. Previous studies of placental pathology and our transcriptomics studies suggest a role for placental maturity in idiopathic spontaneous preterm birth. It is known that placental DNA methylation changes over gestation. We hypothesized that if placental hypermaturity is present in our samples, we would observe a unique idiopathic spontaneous preterm birth DNA methylation profile potentially driving the gene expression differences we previously identified in our placental samples. Our results indicate the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern suggesting hypermaturity. Only seven significant differentially methylated regions fitting the idiopathic spontaneous preterm birth specific (relative to the controls) profile were identified, indicating unusually high similarity in DNA methylation between idiopathic spontaneous preterm birth and term birth samples. We identified an additional 1,718 significantly methylated regions in our gestational age matched controls were the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern, again indicating a striking level of similarity between the idiopathic spontaneous preterm birth and term birth samples. Pathway analysis of these regions revealed differences in genes within the WNT and Cadherin signaling pathways, both of which are essential in placental development and maturation. Taken together, these data demonstrate that the idiopathic spontaneous preterm birth samples are molecularly more mature than expected given their respective gestational age which likely impacts birth timing. Exploratory analysis of multiple traits co-adaptations in the population history During the history of range expansion, the populations encounter with variety of environments. They respond to the local environments by modifying the mutually interacting traits. Therefore, to understand the whole life history of the populations, it is ideal to capture the history of their range expansion with reference to the series of surrounding environments and to infer the coadaptation of the multiple traits. Toward this end, we provide an exploratory analysis based on the features of populations: site frequency spectra of populations, population-specific FST, association between genes and environments, positive selections on traits mapped on the admixture graph, and GWAS results. Correspondence analysis of genes, environments, and traits provides a birds-eye view of the history of population differentiation and range expansion and various types of environmental selections at the times. Principal component analysis of the estimated trait-specific polygenic adaptations mapped on the admixture graph enables to understand the coadaptation of multiple traits. The potential usefulness was confirmed by analyzing a public dataset of wild poplar in northwestern America. In response to the northern cold temperature and longer daylength, the populations increased the photosynthetic activity and nutrient use efficiency at the expense of the risk of pathogen invasion, and in response to warm temperature, they increased the growth. At higher altitude, they shifted the maximum activity to earlier period in spring to reduce the activity in dry summer. The R codes for our representation method and simulations of population colonization used in this study are available as supplementary script. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Differences in gait stability and acceleration characteristics between healthy young and older females Our aim was to evaluate differences in gait acceleration intensity, variability and stability of feet and trunk between older females and young females using inertial sensors. Twenty older females (OF; mean age 68.4, SD 4.1 years) and eighteen young females (YF; mean age 22.3, SD 1.7 years) were asked to walk straight for 100 meters at their preferred speed, while wearing inertial sensors on heels and lower back. We calculated spatiotemporal measures, foot and trunk acceleration characteristics and their variability, as well as trunk stability using the local divergence exponent (LDE). Two-way analysis of variance (including the factors foot and age), Students t-test, and Mann-Whitney U test were used to compare statistical differences of measures between groups. Cohens d effects were calculated for each variable. Foot maximum vertical acceleration and amplitude, trunk-foot vertical acceleration attenuation, as well as their variability were significantly smaller in OF than in YF. In contrast, trunk mediolateral acceleration amplitude, maximum vertical acceleration, and amplitude, as well as their variability were significantly larger in OF than in YF. Moreover, OF showed lower stability (i.e., higher LDE values) in mediolateral acceleration, mediolateral and vertical angular velocity of the trunk. Even though we measured healthy older females, these participants showed lower vertical foot accelerations with higher vertical trunk acceleration, lower trunk-foot vertical acceleration attenuation, less gait stability, and more variability of the trunk, and hence, were more likely to fall. These findings suggest that instrumented gait measurements may help for early detection of changes or impairments in gait performance, even before this can be observed by clinical eye or gait speed. Feed efficiency of lactating Holstein cows was not as repeatable across diets as within diet over subsequent lactation stages BackgroundImproving feed efficiency has become a common target for dairy farmers to meet the requirement of producing more milk with fewer resources. To improve feed efficiency, a prerequisite is to ensure that the cows identified as most or least efficient will remain as such, independently of diet composition. Therefore, the current research analysed the ability of lactating dairy cows to maintain their feed efficiency while changing the energy density of the diet by changing its concentration in starch and fibre. A total of 60 lactating Holstein cows, including 33 primiparous cows, were first fed a high starch diet-low fibre (diet S+F-), then switched over to a low starch diet-high fibre (diet S-F+). To know if diet affect feed efficiency, we compared the ability of feed efficiency to be maintained within a diet over subsequent lactation stages, known as repeatability of feed efficiency, with its ability to be maintained across diets, known as reproducibility of feed efficiency. To do so we used two indicators: the estimation of the error of repeatability/reproducibility, which is commonly used in metrology, and the coefficient of correlation of concordance (CCC), which is used in biology. The effect of diet change could also lead to a change in cows sorting behaviour which could potentially affect feed efficiency if for example the most efficient cows select more concentrate than the least efficient. We therefore analysed the relationship between the differences in individual feed refusals composition and the differences in feed efficiency. To do so, the composition of each feed refusal was described with its near infra-red (NIR) spectroscopy and was performed on each individual feed ingredient, diet and refusals and used as composition variable. The variability of the NIR spectra of the refusals was described with its principal components thanks to a principal component analysis (PCA). The Pearson correlation was estimated to check the relationship between feed efficiency and refusals composition, i.e. sorting behaviour.

ResultsThe error of reproducibility of feed efficiency across diets was 2.95 MJ/d. This error was significantly larger than the errors of repeatability estimated within diet, which were 2.01 MJ/d within diet S-F+ and 2.40 MJ/d within diet S+F-. The CCC was 0.64 between feed efficiency estimated within diet S+F- and feed efficiency estimated within diet S-F+. This CCC was smaller than the one observed for feed efficiency estimated within diet between two subsequent lactation stages (CCC = 0.72 within diet S+F- and 0.85 within diet S-F+). Feed efficiency was poorly correlated to the first two principal components, which explained 90% of the total variability of the NIR spectra of the individual refusals. This suggests that feed sorting behaviour did not explain differences in feed efficiency.

ConclusionsFeed efficiency was significantly less reproducible across diets than repeatable within the same diet over subsequent lactation stages, but cows ranking for feed efficiency was not significantly affected by diet change. This loss in repeatability across diets could be due to a more pronounced feed sorting subsequent to the change in diet composition. However, the differences in sorting behaviour between cows were not associated to feed efficiency differences in this trial neither with the S+F- diet nor with the S-F+ diet. Those results have to be confirmed on diets having different forage to concentrate ratios to ensure that the least and most efficient cows will not change. A multi-state dynamic process confers mechano-adaptation to the bacterial flagellar motor Adaptation is a defining feature of living systems. The bacterial flagellar motor adapts to changes in the external mechanical environment by adding or removing torque-generating stator units. However, the molecular mechanism for mechanosensitive motor remodeling remains unclear. Here, we induced stator disassembly using electrorotation, followed by the time-dependent assembly of the individual stator units into the motor after electrorotation was terminated. From these experiments, we extracted detailed statistics of the dwell times that comprise the stochastic dynamics of the binding and unbinding of stator units. Dwell times reveal multiple timescales, indicating the existence of multiple binding states of the stator units. Based on these results, we propose a minimal model in which the stator unit can occupy four different states - two bound states with very different rates of unbinding, a diffusive unbound state, and a transiently detached state. Our minimal model quantitatively explains multiple features of the experimental data and allows us to determine the transition rates among all four states. Our experiments and modeling suggest a mechanism of mechano-adaptive remodeling of the bacterial flagellar motor in which torque generated by bound stator units controls their effective unbinding rate by modulating the transition between the two bound states. Furthermore, the binding rate of stator units with the motor has a non-monotonic dependence on the number of bound units, likely because of two counteracting effects of motor rotation on the binding process. Identifying and correcting multiple sources of misspecification in GWAS summary statistics for polygenic scores Publicly available genome-wide association studies (GWAS) summary statistics exhibit uneven quality, which can impact the validity of follow-up analyses. First, we present an overview of possible misspecifications that come with GWAS summary statistics. Then, in both simulations and real data analyses, we show that additional information such as imputation INFO scores, allele frequencies, and per-variant sample sizes in GWAS summary statistics can be used to detect possible issues and correct for misspecifications in the GWAS summary statistics. One important motivation for us is to improve the predictive performance of polygenic scores built from these summary statistics. Unfortunately, due to the lack of reporting standards for GWAS summary statistics, this additional information is not systematically reported. We also show that using well-matched LD references can improve model fit and translate into more accurate prediction. Finally, we discuss how to make polygenic score methods such as lassosum and LDpred2 more robust to these misspecifications to improve their predictive power. Lipid membranes modulate the activity of RNA through sequence-dependent interactions RNA is a ubiquitous biomolecule that can serve as both catalyst and information carrier. Understanding how RNA bioactivity is controlled is crucial for elucidating its physiological roles and potential applications in synthetic biology. Here we show that lipid membranes can act as RNA organization platforms, introducing a novel mechanism for ribo-regulation. The activity of R3C ribozyme can be modified by the presence of lipid membranes, with direct RNA-lipid interactions dependent on RNA nucleotide content, base pairing and length. In particular, the presence of guanine in short RNAs is crucial for RNA-lipid interactions, and G-quadruplex formation further promotes lipid binding. Lastly, by artificially modifying the R3C substrate sequence to enhance membrane binding we generated a lipid-sensitive ribozyme reaction with riboswitch-like behavior. These findings introduce RNA-lipid interactions as a tool for developing synthetic riboswitches and novel RNA-based lipid biosensors, and bear significant implications for RNA World scenarios for the origin of life. A small molecule antagonizes jasmonic acid perception and auxin responses in vascular and nonvascular plants The phytohormone JA-Ile regulates many stress responses and developmental processes in plants. A co-receptor complex formed by the F-box protein COI1 (Coronatine Insensitive 1) and a JAZ (Jasmonate ZIM-domain) repressor perceives the hormone. JA-Ile antagonists are invaluable tools for exploring the role of JA-Ile in specific tissues and developmental stages, and for identifying regulatory processes of the signalling pathway. Using two complementary chemical screens, we identified three compounds that exhibit a robust inhibitory effect on both the hormone-mediated COI-JAZ interaction and degradation of JAZ1 and JAZ9 in vivo. One molecule, J4, also restrains specific JA-induced physiological responses in different angiosperm plants, including JA-mediated gene expression, growth inhibition, chlorophyll degradation and anthocyanin accumulation. Interaction experiments with purified proteins indicate that J4 directly interferes with the formation of the Arabidopsis (Arabidopsis thaliana) COI1-JAZ complex otherwise induced by JA. The antagonistic effect of J4 on COI1- JAZ also occurs in the liverwort Marchantia polymorpha, suggesting the mode of action is conserved in land plants. Besides JA signalling, J4 works as an antagonist of the closely-related auxin signalling pathway, preventing TIR1/Aux-IAA interaction and auxin responses in planta, including hormone-mediated degradation of an auxin repressor, gene expression and gravitropic response. However, J4 does not affect other hormonal pathways. Altogether, our results show that this dual antagonist competes with JA-Ile and auxin, preventing the formation of phylogenetically related receptor complexes. J4 may be a useful tool to dissect both the JA-Ile and auxin pathways in particular tissues and developmental stages since it reversibly inhibits these pathways.

One sentence summaryA chemical screen identified a molecule that antagonizes jasmonate perception by directly interfering with receptor complex formation in phylogenetically distant vascular and non-vascular plants. A genomic and morphometric analysis of alpine bumblebees: ongoing reductions in tongue length but no clear genetic component Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyze population structure in these populations. We generated a genome assembly of B. balteatus. We then analyzed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species ("incognitus"). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait likely has a complex genetic and environmental basis. Testing for phylogenetic signal in single-cell RNA-seq data AO_SCPLOWBSTRACTC_SCPLOWPhylogenetic methods are emerging as a useful tool to understand cancer evolutionary dynamics, including tumor structure, heterogeneity, and progression. Most currently used approaches utilize either bulk whole genome sequencing (WGS) or single-cell DNA sequencing (scDNA-seq) and are based on calling copy number alterations and single nucleotide variants (SNVs). Here we explore the potential of single-cell RNA sequencing (scRNA-seq) to reconstruct cancer evolutionary dynamics. scRNA-seq is commonly applied to explore differential gene expression of cancer cells throughout tumor progression. The method exacerbates the single-cell sequencing problem of low yield per cell with uneven expression levels. This accounts for low and uneven sequencing coverage and makes SNV detection and phylogenetic analysis challenging. In this paper, we demonstrate for the first time that scRNA-seq data contains sufficient evolutionary signal and can be utilized in phylogenetic analyses. We explore and compare results of such analyses based on both expression levels and SNVs called from scRNA-seq data. Both techniques are shown to be useful for reconstructing phylogenetic relationships between cells, reflecting the clonal composition of a tumor. Both standardized expression values and SNVs appear to be equally capable of reconstructing a similar pattern of phylogenetic relationship. This pattern is stable even when phylogenetic uncertainty is taken in account. Our results open up a new direction of somatic phylogenetics based on scRNA-seq data. Further research is required to refine and improve these approaches to capture the full picture of somatic evolutionary dynamics in cancer. Emergent dynamics of a three-node regulatory network explain phenotypic switching and heterogeneity: a case study of Th1/Th2/Th17 cell differentiation Naive helper (CD4+) T-cells can differentiate into distinct functional subsets including Th1, Th2, and Th17 phenotypes. Each of these phenotypes has a  master regulator - T-bet (Th1), GATA3 (Th2) and ROR{gamma}T (Th17) - that inhibits the other two master regulators. Such mutual repression among them at a transcriptional level can enable multistability, giving rise to six experimentally observed phenotypes - Th1, Th2, Th17, hybrid Th/Th2, hybrid Th2/Th17 and hybrid Th1/Th17. However, the dynamics of switching among these phenotypes, particularly in the case of epigenetic influence, remains unclear. Here, through mathematical modeling, we investigated the coupled transcription-epigenetic dynamics in a three-node mutually repressing network to elucidate how epigenetic changes mediated by any  master regulator can influence the transition rates among different cellular phenotypes. We show that the degree of plasticity exhibited by one phenotype depends on relative strength and duration of mutual epigenetic repression mediated among the master regulators in a three-node network. Further, our model predictions can offer putative mechanisms underlying relatively higher plasticity of Th17 phenotype as observed in vitro and in vivo. Together, our modeling framework characterizes phenotypic plasticity and heterogeneity as an outcome of emergent dynamics of a three-node regulatory network, such as the one mediated by T-bet/GATA3/ROR{gamma}T. EGFR-driven profibrotic signals in AKI-to-CKD transition are blocked by TNF inhibition: Opportunities for etanercept treatment Proximal-tubule-cell EGFR activation mediates tubule cell proliferation and repair early after kidney injury, while sustained EGFR activation causes kidney fibrosis. Inflammation is a key driver of AKI-to-CKD transition and fibrosis, but mechanisms of EGFR-driven profibrotic responses are not well understood. In a mouse model of AKI-to-CKD transition and CKD progression, we show that EGFR-inhibition significantly reduced kidney expression of many immunoregulatory molecules already by day two after injury, including the potent inflammatory cytokine tumor-necrosis-factor (TNF). Single nuclei RNA-sequencing analysis showed that macrophages were among the main early cellular sources of TNF in the injured kidney. In vitro, EGFR activation in macrophages increased macrophage TNF expression, while EGFR inhibition in vivo reduced kidney macrophage accumulation, as early as two days after injury. Thus, profibrotic EGFR signaling increases kidney TNF levels both directly and indirectly. TNF inhibition did not alter tubule EGFR activation and in contrast to EGFR inhibition did not reduce early macrophage accumulation in the kidney, suggesting that TNF does not promote early infiltration of immune cells in the kidney, but rather regulates profibrotic functions of kidney and/or immune cells. TNF inhibition with etanercept in vivo in AKI-injured mice downregulated a number of cytokines, including TNF itself. This cytokine downregulation overlapped with EGFR inhibition, and additional non-overlapping downregulated cytokines shared the same functions, predicting that TNF inhibition would prevent AKI-to-CKD transition like EGFR inhibition. Indeed, TNF-inhibition with etanercept reduced AKI-induced fibrosis to the same degree as the EGFR-inhibition, while the combination of both treatments showed no additive effect. In conclusion, our results identify TNF as a downstream effector of profibrotic EGFR activation and motivate the examination of TNF pathway inhibition in human AKI or CKD.

Translational statementProximal tubule epidermal-growth-factor (EGFR) activation in mice and likely also in humans drives inflammation in the kidney after acute-kidney-injury (AKI) and causes chronic-kidney-disease (CKD) with fibrosis in a process termed AKI-to-CKD transition. Recent retrospective data shows that patients treated with TNF inhibitors show decreased incidence and progression of CKD. Our current work shows that TNF inhibition in mice is equally as effective as EGFR inhibition in preventing AKI-to-CKD transition and fibrosis. Thus, our results have translational potential and may stimulate examination of short-term TNF inhibition in AKI to prevent AKI-to-CKD transition or possibly of longer-term TNF inhibition in CKD to prevent CKD progression. 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. An image reconstruction framework for characterizing initial visual encoding We developed an image-computable observer model of the initial visual encoding that operates on natural image input, based on the framework of Bayesian image reconstruction from the excitations of the retinal cone mosaic. Our model extends previous work on ideal observer analysis and evaluation of performance beyond psychophysical discrimination, takes into account the statistical regularities of the visual environment, and provides a unifying framework for answering a wide range of questions regarding the visual front end. Using the error in the reconstructions as a metric, we analyzed variations of the number of different photoreceptor types on human retina as an optimal design problem. In addition, the reconstructions allow both visualization and quantification of information loss due to physiological optics and cone mosaic sampling, and how these vary with eccentricity. Furthermore, in simulations of color deficiencies and interferometric experiments, we found that the reconstructed images provide a reasonable proxy for modeling subjects percepts. Lastly, we used the reconstruction-based observer for the analysis of psychophysical threshold, and found notable interactions between spatial frequency and chromatic direction in the resulting spatial contrast sensitivity function. Our method is widely applicable to experiments and applications in which the initial visual encoding plays an important role. Adolescent Alcohol Binge-Drinking induces delayed appearance of behavioral defects in mice Adolescence is a developmental period characterized by significant changes in brain architecture and behavior. The immaturity of the adolescent brain is associated with heightened vulnerability to exogenous agents, including alcohol. Alcohol is the most consumed drug among teenagers, and binge-drinking during adolescence is a major public health concern. Studies have suggested that adolescent alcohol exposure (AAE) may interfere with the maturation of frontal brain regions and lead to long-lasting behavioral consequences. In this study, we used a mouse model of AAE in which adolescent mice reach high blood alcohol concentration after voluntary binge-drinking. In order to assess short- and long-term consequences of AAE, a battery of behavioral tests was performed during late adolescence and during adulthood. We showed that AAE had no short-term effect on young mice behavior but rather increased anxiety- and depressive-like behaviors, as well as alcohol consumption during adulthood. Moreover, alcohol binge-drinking during adolescence dramatically decreased recognition memory performances and behavioral flexibility in both adult males and females. Furthermore, we showed that voluntary consumption of alcohol during adolescence did not trigger any major activation of the innate immune system in the prefrontal cortex (PFC). Together, our data suggest that voluntary alcohol binge-drinking in adolescent mice induces a delayed appearance of behavioral impairments in adulthood. Non-targeted multimodal metabolomics and lipidomics data from ovine rumen fluid fractions BackgroundMetabolomics is a powerful and sensitive approach for investigating low molecular weight metabolite profiles present in rumen biofluids to identify potential roles of metabolites in rumen microbiome and understanding host-level regulatory mechanisms associated with animal production.

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

ConclusionsUntargeted metabolomics provides a detailed snapshot of the ovine ruminal fluid metabolome that can be used as reference for future studies on ovine rumen fluid or as a comparator for other ruminant species. All data and metadata are available for download in the MetaboLights database. O-specific antigen-dependent surface hydrophobicity mediates aggregate assembly type in Pseudomonas aeruginosa Bacteria live in spatially organized aggregates during chronic infections, where they adapt to the host environment, evade immune responses and resist therapeutic interventions. Although it is known that environmental factors such as polymers influence bacterial aggregation, it is not clear how bacterial adaptation during chronic infection impacts the formation and spatial organization of aggregates in the presence of polymers. Here we show that in an in vitro model of cystic fibrosis (CF) containing the polymers eDNA and mucin, O-specific antigen is a major factor in determining the formation of two distinct aggregate assembly types of Pseudomonas aeruginosa due to alterations in cell surface hydrophobicity. Our findings suggest that during chronic infection, interplay between cell surface properties and polymers in the environment may influence the formation and structure of bacterial aggregates, which would shed new light on the fitness costs and benefits of O-antigen production in environments such as CF lungs.

ImportanceDuring chronic infection, several factors contribute to the biogeography of microbial communities. Heterogeneous populations of Pseudomonas aeruginosa form aggregates in cystic fibrosis airways, however, the impact of this population heterogeneity on spatial organization and aggregate assembly is not well understood. In this study we found that changes in O-specific antigen determine the spatial organization of P. aeruginosa cells by altering the relative cell surface hydrophobicity. This finding suggests a role for O-antigen in regulating P. aeruginosa aggregate size and shape in cystic fibrosis airways.