Improved radiation expression profiling in blood by sequential application of sensitive and specific gene signatures PurposeCombinations of expressed genes can discriminate radiation-exposed from normal control blood samples by machine learning based signatures (with 8 to 20% misclassification rates). These signatures can quantify therapeutically-relevant as well as accidental radiation exposures. The prodromal symptoms of Acute Radiation Syndrome (ARS) overlap those present in Influenza and Dengue Fever infections. Surprisingly, these human radiation signatures misclassified gene expression profiles of virally infected samples as false positive exposures. The present study investigates these and other confounders, and then mitigates their impact on signature accuracy.

MethodsThis study investigated recall by previous and novel radiation signatures independently derived from multiple Gene Expression Omnibus datasets on common and rare non-malignant blood disorders and blood-borne infections (thromboembolism, S. aureus bacteremia, malaria, sickle cell disease, polycythemia vera, and aplastic anemia). Normalized expression levels of signature genes are used as input to machine learning-based classifiers to predict radiation exposure in other hematological conditions.

ResultsExcept for aplastic anemia, these blood-borne disorders modify the normal baseline expression values of genes present in radiation signatures, leading to false-positive misclassification of radiation exposures in 8 to 54% of individuals. Shared changes, predominantly in DNA damage response and apoptosis-related gene transcripts in radiation and confounding hematological conditions, compromise the utility of these signatures for radiation assessment. These confounding conditions (sickle cell disease, thromboembolism, S. aureus bacteremia, malaria) induce neutrophil extracellular traps, initiated by chromatin decondensation, DNA damage response and fragmentation followed by programmed cell death. Riboviral infections (for example, Influenza or Dengue fever) have been proposed to bind and deplete host RNA binding proteins, inducing R-loops in chromatin. R-loops that collide with incoming replication forks can result in incompletely repaired DNA damage, inducing apoptosis and releasing mature virus. To mitigate the effects of confounders, we evaluated predicted radiation-positive samples with novel gene expression signatures derived from radiation-responsive transcripts encoding secreted blood plasma proteins whose expression levels are unperturbed by these conditions.

ConclusionsThis approach identifies and eliminates misclassified samples with underlying hematological or infectious conditions, leaving only samples with true radiation exposures. Diagnostic accuracy is significantly improved by selecting genes that maximize both sensitivity and specificity in the appropriate tissue using combinations of the best signatures for each of these classes of signatures.
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.
Long term analysis of social structure: evidence of age-based consistent associations in male Alpine ibex Despite its recognized importance for understanding the evolution of animal sociality as well as for conservation, long term analysis of social networks of animal populations is still relatively uncommon. We investigated social network dynamics in males of a gregarious mountain ungulate (Alpine ibex, Capra ibex) over ten years focusing on groups, sub-groups and individuals, exploring the dynamics of sociality over different scales. Despite the social structure changing between seasons, the Alpine ibex population was highly cohesive: fission-fusion dynamics lead almost every male in the population to associate with each other male at least once. Nevertheless, we found that male Alpine ibex showed preferential associations that were maintained across seasons and years. Age seemed to be the most important factor driving preferential associations while other characteristics, such as social status, appeared less crucial. We also found that centrality measures were influenced by age and were also related to individual physical condition. The multi-scale and long-term frame of our study helped us show that ecological constrains, such as resource availability, may play a role in shaping associations in a gregarious species, but they cannot solely explain sociality and preferential association that are likely also to be driven by life-history linked physiological and social needs. Our results highlight the importance of long-term studies based on individually recognizable subjects to help us build on our understanding of the evolution of animal sociality. TAP: Targeting and analysis pipeline for optimization and verification of coil placement in transcranial magnetic stimulation ObjectiveTranscranial magnetic stimulation (TMS) can modulate brain function via an electric field (E-field) induced in a region of interest (ROI) in the brain. The ROI E-field can be computationally maximized and set to match a specific reference using individualized head models to find the optimal coil placement and stimulus intensity. However, the available software lacks many practical features for prospective planning of TMS interventions and retrospective evaluation of the experimental targeting accuracy.

