A bacterial effector counteracts host autophagy by promoting degradation of an autophagy component Beyond its role in cellular homeostasis, autophagy plays anti- and pro-microbial roles in host-microbe interactions, both in animals and plants. One prominent role of anti-microbial autophagy is to degrade intracellular pathogens or microbial molecules, in a process termed xenophagy. Consequently, microbes evolved mechanisms to hijack or modulate autophagy to escape elimination. Although well-described in animals, the extent to which xenophagy contributes to plant-bacteria interactions remains unknown. Here, we provide evidence that Xanthomonas campestris pv. vesicatoria (Xcv) suppresses host autophagy by utilizing type-III effector XopL. XopL interacts with and degrades the autophagy component SH3P2 via its E3 ligase activity to promote infection. Intriguingly, XopL is targeted for degradation by defense-related selective autophagy mediated by NBR1/Joka2, revealing a complex antagonistic interplay between XopL and the host autophagy machinery. Our results implicate plant antimicrobial autophagy in depletion of a bacterial virulence factor and unravels an unprecedented pathogen strategy to counteract defense-related autophagy.
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.
Increased biting rate and decreased Wolbachia density in irradiated mosquitoes BackgroundReleasing considerable numbers of radiation-sterilized males is a promising strategy to suppress mosquito vectors. However, releases may also include small percentages of biting females which translate to large numbers when releases are large. Currently, the effects of irradiation on the host-seeking and host-biting behaviors have not been exhaustively investigated. Information is also lacking regarding the effects of sterilizing treatment on the endosymbiotic bacterium Wolbachia, which is known to affect the vector competence of infected mosquitos.

MethodsTo ascertain the effects of irradiation on females, the pupae of two Aedes albopictus strains, differing in their natural or artificial Wolbachia infection type, and Ae. aegypti--which is not infected by Wolbachia--were treated with various doses of X-rays and monitored for key fitness parameters and biting behavior over a period of two weeks. The effect of radiation on Wolbachia was investigated by qPCR and FISH analysis.

ResultsPartial Ae. albopictus female sterility was achieved at 28 Gy but the number of weekly bites more than doubled compared to that of the controls. Radiation doses of 35 and 45 Gy completely inhibited progeny production but did not significantly affect the survival or flight ability of Ae. albopictus females and caused a tripling of the number of bites per female per week (compared to untreated controls). These results were also confirmed in Ae. aegypti after treatment at 50 Gy. Wolbachia density decreased significantly in 45-Gray-irradiated females, with the greatest decreases in the early irradiation group (26{+/-}2-hour-old pupae). Wolbachia density also decreased as adults aged. This trend was confirmed in ovaries but not in extra-ovarian tissues. FISH analysis showed a strongly reduced Wolbachia-specific fluorescence in the ovaries of 13{+/-}1-day-old females.

ConclusionsThese results suggest that, under SIT programs, the vector capacity of a target population could increase with the frequency of the irradiated females co-released with the sterile males due to an increased biting rate. In the context of a successful suppression, the related safety issues could be generally negligible, but they should be conservatively evaluated when large scale programs relying on imperfect sexing and high overflooding release ratios are run for long time in areas endemic for arboviral diseases. Also, the effects of irradiation on the vector competence deserve further investigation. Maternal Western-Style Diet Impairs Bone Marrow Development and Drives a Hyperinflammatory Phenotype in Hematopoietic Stem and Progenitor Cells in Fetal Rhesus Macaques BackgroundMaternal obesity adversely impacts the in utero metabolic environment and offsprings health, but its effect on fetal hematopoiesis and immune cell development remains incompletely understood, particularly in models that resemble human development.

MethodsWe studied gestational day 130-135 fetuses derived from rhesus macaque dams chronically exposed to a high-fat Western-style diet (WSD) or a low-fat control diet. Fetal immune cell phenotypes and fetal bone marrow architecture and hematopoietic stem and progenitor cell (FBM HSPC) function were examined using bone computed tomography, histology, flow cytometry, single-cell RNA-sequencing, and HSPC transplantation assays.

