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.
Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression.
NOMA: a high-throughput microchip for robust, sequential measurements of secretions from the same single-cells Despite recent advances in single-cell analysis technologies, lacking simple methods to keep the live single-cells traceable for longitudinal detection reliably poses a significant obstacle in single-cell secretion analysis. Here we reported the high-density NOMA (narrow-opening microwell array) microchip that realized the retention of[&ge;] 97% of trapped single cells in dedicated spatial locations during repetitive detection procedures, verified with both adherent and suspension cells by two researchers independently. We applied it in monitoring single-cell protein secretions sequentially from the same single cells, and we found the digital protein secretion patterns dominate the protein secretion. We also demonstrated the microchip for longitudinally tracking IL-8 and the CD81+EV secretions from the same single-cells over days, which revealed the presence of "super secretors" within the cell population be more persistent to secrete protein or extracellular vesicle for an extended period. The NOMA platform reported here is simple, robust, and easy to operate for tracking sequential measurements from the same single cells, representing a novel and informative tool to inspire new observations in biomedical research.

Table of Contents Graphic

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=66 SRC="FIGDIR/small/427049v4_ufig1.gif" ALT="Figure 1">
View larger version (25K):
org.highwire.dtl.DTLVardef@aa69c7org.highwire.dtl.DTLVardef@6bc79org.highwire.dtl.DTLVardef@3fcbbforg.highwire.dtl.DTLVardef@597f4b_HPS_FORMAT_FIGEXP  M_FIG C_FIG G2 stem cells orchestrate time-directed, long-range coordination of calcium signaling during skin epidermal regeneration Skin epidermal homeostasis is maintained via constant regeneration by stem cells, which must communicate to balance their self-renewal and differentiation. A key molecular pathway, Ca2+ signaling has been implicated as a signal integrator in developing and wounded epithelial tissues[1, 2, 3, 4]. Yet how stem cells carry out this signaling across a regenerative tissue remains unknown due to significant challenges in studying signaling dynamics in live mice, limiting our understanding of the mechanisms of stem cell communication during homeostasis. To interpret high dimensional signals that have complex spatial and temporal patterns, we combined optimized imaging of Ca2+ signaling in thousands of epidermal stem cells in living mice with a new machine learning tool, Geometric Scattering Trajectory Homology (GSTH). Using a combination of signal processing, data geometry, and topology, GSTH captures patterns of signaling at multiple scales, either between direct or distant stem cell neighbors. Here we show that epidermal stem cells display dynamic intercellular Ca2+ signaling among neighborhoods of up to 10 cells that is surprisingly coordinated and directed through time across a pool of thousands of stem cells. We find that this collective coordination is an emergent property of the stem cell compartment, distinct from excitatory quiescent neuronal tissues. We demonstrate that cycling stem cells, specifically G2 cells, govern homeostatic patterns of Ca2+ signaling. Stem cells in different cell cycle stages dynamically regulate localization of the gap junction component Connexin43 (Cx43). Lastly, we uncouple global from local communication and identify Cx43 as the molecular mediator necessary for connectivity between local signaling neighborhoods. This work provides resolution in how stem cells at different stages of the cell cycle communicate and how that diversity of phases is essential for tissue wide communication and signaling flow during epidermal regeneration. Our approach provides a framework to investigate stem cell populations and their signaling dynamics, previously not possible. Chimeric spike mRNA vaccines protect against Sarbecovirus challenge in mice The emergence of SARS-CoV in 2003 and SARS-CoV-2 in 2019 highlights the need to develop universal vaccination strategies against the broader Sarbecovirus subgenus. Using chimeric spike designs, we demonstrate protection against challenge from SARS-CoV, SARS-CoV-2, SARS-CoV-2 B.1.351, bat CoV (Bt-CoV) RsSHC014, and a heterologous Bt-CoV WIV-1 in vulnerable aged mice. Chimeric spike mRNAs induced high levels of broadly protective neutralizing antibodies against high-risk Sarbecoviruses. In contrast, SARS-CoV-2 mRNA vaccination not only showed a marked reduction in neutralizing titers against heterologous Sarbecoviruses, but SARS-CoV and WIV-1 challenge in mice resulted in breakthrough infection. Chimeric spike mRNA vaccines efficiently neutralized D614G, UK B.1.1.7., mink cluster five, and the South African B.1.351 variant of concern. Thus, multiplexed-chimeric spikes can prevent SARS-like zoonotic coronavirus infections with pandemic potential.

