Distribution and molecular evolution of the anti-CRISPR family AcrIF7 Anti-CRISPRs are proteins capable of blocking CRISPR-Cas systems and typically encoded in mobile genetic elements. Since their discovery, numerous anti-CRISPR families have been identified. However, little is known about the distribution and sequence diversity of members within a family, nor how these traits influence the anti-CRISPRs function and evolution. Here we use AcrIF7 to explore the dissemination and molecular evolution of an anti-CRISPR family. We uncovered five sub-clusters and prevalent anti-CRISPR variants within the group. Remarkably, AcrIF7 homologs display high similarity despite their broad geographical, ecological and temporal distribution. Although mainly associated with Pseudomonas aeruginosa, AcrIF7 was identified in distinct genetic backgrounds indicating horizontal dissemination, primarily by phages. Using mutagenesis, we recreated variation observed in databases but also extended the sequence diversity of the group. Characterisation of the variants identified residues key for the anti-CRISPR function and other contributing to its mutational tolerance. Moreover, molecular docking revealed that variants with affected function lose key interactions with its CRISPR-Cas target. Analysis of publicly available data and the generated variants suggests that the dominant AcrIF7 variant corresponds to the minimal and optimal anti-CRISPR selected in the family. Our study provides a blueprint to investigate the molecular evolution of anti-CRISPR families.
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.
Probing effects of the SARS-CoV-2 E protein on membrane curvature and intracellular calcium SARS-CoV-2 contains four structural proteins in its genome. These proteins aid in the assembly and budding of new virions at the ER-Golgi intermediate compartment (ERGIC). Current fundamental research efforts largely focus on one of these proteins - the spike (S) protein. Since successful antiviral therapies are likely to target multiple viral components, there is considerable interest in understanding the biophysical role of its other structural proteins, in particular structural membrane proteins. Here, we have focused our efforts on the characterization of the full-length envelope (E) protein from SARS-CoV-2, combining experimental and computational approaches. Recombinant expression of the full-length E protein from SARS-CoV-2 reveals that this membrane protein is capable of independent multimerization, possibly as a tetrameric or smaller species. Fluorescence microscopy shows that the protein localizes intracellularly, and coarse-grained MD simulations indicate it causes bending of the surrounding lipid bilayer, corroborating a potential role for the E protein in viral budding. Although we did not find robust electrophysiological evidence of ion-channel activity, cells transfected with the E protein exhibited reduced intracellular Ca2+, which may further promote viral replication. However, our atomistic MD simulations revealed that previous NMR structures are relatively unstable, and result in models incapable of ion conduction. Our study highlights the importance of using high-resolution structural data obtained from a full-length protein to gain detailed molecular insights, and eventually permitting virtual drug screening. A nontoxic tetrodotoxin analog, 5,6,11-trideoxytetrodotoxin, as an olfactory chemoattractant for the grass puffer Toxic puffers contain the potent neurotoxin, tetrodotoxin (TTX). Although TTX is considered to serve as a defense substance, previous behavioral studies have demonstrated that TTX acts as an attractive pheromone for some toxic puffers. To elucidate the physiological mechanism of putative pheromonal action of TTX, we examined whether grass puffers Takifugu alboplumbeus can detect TTX. Electroolfactogram (EOG) results suggest that the olfactory epithelium (OE) of grass puffers responded to a type of TTX analog (5,6,11-trideoxyTTX), although it did not respond to TTX. We also examined the attractive action of 5,6,11-trideoxyTTX on grass puffers by recording their swimming behavior under dark conditions. Grass puffers preferred to stay on the side of the aquarium where 5,6,11-trideoxyTTX was administered, and their swimming speed decreased. Additionally, odorant-induced labeling of olfactory sensory neurons by immunohistochemistry against neural activity marker (phosphorylated extracellular signal regulated kinase; pERK) revealed that labeled olfactory sensory neurons were localized in the region surrounding "islets" where there was considered as nonsensory epithelium. 5,6,11-trideoxyTTX has been known to accumulate in grass puffers, but its toxicity is much lower (almost nontoxic) than TTX. Our results suggest that toxic puffers may positively use this unknown function TTX analog as an olfactory chemoattractant. Titers of antibodies the receptor-binding domain (RBD) of ancestral SARS-CoV-2 are predictive for levels of neutralizing antibodies to multiple variants Diagnostic assays currently used to monitor the efficacy of COVID-19 vaccines measure levels of antibodies to the receptor-binding domain of ancestral SARS-CoV-2 (RBDwt). However, the predictive value for protection against new variants of concern (VOCs) has not been firmly established. Here, we used bead-based arrays and flow cytometry to measure binding of antibodies to spike proteins and receptor-binding domains (RBDs) from VOCs in 12,000 sera. Effects of sera on RBD-ACE2 interactions were measured as a proxy for neutralizing antibodies. The samples were obtained from healthy individuals or patients on immunosuppressive therapy who had received two to four doses of COVID-19 vaccines and from COVID-19 convalescents. The results show that anti-RBDwt titers correlate with the levels of binding- and neutralizing antibodies against the Alpha, Beta, Gamma, Delta, Epsilon and Omicron variants. The benefit of multiplexed analysis lies in the ability to measure a wide range of anti-RBD titers using a single dilution of serum for each assay. The reactivity patterns also yield an internal reference for neutralizing activity and binding antibody units per milliliter (BAU/ml). Results obtained with sera from vaccinated healthy individuals and patients confirmed and extended results from previous studies on time-dependent waning of antibody levels and effects of immunosuppressive agents. We conclude that anti-RBDwt titers correlate with levels of neutralizing antibodies against VOCs and propose that our method may be implemented to enhance the precision and throughput of immunomonitoring. Interplay between disordered regions in RNAs and proteins modulates interactions within stress granules and processing bodies Condensation, or liquid-like phase separation, is a phenomenon indispensable for the spatiotemporal regulation of molecules within the cell. Recent studies indicate that the composition and molecular organization of phase-separated organelles such as Stress Granules (SGs) and Processing Bodies (PBs) are highly variable and dynamic. A dense contact network involving both RNAs and proteins controls the formation of SGs and PBs and an intricate molecular architecture, at present poorly understood, guarantees that these assemblies sense and adapt to different stresses and environmental changes. Here, we investigated the physico-chemical properties of SGs and PBs components and studied the architecture of their interaction networks. We found that proteins and RNAs establishing the largest amount of contacts in SGs and PBs have distinct structural properties and intrinsic disorder is enriched in all protein-RNA, protein-protein and RNA-RNA interaction networks. The increase of disorder in proteins is accompanied by an enrichment in single-stranded regions of RNA binding partners. Our results suggest that SGs and PBs quickly assemble and disassemble through fuzzy-like dynamic contacts modulated by unfolded domains of their components.

