Living with high potassium: a balance between nutrient acquisition and stress signaling during K-induced salt stress O_LIHigh potassium (K) in the growth medium is more toxic to plants than Na at similar concentrations. However, the molecular mechanisms underlying plant responses to K-induced salt stress are virtually unknown.
C_LIO_LIWe examined Arabidopsis thaliana and its extremophyte relative Schrenkiella parvula, using a comparative multi-omics approach to identify cellular processes affected by excess K and understand which deterministic regulatory pathways are active to avoid tissue damage while sustaining growth.
C_LIO_LIA. thaliana showed limited capacity to curb excess K accumulation and prevent nutrient depletion contrasting to S. parvula which could limit excess K accumulation without restricting nutrient uptake. Facilitated by a targeted transcriptomic response, promoting nitrogen uptake along with other key nutrients and uninterrupted N assimilation into primary metabolites during excess K-stress allowed S. parvula to boost its antioxidant and osmolyte pools concurrently leading to sustained growth. Antithetically, A. thaliana showed transcriptional responses indicative of a poor balance between stress signaling, increased ROS levels, and reduced photosynthesis, subsequently leading to inhibited growth.
C_LIO_LIThe ability to regulate independent nutrient uptake and a coordinated transcriptomic response to avoid non-specific stress signaling are two main deterministic steps towards building stress resilience to excess K+-induced salt stress.
C_LI
Isoprenoid biosynthesis regulation in poplars by methylerythritol phosphate and mevalonic acid pathways The isoprenoids found in plants are extremely important to survive with various human applications, such as flavoring, fragrance, dye, pharmaceuticals, and biomass used for biofuels. Methylerythritol phosphate (MEP) and mevalonic acid (MVA) pathways are critical in plants, responsible for isoprenoid biosynthesis. 1-deoxy-D-xylulose5-phosphate synthase (DXS) and 1-deoxy-D-xylulose5-phosphate reductoisomerase (DXR) catalyze the rate-limiting steps in the MEP pathway, while 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) catalyzes the rate-limiting step in the MVA pathway. Here, we showed while PtHMGR overexpressors (OEs) exhibited different MEP- and MVA-related gene expressions compared with non-transgenic poplars (NT), the PtDXR-OEs revealed upregulated MEP-related and downregulated MVA-related gene expressions. PtDXR and PtHMGR overexpressions caused changes in MVA-derived trans-zeatin-riboside, isopentenyl adenosine, castasterone, and 6-deoxocastasterone well as MEP-derived carotenoids and gibberellins. In PtHMGR-OEs, the accumulated geranyl diphosphate synthase (GPS) and geranyl pyrophosphate synthase (GPPS) transcript levels in the MEP pathway led to an accumulation of MEP-derived isoprenoids. In contrast, upregulation of farnesyl diphosphate synthase (FPS) expression in the MVA pathway contributed to increased levels of MVA-derived isoprenoids. In addition, PtHMGR-OEs increased MEP-related GPS and GPPS transcript levels, expanded MEP-derived isoprenoid levels, changed FPS transcript levels, and affected MVA-derived isoprenoid yields. These results demonstrate the contribution of MVA and MEP pathways regulating isoprenoid biosynthesis in poplars.
High-plex Multiomic Analysis in FFPE at Subcellular Level by Spatial Molecular Imaging The Spatial Molecular Imaging platform (CosMxTM SMI, NanoString, Seattle, WA) utilizes a high-plex in-situ imaging chemistry for both RNA and protein detection. This automated instrument provides 1000s of plex, at high sensitivity (1 to 2 copies/cell), very low error rate (0.0093 false calls/cell), and background ([~]0.04 counts/cell). The imaging system generates three-dimensional super-resolution localization of analytes (< 50 nm in the X-Y plane), at [~] 1 million cells per sample, four samples per run. Cell segmentation is antibody-based, compatible with FFPE samples. Multiomic data (980 RNAs, 108 proteins) were measured at subcellular resolution using FFPE tissues (non-small cell lung (NSCLC) and breast cancer) and allowed identification of over 18 distinct cell types, 10 unique tumor microenvironments, and 100 pairwise ligand-receptor interactions. Over 800,000 single cells and [~]260 million transcripts data are released into the public domain allowing extended data analysis by the entire spatial biology research community. Metabolic crossfeeding structures the assembly of polysaccharide degrading communities Metabolic processes that fuel the growth of heterotrophic microbial communities are initiated by specialized biopolymer degraders that decompose complex forms of organic matter. It is unclear, however, to what extent degraders control the downstream assembly of the community that follows polymer breakdown. Investigating a model marine microbial community that degrades chitin, we show that chitinases secreted by different degraders produce oligomers of specific chain lengths that not only select for specialized consumers but also influence the metabolites secreted by these consumers into a shared resource pool. Each species participating in the breakdown cascade exhibits unique hierarchical preferences for substrates, which underlies the sequential colonization of metabolically distinct groups as resource availability changes over time. By identifying the metabolic underpinnings of microbial community assembly, we reveal a hierarchical crossfeeding structure that allows biopolymer degraders to shape the dynamics of community assembly.

