Distinct mutations and lineages of SARS-CoV-2 virus in the early phase of COVID-19 pandemic and subsequent global expansion A novel coronavirus, SARS-CoV-2, has caused over 190 million cases and over 4 million deaths worldwide since it occurred in December 2019 in Wuhan, China. Here we conceptualized the temporospatial evolutionary and expansion dynamics of SARS-CoV-2 by taking a series of cross-sectional view of viral genomes from early outbreak in January 2020 in Wuhan to early phase of global ignition in early April, and finally to the subsequent global expansion by late December 2020. Based on the phylogenetic analysis of the early patients in Wuhan, Wuhan/WH04/2020 is supposed to be a more appropriate reference genome of SARS-CoV-2, instead of the first sequenced genome Wuhan-Hu-1. By scrutinizing the cases from the very early outbreak, we found a viral genotype from the Seafood Market in Wuhan featured with two concurrent mutations (i.e. M type) had become the overwhelmingly dominant genotype (95.3%) of the pandemic one year later. By analyzing 4,013 SARS-CoV-2 genomes from different continents by early April, we were able to interrogate the viral genomic composition dynamics of initial phase of global ignition over a timespan of 14-week. 11 major viral genotypes with unique geographic distributions were also identified. WE1 type, a descendant of M and predominantly witnessed in western Europe, consisted a half of all the cases (50.2%) at the time. The mutations of major genotypes at the same hierarchical level were mutually exclusive, which implying that various genotypes bearing the specific mutations were propagated during human-to-human transmission, not by accumulating hot-spot mutations during the replication of individual viral genomes. As the pandemic was unfolding, we also used the same approach to analyze 261,323 SARS-CoV-2 genomes from the world since the outbreak in Wuhan (i.e. including all the publicly available viral genomes) in order to recapitulate our findings over one-year timespan. By 25 December 2020, 95.3% of global cases were M type and 93.0% of M-type cases were WE1. In fact, at present all the four variants of concern (VOC) are the descendants of WE1 type. This study demonstrates the viral genotypes can be utilized as molecular barcodes in combination with epidemiologic data to monitor the spreading routes of the pandemic and evaluate the effectiveness of control measures. Moreover, the dynamics of viral mutational spectrum in the study may help the early identification of new strains in patients to reduce further spread of infection, guide the development of molecular diagnosis and vaccines against COVID-19, and help assess their accuracy and efficacy in real world at real time.
The first sheep graph-based pan-genome reveals the spectrum of structural variations and their effects on tail phenotypes Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here, we assembled 26 haplotype-resolved genome assemblies from 13 genetically diverse sheep using PacBio HiFi sequencing. We constructed a graph-based ovine pan-genome and discovered 142,422 biallelic insertions and deletions, 7,028 divergent alleles and 13,419 multiallelic variations. We then used a graph-based approach to genotype the biallelic SVs in 684 individuals from 45 domestic breeds and two wild species. Integration with RNA-seq data allows to identify candidate expression-associated SVs. We demonstrate a direct link of SVs and phenotypes by localizing the putative causative insertion in HOXB13 gene responsible for the long-tail trait and identifying multiple large SVs associated with the fat-tail. Beyond generating a benchmark resource for ovine structural variants, our study highlights that animal genetic research will greatly benefit from using a pan-genome graph rather than a single reference genome.
Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. ksrates: positioning whole-genome duplications relative to speciation events in KS distributions SummaryWe present ksrates, a user-friendly command-line tool to position ancient whole-genome duplication (WGD) events with respect to speciation events in a phylogeny by comparing paralog and ortholog KS distributions derived from genomic or transcriptomic sequences, while adjusting for substitution rate differences among the lineages involved.

Availability and implementationksrates is implemented in Python 3 and as a Nextflow pipeline. The source code, Singularity and Docker containers, documentation and tutorial are available via https://github.com/VIB-PSB/ksrates.