ApproachThe TMS targeting and analysis pipeline (TAP) software uses an MRI/fMRI-derived brain target to optimize coil placement considering experimental parameters such as the subjects hair thickness and coil placement restrictions. The coil placement optimization is implemented in SimNIBS 3.2, for which an additional graphical user interface (TargetingNavigator) is provided to visualize and adjust procedural parameters. The coil optimization process also computes the E-field at the target, allowing the selection of the TMS device intensity setting to achieve a specific E-field strength. The optimized coil placement information is prepared for neuronavigation software, which supports targeting during the TMS procedure. The neuronavigation system can record the coil placement during the experiment, and these data can be processed in TAP to quantify the accuracy of the experimental TMS coil placement and induced E-field.

Main resultsTAP was demonstrated in a study consisting of three repetitive TMS sessions in five subjects. TMS was delivered by an experienced operator under neuronavigation with the computationally optimized coil placement. Analysis of the experimental accuracy from the recorded neuronavigation data indicated coil location and orientation deviations up to about 2 mm and 2{degrees}, respectively, resulting in an 8% median decrease in the target E-field magnitude compared to the optimal placement.

SignificanceTAP supports navigated TMS with a variety of features for rigorous and reproducible stimulation delivery, including planning and evaluation of coil placement and intensity selection for E-field-based dosing. Frog nest foams exhibit pharmaceutical foam-like properties Foams have frequently been used as systems for the delivery of cosmetic and therapeutic molecules; however, there is high variability in the foamability and long-term stability of synthetic foams. The development of pharmaceutical foams that exhibit desirable foaming properties, delivering appropriate amounts of the active pharmaceutical ingredient (API) and that have excellent biocompatibility is of great interest. The production of stable foams is rare in the natural world; however, certain species of frogs have adopted foam production as a means of providing a protective environment for their eggs and larvae from a predators and parasites, to prevent desiccation, to control gaseous exchange, temperature extremes, and to reduce UV damage. These foams show great stability (up to 10 days in tropical environments) and are highly biocompatible due to the sensitive nature of amphibian skin. This work demonstrates for the first time, that nests of the Tungara frog (Engystomops pustulosus) is stable ex situ with useful physiochemical and biocompatible properties and is capable of encapsulating a range of compounds, including antibiotics. These protein foam mixtures share some properties with pharmaceutical foams and may find utility in a range of pharmaceutical applications such as a topical drug delivery systems (DDS). Design, Mutate, Screen: High-throughput creation of genetic clocks with different period-amplitude characteristics Directed evolution has become an invaluable tool in protein engineering and has also greatly influenced the construction of synthetic gene circuits. The ability to generate diversity at precise targets for directed evolution approaches has improved vastly, allowing researchers to create large, specific mutant libraries with relative ease. Screening approaches for large mutant libraries have similarly come a long way, especially when the desired behavior can easily be tested for with static, single time-point assays. For more complex gene circuits with dynamic phenotypes that change over time, directed evolution approaches to controlling and tuning circuit behavior have been hindered by the lack of sufficiently high-throughput screening methods to isolate variants with desired characteristics. Here we utilize directed mutagenesis and multiplexed microfluidics to develop a workflow for creating, screening and tuning dynamic gene circuits that operate at the population level. Specifically, we create a mutant library of an existing oscillator, the synchronized lysis circuit, and tune its dynamics while uncovering principles regarding its behavior. Lastly, we utilize this directed evolution workflow to construct a new synchronized genetic oscillator that exhibits robust dynamics over long time scales. Unifying Single-Channel Permeability from Rare-Event Sampling and Steady-State Flux Various all-atom molecular dynamics (MD) simulation methods have been developed to compute free energies and crossing rates of ions and small molecules through ion channels. However, a systemic comparison across different methods is scarce. Using a carbon nanotube as a model of small conductance ion channel, we computed the single-channel permeability for potassium ion using umbrella sampling, Markovian milestoning, and steady-state flux under applied voltage. We show that a slightly modified inhomogeneous solubility-diffusion equation yields a single-channel permeability consistent with the mean first passage time (MFPT)-based method. For milestoning, applying cylindrical and spherical bulk boundary conditions yield consistent MFPT if factoring in the effective bulk concentration. The sensitivity of the MFPT to the output frequency of collective variables is highlighted using the convergence and symmetricity of the inward and outward MFPT profiles. The consistent transport kinetic results from all three methods demonstrated the robustness of MD-based methods in computing ion channel permeation. The advantages and disadvantages of each technique are discussed, focusing on the future applications of milestoning in more complex systems. Comparative miRNA Transcriptomics of Mouse and Macaque Reveals MYOC is An Inhibitor for C. neoformans Invasion into Brain Cryptococcal meningoencephalitis is an emerging infection shifted from primarily ART-naive to being ART-experienced HIV/AIDS patients, COVID-19 patients and also in immune competent individuals, mainly caused by the human opportunistic pathogen Cryptococcus neoformans, yet mechanisms of the brain or CNS dissemination remain to elucidate, which is the deadest process for the disease. Meanwhile, illustrations of clinically relevant responses in cryptococcosis were limited, as the low availabilities of clinical samples. In this study, macaque and mouse infection models were employed and miRNA-mRNA transcriptomes were performed and combined, which revealed cytoskeleton, a major feather in HIV/AIDS patients, was a centric pathway regulated in both two infection models. Notably, assays of clinical immune cells confirmed an enhanced "Trojan Horse" in HIV/AIDS patients, which can be shut down by cytoskeleton inhibitors. Furthermore, we identified a novel enhancer for macrophage "Trojan Horse", myocilin, and an enhanced fungal burden was achieved in brains of MYOC transgenic mice. Taking together, this study reveals fundamental roles of cytoskeleton and MYOC in blocking fungal CNS dissemination, which not only helps to understand the high prevalence of cryptococcal meningitis in HIV/AIDS, but also facilitates the development of novel drugs for therapies of meningoencephalitis caused by C. neoformans and other pathogenic microorganisms. High glucose enhances antigen-independent CTL killing via TRAIL Cytotoxic T lymphocytes (CTLs) are involved in development of diabetes. However, the impact of excessive glucose on CTL-mediated antigen-independent killing remains elusive. Here, we report that TNF-related apoptosis inducing ligand (TRAIL) is substantially up- regulated in CTLs in environments with high glucose (HG) both in vitro and in vivo. The PI3K- Akt-NF{kappa}B axis and non-mitochondrial reactive oxygen species are essential in HG-induced TRAIL upregulation in CTLs. TRAILhigh CTLs induce apoptosis of pancreatic beta cell line 1.4E7. Metformin and Vitamin D synergistically reduce HG-enhanced expression of TRAIL in CTLs and coherently protect 1.4E7 cells from TRAIL-mediated apoptosis. Notably, in patients with diabetes, correlation between Vitamin D concentrations in plasma and glucose levels is linked to HG-enhanced TRAIL expression on CTLs. Microarray data reveal that OXCT2, an important enzyme in ketone body catabolism, is a promising target in response to vitamin D. Our work not only reveals a novel mechanism of CTL involvement in progression of diabetes, but also establishes CTLs as a target for combined metformin and vitamin D therapy to protect pancreatic beta cells of diabetic patients. The response network of HSP70 defines vulnerabilities in cancer cells with the inhibited proteasome Human neoplasias are often addicted to the proteasome machinery. However, cancers have evolved efficient response mechanisms to overcome proteasome inhibition with bortezomib and carfilzomib - drugs approved for multiple myeloma treatment. To understand these responses we investigated proteome changes upon the proteasome inhibition with carfilzomib - in multiple myeloma, normal fibroblasts, and cancers of lung, colon, and pancreas. A pathway-oriented siRNA screen based on the proteomics results showed that molecular chaperones, autophagy- and endocytosis-related proteins are cancer-specific vulnerabilities combined with carfilzomib. Targeting of HSPA1A/B (HSP70 family chaperones) most specifically sensitized cancer cells and patient-derived organoids to the proteasome inhibition. A high level of HSPA1A/B mRNA correlated with a low proteasome activity in cancer patient tissues and is a risk factor in cancer patients with a low proteasome expression. Mechanistically, HSPA1A/B governed autophagy, unfolded protein response, endocytic trafficking, and chaperoned the proteasome machinery, suppressing the effect of the proteasome inhibition, but did not control the NRF1/2-driven proteasome subunit transcriptional bounce-back. Consequently, downregulation of NRF1 most specifically decreased the viability of cancer cells with the inhibited proteasome and HSP70. Fatal COVID-19 and non-COVID-19 Acute Respiratory Distress Syndrome is Associated with Incomplete Alveolar Type 1 Epithelial Cell Differentiation from the Transitional State Without Fibrosis ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. AICAR-mediated selection of karyotypically normal cells in Pallister-Killian syndrome patient-derived skin fibroblasts Pallister-Killian Syndrome (PKS) is a mosaic aneuploidy syndrome typically caused by the presence of a supernumerary marker isochromosome containing two extra copies of the short arm of chromosome 12 (iso-chromosome 12p or tetrasomy 12p). This isochromosome is always present in a mosaicism state that has tissue limited distribution. PKS is characterized by severe neurodevelopmental delay, intellectual disability, multisystem involvement and congenital malformations including typical dysmorphic features and skin pigmentation anomalies.

