Numerosity tuning in human association cortices and local image contrast representations in early visual cortex Human early visual cortex response amplitudes monotonically increase with numerosity (object number), regardless of object size and spacing. However, numerosity is typically considered a high-level visual or cognitive feature, while early visual responses follow image contrast in the spatial frequency domain. We found that, at fixed contrast, aggregate Fourier power (at all orientations and spatial frequencies) followed numerosity closely but nonlinearly with little effect of object size, spacing or shape. This would allow straightforward numerosity estimation from spatial frequency domain image representations. Using 7T fMRI, we showed monotonic responses originate in primary visual cortex (V1) at the stimuluss retinotopic location. Responses here and in neural network models followed aggregate Fourier power more closely than numerosity. Truly numerosity tuned responses emerged after lateral occipital cortex and were independent of retinotopic location. We propose numerositys straightforward perception and neural responses may have built on behaviorally beneficial spatial frequency analyses in simpler animals.
Whole-Brain opto-fMRI Map of Mouse VTA Dopaminergic Activation Reflects Structural Projections with Small but Significant Deviations Ascending dopaminergic projections from neurons located in the Ventral Tegmental Area (VTA) are key to the etiology, dysfunction, and control of motivation, learning, and addiction. Due to evolutionary conservation of this nucleus and the extensive use of mice as disease models, establishing an assay for VTA dopaminergic signalling in the mouse brain is crucial for the translational investigation of motivational control as well as of neuronal function phenotypes for diseases and interventions. In this article we use optogenetic stimulation directed at VTA dopaminergic neurons in combination with functional Magnetic Resonance Imaging (fMRI), a method widely used in human deep brain imaging. We present a comprehensive assay producing the first whole-brain opto-fMRI map of dopaminergic activation in the mouse, and show that VTA dopaminergic system function is consistent with its structural VTA projections, diverging only in a few key aspects. While the activation map predominantly highlights target areas according to their relative projection densities (e.g. strong activation of the nucleus accumbens and low activation of the hippocampus), it also includes areas for which a structural connection is not well established (such as the dorsomedial striatum). We further detail the variability of the assay with regard to multiple experimental parameters, including stimulation protocol and implant position, and provide evidence-based recommendations for assay reuse, publishing both reference results and a reference analysis workflow implementation.
Ubiquity and origins of structural maintenance of chromosomes (SMC) proteins in eukaryotes Structural maintenance of chromosomes (SMC) protein complexes are common in Bacteria, Archaea, and Eukaryota. SMC proteins, together with the proteins related to SMC (SMC-related proteins), constitute a superfamily of ATPases. Bacteria/Archaea and Eukaryotes are distinctive from one another in terms of the repertory of SMC proteins. A single type of SMC protein is dimerized in the bacterial and archaeal complexes, whereas eukaryotes possess six distinct SMC subfamilies (SMC1-6), constituting three heterodimeric complexes, namely cohesin, condensin, and SMC5/6 complex. Thus, to bridge the homodimeric SMC complexes in Bacteria and Archaea to the heterodimeric SMC complexes in Eukaryota, we need to invoke multiple duplications of an SMC gene followed by functional divergence. However, to our knowledge, the evolution of the SMC proteins in Eukaryota had not been examined for more than a decade. In this study, we reexamined the ubiquity of SMC1-6 in phylogenetically diverse eukaryotes that cover the major eukaryotic taxonomic groups recognized to date and provide two novel insights into the SMC evolution in eukaryotes. First, multiple secondary losses of SMC5 and SMC6 occurred in the eukaryotic evolution. Second, the SMC proteins constituting cohesin and condensin (i.e., SMC1-4), and SMC5 and SMC6 were derived from closely related but distinct ancestral proteins. Based on the above-mentioned findings, we discuss how SMC1-6 have diverged from the archaeal homologs. The OPAQUE1/DISCORDIA2 myosin XI is required for phragmoplast guidance during asymmetric cell division in maize Asymmetric divisions produce cells with different fates and are critical for development. Here we show that a maize myosin XI protein, OPAQUE1 (O1), is necessary for asymmetric divisions during maize stomatal development. We analyzed stomatal precursor cells prior to and during asymmetric division to determine why o1 mutants have abnormal division planes. Cell polarization and nuclear positioning occur normally in the o1 mutant, and the future site of division is correctly specified. The defect in o1 occurs during late cytokinesis, when the plant-specific phragmoplast -- made of microtubules, actin and other proteins - forms the nascent cell plate. The phragmoplast becomes misguided and does not meet the previously established division site. Initial phragmoplast guidance is correct in o1. However, as phragmoplast expansion continues, phragmoplasts in o1 stomatal precursor cells become misguided and do not meet the cortex at the established division site. To understand how this myosin protein contributes to phragmoplast guidance, we identified O1-interacting proteins. Other myosins, specific actin-binding proteins, and maize kinesins related to the Arabidopsis thaliana division site markers PHRAGMOPLAST ORIENTING KINESINs (POKs) interact with O1. We propose that different myosins are important at multiple steps of phragmoplast expansion, and the O1 actin motor and POK-like microtubule motors work together to ensure correct late-stage phragmoplast guidance. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. Induction of exaggerated cytokine production in human peripheral blood mononuclear cells by a recombinant SARS-CoV-2 spike glycoprotein S1 and its inhibition by dexamethasone An understanding of the pathological inflammatory mechanisms involved in SARS-CoV-2 virus infection is necessary in order to discover new molecular pharmacological targets for SARS-CoV-2 cytokine storm. In this study, the effects of a recombinant SARS-CoV-2 spike glycoprotein S1 was investigated in human peripheral blood mononuclear cells (PBMCs). Stimulation of PBMCs with spike glycoprotein S1 (100 ng/mL) resulted in significant elevation in the production of TNF, IL-6, IL-1{beta} and IL-8. However, pre-treatment with dexamethasone (100 nM) caused significant reduction in the release of these cytokines. Further experiments revealed that S1 stimulation of PBMCs increased phosphorylation of NF-{kappa}B p65 and I{kappa}B, and I{kappa}B degradation. DNA binding of NF-{kappa}B p65 was also significantly increased following stimulation with spike glycoprotein S1. Treatment of PBMCs with dexamethasone (100 nM) or BAY11-7082 (1 {micro}M) resulted in inhibition of spike glycoprotein S1-induced NF-{kappa}B activation. Activation of p38 MAPK by S1 was blocked in the presence of dexamethasone and SKF 86002. CRID3, but not dexamethasone pre-treatment produced significant inhibition of S1-induced activation of NLRP3/caspase-1. Further experiments revealed that S1-induced increase in the production of TNF, IL-6, IL-1{beta} and IL-8 was reduced in the presence of BAY11-7082 and SKF 86002, while CRID3 pre-treatment resulted in the reduction of IL-1{beta} production. These results suggest that SARS-CoV-2 spike glycoprotein S1 stimulated PBMCs to release pro-inflammatory cytokines through mechanisms involving activation of NF-{kappa}B, p38 MAPK and NLRP3 inflammasome. It is proposed that the clinical benefits of dexamethasone in COVID-19 is possibly due to its anti-inflammatory activity in reducing SARS-CoV-2 cytokine storm. Regulation of the PKD2 channel function by TACAN Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in membrane receptor PKD1 or cation channel PKD2. TACAN (also named TMEM120A), recently reported as an ion channel in neuron cells for mechano and pain sensing, is also distributed in diverse non-neuronal tissues such as kidney, heart and intestine, suggesting its involvement in other functions. In this study, we found that TACAN is in complex with PKD2 in native renal cell lines. Using the two-electrode voltage clamp in Xenopus oocytes we found that TACAN inhibited the channel activity of PKD2 gain-of-function mutant F604P. The first and last transmembrane domains of TACAN were found to interact with the PKD2 C-and N-terminal portions, respectively. We showed that the TACAN N-terminus acted as a blocking peptide and that TACAN inhibits the PKD2 function through the PKD2/TACAN binding. By patch clamping in mammalian cells, we found that TACAN inhibits both the single channel conductance and open probability of PKD2 and mutant F604P. PKD2 co-expressed with TACAN, but not PKD2 alone, exhibited pressure sensitivity. Furthermore, we also found that TACAN aggravates PKD2-dependent tail curvature and pronephric cysts in larval zebrafish, in support of the in vitro inhibitory effects of TACAN. In summary, this study revealed that TACAN acts as a PKD2 inhibitor and mediates mechano sensitivity of the PKD2/TACAN channel complex. 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. Building a Knowledge Graph to Enable Precision Medicine Developing personalized diagnostic strategies and targeted treatments requires a deep understanding of disease biology and the ability to dissect the relationship between molecular and genetic factors and their phenotypic consequences. However, such knowledge is fragmented across publications, non-standardized research repositories, and evolving ontologies describing various scales of biological organization between genotypes and clinical phenotypes. Here, we present PrimeKG, a precision medicine-oriented knowledge graph that provides a holistic view of diseases. PrimeKG integrates 20 high-quality resources to describe 17,080 diseases with 4,050,249 relationships representing ten major biological scales, including disease-associated protein perturbations, biological processes and pathways, anatomical and phenotypic scales, and the entire range of approved and experimental drugs with their therapeutic action, considerably expanding previous efforts in disease-rooted knowledge graphs. In addition, PrimeKG supports artificial intelligence analyses of how drugs might target disease-associated molecular perturbations by containing an abundance of 'indications', 'contradictions', and 'off-label use' drug-disease edges lacking in other knowledge graphs. We accompany PrimeKG's graph structure with text descriptions of clinical guidelines to enable multimodal analyses. An epigenomic shift in amygdala marks the transition to maternal behaviors in alloparenting virgin female mice In many species, adults care for young offspring that are not their own, a phenomenon called alloparenting. However, most nonparental adults must be sensitized by repeated or extended exposures to newborns before they will robustly display parental-like behaviors. To capture neurogenomic events underlying the transition to active parental caring behaviors, we analyzed brain gene expression and chromatin profiles of virgin female mice co-housed with mothers during pregnancy and after birth. After an initial display of antagonistic behaviors and a surge of defense-related gene expression, we observed a dramatic shift in the chromatin landscape specifically in amygdala of the pup-exposed virgin females, accompanied by a dampening of anxiety-related gene expression. This epigenetic shift coincided with hypothalamic expression of the oxytocin gene and the emergence of behaviors and gene expression patterns classically associated with maternal care. The results outline a neurogenomic program associated with dramatic behavioral changes and suggest molecular networks relevant to human postpartum mental health. 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. A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. 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. Chromosome breaks in breast cancers occur near herpes tumor virus sequences and explain why cancer drugs stop working One explanation for why breast cancer has a relentless tendency to come back after treatment is that breast cancer phenomes are fluid, so it is difficult for therapy to target new phenotypes that arise. In breast cancer, viral interactions cause chromosome breakage and fusion, producing anaphase bridges containing two centromeres. Cell division then leads to more chromosome breaks, more mutations, and new cancer phenotypes. Analyses of public data from about 2100 breast cancers implicate Epstein-Barr virus (EBV or Human Herpes Virus 4) in the cycles that break chromosomes and produce anaphase bridges. EBV generates breakage cycles with highly fragmented chromosomes (chromothripsis) as an inherent feature. These episodes continually produce new cancer driver mutations. DNA at chromothripsis boundaries in breast cancers consistently matches EBV tumor variants. Breakpoints are often near known EBV binding sites, and breast cancer shares an EBV methylation signature with other known EBV cancers. Breast cancer breakpoints on all chromosomes cluster around the same positions as in cancers definitively associated with EBV infection (nasopharyngeal cancer and Burkitts lymphoma). The status of EBV in the microbiome and the presence of anaphase bridges are phenotypic markers of breast cancer that are not explicitly targeted by chemotherapy. The breast cancer phenome continually supplies new cancer driver genes that foster breast cancer recurrence.