Research HighlightsO_LIWe systematically studied RNA-RNA, protein-protein and RNA-protein interaction networks in stress granules and processing bodies;
C_LIO_LIRNAs enriched in stress granules and processing bodies are more single-stranded and form a large number of contacts with both proteins and RNAs;
C_LIO_LIProteins in stress granules and processing bodies are less structured and contact larger amounts of single-stranded RNAs.
C_LI Here's the twist: How the brain updates the representations of naturalistic events as our understanding of the past changes The brain actively reshapes our understanding of past events in light of new incoming information. In the current study, we ask how the brain supports this updating process during the encoding and recall of naturalistic stimuli. One group of participants watched a movie ("The Sixth Sense") with a cinematic "twist" at the end that dramatically changed the interpretation of previous events. Next, participants were asked to verbally recall the movie events, taking into account the new "twist" information. Most participants updated their recall to incorporate the twist. Two additional groups recalled the movie without having to update their memories during recall: one group never saw the twist; another group was exposed to the twist prior to the beginning of the movie, and thus the twist information was incorporated both during encoding and recall. We found that providing participants with information about the twist beforehand altered neural response patterns during movie-viewing in the default mode network (DMN). Moreover, presenting participants with the twist at the end of the movie changed the neural representation of the previously-encoded information during recall in a subset of DMN regions. Further evidence for this transformation was obtained by comparing the neural activation patterns during encoding and recall and correlating them with behavioral signatures of memory updating. Our results demonstrate that neural representations of past events encoded in the DMN are dynamically integrated with new information that reshapes our understanding in natural contexts. STING activation promotes robust immune response and NK cell-mediated tumor regression in glioblastoma models Immunotherapy has had a tremendous impact on cancer treatment in the past decade, with hitherto unseen responses at advanced and metastatic stages of the disease. However, the aggressive brain tumor glioblastoma (GBM) is highly immunosuppressive and remains largely refractory to current immunotherapeutic approaches. The cGAS-STING cytoplasmic double stranded DNA (dsDNA) sensing pathway has emerged as a next-generation immunotherapy target with potent local immune stimulatory properties.

