Phylogenetic and population structure analyses uncover pervasive misclassification and help assessing the biosafety of Pseudomonas alloputida for biotechnological applications The Pseudomonas putida group comprises strains with biotechnological and clinical relevance. P. alloputida was proposed as a new species and highlighted the misclassification of P. putida. Nevertheless, the population structure of P. alloputida remained unexplored. We retrieved 11,025 Pseudomonas genomes and used P. alloputida Kh7T to delineate the species. The P. alloputida population structure comprises at least 7 clonal complexes (CCs). Clinical isolates are mainly found in CC4 and acquired resistance genes are present at low frequency in plasmids. Virulence profiles support the potential of CC7 members to outcompete other plant or human pathogens through a type VI secretion system. Finally, we found that horizontal gene transfer had an important role in shaping the ability of P. alloputida to bioremediate aromatic compounds such as toluene. Our results provide the grounds to understand P. alloputida genetic diversity and safety for environmental applications.
Phylogenomic evidence for the Origin of Obligately Anaerobic Anammox Bacteria around the Great Oxidation Event The anaerobic ammonium oxidation (anammox) bacteria could transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for up to half of the global nitrogen loss in surface environments. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remains enigmatic. Here, we performed comprehensive phylogenomic analysis and showed a single origin of anammox bacteria within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates, which includes implementing both a traditional cyanobacteria-based and a recently developed mitochondria-based approach, we estimated that anammox bacteria originated at early Proterozoic and most likely around the so-called Great Oxidation Event (GOE; 2.32 to 2.5 billion years ago [Ga]) which fundamentally changed global biogeochemical cycles. We further showed that during the origin of anammox bacteria, genes involved in oxidative stress, bioenergetics and anammox granules formation were recruited, which might have contributed to their survival on an increasingly oxic Earth. Our findings suggest the rising levels of atmospheric oxygen, which made nitrite increasingly available, was a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link the GOE to the evolution of obligate anaerobic bacteria.
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. Innovative GenExpA software for selecting suitable reference genes for reliable normalization of gene expression in melanoma The algorithms commonly used to select the best stable reference gene in RT-qPCR data analysis have their limitations. We showed that simple selection of the reference gene or pair of genes with the lowest stability value from the pool of potential reference genes - a commonly used approach - is not sufficient to accurately and reliably normalize the target gene transcript and can lead to biologically incorrect conclusions. For reliable assessment of changes in a target gene expression level, we propose our innovative GenExpA software, which works in a manner independent of the experimental model and the normalizer used. GenExpA software selects the best reference by combining the NormFinder algorithm with progressive removal of the least stable gene from the candidate genes in a given experimental model and in the set of daughter models assigned to it. The reliability of references is validated based on the consistency of the statistical analyses of normalized target gene expression levels through all models, described by the coherence score (CS). The use of the CS value imparts a new quality to qPCR analysis because it clarifies how low the stability value of reference must be in order for biologically correct conclusions to be drawn. We tested our method on qPCR data for the B4GALT genes family in melanoma, which is characterized by a high mutation rate, and in melanocytes. GenExpA is available at https://github.com/DorotaHojaLukowicz/GenExpA or https://www.sciencemarket.pl/baza-programow-open-source#oferty.

HighlightsGenExpA - next-generation software for normalizer selection and validation

The GenExpA tool defines how low the stability value of the reference should be

The GenExpA tool determines the level of gene expression analysis robustness

