The role of Kabuki Syndrome genes KMT2D and KDM6A in development: Analysis in Human sequencing data and compared to mice and zebrafish KMT2D and KDM6A are epigenetic regulators that have been implicated in Kabuki Syndrome, a rare congenital birth defect with multiple tissue and organ abnormalities, including craniofacial and heart defects. Our previous study identified human families with mutations in the epigenetic modifiers KMT2D and KDM6A, which is implicated in 32% and 10% of Kabuki Syndrome patients respectively. To understand the connection to Kabuki syndrome patients, and the transcriptional targets of KMT2D and KDM6A in humans, we performed RNA sequencing (seq) of lymphoblastoid cells from Kabuki Syndrome patients carrying mutations in KMT2D and KDM6A. We identified 1995 significant changes in transcriptional targets for KMT2D and 1917 for KDM6A, as compared to control. When compared with RNA-seq datasets obtained from other mouse and zebrafish studies, our analysis revealed KMT2D mutations affect the expression of 76 orthologous genes across all three datasets. Similarliy, KDM6A afftects the expprssion of 7 orthologous genes across three datasets. Despite the differences in cell types, stages, and species in the comparison between the transcriptomic datasets, there are common gene expression changes associated with KMT2D and KDM6A mutations. qPCR on novel zebrafish mutants confirmed the differentially expression in KMT2D or KDM6A mutant backgrounds. Taken together, our results show that KMT2D and KDM6A regulate common and unique genes across humans, mice, and zebrafish for early craniofacial and cardiac development and this information contributes to the understanding of epigenetic dysregulation during development of Kabuki syndrome.
Single-cell Transcriptomic Profiling Unveils Cardiac Cell-type Specific Response to Maternal Hyperglycemia Underlying the Risk of Congenital Heart Defects Congenital heart disease (CHD) is the most prevalent structural malformations of the heart affecting [~]1% of live births. To date, both damaging genetic variations and adverse environmental exposure such as maternal diabetes have been found to cause CHD. Clinical studies show [~]fivefold higher risk of CHD in the offspring of mothers with pregestational diabetes. Maternal pregestational diabetes affects the gene regulatory networks key to proper cardiac development in the fetus. However, the cell-type specificity of these gene regulatory responses to maternal diabetes and their association with the observed cardiac defects in the fetuses remains unknown. To uncover the transcriptional responses to maternal diabetes in the early embryonic heart, we used an established murine model of pregestational diabetes. In this model, we have previously demonstrated an increased incidence of CHD. Here, we show maternal hyperglycemia (matHG) elicits diverse cellular responses during heart development by single-cell RNA-sequencing in embryonic hearts exposed to control and matHG environment. Through differential gene-expression and pseudotime trajectory analyses of this data, we identified changes in lineage specifying transcription factors, predominantly affecting Isl1+ second heart field progenitors and Tnnt2+ cardiomyocytes with matHG. Using in vivo cell-lineage tracing studies, we confirmed that matHG exposure leads to impaired second heart field-derived cardiomyocyte differentiation. Finally, this work identifies matHG-mediated transcriptional determinants in cardiac cell lineages elevate CHD risk and show perturbations in Isl1-dependent gene-regulatory network (Isl1-GRN) affect cardiomyocyte differentiation. Functional analysis of this GRN in cardiac progenitor cells will provide further mechanistic insights into matHG-induced severity of CHD associated with diabetic pregnancies.
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. 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. Frameshift and wild-type proteins are highly similar because the genetic code and genomes were optimized for frameshift tolerance Frameshift protein sequences encoded by alternative reading frames of coding genes have been considered meaningless, and frameshift mutations have been considered of little importance for the molecular evolution of coding genes and proteins. However, functional frameshifts have been found widely existing. It was puzzling how a frameshift protein kept its structure and functionality while its amino-acid sequence was changed substantially. Here we show that frame similarities between frameshifts and wild types are higher than random similarities and are defined at the genetic code, gene, and genome levels. In the standard genetic code, frameshift codon substitutions are more conservative than random substitutions. The frameshift tolerability of the standard genetic code ranks in the top 2.0-3.5% of alternative genetic codes, showing that the genetic code is nearly optimal for frameshift tolerance. Furthermore, frameshift-resistant codons (codon pairs) appear more frequently than expected in many genes and certain genomes, showing that the frameshift optimality is reflected not only in the genetic code but more importantly, in its allowance of further optimizing the frameshift tolerance of a particular gene or genome, which shed light on the role of frameshift mutations in molecular and genomic evolution. LY-CoV1404 (bebtelovimab) potently neutralizes SARS-CoV-2 variants SARS-CoV-2 neutralizing monoclonal antibodies (mAbs) can reduce the risk of hospitalization when administered early during COVID-19 disease. However, the emergence of variants of concern has negatively impacted the therapeutic use of some authorized mAbs. Using a high throughput B-cell screening pipeline, we isolated a highly potent SARS-CoV-2 spike glycoprotein receptor binding domain (RBD)-specific antibody called LY-CoV1404 (also known as bebtelovimab). LY-CoV1404 potently neutralizes authentic SARS-CoV-2 virus, including the prototype, B.1.1.7, B.1.351 and B.1.617.2). In pseudovirus neutralization studies, LY-CoV1404 retains potent neutralizing activity against numerous variants including B.1.1.7, B.1.351, B.1.617.2, B.1.427/B.1.429, P.1, B.1.526, B.1.1.529, and the BA.2 subvariant and retains binding to spike proteins with a variety of underlying RBD mutations including K417N, L452R, E484K, and N501Y. Structural analysis reveals that the contact residues of the LY-CoV1404 epitope are highly conserved with the exception of N439 and N501. Notably, the binding and neutralizing activity of LY-CoV1404 is unaffected by the most common mutations at these positions (N439K and N501Y). The breadth of reactivity to amino acid substitutions present among current VOC together with broad and potent neutralizing activity and the relatively conserved epitope suggest that LY-CoV1404 has the potential to be an effective therapeutic agent to treat all known variants causing COVID-19.