Here, we investigated the status of the STING pathway in GBM and the modulation of the brain tumor microenvironment (TME) with the STING agonist ADU-S100. Our data reveal the presence of STING in human GBM specimens, where it stains strongly in the tumor vasculature. We show that human GBM explants can respond to STING agonist treatment by secretion of inflammatory cytokines. In murine GBM models, we show a profound shift in the tumor immune landscape after STING agonist treatment, with massive infiltration of the tumor-bearing hemisphere with innate immune cells including inflammatory macrophages, neutrophils and NK populations. Treatment of established murine intracranial GL261 and CT-2A tumors by biodegradable ADU-S100-loaded intracranial implants demonstrated a significant increase in survival in both models and long-term survival with immune memory in GL261. Responses to treatment were abolished by NK cell depletion. This study reveals therapeutic potential and deep remodeling of the TME by STING activation in GBM and warrants the further examination of STING agonists alone or in combination with other immunotherapies such as cancer vaccines, CAR T cells, NK therapies or immune checkpoint blockade.

Significance statementModulation of the immune microenvironment is critical for immunosuppressive and therapy refractory tumors like glioblastoma. Activation of the STING pathway deeply remodels the brain tumor environment and attracts innate immune cells and natural killer cell populations, producing a robust antitumor effect with long-term immune memory. We further show that human glioblastoma tissue can respond to the therapy and lay the foundations for combined intracranial immunotherapies by using crosslinked biodegradable brain implants. Embryonic exposure to cannabidiol disrupts active neural circuits, an effect increased by Δ9-tetrahydrocannabinol and involving CB1R and CB2R in zebrafish 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. Reliance of neuronal gene expression on cohesin scales with chromatin loop length Cohesin and CTCF are major drivers of 3D genome organization, but their role in neurons is still emerging. Here we show a prominent role for cohesin in the expression of genes that facilitate neuronal maturation and homeostasis. Unexpectedly, we observed two major classes of activity-regulated genes with distinct reliance on cohesin in primary cortical neurons. Immediate early genes remained fully inducible by KCl and BDNF, and short-range enhancer-promoter contacts at the Immediate early gene Fos formed robustly in the absence of cohesin. In contrast, cohesin was required for full expression of a subset of secondary response genes characterised by long-range chromatin contacts. Cohesin-dependence of constitutive neuronal genes with key functions in synaptic transmission and neurotransmitter signaling also scaled with chromatin loop length. Our data demonstrate that key genes required for the maturation and activation of primary cortical neurons depend cohesin for their full expression, and that the degree to which these genes rely on cohesin scales with the genomic distance traversed by their chromatin contacts. A network approach to genetic circuit designs As genetic circuits become more sophisticated, the size and complexity of data about their designs increases. This data captured goes beyond monolithic genetic sequences and towards circuit modularity and functional details, which are beneficial for analyzing circuit performance and establishing design automation techniques. However, the accessibility, visualisation and usability of design data (and metadata) have received relatively little attention to date. Here, we present a method to turn circuit designs into networks and showcase its potential to enhance the utility of design data. Since networks are dynamic structures, initial graphs can be interactively shaped into sub-networks of relevant information based on requirements such as abstraction, hierarchy and protein interactions. Additionally, several visual changes can be applied, such as colouring or clustering nodes based on types (e.g., genes or promoters), resulting in easier comprehension from a user perspective. This approach allows circuit designs to be coupled to other networks, such as metabolic pathways or implementation protocols captured in graph-like formats. Therefore, we advocate using networks to structure, access and improve synthetic biology information. Loss of histone methyltransferase SETD1B in oogenesis results in the redistribution of genomic histone 3 lysine 4 trimethylation Histone 3 lysine 4 trimethylation (H3K4me3) is an epigenetic mark found at gene promoters and CpG islands. H3K4me3 is essential for mammalian development, yet mechanisms underlying its genomic targeting are poorly understood. H3K4me3 methyltransferases SETD1B and MLL2 are essential for oogenesis. We investigated changes in H3K4me3 in Setd1b conditional knockout (cKO) oocytes using ultra-low input ChIP-seq, with comparisons to DNA methylation and gene expression analyses. H3K4me3 was redistributed in Setd1b cKO oocytes showing losses at active gene promoters associated with downregulated gene expression. Remarkably, many regions also gained H3K4me3, in particular those that were DNA hypomethylated, transcriptionally inactive and CpG-rich, which are hallmarks of MLL2 targets. Consequently, loss of SETD1B disrupts the balance between MLL2 and de novo DNA methyltransferases in determining the epigenetic landscape during oogenesis. Our work reveals two distinct, complementary mechanisms of genomic targeting of H3K4me3 in oogenesis, with SETD1B linked to gene expression and MLL2 to CpG content.

