Determining significant correlation between pairs of extant characters in a small parsimony framework When studying the evolutionary relationships among a set of species, the principle of parsimony states that a relationship involving the fewest number of evolutionary events is likely the correct one. Due to its simplicity, this principle was formalized in the context of computational evolutionary biology decades ago by, e.g., Fitch and Sankoff. Because the parsimony framework does not require a model of evolution, unlike maximum likelihood or Bayesian approaches, it is often a good starting point when no reasonable estimate of such a model is available.

In this work, we devise a method for determining if pairs of discrete characters are significantly correlated across all most parsimonious reconstructions, given a set of species on these characters, and an evolutionary tree. The first step of this method is to use Sankoffs algorithm to compute all most parsimonious assignments of ancestral states (of each character) to the internal nodes of the phylogeny. Correlation between a pair of evolutionary events (e.g., absent to present) for a pair of characters is then determined by the (co-) occurrence patterns between the sets of their respective ancestral assignments. The probability of obtaining a correlation this extreme (or more) under a null hypothesis where the events happen randomly on the evolutionary tree is then used to assess the significance of this correlation. We implement this method: parcours (PARsimonious CO-occURrenceS) and use it to identify significantly correlated evolution among vocalizations and morphological characters in the Felidae family.

The parcours tool is freely available at https://github.com/murraypatterson/parcours
Double trouble: having two hosts reduces infection prevalence in vectored trypanosomatids Trypanosomatids are a diverse family of protozoan parasites, some of which cause devastating human and livestock diseases. There are two distinct infection life-cycles in trypanosomatids; some species complete their entire life-cycle in a single host (monoxenous) while others infect two hosts (dixenous). Dixenous trypanosomatids are mostly vectored by insects, and the human trypanosomatid diseases are caused mainly by vectored parasites. While infection prevalence has been described for subsets of hosts and trypanosomatids, little is known about whether monoxenous and dixenous trypanosomatids differ in infection prevalence. Here, we use meta-analyses to synthesise all published evidence of trypanosomatid infection prevalence for the last two decades, encompassing 931 unique host-trypansomatid systems. In examining 581 studies that describe infection prevalence, we find, strikingly, that monoxenous species are two-fold more prevalent than dixenous species across all hosts. We also find that dixenous trypanosomatids have significantly lower infection prevalence in insects than their non-insect hosts. Within monoxenous trypanosomatids, genera infecting bees are most prevalent and infection prevalence does not vary between wild and managed bees. To our knowledge, these results reveal for the first time, a fundamental difference in infection prevalence according to host specificity where vectored species suffer from lower infection prevalence as a result of a  jack of all trades, master of none style trade-off between the vector and subsequent hosts.
Bamboos flower after the return of almost the same sun-moon phasing as at seedling emergence All Asian woody bamboo species of economic importance are semelparous. They remain vegetative during time intervals that are specific to each species and range from three to 120 years, with notable concentrations around a series of values (3, 7-8, 14-17, 29-36, 42-48, 61-64, and 120 years). Then, they flower gregariously within a short period. Entire forests temporarily disappear during these periods, and the dates of these dramatic events have been recorded over the last 200 years.

While searching for a correlated environmental cue, I have found that the concentrations of flowering cycles were highly correlated with the series of successive returns of almost the same sun-moon phasing as at seedling emergence. On basis of knowledge on plant photoperiod sensitivity, I hypothesize that bamboo plants i) run a lunar cellular clock that is set at the full moon, ii) retain in their cellular memory the exact sun-moon phasing of the year of their emergence as seedlings, and iii) inhibit flowering until the occurrence of a unique, species-specific sun-moon phasing that is shifted by a precise amount from the sun-moon phasing at their emergence. Recent evidence of plant responses to lunar cycles supports this hypothesis, for which experimental evidence is now anticipated.