The GenExpA tool increases the speed of analysis and reduces cost Osr40g3 imparts salinity tolerance by regulating GF14e-mediated modulation of gibberellin metabolism to activate expansin protein, EG45 in rice While the role of R40 family of lectin proteins in osmotic stress was already reported, their mechanism of stress tolerance still remained elusive. Among them, the expression of Osr40g3 was most strongly induced in response to salinity. Here, we explored the mechanism of Osr40g3- mediated salinity tolerance in rice. The overexpression of Osr40g3 under control of a constitutive promoter significantly improved salinity tolerance but displayed pollen sterility and impaired seed development. The gene further displayed precise tissue-specific expression which was essential for proper seed development in rice. In salinity-stressed leaves, Osr40g3 interacted with OsEG45 leading to salinity tolerance. Overexpression of OsEG45 significantly improved salinity tolerance while the OsEG45 silencing lines displayed severe sensitivity thus confirming its role in imparting salinity tolerance. On the other hand, the constitutive overexpression of Osr40g3 led to its interaction with OsGF14e in stamens reducing its protein stability. This interaction further interfered with gibberellic acid (GA) metabolism thus leading to pollen sterility. The OsGF14e silencing lines also failed to set seeds and displayed diminished GA content thus supporting its involvement in regulating seed development. Surprisingly, these OsGF14e silencing lines also displayed improved salinity tolerance with higher expression of OsEG45 while the OsGF14e overexpression lines displayed opposite phenotype. This altered expression of OsEG45 was further mediated by GA levels in these transgenic lines. Together, we demonstrate that salt stress induces Osr40g3 which then positively regulates OsEG45 and negatively regulates OsGF14e. Inhibition of OsGF14e reduces GA level which further induces OsEG45 expression thus leading to salinity tolerance in rice. Computational and experimental studies of the breathing motion of a protein loop: implications in PfAMA1-PfRON2 late-stage binding event Protein-protein interactions are important targets for various drug discovery campaigns. One such promising and therapeutically pertinent protein-protein complex is PfAMA1-PfRON2 involved in malarial parasite invasion into human red blood cells. A thorough understanding of the interactions between these macromolecular binding partners is crucial for designing better therapeutics against this age-old disease. Although crystal structures of several PfAMA1-PfRON2 complexes are available, the mechanism of how domain II loop associates with PfRON2 is not clear. The current work investigates how the domain II loop of PfAMA1 exerts its effect on the alpha helix of the PfRON2, thus influencing the overall kinetics of this intricate recognition phenomenon. To this end, we have computationally simulated the dynamics and free energetics of domain II loop closing processes and identified a set of key amino acid residues of PfRON2 helix which are essential for binding. The subsequent evaluation of the binding affinity of Ala-substituted PfRON2 peptide ligands by surface plasmon resonance (SPR) and isothermal titration calorimetry (ITC) validates the relative importance of the residues in context. Together, the combination of computational and experimental investigation reveals that the domain II loop of PfAMA1 is in fact responsible for arresting the PfRON2 molecule from egress, K2027 and D2028 of PfRON2 being the determinant residues for the capturing event. Our study provides a comprehensive understanding of the molecular recognition event between PfAMA1 and PfRON2, specifically in the post binding stage, which could potentially open up new avenues to drug discovery against malaria.