Graphical AbstractIn oogenesis, SETD1B and CXXC1 target H3K4me3 to actively transcribed gene promoters, while MLL2 targets transcriptionally inactive regions based on underlying CpG composition (upper panel). When SETD1B is ablated, H3K4me3 is lost at a subset of active promoters, resulting in downregulation of transcription (lower panel). Loss of SETD1B alters the activity of MLL2, permitting MLL2 to deposit H3K4me3 at CpG-rich regions, many of which should otherwise be DNA methylated. Thus, it is evident that MLL2 and de novo DNMTs compete for genomic occupancy late in oogenesis, and loss of SETD1B disrupts the balance of these mechanisms.

O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY Effect of tag attachment on the flight performance of five raptor species Bio-logging devices play a fundamental and indispensable role in movement ecology studies, particularly in the wild. However, researchers are becoming increasingly aware of the influence that attaching devices can have on animals, particularly on their behaviour, energy expenditure and survival. The way a device is attached to an animals body has also potential consequences for the collected data, and quantifying the type and magnitude of such potential effects is fundamental to enable researchers to combine and compare data from different studies, as much as it is to improve animal welfare.

For over two decades, large terrestrial birds have been in the focus of long-term movement ecology research, employing bio-logging devices attached with different types of harnesses. However, comparative studies investigating the effects of different harness types used on these species are scarce.