SummaryBreaks in breast cancer chromosomes and in fragments of anaphase bridges occur mainly near regions that resemble Epstein-Barr virus (EBV) tumor variants

Extensive breaks are an inherent outcome of anaphase bridges and the reassembly of highly fragmented chromosomes continually generates new cancer driver mutations and new cancer phenotypes.

Treatment does not explicitly target the underlying breakage and reassembly cycles so breast cancer often comes back.

The breast cancer phenome is not fixed and can continuously generate new phenotypes

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=154 HEIGHT=200 SRC="FIGDIR/small/467751v4_ufig1.gif" ALT="Figure 1">
View larger version (38K):
org.highwire.dtl.DTLVardef@6537cforg.highwire.dtl.DTLVardef@132fa50org.highwire.dtl.DTLVardef@1cfcb04org.highwire.dtl.DTLVardef@11979aa_HPS_FORMAT_FIGEXP  M_FIG C_FIG Development of a yeast-based assay for bioavailable phosphorus Preventing eutrophication of inland freshwater ecosystems requires quantifying the phosphorus (P) content of the streams and rivers that feed them. Typical methods for measuring P assess soluble reactive P (SRP) or total P (TP) and require expensive analytical techniques that produce hazardous waste. Here we present a novel method for measuring the more relevant bioavailable P (BAP); this assay utilizes the growth of familiar bakers yeast, avoids production of hazardous waste, and reduces cost relative to measurements of SRP and TP. The yeast BAP (yBAP) assay takes advantage of the observation that yeast density at saturating growth increases linearly with provided P. We show that this relationship can be used to measure P in freshwater in concentration ranges relevant to eutrophication. In addition, we measured yBAP in water containing known amounts of fertilizer and in samples from agricultural waterways. We observed that the majority of yBAP values were between those obtained from standard SRP and TP measurements, demonstrating that the assay is compatible with real-world settings. The cost-effective and nonhazardous nature of the yeast-based assay suggests that it could have utility in a range of settings, offering added insight to identify water systems at risk of eutrophication from excess phosphorus. GPI-anchored FGF directs cytoneme-mediated bidirectional signaling to self-regulate tissue-specific dispersion How signaling proteins generate a multitude of information to organize tissue patterns is critical to understanding morphogenesis. In Drosophila, FGF produced in wing-disc cells regulates the development of the disc-associated air-sac-primordium (ASP). Here, we show that FGF is Glycosylphosphatidylinositol-anchored to the producing cell surface and that this modification both inhibits free FGF secretion and activates target-specific bidirectional FGF-FGFR signaling through cytonemes. FGF-source and recipient ASP cells extend cytonemes that present FGF and FGFR on their surfaces and reciprocally recognize each other over distance by contacting through CAM-like FGF-FGFR binding. Contact-mediated FGF-FGFR binding induces bidirectional signaling, which, in turn, promotes ASP and source cells to polarize cytonemes toward each other and reinforce signaling contacts. Subsequent un-anchoring of FGFR-bound-FGF from the source cell membrane dissociates cytoneme contacts and delivers FGF target-specifically to ASP cytonemes for paracrine functions. Thus, GPI-anchored FGF organizes both source and recipient cells and self-regulates its cytoneme-mediated tissue-specific dispersion and signaling. Improving Martini 3 for disordered and multidomain proteins Coarse-grained molecular dynamics simulations are a useful tool to determine conformational ensembles of proteins. Here, we show that the coarse-grained force field Martini 3 underestimates the global dimensions of intrinsically disordered proteins (IDPs) and multidomain proteins when compared with small angle X-ray scattering (SAXS) data, and that increasing the strength of protein-water interactions favours more expanded conformations. We find that increasing the strength of interactions between protein and water by ca. 10% results in improved agreement with the SAXS data for IDPs and multi-domain proteins. We also show that this correction results in a more accurate description of self-association of IDPs and folded proteins and better agreement with paramagnetic relaxation enhancement data for most IDPs. While simulations with this revised force field still show deviations to experiments for some systems our results suggest that it is overall a substantial improvement for coarse-grained simulations of soluble proteins. A chromosome-level genome assembly of the European Beech (Fagus sylvatica) reveals anomalies for organelle DNA integration, repeat content and distribution of SNPs The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 Mb and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb of the genome in which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome- wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate. An alternative UPF1 isoform drives conditional remodeling of nonsense-mediated mRNA decay The nonsense-mediated mRNA decay (NMD) pathway monitors translation termination to degrade transcripts with premature stop codons and regulate thousands of human genes. Here we show that an alternative mammalian-specific isoform of the core NMD factor UPF1, termed UPF1LL, enables condition-dependent remodeling of NMD specificity. Previous studies indicate that the extension of a conserved regulatory loop in the UPF1LL helicase core confers a decreased propensity to dissociate from RNA upon ATP hydrolysis relative to the major UPF1 isoform, designated UPF1 SL. Using biochemical and transcriptome-wide approaches, we find that UPF1LL overcomes the protective RNA binding proteins PTBP1 and hnRNP L to preferentially bind and down-regulate long 3UTRs normally shielded from NMD. Unexpectedly, UPF1LL supports induction of NMD on new populations of substrate mRNAs in response to activation of the integrated stress response and impaired translation efficiency. Thus, while canonical NMD is abolished by moderate translational repression, UPF1LL activity is enhanced, providing a mechanism to rapidly rewire NMD specificity in response to cellular stress. Probing Patterning in Microbial Consortia with picCASO: a Cellular Automaton for Spatial Organisation Microbial consortia exhibit spatial patterning across diverse environments. Since probing the self-organization of natural microbial communities is limited by their inherent complexity, synthetic models have emerged as attractive alternatives. In this study, we develop novel frameworks of bacterial communication and explore the emergent spatiotemporal organization of microbes. Specifically, we build quorum sensing-mediated models of microbial growth that are utilized to characterize the dynamics of communities from arbitrary initial configurations and establish the effectiveness of our communication strategies in coupling the growth rates of microbes. Our simulations indicate that the behavior of quorum sensing-coupled consortia can be most effectively modulated by the rates of secretion of AHLs. Such a mechanism of control enables the construction of desired relative populations of constituent species in spatially organized populations. Our models accurately recapitulate previous experiments that have investigated pattern formation in synthetic multi-cellular systems. Additionally, our software tool enables the easy implementation and analysis of our frameworks for a variety of initial configurations and simplifies the development of sophisticated gene circuits facilitating distributed computing. Overall, we demonstrate the potential of spatial organization as a tunable parameter in synthetic biology by introducing a communication paradigm based on the location and strength of coupling of microbial strains.

