Microglial refinement of A-fibre projections in the postnatal spinal cord dorsal horn is required for normal maturation of dynamic touch. Sensory systems are shaped in postnatal life by the refinement of synaptic connections. In the dorsal horn of the spinal cord, sensory circuits undergo postnatal activity dependent reorganisation, including the retraction of primary afferent A-fibres from superficial to deeper laminae which is accompanied by decreases in cutaneous sensitivity. Here we show that microglia, the resident immune cells in the CNS, phagocytose A-fibre terminals in superficial laminae in the first weeks of life. Genetic perturbation of microglial engulfment at that time prevents the normal process of A-fibre retraction, resulting in increased sensitivity of dorsal horn cells to dynamic tactile cutaneous stimulation, and behavioural hypersensitivity to dynamic touch. Thus, functional microglia are necessary for normal postnatal development of dorsal horn sensory circuits. In the absence of microglial engulfment, superfluous A-fibre projections remain in the dorsal horn and the balance of sensory connectivity is disrupted, leading to lifelong hypersensitivity to dynamic touch.

Impact statementMicroglia phagocytose superfluous A-fibres in the superficial spinal dorsal horn during normal development, the disruption of which leads to long term aberrant dynamic touch processing and behaviour.
Emergence of task information from dynamic network interactions in the human brain How cognitive task behavior is generated by brain network interactions is a central question in neuroscience. We hypothesized that cognitive task behavior emerges from spatiotemporal activity patterns flowing from cognitive control networks (CCNs) - large-scale brain networks that represent task goals - to the motor network. We tested this hypothesis by first localizing the spatiotemporal signatures of task information in the human brain using MRI-individualized source electroencephalography and dynamic decoding. We then applied a novel network modeling approach to simulate the flow of task-evoked activity over causally interpretable resting-state functional connections (dynamic, lagged, direct and directional) to accurately predict empirical response information dynamics underlying behavior. Simulated network lesions confirmed our hypothesis that CCNs are the primary causal drivers of response information, transitioning from early dorsal attention network-dominated sensory-to-response transformation to later collaborative CCN engagement during task response selection. These results reveal the dynamic network interactions by which CCNs produce cognitive task behavior.
Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. Naming Human Diseases: Ethical Principles of Curating Exclusive Substitute for Inopportune Nosology BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. Circulating SARS-CoV-2 variants B.1.1.7, 501Y.V2, and P.1 have gained ability to utilize rat and mouse Ace2 and altered in vitro sensitivity to neutralizing antibodies and ACE2-Ig The emergence of SARS-CoV-2 variants poses greater challenges to the control of COVID-19 pandemic. Here, we parallelly investigated three important characteristics of seven SARS-CoV-2 variants, including two mink-associated variants, the B.1.617.1 variant, and the four WHO-designated variants of concerns (B.1.1.7, B.1.351, P.1, and B.1.617.2). We first investigated the ability of these variants to bind and use animal ACE2 orthologs as entry receptor. We found that, in contrast to a prototype variant, the B.1.1.7, B.1.351, and P.1 variants had significantly enhanced affinities to cattle, pig, and mouse ACE2 proteins, suggesting increased susceptibility of these species to these SARS-CoV-2 variants. We then evaluated in vitro neutralization sensitivity of these variants to four monoclonal antibodies in clinical use. We observed that all the variants were partially or completely resistant against at least one of the four tested antibodies, with B.1.351 and P.1 showing significant resistance to three of them. As ACE2-Ig is a broad-spectrum anti-SARS-CoV-2 drug candidate, we then evaluated in vitro neutralization sensitivity of these variants to eight ACE2-Ig constructs previously described in three different studies. All the SARS-CoV-2 variants were efficiently neutralized by these ACE2-Ig constructs. Interestingly, compared to the prototype variant, most tested variants including the variants of concern B.1.1.7, B.1.351, P.1, and B.1.617.2 showed significantly increased (up to [~]15-fold) neutralization sensitivity to ACE2-Ig constructs that are not heavily mutated in the spike-binding interface of the soluble ACE2 domain, suggesting that SARS-CoV-2 evolves toward better utilizing ACE2, and that ACE2-Ig is an attractive drug candidate for coping with SARS-CoV-2 mutations. bulkAnalyseR: An accessible, interactive pipeline for analysing and sharing bulk multi-modal sequencing data Bulk sequencing experiments (single- and multi-omics) are essential for exploring wide-ranging biological questions. To facilitate interactive, exploratory tasks, coupled with the sharing of easily accessible information, we present bulkAnalyseR, a package integrating state-of-the-art approaches using an expression matrix as the starting point (pre-processing functions are available as part of the package). Static summary images are replaced with interactive panels illustrating quality-checking, differential expression analysis (with noise detection) and biological interpretation (enrichment analyses, identification expression patterns, followed by inference and comparison of regulatory interactions). bulkAnalyseR can handle different modalities, facilitating robust integration and comparison of cis-, trans- and customised regulatory networks.