In this study, we tested for potential differences in data collected by two commonly used harness types, backpack and leg-loop, on the flight performance of 10 individuals from five raptor species, equipped with high resolution bio-logging devices, in the same area and time. We explored the effect of harness type on vertical speed, horizontal speed, glide ratio, height above sea level, distance travelled, proportion of soaring and flapping behaviour, and VeDBA (a proxy for energy expenditure) between and within individuals, all used as fine-scale measures of flight performance. Birds equipped with leg-loops climbed up to 0.65 ms-1 faster, reached 19% greater heights while soaring, and spent less time in active flight compared to birds equipped with backpacks, suggesting that backpack harnesses, compared to leg-loops, might cause additional drag affecting the birds flight performance. A lower rate of sinking while gliding, a slightly higher glide ratio, higher horizontal speed while soaring, and lower VeDBA, were also indicative of less drag using leg-loops. Our results add to the existing literature highlighting the design-related advantages of leg-loops, and support the use of leg-loops as a better alternative to backpack harnesses for large soaring birds, when possible. Our study also highlights how apparently small changes in device attachment can lead to notable improvements in tagging practice, with implications for animal welfare, data interpretation and comparability. The Skeletal Muscle Circadian Clock Regulates Titin Splicing Through RBM20 Circadian rhythms are maintained by a cell autonomous, transcriptional-translational feedback loop known as the molecular clock. While previous research suggests a role of the molecular clock in regulating skeletal muscle structure and function, no mechanisms have connected the molecular clock to sarcomere filaments. Utilizing inducible, skeletal muscle specific, Bmal1 knockout (iMSBmal1-/-) mice, we showed that knocking out skeletal muscle clock function alters titin isoform expression through both RNAseq and LC-MS. This alteration in titins spring length resulted in sarcomere length heterogeneity. We demonstrate the direct link between altered titin splicing and sarcomere length in vitro using U7 snRNPs that truncate the region of titin altered in iMSBmal1-/- muscle. Finally, we identified a mechanism whereby the skeletal muscle clock regulates titin isoform expression through transcriptional regulation of Rbm20, a potent splicing regulator of titin. Our findings demonstrate the importance of the skeletal muscle circadian clock in maintaining titin isoform through regulation of RBM20 expression. Because circadian rhythm disruption is a feature of many chronic diseases, our results highlight a novel pathway that could be targeted to maintain skeletal muscle structure and function in a range of pathologies. Modular Dynamic Biomolecular Modelling with Bond Graphs: The Unification of Stoichiometry, Thermodynamics, Kinetics and Data. Renewed interest in dynamic simulation models of biomolecular systems has arisen from advances in genome-wide measurement and applications of such models in biotechnology and synthetic biology. In particular, genome-scale models of cellular metabolism beyond the steady state are required in order to represent transient and dynamic regulatory properties of the system. Development of such whole-cell models requires new modelling approaches. Here we propose the energy-based bond graph methodology, which integrates stoichiometric models with thermo-dynamic principles and kinetic modelling. We demonstrate how the bond graph approach intrinsically enforces thermodynamic constraints, provides a modular approach to modelling, and gives a basis for estimation of model parameters leading to dynamic models of biomolecular systems. The approach is illustrated using a well-established stoichiometric model of Escherichia coli (E. coli) and published experimental data. Large Bi-ethnic Study of Plasma Proteome Leads to Comprehensive Mapping of cis-pQTL and Models for Proteome-wide Association Studies Improved understanding of genetic regulation of proteome can facilitate the identification of causal mechanisms for complex traits. We analyzed data on 4,657 plasma proteins from 7,213 European American (EA) and 1,871 African American (AA) individuals from the ARIC study, and further replicated findings on 467 AA individuals from the AASK study. Here we identified 2,004 proteins in EA and 1,618 in AA, with majority overlapping, which showed associations with common variants in cis-regions. Availability of AA samples led to smaller credible sets and significant number of population-specific cis-pQTLs. Elastic-net produced powerful models for protein prediction in both populations. An application of proteome-wide association studies (PWAS) to serum urate and gout, implicated several proteins, including IL1RN, revealing the promise of the drug anakinra to treat acute gout flares. Our study demonstrates the value of large and diverse ancestry study for genetic mechanisms of molecular phenotypes and their relationship with complex traits. 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. 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. Identifying and testing marker-trait associations for growth and phenology in three pine species: implications for genomic prediction In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 SNPs from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments. 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. Accurate prediction of functional states of cis-regulatory modules reveals the common epigenetic rules in humans and mice We proposed a two-step approach for predicting active cis-regulatory modules (CRMs) in a cell/tissue type. We first predict a map of CRM loci in the genome using all available transcription factor binding data in the organism, and then predict functional states of all the putative CRMs in any cell/tissue type using few epigenetic marks. We have recently developed a pipeline dePCRM2 for the first step, and now presented machine-learning methods for the second step. Our approach substantially outperforms existing methods. Our results suggest common epigenetic rules for defining functional states of CRMs in various cell/tissue types in humans and mice. Wolbachia load variation in Drosophila is more likely caused by drift than by host genetic factors Symbiosis is a continuum of long-term interactions ranging from mutualism to parasitism, according to the balance between costs and benefits for the protagonists. The density of endosymbionts is, in both cases, a key factor that determines both the transmission of symbionts and the host extended phenotype, and is thus tightly regulated within hosts. However, the evolutionary and molecular mechanisms underlying bacterial density regulation are currently poorly understood. In this context, the symbiosis between the fruit fly and its intracellular bacteria Wolbachia (wMelPop strain) is particularly interesting to study. Although vertically transmitted, the symbiont is pathogenic, and a positive correlation between virulence and wMelPop density is observed. In addition, the number of repeats of a bacterial genomic region -Octomom-is positively correlated with Wolbachia density, underlying a potential genetic mechanism that controls bacterial density. Interestingly, the number of repeats varies between host individuals, but most likely also within them. Such genetic heterogeneity within the host could promote conflicts between bacteria themselves and with the host, notably by increasing within-host competition between symbiont genotypes through a process analogous to the tragedy of the commons. To characterize the determinisms at play in the regulation of bacterial density, we first introgressed wMelPop in different genetic backgrounds of D. melanogaster, and found different density levels and Octomom copy numbers in each host lineage. To determine whether such variations reflect a host genetic determinism on density regulation through Octomom copy number selection, we replicated the introgressions and performed reciprocal crosses on the two Drosophila populations with the most extreme density levels. In both experiments, we detected an absence of directionality in the patterns of infection, associated with a strong instability of these patterns across generations. Given that bacterial density was highly correlated with Octomom copy numbers in all experiments, these results rather suggest a strong influence of drift and a random increase in the frequency of certain bacterial variants. We then discuss how drift, both on the symbiont population during transmission and on the host population, could limit the efficiency of selection in such a symbiotic system, and the consequences of drift on the regulation of density and composition of bacterial populations. Drug screening at single-organoid resolution via bioprinting and interferometry High-throughput drug screening is an established approach to investigate tumor biology and identify therapeutic leads. Traditional platforms for high-throughput screening use two-dimensional cultures of immortalized cell lines which do not accurately reflect the biology of human tumors. More clinically relevant model systems, such as three-dimensional tumor organoids, can be difficult to screen and scale. For example, manually seeded organoids coupled to destructive endpoint assays allow for the characterization of response to treatment, but do not capture the transitory changes and intra-sample heterogeneity underlying clinically observed resistance to therapy. We therefore developed a pipeline to generate bioprinted tumor organoids linked to label-free, real-time imaging via high-speed live cell interferometry (HSLCI) and machine learning-based quantitation of individual organoids. Bioprinting cells gives rise to 3D organoid structures that preserve tumor histology and gene expression. HSLCI imaging in tandem with machine learning-based image segmentation and organoid classification tools enables accurate, label-free parallel mass measurements for thousands of bioprinted organoids. We demonstrate that our method quantitatively identifies individual organoids as insensitive, transiently sensitive, or persistently sensitive to specific treatments. This opens new avenues for rapid, actionable therapeutic selection using automated tumor organoid screening. Identification of a Golgi-localized altORF peptide reveals a minimal Golgi targeting motif Prior work has identified signal sequences and motifs that are necessary and sufficient to target proteins to specific subcellular regions and organelles such as the plasma membrane, nucleus, Endoplasmic Reticulum, and mitochondria. In contrast, minimal sequence motifs that are sufficient for Golgi localization remain largely elusive. In this work, we identified a 37 amino acid alternative open reading frame (altORF) within the mRNA of the centromere protein, CENP-R. This altORF peptide localizes specifically to the cytoplasmic surface of the Golgi apparatus. Through mutational analysis, we identify a minimal 10 amino acid sequence and a critical cysteine residue that are necessary and sufficient for Golgi localization. Pharmacological perturbations suggest that this peptide is modified through palmitoylation to promote its localization. Together, our work defines a minimal sequence that is sufficient for Golgi targeting and provide a valuable Golgi marker for live cell imaging. Scaling effects of temperature on parasitism from individuals to host-parasite systems Parasitism is expected to change in a warmer future, but whether warming leads to substantial increases in parasitism remains unclear. Understanding how warming effects on parasitism in individual hosts (e.g., parasite load) translate to effects on population-level parasitism (e.g., prevalence, R0) remains a major knowledge gap. We analyzed the temperature dependence of parasitism at both host and population levels in thirteen empirical vector-borne host-parasite systems and found a strong, significant positive correlation between the thermal optima of individual- and population-level parasitism. We also found a significant, positive correlation in eleven environmentally-transmitted parasite systems, though several of these systems exhibited thermal optima >5{o}C apart between individual and population levels. Similarly, parasitism thermal optima were close to host performance thermal optima in vector-borne systems but not in environmentally-transmitted systems. We then adapted and simulated simple models for both transmission modes and found a similar pattern to the empirical systems: thermal optima in vector-borne systems were more strongly correlated across scales compared to environmentally-transmitted systems. Generally, our results suggest that information on the temperature-dependence, and specifically the thermal optimum, at either the individual- or population-level should provide a useful--though not quantitatively exact--baseline for predicting temperature dependence at the other level, especially in vector-borne parasite systems. Environmentally-transmitted parasitism may operate by a different set of rules, in which temperature-dependence is decoupled in some systems, requiring the need for trait-based studies of temperature-dependence at individual and population levels. DeGlyPHER: an ultrasensitive method for analysis of viral spike N-glycoforms Viruses can evade the host immune system by displaying numerous glycans on their surface "spike-proteins" that cover immune epitopes. We have developed an ultrasensitive "single pot" method to assess glycan occupancy and the extent of glycan processing from high-mannose to complex forms at each N-glycosylation site. Though aimed at characterizing glycosylation of viral spike-proteins as potential vaccines, this method is applicable for analysis of site-specific glycosylation of any glycoprotein.