One sentence summarySpecialized biopolymer degraders direct the trajectory of microbial community assembly through interconnected modes of nutrient crossfeeding. The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. A platform of patient-derived microtumors identifies individual treatment responses and therapeutic vulnerabilities in ovarian cancer. BackgroundIn light of the frequent development of therapeutic resistance in cancer treatment, there is a strong need for personalized model systems representing patient tumor heterogeneity, while enabling parallel drug testing and identification of appropriate treatment responses in individual patients. Using ovarian cancer as a prime example of a heterogeneous tumor disease with complex microenvironment and high recurrence rates, we developed a 3D preclinical ovarian cancer model comprised of patient-derived microtumors (PDM) and autologous tumor-infiltrating lymphocytes (TILs) for identification of treatment vulnerabilities and validation of anti-cancer drug efficacy using immunohistochemistry, immune cell phenotyping, functional assays and protein profiling analyses.

MethodsPDM and TILs were isolated from fresh primary ovarian cancer tissue specimen using mechanical disruption and limited enzymatic digestion and were subsequently cultured in suspension in defined media in the absence of serum. The heterogeneous cellular composition of isolated PDM as well as autologous TILs was analyzed by FFPE immunohistochemistry and multi-color flow cytometry, respectively. For in-depth protein profiling of PDM we established Reverse Phase Protein Array (RPPA) analyses of >110 total and phospho-proteins. Treatment efficacy in response to chemotherapeutics as well as immunotherapeutic compounds was assessed in PDM and PDM-TIL co-cultures using a functional viability assay in microplate format.

ResultsThe enzymatic digestion of primary ovarian cancer tissue and suspension culture in defined serum-free media allowed fast and efficient recovery of patient-derived microtumors (PDM). Immunohistochemical analyses demonstrated histopathological comparability of ovarian cancer PDM with corresponding patient tumor tissue. Reverse Phase Protein Array (RPPA)-based analyses of >110 total and phospho-proteins enabled the identification of patient-specific sensitivities to standard, platinum-based therapy and thereby the prediction of potential treatment-responders. Finally, combining PDM and autologous TILs for individual efficacy testing of immune checkpoint inhibitor treatment demonstrated the potential for patient-specific enhancement of cytotoxic TIL activity by this therapeutic approach.