In BriefLY-CoV1404 is a potent SARS-CoV-2-binding antibody that neutralizes all known variants of concern and whose epitope is rarely mutated.

HighlightsO_LILY-CoV1404 potently neutralizes SARS-CoV-2 authentic virus and known variants of concern including the B.1.1.529 (Omicron), the BA.2 Omicron subvariant, and B.1.617.2 (Delta) variants
C_LIO_LINo loss of potency against currently circulating variants
C_LIO_LIBinding epitope on RBD of SARS-CoV-2 is rarely mutated in GISAID database
C_LIO_LIBreadth of neutralizing activity and potency supports clinical development
C_LI Hidden dynamics of proteasome autoregulation discovered by cryo-EM data-driven deep learning The cellular functions are executed by macromolecular complexes in nonequilibrium dynamic processes, which exhibit a vast diversity of conformational states. Solving conformational continuum of important complexes at atomic level remains a key challenge in the methodology of computational biophysics. Here we introduce a deep-learning framework, named AlphaCryo4D, which enables atomic-level cryogenic electron microscopy reconstructions of nonequilibrium conformational continuum. AlphaCryo4D integrates 3D deep residual learning with manifold embedding of energy landscapes, which simultaneously improves 3D classification accuracy and reconstruction resolution via an energy-based particle-voting algorithm. Applications of this methodology to six experimental datasets demonstrated its generality in breaking resolution barrier of visualizing dynamic components of functional complexes, in choreographing continuous inter-subunit motions and in exploring their  hidden conformational space. Importantly, AlphaCryo4D identified 64 conformers of the substrate-bound human 26S proteasome and reconstituted single-nucleotide-exchange dynamics of proteasomal AAA-ATPase motor during translocation initiation, which together illuminates a grand hierarchical allostery for proteasome autoregulation. EPIC: inferring relevant cell types for complex traits by integrating genome-wide association studies and single-cell RNA sequencing More than a decade of genome-wide association studies (GWASs) have identified genetic risk variants that are significantly associated with complex traits. Emerging evidence suggests that the function of trait-associated variants likely acts in a tissue- or cell-type-specific fashion. Yet, it remains challenging to prioritize trait-relevant tissues or cell types to elucidate disease etiology. Here, we present EPIC (cEll tyPe enrIChment), a statistical framework that relates large-scale GWAS summary statistics to cell-type-specific gene expression measurements from single-cell RNA sequencing (scRNA-seq). We derive powerful gene-level test statistics for common and rare variants, separately and jointly, and adopt generalized least squares to prioritize trait-relevant cell types while accounting for the correlation structures both within and between genes. Using enrichment of loci associated with four lipid traits in the liver and enrichment of loci associated with three neurological disorders in the brain as ground truths, we show that EPIC outperforms existing methods. We apply our framework to multiple scRNA-seq datasets from different platforms and identify cell types underlying type 2 diabetes and schizophrenia. The enrichment is replicated using independent GWAS and scRNA-seq datasets and further validated using PubMed search and existing bulk case-control testing results. Critical re-evaluation of the slope factor of the operational model of agonism Background and PurposeAlthough being a relative term, agonist efficacy is a cornerstone in the proper assessment of agonist selectivity and signalling bias. The operational model of agonism (OMA) has become successful in the determination of agonist efficacies and ranking them. In 1985, Black et al. introduced the slope factor to the OMA to make it more flexible and allow for fitting steep as well as flat concentration-response curves.

