Voxel-wise and spatial modelling of binary lesion masks: Comparison of methods with a realistic simulation framework ObjectivesWhite matter lesions are a very common finding on MRI in older adults and their presence increases the risk of stroke and dementia. Accurate and computationally efficient modelling methods are necessary to map the association of lesion incidence with risk factors, such as hypertension. However, there is no consensus in the brain mapping literature whether a voxel-wise modelling approach is better for binary lesion data than a more computationally intensive spatial modelling approach that accounts for voxel dependence.

MethodsWe review three regression approaches for modelling binary lesion masks including massunivariate probit regression modelling with either maximum likelihood estimates, or mean bias-reduced estimates, and spatial Bayesian modelling, where the regression coefficients have a conditional autoregressive model prior to account for local spatial dependence. We design a novel simulation framework of artificial lesion maps to compare the three alternative lesion mapping methods. The age effect on lesion probability estimated from a reference data set (13,680 individuals from the UK Biobank) is used to simulate a realistic voxel-wise distribution of lesions across age. To mimic the real features of lesion masks, we propose matching brain lesion summaries (total lesion volume, average lesion size and lesion count) across the reference data set and the simulated data sets. Thus, we allow for a fair comparison between the modelling approaches, under a realistic simulation setting.

ResultsOur findings suggest that bias-reduced estimates for voxel-wise binary-response generalized linear models (GLMs) overcome the drawbacks of infinite and biased maximum likelihood estimates and scale well for large data sets because voxel-wise estimation can be performed in parallel across voxels. Contrary to the assumption of spatial dependence being key in lesion mapping, our results show that voxel-wise bias-reduction and spatial modelling result in largely similar estimates.