Conclusion3D patient-derived ovarian cancer microtumors represent a preclinical, ex vivo tumor model that reflects intertumoral heterogeneity and represent the cellular complexity of individual patient tumors. Combining protein pathway analysis and anti-cancer drug efficacy testing of PDM enables drug mode-of-action analyses and therapeutic sensitivity prediction within a clinically relevant time frame after surgery. Follow-up studies in larger cohorts are currently under way to further evaluate the applicability of this platform to support clinical decision-making and personalizing cancer treatment. New Segmentation and Feature Extraction Algorithm for Classification of White Blood Cells in Peripheral Smear Images This article addresses a new method for the classification of white blood cells (WBCs) using image processing techniques and machine learning methods. The proposed method consists of three steps: detecting the nucleus and cytoplasm, extracting features, and classification. At first, a new algorithm is designed to segment the nucleus. For the cytoplasm to be detected, only a part of it located inside the convex hull of the nucleus is involved in the process. This attitude helps us overcome the difficulties of segmenting the cytoplasm. In the second phase, three shapes and four novel color features are devised and extracted. Finally, by using an SVM model, the WBCs are classified. The segmentation algorithm can detect the nucleus with a dice similarity coefficient of 0.9675. The proposed method can categorize WBCs in Raabin-WBC, LISC, and BCCD datasets with accuracies of 94.65 %, 92.21 %, and 94.20 %, respectively. It is worth mentioning that the hyperparameters of the classifier are fixed only with the Raabin-WBC dataset, and these parameters are not readjusted for LISC and BCCD datasets. The obtained results demonstrate that the proposed method is robust, fast, and accurate. Genetic targeting of Card19 is linked to disrupted Ninj1 expression, impaired cell lysis, and increased susceptibility to Yersinia infection Cell death plays a critical role in inflammatory responses. During pyroptosis, inflammatory caspases cleave Gasdermin D (GSDMD) to release an N-terminal fragment that generates plasma membrane pores that mediate cell lysis and IL-1 cytokine release. Terminal cell lysis and IL-1{beta} release following caspase activation can be uncoupled in certain cell types or in response to particular stimuli, a state termed hyperactivation. However, the factors and mechanisms that regulate terminal cell lysis downstream of GSDMD cleavage remain poorly understood. In the course of studies to define regulation of pyroptosis during Yersinia infection, we identified a line of Card19-deficient mice (Card19lxcn) whose macrophages were protected from cell lysis and showed reduced apoptosis and pyroptosis, yet had wild-type levels of caspase activation, IL-1 secretion, and GSDMD cleavage. Unexpectedly, CARD19, a mitochondrial CARD-containing protein, was not directly responsible for this, as two independently-generated CRISPR/Cas9 Card19 knockout mice showed no defect in macrophage cell lysis, and expression of CARD19 in Card19lxcn macrophages did not restore cell lysis. Card19 is located on chromosome 13, adjacent to Ninj1, which was recently reported to regulate cell lysis downstream of GSDMD activation. Intriguingly, RNA-seq and western blotting revealed that Card19lxcn BMDMs are hypomorphic for NINJ1 expression, and reconstitution of Ninj1 in Card19lxcn immortalized BMDMs restored cell lysis. Card19lxcn mice exhibited significantly increased susceptibility to Yersinia infection, demonstrating that cell lysis itself plays a key role in protection against bacterial infection. Our findings identify genetic targeting of Card19 being responsible for off-target effects on the adjacent Ninj1 gene, thereby disrupting the ability of macrophages to undergo plasma membrane rupture downstream of gasdermin cleavage and impacting host survival and bacterial control during Yersinia infection.

Author SummaryProgrammed cell death is critical for regulating tissue homeostasis and host defense against infection. Pyroptosis is an inflammatory form of programmed cell death that couples cell lysis with release of inflammatory cytokines. Cell lysis is triggered by activation of particular intracellular pore forming proteins, but how regulation of cell lysis occurs is not well understood. Genetic targeting of Card19 on chromosome 13 resulted in decreased expression of the adjacent gene, Ninj1 which was recently found to regulate terminal lysis events in response to cell death-inducing stimuli. We found that macrophages from Card19-deficient mice were resistant to multiple forms of cell death in response to a variety of inflammatory stimuli, including canonical and non-canonical inflammasome activation, as well as triggers of cell-extrinsic apoptosis. Notably, Card19-deficient mice were more susceptible to Yersinia infection, indicating that cell lysis contributes to control of bacterial infections. Our data provide new insight into the impact of terminal cell lysis on control of bacterial infection and highlight the role of additional factors that regulate lytic cell death downstream of gasdermin cleavage. Effect of tag attachment on the flight performance of five raptor species Bio-logging devices play a fundamental and indispensable role in movement ecology studies, particularly in the wild. However, researchers are becoming increasingly aware of the influence that attaching devices can have on animals, particularly on their behaviour, energy expenditure and survival. The way a device is attached to an animals body has also potential consequences for the collected data, and quantifying the type and magnitude of such potential effects is fundamental to enable researchers to combine and compare data from different studies, as much as it is to improve animal welfare.

For over two decades, large terrestrial birds have been in the focus of long-term movement ecology research, employing bio-logging devices attached with different types of harnesses. However, comparative studies investigating the effects of different harness types used on these species are scarce.

