Theory and Simulations of condensin mediated loop extrusion in DNA Condensation of hundreds of mega-base-pair-long human chromosomes in a small nuclear volume is a spectacular biological phenomenon. This process is driven by the formation of chromosome loops. The ATP consuming motor, condensin, interacts with chromatin segments to actively extrude loops. Motivated by real-time imaging of loop extrusion (LE), we created an analytically solvable model, predicting the LE velocity and step size distribution as a function of external load. The theory fits the available experimental data quantitatively, and suggests that condensin must undergo a large conformational change, induced by ATP binding, bringing distant parts of the motor to proximity. Simulations using a simple model confirm that the motor transitions between an open and a closed state in order to extrude loops by a scrunching mechanism, similar to that proposed in DNA bubble formation during bacterial transcription. Changes in the orientation of the motor domains are transmitted over ~ 50 nm, connecting the motor head and the hinge, thus providing an allosteric basis for LE.
Correlative all-optical quantification of mass density and mechanics of sub-cellular compartments with fluorescence specificity Quantitative measurements of physical parameters become increasingly important for understanding biological processes. Brillouin microscopy (BM) has recently emerged as one technique providing the 3D distribution of viscoelastic properties inside biological samples -- so far relying on the implicit assumption that refractive index (RI) and density can be neglected. Here, we present a novel method (FOB microscopy) combining BM with optical diffraction tomography and epi-fluorescence imaging for explicitly measuring the Brillouin shift, RI and absolute density with specificity to fuorescently labeled structures. We show that neglecting the RI and density might lead to erroneous conclusions. Investigating the nucleoplasm of wild-type HeLa cells, we find that it has lower density but higher longitudinal modulus than the cytoplasm. Thus, the longitudinal modulus is not merely sensitive to the water content of the sample -- a postulate vividly discussed in the field. We demonstrate the further utility of FOB on various biological systems including adipocytes and intracellular membraneless compartments. FOB microscopy can provide unexpected scientific discoveries and shed quantitative light on processes such as phase separation and transition inside living cells.
Cell cycle-independent integration of stress signals promotes Non-G1/G0 quiescence entry Cellular quiescence is a non-proliferative state required for cell survival under stress and during development. In most quiescent cells, proliferation is stopped in a reversible state of low Cdk1 kinase activity; in many organisms, however, quiescent states with high Cdk1 activity can also be established through still uncharacterized stress or developmental mechanisms. Here, we used a microfluidics approach coupled to phenotypic classification by machine learning to identify stress pathways associated with starvation-triggered high-Cdk1 quiescent states in Saccharomyces cerevisiae. We found that low- and high-Cdk1 quiescent states shared a core of stress-associated processes, such as autophagy, protein aggregation, and mitochondrial upregulation, but differed in the nuclear accumulation of the stress transcription factors Xbp1, Gln3, and Sfp1. The decision between low- or high-Cdk1 quiescence was controlled by cell cycle-independent accumulation of Xbp1, which acted as a time-delayed integrator of the duration of stress stimuli. Our results show how cell cycle-independent stress-activated factors promote cellular quiescence outside of G1/G0. Bias correction for inverse variance weighting Mendelian randomization BackgroundInverse-variance weighted two-sample Mendelian randomization (IVW-MR) is the most widely used approach that utilizes genome-wide association studies (GWAS) summary statistics to infer the existence and the strength of the causal effect between an exposure and an outcome. Estimates from this approach can be subject to different biases due to the use of weak instruments and winners curse, which can change as a function of the overlap between the exposure and outcome samples.

MethodsWe developed a method (MRlap) that simultaneously considers weak instrument bias and winners curse, while accounting for potential sample overlap. Assuming spike-and-slab genomic architecture and leveraging LD-score regression and other techniques, we could analytically derive, reliably estimate, and hence correct for the bias of IVW-MR using association summary statistics only.

