Feline Papillomavirus Isolate P20 Assembled from Metagenomic Data Isolated from the Human Skin of a House Cat Owner A feline papillomavirus was assembled from metagenomic sequencing data collected from the human skin of a house cat owner. This circular papillomavirus strain P20 is 8069 bp in length, has a GC content of 54.38%, and displays genome organization typical of feline papillomaviruses with six annotated protein coding regions. The genome exhibits approximately 75% sequence similarity to other feline papillomavirus genomes.
Identification of a novel cationic glycolipid in Streptococcus agalactiae that contributes to brain entry and meningitis Bacterial membrane lipids are critical for membrane bilayer formation, cell division, protein localization, stress responses, and pathogenesis. Despite their critical roles, membrane lipids have not been fully elucidated for many pathogens. Here, we report the discovery of a novel cationic glycolipid, Lysyl-Glucosyl-Diacylglycerol (Lys-Glc-DAG) that is synthesized in high abundance by the bacterium Streptococcus agalactiae (Group B Streptococcus, GBS). To our knowledge, Lys-Glc-DAG is more positively charged than any other known lipids. Lys-Glc-DAG carries two positive net charges per molecule, distinct from the widely described lysylated phospholipid Lysyl-phosphatidylglycerol (Lys-PG) which carries one positive net charge due to the presence of a negatively charged phosphate moiety. We use normal phase liquid chromatography (NPLC) coupled with electrospray ionization (ESI) high-resolution tandem mass spectrometry (HRMS/MS) and genetic approaches to determine that Lys-Glc-DAG is synthesized by the enzyme MprF in GBS, which covalently modifies the neutral glycolipid Glc-DAG with the cationic amino acid lysine. GBS is a leading cause of neonatal meningitis, which requires traversal of the endothelial blood-brain barrier (BBB). We demonstrate that GBS strains lacking mprF exhibit a significant decrease in the ability to invade BBB endothelial cells. Further, mice challenged with a GBS{Delta}mprF mutant developed bacteremia comparably to Wild-Type infected mice yet had less recovered bacteria from brain tissue and a lower incidence of meningitis. Thus, our data suggest that Lys-Glc-DAG may contribute to bacterial uptake into host cells and disease progression. Importantly, our discovery provides a platform for further study of cationic lipids at the host-pathogen interface.
Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Single-Cell Multiomics Reveals Distinct Cell States at the Top of the Human Hematopoietic Hierarchy The advent of single cell (Sc) genomics has challenged the dogma of haematopoiesis as a tree-like structure of stepwise lineage commitment through distinct and increasingly restricted progenitor populations. Instead, analysis of ScRNA-seq has proposed that the earliest events in human hematopoietic stem cell (HSC) differentiation are characterized by only subtle molecular changes, with hematopoietic stem and progenitor cells (HSPCs) existing as a continuum of low-primed cell-states that gradually transition into a specific lineage (CLOUD-HSPCs). Here, we combine ScRNA-seq, ScATAC-seq and cell surface proteomics to dissect the heterogeneity of CLOUD-HSPCs at different stages of human life. Within CLOUD-HSPCs, pseudotime ordering of both mRNA and chromatin data revealed a bifurcation of megakaryocyte/erythroid and lympho/myeloid trajectories immediately downstream a subpopulation with an HSC-specific enhancer signature. Importantly, both HSCs and lineage-restricted progenitor populations could be prospectively isolated based on correlation of their molecular signatures with CD35 and CD11A expression, respectively. Moreover, we describe the changes that occur in this heterogeneity as hematopoiesis develops from neonatal to aged bone marrow, including an increase of HSCs and depletion of lympho-myeloid biased MPPs. Thus, this study dissects the heterogeneity of human CLOUD-HSPCs revealing distinct HSPC-states of relevance in homeostatic settings such as ageing. Stimulated oxygen transport in tissue by magnetic needle acupuncture AimTo understand physiology of magnetic needle acupuncture by investigating O2 transport in tissue during needle intervention.

MethodsO2 transport in tissue is modeled by utilizing COMSOL with magnetic needle inserted into muscle tissue in a 2D porous media. The damaged tissue has been mimicked by an extracted tissue block with 1st order O2 consumption rate. The convection-diffusion O2 transport in the damaged tissue has been further evaluated by varying magnetic flux density B0 of the needle (0-1 T), myoglobin concentration (0-1 mM), O2 tension (5-100 Torr), O2 consumption rates, tissue permeability (10-12-10-6 m2) and porosity (0.1-0.9).

