Changes in Behavior and Salivary Serum Amyloid A in cows with Subclinical Mastitis The aim of this study was to identify detailed changes in behavior, and in salivary serum amyloid A (SAA), associated with subclinical mastitis. This included standard sickness behaviors (which form part of an adaptive response to conserve energy, minimize heat loss and facilitate recovery following infection and inflammation) and lesser-studied behaviors, that are biologically predicted to change at low-level infection, and therefore particularly relevant for detecting subclinical mastitis (here labelled  luxury behaviors). SAA is a physiological marker of systemic inflammation, with levels in milk and serum already known to increase during subclinical mastitis. We investigated whether the same was true for SAA in cow saliva. Data were collected for n = 34 commercial barn-housed dairy cows, tested in 17 matched-pairs. Each pair comprised a cow with subclinical mastitis (SCM) and a healthy control (CTRL), identified using somatic cell count (SCC) (SCM: SCC >200 x1000 cells/ml; CTRL: SCC <100 x1000 cells/ml). SCM cows were selected for study ad-hoc, at which point they were paired with a CTRL cow, based upon parity and calving date; consequently, the full data set was accrued over several months. Data were collected for each pair over 3 days: SCC (l4:00-l5:00h) Day 1; behavior (24h from 00:00h) Day 2; salivary serum amyloid-A (SAA) Day 3. We report, for the first time, that an increase in salivary SAA occurs during subclinical mastitis; SAA was higher in SCM cows and demonstrated a positive (weak) correlation with SCC. The behavioral comparisons revealed that SCM cows had reductions in activity (behavioral transitions and distance moved), social exploration, social reactivity (here: likelihood to move away/be displaced following receipt of agonism), performance of social grooming and head butts, and the receipt of agonistic non-contact challenges. In addition, SCM cows received more head swipes, and spent a greater proportion of time lying with their head on their flank than CTRL cows. SCM cows also displayed an altered feeding pattern; they spent a greater proportion of feeding time in direct contact with two conspecifics, and a lower proportion of feeding time at self-locking feed barriers, than CTRL cows. Behavioral measures were found to correlate, albeit loosely, with serum SAA in a direction consistent with predictions for sickness behavior. These included positive correlations with both lying duration and the receipt of all agonistic behavior, and negative correlations with feeding, drinking, the performance of all social and all agonistic behavior, and social reactivity. We conclude that changes in salivary SAA, social behavior, and activity offer potential in the detection of subclinical mastitis and recommend further investigation to substantiate and refine our findings.
Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level.
Anaerobic Thermophilic Co-digestion of Sugarcane Biorefinery Residues Enhances Biogas Production This is a pioneer study evaluating the methane (CH4) production potential from residues of integrated 1st (vinasse and filter cake) and 2nd (deacetylation pretreatment liquor from straw) generation (1G2G) sugarcane biorefinery, providing a fully chemical characterization of them and their relation with the anaerobic digestion (AD) process. Small-scale assays provided fundamentals for basing the co-digestion optimization by assessing the optimal co-substrates synergistic conditions. Biochemical Methane Potential (BMP) tests showed co-digestion enhanced CH4 yield of isolated substrates, reaching up to 605 NmLCH4 gVS-1. The association of vinasse and deacetylation liquor as co-substrates increased the BMP by ~38% mostly by nutritionally benefiting the methanogenic activity. The kinetic analysis confirmed that the deacetylation liquor was the co-substrate responsible for improving the CH4 production in the co-digestion systems due to the highest CH4 conversion rate. The alkaline characteristic of the liquor (pH~12) also prevented alkalizing from being added to the co-digestion, an input that normally makes the process economically unfeasible to implement on an industrial scale due to the large quantities required for buffering the reactor. The filter cake had the lowest BMP (262 NmLCH4 gVS-1) and digestibility ([&le;] 40%), further limited by the required stirring to improve the mass transfer of biochemical reactions. The present study drives towards more sustainable use of vinasse, the most voluminous waste from the sugarcane industry, and lignin-rich residues derived from pre-treatment alkaline methods, aiming at an energy-efficient utilization, by at least 16% when compared to the traditional vinasse AD. The experimental and modeling elements from this work indicated the lignin-rich liquor is the main responsible for putting the co-digestion as a disruptive technological arrangement within the 1G2G sugarcane biorefineries, reinforcing the biogas production as the hub of the bioeconomy in the agroindustrial sector. Nuclear Progestin Receptor Mediated Linkage of Blood Coagulation and Ovulation Ovulation is a dramatic remodeling process that includes rupture of blood capillaries and clotting, but coagulation is not thought to directly regulate this process. Herein, we report remarkable increases of coagulation factors V (f5, [~]3145-fold) and tissue factor (f3a, [~]120-fold) in zebrafish ovarian follicle cells during ovulation. This increase was mediated through the nuclear progestin receptor (Pgr), which is essential for ovulation in zebrafish, and was totally abolished in ovarian follicular cells from pgr-/- mutants. In addition, promoter activities of f5 and f3a were significantly enhanced by progestin (DHP) via Pgr. Similar regulation of human F5 promoter activity was induced via human PGRB, suggesting a conserved mechanism. Site-directed mutagenesis of the zebrafish f5 promoter further demonstrated a direct regulation of coagulation factors via progestin response elements. Moreover, a stark increase of erythrocytes occurred in capillaries meshed in wildtype preovulatory follicles but was absent in pgr-/- mutants. Interestingly, anticoagulants significantly inhibited ovulation both in vitro and in vivo, respectively. Furthermore, reduced fecundity was observed in f5+/- female zebrafish. Taken together, our study provides plausible evidence for steroid regulation of coagulation factors, and a new hypothesis for blood clotting triggered ovulation in vertebrates. Global membrane tension is independent of polyacrylamide substrate stiffness Most animal cells are surrounded by a cell membrane and an underlying actomyosin cortex. Both structures are linked, and they are under tension. In-plane membrane tension and cortical tension both influence many cellular processes, including cell migration, division, and endocytosis. However, while actomyosin tension is regulated by substrate stiffness, how membrane tension responds to mechanical substrate properties is currently poorly understood. Here, we probed the effective membrane tension of neurons and fibroblasts cultured on glass and polyacrylamide substrates of varying stiffness using optical tweezers. In contrast to actomyosin-based traction forces, both peak forces and steady state tether forces of cells cultured on hydrogels were independent of substrate stiffness and did not change after blocking myosin II activity using blebbistatin, indicating that tether and traction forces are not directly linked. Peak forces in fibroblasts on hydrogels were about twice as high as those in neurons, indicating stronger membrane-cortex adhesion in fibroblasts. Steady state tether forces were generally higher in cells cultured on hydrogels than on glass, which we explain by a mechanical model. Our results provide new insights into the complex regulation of effective membrane tension and pave the way for a deeper understanding of the biological processes it instructs. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Lipid membranes modulate the activity of RNA through sequence-dependent interactions RNA is a ubiquitous biomolecule that can serve as both catalyst and information carrier. Understanding how RNA bioactivity is controlled is crucial for elucidating its physiological roles and potential applications in synthetic biology. Here we show that lipid membranes can act as RNA organization platforms, introducing a novel mechanism for ribo-regulation. The activity of R3C ribozyme can be modified by the presence of lipid membranes, with direct RNA-lipid interactions dependent on RNA nucleotide content, base pairing and length. In particular, the presence of guanine in short RNAs is crucial for RNA-lipid interactions, and G-quadruplex formation further promotes lipid binding. Lastly, by artificially modifying the R3C substrate sequence to enhance membrane binding we generated a lipid-sensitive ribozyme reaction with riboswitch-like behavior. These findings introduce RNA-lipid interactions as a tool for developing synthetic riboswitches and novel RNA-based lipid biosensors, and bear significant implications for RNA World scenarios for the origin of life. SIX1 Reprograms Myogenic Transcription Factors to Maintain the Rhabdomyosarcoma undifferentiated state Rhabdomyosarcoma (RMS) is a pediatric skeletal muscle sarcoma characterized by the expression of the myogenic-lineage transcription factors (TF) MYOD1 and MYOG. Despite high expression of these TFs, RMS cells fail to terminally differentiate, suggesting the presence of factors that alter their function. Here, we demonstrate that the developmental TF, SIX1, is highly expressed in RMS and is critical to maintain a muscle progenitor-like state. SIX1 loss induces terminal differentiation of RMS cells into myotube-like cells and dramatically impedes tumor growth in vivo. We show that SIX1 maintains the RMS undifferentiated state by controlling enhancer activity and MYOD1 occupancy at loci more permissive to tumor growth over terminal muscle differentiation. Finally, we demonstrate that a gene signature derived from SIX1 loss correlates with differentiation status in RMS and predicts RMS progression in human disease. Our findings demonstrate a master regulatory role for SIX1 in the repression of RMS differentiation via genome-wide alterations in MYOD1-mediated transcription.