SentenceChimerized RBD, NTD, and S2 spike mRNA-LNPs protect mice against epidemic, zoonotic, and pandemic SARS-like viruses Mechanical bistability of the mesoderm epithelium facilitates mesoderm invagination during Drosophila gastrulation Apical constriction driven by non-muscle myosin II ("myosin") provides a well-conserved mechanism to mediate epithelial folding. It remains unclear how contractile forces near the apical surface of a cell sheet drive out-of-the-plane bending of the sheet and whether myosin contractility is required throughout folding. By optogenetic-mediated acute inhibition of myosin, we find that during Drosophila mesoderm invagination, myosin contractility is critical to prevent tissue relaxation during the early, "priming" stage of folding but is dispensable for the actual folding step after the tissue passes through a stereotyped transitional configuration. The binary response suggests that the mesoderm is mechanically bistable during gastrulation. Combined modeling analysis and experimental measurements suggest that the observed mechanical bistability may arise from apicobasal shrinkage of the surrounding ectoderm, which promotes mesoderm invagination by facilitating a buckling transition. Our results suggest that Drosophila mesoderm invagination requires a joint action of local myosin contractility and mechanical bistability of the epithelium to trigger epithelial buckling. Multi-omics investigation of tacrolimus nephrotoxicity BackgroundTacrolimus, an immunosuppressive drug prescribed to a majority of transplanted patients is nephrotoxic, through still unclear mechanisms. This study aims to evaluate the impact of tacrolimus on a lineage of proximal tubular cells using a multi-omics approach.

MethodsLLC-PK1 cells were exposed to 5 M of tacrolimus for 24h. Intracellular proteins and metabolites, and extracellular metabolites were extracted and analysed by LC-MS/MS. The transcriptional expression of the dysregulated proteins PCK-1, FBP1 and FBP2 was measured using RT-qPCR.

ResultsIn our cell model, tacrolimus impacted different metabolic pathways including those of arginine (e.g., citrulline, ornithine) (p < 0.0001), amino acids (e.g., valine, isoleucine, aspartic acid) (p < 0.0001) and pyrimidine (p<0.01). In addition, it induced oxidative stress (p < 0.01) as shown by a decrease in total cell glutathione quantity. It impacted cell energy through an increase in Krebs cycle intermediates (e.g., citrate, aconitate, fumarate) (p < 0.01) and down-regulation of PCK-1 (p < 0.05) and FPB1 (p < 0.01), which are key enzymes in gluconeogenesis. Apart from glucose synthesis, gluconeogenesis is an important process in kidney mediated acid-base balance control.