In this study, we tested for potential differences in data collected by two commonly used harness types, backpack and leg-loop, on the flight performance of 10 individuals from five raptor species, equipped with high resolution bio-logging devices, in the same area and time. We explored the effect of harness type on vertical speed, horizontal speed, glide ratio, height above sea level, distance travelled, proportion of soaring and flapping behaviour, and VeDBA (a proxy for energy expenditure) between and within individuals, all used as fine-scale measures of flight performance. Birds equipped with leg-loops climbed up to 0.65 ms-1 faster, reached 19% greater heights while soaring, and spent less time in active flight compared to birds equipped with backpacks, suggesting that backpack harnesses, compared to leg-loops, might cause additional drag affecting the birds flight performance. A lower rate of sinking while gliding, a slightly higher glide ratio, higher horizontal speed while soaring, and lower VeDBA, were also indicative of less drag using leg-loops. Our results add to the existing literature highlighting the design-related advantages of leg-loops, and support the use of leg-loops as a better alternative to backpack harnesses for large soaring birds, when possible. Our study also highlights how apparently small changes in device attachment can lead to notable improvements in tagging practice, with implications for animal welfare, data interpretation and comparability. Sahul's megafauna were vulnerable to plant-community changes due to their position in the trophic network Extinctions stemming from environmental change often trigger trophic cascades and coextinctions. However, it remains unclear whether trophic cascades were a large contributor to the megafauna extinctions that swept across several continents in the Late Pleistocene. The pathways to megafauna extinctions are particularly unclear for Sahul (landmass comprising Australia and New Guinea), where extinctions happened earlier than on other continents. We investigated the role of bottom-up trophic cascades in Late Pleistocene Sahul by constructing pre-extinction (~ 80 ka) trophic network models of the vertebrate community of Naracoorte, south-eastern Australia. These models allowed us to predict vertebrate species vulnerability to cascading extinctions based on their position in the network. We tested whether the observed extinctions could be explained by bottom-up cascades, or if they should be attributed to other external causes. Species that disappeared from the community were more vulnerable, overall, to bottom-up cascades than were species that survived. The position of extinct species in the network - having few or no predators - also suggests they might have been particularly vulnerable to a new predator. These results provide quantitative evidence that trophic cascades and naivety to predators could have contributed to the megafauna extinction event in Sahul. Structural and mechanistic basis of reiterative transcription initiation Reiterative transcription initiation, observed at promoters that contain homopolymeric sequences at the transcription start site, generates RNA products having 5 sequences non-complementary to the DNA template. Here, using crystallography and cryo-EM to define structures, protein-DNA-photocrosslinking to map positions of RNAP leading and trailing edges relative to DNA, and single-molecule DNA nanomanipulation to assess RNAP-dependent DNA unwinding, we show that RNA extension in reiterative transcription initiation: (1) occurs without DNA scrunching, (2) involves a short, 2-3 bp RNA-DNA hybrid, and (3) generates RNA that exits RNAP through the portal by which scrunched nontemplate-strand DNA exits RNAP in standard transcription initiation. The results establish that, whereas RNA extension in standard transcription initiation proceeds through a scrunching mechanism, RNA extension in reiterative transcription initiation proceeds through a slippage mechanism, with slipping of RNA relative to DNA within a short RNA-DNA hybrid, and with extrusion of RNA from RNAP through an alternative RNA exit. Motor cortical output for skilled forelimb movement is selectively distributed across projection neuron classes The interaction of descending neocortical outputs and subcortical premotor circuits is critical for shaping the skilled movements. Two broad classes of motor cortical output projection neurons provide input to many subcortical motor areas: pyramidal tract neurons (PT), which project throughout the neuraxis; and intratelencephalic neurons (IT), which project within cortex and subcortical striatum. It is unclear whether these classes are functionally in series or whether separable components of descending motor control signals are distributed across these distinct classes of projection neurons. Here we combine large-scale neural recordings across all layers of motor cortex with cell-type specific perturbations to study cortically-dependent mouse motor behaviors: kinematically-variable manipulation of a joystick and a kinematically-precise reach-to-grasp. We find that striatum-projecting IT neuron activity preferentially represents amplitude whereas pons-projecting PT neurons preferentially represent the variable direction of forelimb movements. Thus, separable components of descending motor cortical commands are distributed across motor cortical projection cell classes.

