A computational modeling approach for predicting multicell patterns based on signaling-induced differential adhesion Differential adhesion within cell populations enables the emergence of unique patterns in heterogeneous multicellular systems. This process has been previously explored using synthetically engineered heterogenous multicell spheroid systems, in which cell subpopulations engage in bidirectional intercellular signaling to regulate the expression of different cadherins. While engineered cell systems provide excellent experimental tools to observe pattern formation in cell populations, computational models may be leveraged to more systematically explore the key parameters that drive the emergence of different patterns. We developed and validated two- and three-dimensional agent-based models (ABMs) of spheroid patterning for cells engineered with a bidirectional signaling circuit that regulates N- and P-cadherin expression. The model was used to predict how varying initial cell seedings, cadherin induction probabilities, or homotypic adhesion strengths between cells impact spheroid patterning, and unsupervised machine learning techniques were used to map system parameters to unique spheroid patterns. The model was then deployed to design new synthetic cell signaling circuits based on a desired final multicell pattern.
Functional buffering via cell-specific gene expression promotes tissue homeostasis and cancer robustness Functional buffering that ensures biological robustness is critical for maintaining tissue homeostasis, organismal survival, and evolution of novelty. However, the mechanism underlying functional buffering, particularly in multicellular organisms, remains largely elusive. Here, we developed an inference index (C-score) for Cell-specific Expression- BUffering (CEBU), whereby functional buffering is mediated via expression of buffering genes in specific cells and tissues in humans. By computing C-scores across 684 human cell lines using genome-wide CRISPR screens and transcriptomic RNA-seq, we report that C- score-identified putative buffering gene pairs are enriched for members of the same duplicated gene family, pathway, and protein complex. Furthermore, CEBU is especially prevalent in tissues of low regenerative capacity (e.g., bone and neuronal tissues) and is weakest in highly regenerative blood cells, linking functional buffering to tissue regeneration. Clinically, the buffering capacity enabled by CEBU can help predict patient survival for multiple cancers. Our results reveal CEBU as a buffering mechanism contributing to tissue homeostasis and cancer robustness in humans.