ConclusionBias-reduced estimates for voxel-wise GLMs are not only accurate but also computationally efficient, which will become increasingly important as more biobank-scale neuroimaging data sets become available.
Cortical responses to vagus nerve stimulation are modulated by brain state in non-human primates Vagus nerve stimulation (VNS) has been tested as therapy for several brain disorders and as a means to modulate cortical excitability and brain plasticity. Cortical effects of VNS, manifesting as vagal-evoked potentials (VEPs), are thought to arise from activation of ascending cholinergic and noradrenergic systems. However, it is unknown whether those effects are modulated by brain state at the time of stimulation. In 2 freely behaving macaque monkeys, we delivered trains of left cervical VNS at different frequencies (5-300 Hz) while recording local field potentials (LFPs) from sites in contralateral prefrontal, sensorimotor and parietal cortical areas. Brain states were inferred from spectral components of LFPs and the presence of overt movement: active awake, resting awake, REM sleep and NREM sleep. VNS elicited VEPs comprising early (<70 ms), intermediate (70-250 ms) and late (>250 ms) components in all sampled cortical areas. The magnitude of only the intermediate and late components was modulated by brain state and pulsing frequency. These findings have implications for the role of ongoing cortical activity and brain state in shaping cortical responses to peripheral stimuli, for the modulation of vagal interoceptive signaling by cortical states, and for the calibration of VNS therapies.
Endothelial cell signature in muscle stem cells validated by VEGFA-FLT1-AKT1 axis promoting survival of muscle stem cell Endothelial and skeletal muscle lineages arise from common embryonic progenitors. Despite their shared developmental origin, adult endothelial cells (ECs) and muscle stem cells (MuSCs) (satellite cells) have been thought to possess distinct gene signatures and signaling pathways. Here we shift this paradigm by uncovering how adult MuSC behavior is affected by the expression of a subset of EC transcripts. We used several computational analyses including single-cell RNAseq to show that MuSCs express low levels of canonical EC markers. We demonstrate that MuSC survival is regulated by one such prototypic endothelial signaling pathway (VEGFA-FLT1). Using pharmacological and genetic gain- and loss-of-function studies, we identify the FLT1-AKT1 axis as the key effector underlying VEGFA-mediated regulation of MuSC survival. All together, our data support that the VEGFA-FLT1-AKT1 pathway promotes MuSC survival during muscle regeneration, and highlights how the minor expression of select transcripts is sufficient for affecting cell behavior. Nuclear and plastid phylogenomic analyses provide insights into the reticulate evolution, species delimitation and biogeography of the Sino-Japanese disjunctive Diabelia (Caprifoliaceae) Understanding biological diversity and the mechanisms of the Sino-Japanese disjunctions are major challenges in eastern Asia biogeography. The Sino-Japanese flora has been broadly studied as an ideal model for plant phylogeography. Diabelia (Caprifoliaceae) is an East Asian genus, with a disjunctive distribution across the Sino-Japanese region. However, relationships within Diabelia remain elusive. In this study, we reconstructed the phylogeny of Diabelia and inferred historical biogeography and evolutionary patterns based on nuclear and plastid sequences from target enrichment and genome skimming approaches, respectively. We found that the main clades within Diabelia were discordant between nuclear and plastid trees. Both nuclear and plastid phylogenetic analyses supported five main clades: D. serrata, D. tetrasepala, D. sanguinea, D. spathulata var. stenophylla and D. spathulata var. spathulata. Species network analyses revealed that Diabelia tetrasepala is likely the result of a hybridization event. Divergence time estimation and ancestral area reconstructions showed that Diabelia originated in Japan during the early Miocene, with subsequent vicariance and dispersal events between Japan and Korea, and between Japan and China. Overall, our results support the division of Diabelia into five main clades and the recognition of five species in the genus. This research provides new insights in the species delimitation and speciation processes of taxonomically complex lineages such as Diabelia. LADON, a natural antisense transcript of NODAL, promotes an invasive behaviour in melanoma cells. The TGF{beta} family member NODAL, known for its role during embryonic development, has also been associated with tumor progression in several cancers. Some of the evidence supporting its involvement in melanoma appeared contradictory, suggesting that NODAL in this context might rely on a non-canonical signalling mode. We found that NODAL inactivation in a metastatic melanoma cell line prevents the cells from acquiring invasive behaviour. However, we show that this phenotype does not result from the absence of NODAL, but from a defect in the expression of a natural antisense transcript of NODAL, here called LADON. We found that LADON promotes the mesenchymal to amoeboid transition that is critical to melanoma cell invasiveness, and that a WNT/{beta}-CATENIN signalling-dependent increase in LADON expression is required to complete this transition. LADONs downstream effectors include, among others, the proto-oncogene MYCN. These results identify LADON as a player in the regulatory network that governs tumor progression in melanoma, and possibly in other types of cancer. SERCA function is impaired in skeletal and cardiac muscles from young DBA/2J mdx mice The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. Estimating cell-type-specific gene co-expression networks from bulk gene expression data with an application to Alzheimer's disease Inferring and characterizing gene co-expression networks has led to important insights on the molecular mechanisms of complex diseases. Most co-expression analyses to date have been performed on gene expression data collected from bulk tissues with different cell type compositions across samples. As a result, the co-expression estimates only offer an aggregate view of the underlying gene regulations and can be confounded by heterogeneity in cell type compositions, failing to reveal gene coordination that may be distinct across different cell types. In this paper, we describe a flexible framework for estimating cell-type-specific gene co-expression networks from bulk sample data, without making specific assumptions on the distributions of gene expression profiles in different cell types. We develop a novel sparse least squares estimator, referred to as CSNet, that is efficient to implement and has good theoretical properties. Using CSNet, we analyzed the bulk gene expression data from a cohort study on Alzheimers disease and identified previously unknown cell-type-specific co-expressions among Alzheimers disease risk genes, suggesting cell-type-specific disease pathology for Alzheimers disease. The genome of the cereal pest Sitophilus oryzae: a transposable element haven 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. C4-dicarboxylate metabolons: Interaction of C4-dicarboxylate transporters of Escherichia coli with cytosolic enzymes Metabolons represent the structural organization of proteins for metabolic or regulatory pathways. Here the interaction of enzymes fumarase FumB and aspartase AspA with the C4-DC transporters DcuA and DcuB of Escherichia coli was tested by a bacterial two-hybrid (BACTH) assay in situ, or by co-chromatography (mSPINE). DcuB interacted strongly with FumB and AspA, and DcuA with AspA. The fumB-dcuB and the dcuA-aspA genes encoding the respective proteins are known for their colocalization on the genome and the production of co-transcripts. The data consistently suggest the formation of DcuB/FumB, DcuB/AspA and DcuA/AspA metabolons in fumarate respiration for the uptake of L-malate, or L-aspartate, conversion to fumarate and excretion of succinate after reduction. The DcuA/AspA metabolon catalyzes L-Asp uptake and fumarate excretion in concerted action also to provide ammonia for nitrogen assimilation. 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. Changes in Behavior and Salivary Serum Amyloid A in cows with Subclinical Mastitis The aim of this study was to identify detailed changes in behavior, and in salivary serum amyloid A (SAA), associated with subclinical mastitis. This included standard sickness behaviors (which form part of an adaptive response to conserve energy, minimize heat loss and facilitate recovery following infection and inflammation) and lesser-studied behaviors, that are biologically predicted to change at low-level infection, and therefore particularly relevant for detecting subclinical mastitis (here labelled  luxury behaviors). SAA is a physiological marker of systemic inflammation, with levels in milk and serum already known to increase during subclinical mastitis. We investigated whether the same was true for SAA in cow saliva. Data were collected for n = 34 commercial barn-housed dairy cows, tested in 17 matched-pairs. Each pair comprised a cow with subclinical mastitis (SCM) and a healthy control (CTRL), identified using somatic cell count (SCC) (SCM: SCC >200 x1000 cells/ml; CTRL: SCC <100 x1000 cells/ml). SCM cows were selected for study ad-hoc, at which point they were paired with a CTRL cow, based upon parity and calving date; consequently, the full data set was accrued over several months. Data were collected for each pair over 3 days: SCC (l4:00-l5:00h) Day 1; behavior (24h from 00:00h) Day 2; salivary serum amyloid-A (SAA) Day 3. We report, for the first time, that an increase in salivary SAA occurs during subclinical mastitis; SAA was higher in SCM cows and demonstrated a positive (weak) correlation with SCC. The behavioral comparisons revealed that SCM cows had reductions in activity (behavioral transitions and distance moved), social exploration, social reactivity (here: likelihood to move away/be displaced following receipt of agonism), performance of social grooming and head butts, and the receipt of agonistic non-contact challenges. In addition, SCM cows received more head swipes, and spent a greater proportion of time lying with their head on their flank than CTRL cows. SCM cows also displayed an altered feeding pattern; they spent a greater proportion of feeding time in direct contact with two conspecifics, and a lower proportion of feeding time at self-locking feed barriers, than CTRL cows. Behavioral measures were found to correlate, albeit loosely, with serum SAA in a direction consistent with predictions for sickness behavior. These included positive correlations with both lying duration and the receipt of all agonistic behavior, and negative correlations with feeding, drinking, the performance of all social and all agonistic behavior, and social reactivity. We conclude that changes in salivary SAA, social behavior, and activity offer potential in the detection of subclinical mastitis and recommend further investigation to substantiate and refine our findings. The Rise of West Nile Virus in Southern andSoutheastern Europe: a spatial-temporal analysisinvestigating the combined effects of climate, land use and economic changes West Nile Virus (WNV) has recently emerged as a major public health concern in Europe, its recent expansion also coincided with some remarkable socio-economic and environmental changes, including an economic crisis and some of the warmest temperatures on record. Here we empirically investigate the drivers of this phenomenon at a European wide scale by constructing and analyzing a unique spatial-temporal data-set, that includes data on climate, land-use, the economy, and government spending on environmental related sectors. Drivers and risk factors of WNV were identified by building a conceptual framework, and relationships were tested using a Generalized Additive Model (GAM), which could capture complex non-linear relationships and also account for the spatial and temporal auto-correlation. Some of the key risk factors identified in our conceptual framework, such as a higher percentage of wetlands and arable land, climate factors (higher summer rainfall and higher summer temperatures) were positive predictors of WNV infections. Interestingly, winter temperatures of between 2 {degrees}C and 6 {degrees}C were among some of the strongest predictors of annual WNV infections, one possible explanation for this result is that successful overwintering of infected adult mosquitoes (likely Culex pipiens) is key to the intensity of outbreaks for a given year. Furthermore, lower surface water extent over the summer is also associated with more more intense outbreaks, suggesting that drought, which is known to induce positive changes in WNV prevalence in mosquitoes is also contributing to the upward trend in WNV cases in affected regions. Our indicators representing the economic crisis were also strong predictors of WNV infections, suggesting there is an association between austerity and cuts to key sectors, which could have benefited vector species and the virus during this crucial period. These results, taken in the context of recent winter warming due to climate change, and more frequent droughts, may offer an explanation of why the virus has become so prevalent in Europe. Can single-component protein condensates form multiphase architectures? Phase-separated biomolecular condensates that contain multiple coexisting phases are widespread in vitro and in cells. Multiphase condensates emerge readily within multi-component mixtures of biomolecules (e.g. proteins and nucleic acids) when the different components present sufficient physicochemical diversity (e.g. in inter-molecular forces, structure, and chemical composition) to sustain separate coexisting phases. Because such diversity is highly coupled to the solution conditions (e.g. temperature, pH, salt, composition), it can manifest itself immediately from the nucleation and growth stages of condensate formation, develop spontaneously due to external stimuli, or progressively as the condensates age. Here, we investigate thermodynamic factors that can explain the intrinsic transformation of single-component condensates into multiphase architectures during the nonequilibrium process of ageing. We develop a multiscale model that integrates atomistic simulations of proteins, sequence-dependent coarse-grained simulations of condensates, and a minimal model of dynamically ageing condensates with non-conservative inter-molecular forces. Our nonequilibrium simulations of condensate ageing predict that single-component condensates that are initially homogeneous and liquid-like can transform into gel-core/liquid-shell or liquid-core/gel-shell multiphase condensates as they age, due to gradual and irreversible enhancement of inter-protein interactions. The type of multiphase architecture is determined by the ageing mechanism, the molecular organization of the gel and liquid phases, and the chemical make up of the protein. Notably, we predict that inter-protein disorder-to-order transitions within the prion-like domains of intracellular proteins could lead to the required non-conservative enhancement of inter-molecular interactions. Our study, therefore, predicts a potential mechanism