Experimental ApproachFunctional analysis of OMA demonstrates that the slope factor implemented by Black et al. affects relationships among parameters of the OMA, leading to wrong estimates of operational efficacy and affinity of an agonist. In contrast, functional analysis of OMA with Hill coefficient as a slope factor has not this flaw.

Key ResultsFitting of the OMA with Black et al. slope factor to concentration-response curves of several models (homotropic allosteric modulation, substrate inhibition and non-competitive auto-inhibition) gives wrong estimates of operational efficacy and affinity. In contrast, fitting of the OMA modified by the Hill coefficient to concentration-response curves of these several models gives correct estimates of operational efficacy and affinity.

Conclusion and ImplicationsThe OMA modified by the Hill coefficient can be used for correct ranking of agonism and subsequent analysis, like quantification of signalling bias. Predictive coding across the left fronto-temporal hierarchy during language comprehension We used MEG and ERPs to track the time-course and localization of evoked activity produced by expected, unexpected plausible and implausible words during incremental language comprehension. We interpret the full pattern of evoked responses and functional connectivity within a hierarchical predictive coding framework in which evoked activity reflects residual information (error) that is passed up and down the cortical hierarchy. Between 300-500ms, the three conditions produced progressively larger responses within left temporal cortex (lexico-semantic prediction error), while implausible inputs produced a selectively enhanced response within left inferior frontal cortex (event-level prediction error). Between 600-1000ms, unexpected plausible words activated left inferior frontal and middle temporal cortices (feedback activity that induced top-down shifts of event-level and lexico-semantic representations), while highly implausible inputs activated posterior fusiform and medial temporal cortices, supporting orthographic reanalysis and new learning. Therefore, predictive coding may provide a unifying theory that links language comprehension to other domains of cognition. Short-term instantaneous prophylaxis and efficient treatment against SARS-CoV-2 in hACE2 mice conferred by an intranasal nanobody (Nb22) Current COVID-19 vaccines need to take at least one month to complete inoculation and then become effective. Around 51% global population are still not fully vaccinated. Instantaneous protection is an unmet need among those who are not fully vaccinated. In addition, breakthrough infections caused by SARS-CoV-2 are widely reported. All these highlight the unmet needing for short-term instantaneous prophylaxis (STIP) in the communities where SARS-CoV-2 is circulating. Previously, we reported nanobodies isolated from an alpaca immunized with the spike protein, exhibiting ultrahigh potency against SARS-CoV-2 and its variants. Herein, we found that Nb22, among our previously reported nanobodies, exhibited ultrapotent neutralization against Delta variant with an IC50 value of 0.41 ng/ml (5.13 pM). Furthermore, the crystal structural analysis revealed that the binding of Nb22 to WH01 and Delta RBDs both effectively blocked the binding of RBD to hACE2. Additionally, intranasal Nb22 exhibited protection against SARS-CoV-2 Delta variant in the post-exposure prophylaxis (PEP) and pre-exposure prophylaxis (PrEP). Of note, intranasal Nb22 also demonstrated high efficacy against SARS-CoV-2 Delta variant in STIP for seven days administered by single dose and exhibited long-lasting retention in the respiratory system for at least one month administered by four doses, providing a means of instantaneous short-term prophylaxis against SARS-CoV-2. Thus, ultrahigh potency, long-lasting retention in the respiratory system as well as stability at room-temperature make the intranasal or inhaled Nb22 to be a potential therapeutic or STIP agent against SARS-CoV-2.

Brief summaryNb22 exhibits ultrahigh potency against Delta variant in vitro and is exploited by crystal structural analysis; furthermore, animal study demonstrates high effectiveness in the treatment and short-term instantaneous prophylaxis in hACE2 mice via intranasal administration.

HighlightsO_LINb22 exhibits ultrapotent neutralization against Delta variant with an IC50 value of 0.41 ng/ml (5.13 pM).
C_LIO_LIStructural analysis elucidates the ultrapotent neutralization of Nb22 against Delta variant.
C_LIO_LINb22 demonstrates complete protection in the treatment of Delta variant infection in hACE2 transgenic mice.
C_LIO_LIWe complete the proof of concept of STIP against SARS-CoV-2 using intranasal Nb22 with ultrahigh potency and long-lasting retention in respiratory system.
C_LI