ConclusionThe variations found using this multi-omics approach clearly point towards a dysregulation of energy production in epithelial cells of the renal tubule, and potentially of their functions, that may be implicated in tacrolimus nephrotoxicity in the clinics. Mutagenic distinction between the receptor-binding and fusion subunits of the SARS-CoV-2 spike glycoprotein and its upshot We observe that a residue R of the spike glycoprotein of SARS-CoV-2 which has mutated in one or more of the current Variants of Concern or Interest or under Monitoring rarely participates in a backbone hydrogen bond if R lies in the S1 subunit and usually participates in one if R lies in the S2 subunit. A partial explanation for this based upon free energy is explored as a potentially general principle in the mutagenesis of viral glycoproteins. This observation could help target future vaccine cargos for the evolving coronavirus as well as more generally. A study of the Delta and Omicron variants suggests that Delta was an energetically necessary intermediary in the evolution from Wuhan-Hu-1 to Omicron. A multi-dataset evaluation of frame censoring for motion correction in task-based fMRI Subject motion during fMRI can affect our ability to accurately measure signals of interest. In recent years, frame censoring--that is, statistically excluding motion-contaminated data within the general linear model using nuisance regressors--has appeared in several task-based fMRI studies as a mitigation strategy. However, there have been few systematic investigations quantifying its efficacy. In the present study, we compared the performance of frame censoring to several other common motion correction approaches for task-based fMRI using open data and reproducible workflows. We analyzed eight publicly-available datasets representing eleven distinct tasks in child, adolescent, and adult participants. Performance was quantified using maximum t-values in group analyses, and ROI-based mean activation and split-half reliability in single subjects. We compared frame censoring across several thresholds to the use of 6 and 24 canonical motion regressors, wavelet despiking, robust weighted least squares, and untrained ICA-based denoising, for a total of 240 separate analyses. Thresholds used to identify censored frames were based on both motion estimates (FD) and image intensity changes (DVARS). Relative to standard motion regressors, we found consistent improvements for modest amounts of frame censoring (e.g., 1-2% data loss), although these gains were frequently comparable to what could be achieved using other techniques. Importantly, no single approach consistently outperformed the others across all datasets and tasks. These findings suggest that the choice of a motion mitigation strategy depends on both the dataset and the outcome metric of interest. The DC1 domain protein BINUCLEATE POLLEN is required for pollen development in Arabidopsis thaliana Development of the male gametophyte is a tightly regulated process that requires precise control of cell division and gene expression. A relevant aspect to understand the events underlying pollen development regulation constitutes the identification and characterization of the genes required for this process. In this work we showed that the DC1 domain protein BINUCLEATE POLLEN (BNP) is essential for pollen development and germination. Pollen grains carrying the defective BNP allele failed to complete mitosis II and are impaired in pollen germination. By yeast two hybrid screening and bimolecular fluorescence complementation assays, we identified a set of BNP-interacting proteins. Among confirmed interactors we found NAC family transcriptional regulators Vascular plant One-Zinc finger 1 (VOZ1) and VOZ2. VOZ1 localization changes during pollen development, moving to the vegetative nucleus at the tricellular stage. We observed that this relocalization requires BNP, as in their absence in pollen from bnp/BNP plants, VOZ1 nuclear localization resulted impaired. As voz1voz2 double mutants showed the same developmental defect observed in bnp pollen grains, we propose that BNP requirement to complete microgametogenesis could be linked to its interaction with VOZ1/2 proteins. BNP could have a role of scaffold protein, recruiting VOZ1/2 to the endosomal system into assemblies that are required for their further translocation to the nucleus, where they act as transcriptional regulators. 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. 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. Outer hair cell electro-mechanical power conversion Outer hair cells are the cellular motors in the mammalian inner ear responsible for sensitive high-frequency hearing. Motor function requires expression of the protein prestin in the OHC lateral membrane and mechano-electrical transduction in the apical hair bundle. In the present report, electrical power consumption and mechanical power output are determined using previously published voltage clamp data from isolated OHCs and membrane patches. Results reveal that power output by the prestin-motor complex peaks at a best frequency much higher than implied by the low-pass character of nonlinear capacitance, and much higher than the whole-cell resistive-capacitive corner frequency. High frequency power output is enabled by a -90{degrees} shift in the phase of prestin-dependent electrical charge displacement, manifested electrically as emergence of imaginary-valued nonlinear capacitance. Augmenting Recombinant Antibody Production in HEK293E Cells: Optimising Transfection and Culture Parameters BackgroundOptimising recombinant antibody production is important for cost-effective therapeutics and diagnostics. With impact on commercialisation, higher productivity beyond laboratory scales is highly sought, where efficient production can also accelerate antibody characterisations and investigations.

MethodsInvestigating HEK293E cells for mammalian antibody production, various transfection and culture parameters were systematically analysed for antibody light chain production before evaluating them for whole antibody production. Transfection parameters investigated include seeding cell density, the concentration of the transfection reagent and DNA, complexation time, temperature, and volume, as well as culture parameters such as medium replacement, serum deprivation, use of cell maintenance antibiotic, incubation temperature, medium volume, post-transfection harvest day and common nutrient supplements.

ResultsUsing 2 mL adherent HEK293E cell culture transfections with 25 kDa linear Polyethylenimine in the most optimised parameters, we demonstrated a ~2-fold production increase for light chain alone and for whole antibody production reaching 536 and 49 g respectively in a cost-effective manner. With the addition of peptone, {kappa} light chain increased by ~4-fold to 1032 g while whole antibody increased to a lesser extent by ~2.5-fold to 51 g, with benefits potentially for antibodies limited by their light chains in production.