Contactsteven.maere@ugent.vib.be, rolf.lohaus@ugent.vib.be Early functional mismatch between breast cancer cells and their tumour microenvironment suppresses long term growth Cancer cells thrive embedded in a fine-tuned cellular and extracellular environment or tumour microenvironment (TME). There is a general understanding of a co-evolution between cancer cells and their surrounding TME, pointing at a functional connection between cancer cells characteristics and the perturbations induced in their surrounding tissue. However, whether this functional connection needs to be set from the start or if aggressive cancer cells can always be dominating their microenvironment has never been formally proven with a dedicated experimental setting where malignant cells can be challenged to grow in a different TME from the one they would naturally create. Here we generated an experimental setting where we transiently perturb the secretory profile of aggressive breast cancer cells without affecting their intrinsic growth ability. This led to the initial establishment of an atypical TME. Interestingly, even if initially tumours are formed, this atypical TME evolves to impair long term in vivo cancer growth. Using a combination of in vivo transcriptomics, protein arrays and in vitro co-cultures, we found that the atypical TME culminates in the infiltration of macrophages with STAT1high activity. These macrophages show strong anti-tumoural functions which reduce long-term tumour growth, despite lacking canonical M1 markers. Importantly, gene signatures of the mesenchymal compartment of the TME, as well as the anti-tumoural macrophages show striking prognostic power-correlating with less aggressive human breast cancers. Robust and efficient online auditory psychophysics Most human auditory psychophysics research has historically been conducted in carefully controlled environments with calibrated audio-equipment, and over potentially hours of repetitive testing with expert listeners. Here, we operationally define such conditions as having high  auditory hygiene. Hereof, conducting auditory psychophysical paradigms online presents a serious challenge, in that results may hinge on absolute sound presentation level, reliably estimated perceptual thresholds, and sustained motivation and attention. We introduce a set of procedures that address these challenges and facilitate auditory hygiene for online auditory psychophysics. First, we establish a simple means of setting sound presentation levels. Across a set of four level-setting conditions conducted in person, we demonstrate the stability and robustness of this amplitude setting procedure in open air and controlled settings. Second, we test participants tone-in-noise thresholds using widely adopted online experiment platforms and demonstrate that reliable threshold estimates can be derived online in approximately one minute of testing. Third, using these amplitude and threshold setting procedures to establish participant-specific stimulus conditions, we show that an online implementation of the classic probe-signal paradigm can be used to demonstrate frequency-selective attention on an individual participant basis, using a third of the trials used in recent in-lab experiments. Finally, we show how threshold and attentional measures relate to well-validated assays of online participants in-task motivation, fatigue, and confidence. This demonstrates the promise of online auditory psychophysics for asking new auditory neuroscience questions quickly, efficiently, and with more diverse samples. Code for the tests is publicly available through Pavlovia and Gorilla. Design of a radial multi-offset detection pattern for in vivo phase contrast imaging of the retinal ganglion cells in humans Previous work has shown that multi-offset detection in adaptive optics scanning laser ophthalmoscopy (AOSLO) can be used to image retinal ganglion cells (RGCs) in monkeys and humans. However, though images of RGCs in anesthetized monkeys with high light levels produced high contrast images of RGCs, images from humans failed to reach the same contrast due to several drawbacks in the previous dual-wavelength multi-offset approach. Our aim here was to design and build a multi-offset detection pattern for humans at safe light levels that could reveal the retinal ganglion cell layer neurons with a contrast, robustness and acquisition time approaching results only previously obtained in monkeys. Here, we present a new imaging system using only one light source, compared to the previous dual-wavelength used on monkeys. Our single-wavelength solution allows for increased light power and eliminates problematic chromatic aberrations. Then, we demonstrate that a radial multi-offset detection pattern with an offset distance of 8-10 Airy Disk Diameter (ADD) is optimal to detect photons multiply scattered in all directions from RGCs thereby enhancing their contrast. This new setup and image processing pipeline led to improved imaging of retinal ganglion cells using multi-offset imaging in AOSLO. Phylogenomic evidence for the Origin of Obligately Anaerobic Anammox Bacteria around the Great Oxidation Event The anaerobic ammonium oxidation (anammox) bacteria could transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for up to half of the global nitrogen loss in surface environments. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remains enigmatic. Here, we performed comprehensive phylogenomic analysis and showed a single origin of anammox bacteria within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates, which includes implementing both a traditional cyanobacteria-based and a recently developed mitochondria-based approach, we estimated that anammox bacteria originated at early Proterozoic and most likely around the so-called Great Oxidation Event (GOE; 2.32 to 2.5 billion years ago [Ga]) which fundamentally changed global biogeochemical cycles. We further showed that during the origin of anammox bacteria, genes involved in oxidative stress, bioenergetics and anammox granules formation were recruited, which might have contributed to their survival on an increasingly oxic Earth. Our findings suggest the rising levels of atmospheric oxygen, which made nitrite increasingly available, was a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link the GOE to the evolution of obligate anaerobic bacteria. Pathfinder: A gamified measure to integrate general cognitive ability into the biological, medical, and behavioural sciences Genome-wide association (GWA) studies have uncovered DNA variants associated with individual differences in general cognitive ability (g), but these are far from capturing heritability estimates obtained from twin studies. A major barrier to finding more of this  missing heritability is assessment - the use of diverse measures across GWA studies as well as time and cost of assessment. In a series of four studies, we created a 15-minute (40-item), online, gamified measure of g that is highly reliable (alpha = .78; two-week test-retest reliability = .88), psychometrically valid and scalable; we called this new measure Pathfinder. In a fifth study, we administered this measure to 4,751 young adults from the Twins Early Development Study. This novel g measure, which also yields reliable verbal and nonverbal scores, correlated substantially with standard measures of g collected at previous ages (r ranging from .42 at age 7 to .57 at age 16). Pathfinder showed substantial twin heritability (.57, 95% CIs = .43, .68) and SNP heritability (.37, 95% CIs = .04, .70). A polygenic score computed from GWA studies of five cognitive and educational traits accounted for 12% of the variation in g, the strongest DNA-based prediction of g to date. Widespread use of this engaging new measure will advance research not only in genomics but throughout the biological, medical, and behavioural sciences. Increasing triacylglycerol formation and lipid storage by unsaturated lipids protects renal proximal tubules in diabetes In diabetic patients, dyslipidemia frequently contributes to organ damage such as diabetic kidney disease (DKD). DKD is associated with excessive renal deposition of triacylglycerol (TAG) in lipid droplets (LD). Yet, it is unclear whether LDs play a protective or damaging role and how this might be influenced by dietary patterns. By using a diabetes mouse model, we find here that high fat diet enriched in the unsaturated oleic acid (OA) caused more lipid storage in LDs in renal proximal tubular cells (PTC) but less tubular damage than a corresponding butter diet with the saturated palmitic acid (PA). Mechanistically, we identify endoplasmic reticulum (ER) stress as the main cause of PA-induced PTC injury. ER stress is caused by elevated cellular levels of saturated TAG precursors and to higher membrane order in the ER. The resulting cell death is preceded by a transcriptional rewiring of phospholipid metabolism. Simultaneous addition of OA rescues the cytotoxic effects by normalizing membrane order and by increasing the total TAG amount. The latter also stimulates the formation of LDs that in turn can release unsaturated lipids upon demand by lipolysis. Our study thus clarifies mechanisms underlying PA-induced cell stress in PTCs and emphasizes the importance of olive oil for the prevention of DKD. G protein-coupled estrogen receptor activation by Bisphenol-A disrupts protection from apoptosis conferred by estrogen receptors ERα and ERβ in pancreatic beta cells 17{beta}-estradiol protects pancreatic {beta}-cells from apoptosis via the estrogen receptors ER, ER{beta} and GPER. Conversely, the endocrine disruptor Bisphenol-A (BPA), which exerts multiple effects in this cell type via the same estrogen receptors, increased basal apoptosis. The molecular initiated events that trigger these opposite actions have yet to be identified. We demonstrated that combined genetic downregulation and pharmacological blockade of each estrogen receptor increased apoptosis to a different extent. The increase in apoptosis induced by BPA was diminished by the pharmacological blockade or the genetic silencing of GPER, and it was partially reproduced by the GPER agonist G1. BPA and G1-induced apoptosis were abolished upon pharmacological inhibition, silencing of ER and ER{beta}, or in dispersed islet cells from ER{beta} knockout (BERKO) mice. Yet, the ER and ER{beta} agonists, PPT and DPN, respectively, had no effect on beta cell viability. To exert their biological actions, ER and ER{beta} form homodimers and heterodimers. Molecular dynamic simulations together with proximity ligand assay and coimmunoprecipitation experiments indicated that the interaction of BPA with ER and ER{beta} as well as the GPER activation by G1 decreased ER{beta} heterodimers. We propose that ER{beta} heterodimers play an antiapoptotic role in beta cells and that BPA- and G1-induced decrease in ER{beta} heterodimers leads to beta cell apoptosis. Unveiling how different estrogenic chemicals affect the crosstalk among estrogen receptors should help to identify diabetogenic endocrine disruptors.