TOC GRAPHIC

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=95 SRC="FIGDIR/small/463826v2_ufig1.gif" ALT="Figure 1">
View larger version (35K):
org.highwire.dtl.DTLVardef@1c26022org.highwire.dtl.DTLVardef@356357org.highwire.dtl.DTLVardef@10ceb04org.highwire.dtl.DTLVardef@1195633_HPS_FORMAT_FIGEXP  M_FIG C_FIG wnt16 exerts pleiotropic effects on bone and lean mass in zebrafish Osteoporosis, a disease of bone fragility, and sarcopenia, a condition of reduced muscle mass and strength, commonly occur in the same individual--a condition termed osteosarcopenia. Genetic variants at the CPED1-WNT16 locus have pleiotropic effects on human bone mineral density (BMD) and lean tissue mass. However, our current understanding of Wnt16 and its influence on bone fails to explain effects on lean mass at this locus. In this study, we sought to better understand the biological mechanism underlying pleiotropy at the CPED1-WNT16 locus on BMD and lean mass. Using single cell analysis, microCT imaging, and genetic approaches in zebrafish, we show that wnt16 exerts dual influence on muscle and bone, and provide evidence that this occurs through a function in somite development, a period in which precursors for both of these tissues are specified. We show that mutations in wnt16 alter vertebral size and lean mass during axial skeletogenesis. We show that wnt16+ cells transcriptionally resemble cells within the amniote dermomyotome, and provide evidence they function as myogenic precursors during embryonic development. We also provide evidence that WNT16 is a gene of major effect on lean mass at the CPED-WNT16 locus. Our studies support the potential for wnt16 to exert pleiotropic effects on bone and lean mass through a function in myogenic precursors. The voltage and spiking responses of subthreshold resonant neurons to structured and fluctuating inputs: persistence and loss of resonance and variability We systematically investigate the response of neurons to oscillatory currents and synaptic-like inputs and we extend our investigation to non-structured synaptic-like spiking inputs with more realistic distributions of presynaptic spike times. We use two types of chirp-like inputs consisting of (i) a sequence of cycles with discretely increasing frequencies over time, and (ii) a sequence having the same cycles arranged in an arbitrary order. We develop and use a number of frequency-dependent voltage response metrics to capture the different aspects of the voltage response, including the standard impedance (Z) and the peak-to-trough amplitude envelope (VENV) profiles. We show that Z-resonant cells (cells that exhibit subthreshold resonance in response to sinusoidal inputs) also show VENV-resonance in response to sinusoidal inputs, but generally do not (or do it very mildly) in response to square-wave and synaptic-like inputs. In the latter cases the resonant response using Z is not predictive of the preferred frequencies at which the neurons spike when the input amplitude is increased above subthreshold levels. We also show that responses to conductance-based synaptic-like inputs are attenuated as compared to the response to current-based synaptic-like inputs, thus providing an explanation to previous experimental results. These response patterns were strongly dependent on the intrinsic properties of the participating neurons, in particular whether the unperturbed Z-resonant cells had a stable node or a focus. In addition, we show that variability emerges in response to chirp-like inputs with arbitrarily ordered patterns where all signals (trials) in a given protocol have the same frequency content and the only source of uncertainty is the subset of all possible permutations of cycles chosen for a given protocol. This variability is the result of the multiple different ways in which the autonomous transient dynamics is activated across cycles in each signal (different cycle orderings) and across trials. We extend our results to include high-rate Poisson distributed current- and conductance-based synaptic inputs and compare them with similar results using additive Gaussian white noise. We show that the responses to both Poisson-distributed synaptic inputs are attenuated with respect to the responses to Gaussian white noise. For cells that exhibit oscillatory responses to Gaussian white noise (band-pass filters), the response to conductance-based synaptic inputs are low-pass filters, while the response to current-based synaptic inputs may remain band-pass filters, consistent with experimental findings. Our results shed light on the mechanisms of communication of oscillatory activity among neurons in a network via subthreshold oscillations and resonance and the generation of network resonance. Discovery of highly potent pancoronavirus fusion inhibitors that also effectively inhibit COVID-19 variants from the UK (Alpha), South Africa (Beta), and India (Delta) We report here the discovery of several highly potent small molecules that showed low nM potency against SARS-CoV (IC50: as low as 13 nM), SARS-CoV-2 (IC50: as low as 23 nM), and MERS-CoV (IC50: as low as 76 nM) in pseudovirus based assays with excellent selectivity indices (SI: as high as > 5000) demonstrating their pancoronavirus inhibition. Some compounds also show 100% inhibition of CPE (IC100) at 1.25 {micro}M against an authentic SARS-CoV-2 (US_WA-1/2020). Furthermore, the most active inhibitors also potently inhibited variants of concerns (VOCs), such as the UK (B.1.1.7), South Africa (B.1.351), and Delta variant (B.1.617.2), originated in India. We confirmed that one of the potent inhibitors binds to the prefusion spike protein trimer of SARS-CoV-2 by SPR. Besides, we showed that they inhibit virus-mediated cell-cell fusion. The ADME data of one of the most active inhibitors, NBCoV1, show drug-like properties. In vivo PK of NBCoV1 in rats demonstrated excellent half-life (t1/2) of 11.3 h, mean resident time (MRT) of 14.2 h, and oral bioavailability. We expect the lead inhibitors to pave the way for further development to preclinical and clinical candidates. First evidence of in vitro cytotoxic effects of marine microlitter on Merluccius merluccius and Mullus barbatus, two Mediterranean commercial fish species Marine litter is composed mainly of plastics and is recognized as a serious threats to marine ecosystems. Ecotoxicological approaches have started elucidating the potential severity of microplastics (MPs) in controlled laboratory studies with pristine materials but no information exist on marine environmental microlitter as a whole. Here, we characterized the litter in the coastal Northern Tyrrhenian sea and in the stomach of two fish species of socio-economic importance, and exposed primary cell cultures of mucosal and lymphoid organs to marine microlitter for evaluating possible cytotoxic effects. An average of 0.30 {+/-} 0.02 microlitter items m-3 was found in water samples. FT-IR analysis revealed that plastic particles, namely HDPE, polyamide and polypropylene were present in 100% and 83.3% of Merluccius merluccius and Mullus barbatus analyzed, which overall ingested 14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual, respectively. Moreover, microlitter was confirmed as a vector of microorganisms. Lastly, the apical end-point of viability was found to be significantly reduced in splenic cells exposed in vitro to two microlitter conditions. Considering the role of the spleen in the mounting of adaptive immune responses, our results warrant more in-depth investigations for clarifying the actual susceptibility of these two species to anthropogenic microlitter.