HighlightsO_LISIX1 prevents differentiation in RMS while it promotes differentiation during normal development
C_LIO_LIFN-RMS are highly dependent on SIX1 for growth in both zebrafish and mouse xenograft models
C_LIO_LILoss of SIX1 alters the transcriptional landscape of RMS cells, inducing a growth to differentiation switch
C_LIO_LISIX1 knockdown in FN-RMS causes reduced super enhancer-based activity at stem-related genes and enhanced MYOD1 binding to differentiation loci, resulting in the activation of a myogenic differentiation program
C_LIO_LIA gene signature derived from SIX1 loss strongly correlates with myogenic differentiation status and is predictive of advanced RMS.
C_LI A chemical screen based on an interruption of zebrafish gastrulation identifies the HTR2C inhibitor Pizotifen as a suppressor of EMT-mediated metastasis Metastasis is responsible for approximately 90% of cancer-associated mortality but few models exist that allow for rapid and effective screening of anti-metastasis drugs. Current mouse models of metastasis are too expensive and time consuming to use for rapid and high-throughput screening. Therefore, we created a unique screening concept utilizing conserved mechanisms between zebrafish gastrulation and cancer metastasis for identification of potential anti-metastatic drugs. We hypothesized that small chemicals that interrupt zebrafish gastrulation might also suppress metastatic progression of cancer cells and developed a phenotype-based chemical screen to test the hypothesis. The screen used epiboly, the first morphogenetic movement in gastrulation, as a marker and enabled 100 chemicals to be tested in five hours. The screen tested 1280 FDA-approved drugs and identified Pizotifen, an antagonist for serotonin receptor 2C (HTR2C) as an epiboly-interrupting drug. Pharmacologic and genetic inhibition of HTR2C suppressed metastatic progression in a mouse model. Blocking HTR2C with Pizotifen restored epithelial properties to metastatic cells through inhibition of Wnt-signaling. In contrast, HTR2C induced epithelial to mesenchymal transition (EMT) through activation of Wnt-signaling and promoted metastatic dissemination of human cancer cells in a zebrafish xenotransplantation model. Taken together, our concept offers a novel platform for discovery of anti-metastasis drugs. SARS-CoV-2 B.1.617.2 Delta variant emergence, replication and immune evasion The SARS-CoV-2 B.1.617.2 (Delta) variant was first identified in the state of Maharashtra in late 2020 and spread throughout India, outcompeting pre-existing lineages including B.1.617.1 (Kappa) and B.1.1.7 (Alpha). In vitro, B.1.617.2 is 6-fold less sensitive to serum neutralising antibodies from recovered individuals, and 8-fold less sensitive to vaccine-elicited antibodies as compared to wild type Wuhan-1 bearing D614G. Serum neutralising titres against B.1.617.2 were lower in ChAdOx-1 versus BNT162b2 vaccinees. B.1.617.2 spike pseudotyped viruses exhibited compromised sensitivity to monoclonal antibodies against the receptor binding domain (RBD) and N-terminal domain (NTD), in particular to the clinically approved bamlavinimab and imdevimab monoclonal antibodies. B.1.617.2 demonstrated higher replication efficiency in both airway organoid and human airway epithelial systems as compared to B.1.1.7, associated with B.1.617.2 spike being in a predominantly cleaved state compared to B.1.1.7. Additionally we observed that B.1.617.2 had higher replication and spike mediated entry as compared to B.1.617.1, potentially explaining B.1.617.2 dominance. In an analysis of over 130 SARS-CoV-2 infected healthcare workers across three centres in India during a period of mixed lineage circulation, we observed substantially reduced ChAdOx-1 vaccine efficacy against B.1.617.2 relative to non-B.1.617.2. Compromised vaccine efficacy against the highly fit and immune evasive B.1.617.2 Delta variant warrants continued infection control measures in the post-vaccination era. Sarcoplasmic reticulum calcium handling in young DBA/2J vs C57BL/10 mdx models of Duchenne muscular dystrophy The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. TNFα increases Tyrosine Hydroxylase expression in human monocytes Most, if not all, peripheral immune cells in humans and animals express tyrosine hydroxylase (TH), the rate limiting enzyme in catecholamine synthesis. Since TH is typically studied in the context of brain catecholamine signaling, little is known about changes in TH production and function in peripheral immune cells. This knowledge gap is due, in part, to the lack of an adequately sensitive assay to measure TH in immune cells expressing lower TH levels compared to other TH expressing cells. Here, we report the development of a highly sensitive and reproducible Bio-ELISA to quantify picogram levels of TH in multiple model systems. We have applied this assay to monocytes isolated from blood of persons with Parkinsons disease (PD) and to age-matched, healthy controls. Our study unexpectedly revealed that PD patients monocytes express significantly higher levels of TH protein in peripheral monocytes relative to healthy controls. Tumor necrosis factor (TNF{delta}), a pro-inflammatory cytokine, has also been shown to be increased in the brains and peripheral circulation in human PD, as well as in animal models of PD. Therefore, we investigated a possible connection between higher levels of TH protein and the known increase in circulating TNF{delta} in PD. Monocytes isolated from healthy donors were treated with TNF{delta} or with TNF{delta} in the presence of an inhibitor. Tissue plasminogen activator (TPA) was used as a positive control. We observed that TNF{delta} stimulation increased both the number of TH+ monocytes and the quantity of TH per monocyte, without increasing the total numbers of monocytes. These results revealed that TNF{delta} could potentially modify monocytic TH production and serve a regulatory role in peripheral immune function. The development and application of a highly sensitive assay to quantify TH in both human and animal cells will provide a novel tool for further investigating possible PD immune regulatory pathways between brain and periphery. Basal Gp78-dependent mitophagy promotes mitochondrial health and limits mitochondrial ROS production Mitochondria are major sources of cytotoxic reactive oxygen species (ROS) that contribute to cancer progression. Mitophagy, the selective elimination of mitochondria by autophagy, monitors and maintains mitochondrial health and integrity, eliminating ROS-producing mitochondria. However, mechanisms underlying mitophagic control of mitochondrial homeostasis under basal conditions remain poorly understood. Gp78 E3 ubiquitin ligase is an endoplasmic reticulum membrane protein that induces mitochondrial fission and mitophagy of depolarized mitochondria. Here, we report that CRISPR/Cas9 knockout of Gp78 in HT-1080 fibrosarcoma cells increased mitochondrial volume and rendered cells resistant to carbonyl cyanide m-chlorophenyl hydrazone (CCCP)-induced mitophagy. These effects were phenocopied by knockdown of the essential autophagy protein ATG5 in wild-type HT-1080 cells. Use of the mito-Keima mitophagy probe confirmed that Gp78 promoted both basal and damage-induced mitophagy. Application of a spot detection algorithm (SPECHT) to GFP-mRFP tandem fluorescent-tagged LC3 (tfLC3)-positive autophagosomes reported elevated autophagosomal maturation in wild-type HT-1080 cells relative to Gp78 knockout cells, predominantly in proximity to mitochondria. Mitophagy inhibition by either Gp78 knockout or ATG5 knockdown reduced mitochondrial potential and increased mitochondrial ROS. Live cell analysis of tfLC3 in HT-1080 cells showed the preferential association of autophagosomes with mitochondria of reduced potential. Basal Gp78-dependent mitophagic flux is therefore selectively associated with reduced potential mitochondria promoting maintenance of a healthy mitochondrial population and limiting ROS production. Predicting the zoonotic capacity of mammals to transmit SARS-CoV-2 Back and forth transmission of SARS-CoV-2 between humans and animals may lead to wild reservoirs of virus that can endanger efforts toward long-term control of COVID-19 in people, and protecting vulnerable animal populations that are particularly susceptible to lethal disease. Predicting high risk host species is key to targeting field surveillance and lab experiments that validate host zoonotic potential. A major bottleneck to predicting animal hosts is the small number of species with available molecular information about the structure of ACE2, a key cellular receptor required for viral cell entry. We overcome this bottleneck by combining species ecological and biological traits with 3D modeling of virus and host cell protein interactions using machine learning methods. This approach enables predictions about the zoonotic capacity of SARS-CoV-2 for over 5,000 mammals -- an order of magnitude more species than previously possible. The high accuracy predictions achieved by this approach are strongly corroborated by in vivo empirical studies. We identify numerous common mammal species whose predicted zoonotic capacity and close proximity to humans may further enhance the risk of spillover and spillback transmission of SARS-CoV-2. Our results reveal high priority areas of geographic overlap between global COVID-19 hotspots and potential new mammal hosts of SARS-CoV-2. With molecular sequence data available for only a small fraction of potential host species, predictive modeling integrating data across multiple biological scales offers a conceptual advance that may expand our predictive capacity for zoonotic viruses with similarly unknown and potentially broad host ranges. Variation in female songbird state determines signal strength needed to evoke copulation It is the female response to male signals that determines courtship success. In most songbirds, females control reproduction via the copulation solicitation display (CSD), an innate, stereotyped posture produced in direct response to male displays. Because CSD can be elicited in the absence of males by the presentation of recorded song, CSD production enables investigations into the effects of underlying signal features and behavioral state on female mating preferences. Using computer vision to quantify CSD trajectory in female brown-headed cowbirds (Molothrus ater), we show that both song quality and a females internal state predict CSD production, as well as the onset latency and duration of the display. We also show that CSD can be produced in a graded fashion based on both signal strength and internal state. These results emphasize the importance of underlying receiver state in determining behavioral responses and suggest that female responsiveness acts in conjunction with male signal strength to determine the efficacy of male courtship. Vesalius: high-resolution in silico anatomization of Spatial Transcriptomic data using Image Analysis Characterization of tissue architecture promises to deliver insights into development, cell communication and disease. In silico spatial domain retrieval methods have been developed for spatial transcriptomics (ST) data assuming transcriptional similarity of neighboring barcodes. However, domain retrieval approaches with this assumption cannot work in complex tissues composed of multiple cell types. This task becomes especially challenging in cellular resolution ST methods. We developed Vesalius to decipher tissue anatomy from ST data by applying image processing technology. Vesalius uniquely detected territories composed of multiple cell types and successfully recovered tissue structures in high-resolution ST data including in mouse brain, embryo, liver, and colon. Utilizing this tissue architecture, Vesalius identified tissue morphology specific gene expression and regional specific gene expression changes for astrocytes, interneuron, oligodendrocytes, and entorhinal cells in the mouse brain. 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. Bacterial glycyl tRNA synthetase offers glimpses of ancestral protein topologies Aminoacyl tRNA synthetases (aaRSs) are among the proposed proteins present in the Last Universal Common Ancestor (LUCA). There are two types of glycyl tRNA synthetases (GlyRSs), from which the archaeal-eukaryal type is the one suggested to be present in LUCA. Here we solved the crystal structure of a complete bacterial glycyl tRNA synthetase (bacGlyRS) and show that indeed, bacGlyRS carries several structural signals that point it at the origin of all aaRSs. Furthermore, if bacGlyRS is ancestral, it should help to build a reliable Tree of Life (ToL). Given the modular nature of protein evolution, we used only two sub-domain segments with duplicated ancestral topologies, no detected orthologs and an assumed limited horizontal gene transfer (HGT). These motifs correspond to the non-specific RNA binding regions of contemporary bacGlyRS, archaeal CCA-adding enzyme (arch-CCAadd), and eukaryotic rRNA processing enzyme (euk-rRNA). The calculated, rooted bacterial ToL agrees with several phyla relationships unaccounted by the available trees. 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. Target Recognition in Tandem WW Domains: Complex Structures for Parallel and Antiparallel Ligand Orientation in h-FBP21 Tandem WW Protein-protein interactions often rely on specialized recognition domains, such as WW domains, which bind to specific proline-rich sequences. The specificity of these protein-protein interactions can be increased by tandem repeats, i.e. two WW domains connected by a linker. With a flexible linker, the WW domains can move freely with respect to each other. Additionally, the tandem WW domains can bind in two different orientations to their target sequences. This makes the elucidation of complex structures of tandem WW domains extremely challenging. Here, we identify and characterize two complex structures of the tandem WW domain of human formin-binding protein 21 and a peptide sequence from its natural binding partner, the core-splicing protein SmB/B. The two structures differ in the ligand orientation, and consequently also in the relative orientation of the two WW domains. We analyze and probe the interactions in the complexes by molecular simulations and NMR experiments. The workflow to identify the complex structures uses molecular simulations, density-based clustering and peptide docking. It is designed to systematically generate possible complex structures for repeats of recognition domains. These structures will help us to understand the synergistic and multivalency effects that generate the astonishing versatility and specificity of protein-protein interactions. ChIPflow:from raw data to epigenomic dynamics We present ChIPflow, a Snakemake-based pipeline for epigenomic data from the raw fastq files to the differential analysis. It can be applied to any chromatin factor, e.g. histone modification or transcription factor, which can be profiled with ChIP-seq. ChIPflow streamlines critical steps like the quality assessment of the immunoprecipitation using cross-correlation and the replicate comparison for both narrow and broad peaks. For the differential analysis ChIPflow provides linear and nonlinear methods for normalisation between samples as well as conservative and stringent models for estimating the variance and testing the significance of the observed binding/marking differences.