Significance StatementBiomolecular condensates are highly diverse systems spanning not only homogeneous liquid droplets, but also gels, glasses, and even multiphase architectures that contain various coexisting liquid-like and/or gel-like inner phases. Multiphase architectures form when the different biomolecular components in a multi-component condensate establish sufficiently imbalanced inter-molecular forces to sustain different coexisting phases. While such a requirement seems, at first glance, impossible to fulfil for a condensate formed exclusively of chemically-identical proteins (i.e., single-component), our simulations predict conditions under which this may be possible. During condensate ageing, a sufficiently large imbalance in inter-molecular interactions can emerge intrinsically from the accumulation of protein structural transitions--driving even single-component condensates into nonequilibrium liquid-core/gel-shell or gel-core/liquid-shell multiphase architectures. Trade-offs in biosensor optimization for dynamic pathway engineering Recent progress in synthetic biology allows the construction of dynamic control circuits for metabolic engineering. This technology promises to overcome many challenges encountered in traditional pathway engineering, thanks to their ability to self-regulate gene expression in response to bioreactor perturbations. The central components in these control circuits are metabolite biosensors that read out pathway signals and actuate enzyme expression. However, the construction of metabolite biosensors is a major bottleneck for strain design, and a key challenge is to understand the relation between biosensor dose-response curves and pathway performance. Here we employ multiobjective optimization to quantify performance trade-offs that arise in the design and calibration of metabolite biosensors. Our approach reveals strategies for tuning dose-response curves along an optimal trade-off between production flux and the cost of an increased expression burden on the host. We explore properties of control architectures built in the literature, and identify their advantages and caveats in terms of performance and robustness to growth conditions and leaky promoters. We demonstrate the optimality of a control circuit for glucaric acid production in Escherichia coli, which has been shown to increase titer by 2.5-fold as compared to static designs. Our results lay the groundwork for the automated design of control circuits for pathway engineering, with applications in the food, energy and pharmaceutical sectors. Embryonic exposure to cannabidiol disrupts active neural circuits, an effect increased by Δ-9-tetrahydrocannabinol and involving CB1R and CB2R in zebrafish In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. HOW CONCERNING IS A SARS-COV-2 VARIANT OF CONCERN? COMPUTATIONAL PREDICTIONS AND THE VARIANTS LABELING SYSTEM We herein report a study to evaluate the use of computational prediction of SARS-CoV-2 genetic variations in improving the current variants labeling system. First, we reviewed the basis of the system developed by the World Health Organization (WHO) for the labeling of SARS-CoV-2 genetic variants and the adaptations made to it by the United States Center of Diseases Control (CDC). We observed that the labeling system is based upon the virus major attributes. However, we found that the labeling criteria of the SARS-CoV-2 variants derived from these attributes are not accurately defined and are used differently by the two health management agencies. Consequently, discrepancies exist between the labels given by WHO and CDC to same variants. Our observations suggest that giving the VOC label to a new variant is premature and might not be appropriate. Therefore, we carried out a comparative computational study to predict the effects of the mutations on the virus structure and functions of five VOCs. By linking these data to the criteria used by WHO and the CDC for variant labeling, we ascertained that comparative computational predictions of the impact of genetic variations are a better ground for rapid and more accurate labelling of SARS-CoV-2 variants. We propose to label all emergent variants VUM or VBM and to carry out computational predictive studies and thorough variants comparison, upon which more appropriate and informative labels can be attributed. Furthermore, harmonization of the variants labeling system would be globally beneficial to communicate about and fight COVID-19 pandemic. 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. Quantitative analysis of gene expression on the inactive human X chromosome The human X chromosome exists in active (Xa) and inactive (Xi) states. To quantify Xi gene expression and explore its implications, we conducted RNA sequencing in two cell types from individuals with one Xa and zero to three Xis. Linear models fit the data well, indicating equal contributions from each Xi and quantifying significant Xi-driven expression changes for 124 of 406 expressed X-chromosome genes. As revealed by allele-specific expression analyses, these Xi-driven changes reflect both i) Xi transcription and ii) Xi modulation - both positive and negative - of Xa transcript levels. Incorporating population and evolutionary measures of constraint on expression levels, we identified 13 X-chromosome genes most likely to drive the phenotypic consequences of natural variation in Xi copy number. These analytics and the underlying datasets offer a foundation for understanding Turner, Klinefelter, and other sex chromosome aneuploidy syndromes, and sex-chromosome effects on health and disease in euploid males and females. Naturally occurring combinations of receptors from single cell transcriptomics in endothelial cells VEGF inhibitor drugs have been successful, especially in ophthalmology, but not all patients respond to them. Combinations of drugs are likely to be needed for more effective therapies of angiogenesis-related diseases. In this paper we describe naturally occurring combinations of receptors in endothelial cells that might help to understand how cells communicate and to identify targets for drug combinations. We also develop and share a new software tool called DECNEO to identify them.