HighlightThe long-standing enigma of the triggering of the gregarious flowering of bamboos is possibly solved. Flowering would be inhibited until the occurrence of the species-specific sun-moon phasing. Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector. Novel feature selection method via kernel tensor decomposition for improved multi-omics data analysis BackgroundFeature selection of multi-omics data analysis remains challenging owing to the size of omics datasets, comprising approximately 102-105 features. In particular, appropriate methods to weight individual omics datasets are unclear, and the approach adopted has substantial consequences for feature selection. In this study, we extended a recently proposed kernel tensor decomposition (KTD)-based unsupervised feature extraction (FE) method to integrate multi-omics datasets obtained from common samples in a weight-free manner.

MethodKTD-based unsupervised FE was reformatted as the collection of kernelized tensors sharing common samples, which was applied to synthetic and real datasets.

ResultsThe proposed advanced KTD-based unsupervised FE method showed comparative performance to that of the previously proposed KTD method, as well as tensor decomposition-based unsupervised FE, but required reduced memory and central processing unit time. Moreover, this advanced KTD method, specifically designed for multi-omics analysis, attributes P-values to features, which is rare for existing multi-omics-oriented methods.

ConclusionsThe sample R code is available at https://github.com/tagtag/MultiR/ Spatiotemporal transcriptomic atlas of mouse organogenesis using DNA nanoball patterned arrays Spatially resolved transcriptomic technologies are promising tools to study cell fate decisions in a physical microenvironment, which is fundamental for enhancing our knowledge of mammalian development. However, the imbalance between resolution, transcript capture and field of view of current methodologies precludes their systematic application to analyze relatively large and three-dimensional mid- and late-gestation mammalian embryos. Here, we combined DNA nanoball (DNB) patterned arrays and tissue RNA capture to create SpaTial Enhanced REsolution Omics-sequencing (Stereo-seq). This approach allows transcriptomic profiling of large histological sections with high resolution and sensitivity. We have applied Stereo-seq to study the kinetics and directionality of transcriptional variation in a time course of mouse organogenesis. We used this information to gain insight into the molecular basis of regional specification, neuronal migration and differentiation in the developing brain. Furthermore, we mapped the expression of a panel of developmental disease-related loci on our global transcriptomic maps to define the spatiotemporal windows of tissue vulnerability. Our panoramic atlas constitutes an essential resource to investigate longstanding questions concerning normal and abnormal mammalian development. Human genetic analyses of organelles highlight the nucleus, but not the mitochondrion, in age-related trait heritability Most age-related human diseases are accompanied by a decline in cellular organelle integrity, including impaired lysosomal proteostasis and defective mitochondrial oxidative phosphorylation. An open question, however, is the degree to which inherited variation in or near genes encoding each organelle contributes to age-related disease pathogenesis. Here, we evaluate if genetic loci encoding organelle proteomes confer greater-than-expected age-related disease risk. As mitochondrial dysfunction is a "hallmark" of aging, we begin by assessing nuclear and mitochondrial DNA loci near genes encoding the mitochondrial proteome and surprisingly observe a lack of enrichment across 24 age-related traits. Within nine other organelles, we find no enrichment with one exception: the nucleus, where enrichment emanates from nuclear transcription factors. In agreement, we find that genes encoding several organelles tend to be "haplosufficient," while we observe strong purifying selection against heterozygous protein-truncating variants impacting the nucleus. Our work identifies common variation near transcription factors as having outsize influence on age-related trait risk, motivating future efforts to determine if and how this inherited variation then contributes to observed age-related organelle deterioration. FGF8-mediated signaling regulates tooth developmental pace during odontogenesis The developing human and mouse teeth constitute an ideal model system to study the regulatory mechanism underlying organ growth control due to the fact that their teeth share highly conserved and well-characterized developmental processes and their developmental tempo vary notably. In the current study, we manipulated heterogenous recombination between human and mouse dental tissues and demonstrate that the dental mesenchyme dominates the tooth developmental tempo and FGF8 could be a critical player during this developmental process. Forced activation of FGF8 signaling in the dental mesenchyme of mice promoted cell proliferation, prevented cell apoptosis via p38 and perhaps PI3K-Akt intracellular signaling, and impelled the transition of the cell cycle from G1- to S-phase in the tooth germ, resulting in the slowdown of the tooth developmental pace. Our results provide compelling evidence that extrinsic signals can profoundly affect tooth developmental tempo and the dental mesenchymal FGF8 could be a pivotal factor in controlling developmental pace in a non-cell-autonomous manner during mammalian odontogenesis. Heat-Dependent Hairpin Melting Drives Thermo-gated TRPV1 Opening The capsaicin receptor TRPV1 is a bio-thermometer in a primary afferent sensory neuron for noxious heat detection. However, the underlying molecular mechanism is unclear. Here, a hairpin thermodynamic model, together with graph theory, was to examine a role of temperature-dependent hairpin sizes in controlling the temperature threshold and sensitivity and gating states of TRPV1. The results showed that heat-dependent melting of the largest hairpin initiates a temperature threshold for channel opening while the smaller hairpins stabilize the open state and increase the heat efficacy or sensitivity. Therefore, this model may provide a broad structural basis for the thermo-gated TRP channels. Immunogenicity and protective efficacy of a highly thermotolerant, trimeric SARS-CoV-2 receptor binding domain derivative The Receptor Binding Domain (RBD) of SARS-CoV-2 is the primary target of neutralizing antibodies. We designed a trimeric, highly thermotolerant glycan engineered RBD by fusion to a heterologous, poorly immunogenic disulfide linked trimerization domain derived from cartilage matrix protein. The protein expressed at a yield of [~]80-100 mg/liter in transiently transfected Expi293 cells, as well as CHO and HEK293 stable cell lines and formed homogeneous disulfide-linked trimers. When lyophilized, these possessed remarkable functional stability to transient thermal stress of upto 100 {degrees}C and were stable to long term storage of over 4 weeks at 37 {degrees}C unlike an alternative RBD-trimer with a different trimerization domain. Two intramuscular immunizations with a human-compatible SWE adjuvanted formulation, elicited antibodies with pseudoviral neutralizing titers in guinea pigs and mice that were 25-250 fold higher than corresponding values in human convalescent sera. Against the beta (B.1.351) variant of concern (VOC), pseudoviral neutralization titers for RBD trimer were [~] three-fold lower than against wildtype B.1 virus. RBD was also displayed on a designed ferritin-like Msdps2 nanoparticle. This showed decreased yield and immunogenicity relative to trimeric RBD. Replicative virus neutralization assays using mouse sera demonstrated that antibodies induced by the trimers neutralized all four VOC to date, namely B.1.1.7, B.1.351, P.1 and B.1.617.2 without significant differences. Trimeric RBD immunized hamsters were protected from viral challenge. The excellent immunogenicity, thermotolerance, and high yield of these immunogens suggest that they are a promising modality to combat COVID-19, including all SARS-CoV-2 VOC to date. Localization of signaling receptors maximizes cellular informationacquisition in spatially-structured natural environments Cells in natural environments like tissue or soil sense and respond to extracellular ligands with intricately structured and non-monotonic spatial distributions, sculpted by processes such as fluid flow and substrate adhesion. In this work, we show that spatial sensing and navigation can be optimized by adapting the spatial organization of signaling pathways to the spatial structure of the environment. We develop an information-theoretic framework for computing the optimal spatial organization of a sensing system for a given signaling environment. We find that receptor localization maximizes information acquisition in simulated natural contexts, including tissue and soil. Receptor localization extends naturally to produce a dynamic protocol for continuously redistributing signaling receptors, which when implemented using simple feedback, boosts cell navigation efficiency by 30-fold. Broadly, our work shows how cells can maximize the fidelity of information transfer by adapting the spatial organization of signaling molecules to the spatial structure of the environment. Sleep-dependent upscaled excitability, saturated neuroplasticity, and modulated cognition in the human brain Sleep strongly affects synaptic strength, making it critical for cognition, especially learning, and memory formation. Whether and how sleep deprivation modulates human brain physiology and cognition is not well understood. Here we examined how overnight sleep deprivation vs sufficient sleep affects (a) cortical excitability, as measured by transcranial magnetic stimulation (TMS), (b) inducibility of LTP-and-LTD-like plasticity via transcranial direct current stimulation (tDCS), and (c) learning, memory and attention. The results suggest that sleep deprivation upscales cortical excitability due to enhanced glutamate-related cortical facilitation and decreased and/or reversed GABAergic cortical inhibition. Furthermore, tDCS-induced LTP-like plasticity abolishes while LTD-like plasticity converts to excitatory LTP-like plasticity under sleep deprivation. This is associated with increased EEG theta activity, which is related to sleep pressure. Finally, we show that learning and memory formation, behavioral counterparts of plasticity, and working memory and attention, which rely on cortical excitability, are impaired during sleep deprivation. Our data suggest that upscaled brain excitability, and altered plasticity, due to sleep deprivation, are associated with impaired cognitive performance. Allopurinol Blocks The Formation And Progression Of Aortic Aneurysm In A Mouse Model Of Marfan Syndrome BackgroundThe pathogenesis and progression of aortic aneurysm in Marfan syndrome (MFS) involves dysregulated TGF-{beta} and nitric oxide signaling, altered hemodynamics, and biomechanical forces. Increasing evidence indicates that redox stress participates in MFS aortopathy, though its contribution is not well established. We reported elevated reactive oxygen species (ROS) formation and NADPH oxidase NOX4 upregulation in MFS patients and mice aortic samples. Here we address the contribution of xanthine oxidoreductase (XOR), which is transcribed as dehydrogenase (XDH) and then post-translationally modified to oxidase (XO), catabolizing mainly purines into uric acid and ROS.