FindingsMaternal WSD induced premature FBM cavity opening and a codominant increase in the number of FBM adipocytes. Furthermore, a maternal WSD induced a proinflammatory transcriptional response in FBM HSPCs. FBM macrophages from the WSD group exhibited heightened proinflammatory responses to toll-like receptor agonist stimulation. Maternal WSD exposure suppressed the expression of genes required for B-cell development and decreased the frequencies of FBM B-cells. Finally, maternal WSD led to poor engraftment of FBM HSPCs in nonlethally irradiated immunodeficient NOD/SCID/IL2r{gamma}-/- mice.

InterpretationsMaternal WSD impairs FBM development, drives a hyperinflammatory phenotype, and induces functional and differentiation impairment in FBM HSPCs in a translationally relevant nonhuman primate model.

FundingNational Institute of Health

RESEARCH IN CONTEXTO_ST_ABSEvidence before this studyC_ST_ABSMaternal obesity is associated with increased risk of infections and proinflammatory disease in offspring. The translationally-relevant rhesus macaque model was utilized to address the effects of maternal obesogenic diet on fetal hematopoietic and immune cell development.

Added value of this studyWe assessed changes in fetal immune cell phenotypes and fetal hematopoietic stem and progenitor cell function using immunohistochemistry, flow cytometry, single-cell RNA sequencing, and transplantation assays. We determined that chronic consumption of a maternal obesogenic diet induced the development of adipogenic and proinflammatory environments in the fetal bone marrow. Additionally, we detected the impairment in B-cell differentiation program in fetal hematopoietic stem and progenitor cells.

Implications of all the available evidenceThese data demonstrate that maternal obesogenic diet modulates fetal hematopoietic development and could impact the offsprings immune system, including proinflammatory phenotype and a decline in B-cell function. Horizontal transposon transfer and its implications for the ancestral ecology of hydrophiine snakes Transposable elements (TEs), also known as jumping genes, are sequences able to move or copy themselves within a genome. As TEs move throughout genomes they often act as a source of genetic novelty, hence understanding TE evolution within lineages may help in understanding environmental adaptation. Studies into the TE content of lineages of mammals such as bats have uncovered horizontal transposon transfer (HTT) into these lineages, with squamates often also containing the same TEs. Despite the repeated finding of HTT into squamates, little comparative research has examined the evolution of TEs within squamates. Here we examine a diverse family of Australo-Melanesian snakes (Hydrophiinae) to examine if the previously identified, order-wide pattern of variable TE content and activity holds true on a smaller scale. Hydrophiinae diverged from Asian elapids [~]30 Mya and have since rapidly diversified into six amphibious, [~]60 marine and [~]100 terrestrial species which fill a broad range of ecological niches. We find TE diversity and expansion differs between hydrophiines and their Asian relatives and identify multiple HTTs into Hydrophiinae, including three likely transferred into the ancestral hydrophiine from fish. These HTT events provide the first tangible evidence that Hydrophiinae reached Australia from Asia via a marine route. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Fungicide resistance characterised across seven modes of action in a Botrytis cinerea isolated from Australian vineyards BACKGROUNDBotrytis bunch rot, caused by Botrytis cinerea, is an economically important disease of grapes in Australia and across grape growing regions worldwide. Control of this disease relies on canopy management and the application of fungicides. Fungicide application can lead to the selection of resistant B. cinerea populations, which has an adverse effect on management of the disease. Characterising the distribution and severity of resistant B. cinerea populations is needed to inform resistance management strategies.

RESULTSIn this study, 725 isolates were sampled from 75 Australian vineyards during 2013 - 2016 and were screened against seven fungicides with different modes of action (MOAs). The resistance frequencies for azoxystrobin, boscalid, fenhexamid, fludioxonil, iprodione, pyrimethanil and tebuconazole were 5, 2.8, 2.1, 6.2, 11.6, 7.7 and 2.9% respectively. Nearly half of the resistant isolates (43.8%) were resistant to more than one of the fungicides tested. The frequency of vineyards with at least one isolate simultaneously resistant to 1, 2, 3, 4 or 5 fungicides was 19.5, 7.8, 6.5, 10.4 and 2.6%. Resistance was associated with previously published genotypes in CytB (G143A), SdhB (H272R/Y), Erg27 (F412S), Mrr1 (D354Y), Bos1 (I365S, N373S + Q369P, I365S + D757N) and Pos5 (V273I, P319A, L412F/V). Novel genotypes were also described in Mrr1 (S611N, D616G) Pos5 (V273L) and Cyp51 (P347S). Expression analysis was used to characterise fludioxonil resistant isolates exhibiting overexpression (6.3-9.6-fold) of the ABC transporter gene AtrB (MDR1 phenotype).