Single-cell gene expression data are used to identify a set of co-expressed endothelial cell receptors, conserved among species (mice and humans) and enriched, within a network, of connections to up-regulated genes. This set includes several receptors previously shown to play a role in angiogenesis. Multiple statistical tests from large datasets, including an independent validation set, support the reproducibility, evolutionary conservation and role in angiogenesis of these naturally occurring combinations of receptors. We also show tissue-specific combinations and, in the case of choroid endothelial cells, consistency with both well-established and recent experimental findings, presented in a separate paper.

The results and methods presented here advance the understanding of signaling to endothelial cells. The methods are generally applicable to the decoding of intercellular combinations of signals. Extracellular Matrix Hydrogels Promote Expression of Muscle-Tendon Junction Proteins Muscle and tendon injuries are prevalent and range from minor sprains and strains to traumatic, debilitating injuries. However, the interactions between these tissues during injury and recovery remain unclear. Three-dimensional tissue models that incorporate both tissues and a physiologically relevant junction between muscle and tendon may help understand how the two tissues interact. Here, we use tissue specific extracellular matrix (ECM) derived from muscle and tendon to determine how cells of each tissue interact with the microenvironment of the opposite tissue resulting in junction specific features. ECM materials were derived from the Achilles tendon and gastrocnemius muscle, decellularized, and processed to form tissue specific pre-hydrogel digests. ECM materials were unique in respect to protein composition and included many types of ECM proteins, not just collagens. After digestion and gelation, ECM hydrogels had similar complex viscosities which were less than type I collagen hydrogels at the same concentration. C2C12 myoblasts and tendon fibroblasts were cultured in tissuespecific ECM conditioned media or encapsulated in tissue-specific ECM hydrogels to determine cell-matrix interactions and the effects on a muscle-tendon junction marker, paxillin. ECM conditioned media had only a minor effect on upregulation of paxillin in cells cultured in monolayer. However, cells cultured within ECM hydrogels had 50-70% higher paxillin expression than cells cultured in type I collagen hydrogels. Contraction of the ECM hydrogels varied by the type of ECM used. Subsequent experiments with varying density of type I collagen (and thus contraction) showed no correlation between paxillin expression and the amount of gel contraction, suggesting that a constituent of the ECM was the driver of paxillin expression in the ECM hydrogels. In addition, the extracellular matrix protein type XXII collagen had similar expression patterns as paxillin, with smaller effect sizes. Using tissue specific ECM allowed for the de-construction of the cell-matrix interactions similar to muscletendon junctions to study the expression of MTJ specific proteins.