Summary blurbWe unveil a genome-wide functional buffering mechanism, termed Cell-specific Expression Buffering (CEBU), whereby gene expression contributes to functional buffering in specific cell types and tissues. We link CEBU to genetic interactions, tissue homeostasis and cancer robustness.
Hypertrophic Chondrocytes Serve as a Reservoir for Marrow Associated Skeletal Stem and Progenitor Cells, Osteoblasts, and Adipocytes During Skeletal Development Hypertrophic chondrocytes give rise to osteoblasts during skeletal development; however, the process by which these non-mitotic cells make this transition is not well understood. Prior studies have also suggested that skeletal stem and progenitor cells (SSPCs) localize to the surrounding periosteum and serve as a major source of marrow associated SSPCs, osteoblasts, osteocytes, and adipocytes during skeletal development. To further understand the cell transition process by which hypertrophic chondrocytes contribute to osteoblasts or other marrow associated cells, we utilized inducible and constitutive hypertrophic chondrocyte lineage tracing and reporter mouse models (Col10a1CreERT2; R26-tdTomatof/+ and Col10a1Cre; R26-tdTomatof/+) in combination with a PDGFRa-H2B-GFP transgenic line, single cell RNA-sequencing, bulk RNA-sequencing, immunofluorescence staining, and cell transplantation assays. Our data demonstrate that hypertrophic chondrocytes undergo a process of dedifferentiation to generate marrow associated SSPCs that serve as a primary source of osteoblasts during skeletal development. These hypertrophic chondrocyte derived SSPCs commit to a CXCL12-abundant reticular (CAR) cell phenotype during skeletal development and demonstrate unique abilities to recruit vasculature and promote bone marrow establishment, while also contributing to the adipogenic lineage. A generic hidden Markov model for multi-parent populations A common step in the analysis of multi-parent populations is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of multi-parent populations. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly-available data on Diversity Outbred and Collaborative Cross mice Genetic diversity loss in the Anthropocene More species than ever before are at risk of extinction due to anthropogenic habitat loss and climate change. But even species that are not threatened have seen reductions in their populations and geographic ranges, likely impacting their genetic diversity. Although preserving genetic diversity is key to maintaining adaptability of species, we lack predictive tools and global estimates of genetic diversity loss across ecosystems. By bridging theories of biodiversity and population genetics, we introduce a mathematical framework to understand the loss of naturally occurring DNA mutations within decreasing habitat within a species. Analysing genome-wide variation data of 10,095 geo-referenced individuals from 20 plant and animal species, we show that genome-wide diversity follows a power law with geographic area (the mutations-area relationship), which can predict genetic diversity loss in spatial computer simulations of local population extinctions. Given pre-21st century values of ecosystem transformations, we estimate that over 10% of genetic diversity may already be lost, surpassing the United Nations targets for genetic preservation. These estimated losses could rapidly accelerate with advancing climate change and habitat destruction, highlighting the need for forecasting tools that facilitate implementation of policies to protect genetic resources globally. Parallel processing in speech perception with local and global representations of linguistic context Speech processing is highly incremental. It is widely accepted that human listeners continuously use the linguistic context to anticipate upcoming concepts, words, and phonemes. However, previous evidence supports two seemingly contradictory models of how a predictive context is integrated with the bottom-up sensory input: Classic psycholinguistic paradigms suggest a two-stage process, in which acoustic input initially leads to local, context-independent representations, which are then quickly integrated with contextual constraints. This contrasts with the view that the brain constructs a single coherent, unified interpretation of the input, which fully integrates available information across representational hierarchies, and thus uses contextual constraints to modulate even the earliest sensory representations. To distinguish these hypotheses, we tested magnetoencephalography responses to continuous narrative speech for signatures of local and unified predictive models. Results provide evidence that listeners employ both types of models in parallel. Two local context models uniquely predict some part of early neural responses, one based on sublexical phoneme sequences, and one based on the phonemes in the current word alone; at the same time, even early responses to phonemes also reflect a unified model that incorporates sentence level constraints to predict upcoming phonemes. Neural source localization places the anatomical origins of the different predictive models in non-identical parts of the superior temporal lobes bilaterally, with the right hemisphere showing a relative preference for more local models. These results suggest that speech processing recruits both local and unified predictive models in parallel, reconciling previous disparate findings. Parallel models might make the perceptual system more robust, facilitate processing of unexpected inputs, and serve a function in language acquisition.