Graphic Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=199 SRC="FIGDIR/small/459055v2_ufig1.gif" ALT="Figure 1">
View larger version (44K):
org.highwire.dtl.DTLVardef@144516corg.highwire.dtl.DTLVardef@3dc17forg.highwire.dtl.DTLVardef@6a8962org.highwire.dtl.DTLVardef@619cd7_HPS_FORMAT_FIGEXP  M_FIG C_FIG Network centrality is the most important variable determining gene sequence evolution rates in Schizosaccharomyces pombe Which variables determine the rate of gene and protein sequence evolution is a central question in evolutionary genomics. In the model organism fission yeast (Schizosaccharomyces pombe), the determinants of the rate of sequence evolution have yet to be determined. Previous studies in other organisms have typically found gene expression levels to be most significant, with numerous other variables identified as having a smaller impact. Here, partial least squares regression (PLS) and partial correlation analysis are used to model sequence evolution rates in the fission yeast genome by a range of variables. Variable importance in projection (VIP) scores as well as partial correlation coefficients are calculated for each variable, and used as estimates of the influence of each independent variable on sequence evolution rate. Unlike many previous studies in other organisms, centrality in the PPI network is shown to be the most important variable, and gene expression found to be less influential. Considerable heterogeneity is found in the influence of different gene ontology terms as well as amino acid composition. However, the majority of variance in constraint in fission yeast remains unexplained by this study, indicating that variables not yet considered and stochastics probably have considerable impact on the rate of molecular evolution. The impact of housing conditions on mesenchymal stromal/stem cell populations differ between adipose tissue and skeletal muscle BackgroundIn pigs, the ratio between lean mass and fat mass in the carcass determines production efficiency and is strongly influenced by the number and size of cells in tissues. During growth, the increase in the number of cells results from the recruitment of different populations of multipotent mesenchymal stromal/stem cells (MSCs) residing in the tissues. We hypothesized that the impact of hygiene of housing conditions during growth on the proportions of MSCs in adipose tissue and skeletal muscle may differ between pigs with different residual feed intake (RFI), a measure of feed efficiency.

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