ResultsWe tested our approach using simulated data for a wide range of realistic settings. In all the explored scenarios, our correction reduced the bias, in some situations by as much as 30 folds. When applied to real data on obesity-related exposures, we observed significant differences between IVW-based and corrected effects, both for non-overlapping and fully overlapping samples. Additionally, our results are consistent with the fact that the strength of the biases will decrease as the sample size increases. Finally, we also showed that the overall bias is also dependent on the genetic architecture of the exposure, and traits with low heritability and/or high polygenicity are more strongly affected.

ConclusionsOur method not only reduces bias in causal effect estimation but also enables the use of much larger GWAS sample sizes, by allowing for potentially overlapping samples. Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. Use of lignocellulosic residue from second-generation ethanol production to enhance methane production through co-digestion 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. Genetic coupling of life-history and aerobic performance in Atlantic salmon indicates potential constrains on life-history evolution A better understanding of the genetic and phenotypic architecture underlying life-history variation is a longstanding aim in biology. Theories suggest energy metabolism determines life-history variation by modulating resource acquisition and allocation trade-offs, but the genetic underpinnings of the relationship and its dependence on ecological conditions have rarely been demonstrated. The strong genetic determination of age-at-maturity by two unlinked genomic regions (vgll3 and six6) makes Atlantic salmon (Salmo salar) an ideal model to address these questions. Using more than 250 juveniles in common garden conditions, we quantified the covariation between metabolic phenotypes -standard and maximum metabolic rates (SMR and MMR), and aerobic scope (AS) - and the life-history genomic regions and tested if food availability modulates the relationships. We found that the early maturation genotype in vgll3 was associated with higher MMR and consequently AS. Additionally, MMR exhibited physiological epistasis; it was decreased when late maturation genotypes co-occurred in both genomic regions. Contrary to our expectation, the life-history genotypes had no effects on SMR. Further, food availability had no effect on the genetic covariation, suggesting a lack of genotype-by-environment interactions. Our results provide insights on the key organismal processes that link energy use at the juvenile-stage to age-at-maturity, indicating potential mechanisms by which metabolism and life-history can coevolve. Structure-Functional-Selectivity Relationship Studies on A-86929 Analogs and Small biaryl Fragments toward Discovery of biased D1 agonists Dopamine regulates normal functions such as movement, reinforcement learning, and cognition, and its dysfunction has been implicated in multiple psychiatric and neurological disorders. Dopamine acts through the D1- (D1R and D5R) and D2-class (D2R, D3R and D4R) of seven transmembrane receptors, and activates both G-protein- and {beta}-arrestin-dependent signaling pathways, to mediate its physiological effects. Current dopamine receptor-based therapies are used to ameliorate motor deficits in Parkinsons disease, or as antipsychotic medications for schizophrenia. These drugs show efficacy for ameliorating only some symptoms caused by dopamine dysfunction and are plagued by debilitating side-effects. Studies in primates and rodents have shown that shifting the balance of dopamine receptor signaling towards the arrestin pathway can be beneficial for inducing normal movement, while reducing motor side-effects such as dyskinesias, and can be efficacious at enhancing cognitive function compared to balanced agonists. Several structure-activity-relationship (SAR) studies have embarked on discovering {beta}-arrestin-biased dopamine agonists, focused on D2 partial agonists, non-catechol D1 agonists, and mixed D1/D2R dopamine receptor agonists. Here, we describe an SAR study to identify novel D1R {beta}-arrestin biased ligands using A-86929, a high-affinity D1R catechol agonist, as a core scaffold. Previously described and novel analogs of A-86929 were synthesized and screened in vitro for structure-functional-selectivity relationships (SFSR) studies to identify chemical motifs responsible for {beta}-arrestin biased activity at both D1 and D2Rs. Most of the A-86929 analogs screened were G protein biased but none of them were exclusively arrestin-biased. Additionally, various catechol aryl fragments were designed and synthesized. Other compounds surveyed included hydroxyl and chloro analogs of dopamine to test for hydrogen bonding and ionic interactions. Some of these small molecular probes displayed weak bias towards the {beta}-arrestin pathway. Continued in-depth SFSR studies informed by structure determination, molecular modeling, and mutagenesis studies will facilitate discovery of potent and efficacious arrestin-biased dopamine receptor ligands. Assigning metabolic rate measurements to torpor and euthermy in heterothermic endotherms: 'torpor', a new package for R. Torpor is a state of controlled reduction of metabolic rate (M) in endotherms. Assigning measurements of M to torpor or euthermy can be challenging, especially when the difference between euthermic M and torpid M is small, in species defending a high minimal body temperature in torpor, in thermolabile species, and slightly below the thermoneutral zone (TNZ). Here, we propose a novel method for distinguishing torpor from euthermy. We use the variation in M measured during euthermic rest and torpor at varying ambient temperatures (Ta) to objectively estimate the lower critical temperature (Tlc) of the TNZ and to assign measurements to torpor, euthermic rest or rest within TNZ. In addition, this method allows the prediction of M during euthermic rest and torpor at varying Ta, including resting M within the TNZ. The present method has shown highly satisfactory results using 28 published sets of metabolic data obtained by respirometry on 26 species of mammals. Ultimately, this novel method aims to facilitate analysis of respirometry data in heterothermic endotherms. Finally, the development of the associated R-package (torpor) will enable widespread use of the method amongst biologists.