ChIPflow can process in parallel multiple chromatin factors with different experimental designs, number of biological replicates and/or conditions. It also facilitates the specific parametrisation of each dataset allowing both narrow or broad peak calling, as well as comparisons between the conditions using multiple statistical settings. Finally, complete reports are produced at the end of the bioinformatic and the statistical part of the analysis, which facilitate the data quality control and the interpretation of the results.

We explored the discriminative power of the statistical settings for the differential analysis, using a published dataset of three histone marks (H3K4me3, H3K27ac and H3K4me1) and two transcription factors (Oct4 and Klf4) profiled with ChIP-seq in two biological conditions (shControl and shUbc9). We show that distinct results are obtained depending on the sources of ChIP-seq variability and the dynamics of the chromatin factor under study. We propose that ChIPflow can be used to measure the richness of the epigenomic landscape underlying a biological process by identifying diverse regulatory regimes and the associated genes sets. An mRNA SARS-CoV-2 vaccine employing Charge-Altering Releasable Transporters with a TLR-9 agonist induces neutralizing antibodies and T cell memory The SARS-CoV-2 pandemic has necessitated the rapid development of prophylactic vaccines. Two mRNA vaccines have been approved for emergency use by the FDA and have demonstrated extraordinary effectiveness. The success of these mRNA vaccines establishes the speed of development and therapeutic potential of mRNA. These authorized vaccines encode full-length versions of the SARS-CoV-2 spike protein. They are formulated with Lipid Nanoparticle (LNP) delivery vehicles that have inherent immunostimulatory properties. Different vaccination strategies and alternative mRNA delivery vehicles would be desirable to ensure flexibility of future generations of SARS-CoV-2 vaccines and the development of mRNA vaccines in general.