CONCLUSIONResistance frequencies were lower when compared to most previously published surveys of B. cinerea resistance in grape and other crops. Nevertheless, continued monitoring of critical MOAs used in Australian vineyards is recommended. 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. Sgs1BLM independent role of Dna2DNA2 nuclease at DNA double strand break is inhibited by Nej1XLF The two major pathways of DNA double strand break (DSB) repair, non-homologous end-joining (NHEJ) and homologous recombination (HR), are highly conserved from yeast to mammals. The regulation of 5 DNA resection controls repair pathway choice and influences repair outcomes. Nej1 was first identified as a canonical NHEJ factor involved in stimulating the ligation of broken DNA ends and more recently, it was shown to be important for DNA end-bridging and inhibiting 5 resection mediated by Dna2-Sgs1. Nej1 interacts with Sae2 and this impacts DSB repair in three ways. First, Nej1 inhibits MRX-Sae2 interactions and Sae2 localization to a DSB. Second, Nej1 inhibits Sae2-dependent recruitment of Dna2 in the absence of Sgs1. Third, NEJ1 and SAE2 showed an epistatic relationship for DNA end-bridging, an event that restrains the broken ends and reduces the frequency of genomic deletions from developing at the DSB. Deletion of NEJ1 suppressed the synthetic lethality of sae2{Delta} sgs1{Delta} and was dependent on the nuclease activity of Dna2. These Nej1 functions promote end-joining DSB repair, but could also be relevant for controlling resection initiation during HR repair.

HighlightsO_LINej1 physically interacts with Sae2 and inhibits end-resection at a DSB.
C_LIO_LINej1 inhibits Sae2 interactions with the MRX complex.
C_LIO_LINej1 inhibits Sae2-dependent recruitment of Dna2 to a DSB.
C_LIO_LINEJ1 and SAE2 are epistatic for DNA end-bridging.
C_LIO_LIDeletion of NEJ1 suppresses the synthetic lethality of sae2{Delta} sgs1{Delta}, which is dependent on Dna2 nuclease activity.
C_LI Improved prediction of protein-protein interactions using AlphaFold2 Predicting the structure of interacting protein chains is a fundamental step towards understanding protein function. Unfortunately, no computational method can produce accurate structures of protein complexes. AlphaFold2, has shown unprecedented levels of accuracy in modelling single chain protein structures. Here, we apply AlphaFold2 for the prediction of heterodimeric protein complexes. We find that the AlphaFold2 protocol together with optimized multiple sequence alignments, generate models with acceptable quality (DockQ[&ge;]0.23) for 63% of the dimers. From the predicted interfaces we create a simple function to predict the DockQ score which distinguishes acceptable from incorrect models as well as interacting from non-interacting proteins with state-of-art accuracy. We find that, using the predicted DockQ scores, we can identify 51% of all interacting pairs at 1% FPR. The protocol can be found at: https://gitlab.com/ElofssonLab/FoldDock. Deep inference of seabird dives from GPS-only records: performance and generalization properties AO_SCPLOWBSTRACTC_SCPLOWAt-sea behaviour of seabirds have received significant attention in ecology over the last decades as it is a key process in the ecology and fate of these populations. It is also, through the position of top predator that these species often occupy, a relevant and integrative indicator of the dynamics of the marine ecosystems they rely on. Seabird trajectories are recorded through the deployment of GPS, and a variety of statistical approaches have been tested to infer probable behaviours from these location data. Recently, deep learning tools have shown promising results for the segmentation and classification of animal behaviour from trajectory data. Yet, these approaches have not been widely used and investigation is still needed to identify optimal network architecture and to demonstrate their generalization properties. From a database of about 300 foraging trajectories derived from GPS data deployed simultaneously with pressure sensors for the identification of dives, this work has benchmarked deep neural network architectures trained in a supervised manner for the prediction of dives from trajectory data. It first confirms that deep learning allows better dive prediction than usual methods such as Hidden Markov Models. It also demonstrates the generalization properties of the trained networks for inferring dives distribution for seabirds from other colonies and ecosystems. In particular, convolutional networks trained on Peruvian boobies from a specific colony show great ability to predict dives of boobies from other colonies and from distinct ecosystems. We further investigate accross-species generalization using a transfer learning strategy known as  fine-tuning. Starting from a convolutional network pre-trained on Guanay cormorant data reduced by two the size of the dataset needed to accurately predict dives in a tropical booby from Brazil. We believe that the networks trained in this study will provide relevant starting point for future fine-tuning works for seabird trajectory segmentation. To make easier such applications, we share online the different models we trained on the considered datasets (https://github.com/AmedeeRoy/BirdDL/models).