Impact statementSpeech processing engages multiple predictive models, using sublexical, word- and sentence contexts in parallel to anticipate upcoming phonemes. Differential regulation of alternative promoters emerges from unified kinetics of enhancer-promoter interaction Many eukaryotic genes contain alternative promoters with distinct expression patterns. How these promoters are differentially regulated remains elusive. Here, we apply single-molecule imaging to quantify the transcriptional regulation of two alternative promoters (P1 and P2) of the Bicoid (Bcd) target gene hunchback in syncytial blastoderm Drosophila embryos. Contrary to the previous notion that Bcd only activates P2, we find that Bcd activates both promoters via the same two enhancers. P1 activation is less frequent and requires binding of more Bcd molecules than P2 activation. Using a theoretical model to relate promoter activity to enhancer states, we show that the two promoters follow common transcription kinetics driven by sequential Bcd binding at the two enhancers. Bcd binding at either enhancer primarily activates P2, while P1 activation relies more on Bcd binding at both enhancers. These results provide a quantitative framework for understanding the kinetic mechanisms of complex eukaryotic gene regulation. Antibacterial potency of Type VI amidase effector toxins is dependent on substrate topology and cellular context Members of the bacterial T6SS amidase effector (Tae) superfamily of toxins are delivered between competing bacteria to degrade cell wall peptidoglycan. Although Taes share a common substrate, they exhibit distinct antimicrobial potency across different competitor species. To investigate the molecular basis governing these differences, we quantitatively defined the functional determinants of Tae1 from Pseudomonas aeruginosa PAO1 using a combination of nuclear magnetic resonance (NMR) and a high-throughput in vivo genetic approach called deep mutational scanning (DMS). As expected, combined analyses confirmed the role of critical residues near the Tae1 catalytic center. Unexpectedly, DMS revealed substantial contributions to enzymatic activity from a much larger, ring-like functional hot spot extending around the entire circumference of the enzyme. Comparative DMS across distinct growth conditions highlighted how functional contribution of different surfaces is highly context-dependent, varying alongside composition of targeted cell walls. These observations suggest that Tae1 engages with the intact cell wall network through a more distributed three-dimensional interaction interface than previously appreciated, providing an explanation for observed differences in antimicrobial potency across divergent Gram-negative competitors. Further binding studies of several Tae1 variants with their cognate immunity protein demonstrate that requirements to maintain protection from Tae activity may be a significant constraint on the mutational landscape of tae1 toxicity in the wild. In total, our work reveals that Tae diversification has likely been shaped by multiple independent pressures to maintain interactions with binding partners that vary across bacterial species and conditions. Pathogenic, but not non-pathogenic, Rickettsia evade inflammasome-dependent IL-1 responses to establish an intracytosolic replication niche. Rickettsia species (spp.) are strict obligate intracellular bacteria, with some being pathogenic in their mammalian host, including humans. One critical feature of these stealthy group of pathogens is their ability to manipulate hostile cytosolic environments to their benefits. Although our understanding of Rickettsia cell biology and pathogenesis are evolving, the mechanisms by which pathogenic Rickettsia spp. evade host innate immune detection remains elusive. Here, we showed that disease severity in wild-type (WT) C57BL/6J mice infected with R. typhi (etiologic agent of murine typhus) and R. rickettsii (etiologic agent of Rocky Mountain Spotted Fever), but not with non-pathogenic R. montanensis, correlated with levels of bacterial burden as detected in the spleens, as well as the serum concentrations of pro-inflammatory cytokine IL-1 and to a lesser extent IL- 1{beta}. Antibody-mediated neutralization of IL-1 confirmed a key role in controlling mortality rates and bacterial burdens of rickettsiae-infected WT mice. As macrophages are a primary source of both IL-1 and IL-1{beta} cytokines, we determined the mechanism of the anti-rickettsial activities using bone-marrow-derived macrophages. We found that pathogenic R. typhi and R. rickettsii, but not non-pathogenic R. montanensis, eluded pro- IL-1 induction and benefited pre-dominantly from the reduced IL-1 secretion, via a Caspase-11-Gsdmd-dependent pathway, to facilitate intracytosolic replication. Adoptative transfer experiments identified that IL-1 secretion by macrophages was critical for controlling rickettsiosis in WT mice. In sum, we identified a previously unappreciated pathway by which pathogenic, unlike non-pathogenic, rickettsiae preferentially target the Caspase-11-Gsdmd-IL-1 signaling axis in macrophages thus supporting their replication within the host.