Results1) Active O2 transport carried by interstitial flow is enhanced with the intervention of a magnetic needle by generating a high gradient magnetic field around the tip, which exerts a strong force (104 N/m3) on the diamagnetic interstitial fluid to accelerate the flow. 2) This interstitial flow can reach 30 m/s at B0 = 1T and strongly correlates to the needle B0 and tissue permeability. 3) The needle stirs the interstitial flow can pump O2 flux by 1-2 orders of magnitude compared to that without magnetic field. 4) The enhancement of active O2 transport by magnetic needle is site-specific to the tissue in the vicinity of the tip. This enhancement is more effective in edema condition with a high tissue permeability (>10-9 m2).

ConclusionsThe dramatic enhancement of O2 transport to restore the O2 mitochondria metabolism for dysfunctional muscle tissues is the fundamental physiological mechanism of magnetic needle acupuncture. Citation needed?Wikipedia bibliometrics during the first waveof the COVID-19 pandemic. BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

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

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. FGFR2 fusion protein-driven mouse models of intrahepatic cholangiocarcinoma unveil a necessary role for Erk signaling Background and aimsAbout 15% of intrahepatic cholangiocarcinoma (iCCA) express fibroblast growth factor receptor 2 (FGFR2) fusion proteins (FFs), most often in concert with mutationally inactivated TP53, CDKN2A or BAP1. FFs span residues 1-768 of FGFR2 fused to sequences encoded by any of a long list (>60) of partner genes, a configuration sufficient to ignite oncogenic FF activation. In line, FGFR-specific tyrosine kinase inhibitors (F-TKI) were shown to provide clinical benefit in FF+ iCCA, although responses were partial and/or limited by resistance mechanisms, including the FF V565F gatekeeper mutation. Herein we present an FF-driven murine iCCA model and exploit its potential for pre-clinical studies on FF therapeutic targeting.

MethodsFour iCCA FFs carrying different fusion sequences were expressed in Tp53-/- mouse liver organoids. Tumorigenic properties of genetically modified liver organoids were assessed by intrahepatic/subcutaneous transplantation in immuno-deficient mice. Cellular models derived from neoplastic lesions were exploited for pre-clinical studies.

ResultsTumors diagnosed as CCA were obtained upon transplantation of FF-expressing liver organoids. The penetrance of this tumorigenic phenotype was influenced by FF identity. Tumor organoids and 2D cell lines derived from CCA lesions were addicted to FF signaling via Ras-Erk, regardless of FF identity or presence of V565F mutation. Double blockade of FF-Ras-Erk pathway by concomitant pharmacological inhibition of FFs and Mek1/2 provided greater therapeutic efficacy than single agent F-TKI in vitro and in vivo.

ConclusionsFF-driven iCCA pathogenesis was successfully modelled in murine Tp53-/- background. This model revealed biological heterogeneity among structurally different FFs. Double blockade of FF-Erk signaling deserves consideration for improving precision-based approaches against human FF+ iCCA.