Author SummaryOver the last decades, the use of miniaturized electronic devices enabled the tracking of many wide-ranging animal species. The deployment of GPS has notably informed on migratory, habitat and foraging strategies of numerous seabird species. A key challenge in movement ecology is to identify specific behavioural patterns (e.g. travelling, resting, foraging) through the observed movement data. In this work, we address the inference of seabird diving behaviour from GPS data using deep learning methods. We demonstrate the performance of deep networks to accurately identify movement patterns from GPS data over state-of-the-art tools, and we illustrate their great accross-species generalization properties (i.e. the ability to generalize prediction from one seabird species to aother). Our results further supports the relevance of deep learning schemes as  ready-to-use tools which could be used by ecologists to segmentate animal trajectories on new (small) datasets, including when these datasets do not include groundtruthed labelled data for a supervised training. Mechanistic basis for understanding the dual activities of the bifunctional Azotobacter vinelandii mannuronan C-5 epimerase and alginate lyase AlgE7 The functional properties of alginates are dictated by the monomer composition and molecular weight distribution. Mannuronan C-5 epimerases determine the former by epimerizing {beta}-D-mannuronic acid residues (M) into -L-guluronic acid residues (G). The molecular weight is affected by alginate lyases, which cleave alginate chains through {beta}-elimination. The reaction mechanisms for the epimerization and cleavage are similar and some enzymes can perform both. These dualistic enzymes share high sequence identity with mannuronan C-5 epimerases without lyase activity, and the mechanism behind their activity as well as the amino acids responsible for it are still unknown. In this study, we investigate mechanistic determinants of the bifunctional epimerase and lyase activity of AlgE7 from Azotobacter vinelandii. Based on sequence analyses, a range of AlgE7 variants were constructed and subjected to activity assays and product characterization by NMR. Our results show that the lyase activity of AlgE7 is regulated by the type of ion present: Calcium promotes it, whereas NaCl reduces it. By using defined poly-M and poly-MG substrates, the preferred cleavage sites of AlgE7 were found to be M{downarrow}XM and G{downarrow}XM, where X can be either M or G. By studying AlgE7 mutants, R148 was identified as an important residue for the lyase activity, and the point mutant R148G resulted in an enzyme with only epimerase activity. Based on the results obtained in the present study we suggest a unified catalytic reaction mechanism for both epimerase and lyase activities where H154 functions as the catalytic base and Y149 as the catalytic acid.

ImportancePost-harvest valorisation and upgrading of algal constituents is a promising strategy in the development of a sustainable bioeconomy based on algal biomass. In this respect, alginate epimerases and lyases are valuable enzymes for tailoring of the functional properties of alginate, a polysaccharide extracted from brown seaweed with numerous applications in food, medicine and material industries. By providing a better understanding of the reaction mechanism and of how the two enzyme reactions can be altered by changes in reaction conditions, this study opens for further applications of bacterial epimerases and lyases in enzymatic tailoring of alginate polymers. Polycystins are upregulated in fibrotic kidneys and promote renal fibrosis Autosomal polycystic kidney disease (ADPKD) is caused by mutations in PKD1or PKD2 gene, which encodes polycystin-1 or polycystin-2 respectively. Paradoxically, polycystins are up-regulated in human ADPKD kidneys and overexpression of polycystin-1 even induced renal cyst formation in mice, suggesting that up-regulation of polycystin is pathogenic. In this study, we determined the role of polycystin in renal fibrosis and explored its underlying mechanisms. We found that the expression of polycystin-1 or polycystin-2 was positively correlated with the expression of collagen-I in the kidney of three different fibrotic mouse models. Pharmaceutical inhibition or genetic knockout of polycystin-2 reduced the expression of pro-fibrotic markers in fibrotic kidneys of two animal models. Moreover, conditional knockout of Pkd1 gene attenuated renal fibrosis in two mouse models. Enhancer of Zeste homolog 2 (EZH2) is a methyltransferase which plays important role in renal fibrosis. We found that EZH2 was positively correlated with the expression of polycystins in fibrotic kidneys. Conditional knockout of the Ezh2 attenuated the anti-fibrotic responses induced by Pkd1 deletion in fibrotic kidneys.