IMPORTANCECurrently, no vaccines are available to prevent rickettsioses, while vector-borne rickettsial infections in humans are on the rise globally. In fact, the insufficient understanding of how pathogenic Rickettsia species circumvent host immune defense mechanisms has significantly hindered the development of more effective therapeutics. Here, we identified a previously unappreciated role for the Caspase-11-Gsdmd-IL-1 signaling axis, to limiting the replication of pathogenic R. rickettsia and R. typhi species in murine macrophages and wild-type (WT) C57BL/6J mice. Adoptative transfer studies further identified IL-1-secreting macrophages as critical mediators in controlling rickettsial infection in WT mice. Collectively, these findings provide insight into the potential mechanism of how pathogenic, but not non-pathogenic Rickettsia spp., benefit from a reduction in the Caspase-11-Gsdmd-mediated release of IL-1 to support host colonization. Deficits in seizure threshold and other behaviors in adult mice without gross neuroanatomic injury after late gestation transient prenatal hypoxia Intrauterine hypoxia is a common cause of brain injury in children resulting in a broad spectrum of long-term neurodevelopmental sequela, including life-long disabilities that can occur even in the absence of severe neuroanatomic damage. Postnatal hypoxia-ischemia rodent models are commonly used to understand the effects of ischemia and transient hypoxia on the developing brain. Postnatal models, however, have some limitations. First, they do not test the impact of placental pathologies on outcomes from hypoxia. Second, they primarily recapitulate severe injury because they provoke substantial cell death, which is not seen in children with mild hypoxic injury. Lastly, they do not model preterm hypoxic injury. Prenatal models of hypoxia in mice may allow us to address some of these limitations to expand our understanding of developmental brain injury. The published rodent models of prenatal hypoxia employ multiple days of hypoxic exposure or complicated surgical procedures, making these models challenging to perform consistently in mice. Furthermore, large animal models suggest that transient prenatal hypoxia without ischemia is sufficient to lead to significant functional impairment to the developing brain. However, these large animal studies are resource-intensive and not readily amenable to mechanistic molecular studies. Therefore, here we characterized the effect of late gestation (embryonic day 17.5) transient prenatal hypoxia (5% inspired oxygen) on long-term anatomical and neurodevelopmental outcomes in mice. Late gestation transient prenatal hypoxia increased hypoxia-inducible factor 1 alpha protein levels (a marker of hypoxic exposure) in the fetal brain. Hypoxia exposure predisposed animals to decreased weight at postnatal day 2, which normalized by day 8. However, hypoxia did not affect gestational age at birth, litter size at birth, or pup survival. No differences in fetal brain cell death or long-term gray or white matter changes resulted from hypoxia. Animals exposed to prenatal hypoxia did have several long-term functional consequences, including sex-dichotomous changes. Hypoxia exposure was associated with a decreased seizure threshold and abnormalities in hindlimb strength and repetitive behaviors in males and females. Males exposed to hypoxia had increased anxiety-related deficits, whereas females had deficits in social interaction. Neither sex developed any motor or visual learning deficits. This study demonstrates that late gestation transient prenatal hypoxia in mice is a simple, clinically relevant paradigm for studying putative environmental and genetic modulators of the long-term effects of hypoxia on the developing brain. Contingency and chance erase necessity in the experimental evolution of ancestral proteins To understand why evolution produced the biological systems that exist today, we must know how important chance, contingency, and necessity were during history. Previous observations suggest that each of these modes of causality affects evolution in various settings, but their relative roles and interactions are not well characterized because they have never been systematically assessed in a single system or on a timescale relevant to evolutionary history. To this end, we reconstructed ancestral B-cell-lymphoma-2-family proteins and developed a continuous evolution method to select for defined protein-protein interaction specificities. By repeatedly evolving a series of ancestral proteins to acquire specificities that occurred during history, we show that contingency steadily overwhelms chance and erases necessity as the primary cause of sequence variation in proteins over long phylogenetic timescales. As a result, evolutionary trajectories launched from distant starting points are essentially unpredictable, even under strong and identical selection pressures. Genetic dissection of the outcomes shows that chance arises because numerous sets of mutations can alter specificity at any point in time, while contingency arises because historical substitutions change these sets. Patterns of variation in extant protein sequences are therefore largely the idiosyncratic product of a particular course of unpredictable historical events. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. Spontaneous differentiation leads to emergence of hybrid cell states relate to poor prognosis in oral cancer PurposeCellular dynamics between phenotypically heterogeneous subpopulations of cancer cells within individual tumor is shown to be responsible for drug tolerance and overall poor prognosis; however, evidences were largely missing in oral cancer. Therefore, this study was undertaken to describe the dynamic phenotypic states among oral cancer cells, its influence on transcriptomic heterogeneity as well as its clinical significance.