Here, we report on the development of an alternative mRNA vaccine approach using a delivery vehicle called Charge-Altering Releasable Transporters (CARTs). Using these inherently nonimmunogenic vehicles we can tailor the vaccine immunogenicity by inclusion of co-formulated adjuvants such as oligodeoxynucleotides with CpG motifs (CpG-ODN). Mice vaccinated with the mRNA-CART vaccine developed therapeutically relevant levels of RBD-specific neutralizing antibodies in both the circulation and in the lung bronchial fluids. In addition, vaccination elicited strong and long lasting RBD-specific TH1 T cell responses including CD4+ and CD8+ T cell memory. Intramolecular 13C signals inform about carbon allocation in plants Stable isotopes at natural abundance are key tools to study physiological processes occurring outside the temporal scope of manipulation and monitoring experiments. Whole-molecule carbon isotope ratios (13C/12C) enable assessments of plant carbon uptake yet conceal information about carbon allocation. Here, we identify an intramolecular 13C/12C signal at treering glucose C-5 and C-6 and develop experimentally testable theories on its origin. More specifically, we assess the potential of processes within C3 metabolism for signal introduction based (inter alia) on constraints on signal propagation posed by metabolic networks. We propose that the intramolecular signal reports carbon allocation into major metabolic pathways in actively photosynthesising leaf cells including the anaplerotic, shikimate, and non-mevalonate pathway. We support our theoretical framework by linking it to previously reported whole-molecule 13C/12C increases in cellulose of ozone-treated Betula pendula and a highly significant relationship between the intramolecular signal and tropospheric ozone concentration. Our theory postulates a pronounced preference of leaf-cytosolic triose-phosphate isomerase to catalyse the forward reaction in vivo (dihydroxyacetone phosphate to glyceraldehyde 3-phosphate). In conclusion, intramolecular 13C/12C analysis resolves information about carbon uptake and allocation enabling more comprehensive assessments of carbon metabolism than whole-molecule 13C/12C analysis.