ConclusionsOur optimised findings show promise for a more efficient and convenient antibody production method through transfection and culture optimisations that can be incorporated to scale up processes and with potential transferability to other mammalian-based recombinant protein production using HEK293E cells.

Statement of SignificanceRecombinant antibody production is crucial for antibody research and development. Systematically investigating transfection and culture parameters such as PEI/DNA concentrations, complexation time, volume, and temperature, supplements, etc., we demonstrated a ~4-fold light chain alone production increase to 1032 g and a 2.5-fold whole antibody production increase to 51 g from 2 mL transfections. A variable-gain stochastic pooling motif mediates information transfer from receptor assemblies into NF-κB A myriad of inflammatory cytokines regulate signaling pathways to maintain cellular homeostasis. The IKK complex is an integration hub for cytokines that govern NF-{kappa}B signaling. In response to inflammation, IKK is activated through recruitment to receptor-associated protein assemblies. How and what information IKK complexes transmit about the milieu are open questions. Here we track dynamics of IKK complexes and nuclear NF-{kappa}B to identify upstream signaling features that determine same-cell responses. Experiments and modeling of single complexes reveals their size, number, and timing relays cytokine-specific information with feedback control that is independent of transcription. Our results provide evidence for variable-gain stochastic pooling, a noise-reducing motif that enables parsimonious and cytokine-specific information transfer. We propose that emergent properties of stochastic pooling are general principles of receptor signaling.

One Sentence SummaryA variable-gain stochastic pooling motif mediates robust and tunable information transmission from the extracellular milieu into the cell. BRPF3 knockdown or inhibition moderately reverses olaparib resistance in high grade serous ovarian carcinoma BackgroundPARP inhibitors (PARPi) kill cancer cells by stalling DNA replication and preventing DNA repair, resulting in a critical accumulation of DNA damage. Resistance to PARPi is a growing clinical problem in the treatment of high grade serous ovarian carcinoma (HGSOC). Acetylation of histone H3 lysine 14 (H3K14ac) and associated histone acetyltransferases (HATs) have known functions in DNA repair and replication, but their expression and activities have not been examined in the context of PARPi-resistant HGSOC.

ResultsUsing mass spectrometry profiling of histone modifications, we observed altered H3K14ac enrichment in PARPi-resistant HGSOC cells relative to isogenic PARPi-sensitive lines. By RT-qPCR and RNA-Seq, we also observed altered expression of numerous HATs in PARPi-resistant HGSOC cells and a PARPi-resistant PDX model. Knockdown of HATs only modestly altered PARPi response, although knockdown and inhibition of PCAF significantly increased resistance. Pharmacologic inhibition of HBO1 severely depleted H3K14ac but did not affect PARPi response. However, knockdown and inhibition of BRPF3, which is known to interact in a complex with HBO1, did reduce PARPi resistance.