ResultsAdipose tissue and muscle shared some common MSC populations although MSC diversity was higher in muscle than in adipose tissue. In muscle, the CD45-CD56+CD34-CD140a+ and CD45-CD56+CD34+CD140a+ cell populations were abundant. Of these two cell populations, only the proportions of CD45-CD56+CD34+CD140a+ cells increased (P < 0.05) in pigs housed in poor hygiene as compared with pigs in good hygiene conditions. For the CD45-CD56-CD34- cell population, present in low proportion, there was an interaction between hygiene condition and genetic line (P < 0.05) with a decrease in low RFI pigs housed in poor hygiene conditions. In adipose tissue, the two abundant MSC populations were CD45-CD56-CD34- and CD45-CD56+CD34-. The proportion of CD45-CD56-CD34- cells increased (P < 0.05) whereas the proportion of CD45-CD56+CD34- tended to decrease (P < 0.1) in pigs housed in poor conditions. This study shows that the proportions of some MSC populations were affected by hygiene of housing conditions in a tissue-dependent manner in pigs of both RFI lines. It suggests that MSCs may play a significant role in adipose tissue and skeletal muscle homeostasis and may influence later growth and body composition in growing animals. Re-convolving the compositional landscape of primary and recurrent glioblastoma reveals prognostic andtargetable tissue states Glioblastoma is an aggressive diffusely infiltrating neoplasm that spreads beyond surgical resection margins, where it intermingles with non-neoplastic brain cells. This complex microenvironment harbouring infiltrating glioma and non-neoplastic brain cells is the origin of tumor recurrence. Thus, understanding the cellular and molecular features of the glioma microenvironment is therapeutically and prognostically important. Here, we used single-nucleus RNA sequencing (snRNAseq) of primary and recurrent glioma to define three compositional  tissue-states rooted in cohabitation of cell-types and transcriptional states. These comprise a state featuring A) abundant normal brain microenvironment cells, B) reactive/inflamed brain microenvironment, and C) cellular/proliferative tumor. All these states exhibited variable degrees of infiltration by glioma cells. Spatial transcriptomics confirmed that the cell-types and transcriptomics states which compositionally cohabitate tissue states also colocalize in space. Tissue states are clinically significant because they correlate with radiographic, histopathologic, and prognostic features. We found that enrichment of tissue state B signature correlated with shorter survival. Importantly, we found that our compositionally defined tissue states are enriched in distinct metabolic pathways. One such pathway is fatty acid biosynthesis, which was enriched in tissue state B - a state enriched in recurrent glioblastoma and composed of astrocyte-like/mesenchymal glioma cells, reactive astrocytes resembling those seen in neurodegeneration, and monocyte-like myeloid cells. We showed that treating acute slices of GBM with a fatty acid synthesis inhibitor is sufficient to deplete the transcriptional signature of tissue state B. Our findings define a novel compositional approach to the glioma-infiltrated tissue which allows us to discover prognostic and targetable features, paving the way to new mechanistic and therapeutic discoveries. Interfacial tension driven open droplet microfluidics We developed an open channel droplet microfluidic system that autonomously generates droplets at low Ca (~10-4-10-3) by leveraging competing hydrostatic and capillary pressure. With only our open channel polytetrafluoroethylene (PTFE) device, pipettes, and commercially available carrier fluid, we produce hundreds of microliter droplets; tubing, electronics, or pumps are not required, making droplet technology feasible for research labs without external flow generators. Furthermore, we demonstrated conceptual applications that showcase the process of droplet generation, splitting, transport, incubation, mixing, and sorting in our system. Unlike conventional droplet microfluidics, the open nature of the device enables the use of physical tools such as tweezers and styli to directly access the system; with this, we developed a new method of droplet sorting and transfer that capitalizes on the Cheerios effect, the aggregation of buoyant objects along a liquid interface. Our platform offers enhanced usability, direct access to the droplet contents, easy manufacturability, compact footprint, and high customizability. This design is a first step in exploring the space of power-free open droplet microfluidic systems and provide design rules for similar channel designs. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. Heterologous autoimmunity and prokaryotic immune defense Some prokaryotes possess CRISPR-Cas systems that provide adaptive immunity to viruses guided by DNA segments called spacers acquired from invading phage. However, the patchy incidence and limited memory breadth of CRISPR-Cas systems suggest that their fitness benefits are offset by costs. Here, we propose that cross-reactive CRISPR targeting can lead to heterologous autoimmunity, whereby foreign spacers guide self-targeting in a spacer-length dependent fashion. Balancing antiviral defense against autoimmunity predicts a scaling relation between spacer length and CRISPR repertoire size. We find evidence for this scaling through comparative analysis of sequenced prokaryotic genomes, and show that this association also holds at the level of CRISPR types. In contrast, the scaling is absent in strains with nonfunctional CRISPR loci. Finally, we demonstrate that stochastic spacer loss can explain variations around the scaling relation, even between strains of the same species. Our results suggest that heterologous autoimmunity is a selective factor shaping the evolution of CRISPR-Cas systems. PRSS2 stimulates tumor growth by remodeling the TME via repression of Tsp1 In the earliest stages of tumor development, epithelial tumors (carcinomas) are physically confined to the area of the tissue in which they form. These nascent lesions (carcinomas in situ) are sequestered from the tissue parenchyma by the basement membrane. Within the tissue parenchyma lie a myriad of cell types comprised of fibroblasts, immune and inflammatory cells and endothelial cells. Upon invasion across the basement membrane and into the tissue parenchyma, tumors must manipulate the expression of pro- and anti-tumorigenic proteins such that pro-tumorigenic factors are produced in excess to anti-tumorigenic proteins. One such anti-tumorigenic protein is Thrombospondin-1 (Tsp-1). We have previously demonstrated that stimulation of Tsp-1 in the tumor microenvironment (TME) potently inhibits tumor growth and progression and in some cases induces tumor regression. Here, we identify a novel tumor-mediated mechanism to repress the expression of Tsp-1 in the TME via secretion of the serine protease PRSS2. We demonstrate that PRSS2 represses Tsp-1, not via its enzymatic activity, but by binding to low-density lipoprotein receptor-related protein 1 (LRP1). These findings describe a novel activity for PRSS2 through binding to LRP1 and represent a potential therapeutic strategy to treat cancer by blocking the PRSS2-mediated repression of Tsp-1. Quantifying the Morphology and Mechanisms of Cancer Progression in 3D in-vitro environments: Integrating Experiments and Multiscale Models Mathematical models of cancer growth have become increasingly more accurate both in the space and time domains. However, the limited amount of data typically available has resulted in a larger number of qualitative rather than quantitative studies. In the present study, we provide an integrated experimental-computational framework for the quantification of the morphological characteristics and the mechanistic modelling of cancer progression in 3D environments. The proposed framework allows for the calibration of multiscale, spatiotemporal models of cancer growth using state-of-the-art 3D cell culture data, and their validation based on the resulting experimental morphological patterns. Its implementation enabled us to pursue two goals; first, the quantitative description of the morphology of cancer progression in 3D cultures, and second, the relation of tumour morphology with underlying biophysical mechanisms that govern cancer growth and migration. We applied this framework to the study of the spatiotemporal progression of Triple Negative Breast Cancer cells cultured in Matrigel scaffolds, and validated the hypothesis of chemotactic migration using a multiscale Keller-Segel model. The results revealed transient, non-random spatial distributions of cancer cells that consist of clustered, and dispersion patterns. The proposed model was able to describe the general characteristics of the experimental observations and suggests that cancer cells exhibited chemotactic migration and accumulation, as well as random motion during the examined time period of development. To our knowledge, this is the first time that a multiscale model is used to quantify the relationship between the spatial patterns and the underlying mechanisms of cancer growth in 3D environments. A Wnt/β-catenin signaling driven adipocyte population required for beiging and glucose homeostasis Wnt/{beta}-catenin signaling has been well established as a potent inhibitor of adipogenesis. Here, we identified a population of adipocytes that exhibit persistent activity of Wnt/{beta}-catenin signaling, as revealed by the Tcf/Lef-Gfp reporter allele, in embryonic and adult mouse fat depots, named as Wnt+ adipocytes. We showed that the {beta}-catenin-mediated signaling activation in these cells is Wnt ligand- and receptor-independent but relies on AKT/mTOR pathway and is essential for cell survival. Such adipocytes are distinct from classical ones in transcriptomic and genomic signatures and can be induced from various sources of mesenchymal stromal cells including human cells. Genetic lineage-tracing and targeted cell ablation studies revealed that these adipocytes convert into beige adipocytes directly and are also required for beige fat recruitment under thermal challenge, demonstrating both cell autonomous and non-cell autonomous roles in adaptive thermogenesis. Furthermore, mice bearing targeted ablation of these adipocytes exhibited glucose intolerance, while mice receiving exogenously supplied such cells manifested enhanced glucose utilization. Our studies uncover a unique adipocyte population in regulating beiging in adipose tissues and systemic glucose homeostasis. Chromatin-contact atlas reveals disorder-mediated protein interactions and moonlighting chromatin-associated RBPs RNA-binding proteins (RBPs) play diverse roles in regulating co-transcriptional RNA-processing and chromatin functions, but our knowledge of the repertoire of chromatin-associated RBPs (caRBPs) and their interactions with chromatin remains limited. Here, we developed SPACE (Silica Particle Assisted Chromatin Enrichment) to isolate global and regional chromatin components with high specificity and sensitivity, and SPACEmap to identify the chromatin-contact regions in proteins. Applied to mouse embryonic stem cells, SPACE identified 1,459 chromatin-associated proteins, [~]48% of which are annotated as RBPs, indicating their dual roles in chromatin and RNA-binding. Additionally, SPACEmap stringently verified chromatin-binding of 404 RBPs and identified their chromatin-contact regions. Notably, SPACEmap showed that about half of the caRBPs bind to chromatin by intrinsically disordered regions (IDRs). Studying SPACE and total proteome dynamics from mES cells grown in 2iL and serum medium indicates significant correlation (R = 0.62). One of the most dynamic caRBPs is Dazl, which we find co-localized with PRC2 at transcription start sites of genes that are distinct from Dazl mRNA binding. Dazl and other PRC2-colocalised caRBPs are rich in intrinsically disordered regions (IDRs), which could contribute to the formation and regulation of phase-separated PRC condensates. Together, our approach provides an unprecedented insight into IDR-mediated interactions and caRBPs with moonlighting functions in native chromatin. Artificial seawater based long-term culture of colonial ascidians Tunicates are highly diverse marine invertebrate filter-feeders that are vertebrates closest relatives. These organisms, despite a drastically different body plan during their adulthood, have a tissue complexity related to that of vertebrates. Ascidians, which compose most of the Tunicata, are benthic sessile hermaphrodites that reproduce sexually through a motile tadpole larval stage. Over half of the known ascidians species are able to reproduce asexually by budding, typically leading to the formation of colonies where animals, called zooids, are interconnected through an external vascular system. In addition, colonial ascidians are established models for important biological processes including allorecognition, immunobiology, aging, angiogenesis and whole-body regeneration. However, the current paucity in breeding infrastructures limits the study of these animals to coastal regions.