In conclusion, polycystins are up-regulated in injured kidneys and polycystins promote renal fibrosis through up-regulation of EZH2 in fibrotic kidneys, suggesting that renal fibrosis in polycystic kidneys is not the direct effect of PKD gene mutations.

SignificanceMutations of polycystins cause autosomal polycystic kidney disease (ADPKD) which is associated with increased renal fibrosis in cystic kidneys. However, the role of polycystin in renal fibrosis remains unclear. Here, we showed that the expression of polycystins is positively correlated with the development of renal fibrosis in injured kidneys. Pharmaceutical inhibition or genetic deletion of polycystin-2 attenuated renal fibrosis in fibrotic kidneys. Moreover, genetic deletion of polycystin-1 inhibited renal fibrosis in two different animal models. Finally, we showed that ablation of the methyltransferase EZH2 mitigated the anti-fibrotic effect of polycystin-1 deletion in fibrotic kidneys. We conclude that polycystin promotes through EZH2 and renal fibrosis in cystic kidney is not the direct effect of PKD gene mutations. Agonist efficiency from concentration-response curves: structural implications and applications Agonists are evaluated by a concentration-response curve (CRC), with a midpoint (EC50) that indicates potency, a high-concentration asymptote that indicates efficacy and a low-concentration asymptote that indicates constitutive activity. A third agonist attribute, efficiency ({eta}), is the fraction of binding energy that is applied to the conformational change that activates the receptor. We show that {eta} can be calculated from EC50 and the asymptotes of a CRC derived from either single-channel or whole-cell responses. For 20 agonists of skeletal muscle nicotinic receptors, the distribution of {eta} values is bimodal with population means at 51% (including ACh, nornicotine and DMPP) and 40% (including epibatidine, varenicline and cytisine). The value of {eta} is related inversely to the size of the agonists head-group, with high-versus low-efficiency ligands having an average volume of 70 [A]3 versus 102 [A]3. Most binding site mutations have only a small effect on ACh efficiency except for Y190A (35%), W149A (60%) and those at G153 (42%). If {eta} is known, the midpoint and high-concentration asymptote can be calculated from each other. Hence, an entire CRC can be estimated from the response to a single agonist concentration, and efficacy can be estimated from EC50 of a CRC that has been normalized to 1. Given {eta}, the level of constitutive activity can be estimated from a single CRC.