Methods and ResultsWe report that in MFS patients, XOR protein levels increased in aortic samples. In MFS mice (Fbn1C1041G/+), XOR mRNA transcripts augmented in the aorta of young but not older MFS mice. Moreover, the enzymatic activity of XO form increased with respect to XDH. The administration of the XOR inhibitor allopurinol (ALO) blocked the progression of aortic root aneurysm in MFS mice. ALO was also protective when administrated before the appearance of aneurysm. MFS-induced elastic fiber fragmentation, fibrotic remodeling, nuclear translocation of pNRF2, and increased 3-nitrotyrosine levels in the aortic tunica media as well as endothelial dysfunction were all prevented by allopurinol. The large overproduction of H2O2 in MFS aorta did not occur in allopurinol-treated mice whose uric acid plasma levels remained unaltered.

ConclusionsThis study strengthens the concept that redox stress is an important determinant of aortic aneurysm formation and progression in MFS and supports a clinical trial for ALO in the pharmacological treatment of MFS aortopathy.

GRAPHICAL ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=146 SRC="FIGDIR/small/464182v3_ufig1.gif" ALT="Figure 1">
View larger version (31K):
org.highwire.dtl.DTLVardef@eab02aorg.highwire.dtl.DTLVardef@1d3b39eorg.highwire.dtl.DTLVardef@1b54d0forg.highwire.dtl.DTLVardef@1d61ad0_HPS_FORMAT_FIGEXP  M_FIG C_FIG Risky choice: probability weighting explains Independence Axiom violations in monkeys Expected Utility Theory (EUT) provides axioms for maximizing utility in risky choice. The Independence Axiom (IA) is its most demanding axiom: preferences between two options should not change when altering both options equally by mixing them with a common gamble. We tested common consequence (CC) and common ratio (CR) violations of the IA over several months in thousands of stochastic choices using a large variety of binary option sets. Three monkeys showed consistently few outright Preference Reversals (8%) but substantial graded Preference Changes (46%) between the initial preferred gamble and the corresponding altered gamble. Linear Discriminant Analysis (LDA) indicated that gamble probabilities predicted most Preference Changes in CC (72%) and CR (88%) tests. The Akaike Information Criterion indicated that probability weighting within Cumulative Prospect Theory (CPT) explained choices better than models using Expected Value (EV) or EUT. Fitting by utility and probability weighting functions of CPT resulted in nonlinear and non-parallel indifference curves (IC) in the Marschak-Machina triangle and suggested IA non-compliance of models using EV or EUT. Indeed, CPT models predicted Preference Changes better than EV and EUT models. Indifference points in out-of-sample tests were closer to CPT-estimated ICs than EV and EUT ICs. Finally, while the few outright Preference Reversals may reflect the long experience of our monkeys, their more graded Preference Changes corresponded to those reported for humans. In benefitting from the wide testing possibilities in monkeys, our stringent axiomatic tests contribute critical information about risky decision-making and serve as basis for investigating neuronal decision mechanisms. The cardiac TRPA1 channel drives calcium-mediated mechano-arrhythmogenesis SUMMARY PARAGRAPHMaintenance of cardiac function involves a regulatory loop in which electrical excitation causes the heart to contract through excitation-contraction coupling (ECC),1 and the mechanical state of the heart directly affects its electrical activity through mechano-electric coupling (MEC).2 However, in pathological states such as acute ischaemia that alter early or late electro-mechanical coordination (i.e., disturbances in ECC or repolarisation-relaxation coupling, RRC), MEC may contribute to the initiation and / or sustenance of arrhythmias (mechano-arrhythmogenesis).3 The molecular identity of specific factor(s) underlying mechano-arrhythmogenesis in acute ischaemia, however, remain undefined.4 By rapid stretch of rabbit single left ventricular cardiomyocytes, we show that upon ATP-sensitive potassium channel-induced alterations of RRC, overall vulnerability to mechano-arrhythmogenesis is increased, with mechano-sensitive5-11 transient receptor potential kinase ankyrin 1 (TRPA1) channels12 acting as the molecular driver through a Ca2+-mediated mechanism. Specifically, TRPA1 activation drives stretch-induced excitation and creates a substrate for self-sustained arrhythmias, which are maintained by increased cytosolic free Ca2+ concentration ([Ca2+]i) and spontaneous [Ca2+]i oscillations. This TRPA1-dependent mechano-arrhythmogenesis involves microtubules, and can be prevented by block of TRPA1 or buffering of [Ca2+]i. Thus, in cardiac pathologies with disturbed RRC dynamics and / or augmented TRPA1 activity, TRPA1 may represent an anti-arrhythmic target with untapped therapeutic potential.13-17 A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation. Fiber orientation downsampling compromises the computation of white matter tract-related deformation Incorporating neuroimaging-revealed structural details into finite element (FE) head models opens vast new opportunities to better understand brain injury mechanisms. Recently, growing efforts have been made to integrate fiber orientation from diffusion tensor imaging (DTI) into FE models to predict white matter (WM) tract-related deformation that is biomechanically characterized by tract-related strains. Commonly used approaches often downsample the spatially enriched fiber orientation to match the FE resolution with one orientation per element (i.e., element-wise orientation implementation). However, the validity of such downsampling operation and corresponding influences on the computed tract-related strains remain elusive. To address this, the current study proposed a new approach to integrate voxel-wise fiber orientation from one DTI atlas (isotropic resolution of 1 mm3) into FE models by embedding orientations from multiple voxels within one element (i.e., voxel-wise orientation implementation). By setting the responses revealed by the newly proposed voxel-wise orientation implementation as the reference, we evaluated the reliability of two previous downsampling approaches by examining the downsampled fiber orientation and the computationally predicted tract-related strains secondary to one concussive impact. Two FE models with varying element sizes (i.e., 6.37 {+/-} 1.60 mm and 1.28 {+/-} 0.55 mm, respectively) were incorporated. The results showed that, for the model with a large voxel-mesh resolution mismatch, the downsampled element-wise fiber orientation, with respect to its voxel-wise counterpart, exhibited an absolute deviation over 30{degrees} across the WM/gray matter interface and the pons regions. Accordingly, this orientation deviation compromised the computation of tract-related strains with normalized root-mean-square errors up to 30% and underestimated the peak tract-related strains up to 10%. For the other FE model with finer meshes, the downsampling-induced effects were lower, both on the fiber orientation and tract-related strains. Taken together, the voxel-wise orientation implementation is recommended in future studies as it leverages the DTI-delineated fiber orientation to a larger extent than the element-wise orientation implementation. Thus, this study yields novel insights on integrating neuroimaging-revealed fiber orientation into FE models and may better inform the computation of WM tract-related deformation, which are crucial for advancing the etiological understanding and computational predictability of brain injury.