Abbreviations used in this paper: ANOVA, analysis of variance; Bap1, BRCA1-Associated-Protein 1; Cdkn2a, cyclin-dependent kinase inhibitor 2A; Cftr, cystic fibrosis transmembrane conductance regulator; Ck19, cytokeratin 19; Cyp3A, cytochrome P450, family 3, subfamily A; EGF, Epidermal growth factor; EGFR, Epidermal growth factor receptor; EpCAM, epithelial cell adhesion molecule; Erk, extracellular signal-regulated kinase; FGFR2, fibroblast growth factor receptor 2; FRS2, fibroblast growth factor receptor substrate 2; GRB2, growth factor receptor-bound 2; GSEA, gene set enrichment analysis; GSVA, gene set variation analysis; H&E, hematoxylin and eosin; HepPar1, hepatocyte Paraffin 1; Hnf4-7, hepatocyte nuclear factor 4 alpha; Hprt, hypoxanthine-guanine phosphoribosyl transferase; LGR5, leucine-rich repeat-containing G-protein coupled receptor 5; NTRK, neurotrophic Tyrosine Kinase; Parp, Poly (ADP-ribose) polymerase; RECIST, response evaluation criteria in solid tumors; SHP2, Src homology phosphatase 2; Ttr, transthyretin. 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. OrganoID: a versatile deep learning platform for tracking and analysis of single-organoid dynamics Organoids have immense potential as ex vivo disease models for drug discovery and personalized drug screening. Dynamic changes in individual organoid morphology, number, and size can indicate important drug responses, however these metrics are difficult and labor-intensive to obtain for high-throughput image datasets. Here, we present OrganoID, a robust image analysis platform that automatically recognizes, labels, and tracks single organoids, pixel-by-pixel, in brightfield and phase-contrast microscopy experiments. The platform was trained on images of pancreatic cancer organoids and validated on separate images of pancreatic, lung, colon, and adenoid cystic carcinoma organoids, which showed excellent agreement with manual measurements of organoid count (96%) and size (95%) without any parameter adjustments. Single-organoid tracking accuracy remained above 89% over a four-day time-lapse microscopy study. Automated single-organoid morphology analysis of a chemotherapy dose-response experiment identified decreased organoid circularity as an important morphological feature reflecting drug response. OrganoID enables straightforward, detailed, and accurate image analysis to accelerate the use of organoids in high-throughput, data-intensive biomedical applications. Hidden dynamics of proteasome autoregulation discovered by cryo-EM data-driven deep learning The cellular functions are executed by macromolecular complexes in nonequilibrium dynamic processes, which exhibit a vast diversity of conformational states. Solving conformational continuum of important complexes at atomic level remains a key challenge in the methodology of computational biophysics. Here we introduce a deep-learning framework, named AlphaCryo4D, which enables atomic-level cryogenic electron microscopy reconstructions of nonequilibrium conformational continuum. AlphaCryo4D integrates 3D deep residual learning with manifold embedding of energy landscapes, which simultaneously improves 3D classification accuracy and reconstruction resolution via an energy-based particle-voting algorithm. Applications of this methodology to six experimental datasets demonstrated its generality in breaking resolution barrier of visualizing dynamic components of functional complexes, in choreographing continuous inter-subunit motions and in exploring their  hidden conformational space. Importantly, AlphaCryo4D identified 64 conformers of the substrate-bound human 26S proteasome and reconstituted single-nucleotide-exchange dynamics of proteasomal AAA-ATPase motor during translocation initiation, which together illuminates a grand hierarchical allostery for proteasome autoregulation. Higher-order epistasis and phenotypic prediction Contemporary high-throughput mutagenesis experiments are providing an increasingly detailed view of the complex patterns of genetic interaction that occur between multiple mutations within a single protein or regulatory element. By simultaneously measuring the effects of thousands of combinations of mutations, these experiments have revealed that the genotype-phenotype relationship typically reflects genetic interactions not only between pairs of sites, but also higher-order interactions among larger numbers of sites. However, modeling and understanding these higher-order interactions remains challenging. Here, we present a method for reconstructing sequence-to-function mappings from partially observed data that can accommodate all orders of genetic interaction. The main idea is to make predictions for unobserved genotypes that match the type and extent of epistasis found in the observed data. This information on the type and extent of epistasis can be extracted by considering how phenotypic correlations change as a function of mutational distance, which is equivalent to estimating the fraction of phenotypic variance due to each order of genetic interaction (additive, pairwise, three-way, etc.). Using these estimated variance components, we then define an empirical Bayes prior that in expectation matches the observed pattern of epistasis, and reconstruct the genotype-phenotype mapping by conducting Gaussian process regression under this prior. To demonstrate the power of this approach, we present an application to the antibody-binding domain GB1 and also provide a detailed exploration of a dataset consisting of high-throughput measurements for the splicing efficiency of human pre-mRNA 5' splice sites, for which we also validate our model predictions via additional low-throughput experiments. Possible Cross Reactivity of Feline and White-tailed Deer Antibodies Against the SARS-CoV-2 Receptor Binding Domain In late 2019, a novel coronavirus began circulating within humans in central China. It was designated SARS-CoV-2 because of its genetic similarities to the 2003 SARS coronavirus (SARS-CoV). Now that SARS-CoV-2 has spread worldwide, there is a risk of it establishing new animal reservoirs and recombination with native circulating coronaviruses. To screen local animal populations in the United States for exposure to SARS-like coronaviruses, we developed a serological assay using the receptor binding domain (RBD) from SARS-CoV-2. SARS-CoV-2s RBD is antigenically distinct from common human and animal coronaviruses allowing us to identify animals previously infected with SARS-CoV or SARS-CoV-2. Using an indirect ELISA for SARS-CoV-2s RBD, we screened serum from wild and domestic animals for the presence of antibodies against SARS-CoV-2s RBD. Surprisingly pre-pandemic feline serum samples submitted to the University of Tennessee Veterinary Hospital were [~]50% positive for anti-SARS RBD antibodies. Some of these samples were serologically negative for feline coronavirus (FCoV), raising the question of the etiological agent generating anti-SARS-CoV-2 RBD cross-reactivity. We also identified several white-tailed deer from South Carolina with anti-SARS-CoV-2 antibodies. These results are intriguing as cross-reactive antibodies towards SARS-CoV-2 RBD have not been reported to date. The etiological agent responsible for seropositivity was not readily apparent, but finding seropositive cats prior to the current SARS-CoV-2 pandemic highlights our lack of information about circulating coronaviruses in other species.