HighlightsO_LI0.30 {+/-} 0.02 microlitter items m-3 were found at the surface of coastal Northern Tyrrhenian sea
C_LIO_LI14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual were retrieved from the stomach of hakes and mullets
C_LIO_LIThe ingested microlitter contained plastic items
C_LIO_LIMicrolitter was validated as a carrier of bacteria, fungi and flagellates
C_LIO_LISplenic cells exposed to two microlitter conditions for 72 hours suffered cytotoxicity
C_LI Sensing their plasma membrane curvature allows migrating cells to circumvent obstacles Cell migration is a hallmark out-of-equilibrium process in biology. In addition to persistent self-propelled motion, many cells display remarkable adaptive behaviors when they navigate complex environments within the body. Combining theory and experiments, we identify a curvature-sensing mechanism underlying obstacle avoidance in immune-like cells. The genetic perturbation of this machinery leads to a reduced capacity to evade obstructions combined with faster and more persistent cell migration in obstacle-free environments. We propose that the active polymerization of the actin cytoskeleton at the advancing edge of migrating cells is locally inhibited by the curvature-sensitive BAR protein Snx33 in regions with inward plasma membrane curvature. This coupling between actin and membrane dynamics leads to a mechanochemical instability that generates complex protrusive patterns at the cellular front. Adaptive motility thus arises from two simultaneous curvature-dependent effects, i) the specific reduction of propulsion in regions where external objects deform the plasma membrane and ii) the intrinsic patterning capacity due to the membrane-actin coupling that promotes spontaneous changes in the cells protrusions. Our results show how cells utilize actin- and plasma membrane biophysics to sense their environment, allowing them to adaptively decide if they should move ahead or turn away. On the basis of our findings, we propose that the natural diversity of BAR proteins may allow cells to tune their curvature sensing machinery to match the shape characteristics in their environment. Nuclear alpha-synuclein is present in the human brain and is modified in dementia with Lewy bodies Dementia with Lewy bodies is pathologically defined by the cytoplasmic accumulation of alpha-synuclein within neuronal cells in the brain. Alpha-synuclein is predominately pre-synaptic, but has been reported present in various subcellular compartments in cell and animal models. In particular, nuclear alpha-synuclein is evident in-vitro and in disease models and has been associated with altered DNA integrity, gene transcription, nuclear homeostasis. However, owing to various factors, the presence of alpha-synuclein in the nuclei of human brain cells remains controversial, as does its role in synucleinopathies.

Here, we close this gap and provide a unique demonstration confirming the presence of nuclear alpha-synuclein in post-mortem brain tissue obtained from cases of dementia with Lewy bodies as well as from controls via immunohistochemistry, immunoblot, and label-free mass-spectrometry.