Experimental DesignWe multiplexed phenotypic markers of putative oral-stem-like cancer cells (SLCCs) and characterized diversity among CD44-positive oral cancer cell subpopulations with respect to distinct expression of CD24 and aldehyde dehydrogenase (ALDH)-activity in multiple cell lines. Population trajectories were characterized by Markov model and cell states were defined based on the population specific RNA sequencing (RNAseq). ssGSEA based gene expression signatures were explore for prognostic significance.

ResultsOral cancer cells followed two distinct patterns of spontaneous repopulation dynamics with stochastic inter-conversions on  ALDH-axis, however a strict non-interconvertible transition on  CD24-axis. Interestingly, plastic  ALDH-axis was harnessed to enrich ALDHHigh subpopulations in response to Cisplatin treatment, to adapt a drug tolerant state. Phenotype-specific RNAseq results suggested the organization of subpopulations into hierarchical structure with possible maintenance of intermediate states of stemness within the differentiating oral cancer cells. Further, survival analysis with each subpopulation-specific gene signature strongly suggested that the cell-state dynamics may act as possible mechanism to drive ITH, resulting in poor prognosis in patient.

ConclusionsOur results emphasized the prognostic power of the population dynamics in oral cancer. Importantly, we have described the phenotypic-composition of heterogeneous subpopulations critical for global tumor behaviour in oral cancer; which is a prerequisite knowledge important for precision treatment, however largely lacking for most solid tumors.

Graphical AbstractWe have characterized diversity among CD44-positive oral cancer cells lines with respect to distinct expression of CD24 and ALDH-activity. Subpopulations showed stochastic inter-conversions on ALDH-axis but a strict non-interconvertible transition of CD24Low to CD24High phenotype, even in response to chemotherapy-induced stress. RNAseq study suggested the organization of subpopulations into hierarchical structure with possible maintenance of intermediate alternate states of stemness within the differentiating oral cancer cells. The described population dynamics demonstrtaed influence tumor behaviour possibly by increasing intratumoral heterogeneity in aggressive oral tumors.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=112 SRC="FIGDIR/small/457509v3_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@1dc6a7org.highwire.dtl.DTLVardef@dd0948org.highwire.dtl.DTLVardef@18c9603org.highwire.dtl.DTLVardef@cca8fb_HPS_FORMAT_FIGEXP  M_FIG C_FIG Translational relevanceIntratumoral heterogeneity (ITH) has been the clinically important factor, impacting aggressive cancer behaviour, drug tolerance and overall poor prognosis. Recent high-throughput studies have provided better cellular and molecular resolution of ITH; however, the prerequisite knowledge which defines the composition of subpopulations critical for global tumor behaviour is majorly lacking for most of the solid tumors. By combining phenotypic markers, we have defined four subpopulations of oral cancer cells. These subpopulations showed stochastic inter-conversions as well as a strict non-interconvertible transition among them to acheive heterogeneity. Importantly, transcriptional states of each subpopulations indicated a clinically relevant signatures for patient prognosis. Also, we observed interconversions of these subpopulations in response to Cisplatin to accumulate drug-tolerant cell state, as rapid and reversible strategy to respond to chemotherapy induced stress. Thus, the characteristics of described phenotypic subgroups may be translated to the clinic for estimating the extent of intratumoral heterogeneity in oral cancer patients. 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. Design, Mutate, Screen: High-throughput creation of genetic clocks with different period-amplitude characteristics Directed evolution has become an invaluable tool in protein engineering and has also greatly influenced the construction of synthetic gene circuits. The ability to generate diversity at precise targets for directed evolution approaches has improved vastly, allowing researchers to create large, specific mutant libraries with relative ease. Screening approaches for large mutant libraries have similarly come a long way, especially when the desired behavior can easily be tested for with static, single time-point assays. For more complex gene circuits with dynamic phenotypes that change over time, directed evolution approaches to controlling and tuning circuit behavior have been hindered by the lack of sufficiently high-throughput screening methods to isolate variants with desired characteristics. Here we utilize directed mutagenesis and multiplexed microfluidics to develop a workflow for creating, screening and tuning dynamic gene circuits that operate at the population level. Specifically, we create a mutant library of an existing oscillator, the synchronized lysis circuit, and tune its dynamics while uncovering principles regarding its behavior. Lastly, we utilize this directed evolution workflow to construct a new synchronized genetic oscillator that exhibits robust dynamics over long time scales. 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. Identification of differentially distributed gene expression and distinct sets of cancer-related genes identified by changes in mean and variability There is increasing evidence that changes in the variability or overall distribution of gene expression are important both in normal biology and in diseases, particularly cancer. Genes whose expression differs in variability or distribution without a difference in mean are ignored by traditional differential expression-based analyses. Using a Bayesian hierarchical model that provides tests for both differential variability and differential distribution for bulk RNA-seq data, we report here an investigation into differential variability and distribution in cancer. Analysis of eight paired tumour-normal datasets from The Cancer Genome Atlas confirms that differential variability and distribution are able to identify cancer-related genes. We further demonstrate that differential variability identifies cancer-related genes that are missed by differential expression analysis, and that differential expression and differential variability identify functionally distinct sets of genes. These results suggest that differential variability analysis may provide insights into genetic aspects of cancer that would not be revealed by differential expression, and that differential distribution analysis may allow for more comprehensive identification of cancer-related genes than analyses based on changes in mean or variability alone. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Fungicide resistance characterised across seven modes of action in a Botrytis cinerea isolated from Australian vineyards BACKGROUNDBotrytis bunch rot, caused by Botrytis cinerea, is an economically important disease of grapes in Australia and across grape growing regions worldwide. Control of this disease relies on canopy management and the application of fungicides. Fungicide application can lead to the selection of resistant B. cinerea populations, which has an adverse effect on management of the disease. Characterising the distribution and severity of resistant B. cinerea populations is needed to inform resistance management strategies.