ImportanceWe report cross-reactive antibodies from pre-pandemic cats and post-pandemic South Carolina white-tailed deer that are specific for that SARS-CoV RBD. There are several potential explanations for this cross-reactivity, each with important implications to coronavirus disease surveillance. Perhaps the most intriguing possibility is the existence and transmission of an etiological agent (such as another coronavirus) with similarity to SARS-CoV-2s RBD region. However, we lack conclusive evidence of pre-pandemic transmission of a SARS-like virus. Our findings provide impetus for the adoption of a One Health Initiative focusing on infectious disease surveillance of multiple animal species to predict the next zoonotic transmission to humans and future pandemics. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. Single-cell Transcriptomic Profiling Unveils Cardiac Cell-type Specific Response to Maternal Hyperglycemia Underlying the Risk of Congenital Heart Defects Congenital heart disease (CHD) is the most prevalent structural malformations of the heart affecting [~]1% of live births. To date, both damaging genetic variations and adverse environmental exposure such as maternal diabetes have been found to cause CHD. Clinical studies show [~]fivefold higher risk of CHD in the offspring of mothers with pregestational diabetes. Maternal pregestational diabetes affects the gene regulatory networks key to proper cardiac development in the fetus. However, the cell-type specificity of these gene regulatory responses to maternal diabetes and their association with the observed cardiac defects in the fetuses remains unknown. To uncover the transcriptional responses to maternal diabetes in the early embryonic heart, we used an established murine model of pregestational diabetes. In this model, we have previously demonstrated an increased incidence of CHD. Here, we show maternal hyperglycemia (matHG) elicits diverse cellular responses during heart development by single-cell RNA-sequencing in embryonic hearts exposed to control and matHG environment. Through differential gene-expression and pseudotime trajectory analyses of this data, we identified changes in lineage specifying transcription factors, predominantly affecting Isl1+ second heart field progenitors and Tnnt2+ cardiomyocytes with matHG. Using in vivo cell-lineage tracing studies, we confirmed that matHG exposure leads to impaired second heart field-derived cardiomyocyte differentiation. Finally, this work identifies matHG-mediated transcriptional determinants in cardiac cell lineages elevate CHD risk and show perturbations in Isl1-dependent gene-regulatory network (Isl1-GRN) affect cardiomyocyte differentiation. Functional analysis of this GRN in cardiac progenitor cells will provide further mechanistic insights into matHG-induced severity of CHD associated with diabetic pregnancies. An exon-intron split framework to prioritize miRNA-driven post-transcriptional regulatory signals and its application to study energy homeostasis in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Bayesian Additive Regression Trees for Genotype by Environment Interaction Models We propose a new class of models for the estimation of genotype by environment (GxE) interactions in plant-based genetics. Our approach, named AMBARTI, uses semi-parametric Bayesian additive regression trees to accurately capture marginal genotypic and environment effects along with their interaction in a fully Bayesian model. We demonstrate that our approach is competitive or superior to similar models widely used in the literature via both simulation and a real world data set. Furthermore, we introduce new types of visualisation to properly assess both the marginal and interactive predictions from the model. An R package that implements our approach is available at https://github.com/ebprado/ambarti. 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. Artificial reefs geographical location matters more than shape, age and depth for sessile invertebrate colonization in the Gulf of Lion (NorthWestern Mediterranean Sea) AO_SCPLOWBSTRACTC_SCPLOWArtificial reefs (ARs) have been used to support fishing activities. Sessile invertebrates are essential components of trophic networks within ARs, supporting fish productivity. However, colonization by sessile invertebrates is possible only after effective larval dispersal from source populations, usually in natural habitat. While most studies focused on short term colonization by pioneer species, we propose to test the relevance of geographic location, shape, age and depth of immersion on the ARs long term colonization by species found in natural stable communities in the Gulf of Lion. We recorded the presence of five sessile invertebrates species, with contrasting life history traits and regional distribution in the natural rocky habitat, on ARs with different shapes deployed during two immersion time periods (1985 and the 2000s) and in two depth ranges (<20m and >20m). At the local level ([~]5kms), neither shape, depth nor immersion duration differentiated ARs assemblages. At the regional scale (>30kms), colonization patterns differed between species, resulting in diverse assemblages. This study highlights the primacy of geographical positioning over shape, immersion duration and depth in ARs colonization, suggesting it should be accounted for in maritime spatial planning. First detection of herpesvirus and prevalence of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Gradient-mixing LEGO robots for purifying DNA origami nanostructures of multiple components by rate-zonal centrifugation DNA origami purification is critical in emerging applications of functionalized DNA nanostructures from basic fundamental biophysics, nanorobots to therapeutics. Advances in DNA origami purification have led to the establishment of rate-zonal centrifugation (RZC) as a scalable, high-yield, and contamination-free approach to purifying DNA origami nanostructures. In RZC purification, a linear density gradient is created using viscous agents, such as glycerol and sucrose, to separate molecules based on their mass and shape during high-rpm centrifugation. However, current methods for creating density gradients are typically time-consuming because of their reliance on slow passive diffusion. Here, we built a LEGO gradient mixer to rapidly create a quasi-continuous density gradient with minimal layering of concentrations using simple rotational motion. We found that rotating two layers of different concentrations at an angle can reduce the diffusion time from a few hours to mere minutes. The instrument needed to perform the movement can be constructed from low-cost components, such as Arduino and LEGO Mindstorms pieces, and has comparable efficacy to commercial gradient mixers currently available. Our results demonstrate that the creation of a linear density gradient can be achieved with minimal labor, time, and cost with this machine. With the recent advances in DNA origami production, we anticipate our findings to further improve the viability of scaling up DNA origami purification in grams quantities. Our simple process enables automated large-scale purification of functionalized DNA origami more feasible in resource-constrained settings. From GWAS Variant to Function: a Study of ~148,000 Variants for Blood Cell Traits. Genome-wide association studies (GWAS) have identified hundreds of thousands of genetic variants associated with complex diseases and traits. However, most variants are noncoding and not clearly linked to genes, making it challenging to interpret these GWAS signals. We present a systematic variant-to-function study, prioritizing the most likely functional elements of the genome for experimental follow-up, for >148,000 variants identified for hematological traits. Specifically, we developed VAMPIRE: Variant Annotation Method Pointing to Interesting Regulatory Effects, an interactive web application implemented in R Shiny (http://shiny.bios.unc.edu/vampire/). This tool efficiently integrates and displays information from multiple complementary sources, including epigenomic signatures from blood cell relevant tissues or cells, functional and conservation summary scores, variant impact on protein and gene expression, chromatin conformation information, as well as publicly available GWAS and phenome-wide association study (PheWAS) results. Leveraging data generated from independently performed functional validation experiments, we demonstrate that our prioritized variants, genes, or variant-gene links are significantly more likely to be experimentally validated. This study not only has important implications for systematic and efficient revelation of functional mechanisms underlying GWAS variants for hematological traits, but also provides a prototype that can be adapted to many other complex traits, paving the path for efficient variant to function (V2F) analyses. 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. Conformist social learning leads to self-organised prevention against adverse bias in risky decision making Given the ubiquity of potentially adverse behavioural bias owing to myopic trial-and-error learning, it seems paradoxical that improvements in decision-making performance through conformist social learning, a process widely considered to be bias amplification, still prevail in animal collective behaviour. Here we show, through model analyses and large-scale interactive behavioural experiments with 585 human subjects, that conformist influence can indeed promote favourable risk taking in repeated experience-based decision making, even though many individuals are systematically biased towards adverse risk aversion. Although strong positive feedback conferred by copying the majoritys behaviour could result in unfavourable informational cascades, our differential equation model of collective behavioural dynamics identified a key role for increasing exploration by negative feedback arising when a weak minority influence undermines the inherent behavioural bias. This  collective behavioural rescue, emerging through coordination of positive and negative feedback, highlights a benefit of collective learning in a broader range of environmental conditions than previously assumed and resolves the ostensible paradox of adaptive collective behavioural flexibility under conformist influences. Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector. Pupil diameter tracked during motor adaptation in humans Pupil diameter, under constant illumination, is known to reflect individuals internal states, such as surprise about observation and environmental uncertainty. Despite the growing use of pupillometry in cognitive learning studies as an additional measure for examining internal states, few studies have used pupillometry in human motor learning studies. Here we provide the first detailed characterization of pupil diameter changes in a short-term reach adaptation paradigm. We measured pupil changes in 91 human participants while they adapted to abrupt, gradual, or switching force field conditions. Sudden increases in movement error caused by the introduction/reversal of the force field resulted in strong phasic pupil dilation during movement accompanied by a transient increase in tonic pre-movement baseline pupil diameter in subsequent trials. In contrast, clear changes in pupil responses were absent when the force field was gradually introduced, indicating that error drove the changes in pupil responses. Nevertheless, we found an association between baseline pupil diameter and awareness of the gradually-introduced perturbation assessed post-experimentally. In all experiments, we also found a strong co-occurrence of larger baseline pupil diameter and slower reaction and movement time after each set break. Interestingly, error-induced pupil responses gradually became insensitive after experiencing multiple reversals. Collectively, these results suggest that tonic baseline pupil diameter reflects ones belief about environmental uncertainty, whereas phasic pupil dilation during movement reflects surprise about a sensory outcome (i.e., movement error), and both effects are modulated by novelty. Our results provide a new approach for non-verbally assessing participants internal states during motor learning. Phytoene synthase 2 in tomato fruits remains functional and contributes to abscisic acid formation In ripening tomato fruits, the leaf-specific carotenoids biosynthesis mediated by phytoene synthase 2 (PSY2) is replaced by a fruit-specific pathway by the expression of two chromoplast-specific genes: phytoene synthase 1 (PSY1) and lycopene-{beta}-cyclase (CYCB). Though both PSY1 and PSY2 genes express in tomato fruits, the functional role of PSY2 is not known. To decipher whether PSY2-mediated carotenogenesis operates in ripening fruits, we blocked the in vivo activity of lycopene-{beta}-cyclases in fruits of several carotenoids and ripening mutants by CPTA (2-(4-Chlorophenylthio)triethylamine hydrochloride), an inhibitor of lycopene-{beta}-cyclases. The CPTA-treatment induced accumulation of lycopene in leaves, immature-green and ripening fruits. Even, in psy1 mutants V7 and r that are deficient in fruit-specific carotenoid biosynthesis, CPTA triggered lycopene accumulation but lowered the abscisic acid level. Differing from fruit-specific carotenogenesis, CPTA-treated V7 and r mutant fruits accumulated lycopene but not phytoene and phytofluene. The lack of phytoene and phytofluene accumulation was reminiscent of PSY2-mediated leaf-like carotenogenesis, where phytoene and phytofluene accumulation is never seen. The lycopene accumulation was associated with the partial transformation of chloroplasts to chromoplasts bearing thread-like crystalline structures. Our study uncovers the operation of a parallel carotenogenesis pathway mediated by PSY2 that provides precursors for abscisic acid biosynthesis in ripening tomato fruits.

Significance statementIt is believed that in ripening tomato fruits phytoene synthase 2 that drives carotenoid biosynthesis in leaves is redundant. Contrary to this, we show that in phytoene synthase 1 mutant fruit that is bereft of lycopene, the chemical inhibition of lycopene {beta}-cyclases triggers lycopene accumulation. Our results uncover that phytoene synthase 2 remains functional in ripening fruits and provides precursors for abscisic acid formation in fruits.