Aneuploid cells, irrespective of the identity of the supernumerary chromosome, including cancer cells, yeast cells and mouse embryonic fibroblasts (MEFs), have been demonstrated to present a disruption of protein homeostasis and increased basal stress levels; resulting in a greater sensitivity to chemical compounds inducing cellular energy stress compared to euploid cell lines. The burden of trisomy 21 has also been recently shown to impair the proteostasis network in lymphoblastoid cell lines and fibroblasts obtained from individuals with Down syndrome.

In this study, we demonstrate that AICAR, 5-Aminoimidazole-4-carboxamide ribonucleotide, a known energy stress inducing drug with antiproliferative effects on aneuploidy cancer cells and MEFs, is also able to selectively eliminate cells carrying the isochromosme12p in PKS clones in a time and dosage dependent manner. Collectively, our results indirectly provide evidence of increased basal energy and proteotoxic stress in PKS cells carrying isochromosome 12p, and suggest a potential therapeutic drug-based strategy that, selectively acting as a stressor for aneuploid cells, may establish the euploid state in PKS and a broader spectrum of human mosaic disorders. Critical re-evaluation of the slope factor of the operational model of agonism Background and PurposeAlthough being a relative term, agonist efficacy is a cornerstone in the proper assessment of agonist selectivity and signalling bias. The operational model of agonism (OMA) has become successful in the determination of agonist efficacies and ranking them. In 1985, Black et al. introduced the slope factor to the OMA to make it more flexible and allow for fitting steep as well as flat concentration-response curves.

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

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