Summary statementThe presented method and its associated R-package (torpor) enable the assignment of metabolic rate measurements to torpor or euthermy, ultimately improving the standardization of respirometry analyses in heterotherms. 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. Simple and practical sialoglycan encoding system reveals vast diversity in nature and identifies a universal sialoglycan-recognizing probe derived from AB5 toxin B subunits Vertebrate sialic acids (Sias) display much diversity in modifications, linkages and underlying glycans. Slide microarrays allow high-throughput explorations of sialoglycan-protein interactions. A microarray presenting [~]150 structurally-defined sialyltrisaccharides with various Sias linkages and modifications still poses challenges in planning, data sorting, visualization and analysis. To address these issues, we devised a simple 9-digit code for sialyltrisaccharides with terminal Sias and underlying two monosaccharides assigned from the non-reducing end, with three digits assigning a monosaccharide, its modifications, and linkage. Calculations based on the encoding system reveals >113,000 likely linear sialyltrisaccharides in nature. Notably a biantennary N-glycan with two terminal sialyltrisaccharides could thus have >1010 potential combinations and a triantennary N-glycan with three terminal sequences, >1015 potential combinations. While all possibilities likely do not exist in nature, sialoglycans encode enormous diversity. While glycomic approaches are used to probe such diverse sialomes, naturally-occurring bacterial AB5 toxin B subunits are simpler tools to track the dynamic sialome in biological systems. Sialoglycan microarray was utilized to compare sialoglycan-recognizing bacterial toxin B subunits. Unlike the poor correlation between B subunits and species phylogeny, there is stronger correlation with Sia-epitope preferences. Further supporting this pattern, we report a B subunit (YenB) from Yersinia enterocolitica (broad host range) recognizing almost all sialoglycans in the microarray, including 4-O-acetylated-Sias not recognized by a Y. pestis orthologue (YpeB). Differential Sia-binding patterns were also observed with phylogenetically-related B subunits from Escherichia coli (SubB), Salmonella Typhi (PltB), S. Typhimurium (ArtB), extra-intestinal E.coli (EcPltB), Vibrio cholera (CtxB), and cholera family homologue of E. coli (EcxB). Cast away on Mindoro island: could density dependence threaten tamaraw (Bubalus mindorensis) on the long term? Endangered species, despite often living at low population densities, may undergo unexpected density-dependent feedbacks in the case of successful recovery or marked reduction in range. Because density-dependence dynamics can increase risk of extinction, these effects can hamper conservation efforts. In this study, we analyse the dynamics of the largest population of the tamaraw (Bubalus mindorensis), a critically endangered species endemic to the Mindoro island, Philippines. The population is located within a < 3,000 ha area in Mounts Iglit-Baco Natural Park, with limited expansion possibilities. We took advantage of a 21 year time series of tamaraw counts to estimate annual population growth rate and possible density-dependence, accounting for sampling errors in the counts. The tamaraw population has been increasing in size at an average rate of +5% per year, as would be expected given its protected status by law. Population growth showed strong spatial structuring, with a population growth close to +10% in the core area of protection, and a reduction of abundance of - 5% at the periphery of its range, inside the protected area. This range constriction is concerning because our best population dynamic model suggests significant negative density-dependence (Bayes factor = 0.9). The contraction of tamaraw range is likely caused by anthropogenic pressures forcing herbivores to live at relatively high densities in the core zone where protection is most effective. Our study highlights the fact that, despite the continuous population growth over the last two decades, the long-term viability of the Mounts Iglit-Baco Natural Park tamaraw population remains uncertain. A Multi-omics Method Enabled by Sequential Metabolomics and Proteomics for Human Pluripotent Stem Cell-derived Cardiomyocytes Human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) show immense promise for patient-specific disease modeling, cardiotoxicity screening, and regenerative therapy development. However, hPSC-CMs in culture have not recapitulated the structural or functional properties of adult CMs in vivo thus far. To gain global insight into hPSC-CM biology, we established a multi-omics method for analyzing the hPSC-CM metabolome and proteome from the same cell culture, creating multi-dimensional profiles of hPSC-CMs. Specifically, we developed a sequential extraction to capture metabolites and proteins from the same hPSC-CM monolayer cultures, and analyzed these extracts using high-resolution mass spectrometry (MS). Using this method, we annotated 205 metabolites/lipids and 4,008 proteins from 106 cells with high reproducibility. We further integrated the proteome and metabolome measurements to create network profiles of molecular phenotypes for hPSC-CMs. Out of 310 pathways identified using metabolomics and proteomics, 40 pathways were considered significantly overrepresented (FDR-corrected p [&le;] 0.05). Highly populated pathways included those involved in protein synthesis (ribosome, spliceosome), ATP generation (oxidative phosphorylation), and cardiac muscle contraction. This multi-omics method achieves deep coverage of metabolites and proteins, creating a multidimensional view of the hPSC-CM phenotype, which provides a strong technological foundation to advance the understanding of hPSC-CM biology. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. T606-phosphorylation deprives the function of Kaiso as a transcription factor and cancer repressor It is well known that Kaiso protein encoded by ZBTB33 gene is a transcription repressor and that Kaiso-P120ctn interaction increases the shift of Kaiso from the nucleus into the cytoplasm. However, the regulatory mechanisms of Kaiso compartmentalization are far from clear. Here, we reported that AKT1 could phosphorylate 606-threonine residue (T606) within the RSSTIP motif of Kaiso in the cytoplasm. The T606-phosphorylated Kaiso (pT606-Kaiso) could directly bind to 14-3-3 family proteins and the depletion of T606 phosphorylation by T606A mutation abolished most of the Kaiso-14-3-3 binding. In addition, the Kaiso-P120ctn interaction was essential for the pT606-Kaiso accumulation in the cytoplasm. Notably, enforced 14-3-3{sigma} (SFN) overexpression could increase the pT606-Kaiso accumulation in the cytoplasm and de-repress the transcription of Kaiso target gene CDH1. Decreased amounts of both pT606-Kaiso and CDH1 proteins were frequently observed in human gastric cancer tissues relative to paired normal controls. The mRNA levels of 14-3-3{sigma} and Kaiso target gene CDH1 were positively and significantly correlated with each other in bioinformatics analyses using publicly available RNA-seq datasets for human normal tissues (n=11688, r=0.60, p<0.001) in the GTEx project and for cancer cell lines (n=1156, r=0.41, p<0.001) in the CCLE project. Furthermore, Kaiso T606A mutant (unable to be phosphorylated) significantly increased the migration and invasion of cancer cells in vitro as well as boosted the growth of these cells in vivo. In conclusion, Kaiso could be phosphorylated by AKT1 at the T606 and the pT606-Kaiso accumulates in the cytoplasm through binding to 14-3-3/P120ctn that de-represses the expression of Kaiso target gene CDH1 in normal tissues. Decreased Kaiso phosphorylation may contribute to the development of gastrointestinal cancer. The status of Kaiso phosphorylation is a determinant factor for the role of Kaiso in the development of cancer. Organ- and function-specific anatomical organization and bioelectronic modulation of the vagus nerve Afferent and efferent fibers in the vagus nerve travel inside fascicles and form branches to innervate organs and regulate organ functions. The spatial organization of fibers and fascicles, with respect to the functions they mediate and the organs they innervate, is unknown. Accordingly, it is unknown whether such organization can be leveraged by bioelectronic devices for function- and organ-specific neurostimulation. To characterize the functional microscopic anatomy of the vagus, we developed a pipeline consisting of micro-computed tomography imaging of fascicles and quantitative immunohistochemistry of single fibers. We found that, in swine, the vagus contains clusters of fascicles, organized along 2 overlapping functional axes: one axis specific to sensory vs. motor somatic and visceral functions, with clusters increasingly separating in the cephalad direction, and a second axis specific to innervated organs, including larynx, lungs and heart, with increasing separation in the caudad direction. Large myelinated afferents and efferents, and small unmyelinated efferents occupy large areas of specific clusters of fascicles; small unmyelinated afferents are numerous and occupy small areas across every fascicle. To test whether fiber populations can be selectively modulated, we developed a multi-contact cuff electrode accommodating the nerve fascicular structure. Stimulation targeting separate nerve sections evokes compound action potentials from fibers of distinct functional types and elicits differential organ responses, including laryngeal muscle contraction, cough reflex, and changes in breathing and heart rate. Our results indicate that vagus fibers are anatomically organized according to functions they mediate and organs they innervate and can be differentially modulated by spatially selective nerve stimulation. Simpati: patient classifier identifies signature pathways based on similarity networks for the disease prediction BACKGROUNDPathway-based patient classification is a supervised learning task which implies a model learning pathways as features to predict the classes of patients. The counterpart of enrichment tools for the pathway analysis are fundamental methods for clinicians and biomedical scientists. They allow to find signature cellular functions which help to define and annotate a disease phenotype. They provide results which lead human experts to manually classify patients. It is a paradox that pathwaybased classifiers which natively resolve this objective are not strongly developed. They could simulate the human way of thinking, decipher hidden multivariate relationships between the deregulated pathways and the disease phenotype, and provide more information than a probability value. Instead, there are currently only two classifiers of such kind, they require a nontrivial hyperparameter tuning, are difficult to interpret and lack in providing new insights. There is the need of new classifiers which can provide novel perspectives about pathways, be easy to apply with different biological omics and produce new data enabling a further analysis of the patients.