Discrete intra-nuclear alpha-synuclein puncta reactive against phosphorylated serine 129-alpha-synuclein and pan-alpha-synuclein antibodies were observed in cortical neurons and non-neuronal cells in fixed brain sections and in isolated nuclear preparations from Dementia with Lewy bodies cases and matched controls. Subsequent biochemical analysis of subcellular fractionated tissue confirmed alpha-synuclein as present in a nuclear fraction at levels ~ 10-fold lower than in the cytoplasm. Critically, however, an increase in monomeric nuclear alpha-synuclein phosphorylated as serine 129 was observed in cases of dementia with Lewy bodies alongside higher molecular weight pan- and phosphorylation reactive alpha-synuclein species, consistent with the formation of intranuclear phosphorylated alpha-synuclein oligomers. Furthermore, the presence of nuclear alpha-synuclein was confirmed via label free mass spectrometry, as 6 unique alpha-synuclein derived peptide sequences were identified in nuclear fractions (71.4% sequence coverage).

Collectively, our data confirm the presence of nuclear alpha-synuclein in human brain tissue and describe nuclear pathology associated with dementia with Lewy bodies. These findings address a major controversy in the synucleinopathy field by confirming the presence of nuclear alpha-synuclein in autoptic human brain tissue and, for the first time, identify that alpha-synuclein is aggregated into novel and potentially pathological assemblies in the nucleus as part of the disease process associated with dementia with Lewy bodies and thus may contribute to the disease phenotype. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Computation of the mitochondrial age distribution along the axon length We describe a compartmental model of mitochondrial transport in axons, which we apply to compute mitochondrial age at different distances from the soma. The model predicts that the average mitochondrial age at the tip of an axon with a length of 1 cm is approximately 22 hours. The mitochondria are youngest closest to the soma and their age scales approximately linearly with distance from the soma. To the best of the authors knowledge, this is the first attempt to predict the spatial distribution of mitochondrial age within an axon. A sensitivity study of the mean age of mitochondria to various model parameters is also presented. Towards a unified model of naive T cell dynamics across the lifespan Naive CD4 and CD8 T cells are cornerstones of adaptive immunity, but the dynamics of their establishment early in life and how their kinetics change as they mature following release from the thymus are poorly understood. Further, due to the diverse signals implicated in naive T cell survival, it has been a long-held and conceptually attractive view that they are sustained by active homeostatic control as thymic activity wanes. Here we employ multiple experimental systems to identify a unified model of naive CD4 and CD8 T cell population dynamics across the mouse lifespan. We infer that both subsets divide rarely and progressively increase their survival capacity with cell age. Strikingly, this simple model captures naive CD4 T cell dynamics throughout life. In contrast, we find that newly generated naive CD8 T cells are lost more rapidly during the first 3-4 weeks of life, likely due to increased recruitment into memory. We find no evidence for elevated division rates in neonates, or for feedback regulation of naive T cell numbers at any age. We show how confronting mathematical models with diverse datasets can reveal a quantitative and remarkably simple picture of naive T cell dynamics from birth into old age. The structurome of a Clostridium difficile phage and the remarkable accurate prediction of its novel phage receptor-binding protein As natural bacterial predators, bacteriophages have the potential to be developed to tackle antimicrobial resistance, but our exploitation of them is limited by understanding their vast uncharacterised genetic diversity1,2. Fascinatingly, this genetic diversity reflects many ways that phages can make proteins, performing similar functions that together form the familiar phage particle. Critical to infection are phage receptor-binding proteins (RBPs) that bind bacterial  receptors and initiate bacterial entry3. Here we identified and characterised Gp22, a novel RBP for phage CDHS-1 that infects pathogenic C. difficile, but that had no recognisable RBPs. We showed that Gp22 antibodies neutralised CDHS-1 infection and used immunogold-labelling and transmission electron microscopy to identify their location on the capsid. The Gp22 three-dimensional structure was resolved by X-ray crystallography revealing a new RBP class with an N-terminal L-shaped -helical superhelix domain and a C-terminal Mg2+-binding domain. The findings provide novel insights into C. difficile phage biology and phage-host interactions. This will facilitate optimal phage development and future engineering strategies4,5. Furthermore, the AlphaFold2-predicted Gp22 structure, which was strikingly accurate, paves the way for a structurome based transformation and guidance of future phage studies where many proteins lack sequence homology but have recognisable protein structures. Characterization of methylation profiles in spontaneous preterm birth placental villous tissue Preterm birth is a global public health crisis which results in significant neonatal and maternal mortality. Yet little is known regarding the molecular mechanisms of idiopathic spontaneous preterm birth, and we have few diagnostic markers for adequate assessment of placental development and function. Previous studies of placental pathology and our transcriptomics studies suggest a role for placental maturity in idiopathic spontaneous preterm birth. It is known that placental DNA methylation changes over gestation. We hypothesized that if placental hypermaturity is present in our samples, we would observe a unique idiopathic spontaneous preterm birth DNA methylation profile potentially driving the gene expression differences we previously identified in our placental samples. Our results indicate the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern suggesting hypermaturity. Only seven significant differentially methylated regions fitting the idiopathic spontaneous preterm birth specific (relative to the controls) profile were identified, indicating unusually high similarity in DNA methylation between idiopathic spontaneous preterm birth and term birth samples. We identified an additional 1,718 significantly methylated regions in our gestational age matched controls were the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern, again indicating a striking level of similarity between the idiopathic spontaneous preterm birth and term birth samples. Pathway analysis of these regions revealed differences in genes within the WNT and Cadherin signaling pathways, both of which are essential in placental development and maturation. Taken together, these data demonstrate that the idiopathic spontaneous preterm birth samples are molecularly more mature than expected given their respective gestational age which likely impacts birth timing. Bone marrow haematopoietic stem cells influence liver homeostatic networks and cancer development In fatty liver disease, systemic homeostasis is perturbed. While pre-clinical models are used to understand its pathogenesis, translating this knowledge to patients is difficult. However, by focusing on the most preserved homeostasis systems between species and models, novel disease dimensions can be unearthed. We interrogated core liver gene co-expression networks in a mouse model of liver cancer following dietary challenge. Immunometabolic modules showed temporal changes under the influence of diet duration and aging. The behaviour of immune network in tumours mirrored their counterparts in non-tumour liver. A high immune response network was associated with a lower tumour burden in mice and humans. This module in mice was enriched for genes related to haematopoietic cell differentiation. Consistently, the bone marrow haematopoietic stem and progenitor cells response was reflective of the liver immune response. Linking haematopoiesis to hepatic homeostasis uncovers a hitherto unexplored dimension of tissue crosstalk that can inform pathogenesis. Evaluation of methods for estimating coalescence times using ancestral recombination graphs The ancestral recombination graph (ARG) is a structure that describes the joint genealogies of sampled DNA sequences along the genome. Recent computational methods have made impressive progress towards scalably estimating whole-genome genealogies. In addition to inferring the ARG, some of these methods can also provide ARGs sampled from a defined posterior distribution. Obtaining good samples of ARGs is crucial for quantifying statistical uncertainty and for estimating population genetic parameters such as effective population size, mutation rate, and allele age. Here, we use standard neutral coalescent simulations to benchmark the estimates of pairwise coalescence times from three popular ARG inference programs: ARGweaver, Relate, and tsinfer+tsdate. We compare 1) the true coalescence times to the inferred times at each locus; 2) the distribution of coalescence times across all loci to the expected exponential distribution; 3) whether the sampled coalescence times have the properties expected of a valid posterior distribution. We find that inferred coalescence times at each locus are most accurate in ARGweaver, and often more accurate in Relate than in tsinfer+tsdate. However, all three methods tend to overestimate small coalescence times and underestimate large ones. Lastly, the posterior distribution of ARGweaver is closer to the expected posterior distribution than Relates, but this higher accuracy comes at a substantial trade-off in scalability. The best choice of method will depend on the number and length of input sequences and on the goal of downstream analyses, and we provide guidelines for the best practices. Failure to account for behavioral variability significantly compromises accuracy in indirect population monitoring Indirect wildlife population surveying largely depends upon counts of artefacts of behavior (e.g., nests or dung). Likelihood to encounter these artefacts is derived from both artefact production and decay, and variability in artefact production behavior is considered to contribute minimally to inaccuracy in wildlife estimation. Here we demonstrate how ignoring behavioral variability contributes to significant population misestimation, using an example of an endangered ape, the bonobo (Pan paniscus). Until now, a single estimate of (i.e., nest construction) has been used to extrapolate bonobo densities, assumed to be broadly representative of bonobo sign production behavior. We estimated nest construction rate across seasons and social groups at the Kokolopori Bonobo Reserve, DRC and find that nest construction rates in bonobos to be highly variable across populations as well as seasonal. Failure to account for this variability led to potentially severe degradation in accuracy of bonobo population estimates of abundance, accounting for a likely overestimation of bonobo numbers by 34%, and at worst as high as 80%. With this example, we demonstrate that failure to account for inter- and intra-population behavioral variation compromises the ability to estimate both relative and absolute wildlife abundances. We argue that variation in sign production is but one of several potential ways that behavioral variability can affect conservation monitoring, should be measured across contexts whenever possible, and must be considered in population estimation confidence intervals. With increasing attention to behavioral variability as a potential tool for conservation, conservationists must also account for the impact that behavioral variability can play upon wildlife population estimation. Our results underline the importance of observational research to wildlife monitoring schemes as a critical but overlooked component of conservation management. We extensively discuss the various avenues through which behavioral variability is likely to impact wildlife monitoring accuracy and precision and propose potential avenues for accounting for behavioral variability in wildlife monitoring. A network approach to genetic circuit designs As genetic circuits become more sophisticated, the size and complexity of data about their designs increases. This data captured goes beyond monolithic genetic sequences and towards circuit modularity and functional details, which are beneficial for analyzing circuit performance and establishing design automation techniques. However, the accessibility, visualisation and usability of design data (and metadata) have received relatively little attention to date. Here, we present a method to turn circuit designs into networks and showcase its potential to enhance the utility of design data. Since networks are dynamic structures, initial graphs can be interactively shaped into sub-networks of relevant information based on requirements such as abstraction, hierarchy and protein interactions. Additionally, several visual changes can be applied, such as colouring or clustering nodes based on types (e.g., genes or promoters), resulting in easier comprehension from a user perspective. This approach allows circuit designs to be coupled to other networks, such as metabolic pathways or implementation protocols captured in graph-like formats. Therefore, we advocate using networks to structure, access and improve synthetic biology information. The role of Kabuki Syndrome genes KMT2D and KDM6A in development: Analysis in Human sequencing data and compared to mice and zebrafish KMT2D and KDM6A are epigenetic regulators that have been implicated in Kabuki Syndrome, a rare congenital birth defect with multiple tissue and organ abnormalities, including craniofacial and heart defects. Our previous study identified human families with mutations in the epigenetic modifiers KMT2D and KDM6A, which is implicated in 32% and 10% of Kabuki Syndrome patients respectively. To understand the connection to Kabuki syndrome patients, and the transcriptional targets of KMT2D and KDM6A in humans, we performed RNA sequencing (seq) of lymphoblastoid cells from Kabuki Syndrome patients carrying mutations in KMT2D and KDM6A. We identified 1995 significant changes in transcriptional targets for KMT2D and 1917 for KDM6A, as compared to control. When compared with RNA-seq datasets obtained from other mouse and zebrafish studies, our analysis revealed KMT2D mutations affect the expression of 76 orthologous genes across all three datasets. Similarliy, KDM6A afftects the expprssion of 7 orthologous genes across three datasets. Despite the differences in cell types, stages, and species in the comparison between the transcriptomic datasets, there are common gene expression changes associated with KMT2D and KDM6A mutations. qPCR on novel zebrafish mutants confirmed the differentially expression in KMT2D or KDM6A mutant backgrounds. Taken together, our results show that KMT2D and KDM6A regulate common and unique genes across humans, mice, and zebrafish for early craniofacial and cardiac development and this information contributes to the understanding of epigenetic dysregulation during development of Kabuki syndrome. A comparison of high resolution image reconstruction methods from confocal array detector with a new deconvolution algorithm Array detector allows a resolution gain for confocal microscopy by combining images sensed by a set of photomultipliers tubes (or sub-detectors). Several methods have been proposed to reconstruct a high resolution image by linearly combining sub-detector images. To overcome the limitations of these techniques, we propose a new reconstruction method that takes the full stack of spatially reassigned detector signals as input. We show on both calibration slides and real data that our deconvolution method allows to achieve a better reconstruction performance in terms of resolution, image contrast, and spatial intensity homogeneity. The tested algorithms are available in an open source software. Asymmetrical dose-responses shape the evolutionary trade-off between antifungal resistance and nutrient use Antimicrobial resistance is an emerging threat for public health. The success of resistance mutations depends on the trade-off between the benefits and costs they incur. This trade-off is largely unknown and uncharacterized for antifungals. Here, we systematically measure the effect of all amino acid substitutions in the yeast cytosine deaminase Fcy1, the target of the antifungal 5-FC (flucytosine). We identify over 900 missense mutations granting resistance to 5-FC, a large fraction of which appear to act through destabilisation of the protein. The relationship between 5-FC resistance and growth sustained by cytosine deamination is characterized by a sharp trade-off, such that small gains in resistance universally lead to large losses in canonical enzyme function. We show that this steep relationship can be explained by differences in the dose-response functions of 5-FC and cytosine. Finally, we observe the same trade-off shape for the ortholog of FCY1 in Cryptoccocus neoformans, a human pathogen. Our results provide a powerful resource and platform for interpreting drug target variants in fungal pathogens as well as unprecedented insights into resistance-function trade-offs. Specific sequence of arrival promotes coexistence via spatial niche preemption by the weak competitor Historical contingency, such as the order of species arrival, can modify competitive outcomes via niche modification or preemption. However, how these mechanisms ultimately modify stabilising niche and average fitness differences remains largely unknown. By experimentally assembling two congeneric spider mite species feeding on tomato plants during two generations, we show that order of arrival affects species competitive ability and changes the outcome of competition. Contrary to expectations, order of arrival did not cause positive frequency dependent priority effects. Instead, coexistence was predicted when the inferior competitor (Tetranychus urticae) arrived first. In that case, T. urticae colonised the preferred feeding stratum (leaves) of T. evansi leading to spatial niche preemption, which equalised fitness and reduced niche differences, driving community assembly to a close-to-neutrality scenario. Our study demonstrates how the order of species arrival and the spatial context of competitive interactions can jointly determine whether species can coexist. Convergent evolution of multiple mutations improves the viral fitness of SARS-CoV-2 variants by balancing positive and negative selection Multiple mutations have been seen to undergo convergent evolution in SARS-CoV-2 variants of concern. One such evolution occurs in Beta, Gamma, and Omicron variants at three amino acid positions K417, E484, and N501 in the receptor binding domain of the spike protein. We examined the physical mechanisms underlying the convergent evolution of three mutations K417T/E484K/N501Y by delineating the individual and collective effects of mutations on binding to angiotensin converting enzyme 2 receptor, immune escape from neutralizing antibodies, protein stability and expression. Our results show that each mutation serves a distinct function that improves virus fitness supporting its positive selection, even though individual mutations have deleterious effects that make them prone to negative selection. Compared to the wild-type, K417T escapes Class 1 antibodies, has increased stability and expression; however, it has decreased receptor binding. E484K escapes Class 2 antibodies; however, it has decreased receptor binding, stability and expression. N501Y increases receptor binding; however, has decreased stability and expression. When these mutations come together, the deleterious effects are mitigated due to the presence of compensatory effects. Triple mutant K417T/E484K/N501Y has increased receptor binding, escapes both Class 1 and Class 2 antibodies, and has similar stability and expression as that of the wild-type. These results show the implications of presence of multiple mutations on virus evolution that enhance viral fitness on different fronts by balancing both positive and negative selection and improves the chances of selection of mutations together.