Author SummaryInteracting microorganisms that coexist in a given environment tend to have well-defined spatial arrangements. While the emergence of such organization is seen across different microbiomes in nature, it is hitherto not well understood. Decoding the inherent spatial patterning of microbes is constrained by the complexity of their natural habitats. Here, we take advantage of synthetic models of microbial communities to study the dynamics of emergent spatial organization. Our framework of bacterial communication utilizes modular synthetic devices to couple growth rates. In addition to uncovering potential principles of spatial organization, this work enables the construction of complex genetic circuits distributed across communicating strains. Additionally, we have developed a software tool, picCASO, that streamlines the investigation of microbial systems communicating through such frameworks. Discovery of an Antarctic ascidian-associated uncultivated Verrucomicrobia with antimelanoma palmerolide biosynthetic potential The Antarctic marine ecosystem harbors a wealth of biological and chemical innovation that has risen in concert over millennia since the isolation of the continent and formation of the Antarctic circumpolar current. Scientific inquiry into the novelty of marine natural products produced by Antarctic benthic invertebrates led to the discovery of a bioactive macrolide, palmerolide A, that has specific activity against melanoma and holds considerable promise as an anticancer therapeutic. While this compound was isolated from the Antarctic ascidian Synoicum adareanum, its biosynthesis has since been hypothesized to be microbially mediated, given structural similarities to microbially-produced hybrid non-ribosomal peptide-polyketide macrolides. Here, we describe a metagenome-enabled investigation aimed at identifying the biosynthetic gene cluster (BGC) and palmerolide A-producing organism. A 74 Kbp candidate BGC encoding the multi-modular enzymatic machinery (hybrid Type I-trans-AT polyketide synthase-non-ribosomal peptide synthetase and tailoring functional domains) was identified and found to harbor key features predicted as necessary for palmerolide A biosynthesis. Surveys of ascidian microbiome samples targeting the candidate BGC revealed a high correlation between palmerolide-gene targets and a single 16S rRNA gene variant (R=0.83 - 0.99). Through repeated rounds of metagenome sequencing followed by binning contigs into metagenome-assembled genomes, we were able to retrieve a near-complete genome (10 contigs) of the BGC-producing organism, a novel verrucomicrobium within the Opitutaceae family that we propose here as Candidatus Synoicihabitans palmerolidicus. The refined genome assembly harbors five highly similar BGC copies, along with structural and functional features that shed light on the host-associated nature of this unique bacterium.