One-sentence summarySeparable components of cortical motor commands are distributed across distinct glutamatergic projection neuron cell-types. Standardised method for cardiomyocyte isolation and purification from individual murine neonatal, infant, and adult hearts Primary cardiomyocytes are invaluable for understanding postnatal heart development. However, a universal method to obtain freshly purified cardiomyocytes without using different agedependent isolation procedures and cell culture, is lacking. Here, we report the development of a standardised method that allows rapid isolation and purification of high-quality cardiomyocytes from individual neonatal through to adult C57BL/6J murine hearts. Langendorff retrograde perfusion, which is currently limited to adult hearts, was adapted for use in neonatal and infant hearts by developing an easier in situ aortic cannulation technique. Tissue digestion conditions were optimised to achieve efficient digestion of hearts of all ages in a comparable timeframe (<14 min). This resulted in a high yield (1.56-2.2x106 cells/heart) and viability (~70-100%) of cardiomyocytes post-isolation. An immunomagnetic cell separation step was then applied to yield highly purified cardiomyocytes (~95%) as confirmed by immunocytochemistry, flow cytometry, and qRT-PCR. For cell-type specific studies, cardiomyocyte DNA, RNA, and protein could be extracted in sufficient yields to conduct molecular experiments. We generated transcriptomic datasets for neonatal cardiomyocytes from individual hearts, for the first time, which revealed nine sex-specific genes (FDR<0.05) encoded on the sex chromosomes. Finally, we also developed an in situ fixation protocol that preserved the native cytoarchitecture of cardiomyocytes (~94% rodshaped post-isolation), and used it to evaluate cell morphology during cardiomyocyte maturation, as well as capture spindle-shaped neonatal cells undergoing cytokinesis. Together, these procedures allow molecular and morphological profiling of high-quality cardiomyocytes from individual hearts of any postnatal age. 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. A Wnt/β-catenin signaling driven adipocyte population required for beiging and glucose homeostasis Wnt/{beta}-catenin signaling has been well established as a potent inhibitor of adipogenesis. Here, we identified a population of adipocytes that exhibit persistent activity of Wnt/{beta}-catenin signaling, as revealed by the Tcf/Lef-Gfp reporter allele, in embryonic and adult mouse fat depots, named as Wnt+ adipocytes. We showed that the {beta}-catenin-mediated signaling activation in these cells is Wnt ligand- and receptor-independent but relies on AKT/mTOR pathway and is essential for cell survival. Such adipocytes are distinct from classical ones in transcriptomic and genomic signatures and can be induced from various sources of mesenchymal stromal cells including human cells. Genetic lineage-tracing and targeted cell ablation studies revealed that these adipocytes convert into beige adipocytes directly and are also required for beige fat recruitment under thermal challenge, demonstrating both cell autonomous and non-cell autonomous roles in adaptive thermogenesis. Furthermore, mice bearing targeted ablation of these adipocytes exhibited glucose intolerance, while mice receiving exogenously supplied such cells manifested enhanced glucose utilization. Our studies uncover a unique adipocyte population in regulating beiging in adipose tissues and systemic glucose homeostasis. CEP162 deficiency causes human retinal degeneration and reveals a dual role inciliogenesis and neurogenesis Defects in primary or motile cilia result in a variety of human pathologies, and retinal degeneration is frequently associated with these so-called ciliopathies. We show that homozygosity for a truncating variant in CEP162, a centrosome and microtubule-associated protein required for transition zone (TZ) assembly during ciliogenesis and neuronal differentiation in the retina, causes late-onset retinitis pigmentosa in 2 unrelated families. The mutant CEP162-E646R*5 protein is expressed and properly localized to the mitotic spindle but missing from the basal body in primary and photoreceptor cilia. This impairs recruitment of TZ components to the basal body and corresponds to complete loss of CEP162 function at the ciliary compartment, reflected by delayed formation of dysmorphic cilia. In contrast, rescue of increased cell death in the developing mouse retina after shRNA knockdown of Cep162 by expression of CEP162-E646R*5 indicates that the mutant retains its role for retinal neurogenesis. Human retinal degeneration thus results from specific loss of ciliary CEP162 function.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=134 HEIGHT=200 SRC="FIGDIR/small/469779v2_ufig1.gif" ALT="Figure 1">
View larger version (50K):
org.highwire.dtl.DTLVardef@83457eorg.highwire.dtl.DTLVardef@47775corg.highwire.dtl.DTLVardef@1ed643eorg.highwire.dtl.DTLVardef@5032d2_HPS_FORMAT_FIGEXP  M_FIG C_FIG Brominated Lipid Probes Expose Structural Asymmetries in Constricted Membranes Lipids in biological membranes are thought to be functionally organized, but few experimental tools can probe nanoscale membrane structure. Using brominated lipids as contrast probes for cryo-EM and a model ESCRT-III membrane remodeling system, we observed leaflet-level and protein-localized lipid structural patterns within highly constricted and thinned membrane nanotubes. These nanotubes differed markedly from protein-free, flat bilayers in leaflet thickness, lipid diffusion rates, and lipid compositional and conformational asymmetries. Simulations and cryo-EM imaging of brominated SDPC showed how a pair of phenylalanine residues scored the outer leaflet with a helical hydrophobic defect where polyunsaturated docosahexaenoyl (DHA) tails accumulated at the bilayer surface. Combining cryo-EM of halogenated lipids with molecular dynamics thus enables new characterizations of the composition and structure of membranes on molecular length scales.