ConclusionsThis study demonstrates that severe depletion of H3K14ac does not affect PARPi response in HGSOC. Our data suggest that bromodomain functions of HAT proteins such as PCAF, or accessory proteins such as BRPF3, may play a greater role in PARPi response than acetyltransferase functions. Failure to account for behavioral variability significantly compromises accuracy in indirect population monitoring Indirect wildlife population surveying largely depends upon counts of artefacts of behavior (e.g., nests or dung). Likelihood to encounter these artefacts is derived from both artefact production and decay, and variability in artefact production behavior is considered to contribute minimally to inaccuracy in wildlife estimation. Here we demonstrate how ignoring behavioral variability contributes to significant population misestimation, using an example of an endangered ape, the bonobo (Pan paniscus). Until now, a single estimate of (i.e., nest construction) has been used to extrapolate bonobo densities, assumed to be broadly representative of bonobo sign production behavior. We estimated nest construction rate across seasons and social groups at the Kokolopori Bonobo Reserve, DRC and find that nest construction rates in bonobos to be highly variable across populations as well as seasonal. Failure to account for this variability led to potentially severe degradation in accuracy of bonobo population estimates of abundance, accounting for a likely overestimation of bonobo numbers by 34%, and at worst as high as 80%. With this example, we demonstrate that failure to account for inter- and intra-population behavioral variation compromises the ability to estimate both relative and absolute wildlife abundances. We argue that variation in sign production is but one of several potential ways that behavioral variability can affect conservation monitoring, should be measured across contexts whenever possible, and must be considered in population estimation confidence intervals. With increasing attention to behavioral variability as a potential tool for conservation, conservationists must also account for the impact that behavioral variability can play upon wildlife population estimation. Our results underline the importance of observational research to wildlife monitoring schemes as a critical but overlooked component of conservation management. We extensively discuss the various avenues through which behavioral variability is likely to impact wildlife monitoring accuracy and precision and propose potential avenues for accounting for behavioral variability in wildlife monitoring. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. Natural selection promotes the evolution of recombination 2: during the process of natural selection The ubiquity of sex and recombination in nature has eluded unified explanation since the time of Darwin. Conditions that promote the evolution of recombination, broadly defined as any form of genetic mixing, are fairly well understood: it is favored when genomes tend to contain more selectively mismatched combinations of alleles than can be explained by chance alone. Yet, while a variety of theoretical approaches have been put forth to explain why such conditions would prevail in natural populations, each has turned out to be of limited scope and applicability. Here, we show, simply and surprisingly, that natural selection acting on standing heritable variation always creates conditions favoring the evolution of recombination, in expectation. Specifically, we find that, in expectation: 1) the mean selective advantage of recombinants is non-negative, 2) the mean selective advantage of a recombination-competent modifier is non-negative, and 3) the asymptotic frequency of a recombination-competent modifier is close to one and is independent of the strength of selection. Remarkably, these findings are independent of the distribution of genic fitnesses in the standing heritable variation upon which natural selection acts, implying that the source of this variation is immaterial. Taken together, our findings indicate that: 1) the evolution of recombination should be promoted in expectation wherever natural selection is operating, and 2) sex and recombination may have evolved more as a byproduct than as a catalyst of natural selection. Multivariate genome-wide association study on tissue-sensitive diffusion metrics identifies key molecular pathways for axonal growth, synaptogenesis, and astrocyte-mediated neuroinflammation The molecular determinants of tissue composition of the human brain remain largely unknown. Recent genome-wide association studies (GWAS) on this topic have had limited success due to methodological constraints. Here, we apply advanced whole-brain analyses on multi-shell diffusion imaging data and multivariate GWAS to two large scale imaging genetic datasets (UK Biobank and the Adolescent Brain Cognitive Development study) to identify and validate genetic association signals. We discovered 503 unique genetic loci that have impact on multiple regions of human brain. Among them, more than 79% were validated in either of two large-scale independent imaging datasets. Key molecular pathways involved in axonal growth, astrocyte-mediated neuroinflammation, and synaptogenesis during development were found to significantly impact the measured variations in tissue-specific imaging features. Our results shed new light on the biological determinants of brain tissue composition and their potential overlap with the genetic basis of neuropsychiatric disorders. First detection of herpesvirus and prevalence of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Development of a Novel PROTAC using the Nucleic Acid Aptamer as a Targeting Ligand for Tumor Selective Degradation of Nucleolin PROteolysis TArgeting Chimeras (PROTACs) induce targeted protein degradation by hijacking the intracellular ubiquitin proteasome system, thus emerging as a new strategy for drug development. However, most PROTACs generated lack cell-type selectivity and have poorly soluble in water. To address this drawback, we developed a novel PROTAC ZL216 using aptamer AS1411 as a targeting ligand of nucleolin to conjugate with a small molecule ligand of E3 ligase VHL, which shows high aqueous solubility and serum stability. Based on the differential expression of nucleolin on the cell surface, ZL216 could bind to and internalize into breast cancer cells, but not normal breast cells. Furthermore, we revealed that ZL216 promoted the formation of a nucleolin-ZL216-VHL ternary complex in breast cancer cells and potently induced nucleolin degradation in vitro and in vivo. As a result, ZL216 inhibited the proliferation and migration of breast cancer cells. These studies demonstrate that in addition to peptides and small molecule compounds, nuclei acid aptamers can also be used to generate PROTACs, which broadens the toolbox constructing PROTACs and provides a promising strategy for development of tumor-selective PROTACs. 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. Thermodynamic constraints on the assembly and diversity ofmicrobial ecosystems are different near to and far fromequilibrium Non-equilibrium thermodynamics has long been an area of substantial interest to ecologists because most fundamental biological processes, such as protein synthesis and respiration, are inherently energy-consuming. However, most of this interest has focused on developing coarse ecosystem-level maximisation principles, providing little insight into underlying mechanisms that lead to such emergent constraints. Microbial communities are a natural system to decipher this mechanistic basis because their interactions in the form of substrate consumption, metabolite production, and cross-feeding can be described explicitly in thermodynamic terms. Previous work has considered how thermodynamic constraints impact competition between pairs of species, but restrained from analysing how this manifests in complex dynamical systems. To address this gap, we develop a thermodynamic microbial community model with fully reversible reaction kinetics, which allows direct consideration of free-energy dissipation. This also allows species to interact via products rather than just substrates, increasing the dynamical complexity, and allowing a more nuanced classification of interaction types to emerge. Using this model, we find that community diversity increases with substrate lability, because greater free-energy availability allows for faster generation of niches. Thus, more niches are generated in the time frame of community establishment, leading to higher final species diversity. We also find that allowing species to make use of near-to-equilibrium reactions increases diversity in a low free-energy regime. In such a regime, two new thermodynamic interaction types that we identify here reach comparable strengths to the conventional (competition and facilitation) types, emphasising the key role that thermodynamics plays in community dynamics. Our results suggest that accounting for realistic thermodynamic constraints is vital for understanding the dynamics of real-world microbial communities.