ImportancePalmerolide A has potential as a chemotherapeutic agent to target melanoma. We interrogated the microbiome of the Antarctic ascidian, Synoicum adareanum, using a cultivation-independent high-throughput sequencing and bioinformatic strategy. The metagenome-encoded biosynthetic machinery predicted to produce palmerolide A was found to be associated with the genome of a member of the S. adareanum core microbiome. Phylogenomic analysis suggests the organism represents a new deeply-branching genus, Candidatus Synoicihabitans palmerolidicus, in the Opitutaceae family of the Verrucomicrobia phylum. The Ca. S. palmerolidicus 4.29 Mb genome encodes a repertoire of carbohydrate-utilizing and transport pathways enabling its ascidian-associated lifestyle. The palmerolide-producers genome also contains five distinct copies of the large palmerolide biosynthetic gene cluster that may provide structural complexity of palmerolide variants. Genetic analysis of mitochondrial DNA copy number and associated traits identifies loci implicated in nucleotide metabolism, platelet activation, and megakaryocyte proliferation, and reveals a causal association of mitochondrial function with mortality Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses were used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p = 3.09 x 10-15) and the gene ontology (GO) terms for mtDNA metabolism (p = 1.43 x 10-8) and mtDNA replication (p = 1.2 x 10-7). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell related traits, kidney function, liver function and overall (p = 0.044) and non-cancer mortality (p = 6.56 x 10-4). Common Dynamic Determinants Govern Quorum Quenching Activity in N-terminal Serine Hydrolases Growing concerns about microbial antibiotic resistance have motivated extensive research into ways of overcoming antibiotic resistance. Quorum quenching (QQ) processes disrupt bacterial communication via quorum sensing, which enables bacteria to sense the surrounding bacterial cell density and markedly affects their virulence. Due to its indirect mode of action, QQ is believed to exert limited pressure on essential bacterial functions and may thus avoid inducing resistance. Although many enzymes display QQ activity against various bacterial signaling molecules, their mechanisms of action are poorly understood, limiting their potential optimization as QQ agents. Here we evaluate the capacity of three N-terminal serine hydrolases to degrade N-acyl homoserine lactones that serve as signaling compounds for Gram-negative bacteria. Using molecular dynamics simulations of the free enzymes and their complexes with two signaling molecules of different lengths, followed by quantum mechanics/molecular mechanics molecular dynamics simulations of their initial catalytic steps, we clarify the molecular processes underpinning their QQ activity. We conclude that all three enzymes degrade bacterial signaling molecules via similar reaction mechanisms. Moreover, we experimentally confirmed the activity of two penicillin G acylases from Escherichia coli (ecPGA) and Achromobacter spp. (aPGA), adding these biotechnologically well-optimized enzymes to the QQ toolbox. We also observed enzyme- and substrate-dependent differences in the catalytic actions of these enzymes, arising primarily from the distinct structures of their acyl-binding cavities and the dynamics of their molecular gates. As a consequence, the first reaction step catalyzed by ecPGA with a longer substrate had an elevated energy barrier because its shallow acyl binding site could not accommodate a productive substrate-binding configuration. Conversely, aPGA in complex with both substrates exhibited unfavorable energetics in both reaction steps due to the dynamics of the residues gating the acyl binding cavity entrance. Finally, the energy barriers of the second reaction step catalyzed by Pseudomonas aeruginosa acyl-homoserine lactone acylase with both substrates were higher than in the other two enzymes due to the unique positioning of Arg297{beta} in this enzyme. The discovery of these dynamic determinants will guide future efforts to design robust QQ agents capable of selectively controlling virulence in resistant bacterial species. A novel method for producing functionalized vesicles that efficiently deliver oligonucleotides in vitro in cancer cells and in vivo in mice Nano-based delivery systems have enhanced our ability to administer and target drugs and macromolecules to their targets. Oligonucleotide drugs have great therapeutic potential but often have off-target effects and stability issues. Therefore, they are often encapsulated in vesicles with targeting ligands such as antibodies (Ab) to deliver their cargo. Herein, we describe a novel, scalable and straightforward approach to producing functionalized vesicles called the "Functionalized Lipid Insertion Method," which differs from the older "Detergent-Dialysis Method." The latter method required excess detergent and extensive dialysis over many hours to produce the functionalized vesicles. With our method, only the functionalized lipid is detergent-solubilized during the engineering of the vesicle. The approach reduces the dialysis time, keeps the vesicles intact while orienting the targeting moieties of the functionalized lipid toward the outside of the vesicle. Pilot in vitro and in vivo experiments was performed to show the feasibility of our method. Dynamic light scattering (DLS) experiments suggested that the original vesicular structure was relatively unperturbed, and the functionalized lipid was inserted externally. Our approach efficiently delivered oligonucleotides and affected the function of liver cancer HepG2 cells. Furthermore, functionalized vesicles achieved targeted delivery of oligonucleotides in mice without inducing a significant innate immune response. The industrial and therapeutic significance and implications of functionalized vesicles produced by our method are also discussed. Additional experiments and analyses are recommended to bring out the full potential of this molecular delivery technology. BRD2 inhibition blocks SARS-CoV-2 infection by reducing transcription of the host cell receptor ACE2 SARS-CoV-2 infection of human cells is initiated by the binding of the viral Spike protein to its cell-surface receptor ACE2. We conducted a targeted CRISPRi screen to uncover druggable pathways controlling Spike protein binding to human cells. We found that the protein BRD2 is required for ACE2 transcription in human lung epithelial cells and cardiomyocytes, and BRD2 inhibitors currently evaluated in clinical trials potently block endogenous ACE2 expression and SARS-CoV-2 infection of human cells, including those of human nasal epithelia. Moreover, pharmacological BRD2 inhibition with the drug ABBV-744 inhibited SARS-CoV-2 replication in Syrian hamsters. We also found that BRD2 controls transcription of several other genes induced upon SARS-CoV-2 infection, including the interferon response, which in turn regulates the antiviral response. Together, our results pinpoint BRD2 as a potent and essential regulator of the host response to SARS-CoV-2 infection and highlight the potential of BRD2 as a novel therapeutic target for COVID-19. Chimeric mRNA based COVID-19 vaccine induces protective immunity against Omicron and Delta The emerging SARS-CoV-2 variants of concern (VOCs) exhibit enhanced transmission and immune escape, reducing the efficacy and effectiveness of the two FDA-approved mRNA vaccines. Here, we explored various strategies to develop novel mRNAs vaccines to achieve safer and wider coverage of VOCs. Firstly, we constructed a cohort of mRNAs that feature a furin cleavage mutation in the spike (S) protein of predominant VOCs, including Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1) and Delta (B.1.617.2). Not present in the mRNA vaccines currently in use, the mutation abolished the cleavage between the S1 and S2 subunits, potentially enhancing the safety profile of the immunogen. Secondly, we systematically evaluated the induction of neutralizing antibodies (nAb) in vaccinated mice, and discovered that individual VOC mRNAs elicited strong neutralizing activity in a VOC-specific manner. Thirdly, the IgG produced in mice immunized with Beta-Furin and Washington (WA)-Furin mRNAs showed potent cross-reactivity with other VOCs, which was further corroborated by challenging vaccinated mice with the live virus of VOCs. However, neither WA-Furin nor Beta-Furin mRNA elicited strong neutralizing activity against the Omicron variant. Hence, we further developed an Omicron-specific mRNA vaccine that restored protection against the original and the sublineages of Omicron variant. Finally, to broaden the protection spectrum of the new Omicron mRNA vaccine, we tested the concept of bivalent immunogen. Instead of just fusing two RBDs head-to-tail, we for the first time constructed an mRNA-based chimeric immunogen by introducing the RBD of Delta variant into the entire S antigen of Omicron. The resultant chimeric mRNA was capable of inducing potent and broadly acting nAb against Omicron (both BA.1 and BA.2) and Delta, which paves the way to develop new vaccine candidate to target emerging variants in the future.