bulkAnalyseR is available on CRAN and GitHub, with extensive documentation and usage examples (https://github.com/Core-Bioinformatics/bulkAnalyseR, https://cran.r-project.org/web/packages/bulkAnalyseR/)

ContactIrina Mohorianu iim22@cam.ac.uk Distributed brain co-processor for tracking electrophysiology and behavior during electrical brain stimulation Early implantable epilepsy therapy devices provided open-loop electrical stimulation without brain sensing, computing, or an interface for synchronized behavioral inputs from patients. Recent epilepsy stimulation devices provide brain sensing but have not yet developed analytics for accurately tracking and quantifying behavior and seizures. Here we describe a distributed brain co-processor providing an intuitive bi-directional interface between patient, implanted neural stimulation and sensing device, and local and distributed computing resources. Automated analysis of continuous streaming electrophysiology is synchronized with patient reports using a hand-held device and integrated with distributed cloud computing resources for quantifying seizures, interictal epileptiform spikes, and patient symptoms during therapeutic electrical brain stimulation. The classification algorithms for interictal epileptiform spikes and seizures were developed and parameterized using long-term ambulatory data from 9 humans and 8 canines with epilepsy, and then implemented prospectively in out-of-sample testing in 2 pet canines and 4 humans with drug resistant epilepsy living in their natural environments. Phylogenomic analyses of the East Asian endemic Abelia (Caprifoliaceae) shed insights into the temporal and spatial diversification history with widespread hybridizations Background and AimsAbelia (Caprifoliaceae) is a small genus with five species, including one man-made hybrid and several natural hybrids. The genus has a discontinuous distribution in mainland China, Taiwan Island, and the Ryukyu islands, providing a model system to explore mechanisms of species dispersal in the East Asian flora. However, the current phylogenetic relationships within Abelia remain uncertain.

MethodsIn this study, we reconstructed phylogenetic relationships within Abelia using nuclear loci generated by target enrichment and plastomes from genome skimming. Divergence time estimation, ancestral area reconstruction, and ecological niche modelling (ENM) were used to examine the diversification history of Abelia.

Key ResultsWe found extensive cytonuclear discordance across the genus. By integrating lines of evidence from molecular phylogenies, divergence times, and morphology, we propose to merge A. macrotera var. zabelioides into A. uniflora. Network analyses suggested that there have been widespread and multiple hybridization events among Abelia species. These hybridization events may have contributed to the speciation mechanism and resulted in a high observed morphological diversity. The diversification of Abelia began in the early Eocene, followed by A. chinensis var. ionandra colonizing the island of Taiwan in the Middle Miocene. The ENM results suggested an expansion of climatically suitable areas during the Last Glacial Maximum and range contraction during the Last Interglacial. Disjunction between the Himalayan-Hengduan Mountain region (HHM) and the island of Taiwan is most likely the consequence of topographic isolation and postglacial contraction.

ConclusionsWe used genomic data to reconstruct the phylogeny of Abelia and found a clear pattern of reticulate evolution in the group. In addition, our results support shrinkage of postglacial range and the heterogeneity of the terrain have led to the disjunction of the mainland China-Taiwan island. This study provides important new insights into the speciation process and taxonomy of Abelia. Placental uptake and metabolism of 25(OH)Vitamin D determines its activity within the fetoplacental unit Pregnancy 25-hydroxyvitamin D (25(OH)D) concentrations are associated with maternal and fetal health outcomes. Using physiological human placental perfusion and villous explants, we investigate the role of the placenta in regulating the relationships between maternal 25(OH)D and fetal physiology. We demonstrate active placental uptake of 25(OH)D3 by endocytosis, placental metabolism of 25(OH)D3 into 24,25-dihydroxyvitamin D3 and active 1,25-dihydroxyvitamin D [1,25(OH)2D3], with subsequent release of these metabolites into both the maternal and fetal circulations. Active placental transport of 25(OH)D3 and synthesis of 1,25(OH)2D3 demonstrate that fetal supply is dependent on placental function rather than simply the availability of maternal 25(OH)D3. We demonstrate that 25(OH)D3 exposure induces rapid effects on the placental transcriptome and proteome. These map to multiple pathways central to placental function and thereby fetal development, independent of vitamin D transfer. Our data suggest that the underlying epigenetic landscape helps dictate the transcriptional response to vitamin D treatment. This is the first quantitative study demonstrating vitamin D transfer and metabolism by the human placenta, with widespread effects on the placenta itself. These data demonstrate a complex interplay between vitamin D and the placenta and will inform future interventions using vitamin D to support fetal development and maternal adaptations to pregnancy. Development and Characterization of PAR-Trackers: New Tools for Detecting Poly(ADP-ribose) In Vitro and In Vivo AO_SCPLOWBSTRACTC_SCPLOWADP-ribosylation (ADPRylation) is a reversible post-translation modification resulting in the covalent attachment of ADP-ribose (ADPR) moieties on substrate proteins. Naturally-occurring protein motifs and domains, including WWEs, PBZs, and macrodomains, act as "readers" for protein-linked ADPR. Although recombinant, antibody-like ADPR detection reagents containing these readers have facilitated the detection of ADPR, they are limited in their ability to capture the dynamic nature of ADPRylation. Herein, we describe and characterize a set of poly(ADP-ribose) (PAR) Trackers (PAR-Ts) - optimized dimerization-dependent or split-protein reassembly PAR sensors in which a naturally occurring PAR binding domain, WWE, was fused to both halves of dimerization-dependent GFP (ddGFP) or split Nano Luciferase (NanoLuc), respectively. We demonstrate that these new tools allow the detection and quantification of PAR levels in extracts, living cells, and living tissues with greater sensitivity, as well as temporal and spatial precision. Importantly, these sensors detect changes in cellular ADPR levels in response to physiological cues (e.g., hormone-dependent induction of adipogenesis without DNA damage), as well as xenograft tumor tissues in living mice. Our results indicate that PAR Trackers have broad utility for detecting ADPR in many different experimental and biological systems. DNA Strand Breaks and Gaps Target Retroviral Binding and Integration Integration into a host genome is essential for retrovirus infection and is catalyzed by a nucleoprotein complex (Intasome) containing the viral integrase (IN) and reverse transcribed (RT) copy DNA (cDNA). Previous studies demonstrated DNA site recognition limited intasome integration. Using single molecule Forster resonance energy transfer (smFRET), we show Prototype Foamy Virus (PFV) intasomes pause at DNA strand breaks and gaps. The break/gap discontinuities are similar to base excision repair (BER) lesion-processing intermediates, which affect retrovirus integration in vivo. Pausing targeted site-directed integration at the break/gap without inducing intasome conformational alterations. An 8-oxo-guanine lesion normally processes by BER and a G/T mismatch or a +T nucleotide insertion that induce flexibility or a bend in the DNA backbone did not promote intasome pausing or targeted integration. These results suggest that repair intermediates can modulate dynamic intasome-DNA interactions which target retroviral integration. Sarcoplasmic reticulum calcium handling in young DBA/2J vs C57BL/10 mdx models of Duchenne muscular dystrophy The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. HuR-dependent SOD2 protein synthesis is an early adaptation to anchorage independence During metastasis, cancer cells must adapt to survive a loss of anchorage and evade anoikis. An important pro-survival adaptation is the ability of metastatic tumor cells to increase their antioxidant capacity and restore cellular redox balance. Although much is known about the transcriptional regulation of antioxidant enzymes in response to stress, how cells rapidly adapt to alter antioxidant enzyme levels is less well understood. Using ovarian cancer cells as a model, we demonstrate that an increase in protein expression of the mitochondrial superoxide dismutase SOD2 is a very early event initiated in response to cellular detachment. SOD2 protein synthesis is rapidly induced within 0.5-2 hours of matrix detachment, and polyribosome profiling demonstrates an increase in the number of ribosomes bound to SOD2 mRNA, indicating an increase in SOD2 translation in response to anchorage-independence. Mechanistically, we find that anchorage-independence specifically induces cytosolic accumulation of the RNA binding protein HuR/ELAVL1 and leads to increased HuR binding to SOD2 mRNA. Using HuR siRNA-mediated knock-down, we show that the presence of HuR is necessary for the increase in SOD2 mRNA association with the heavy polyribosome fraction and SOD2 protein synthesis observed in anchorage-independence. Cellular detachment activates the stress-response protein kinase p38 MAPK, which is necessary for HuR-SOD2 mRNA binding and optimal increases in SOD2 protein expression. These findings illustrate a novel post-transcriptional regulatory mechanisms of SOD2, enabling cells to rapidly increase their mitochondrial antioxidant capacity as an acute response to anchorage-independence. Identifying and correcting for misspecifications in GWAS summary statistics and polygenic scores Publicly available genome-wide association studies (GWAS) summary statistics exhibit uneven quality, which can impact the validity of follow-up analyses. First, we present an overview of possible misspecifications that come with GWAS summary statistics. Then, in both simulations and real data analyses, we show that additional information such as imputation INFO scores, allele frequencies, and per-variant sample sizes in GWAS summary statistics can be used to detect possible issues and correct for misspecifications in the GWAS summary statistics. One important motivation for us is to improve the predictive performance of polygenic scores built from these summary statistics. Unfortunately, due to the lack of reporting standards for GWAS summary statistics, this additional information is not systematically reported. We also show that using well-matched LD references can improve model fit and translate into more accurate prediction. Finally, we discuss how to make polygenic score methods such as lassosum and LDpred2 more robust to these misspecifications to improve their predictive power. Presence of a SARS-COV-2 protein enhances Amyloid Formation of Serum Amyloid A A marker for the severeness and disease progress of COVID-19 is overexpression of serum amyloid A (SAA) to levels that in other diseases are associated with a risk for SAA amyloidosis. In order to understand whether SAA amyloidosis could also be a long-term risk of SARS-COV-2 infections we have used long all-atom molecular dynamic simulations to study the effect of a SARS-COV-2 protein segment on SAA amyloid formation. Sampling over 40 {micro}s we find that presence of the nine-residue segment SK9, located at the C-terminus of the Envelope protein, increases the propensity for SAA fibril formation by three mechanisms: it reduces the stability of the lipid-transporting hexamer shifting the equilibrium toward monomers, it increases the frequency of aggregation-prone configurations in the resulting chains, and it raises the stability of SAA fibrils. Our results therefore suggest that SAA amyloidosis and related pathologies may be a long-term risk of SARS-COV-2 infections. 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. Actomyosin activity-dependent apical targeting of Rab11 vesicles reinforces apical constriction During tissue morphogenesis, cell shape changes resulting from cell-generated forces often require active regulation of intracellular trafficking. How mechanical stimuli influence intracellular trafficking and how such regulation impacts tissue mechanics are not fully understood. In this study, we identify an actomyosin dependent mechanism involving Rab11- mediated trafficking in regulating apical constriction in the Drosophila embryo. During Drosophila mesoderm invagination, apical actin and Myosin II (actomyosin) contractility induces apical accumulation of Rab11-marked vesicle-like structures ("Rab11 vesicles") by promoting a directional bias in dynein mediated vesicle transport. At the apical domain, Rab11 vesicles are enriched near the adherens junctions (AJs). The apical accumulation of Rab11 vesicles is essential to prevent fragmented apical AJs, breaks in the supracellular actomyosin network and a reduction in the apical constriction rate. This Rab11 function is separate from its role in promoting apical Myosin II accumulation. These findings suggest a feedback mechanism between actomyosin activity and Rab11-mediated intracellular trafficking that regulates the force generation machinery during tissue folding. Poly ADP-ribosylation of SET8 leads to aberrant H4K20 methylation domains in mammalian cells In mammalian cells, SET8 mediated Histone H4 Lys 20 monomethylation (H4K20me1) has been implicated in regulating mitotic condensation, DNA replication, DNA damage response, and gene expression. Here we show SET8, the only known enzyme for H4K20me1 is post-translationally poly ADP-ribosylated by PARP1 on lysine residues. PARP1 interacts with SET8 in a cell cycle-dependent manner. Poly ADP-ribosylation on SET8 renders it catalytically compromised, and degradation via ubiquitylation pathway. Knockdown of PARP1 led to an increase of SET8 protein levels, leading to aberrant H4K20me1 and H4K20me3 domains genome-wide. H4K20me1 is associated with higher gene transcription levels while the increase of H4K20me3 levels was predominant in DNA repeat elements correlating with modulation of chromatin states. Hence, SET8 mediated chromatin remodeling in mammalian cells are modulated by poly ADP-ribosylation by PARP1. 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. Implicit Laplacian of Enhanced Edge: An Unguided Algorithm for Accurate and Automated Quantitative Analysis of Cytoskeletal Images The eukaryotic cytoskeleton plays essential roles in cell signaling and trafficking, which is broadly associated with immunity and diseases of human and plants. To date, most analyses aiming at defining the temporal and spatial dynamics of the cytoskeleton have relied on qualitative and quantitative analysis of fluorescence images to describe cytoskeletal function. While state-of-the-art, these approaches have limitations: the diverse shape and brightness of the cytoskeleton cause considerable inaccuracy in both human-driven and automated approaches, and the widely adopted image projection process (3D to 2D) leads to bias and information loss. Here, we describe the development and application of Implicit Laplacian of Enhanced Edge (ILEE), an unguided approach that uses a 2D/3D-compatible local thresholding algorithm for the quantitative evaluation of cytoskeletal status and organization at high performance. Using ILEE, we constructed a Python library to enable automated cytoskeletal image analysis, providing numerous biologically-interpretable indices measuring the density, bundling, severing, branching, and directionality of the cytoskeleton. The data presented herein demonstrate that ILEE resolves the defects of classic cytoskeleton analysis approaches, enables the measurement of novel cytoskeletal features, and yields quantitatively descriptive data with superior accuracy, stability, and robustness. We released the ILEE algorithm as an open-source library and further developed a Google Colab interface as a community resource. Predictive uncertainty in mechanistic models of cellular processes calibrated to experimental data Mathematical models are often used to explore network-driven cellular processes from a systems perspective. However, a dearth of quantitative data suitable for model calibration leads to models with parameter unidentifiability and questionable predictive power. Here we introduce a Bayesian and Machine-Learning based Measurement Model approach to explore how quantitative and non-quantitative data constrain models of apoptosis execution within a missing data context. We find two orders of magnitude more ordinal (e.g. immunoblot) data are necessary to achieve accuracy comparable to quantitative (e.g. fluorescence) data. Notably, ordinal and nominal (e.g. immunostain) non-quantitative data synergize to reduce model uncertainty and improve accuracy. Further, model prediction accuracy and certainty strongly depend on rigorous data-driven formulations of the measurement, and the size and make-up of the datasets. Finally, we demonstrate the potential of a data-driven Measurement Model approach to identify model features that could lead to informative experimental measurements and improve model predictive power. Modeling miRNA-driven post-transcriptional regulatory signals by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Ensovibep, a novel trispecific DARPin candidate that protects against SARS-CoV-2 variants 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. An adaptive teosinte mexicana introgression modulates phosphatidylcholine levels and is associated with maize flowering time Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp.parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability, and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern US, Canada and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time. Outsourced hearing in an orb-weaving spider that uses its web as an auditory sensor Hearing is a fundamental sense of many animals, including all mammals, birds, some reptiles, amphibians, fish, and arthropods. The auditory organs of these animals are extremely diverse in anatomy after hundreds of millions of years of evolution, yet all are made up of cellular tissues and are morphologically part of bodies of animals. Here we show hearing in the orb-weaving spider, Larinioides sclopetarius is not constrained by the organisms body but is extended through outsourcing hearing to its extended phenotype, the proteinaceous, self-manufactured orb-web. We find the wispy, wheel-shaped orb-web acts as a hyperacute acoustic "antenna" to capture the sound-induced air particle movements that approach the maximum physical efficiency, better than the acoustic responsivity of all previously known eardrums. By sensing the motion of web threads, the spider remotely detects and localizes the source of an incoming airborne acoustic wave such as those emitted by approaching prey or predators. By outsourcing its acoustic sensors to its web, the spider is released from body size constraints and permits the araneid spider to increase its sound-sensitive surface area enormously, up to 10,000 times greater than the spider itself. The spider also enables the flexibility to functionally adjust and regularly regenerate its "external eardrum" according to its needs. The "outsourcing" and "supersizing" of auditory function in spiders provides unique features for studying extended and regenerative sensing, and designing novel acoustic flow detectors for precise fluid dynamic measurement and manipulation. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development.