Graphic abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=142 SRC="FIGDIR/small/471622v2_ufig1.gif" ALT="Figure 1">
View larger version (55K):
org.highwire.dtl.DTLVardef@12ebda9org.highwire.dtl.DTLVardef@1a5d624org.highwire.dtl.DTLVardef@8e068borg.highwire.dtl.DTLVardef@4f8ff5_HPS_FORMAT_FIGEXP  M_FIG C_FIG F1F0 ATP Hydrolysis is a Determinant of Metabolic Rate, a Correlate of Lifespan, and a Weakness of Cancer Dramatic rectal temperature decrease, in mice administered with a drug that specifically inhibits F1F0 ATP hydrolysis, suggests that F1F0 ATP hydrolysis is a major determinant of metabolic rate in vivo. Across twelve investigated species, less F1F0 ATP hydrolysis correlates with greater maximal lifespan. Specific drug inhibition of F1F0 ATP hydrolysis exerts potent anti-cancer activity in vitro. Drug development of an affinity enhanced, broadly neutralizing heavy chain-only antibody that restricts SARS-CoV-2 in rodents We have identified camelid single-domain antibodies (VHHs) that cross-neutralize SARS-CoV-1 and -2, such as VHH72, which binds to a unique highly conserved epitope in the viral receptor-binding domain (RBD) that is difficult to access for human antibodies. Here, we establish a protein engineering path for how a stable, long-acting drug candidate can be generated out of such a VHH building block. When fused to human IgG1-Fc, the prototype VHH72 molecule prophylactically protects hamsters from SARS-CoV-2. In addition, we demonstrate that both systemic and intranasal application protects hACE-2-transgenic mice from SARS-CoV-2 induced lethal disease progression. To boost potency of the lead, we used structure-guided molecular modeling combined with rapid yeast-based Fc-fusion prototyping, resulting in the affinity-matured VHH72_S56A-Fc, with subnanomolar SARS-CoV-1 and -2 neutralizing potency. Upon humanization, VHH72_S56A was fused to a human IgG1 Fc with optimized manufacturing homogeneity and silenced effector functions for enhanced safety, and its stability as well as lack of off-target binding was extensively characterized. Therapeutic systemic administration of a low dose of VHH72_S56A-Fc antibodies strongly restricted replication of both original and D614G mutant variants of SARS-CoV-2 virus in hamsters, and minimized the development of lung damage. This work led to the selection of XVR011 for clinical development, a highly stable anti-COVID-19 biologic with excellent manufacturability. Additionally, we show that XVR011 is unaffected in its neutralizing capacity of currently rapidly spreading SARS-CoV-2 variants, and demonstrate its unique, wide scope of binding across the Sarbecovirus clades. Feasibility of Administering Human Pancreatic Cancer Chemotherapy in a Spontaneous Pancreatic Cancer Mouse Model BackgroundBoth modified FOLFIRINOX (mFFX) and gemcitabine/nab-paclitaxel chemotherapy regimens have been shown to improve clinical outcomes in patients with pancreatic cancer, and are often used interchangeably as the standard of care. Preclinical studies often do not use these regimens, since administering these multiagent approaches can be difficult. In this study, we assessed the feasibility of administering these two chemotherapy regimens in spontaneous pancreatic tumors using KPC mice with the ultimate goal of advancing preclinical studies.