HighlightIntramolecular 13C/12C analysis resolves information about carbon uptake and allocation (and associated environmental controls) enabling more comprehensive assessments of carbon metabolism, plant-environment interactions, and environmental variability than whole-molecule 13C/12C analysis. Multivariate Genomewide Association Analysis by Iterative Hard Thresholding 1In genome-wide association studies (GWAS), analyzing multiple correlated traits is potentially superior to conducting multiple univariate analyses. Standard methods for multivariate GWAS operate marker-by-marker and are computationally intensive. We present a penalized regression algorithm for multivariate GWAS based on iterative hard thresholding (IHT) and implement it in a convenient Julia package MendelIHT.jl (https://github.com/OpenMendel/MendelIHT.jl). In simulation studies with up to 100 traits, IHT exhibits similar true positive rates, smaller false positive rates, and faster execution times than GEMMAs linear mixed models and mv-PLINKs canonical correlation analysis. On UK Biobank data, our IHT software completed a 3-trait joint analysis in 20 hours and an 18-trait joint analysis in 53 hours, requiring up to 80GB of computer memory. In short, our software enables geneticists to fit a single regression model that simultaneously considers the effect of all SNPs and dozens of traits. Synthetic Cooling Agents in US-marketed E-cigarette Refill Liquids and Popular Disposable E-cigarettes: Chemical Analysis and Risk Assessment. BackgroundMenthol, through its cooling and pleasant sensory effects, facilitates smoking and tobacco product initiation, resulting in the high popularity of mint/menthol-flavored E-cigarettes. More recently, E-cigarette vendors started marketing synthetic cooling agents as additives that impart a cooling effect but lack a characteristic minty odor. Knowledge about content of synthetic coolants in US-marketed E-cigarette products and associated health risks is limited.

MethodsE-liquid vendor sites were searched with the terms "koolada", "kool/cool", "ice", or WS-3/WS-23, denoting individual cooling agents, and relevant refill E-liquids were purchased. "Ice" flavor varieties of Puffbar, the most popular disposable E-cigarette brand, were compared with non-"Ice" varieties. E-liquids were characterized, and synthetic coolants quantified using GC/MS. Margin of exposure (MOE), a risk assessment parameter, was calculated to assess the risk associated with synthetic coolant exposure from E-cigarette use.

ResultsWS-3 was detected in 24/25 refill E-liquids analyzed. All Puffbar flavor varieties contained either WS-23 (13/14) or WS-3 (5/14), in both "Ice"- and non-"Ice" flavors. Modeling consumption of WS-3 from vaped E-liquids, resulted in MOEs below the safe margin of 100 for most daily use scenarios. MOEs for WS-23 were <100 for 10/13 Puffbar flavors in all use scenarios. Puffbar power specifications are identical to Juul devices.

ConclusionsSynthetic cooling agents (WS-3/WS-23) were present in US-marketed E-cigarettes, at levels that may result in consumer exposures exceeding safety thresholds set by regulatory agencies. Synthetic coolants are not only found in mint-or menthol-flavored products, but also in fruit- and candy-flavored products, including popular disposable E-cigarette products such as Puffbar.

ImplicationsSynthetic cooling agents are widely used in "kool/cool"- and "ice"-flavored E-liquids and in E-liquids without these labels, both as a potential replacement for menthol or to add cooling  notes to non-menthol flavors. These agents may be used to bypass current and future regulatory limits on menthol content in tobacco products, and not just E-cigarettes. Since synthetic cooling agents are odorless, they may not fall under the category of "characterizing flavor", potentially circumventing regulatory measures based on this concept. Regulators need to consider the additional health risks associated with exposure to synthetic cooling agents. Integrative Modelling of Signalling Network Dynamics Identifies Cell Type-selective Therapeutic Strategies for FGFR4-driven Cancers Oncogenic FGFR4 signalling represents a potential therapeutic target in many cancer types, including triple negative breast cancer (TNBC) and hepatocellular carcinoma (HCC). However, resistance to single-agent therapy directed at FGFR4 remains a major challenge, prompting the need to identify more effective combinatorial therapeutic strategies. Here, we integrated computational network modelling and experimental validation to characterise dynamic reprogramming of the FGFR4 signalling network in TNBC following FGFR4 kinase inhibition. We found that AKT, which signals downstream of FGFR4, displayed a rapid and potent reactivation following FGFR4 targeting. Through model-based simulation and systematic prediction of the effect of co-targeting specific network nodes, we predicted, and validated experimentally, strong synergism of co-targeting FGFR4 and particular ErbB kinases or AKT, but not the upstream kinase PI3K. Further, incorporation of protein expression data from hundreds of cancer cell types enabled us to adapt our model to other diverse cellular contexts, leading to the prediction that while AKT rebound occurs frequently, it is not a general phenomenon. Instead, ERK is reactivated in a subset of cell types, including the FGFR4-driven HCC cell line Hep3B. This was subsequently corroborated, and moreover, co-targeting FGFR4 and MEK in Hep3B cells markedly enhanced inhibition of cell proliferation. Overall, these findings provide novel insights into the dynamics of drug-induced network remodelling in cancer cells, highlight the impact of protein expression heterogeneity on network response to targeted therapy and identify candidate cell type-selective combination treatments for FGFR4-driven cancer.