RESULTSWe propose Simpati, an innovative and interpretable patient classifier based on pathway-specific patient similarity networks. The first classifier to adopt ad-hoc novel algorithms for such graph type. It standardizes the biological high-throughput dataset of patients profiles with a propagation algorithm that considers the interconnected nature of the cells molecules for inferring a new activity score. This allows Simpati to classify with dense, sparse, and non-homogenous omic data. Simpati organizes patients molecules in pathways represented by patient similarity networks for being interpretable, handling missing data and preserving the patient privacy. A network represents patients as nodes and a novel similarity measure determines how much every pair act co-ordinately in a pathway. Simpati detects signature biological processes based on how much the topological properties of the related networks separate the patient classes. In this step, it includes a new cohesive subgroup detection algorithm to handle patients not showing the same pathway activity as the other class members. An unknown patient is then classified by a unique recommender system which considers how much is similar to known patients and distant from being an outlier. Simpati outperforms previously published classifiers on five cancer datasets described with two biological omics, classifies well with sparse data, identifies more relevant pathways associated to the patients disease than the competitors and has the lowest computational requirements.

CONCLUSIONSimpati can serve as generic-purpose pathway-based classifier of patient classes. It provides signature pathways to unveil the altered biological mechanisms of a disease phenotype and to classify patients according to the learnt pathway-specific similarities. The signature condition and patient prediction can be deciphered considering the patient similarity networks which must reveal the members of a patient class more cohesive and similar than the non-members. Simpati divides the pathways in up and downinvolved. Upinvolved when the signaling cascades generated by the altered molecules of the disease patients impact stronger the pathway than the ones of the control class. We provide an R implementation, a graphical user interface and a visualization function for the patient similarity networks. The software is available at: https://github.com/LucaGiudice/Simpati Cell-specific Bioorthogonal Tagging of Glycoproteins Altered glycosylation is an undisputed corollary of cancer development. Understanding these alterations is paramount but hampered by limitations underlying cellular model systems. For instance, the intricate interactions between tumour and host cannot be adequately recapitulated in monoculture of tumour-derived cell lines. More complex co-culture models usually rely on sorting procedures for proteome analyses and rarely capture the details of protein glycosylation. Here, we report a strategy termed Bio-Orthogonal Cell line-specific Tagging of Glycoproteins (BOCTAG). Cells are equipped by transfection with an artificial biosynthetic pathway that transforms bioorthogonally tagged sugars into the corresponding nucleotide-sugars. Only transfected cells incorporate bioorthogonal tags into glycoproteins in the presence of non-transfected cells. We employ BOCTAG as an imaging technique and to annotate cell-specific glycosylation sites in mass spectrometry-glycoproteomics. We demonstrate application in co-culture and mouse models, allowing for profiling of the glycoproteome as an important modulator of cellular function. The transposable element-rich genome of the cereal pest Sitophilus oryzae BackgroundThe rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions.