RESULTSIn this study, 725 isolates were sampled from 75 Australian vineyards during 2013 - 2016 and were screened against seven fungicides with different modes of action (MOAs). The resistance frequencies for azoxystrobin, boscalid, fenhexamid, fludioxonil, iprodione, pyrimethanil and tebuconazole were 5, 2.8, 2.1, 6.2, 11.6, 7.7 and 2.9% respectively. Nearly half of the resistant isolates (43.8%) were resistant to more than one of the fungicides tested. The frequency of vineyards with at least one isolate simultaneously resistant to 1, 2, 3, 4 or 5 fungicides was 19.5, 7.8, 6.5, 10.4 and 2.6%. Resistance was associated with previously published genotypes in CytB (G143A), SdhB (H272R/Y), Erg27 (F412S), Mrr1 (D354Y), Bos1 (I365S, N373S + Q369P, I365S + D757N) and Pos5 (V273I, P319A, L412F/V). Novel genotypes were also described in Mrr1 (S611N, D616G) Pos5 (V273L) and Cyp51 (P347S). Expression analysis was used to characterise fludioxonil resistant isolates exhibiting overexpression (6.3-9.6-fold) of the ABC transporter gene AtrB (MDR1 phenotype).

CONCLUSIONResistance frequencies were lower when compared to most previously published surveys of B. cinerea resistance in grape and other crops. Nevertheless, continued monitoring of critical MOAs used in Australian vineyards is recommended. A cell type-specific multiomics uncovers a guard cell-specific RAF15-SnRK2.6/OST1 kinase cascade Multicellular organisms such as plants contain different cell types with specialized functions. Analyzing the characteristics of each cell type reveals specific cell functions and enhances understanding of organization and function at the organismal level. Here we report a highly-sensitive and efficient cell type-specific multiomics pipeline, combining simplified flow cytometry-based cell sorting for fluorescent protoplasts and optimized nanoscale proteomics and metabolomics methods, which allow in-depth analysis of the proteomes and metabolomes of a particular cell type. Using this method, we quantitatively compared the proteomes and metabolomes of guard cells and mesophyll cells and revealed that the enrichment of signal transduction-related proteins enables guard cells to respond rapidly to various environmental stimuli. We uncovered a guard cell-specific kinase cascade whereby RAF15 and OST1 mediate ABA-induced stomatal closure. This pipeline can be applied to various cell types in plant or non-plant systems to learn how cells function in highly organized multicellular organisms. Common coupling map advances GPCR-G protein selectivity Two-thirds of human hormones and one-third of clinical drugs act on membrane receptors that couple to G proteins to achieve appropriate functional responses. While G protein transducers from literature are annotated in the Guide to Pharmacology database, two recent large-scale datasets now expand the receptor-G protein  couplome. However, these three datasets differ in scope and reported G protein couplings giving different coverage and conclusions on GPCR-G protein signaling. Here, we report a common coupling map uncovering novel couplings supported by both large-scale studies, the selectivity/promiscuity of GPCRs and G proteins, and how the co-coupling and co-expression of G proteins compare to the families from phylogenetic relationships. The coupling map and insights on GPCR-G protein selectivity will catalyze advances in receptor research and cellular signaling towards the exploitation of G protein signaling bias in design of safer drugs. 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. Evidence for Deleterious Antigenic Imprinting in SARS-CoV-2 Immune Response A previous report demonstrated the strong association between the presence of antibodies binding to an epitope region from SARS-CoV-2 nucleocapsid, termed Ep9, and COVID-19 disease severity. Patients with anti-Ep9 antibodies (Abs) had hallmarks of antigenic imprinting (AIM), including early IgG upregulation and cytokine-associated injury. Thus, the immunological memory of a previous infection was hypothesized to drive formation of suboptimal anti-Ep9 Abs in severe COVID-19 infections. This study identifies a putative primary antigen capable of stimulating production of cross-reactive, anti-Ep9 Abs. Binding assays with patient blood samples directly show cross-reactivity between Abs binding to Ep9 and only one bioinformatics-derived, homologous potential antigen, a sequence derived from the neuraminidase protein of H3N2 Influenza A virus. This cross-reactive binding is highly influenza strain specific and sensitive to even single amino acid changes in epitope sequence. The neuraminidase protein is not present in the influenza vaccine, and the anti-Ep9 Abs likely resulted from the widespread influenza infection in 2014. Therefore, AIM from a previous infection could underlie some cases of COVID-19 disease severity.