HighlightsO_LIPharmacological blockade and gene silencing of estrogen receptors ER, ER{beta} and GPER indicate that they are antiapoptotic in basal conditions.
C_LIO_LIGPER activation by G1 and BPA triggered apoptosis via a crosstalk with ER and ER{beta}.
C_LIO_LIBPA interaction with ER and ER{beta} as well as GPER activation decreased ER{beta} heterodimers, which was associated to increased apoptosis.
C_LIO_LIThis pathway represents a novel molecular initiating event underlying the pro-apoptotic effect of BPA
C_LIO_LIThe EndoC-{beta}H1 cell line may be a valid model of human {beta}-cells for identifying diabetogenic pollutants.
C_LI A novel nongenetic murine model of hyperglycemia and hyperlipidemia-associated aggravated atherosclerosis ObjectiveAtherosclerosis, the main pathology underlying cardiovascular diseases is accelerated in diabetic patients. Genetic mouse models require breeding efforts which are time-consuming and costly. Our aim was to establish a new nongenetic model of inducible metabolic risk factors that mimics hyperlipidemia, hyperglycemia, or both and allows the detection of phenotypic differences dependent on the metabolic stressor(s).

Methods and ResultsWild-type mice were injected with gain-of-function PCSK9D377Y (proprotein convertase subtilisin/kexin type 9) mutant adeno-associated viral particles (AAV) and streptozotocin and fed either a high-fat diet (HFD) for 12 or 20 weeks or a high-cholesterol/high-fat diet (Paigen diet, PD) for 8 weeks. To evaluate atherosclerosis, two different vascular sites (aortic sinus and the truncus of the brachiocephalic artery) were examined in the mice. Combined hyperlipidemic and hyperglycemic (HGHCi) mice fed a HFD or PD displayed characteristic features of aggravated atherosclerosis when compared to hyperlipidemia (HCi HFD or PD) mice alone. Atherosclerotic plaques of HGHCi HFD animals were larger, showed a less stable phenotype (measured by the increased necrotic core area, reduced fibrous cap thickness, and less -SMA-positive area) and had more inflammation (increased plasma IL-1{beta} level, aortic pro-inflammatory gene expression, and MOMA-2-positive cells in the BCA) after 20 weeks of HFD. Differences between the HGHCi and HCi HFD models were confirmed using RNA-seq analysis of aortic tissue, revealing that significantly more genes were dysregulated in mice with combined hyperlipidemia and hyperglycemia than in the hyperlipidemia-only group. The HGHCi-associated genes were related to pathways regulating inflammation (increased Cd68, iNos, and Tnfa expression) and extracellular matrix degradation (Adamts4 and Mmp14). When comparing HFD with PD, the PD aggravated atherosclerosis to a greater extent in mice and showed plaque formation after 8 weeks. Hyperlipidemic and hyperglycemic mice fed a PD (HGHCi PD) showed less collagen (Sirius red) and increased inflammation (CD68-positive cells) within aortic plaques than hyperlipidemic mice (HCi PD). HGHCi-PD mice represent a directly inducible hyperglycemic atherosclerosis model compared with HFD-fed mice, in which atherosclerosis is severe by 8 weeks.