Statement of significanceReceptors are molecular machines that convert chemical energy from agonist binding into mechanical energy of a global conformational change that generates a cell response. Agonists are distinguished by their potency (proportional to affinity) and efficacy, but also by the efficiency at which their binding energy is applied to receptor activation. Here we show that agonist efficiency can be estimated from a single concentrationresponse curve, and estimate efficiencies of 20 nicotinic receptor agonists. These have a bimodal distribution with larger agonists belonging to the lower-efficiency population. We further show that mutations of some binding site residues alter efficiency, and that knowledge of agonist efficiency simplifies and extends dose-response curve analysis. Effects of substrate on shrimp growth, water quality and bacterial community in the biofloc system nursing Penaeus vannamei postlarvae This study aimed to investigate the effects of substrate on water quality, shrimp growth and bacterial community in the biofloc system with a salinity of 5{per thousand}. Two treatments, biofloc system with (sB) or without (nB) addition of elastic solid packing filler (nylon) as substrate, were set up. Penaeus vannamei postlarvae (PL, [~] stage 15) were stocked at a density of 4000 PL m-3 in each treatment with triplicates for a 28-days culture experiment, taking glucose as carbon source (C:N 15:1). Results showed that the survival rate (96.3{+/-}3.6%), FCR (0.76{+/-}0.06) and productivity (1.54{+/-}0.12 kg m-3) in sB treatment were significantly better than those in nB treatment (81.0{+/-}7.1%, 0.98{+/-}0.08 and 1.14{+/-}0.09 kg m-3, P<0.05). All water parameters were in the recommended ranges. Substrate showed significant effect on TAN, TSS, turbidity, biofloc volume, pH and carbonate alkalinity (P < 0.05). Actinobacteria (4.0-22.7%), Bacteroidetes (10.4-33.5%), Firmicutes (0.2-11.2%), Planctomycetes (4.0-14.9%) and Proteobacteria (29.4-59.0%) were the most dominant phyla for both treatments. However, the bacterial community in sB treatment showed to be significantly different from that in nB treatment (Jaccard distance 0.94{+/-}0.01, P=0.001). Substrate showed significant effects on Shannon, Heip, Pielou and Simpson index, as well as relative abundances of Actinobacteria, Bacteroidetes and Planctomycetes (P < 0.05). The results suggested that addition of substrate affected the shrimp growth, water quality and bacterial community in the biofloc system nursing P. vannamei PL with a 5{per thousand} salinity. CTP and parS coordinate ParB partition complex dynamics and ParA-ATPase activation for ParABS-mediated DNA partitioning ParABS partition systems, comprising the centromere-like DNA sequence parS, the parS-binding ParB-CTPase and the nucleoid-binding ParA-ATPase, ensure faithful segregation of bacterial chromosomes and low-copy-number plasmids. F-plasmid partition complexes containing ParBF and parSF move by generating and following a local concentration gradient of nucleoid-bound ParAF. However, the process through which ParBF activates ParAF-ATPase has not been defined. We studied CTP- and parSF-modulated ParAF--ParBF complex assembly, in which DNA-bound ParAF-ATP dimers are activated for ATP hydrolysis by interacting with two ParBF N-terminal domains. CTP or parSF enhances the ATPase rate without significantly accelerating ParAF--ParBF complex assembly. Together, parSF and CTP accelerate ParAF--ParBF assembly without further significant increase in ATPase rate. Magnetic-tweezers experiments showed that CTP promotes multiple ParBF loading onto parSF-containing DNA, generating condensed partition complex-like assemblies. We propose that ParBF in the partition complex adopts a conformation that enhances ParBF--ParBF and ParAF--ParBF interactions promoting efficient partitioning. A General Biphasic-Bodyweight Model for Scaling Basal Metabolic Rate, Glomerular Filtration Rate and Drug Clearance from Birth to Adulthood Understanding the maturation process of human physiology and metabolism has broad medical and pharmaceutical implications. Age and bodyweight are frequently considered as separate variables in modeling the dynamical changes of human organ functions and of drug clearance from birth to adulthood. The objective of this study is to propose a unified, continuous and bodyweight-only equation to quantify the changes of human basal metabolic rate (BMR), glomerular filtration rate (GFR) and drug clearance (CL) from infancy to adulthood. The BMR datasets were retrieved from a comprehensive historical database of male and female subjects (0.02 to 64 years). The CL datasets for 17 drugs and the GFR dataset were generated from age-incorporated maturation-and-growth models with reported parameter values. The model used in the simulation is independent of the proposed model. A statistical approach was used to simulate the model generated CL and GFR data for a hypothetical population with 26 age groups (ranging from 0 to 20 years). Besides, individual CL data for one drug, and sparse PBPK-modeled CL values for two drugs were also included for analysis. A 4-parameter, mixed-allometry equation with two power-law functions of bodyweight was proposed and evaluated as a general model using nonlinear regression and dimensionless analysis. All datasets universally reveal biphasic curves with two distinct linear segments on log-log plots. Compared with simple allometry, the biphasic model fits satisfactorily to all datasets (based on Akaikes Information Criterion and residual plots). The biphasic equation consists of two reciprocal allometric terms that asymptotically determine the overall curvature. The fitting results show a superlinear scaling phase (slope >1; ca. 1.5 - 3.5) below the characteristic bodyweight at the phase transition; and above which, a sublinear scaling phase (slope <1; ca. 0.5 - 0.7) is evident. The phase-transition bodyweight is ranging from 5 to 20 kg (corresponding to 0.5 - 9 years) and the mean value is around 10 kg ([~]2 years) for all data sets. The dimensionless analysis generalizes, and offers quantitative realization of, the maturation and growth process. In conclusion, the proposed mixed-allometry equation is a generic model that quantitatively describes the phase transition occurring in the human maturation process of BMR, GFR and drug CL. 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. Patterned apoptosis has an instructive role for local growth and tissue shape in a fast-growing epithelium What regulates organ size and shape remains one of the fundamental mysteries of modern biology. So far, research in this area has primarily focused on deciphering the regulation in time and space of growth and cell division, while the contribution of cell death has been much more neglected. This includes studies of the Drosophila wing imaginal disc, the prospective fly wing which undergoes massive growth during larval stage, and represents one of the best characterised systems for the study of growth and patterning. So far, it has been assumed that cell death was relatively neglectable in this tissue and as a result the pattern of growth was usually attributed to the distribution of cell division. Here, using systematic mapping and registration combined with quantitative assessment of clone size and disappearance, we show for the first time that cell death is not neglectable, and outline a persistent pattern of cell death and clone elimination in the disc. Local variation of cell death is associated with local variation of clone size, pointing to an impact of cell death on local growth which is not fully compensated by proliferation. Using morphometric analyses of adult wing shape and genetic perturbations, we provide evidence that patterned death affects locally and globally adult wing shape and size. This study describes a roadmap for accurate assessment of the contribution of cell death to tissue shape, and outlines for the first time an important instructive role of cell death in modulating quantitatively local growth and the morphogenesis of a fast-growing tissue. Genetic Loci and Metabolic States Associated With Murine Epigenetic Aging Changes in DNA methylation (DNAm) are linked to aging. Here, we profile highly conserved CpGs in 339 predominantly female mice belonging to the BXD family for which we have deep longevity and genomic data. We use a  pan-mammalian microarray that provides a common platform for assaying the methylome across mammalian clades. We computed epigenetic clocks and tested associations with DNAm entropy, diet, weight, metabolic traits, and genetic variation. We describe the multifactorial variance of methylation at these CpGs, and show that high fat diet augments the age-associated changes. Entropy increases with age. The progression to disorder, particularly at CpGs that gain methylation over time, was predictive of genotype-dependent life expectancy. The longer-lived BXD strains had comparatively lower entropy at a given age. We identified two genetic loci that modulate rates of epigenetic age acceleration (EAA): one on chromosome (Chr) 11 that encompasses the Erbb2/Her2 oncogenic region, and a second on Chr19 that contains a cytochrome P450 cluster. Both loci harbor genes associated with EAA in humans including STXBP4, NKX2-3, and CUTC. Transcriptome and proteome analyses revealed associations with oxidation-reduction, metabolic, and immune response pathways. Our results highlight concordant loci for EAA in humans and mice, and demonstrate a tight coupling between the metabolic state and epigenetic aging. 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. High performance single-cell gene regulatory network inference atscale: The Inferelator 3.0 MotivationGene regulatory networks define regulatory relationships between transcription factors and target genes within a biological system, and reconstructing them is essential for understanding cellular growth and function. Methods for inferring and reconstructing networks from genomics data have evolved rapidly over the last decade in response to advances in sequencing technology and machine learning. The scale of data collection has increased dramatically; the largest genome-wide gene expression datasets have grown from thousands of measurements to millions of single cells, and new technologies are on the horizon to increase to tens of millions of cells and above.