Author summaryThere is a growing interest in microbial communities due to their important role in biogeochemical cycling as well as plant and animal health. Although our understanding of thermodynamic constraints on individual cells is rapidly improving, the impact of these constraints on complex microbial communities remains largely unexplored theoretically and empirically. Here, we develop a new microbial community model which allows thermodynamic efficiency and entropy production to be calculated directly. We find that availability of substrates with greater free-energy allows for a faster rate of niche generation, leading to higher final species diversity. We also show that when the free-energy availability is low, species with reactions close to thermodynamic equilibrium are favoured, leading to more diverse and efficient communities. In addition to the conventional interaction types (competition and facilitation), our model reveals the existence of two novel interaction types mediated by products rather than substrates. Though the conventional interactions are generally the strongest, the novel interaction types are significant when free-energy availability is low. Our results suggest that non-equilibrium thermodynamics need to be considered when studying microbial community dynamics. 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. Construction of a new chromosome-scale, long-read reference genome assembly for the Syrian hamster, Mesocricetus auratus BackgroundThe Syrian hamster (Mesocricetus auratus) has been suggested as a useful mammalian model for a variety of diseases and infections, including infection with respiratory viruses such as SARS-CoV-2. The MesAur1.0 genome assembly was generated in 2013 using whole-genome shotgun sequencing with short-read sequence data. Current more advanced sequencing technologies and assembly methods now permit the generation of near-complete genome assemblies with higher quality and greater continuity.

FindingsHere, we report an improved assembly of the M. auratus genome (BCM_Maur_2.0) using Oxford Nanopore Technologies long-read sequencing to produce a chromosome-scale assembly. The total length of the new assembly is 2.46 Gbp, similar to the 2.50 Gbp length of a previous assembly of this genome, MesAur1.0. BCM_Maur_2.0 exhibits significantly improved continuity with a scaffold N50 that is 6.7 times greater than MesAur1.0. Furthermore, 21,616 protein coding genes and 10,459 noncoding genes are annotated in BCM_Maur_2.0 compared to 20,495 protein coding genes and 4,168 noncoding genes in MesAur1.0. This new assembly also improves the unresolved regions as measured by nucleotide ambiguities, where approximately 17.11% of bases in MesAur1.0 were unresolved compared to BCM_Maur_2.0 in which the number of unresolved bases is reduced to 3.00%.

ConclusionsAccess to a more complete reference genome with improved accuracy and continuity will facilitate more detailed, comprehensive, and meaningful research results for a wide variety of future studies using Syrian hamsters as models.