ConclusionWe established a nongenetically inducible mouse model allowing comparative analyses of atherosclerosis in HCi and HGHCi conditions and its modification by diet, allowing analyses of multiple metabolic hits in mice. Pan-cancer HLA Gene-mediated Tumor Immunogenicity and Immune Evasion Human Leukocyte Antigen (HLA) expression contributes to the activation of anti-tumor immunity through interactions with T cell receptors. However, pan-cancer HLA-mediated immunogenicity and immunoediting mechanisms has not been systematically studied. In a retrospective analysis of 33 tumor types from the Cancer Genome Atlas, we uncovered HLA class I and class II differential expression, which outperformed traditional clinical metrics in predicting patient survival. We also characterized the distribution of HLA supertypes across cancers and showed that patients with high HLA allelic diversity and gene expression had better prognosis. Immune microenvironments with varied survival outcomes could be predicted using a neural network model trained on HLA expression data. Furthermore, we identified a subset of tumors which upregulated HLA class I but not class II genes and exploited HLA-mediated escape strategies. Our results suggest the potential of using HLA expression data to predict patient prognosis. Taken together, we emphasize the crucial role of HLA upregulation in shaping prolonged anti-tumor immunity.

SynopsisIn a retrospective analysis of 11080 patients of 33 TCGA cancer types, we showed HLA class I and class II differential expression shape various immune microenvironments and tumor immunoediting mechanisms, predicting tumor immunogenicity and patient survival. Estimating animal density with camera traps: a practitioner's guide of the REST model Camera traps are increasingly popular in wildlife research and have the potential to be reliable and cost-effective for estimating animal density. Although density estimation using this automatic technique has long been restricted to species with individually recognizable markings, several analytical approaches have been proposed to target animals lacking such markings (i.e., unmarked populations). Among these approaches, the random encounter and staying time (REST) model may be an efficient and cost-effective approach, even though the procedures for the implementations have not yet been shared with researchers. This paper presents a working protocol for implementing the REST model our research group has been developing. We also present the R code to perform parameter estimation with a maximum likelihood and Bayesian approach. We suggest that this model has potential for further development. We strongly hope that this paper will encourage many researchers to use the REST for density estimation in a wide variety of species in various habitats and make a significant contribution to advancing wildlife conservation and management. Phenomics data processing: Extracting dose-response curve parameters from high-resolution temperature courses and repeated field-based wheat height measurements Temperature is a main driver of plant growth and development. New phenotyping tools enable quantifying the temperature response of hundreds of genotypes. Yet, for field-derived data, temperature response modeling bears flaws and pitfalls concerning the interpretation of derived parameters. In this study, climate data from five growing seasons with differing temperature distributions served as starting point for a growth simulation of wheat stem elongation, based on a four-parametric temperature response function (Wang-Engel) including all cardinal temperatures. In a novel approach, we re-extracted dose-responses from the simulation by combining high-resolution (hours) temperature courses with low-resolution (days) height data. The collection of such data is common in field phenotyping platforms. To take advantage of the lack of supra-optimal temperatures during the stem elongation, simpler (linear and asymptotic) models to predict temperature-response parameters were investigated. The asymptotic model extracted the base temperature of growth and the maximum absolute growth rate with high precision, whereas simpler, linear models failed to do so. Additionally, the asymptotic model provided a proxy estimate for the optimum temperature. However, when including seasonally changing cardinal temperatures, the prediction accuracy of the asymptotic model was strongly reduced. In a field study with three winter wheat varieties, significant differences were found for all three asymptotic dose-response curve parameters. We conclude that the asymptotic model based on high-resolution temperature courses is suitable to extract meaningful parameters from field-based data. 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. Rapid expansion of SARS-CoV-2 variants of concern is a result of adaptive epistasis The SARS-CoV-2 pandemic recently entered an alarming new phase with the emergence of the variants of concern (VOC) and understanding their biology is paramount to predicting future ones. Current efforts mainly focus on mutations in the spike glycoprotein (S), but changes in other regions of the viral proteome are likely key. We analyzed more than 900,000 SARS-CoV-2 genomes with a computational systems biology approach including a haplotype network and protein structural analyses to reveal lineage-defining mutations and their critical functional attributes. Our results indicate that increased transmission is promoted by epistasis, i.e., combinations of mutations in S and other viral proteins. Mutations in the non-S proteins involve immune-antagonism and replication performance, suggesting convergent evolution. Furthermore, adaptive mutations appear in geographically disparate locations, suggesting that either independent, repeat mutation events or recombination among different strains are generating VOC. We demonstrate that recombination is a stronger hypothesis, and may be accelerating the emergence of VOC by bringing together cooperative mutations. This emphasizes the importance of a global response to stop the COVID-19 pandemic. Mouse CHD4-NURD is required for spermatogonia survival and normal meiotic progression Testis development and sustained germ cell production in adults rely on the establishment and maintenance of spermatogonia stem cells and their proper differentiation into spermatocytes. Chromatin remodeling complexes regulate critical processes during gamete development by restricting or promoting accessibility of DNA repair and gene expression machineries to the chromatin. Here, we investigated the role of CHD4 and CHD3 catalytic subunits of the NURD complex during spermatogenesis. Germ cell-specific deletion of Chd4 early in gametogenesis, but not Chd3, resulted in arrested early gamete development due to failed cell survival of neonate undifferentiated spermatogonia stem cell population. Candidate assessment revealed that CHD4 controls expression of Dmrt1 and its downstream target Plzf, both described as prominent regulators of spermatogonia stem cell maintenance. Our results show the requirement of CHD4 in mammalian gametogenesis pointing to functions in gene expression early in the process. Decoding silence in free recall In experiments on free recall from lists of items, not all memory retrievals are necessarily reported. Previous studies investigated unreported retrievals by attempting to induce their externalization. We show that, without any intervention, their statistics may be directly estimated through a model-free analysis of inter-response times - the silent intervals between recalls. A delay attributable to unreported recalls emerges in three situations: if the final item was already recalled ("silent recency effect"); if the item that, within the list, follows the latest recalled item was already recalled ("silent contiguity effect"); and in sequential recalls within highly performing trials ("sequential slowdown"). We then turn to reproducing all these effects by a minimal model where the discarding of memories ("bouncing") occurs either if they are repetitious or, in strategically organized trials, if they are not sequential. Based on our findings, we propose various approaches to further probing the submerged dynamics of memory retrieval. 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. Differential regulation of alternative promoters emerges from unified kinetics of enhancer-promoter interaction Many eukaryotic genes contain alternative promoters with distinct expression patterns. How these promoters are differentially regulated remains elusive. Here, we apply single-molecule imaging to quantify the transcriptional regulation of two alternative promoters (P1 and P2) of the Bicoid (Bcd) target gene hunchback in syncytial blastoderm Drosophila embryos. Contrary to the previous notion that Bcd only activates P2, we find that Bcd activates both promoters via the same two enhancers. P1 activation is less frequent and requires binding of more Bcd molecules than P2 activation. Using a theoretical model to relate promoter activity to enhancer states, we show that the two promoters follow common transcription kinetics driven by sequential Bcd binding at the two enhancers. Bcd binding at either enhancer primarily activates P2, while P1 activation relies more on Bcd binding at both enhancers. These results provide a quantitative framework for understanding the kinetic mechanisms of complex eukaryotic gene regulation. 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. Chimeric mRNA based COVID-19 vaccine induces protective immunity against Omicron and Delta The emerging SARS-CoV-2 variants of concern (VOCs) exhibit enhanced transmission and immune escape, reducing the efficacy and effectiveness of the two FDA-approved mRNA vaccines. Here, we explored various strategies to develop novel mRNAs vaccines to achieve safer and wider coverage of VOCs. Firstly, we constructed a cohort of mRNAs that feature a furin cleavage mutation in the spike (S) protein of predominant VOCs, including Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1) and Delta (B.1.617.2). Not present in the mRNA vaccines currently in use, the mutation abolished the cleavage between the S1 and S2 subunits, potentially enhancing the safety profile of the immunogen. Secondly, we systematically evaluated the induction of neutralizing antibodies (nAb) in vaccinated mice, and discovered that individual VOC mRNAs elicited strong neutralizing activity in a VOC-specific manner. Thirdly, the IgG produced in mice immunized with Beta-Furin and Washington (WA)-Furin mRNAs showed potent cross-reactivity with other VOCs, which was further corroborated by challenging vaccinated mice with the live virus of VOCs. However, neither WA-Furin nor Beta-Furin mRNA elicited strong neutralizing activity against the Omicron variant. Hence, we further developed an Omicron-specific mRNA vaccine that restored protection against the original and the sublineages of Omicron variant. Finally, to broaden the protection spectrum of the new Omicron mRNA vaccine, we tested the concept of bivalent immunogen. Instead of just fusing two RBDs head-to-tail, we for the first time constructed an mRNA-based chimeric immunogen by introducing the RBD of Delta variant into the entire S antigen of Omicron. The resultant chimeric mRNA was capable of inducing potent and broadly acting nAb against Omicron (both BA.1 and BA.2) and Delta, which paves the way to develop new vaccine candidate to target emerging variants in the future. Feline tooth resorption findings from a questionnaire survey on 8115 Finnish cats and a case-control study for a more defined population drawn from the survey Tooth resorption (TR) is one of the most common dental diseases in cats. It is a painful disease characterised by progressive dental destruction, which eventually results in loss of teeth. The aetiology of the disease remains unclear, but associations with old age, breed, other oral and dental diseases, and certain environmental factors have been suspected. We collected health and environmental information of 8115 Finnish cats from 41 breeds through an online survey for cat owners to study TR in Finland with an observational study. The frequency of TR that the owner reported as veterinarian-diagnosed was 3.9% in the entire data (316/8115) and 15% in cats that the owner reported as having veterinarian-diagnosed oral or dental disease (316/2070). We determined the risk factors and breed variation of feline TR by logistic regression with the entire data and a case-control study in a more defined population drawn from the entire data. In the entire data, we defined those cats that the owners had reported as being diagnosed with TR by a veterinarian as TR-cats and compared them to all other cats. In the case-control study, we included only cats with owner-reported veterinarian-diagnosed oral and dental diseases. We defined cats with TR-diagnosis made under sedation as TR-cases. Other veterinarian-diagnosed TR-cats were excluded. Cats with other oral and dental diseases were used as controls. The model with the entire population included 308 TR-cats compared with 7508 cats, and the case-control study included 198 TR-cats and 1683 controls. Although the diagnoses were reportedly made by a veterinarian, the owner answered the survey and might have understood, remembered, or interpreted the diagnoses in a way that could affect the results. The results indicated an increasing risk of TR with age. Dental calculus, gingivitis, and periodontitis were potential risk factors for TR. These findings and the interaction between dental calculus and gingivitis suggest that inflammatory changes caused by dental calculus increase the risk of TR. Keeping food available constantly was a potential protective factor. We found that Cornish Rex, European, and Ragdoll breeds were at higher risk for TR. The observed differences between breeds highlight a genetic contribution. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. Protein-protein interactions on membrane surfaces analysed using pull-downs with supported bilayers on silica beads Discovery-based proteomics workflows that identify novel interactors rely on immunoprecipitations or pull-downs with genetically-tagged bait proteins immobilized on appropriate matrices. But strategies to analyse protein interactions on a diffusible membrane surface combined with the practical ease of pull-downs remain unavailable. Such strategies are important to analyse protein complexes that mature in composition and stability because of diffusion-based encounter between participant proteins. Here, we describe a generic pull-down strategy to analyse such complexes using chelating lipid-containing supported bilayers formed on silica beads. These templates can display desired His-tagged bait proteins on a diffusible membrane surface. Using clathrin-mediated endocytosis as a paradigm, we find that the clathrin-binding adaptor protein epsin1 displayed as bait on these templates pulls down significantly higher amounts of clathrin from brain lysates than when immobilized on conventional matrices. Together, our results establish the potential of such templates as superior matrices for analysing protein-protein interactions and resultant complexes formed on membrane surfaces.