ImportanceInfections with SARS-COV-2 result in diverse disease outcomes, ranging from asymptomatic to fatal. The mechanisms underlying different disease outcomes remain largely unexplained. Previously, our laboratory identified a strong association between the presence of an antibody and increased disease severity in a subset of COVID-19 patients. Here, we report that this severity-associated antibody cross-reacts with viral proteins from an influenza A viral strain from 2014. Therefore, we speculate that antibodies generated against previous infections, like the 2014 influenza A, play a significant role in directing some peoples immune responses against SARS-COV-2. Such understanding of the sources and drivers of COVID-19 disease severity can help early identification and pre-emptive treatment. B1 SINE-binding ZFP266 impedes reprogramming through suppression of chromatin opening mediated by pioneering factors Induced pluripotent stem cell reprogramming is inherently inefficient and understanding the molecular mechanisms underlying this inefficiency holds the key to successfully control cellular identity. Here, we report 16 novel reprogramming roadblock genes identified by CRISPR/Cas9-mediated genome-wide knockout (KO) screening. Of these, depletion of the predicted KRAB zinc finger protein (KRAB-ZFP) Zfp266 strongly and consistently enhanced iPSC generation in several iPSC reprogramming settings, emerging as the most robust roadblock. Further analyses revealed that ZFP266 binds Short Interspersed Nuclear Elements (SINEs) adjacent to binding sites of pioneering factors, OCT4 (POU5F1), SOX2 and KLF4, and impedes chromatin opening. Replacing the KRAB co-suppressor with a co-activator domain converted ZFP266 from a reprogramming inhibitor to a potent reprogramming facilitator. This work proposes SINE-KRAB-ZFP interaction to be a critical regulator of chromatin accessibility at enhancers for efficient cellular identity changes and also serves as a resource to further illuminate molecular mechanisms hindering reprogramming. Intermittent compression increases flow through a gel: a mechanical model perhaps relevant to flow of cerebrospinal fluid Exchange of molecules between cerebrospinal fluid (CSF) and brain cells contributes to brain function and protection from dementia. Despite widespread acceptance of the  glymphatic theory, the route by which CSF is brought close enough to the neural tissue for solutes to be exchanged by extracellular diffusion is not entirely clear. Exogenous molecules injected into CSF are observed to reach the basement lamina that surrounds the dense capillary network. Transport of solutes by diffusion along the basement lamina, a gel of macromolcules about 100 nm thick, would be too slow; bulk flow in a static geometry would require unphysiologically high pressures. However, it is known that the pulsation of blood aids transport of CSF, and we hypothesized that this is because the pulsation intermittently squeezes the pericapillary lamina. In a primitive mimicry, we have tested whether intermittent squeezing increases flow through an agar gel. In all but one of 216 tests, pulsation caused a reversible increase, sometimes by a factor of 100 or more. The enhancement was greatest for frequencies 5-11 Hz and, over the tested range of pressure heads (20 - 50 cmH2O), was greatest for the lowest pressure. The results suggest a reason why exercise slows the aging of the brain. Parameterized syncmer schemes improve long-read mapping MotivationSequencing long reads presents novel challenges to mapping. One such challenge is low sequence similarity between the reads and the reference, due to high sequencing error and mutation rates. This occurs, e.g., in a cancer tumor, or due to differences between strains of viruses or bacteria. A key idea in mapping algorithms is to sketch sequences with their minimizers. Recently, syncmers were introduced as an alternative sketching method that is more robust to mutations and sequencing errors.