One-Sentence SummaryCryo-EM imaging of halogenated lipids and MD simulations provide molecular-scale insight into constricted bilayer structure. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. A thermodynamic model for water activity and redox potential in evolution and development Reactions involving water and oxygen are basic features of geological and biological processes. To understand how life interacts with its environment requires monitoring interactions with H2O and O2 at timescales relevant to not only organismal growth but also billions of years of geobiological evolution. Here, chemical transformations intrinsic to evolution and development were characterized by analyzing data from recent phylostratigraphic and proteomic studies. This two-stage analysis involves obtaining chemical metrics (carbon oxidation state and stoichiometric hydration state) from the elemental compositions of proteins followed by modeling the relative stabilities of target proteins against a proteomic background to infer thermodynamic parameters (oxygen fugacity, water activity, and virtual redox potential (Eh)). The main results of this study are a rise in carbon oxidation state of proteins spanning the time of the Great Oxidation Event, a rise in virtual redox potential that coincides with the likely emergence of aerobic metabolism, a rise in carbon oxidation state of proteins inferred from the transcriptome in late stages of Bacillus subtilis biofilm growth, and a drop in stoichiometric hydration state of the fruit fly developmental proteome at the same time as a drop in organismal water content. Stoichiometric hydration state also decreases for proteins with more recent gene ages and through stages of biofilm development, leading to predictions of higher hydration potentials at earlier time points. By building chemical representations of protein evolution and developmental proteomes, exciting new types of geobiochemical proxies can be developed. Network Inference from Perturbation Time Course Data Networks underlie much of biology from subcellular to ecological scales. Yet, understanding what experimental data are needed and how to use them for unambiguously identifying the structure of even small networks remains a broad challenge. Here, we integrate a dynamic least squares framework into established modular response analysis (DL-MRA), that specifies sufficient experimental perturbation time course data to robustly infer arbitrary two and three node networks. DL-MRA considers important network properties that current methods often struggle to capture: (i) edge sign and directionality; (ii) cycles with feedback or feedforward loops including self-regulation; (iii) dynamic network behavior; (iv) edges external to the network; and (v) robust performance with experimental noise. We evaluate the performance of and the extent to which the approach applies to cell state transition networks, intracellular signaling networks, and gene regulatory networks. Although signaling networks are often an application of network reconstruction methods, the results suggest that only under quite restricted conditions can they be robustly inferred. For gene regulatory networks, the results suggest that incomplete knockdown is often more informative than full knockout perturbation, which may change experimental strategies for gene regulatory network reconstruction. Overall, the results give a rational basis to experimental data requirements for network reconstruction and can be applied to any such problem where perturbation time course experiments are possible. High resolution, serial imaging of early mouse and human liver bud morphogenesis in three dimensions BackgroundLiver organogenesis has thus far served as a paradigm for solid organ formation. The developing liver bud is a well established model of organogenesis, and murine genetic studies demonstrate key molecules involved in key morphogenetic changes. However, the analysis of the liver bud is typically limited to 2D tissue sections, which precludes extensive visualization, quantitation, and analysis. Further, the lack of human liver bud data has further hindered our understanding of human liver organogenesis. Therefore, new analytical and visualization approaches are needed to elicit further morphogenetic details of liver organogenesis and to elucidate differences between mouse and human liver bud growth.