MethodsKPC mice were created by breeding KrasLSL-G12D/+ to Trp53fl/fl;Ptf1Cre/+, resulting in KrasLSL-G12D/+;p53fl/+;Ptf1Cre/+ mice. At 14 weeks of age, mice were palpated for spontaneous tumor growth that was verified using ultrasounds. Mice with tumors under 15 mm in diameter were used. The mice were assigned to one of seven treatment regimens: 1 cycle of mFOLFIRINOX (FFX X1), 2 cycles of mFOLFIRINOX (FFX X2), 1 cycle of mFOLFIRINOX with 40Gy SBRT (FFX SBRT), 1 cycle of gemcitabine/nab-paclitaxel (GEM/AB X1), 2 cycles of gemcitabine/nab-paclitaxel (GEM/AB X2), 2 cycles of gemcitabine/nab-paclitaxel with 40Gy SBRT (GEM/AB SBRT), or no chemotherapy treatment (control).

ResultsIn total, 92 mice were included. The median OS in the FFX X2 group was slightly longer that the median OS in the FFX X1 group (15 days vs 11 days, P=0.003). Mice in the GEM/AB X2 group had longer OS when compared to mice in the GEM/AB X1 group (33.5 vs 13 days, P=0.001). Mice treated with chemotherapy survived longer than untreated control animals (median OS: 6.5 days, P <0.001). Moreover, in mice treated with chemotherapy, mice that received 2 cycles of GEM/AB X2 had the longest survival, while the FFX X1 group had the poorest OS (P <0.001). Lastly, chemotherapy followed by consolidative SBRT trended towards increased local control and survival.