ResultsWe sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 Kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino-acids required for insect development and cuticle biosynthesis.

ConclusionsHere we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes. 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. An epigenomic shift in amygdala marks the transition to maternal behaviors in alloparenting virgin female mice In many species, adults care for young offspring that are not their own, a phenomenon called alloparenting. However, most nonparental adults must be sensitized by repeated or extended exposures to newborns before they will robustly display parental-like behaviors. To capture neurogenomic events underlying the transition to active parental caring behaviors, we analyzed brain gene expression and chromatin profiles of virgin female mice co-housed with mothers during pregnancy and after birth. After an initial display of antagonistic behaviors and a surge of defense-related gene expression, we observed a dramatic shift in the chromatin landscape specifically in amygdala of the pup-exposed virgin females, accompanied by a dampening of anxiety-related gene expression. This epigenetic shift coincided with hypothalamic expression of the oxytocin gene and the emergence of behaviors and gene expression patterns classically associated with maternal care. The results outline a neurogenomic program associated with dramatic behavioral changes and suggest molecular networks relevant to human postpartum mental health. A protein microarray-based in vitro transglutaminase assay platform for epitope mapping and peptide discovery: implication for transglutaminase-mediated immunodominant determination Transglutaminases (TGs) are a family of crosslinking enzymes catalyzing the formation of intra- and inter-molecular glutamine-lysine isopeptide bonds in a calcium dependent manner. Activation of transglutaminases is pathogenically associated with severe human diseases including neurodegenerations, cardiovascular diseases, and autoimmune diseases. Although continuous efforts determining the enzymes substrate preference have been witnessed, a high-throughput assay platform with the "omic" efficiency is still missing for the global identification of substrate-specific TG modification sites. In this study we report a protein microarray-based in vitro TG assay platform for rapid and large-scale (up to 30000 reactions per chip) determination of the glutamine (Q)-bearing TG modification motifs. With this platform we identified the Q16 in superoxide dismutase 1 and Q109 in alpha-synuclein as the primary modification sites for tissue transglutaminase (TG2), the most ubiquitous member of the enzyme family. Of particular interest, we found a close match between the modification motifs and published vaccine epitope sequences in alpha-synuclein, implying an essential and intrinsic role transglutaminase might play in the determination of immunodominant epitopes. Our data collectively suggest the glutamine and its follow-up five residues on the C terminal of a protein compose a minimal determinant motif for TG2 modification and the TG2 modification motifs determined by our platform could finally correspond to the substrates epitope sequences in antigen processing. To screen for site-specific interfering peptides and assist gene editing-based protein engineering, we optimized through onchip amino-acid scanning the TG2 modification motif QQIV in the extracellular matrix protein fibronectin and obtained 14 variants with significantly higher TG2 reactivity that might serve as the competitive inhibitor peptides and 1 with lower reactivity. We further confirmed this optimization approach with 12-mer peptides, the longest ones that could be synthesized on the chip. Taken together, our synthetic transglutaminase assay platform might be able to deliver a precise epitope blueprint for immunotherapeutic targeting and provide proof-of-concept and directional studies for TG-based peptide discovery and protein design. Evaluating the Anti-aging Effects of Chinese Herbal Medicine Rhodiola rosea L. in Cultured Keratinocytes Among the many health benefits Chinese herbal medicine (CHM) presents, anti-aging is of special interest (Shen, Jiang, Yang, Wang, & Zhu, 2017; Zhao & Luo, 2017). Reported to possess anti-aging effects, the CHM Rhodiola rosea L., known colloquially as the "golden root", has been widely incorporated in various drinks, daily supplements, and even cosmetics. This study investigates the effects of commercial R. rosea extracts and natural R. rosea roots on preventing UV-induced photoaging of the skin and correlates such effects with the composition of active ingredients in the extracts. To simulate the photoaging process, drug treated HaCaT cells were exposed to UVA and UVB radiation. The pharmacological anti-aging effects of R. rosea extracts were evaluated qualitatively and quantitatively through confocal immunofluorescence images with {gamma}-H2A.X marker and telomerase activity assay (Telo TAGGG Telomerase PCR-ELISA). Preparatory thin layer chromatography and high-performance liquid chromatography were performed to isolate and quantify active ingredients. Cultured HaCaT cells showed morphological change after exposure to both UVA (>15.0 J) and UVB (>2000 mJ). The photoaging of keratinocytes was rescued by pretreating cells with R. rosea extracts as well as salidroside and rosavin active ingredients (P < 0.05). R. rosea-treated cells were characterized by increased telomerase activity and fewer {gamma}-H2A.X foci compared to that of the control. Extracts with better preventative effects contained higher salidroside and rosavin content. The findings in this study reaffirm R. roseas efficacy as an anti-aging remedy and provide a basis for CHMs integration into the mainstream of global healthcare. Fatal COVID-19 and non-COVID-19 Acute Respiratory Distress Syndrome is Associated with Incomplete Alveolar Type 1 Epithelial Cell Differentiation from the Transitional State Without Fibrosis ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. Species-specific developmental timing dictates expansion of the avian wing skeletal pattern A fundamental question in biology is how species-specific rates of embryonic development are controlled in differently sized species. To address this problem, we compared wing development in the quail and the larger chick. We reveal that pattern formation (humerus to digits) is faster in the quail wing as determined by the earlier activation of 5 Hox genes, termination of developmental organisers (Shh and Fgf8) and the laying down of the skeleton (Sox9). Using interspecies tissue grafts, we show that developmental timing can be reset during a critical window in which retinoic acid signalling is active. Accordingly, retinoic acid can switch developmental timing between the quail and chick, and also between the chick and the larger turkey. We show that the incremental growth rate is equivalent between all three species, suggesting that differences in the expansion of the skeletal pattern are governed primarily by the tempo of development. Our findings have implications for how development is timed both within and between different species. 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.