ResultsTo address this need, we focused on high resolution imaging, visualization, and analysis of early liver growth by using available online databases for both mouse (EMAP, Toronto Phenogenomics center) and human (3D Atlas of Human Embryology), noninvasive multimodality imaging studies of the murine embryo, and mouse/human liver weight data. First, we performed three-dimensional (3D reconstructions) of stacked, digital tissue sections that had been initially segmented for the liver epithelium and the septum transversum mesenchyme (STM). 3D reconstruction of both mouse and human data sets enabled visualization and analysis of the dynamics of liver bud morphogenesis, including hepatic cord formation and remodeling, mechanisms of growth, and liver-epithelial STM interactions. These studies demonstrated potentially under-appreciated mechanisms of growth, including rapid exponential growth that is matched at the earliest stages by STM growth, and unique differences between mouse and human liver bud growth. To gain further insight into the exponential liver bud growth that was observed, we plotted volumetric data from 3D reconstruction together with fetal liver growth data from multimodality (optical projection tomography, magnetic resonance imaging, micro-CT) and liver weight data to compose complete growth curves during mouse (E8.5-E18) and human (day 25-300) liver development. For further analysis, we performed curve fitting and parameter estimation, using Gompertzian models, which enables the comparison between mouse and human liver bud growth, as well as comparisons to processes like liver regeneration. To demonstrate the importance of mesenchyme in rapid liver bud growth and morphogenesis in the human liver bud, we performed functional analysis in which human pluripotent stem cell (hPSC)-derived hepatic organoids were used to model collective migration that occurs in vivo, demonstrating that migration is strongly dependent upon mesenchyme.

DiscussionOur data demonstrates improved visualization with 3D images, under-appreciated and potentially new mechanisms of growth, complete liver growth curves with quantitative analysis through embryonic and fetal stages, and a new functional human stem cell-derived liver organoid assay demonstrating mesenchyme-driven collective migration. These data enhance our understanding of liver organogenesis. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. EPIC: inferring relevant cell types for complex traits by integrating genome-wide association studies and single-cell RNA sequencing More than a decade of genome-wide association studies (GWASs) have identified genetic risk variants that are significantly associated with complex traits. Emerging evidence suggests that the function of trait-associated variants likely acts in a tissue- or cell-type-specific fashion. Yet, it remains challenging to prioritize trait-relevant tissues or cell types to elucidate disease etiology. Here, we present EPIC (cEll tyPe enrIChment), a statistical framework that relates large-scale GWAS summary statistics to cell-type-specific gene expression measurements from single-cell RNA sequencing (scRNA-seq). We derive powerful gene-level test statistics for common and rare variants, separately and jointly, and adopt generalized least squares to prioritize trait-relevant cell types while accounting for the correlation structures both within and between genes. Using enrichment of loci associated with four lipid traits in the liver and enrichment of loci associated with three neurological disorders in the brain as ground truths, we show that EPIC outperforms existing methods. We apply our framework to multiple scRNA-seq datasets from different platforms and identify cell types underlying type 2 diabetes and schizophrenia. The enrichment is replicated using independent GWAS and scRNA-seq datasets and further validated using PubMed search and existing bulk case-control testing results. Ghost lineages highly influence the interpretation of introgression tests Most species are extinct; those that are not are often unknown. Sequenced and sampled species are often a minority of known ones. Past evolutionary events involving horizontal gene flow, such as horizontal gene transfer, hybridization, introgression and admixture, are therefore likely to involve "ghosts", i.e. extinct, unknown or unsampled lineages. The existence of these ghost lineages is widely acknowledged, but their possible impact on the detection of gene flow and on the identification of the species involved is largely overlooked. It is generally considered as a possible source of error that, with reasonable approximation, can be ignored. We explore the possible influence of absent species on an evolutionary study by quantifying the effect of ghost lineages on introgression as detected by the popular D-statistic method. We show from simulated data that under certain frequently encountered conditions, the donors and recipients of horizontal gene flow can be wrongly identified if ghost lineages are not taken into account. In particular, having a distant outgroup, which is usually recommended, leads to an increase in the error probability and to false interpretations in most cases. We conclude that introgression from ghost lineages should be systematically considered as an alternative possible, even probable, scenario. 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. 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.