Conclusion and ImplicationsThe OMA modified by the Hill coefficient can be used for correct ranking of agonism and subsequent analysis, like quantification of signalling bias. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. Self-learning algorithm for denoising and advancing the integration of scRNA-seq datasets improves the identification of resilient and susceptible retinal ganglion cell types The clustering of single cell RNA-sequencing (scRNA-seq) data enables the classification of cell types, and the development of integration algorithms has enabled the tracing of altered-from-baseline transcriptomes to their respective cell type origins. Here, we developed an algorithm that removes sources of noise from scRNA-seq datasets and improves the accuracy of their integration. We found that this improvement of integration is further increased by an algorithm we developed for auto-determination of optimal weight-assignment to anchors. As a proof-of-concept, we show that our algorithm determined the type-origin of the 17% of injured retinal ganglion cells (RGCs) that an existing integration algorithm could not identify. As we found that the majority of the originally unassigned cells belonged to only several out of 45 RGC types, a consequent change in the proportions of the surviving types resulted in a moderately different ranking of resiliency to injury. This updated standard for RGC type resilience ranking improved the understanding of RGC type biology and the ranking-dependent prediction of genes for therapeutic neuroprotective targeting. Additional bioinformatic analyses revealed single-gene cluster-markers for RGC types and contributed new insights into the global characteristics of RGC types and how axonal injury affects them, showing how the extent of variance between transcriptomes of RGC types increases differentially during maturation and after injury. We also show, for the first time, that specific global properties of the transcriptome can predict the resilience to injury of certain cell types. Furthermore, we developed a website for cluster-by-cluster comparison of gene expression between uninjured and injured RGC types. The publicly available R library package, CellTools, for denoising and improving the integration of scRNA-seq datasets, as well as for associating cell types with the global properties of the transcriptome, will assist scRNA-seq studies of developing, injured, or diseased cell types across biological fields. Vectofusin-1 based T-cell transduction approach for generating murine chimeric antigen receptor T-cells for cancer. Gene transfer into human and murine T-cells using viral-based approaches has several promising therapeutic applications including the production of chimeric antigen receptor T-cell (CAR-T) therapy. The generation of murine CAR-T is paramount to test and validate immunocompetent mouse models for CAR-T therapy. Several viral transduction enhancers already exist for gene therapy with few limitations. In this study, we tested vectofusin-1, a short cationic peptide, as a soluble transduction enhancer for gammaretroviral transduction for the generation of anti-CD19 murine CAR-T. We found that in comparison to Retronectin, Vectofusin-1 is an equally optimal transduction enhancer for the generation of murine CAR-T cells. Protein scaffold-based multimerization of soluble ACE2 efficiently blocks SARS-CoV-2 infection in vitro and in vivo Soluble ACE2 (sACE2) decoy receptors are promising agents to inhibit SARS-CoV-2, as their efficiency is less likely to be affected by common escape mutations in viral proteins. However, their success may be limited by their relatively poor potency. To address this challenge, we developed a large decoy library of sACE2 fusion proteins, generated with several protease inhibitors or multimerization tags. Among these decoys, multimeric sACE2 consisting of SunTag or MoonTag systems, which were originally utilized for signal amplification or gene activation systems, were extremely effective in neutralizing SARS-CoV-2 in pseudoviral systems and in clinical isolates. These novel sACE2 fusion proteins exhibited greater than 100-fold SARS-CoV-2 neutralization efficiency, compared to monomeric sACE2. SunTag or MoonTag in combination with a more potent version of sACE2, which has multiple point mutations for greater binding (v1), achieved near complete neutralization at a sub-nanomolar range, comparable with clinical monoclonal antibodies. Pseudoviruses bearing mutant versions of Spike, alpha, beta, gamma or delta variants, were also neutralized efficiently with SunTag or MoonTag fused sACE2(v1). Finally, therapeutic treatment of sACE2(v1)-MoonTag provided protection against SARS-CoV-2 infection in an in vivo mouse model. Overall, we suggest that the superior activity of the sACE2-SunTag or sACE2-MoonTag fusions is due to the greater occupancy of the multimeric sACE2 receptors on Spike protein as compared to monomeric sACE2. Therefore, these highly potent multimeric sACE2 decoy receptors may offer a promising treatment approach against SARS-CoV-2 infections.