To promote a wider scientific spreading and popularity of colonial ascidians, we have developed a flexible recirculating husbandry setup for their long-term in-lab culture. Our system is inspired both by the flow-through aquariums used by coastal ascidian labs, as well as by the recirculating in-lab systems used for zebrafish research. Our hybrid system thus combines colony breeding, water filtering and food culturing in a semi-automated system where specimens develop on hanging microscopy glass slides. Temperature, light/dark cycles, flow speed and feeding rates can be controlled independently in four different breeding environments to provide room for species-specific optimization as well as for running experiments. This setup is complemented with a quarantine for the acclimatization of wild isolates.

Herein we present our success in breeding Botrylloides diegensis, a species of colonial ascidians, for more than 3 years in recirculating artificial seawater over 600 km away from their natural habitat. We show that colonies adapt well to in-lab culturing provided that a suitable marine microbiome is present, and that a specific strain can be isolated, propagated and efficiently used for research over prolonged periods of time. The flexible and modular structure of our system can be scaled and adapted to the needs of specific species, such as Botryllus schlosseri, as well as of particular laboratory spaces. Overall, we show that Botrylloides diegensis can be proficiently bred in-land and suggest that our results can be extended to other species of colonial ascidians to promote research on these fascinating animals.

HighlightsO_LIFirst artificial seawater based recirculating aquaculture for colonial ascidians
C_LIO_LIOver 3 years of continuous breeding
C_LIO_LISemi-automated setup with minimized maintenance
C_LIO_LIGood biomass production for strain propagation
C_LIO_LI4 different culture conditions for optimized breeding for species of interest
C_LI DOMINANT AWN INHIBITOR encodes the ALOG protein originating from gene duplication and inhibits awn elongation by suppressing cell proliferation and elongation in sorghum The awn, a needle-like structure extending from the tip of the lemma in grass species, plays a role in environmental adaptation and fitness. In some crops, awns appear to have been eliminated during domestication. Although numerous genes involved in awn development have been identified, several dominant genes that eliminate awns are also known to exist. For example, in sorghum (Sorghum bicolor), the dominant awn inhibiting gene has been known since 1921; however, its molecular features remain uncharacterized. In this study, we conducted quantitative trait locus analysis and a genome-wide association study of awn-related traits in sorghum and identified DOMINANT AWN INHIBITOR (DAI), which encodes the ALOG family protein on chromosome 3. DAI appeared to be present in most awnless sorghum cultivars, likely because of its effectiveness. Detailed analysis of the ALOG protein family in cereals revealed that DAI originated from duplication of its twin paralog (DAIori) on chromosome 10. Observations of immature awns in near-isogenic lines revealed that DAI inhibits awn elongation by suppressing both cell proliferation and elongation. We also found that only DAI gained a novel function to inhibit awn elongation through an awn-specific expression pattern distinct from that of DAIori. Interestingly, heterologous expression of DAI with its own promoter in rice inhibited awn elongation in the awned cultivar Kasalath. We found that DAI originated from gene duplication, providing an interesting example of gain-of-function that occurs only in sorghum but shares its functionality with rice and sorghum. Including imprecisely georeferenced specimens improves accuracy of species distribution models and estimates of niche breadth: Don't let the perfect be the enemy of the good Museum and herbarium specimen records are frequently used to assess species conservation status and responses to climate change. Typically, records of occurrence with imprecise geolocality information are discarded because they cannot be matched confidently to environmental conditions, and are thus expected to increase uncertainty in downstream analyses. However, using only those records that can be precisely georeferenced risks undersampling of species environmental and geographic distributions. Using simulated and real species, we compared the effect of discarding versus retaining imprecise records on the accuracy of ecological niche models (ENMs), and estimates of niche breadth and extent of occurrence. Imprecise records were assigned to locations or climate values using a conservative approach. For simulated species, including imprecisely-georeferenced records alongside precisely-georeferenced records improved the accuracy of ENMs projected to the present and the future. Improvements were especially notable for species with <~20 precise occurrences, but accuracy continued to improve and eventually matched the accuracy of ENMs using all records with no georeferencing error. Using only precise records underestimated loss in suitable habitat and overestimated the amount of suitable habitat in both the present and future. Including imprecise records also improved univariate and multivariate estimates of niche breadth, and of species extent of occurrence. An analysis of 44 species in the genus Asclepias (milkweeds) revealed similarly large differences between cases using only precise records and those using precise plus imprecise records. Although the inclusion of imprecisely-georeferenced occurrence records may not always be appropriate, careful consideration must be made regarding the trade-offs between fidelity of the match between occurrences and the environment, versus sample size and how well precisely-geolocated records sample species environmental and geographic distributions. The preponderance of imprecisely georeferences specimens in museums and herbaria could be gainfully employed to address known shortfalls in sampling of species distributions and climatic niche tolerances. Population bottlenecks constrain host microbiome diversity and genetic variation impeding fitness It is becoming increasingly clear that microbial symbionts influence key aspects of their hosts fitness, and vice versa. This may fundamentally change our thinking about how microbes and hosts interact in influencing fitness and adaptation to changing environments. Here we explore how reductions in population size commonly experienced by threatened species influence microbiome diversity. Consequences of such reductions are normally interpreted in terms of a loss of genetic variation, increased inbreeding and associated inbreeding depression. However, fitness effects of population bottlenecks might also be mediated through microbiome diversity, such as through loss of functionally important microbes. Here we utilise 50 Drosophila melanogaster lines with different histories of population bottlenecks to explore these questions. The lines were phenotyped for egg-to-adult viability and their genomes sequenced to estimate genetic variation. The bacterial 16S rRNA gene was amplified in these lines to investigate microbial diversity. We found that 1) host population bottlenecks constrained microbiome richness and diversity, 2) core microbiomes of hosts with low genetic variation were constituted from subsets of microbiomes found in flies with higher genetic variation, 3) both microbiome diversity and host genetic variation contributed to host population fitness, 4) connectivity and robustness of bacterial networks was low in the inbred lines regardless of host genetic variation, 5) reduced microbial diversity was associated with weaker evolutionary responses of hosts in stressful environments, and 6) these effects were unrelated to Wolbachia density. These findings suggest that population bottlenecks reduce hologenomic variation (combined host and microbial genetic variation). Thus, while the current biodiversity crisis focuses on population sizes and genetic variation of eukaryotes, an additional focal point should be the microbial diversity carried by the eukaryotes, which in turn may influence host fitness and adaptability with consequences for the persistence of populations.