ResultsWe introduce parameterized syncmer schemes, a generalization of syncmers, and provide a theoretical analysis for multi-parameter schemes. By combining these schemes with downsampling or minimizers we can achieve any desired compression and window guarantee. We implemented the use of parameterized syncmer schemes in the popular minimap2 and Winnowmap2 mappers. In tests on simulated and real long read data from a variety of genomes, the syncmer-based algorithms, with scheme parameters selected on the basis of the theoretical analysis, reduced unmapped reads by 20-60% at high compression while usually using less memory. The advantage was more pronounced at low sequence identity. At sequence identity of 75% and medium compression, syncmer-minimap had only 37% as many unmapped reads, and 8% fewer of the reads that did map were incorrectly mapped. Even at lower compression and error rates, parameterized syncmer based mapping mapped more reads than the original minimizer-based mappers as well as mappers using the original syncmer schemes. We conclude that using parameterized syncmer schemes can improve mapping of long reads in a wide range of settings.

Availabilityhttps://github.com/Shamir-Lab/syncmer_mapping

Supplementary informationSupplementary data are available at https://github.com/Shamir-Lab/syncmer_mapping.

Author summaryPopular long read mappers use minimizers, the minimal hashed k-mers from overlapping windows, as alignment seeds. Recent work showed that syncmers, which select a fixed set of k-mers as seeds, are more likely to be conserved under errors or mutations than minimizers, making them potentially useful for mapping error-prone long reads. We introduce a framework for creating syncmers, that we call parameterized syncmer schemes, which generalize those introduced so far, and provide a theoretical analysis of their properties. We implemented parameterized syncmer schemes in the minimap2 and Winnowmap2 long read mappers. Using parameters selected on the basis of our theoretical analysis we demonstrate improved mapping performance, with fewer unmapped and incorrectly mapped reads on a variety of simulated and real datasets. The improvements are consistent across a broad range of compression rates and sequence identities, with the most significant improvements for lower sequence identity (high error or mutation rates) and high compression.