One Sentence SummaryMultimerization of sACE2 markedly enhanced the neutralization of SARS-CoV-2 by blocking multiple viral spike proteins simultaneously. Text has been revised, image definition has been revised. Maize is a staple food all over the world. Models for climate change suggest that, in the future, cloud formation might be reduced in the tropics increasing the exposure to Ultraviolet-B (UV-B) radiation, a DNA-damaging agent. UV-B (290 to 320 nm) has been shown to affect yield in maize. In this project we have determined the differences in DNA repair efficiencies between U.S. inbred lines B73 and Mo17, and Central American purple landraces from Guatemala and Costa Rica. Our results from single cell electrophoresis experiments (Comet Assay) suggest that the landrace Pujagua Santa Cruz (P1, Costa Rican) was resistant to damage caused by the radiomimetic agent zeocin (24h/100 micrograms per mL), while landrace Pujagua La Cruz (P2, Costa Rican) was able to repair DNA damage after one hour. On the other hand, line Mo17 (Missouri, USA) was unable to repair the damage, while B73 (Iowa, USA) and the landraces Jocopilas (Guatemalan), Orotina Congo (Costa Rican) and Talamanca (Costa Rican) were partially able to repair DNA damage. High Resolution Melting (HRM) curve analysis of putative homologous DNA repair gene ZeaATM1 showed that both P1 and P2 had differences in the melting temperatures for this gene compared to B73 and Mo17, while P1 showed additional differences in ZeaSOG1, ZeaRAD51 and ZeaBRCA1, suggesting that in this landrace the presence of polymorphisms may be common among key genes for this pathway. Taken together our results suggest that key adaptive differences in DNA repair efficiencies exist between inbred lines and landraces of maize and that some Central America landraces could be used as a valuable pool of alleles for plant breeding aiming to increase tolerance to radiation. Sahul's megafauna were vulnerable to plant-community changes due to their position in the trophic network Extinctions stemming from environmental change often trigger trophic cascades and coextinctions. However, it remains unclear whether trophic cascades were a large contributor to the megafauna extinctions that swept across several continents in the Late Pleistocene. The pathways to megafauna extinctions are particularly unclear for Sahul (landmass comprising Australia and New Guinea), where extinctions happened earlier than on other continents. We investigated the role of bottom-up trophic cascades in Late Pleistocene Sahul by constructing pre-extinction (~ 80 ka) trophic network models of the vertebrate community of Naracoorte, south-eastern Australia. These models allowed us to predict vertebrate species vulnerability to cascading extinctions based on their position in the network. We tested whether the observed extinctions could be explained by bottom-up cascades, or if they should be attributed to other external causes. Species that disappeared from the community were more vulnerable, overall, to bottom-up cascades than were species that survived. The position of extinct species in the network - having few or no predators - also suggests they might have been particularly vulnerable to a new predator. These results provide quantitative evidence that trophic cascades and naivety to predators could have contributed to the megafauna extinction event in Sahul. Comprehensive O-GlcNAc glycoproteomics on NOTCH1 EGF repeats implicated unique Lewis X epitopes in mammals The O-GlcNAc modification of Notch receptors regulates Notch ligand interactions in a manner distinct from other forms of O-glycans on epidermal growth factor-like (EGF) repeats of Notch receptors. Although many proteins, besides Notch receptors, are expected to be O-GlcNAcylated by EGF domain-specific O-GlcNAc transferase (EOGT), only a small number of proteins have been reported to be modified in vivo, and elongated O-GlcNAc glycans have not been extensively explored. To extend our view of the specificity and variety of the glycan modification, we conducted a comprehensive analysis of O-GlcNAc glycans on NOTCH1 in mammals. Mass spectrometric analysis of NOTCH1 fragments expressed in HEK293T cells revealed that several EGF domains with putative O-GlcNAcylation sites were hardly modified with O-GlcNAc. Although amino acid residues before the modification site are preferentially occupied with aromatic residues, Phe and Tyr are preferable to Trp for the apparent modification with O-GlcNAc. Furthermore, a minor form of fucosylated O-GlcNAc glycans was detected in a subset of EGF domains. Fucosylation of O-GlcNAc glycans was enhanced by FUT1, FUT2, or FUT9 expression. The FUT9-dependent Lewis X epitope was confirmed by immunoblotting using an anti-Lewis X antibody. As expected from the similarity in the glycan structures, the Lexis X antigen was detected on O-fucose glycans. Our results refined the putative consensus sequence for the EOGT-dependent extracellular O-GlcNAc modification in mammals and revealed the structural diversity of functional Notch O-glycans. Single cell genomics reveals plastid-lacking Picozoa are close relatives of red algae The endosymbiotic origin of plastids from cyanobacteria gave eukaryotes photosynthetic capabilities and launched the diversification of countless forms of algae. These primary plastids are found in members of the eukaryotic supergroup Archaeplastida. All known archaeplastids still retain some form of primary plastids, which are widely assumed to have a single origin. Here, we used single-cell genomics from natural samples combined with phylogenomics to infer the evolutionary origin of the phylum Picozoa, a globally distributed but seemingly rare group of marine microbial heterotrophic eukaryotes. Strikingly, the analysis of 43 single-cell genomes shows that Picozoa belong to Archaeplastida, specifically related to red algae and the phagotrophic rhodelphids. These picozoan genomes support the hypothesis that Picozoa lack a plastid, and further reveal no evidence of an early cryptic endosymbiosis with cyanobacteria. These findings change our understanding of plastid evolution as they either represent the first complete plastid loss in a free-living taxon, or indicate that red algae and rhodelphids obtained their plastids independently of other archaeplastids. 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. High-quality genome assembly of Cinnamomum burmannii (chvar. Borneol) provides insights into the natural borneol biosynthesis Cinnamomum burmannii (chvar. Borneol) is a well-known medicinal and industrial plant cultivated in the Lingnan region of China. It is the key source from organism of natural borneol (D-borneol), one of the precious and widely used Chinese herbal medicines with a variety of medicinal effects. Here, we report a high-quality chromosome-scale genome assembly of C. burmannii (chvar. Borneol) using Pacbio single-molecule sequencing and Hi-C technology. The assembled genome size was 1.14 GB with a scaffold N50 of 94.30 Mb, while 98.77% of the assembled sequences were anchored on 12 pseudochromosomes including 41549 protein-coding genes. Genomic evolution analysis revealed C. burmannii and C. micranthum shared two Lauraceae unique ancestral whole-genome duplication (WGD) events. Likewise, comparative genomic analysis showed strong collinearity between these two species. Besides, the analysis for Long Terminal Repeat Retrotransposons (LTR-RTs) indicated the outbreak of LTR-RTs insertion made a great contribution to the size difference of genomes between C. burmannii and C. micranthum. Furthermore, the candidate genes in pathway associated with natural borneol synthesis were identified on the genome and their differential expressions were analyzed in various biological tissues. We considered that several of genes in Mevalonate (MVA) Methylerythritol Phosphate (MEP) pathways or in downstream pathway have the potential to be the key factors in the biosynthesis of D-borneol. We also constructed the genome database (CAMD; http://www.cinnamomumdatabase.com/) of Cinnamomum species for a better data utilization in the future. All these results will enrich the genomic data of Lauraceae plants and facilitate genetic improvement of this commercially important plant. 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. 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. TRIM28-dependent SUMOylation protects the adult ovary from the testicular pathway Gonadal sexual fate in mammals is determined during embryonic development and must be actively maintained in adulthood. In the mouse ovary, oestrogen receptors and FOXL2 protect ovarian granulosa cells from transdifferentiation into Sertoli cells, their testicular counterpart. However, the mechanism underlying their protective effect is unknown. Here, we show that TRIM28 is required to prevent female-to-male sex reversal of the mouse ovary after birth. We found that upon loss of Trim28, ovarian granulosa cells transdifferentiate to Sertoli cells through an intermediate cell type, different from gonadal embryonic progenitors. TRIM28 is recruited on chromatin in the proximity of FOXL2 to maintain the ovarian pathway and to repress testicular-specific genes. The role of TRIM28 in ovarian maintenance depends on its E3-SUMO ligase activity that regulates the sex-specific SUMOylation profile of ovarian-specific genes. Our study identifies TRIM28 as a key factor in protecting the adult ovary from the testicular pathway.