Author summaryIt is becoming increasingly clear that organisms and the microbes that live on or in them - their microbiome - affect each other in profound ways that we are just beginning to understand. For instance, a diverse microbiome can help maintain metabolic functions or fight pathogens causing diseases. A disrupted microbiome may be especially critical for animals and plants that occur in low numbers because of threats from e.g. human exploitation or climate change, as they may already suffer from genetic challenges such as inbreeding and reduced evolutionary potential. The importance of such a reduction in population size, called a bottleneck, on the microbial diversity and the potential interactive effects on host health remains unexplored. Here we experimentally test these associations by investigating the microbiomes of 50 inbred or non-inbred populations of vinegar flies. We found that restricting the population size constrain the hosts genetic variation and simultaneously decreases the diversity of the microbiome that they harbor, and that both effects were detrimental to host fitness. The microbial communities in inbred host populations were less robust than in their non-inbred counterparts, suggesting that we should increasingly consider the microbiome diversity, which may ultimately influence the health and persistence of threatened species. Different contribution of the monkey prefrontal and premotor dorsal cortex in decision making during a transitive inference task Several studies have reported similar neural modulations between brain areas of the frontal cortex, such as the dorsolateral prefrontal (DLPFC) and the premotor dorsal (PMd) cortex, in tasks requiring encoding of the abstract rules for selecting the proper action. Here we compared the neuronal modulation of the DLPFC and PMd of monkeys trained to choose the higher rank from a pair of abstract images (target item), selected from an arbitrarily rank-ordered set (A>B>C>D>E>F) in the context of a transitive inference task. Once acquired by trial-and-error, the ordinal relationship between pairs of adjacent images (i.e., A>B; B>C; C>D; D>E; E>F), monkeys were tested in indicating the ordinal relation between items of the list not paired during learning. During these decisions, we observed that the choice accuracy increased and the reaction time decreased as the rank difference between the compared items enhanced. This result is in line with the hypothesis that after learning, the monkeys built an abstract mental representation of the ranked items, where rank comparisons correspond to the items position comparison on this representation. In both brain areas, we observed higher neuronal activity when the target item appeared in a specific location on the screen with respect to the opposite position and that this difference was particularly enhanced at lower degrees of difficulty. By comparing the time evolution of the activity of the two areas, we observed that the neural encoding of target item spatial position occurred earlier in the DLPFC than in the PMd.