ConclusionsWe demonstrate the utility and feasibility of clinicall relevant mFOLFIRINOX and gemcitabine/nab-paclitaxel in preclinical models of pancreatic cancer. Patterns of piRNA regulation in Drosophila revealed through transposable element clade inference Transposable elements (TEs) are self-replicating "genetic parasites" ubiquitous to eukaryotic genomes. In addition to conflict between TEs and their host genomes, TEs of the same family are in competition with each other. They compete for the same genomic niches while experiencing the same regime of copy-number selection. This suggests that competition among TEs may favor the emergence of new variants that can outcompete their ancestral forms. To investigate the sequence evolution of TEs, we developed a method to infer clades: collections of TEs that share SNP variants and represent distinct TE family lineages. We applied this method to a panel of 85 Drosophila melanogaster genomes and found that the genetic variation of several TE families shows significant population structure that arises from the population- specific expansions of single clades. We used population genetic theory to classify these clades into younger versus older clades and found that younger clades are associated with a greater abundance of sense and antisense piRNAs per copy than older ones. Further, we find that the abundance of younger, but not older clades, is positively correlated with antisense piRNA production, suggesting a general pattern where hosts preferentially produce antisense piRNAs from recently active TE variants. Together these findings suggest a pattern whereby new TE variants arise by mutation and then increase in copy number, followed by the host producing antisense piRNAs that may be used to silence these emerging variants. Estimating the fundamental niche: accounting for the uneven availability of existing climates Studies that question important conceptual and methodological aspects of the field of ecological niche modeling (and species distribution modeling) have cast doubts on whether it is possible to estimate the fundamental niche of a species using presence-only data. The main limitation in niche estimation is that presence data come from the realized niche, which is only a subset of the fundamental niche. Most of the existing methods lack the ability to overcome this limitation and therefore fit niches that resemble the realized niche. To obtain a more accurate estimate of the fundamental niche, we propose using the geographic region that is accessible to a species (based on its dispersal ability) to determine a sampling distribution, in environmental space, from which we can quantify the likelihood of observing a particular environmental combination in a sample of presence points. We incorporate this sampling distribution into a multivariate normal model (i.e., a Mahalanobis distance model) by creating a weight function that modifies the probabilities of observing different environmental combinations in a sample of presences. This modification accounts for the uneven availability of environmental conditions. We show that the parameters of this weighted-normal model can be approximated with a maximum likelihood estimation approach; and then used to draw ellipsoids (confidence regions) that represent the fundamental niche of the species. We illustrate the application of our model with two worked examples. First, we use presence data for an invasive species and an accessible area that includes only its native range to evaluate whether the fitted model predicts confirmed establishments of the species outside its native range. Second, we use presence data for closely related species with known accessible areas to demonstrate how the different dispersal abilities of the species constrain a classic Mahalanobis distance model. Overall, we show that accounting for the distribution of environmental conditions that are accessible to a species indeed affects the estimation of the ellipsoids used to model its fundamental niche. Multispecific DARPin(R) therapeutics demonstrate very high potency against SARS-CoV-2 variants in vitro SARS-CoV-2 has infected millions of people globally and continues to undergo evolution. Emerging variants can be partially resistant to vaccine induced immunity and therapeutic antibodies, emphasizing the urgent need for accessible, broad-spectrum therapeutics. Here, we report a comprehensive study of ensovibep, the first trispecific clinical DARPin candidate, that can simultaneously engage all three units of the spike protein trimer to potently inhibit ACE2 interaction, as revealed by structural analyses. The cooperative binding of the individual modules enables ensovibep to retain inhibitory potency against all frequent SARS-CoV-2 variants, including Omicron BA.1 and BA.2, as of February 2022. Moreover, viral passaging experiments show that ensovibep, when used as a single agent, can prevent development of escape mutations comparably to a cocktail of monoclonal antibodies (mAb). Finally, we demonstrate that the very high in vitro antiviral potency also translates into significant therapeutic protection and reduction of pathogenesis in Roborovski dwarf hamsters infected with either the SARS-CoV-2 wild-type or the Alpha variant. In this model, ensovibep prevents fatality and provides substantial protection equivalent to the standard of care mAb cocktail. These results support further clinical evaluation and indicate that ensovibep could be a valuable alternative to mAb cocktails and other treatments for COVID-19. 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 thermodynamic model for water activity and redox potential in evolution and development Reactions involving water and oxygen are basic features of geological and biological processes. To understand how life interacts with its environment requires monitoring interactions with H2O and O2 at timescales relevant to not only organismal growth but also billions of years of geobiological evolution. Here, chemical transformations intrinsic to evolution and development were characterized by analyzing data from recent phylostratigraphic and proteomic studies. This two-stage analysis involves obtaining chemical metrics (carbon oxidation state and stoichiometric hydration state) from the elemental compositions of proteins followed by modeling the relative stabilities of target proteins against a proteomic background to infer thermodynamic parameters (oxygen fugacity, water activity, and virtual redox potential (Eh)). The main results of this study are a rise in carbon oxidation state of proteins spanning the time of the Great Oxidation Event, a rise in virtual redox potential that coincides with the likely emergence of aerobic metabolism, a rise in carbon oxidation state of proteins inferred from the transcriptome in late stages of Bacillus subtilis biofilm growth, and a drop in stoichiometric hydration state of the fruit fly developmental proteome at the same time as a drop in organismal water content. Stoichiometric hydration state also decreases for proteins with more recent gene ages and through stages of biofilm development, leading to predictions of higher hydration potentials at earlier time points. By building chemical representations of protein evolution and developmental proteomes, exciting new types of geobiochemical proxies can be developed.