ResultsIn this work, we present the Inferelator 3.0, which has been significantly updated to integrate data from distinct cell types to learn context-specific regulatory networks and aggregate them into a shared regulatory network, while retaining the functionality of the previous versions. The Inferelator is able to integrate the largest single-cell datasets and learn cell-type specific gene regulatory networks. Compared to other network inference methods, the Inferelator learns new and informative Saccharomyces cerevisiae networks from single-cell gene expression data, measured by recovery of a known gold standard. We demonstrate its scaling capabilities by learning networks for multiple distinct neuronal and glial cell types in the developing Mus musculus brain at E18 from a large (1.3 million) single-cell gene expression dataset with paired single-cell chromatin accessibility data.

AvailabilityThe inferelator software is available on GitHub (https://github.com/flatironinstitute/inferelator) under the MIT license and has been released as python packages with associated documentation (https://inferelator.readthedocs.io/). Steady-state measures of visual suppression In the early visual system, suppression occurs between neurons representing different stimulus properties. This includes features such as orientation (cross-orientation suppression), eye-of-origin (interocular suppression) and spatial location (surround suppression), which are thought to involve distinct anatomical pathways. We asked if these separate routes to suppression can be differentiated by their pattern of gain control on the contrast response function measured in human participants using steady-state electroencephalography. Changes in contrast gain shift the contrast response function laterally, whereas changes in response gain scale the function vertically. We used a Bayesian hierarchical model to summarise the evidence for each type of gain control. A computational meta-analysis of 16 previous studies found the most evidence for contrast gain effects with overlaid masks, but no clear evidence favouring either response gain or contrast gain for other mask types. We then conducted two new experiments, comparing suppression from four mask types (monocular and dichoptic overlay masks, and aligned and orthogonal surround masks) on responses to sine wave grating patches flickering at 5Hz. At the occipital pole, there was strong evidence for contrast gain effects in all four mask types at the first harmonic frequency (5Hz). Suppression generally became stronger at more lateral electrode sites, but there was little evidence of response gain effects. At the second harmonic frequency (10Hz) suppression was stronger overall, and involved both contrast and response gain effects. Although suppression from different mask types involves distinct anatomical pathways, gain control processes appear to serve a common purpose, which we suggest might be to suppress less reliable inputs. Extending Bayesian Elo-rating to quantify the steepness of dominance hierarchies The steepness of dominance hierarchies provides information about the degree of competition within animal social groups and is thus an important concept in socioecology. The currently most widely-used metrics to quantify steepness are based on Davids scores (DS) derived from dominance interaction networks. One serious drawback of these DS-based metrics is that they are biased, i.e., network density systematically decreases steepness values. Here, we provide a novel approach to estimate steepness based on Elo-ratings, implemented in a Bayesian framework (STEER: Steepness estimation with Elo-rating). Our new metric has two key advantages. First, STEER is unbiased, precise and more robust to data density than DS-based steepness. Second, it provides explicit probability distributions of the estimated steepness coefficient, which allows uncertainty assessment. In addition, it relies on the same underlying concept and is on the same scale as the original measure, and thus allows comparison to existing published results. We evaluate and validate performance of STEER by means of experimentation on empirical and artificial data sets and compare its performance to that of several other steepness estimators. Our results suggest that STEER provides a considerable improvement over existing methods to estimate dominance hierarchy steepness. We provide an R package EloSteepness to calculate the new steepness measure, and also show an example of using steepness in a comparative analysis. Intermittent compression increases flow through a gel: a mechanical model perhaps relevant to flow of cerebrospinal fluid Exchange of molecules between cerebrospinal fluid (CSF) and brain cells contributes to brain function and protection from dementia. Despite widespread acceptance of the  glymphatic theory, the route by which CSF is brought close enough to the neural tissue for solutes to be exchanged by extracellular diffusion is not entirely clear. Exogenous molecules injected into CSF are observed to reach the basement lamina that surrounds the dense capillary network. Transport of solutes by diffusion along the basement lamina, a gel of macromolcules about 100 nm thick, would be too slow; bulk flow in a static geometry would require unphysiologically high pressures. However, it is known that the pulsation of blood aids transport of CSF, and we hypothesized that this is because the pulsation intermittently squeezes the pericapillary lamina. In a primitive mimicry, we have tested whether intermittent squeezing increases flow through an agar gel. In all but one of 216 tests, pulsation caused a reversible increase, sometimes by a factor of 100 or more. The enhancement was greatest for frequencies 5-11 Hz and, over the tested range of pressure heads (20 - 50 cmH2O), was greatest for the lowest pressure. The results suggest a reason why exercise slows the aging of the brain. 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.