Impact StatementThe muscle-tendon junction is an important feature of muscle-tendon units; however, despite crosstalk between the two tissue types, it is overlooked in current research. Deconstructing the cell-matrix interactions will provide the opportunity to study significant junction specific features and markers that should be included in tissue models of the muscletendon unit, while gaining a deeper understanding of the natural junction. This research aims to inform future methods to engineer a more relevant multi-tissue platform to study the muscletendon unit. NUDT18 catalyzes the hydrolysis of active metabolites of the antivirals Remdesivir, Ribavirin and Molnupiravir Remdesivir and Molnupiravir have gained considerable interest due to their activity against SARS-CoV-2. Cellular hydrolysis of their active triphosphate forms, Remdesivir-TP and Molnupiravir-TP, would decrease drug efficiency. We therefore tested Remdesivir-TP as a substrate against a panel of human hydrolases and found that NUDT18 catalyzes the hydrolysis of Remdesivir-TP. The kcat value of NUDT18 for Remdesivir-TP was determined to 2.6 s-1 and the Km value was 156 M, suggesting that NUDT18 catalyzed hydrolysis of Remdesivir-TP occurs in cells. We demonstrate that the triphosphates of the antivirals Ribavirin and Molnupiravir are hydrolyzed by NUDT18, albeit with a lower efficiency compared to Remdesivir-TP. NUDT18 also hydrolyses the triphosphates of Sofosbuvir and Aciclovir although with significantly lower activity. These results suggest that NUDT18 can act as a cellular sanitizer of modified nucleotides and may influence the antiviral efficacy of Remdesivir, Molnupiravir and Ribavirin. NUDT18 is expressed in respiratory epithelial cells and may limit the antiviral efficacy of Remdesivir and Molnupiravir against SARS-CoV2 replication by decreasing the intracellular concentration of their active metabolites at their intended site of action. Bovine neutrophils release extracellular traps and cooperate with macrophages in Mycobacterium avium subsp. paratuberculosis clearance in vitro Mycobacterium avium subsp. paratuberculosis (Map) is the underlying pathogen causing bovine paratuberculosis (PTB), an enteric granulomatous disease that mainly affects ruminants and for which an effective treatment is needed. Macrophages are the primary target cells for Map, which survives and replicates intracellularly by inhibiting phagosome maturation. Neutrophils are present at disease sites during the early stages of the infection, but seem to be absent in the late stage, in contrast to healthy tissue. Although neutrophil activity has been reported to be impaired following Map infection, their role in PTB pathogenesis has not been fully defined. Neutrophils are capable of releasing extracellular traps consisting of extruded DNA and proteins that immobilize and kill microorganisms, but this mechanism has not been evaluated against Map. Our main objective was to study the interaction of neutrophils with macrophages during an in vitro mycobacterial infection. For this purpose, neutrophils and macrophages from the same animal were cultured alone or together in the presence of Map or Mycobacterium bovis Bacillus-Calmette-Guerin (BCG). Extracellular trap release, mycobacteria killing as well as IL-1{beta} and IL-8 release were assessed. Neutrophils released extracellular traps against mycobacteria when cultured alone and in the presence of macrophages without direct cell contact, but resulted inhibited in direct contact. Macrophages were extremely efficient at killing BCG, but ineffective at killing Map. In contrast, neutrophils showed similar killing rates for both mycobacteria. Co-cultures infected with Map showed the expected killing effect of combining both cell types, whereas co-cultures infected with BCG showed a potentiated killing effect beyond the expected one, indicating a potential synergistic cooperation. In both cases, IL-1{beta} and IL-8 levels were lower in co-cultures, suggestive of a reduced inflammatory reaction. These data indicate that cooperation of both cell types can be beneficial in terms of decreasing the inflammatory reaction while the effective elimination of Map can be compromised. These results suggest that neutrophils are effective at Map killing and can exert protective mechanisms against Map that seem to fail during PTB disease after the arrival of macrophages at the infection site. Transcripts switched off at the stop of phloem unloading highlight the energy efficiency of sugar import in the ripening V. vinifera fruit. Transcriptomic changes at the cessation of sugar accumulation in the pericarp of Vitis vinifera were addressed on single berries re-synchronized according to their individual growth patterns. The net rates of water, sugars and K+ accumulation inferred from individual growth and solute concentration confirmed that these inflows stopped simultaneously in the ripe berry, while the small amount of malic acid remaining at this stage was still being oxidized at low rate. Re-synchronized individual berries displayed negligible variations in gene expression among triplicates. RNA-Seq studies revealed sharp reprogramming of cell wall enzymes and structural proteins at the stop of phloem unloading, associated with an 80% repression of multiple sugar transporters and aquaporins on the plasma or tonoplast membranes, with the noticeable exception of H+/sugar symporters, that were rather weakly and constitutively expressed. This was verified in three genotypes placed in contrasted thermo-hydric conditions. The prevalence of SWEET suggests that electrogenic transporters would play a minor role on the plasma membranes of SE/CC complex and the one of the flesh, while sucrose/H+ exchangers dominate on its tonoplast. Cis-regulatory elements present in their promoters allowed to sort these transporters in different groups, also including specific TIPs and PIPs paralogs, and cohorts of cell wall related genes. Together with simple thermodynamic considerations, these results lead to propose that H+/sugar exchangers at the tonoplast, associated with a considerably acidic vacuolar pH, may exhaust cytosolic sugars in the flesh and alleviate the need for supplementary energization of sugar transport at the plasma membrane. The initial response of females towards congeneric males matches the propensity to hybridise in Ophthalmotilapia. AO_SCPLOWBSTRACTC_SCPLOWCichlid radiations often harbour closely related species with overlapping niches and distribution ranges. Such species sometimes hybridise in nature, which raises the question how can they coexist. This also holds for the Tanganyika mouthbrooders Ophthalmotilapia ventralis and O. nasuta. Earlier studies found indications of asymmetrical hybridisation with females of O. ventralis accepting males of O. nasuta, but not the other way around. We hypothesised that this was due to differences in the capacity for species recognition. Given the higher propensity of O. ventralis females towards hybridisation, we expect a reduced ability for species recognition in O. ventralis females, compared to O. nasuta females. We staged two experiments, one focusing on 22 female O. nasuta and one on 21 female O. ventralis. These fish were placed in one half of a tank and briefly exposed to a conspecific or a heterospecific male, a conspecific female, or nothing (control). Female response was evaluated by scoring six tracking parameters and by noting the occurrence of ten discrete behaviours before and during the encounter. Females always responded to the presence of another fish by approaching it. Remarkably, for both O. nasuta and O. ventralis, we did not find a different response between encounters with conspecific males and females. However, in agreement with our hypothesis, O. nasuta females behaved differently towards conspecific or heterospecific males, whereas O. ventralis females did not. When presented with a heterospecific male, O. nasuta females performed a lower number of  ram behaviours. Additionally, they never displayed the  flee behaviour, a component of the species mating repertoire that was seen in all but one of the presentations with a conspecific male. Our findings show that differences in species recognition at first encounter predict to a large degree the outcome of the mating process, even in the absence of mating behaviour. G protein-coupled receptor signaling regulates ER-mitochondria contacts Interactions between the endoplasmic reticulum (ER) and mitochondria (Mito) are crucial for many cellular functions, and their interaction levels change dynamically depending on the cellular environment. Little is known about how the interactions between these organelles are regulated within the cell. Here we screened a compound library to identify chemical modulators for ER-Mito contacts. Multiple agonists of G-protein coupled receptors (GPCRs), beta-adrenergic receptors ({beta}-ARs) in particular, scored in this screen. Analyses in multiple orthogonal assays validated that GPCR activation promotes physical and functional interactions between the two organelles. Furthermore, we have elucidated potential downstream effectors mediating GPCR-induced ER-Mito contacts. Given the role of GPCRs in sensing various extracellular signals and subsequently eliciting appropriate cellular responses, our data provide significant insights into the roles of ER-Mito contacts in maintaining cellular homeostasis and in responding to physiological demands or stresses. 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.