Pathways of savannization in a mesic African savanna-forest mosaic following an extreme fire O_LIFires in savannas limit tree cover, thereby promoting flammable grass accumulation and fuelling further frequent fires. Meanwhile, forests and thickets form dense canopies that reduce C4-grass fuel loads and creating a humid microclimate, thereby excluding fires under typical climatic conditions.
C_LIO_LIHowever, extreme fires occasionally burn into these closed-canopy systems. Although these rare fires cause substantial tree mortality and can make repeat fires more likely, the long-term consequences of an extreme fire for closed canopy vegetation structure and potential to convert to savanna (hereafter "savannization") remain largely unknown.
C_LIO_LIHere, we analysed whether an extreme fire could, alone, alter species composition, vegetation structure, and fire regimes of closed-canopy ecosystems in an intact savanna-forest-thicket mosaic, or whether successive fires after an initial extreme fire were necessary to trigger a biome transition between from forest to savanna.
C_LIO_LIWe found that forests that only burned once recovered, whereas those that burned again following an initial extreme fire transitioned from closed-canopy forests towards open, grassy savannas.
C_LIO_LIWhile thickets had less tree mortality in fires than forests, repeat fires nonetheless precipitated a transition towards savannas.
C_LIO_LIColonization of the savanna tree community lagged behind the grass community, but also began to transition.
C_LI

SynthesisOur results suggest that rare extreme fires, followed by repeated burning can indeed result in savannization in places where savanna and forest represent alternative stable states.
From adverse to beneficial - contrasting dietary effects of freshwater mixotrophs on zooplankton O_LIThe importance of mixotrophic algae as key bacterivores in microbial food webs is increasingly acknowledged, but their effects on the next trophic level remain poorly understood. Their high stoichiometric food quality is contrasted by anti-grazing strategies.
C_LIO_LIWe tested the quality of freshwater mixotrophs as prey for zooplankton, using four non-colonial chrysophyte species and a cryptophyte as a high quality reference food. We (1) analyzed the stoichiometric and biochemical (fatty acid) composition of the mixotrophs, and (2) quantified their dietary effects on Daphnia longispina survival.
C_LIO_LISurvival of D. longispina significantly depended on the identity of species provided as food, ranging from higher to lower as compared to starvation. This was not reflected in differences in cellular stoichiometry or fatty acid profiles of the mixotrophs. We suggest that toxicity may be the driver for the observed differences.
C_LIO_LIGeneralization of the dietary effects of mixotrophic chrysophytes does not appear straightforward. Besides fundamental species-specific differences, potential toxic effects may vary depending on environmental cues or physiological strategies. Notably in our study, Ochromonas tuberculata, a species previously reported to be deleterious turned out to be a beneficial food source, in terms of enabling high survival of D. longispina.
C_LIO_LIWe challenge the generality of the assumption that chrysophytes are of low value as food for zooplankton. We recommend that future studies test how environmental conditions and physiological strategies shape the quality of mixotrophs as food for consumers at higher trophic levels, specifically focusing on effects of dietary toxicity.
C_LI
A phenomics approach for in vitro antiviral drug discovery BackgroundThe current COVID-19 pandemic has highlighted the need for new and fast methods to identify novel or repurposed therapeutic drugs. Here we present a method for untargeted phenotypic drug screening of virus-infected cells, combining Cell Painting with antibody-based detection of viral infection in a single assay. We designed an image analysis pipeline for segmentation and classification of virus-infected and non-infected cells, followed by extraction of morphological properties.

ResultsWe show that the methodology can successfully capture virus-induced phenotypic signatures of MRC-5 human lung fibroblasts infected with Human coronavirus 229E (CoV-229E). Moreover, we demonstrate that our method can be used in phenotypic drug screening using a panel of nine host- and virus-targeting antivirals. Treatment with effective antiviral compounds reversed the morphological profile of the host cells towards a non-infected state.

ConclusionsThe method can be used in drug discovery for morphological profiling of novel antiviral compounds on both infected and non-infected cells. Alkaline phosphatase activity of serum affects osteogenic differentiation cultures Fetal bovine serum (FBS) is a widely used supplement in cell culture medium, despite its known variability in composition which greatly affects cellular function and consequently the outcome of studies. In bone tissue engineering, the deposited mineralized matrix is one of the main outcome parameters, but using different brands of FBS can result in large variations. Alkaline phosphatase (ALP) is present in FBS. Not only is ALP used to judge the osteogenic differentiation of bone cells, it may affect deposition of mineralized matrix. The present study focused on the enzymatic activity of ALP in FBS of different suppliers and its contribution to mineralization in osteogenic differentiation cultures. It was hypothesized that culturing cells in a medium with high intrinsic ALP activity of FBS will lead to higher mineral deposition compared to media with lower ALP activity. The used FBS types were shown to have significant differences in enzymatic ALP activity. Our results indicate that the ALP activity of the medium not only affected the deposited mineralized matrix but also the osteogenic differentiation of cells as measured by a changed cellular ALP activity of human bone marrow derived mesenchymal stromal cells (hBMSC). In media with low inherent ALP activity, the cellular ALP activity was increased and played the major role in the mineralization process; while, in media with high intrinsic ALP activity contribution from the serum, less cellular ALP activity was measured and the ALP activity of the medium also contributed to mineral formation substantially. Our results highlight the diverse effects of ALP activity intrinsic to FBS on osteogenic differentiation and matrix mineralization and how FBS can determine the experimental outcomes, in particular for studies investigating matrix mineralization. Once again, the need to replace FBS with more controlled and known additives is highlighted.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=55 SRC="FIGDIR/small/462693v2_ufig1.gif" ALT="Figure 1">
View larger version (17K):
org.highwire.dtl.DTLVardef@902078org.highwire.dtl.DTLVardef@43b6fborg.highwire.dtl.DTLVardef@1a244fborg.highwire.dtl.DTLVardef@5ef819_HPS_FORMAT_FIGEXP  M_FIG C_FIG Do males know? Evidence-driven rainmaking mating strategy by male bonobos to meet the fertile window of females Female bonobos exhibit prolonged receptivity, one suggested function of which is to lower male mating competition. However, it is questionable whether easier access to receptive females can reduce male-male competition, given the exclusive nature of male reproductive success. We tested whether males could determine the fertile phase of females. We found that ovulation probability predicted male mating effort. High-ranking males copulated with females with higher fertility, and male-male agonistic interactions increased when there were fertile females in the party. When there were multiple females with maximal swelling, males concentrated their mating effort on females with an older infant whose maximal swelling started earlier, and they continued mating efforts until detumescence (rainmaking). These findings suggest that male bonobos distinguish between fertile and non-fertile phases of females and that having more receptive females in the party does not reduce male-male competition for fertile females.

TeaserMales use the rainmakers rule to meet the periovulatory phase of the female bonobo for better reproductive success. Adversarial domain translation networks for fast and accurate integration of large-scale atlas-level single-cell datasets The rapid emergence of large-scale atlas-level single-cell RNA-seq datasets presents remarkable opportunities for broad and deep biological investigations through integrative analyses. However, harmonizing such datasets requires integration approaches to be not only computationally scalable, but also capable of preserving a wide range of fine-grained cell populations. We created Portal, a unified framework of adversarial domain translation to learn harmonized representations of datasets. With innovation in model and algorithm designs, Portal achieves superior performance in preserving biological variation during integration, while achieving integration of millions of cells in minutes with low memory consumption. We show that Portal is widely applicable to integrating datasets across samples, platforms and data types (including scRNA-seq, snRNA-seq and scATAC-seq). Finally, we demonstrate the power of Portal by applying it to the integration of cross-species datasets with limited shared information among them, elucidating biological insights into the similarities and divergences in the spermatogenesis process among mouse, macaque and human. ACE2-independent SARS-CoV-2 infection and mouse adaption emerge after passage in cells expressing human and mouse ACE2 Human ACE2 (hACE2) is the key cell attachment and entry receptor for SARS-CoV-2, with the original SARS-CoV-2 isolates unable to use mouse ACE2 (mACE2). Herein we describe a new system for generating mouse-adapted SARS-CoV-2 in vitro by serial passaging virus in co-cultures of cell lines expressing hACE2 and mACE2. Mouse-adapted viruses emerged with up to five amino acid changes in the spike protein, all of which have been seen in human isolates. Mouse-adapted viruses replicated to high titers in C57BL/6J mouse lungs and nasal turbinates, and caused severe lung histopathology. One mouse-adapted virus was also able to replicate efficiently in ACE2-negative cell lines, with ACE2-independent entry by SARS-CoV-2 representing a new biology for SARS-CoV-2 that has potential widespread implications for disease and intervention development. A burst of regulatory and protein innovation at the origin of placental mammals drove the emergence of placenta and underpins divergent early pregnancy strategies in modern mammals. Summary paragraphThe origin of live birth in mammals [~]148 million years ago was a dramatic shift in reproductive strategy, yet the molecular changes that established mammal viviparity are largely unknown. Although progesterone receptor signalling predates the origin of mammals andis highly conserved in, and critical for, successful mammal pregnancy, it alone cannot explain the origin and subsequent diversity of implantation strategies throughout the placental mammal radiation. MiRNAs are known to be flexible and dynamic regulators with a well established role in the pathophysiology of mammal placenta. We propose that a dynamic core microRNA (miRNA) network originated early in placental mammal evolution, responds to conserved mammal pregnancy cues (e.g. progesterone), and facilitates species-specific responses. Here we identify 13 miRNAs that arose at the origin of placental mammals and were subsequently retained in all descendent lineages. The expression of these 13 miRNAs in response to early pregnancy molecules is regulated in a species-specific manner in endometrial epithelia of species with extremes of implantation strategies. Furthermore, these 13 miRNAs preferentially target 84 proteins under positive selective pressure on ancestral eutherian. Discovery of this core "live-birth" toolkit and specifically adapted proteins helps explain the origin and evolution of the placenta in mammals. An exon-intron split framework to prioritize miRNA-driven post-transcriptional regulatory signals and its application to study energy homeostasis in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Reconstitution of minimal motility system based on Spiroplasma swimming by expressing two bacterial actins in synthetic minimal bacterium Motility is one of the most important features of life, but its evolutionary origin remains unknown. In this study, we focused on Spiroplasma, commensal, or parasitic bacteria. They swim by switching the helicity of a ribbon-like cytoskeleton that comprises six proteins, each of which evolved from a nucleosidase and bacterial actin called MreB. We expressed these proteins in a synthetic, non-motile minimal bacterium, JCVI-syn3.0B, whose reduced genome was computer-designed and chemically synthesized. The synthetic bacterium exhibited swimming motility with features characteristic of Spiroplasma swimming. Moreover, some combinations of the two proteins produced a helical cell shape and swimming, suggesting that the swimming originated from the differentiation and coupling of bacterial actins, and we also obtained a minimal system for motility of the synthetic bacterium.

One-Sentence SummaryThe minimal system comprised two bacterial actins that provided cell helicity and swimming to the synthetic minimal bacterium. Membrane insertion mechanism of the caveolae coat protein Cavin1 Caveolae are small plasma membrane invaginations, important for control of membrane tension, signaling cascades and lipid sorting. The caveolae coat protein Cavin1 is essential for shaping such high curvature membrane structures. Yet, a mechanistic understanding of how Cavin1 assembles at the membrane interface is lacking. Here, we used model membranes combined with biophysical dissection and computational modelling to show that Cavin1 inserts into membranes. We establish that initial PI(4,5)P2-dependent membrane adsorption of the trimeric helical region 1 (HR1) of Cavin1 mediates the subsequent partial separation and membrane insertion of the individual helices. Insertion kinetics of the HR1 is further enhanced by the presence of flanking negatively charged disordered regions, which was found important for the co-assembly of Cavin1 with Caveolin1 in living cells. We propose that this intricate mechanism potentiates membrane curvature generation and facilitates dynamic rounds of assembly and disassembly of Cavin1 at the membrane.

Significance statementCaveolae are cholesterol enriched membrane invaginations coupled to severe muscle and lipid disorders. Their formation is dependent on assembly of the protein Cavin1 at the lipid membrane interface driving membrane curvature. In this work, we dissect the mechanism for how Cavin1 binds and inserts into membranes using a combination of biochemical and biophysical characterization as well as computational modelling. The proposed model for membrane assembly potentiates dynamic switching between shielded and exposed hydrophobic helices used for membrane insertion and clarifies how Cavin1 can drive membrane curvature and the formation of caveolae. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Liver Endothelium Promotes HER3-mediated Cell Survival in KRAS Wild-type and Mutant Colorectal Cancer Cells We previously showed that liver endothelial cells (ECs) secreted soluble factors in a paracrine fashion and activated human epidermal growth factor receptor 3 (HER3, also known as ERBB3) for promoting colorectal cancer (CRC) growth and chemoresistance. However, RAS proteins play a critical role in receptor tyrosine kinase signaling pathways, and KRAS mutations mediate CRC resistance to therapies targeting EGFR, another HER protein. Therefore, the role of KRAS mutation status in EC-induced HER3 activation and CRC survival was investigated as it has therapeutic implications. We used CRC cell lines and patient-derived xenografts harboring KRAS wild-type or mutant genes and demonstrated that liver EC-secreted factors promoted HER3-mediated CRC cell growth independent of KRAS mutation status. Also, blocking HER3 in CRC cells by siRNAs or a HER3 antibody seribantumab blocked EC-induced CRC survival. Our findings highlight the potential of utilizing HER3 targeted therapies for treating patients with mCRC independent of RAS mutational status. Essential Regression - a generalizable framework for inferring causal latent factors from multi-omic human datasets High-dimensional cellular and molecular profiling of human samples highlights the need for analytical approaches that can integrate multi-omic datasets to generate predictive biomarkers and prioritized causal inferences. Current methods are limited by high dimensionality of the combined datasets, the differences in their data distributions and their integration to infer causal relationships. Here we present Essential Regression (ER), an interpretable machine learning approach for high-dimensional multi-omic datasets, that addresses these problems by identifying latent factors and their likely cause-effect relationships with the system-wide outcome/properties of interest. ER is a novel data-distribution-free latent-factor regression model that integrates multi-omic datasets and identifies latent factors significantly associated with an outcome. ER outperforms a range of state-of-the-art methods in terms of prediction performance on simulated datasets. ER can be coupled with probabilistic graphical modeling thereby strengthening the causal inferences. ER generates novel cellular and molecular predictions, using multi-omic human systems immunology datasets, pertaining to immunosenescence and immune dysregulation. The Inflammatory Role of Milk Fat Globule Epidermal Growth Factor VIII in Angiotensin II Induced Arterial Remodeling BackgroundAngiotensin II (Ang II) and milk fat globule-epidermal growth factor VIII (MFG-E8) are involved in age-associated arterial remodeling; however, the inflammatory role of MFG-E8 in Ang II associated arterial remodeling with aging remains unknown.

Methods and ResultsIn this study, 30-week-old MFG-E8 knock out (KO) and age-matched wild-type (WT) mice were infused with Ang II or saline. After infusion the with Ang II, the aortic molecular, cellular, and structural remodeling were observed in mice and compared to those infused with saline, but these effects were dependent on the expression of MFG-E8: (1) In the WT mice, Ang II infusion substantially increased intimal-medial thickness, elastic lamina degradation, collagen deposition, and the proliferation of VSMCs; in contrast, in the KO mice, these effects were significantly reduced; (2) In the WT mice, Ang II treatment significantly increased the activation and expression of MMP2, TGF-{beta}1, and its downstream signaling molecule p-SMAD2, and collagen type I production, however, in the KO mice, these molecular effects were significantly reduced; (3) In the WT mice, Ang II treatment increased inflammatory p-NF-{kappa}B p65, MCP1, TNF-, ICAM1, and VCAM1 molecular expression, while conversely, in the KO mice, no significant inflammatory changes were found; (4) Importantly, compared to untreated control mice with a wide range of age from 4-96 weeks, Ang II infused "younger" mice produced an "older" arterial inflammatory phenotype, which was alleviated by MFG-E8 deficiency.

ConclusionsMFG-E8 mediates Ang II associated arterial inflammatory remodeling. Targeting MFG-E8 is a novel molecular approach to curb adverse arterial remodeling during aging and hypertension.

CLINICAL PERSPECTIVEO_ST_ABSWhat Is New?C_ST_ABS* Both Ang II and MFG-E8 increases are involved in proinflammatory arterial remodeling mediating the molecular, cellular and tissue events in aging and hypertension.
* MFG-E8 is essential for Ang II induced and age-associated adverse arterial remodeling via the increase of proinflammation, intimal medial thickening, elastin fragmentation, collagen deposition, and VSMC proliferation.


What Are the Clinical Implications?Since MFG-E8 mediates Ang II induced proinflammation in arterial wall remodeling in aging and hypertension, targeting MFG-E8 is a potential molecular approach to curb inflammatory arterial remodeling, maintaining the health of the vascular system during aging and hypertension. Secretory cells are the primary source of pIgR in small airways BackgroundLoss of secretory immunoglobulin A (SIgA) is common in COPD small airways and likely contributes to disease progression. We hypothesized loss of SIgA results from reduced expression of pIgR, a chaperone protein needed for SIgA transcytosis, in the COPD small airway epithelium.

MethodspIgR-expressing cells were defined and quantified at single-cell resolution in human airways using RNA in-situ hybridization, immunostaining, and single-cell RNA sequencing. Complementary studies in mice utilized immunostaining, primary murine tracheal epithelial cell (MTEC) culture, and transgenic mice with secretory or ciliated cell-specific knockout of pIgR. SIgA degradation by human neutrophil elastase or secreted bacterial proteases from non-typeable Haemophilus influenzae (NTHi) was evaluated in vitro.

ResultsWe found that secretory cells are the predominant cell type responsible for pIgR expression in human and murine airways. Loss of SIgA in small airways was not associated with a reduction in secretory cells but rather a reduction in pIgR protein expression despite intact PIGR mRNA expression. Neutrophil elastase and NTHi-secreted proteases are both capable of degrading SIgA in vitro and may also contribute to a deficient SIgA immunobarrier in COPD.

InterpretationLoss of the SIgA immunobarrier in small airways of patients severe COPD is complex and likely results from both pIgR-dependent defects in IgA transcytosis and SIgA degradation.

Key MessagesWhat is the key question? Localized SIgA deficiency in small airways is an established driver of COPD pathogenesis, but the mechanism of loss remains unclear. We hypothesized loss of SIgA is due to reduced numbers of pIgR-expressing cells in SIgA-deficient small airways.

What is the bottom line? pIgR is primarily expressed by secretory cells in human and murine airways. Although numbers of secretory cells are similar between SIgA-deficient and SIgA-replete airways in COPD, there is reduced expression of pIgR protein, but not mRNA, in SIgA-deficient airways. Additionally, host and bacterial proteases degrade SIgA in vitro, suggesting loss of SIgA may relate to both impaired transcytosis and increased degradation.

Why read on? This study highlights the complexity of SIgA immunobarrier maintenance and suggests that strategies aimed at restoring the SIgA immunobarrier will need to account for both impaired transcytosis and degradation by host and/or bacterial proteases. PIN-FORMED1 polarity in the shoot is insensitive to the polarity of neighbouring cells In plants, epidermal cells are planar-polarized along an axis marked by the asymmetric localization patterns of several proteins including PIN-FORMED1 (PIN1), which facilitates the directional efflux of the plant hormone auxin to pattern phyllotaxis (Heisler et al., 2010; Mansfield et al., 2018; Reinhardt et al., 2003). While PIN1 polarity is known to be regulated non-cell autonomously via the MONOPTEROS (MP) transcription factor, how this occurs has not been determined (Bhatia et al., 2016). Here we use mosaic expression of the serine threonine kinase PINOID (PID) to test whether PIN1 polarizes according to the polarity of neighbouring cells. Our findings reveal that PIN1 is insensitive to the polarity of PIN1 in neighbouring cells arguing against auxin flux or extracellular auxin concentrations acting as a polarity cue, in contrast to previous model proposals (Abley et al., 2016; Stoma et al., 2008). Classification of Pacemaker Dynamics in the Mouse Intestine by Field Potential Microimaging The flexible and sophisticated movement of the gastrointestinal (GI) tract implies the involvement of mechanisms other than enteric neural circuits, to coordinate excitation in microregions. We thus performed microimaging of pacemaker dynamics in the small intestine of mice since it contains typical network-forming pacemaker cells. A dialysis membrane-reinforced low-impedance microelectrode array (MEA) enabled field potentials over a wide frequency range to be stably measured in microregions. The pacemaker dynamics were classified into basic patterns despite large variations. In the developmental process, pacemaker activity was categorized as either an  expanding or a  migrating pattern that was initiated in or propagated to the MEA sensing area, respectively. The intercellular current of the volume conductor complicated the waveform of both activities. The existence of  expanding and  migrating patterns was attributable to duplicated pacemaker systems such as intracellular Ca2+ oscillation-activated and voltage-gated mechanisms. Additionally, from the spatio-temporal feature during the period of pacemaker events, the  bumpy/aberrant pattern was defined by aberrant, incoherent propagation, and associated with local impairment of excitability, while the  colliding/converging pattern involved the interaction of multiple activities in the MEA area. Interconversion between the four micro-coordination patterns occurred in the same microregion. 5-Hydroxytryptamine (5-HT) promoted  migrating activity, implying an improvement or restoration of spatial conductivity. These results agree well with the action of 5-HT to change GI movement toward propulsion. In conclusion, our MEA method of microimaging classification enables the quantitative assessment of spatio-temporal electric coordination underlying GI motility, suggesting its application to small model animals. Characterization of regulatory transcriptional mechanisms in hepatocyte lipotoxicity Non-alcoholic fatty liver disease is a continuum of disorders among which non-alcoholic steatohepatitis (NASH) is particularly associated with a negative prognosis. Hepatocyte lipotoxicity is one of the main pathogenic factors of liver fibrosis and NASH. However, the molecular mechanisms regulating this process are poorly understood. The main aim of this study was to dissect transcriptional mechanisms regulated by lipotoxicity in hepatocytes. We achieved this aim by combining transcriptomic, proteomic and chromatin accessibility analyses from human liver and mouse hepatocytes. This integrative approach revealed several transcription factor networks deregulated by NASH and lipotoxicity. To validate these predictions, genetic deletion of the transcription factors MAFK and TCF4 was performed, resulting in hepatocytes that were better protected against saturated fatty acid oversupply. MAFK- and TCF4-regulated gene expression profiles suggest a mitigating effect against cell stress, while promoting cell survival and growth. Moreover, in the context of lipotoxicity, some MAFK and TCF4 target genes were to the corresponding differentially regulated transcripts in human liver fibrosis. Collectively, our findings comprehensively profile the transcriptional response to lipotoxicity in hepatocytes, revealing new molecular insights and providing a valuable resource for future endeavours to tackle the molecular mechanisms of NASH. 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. 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. Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin BackgroundEvolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation at many enhancers, distal regulatory elements that control spatio-temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific, code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

ResultsWe present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that is able to predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Then, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity and liver enhancers associated with liver regeneration have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

ConclusionThe framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements. HYBRID REACTION-DIFFUSION AND CLOCK-AND-WAVEFRONT MODEL FOR THE ARREST OF OSCILLATIONS IN THE SOMITOGENESIS SEGMENTATION CLOCK AO_SCPLOWBSTRACTC_SCPLOWThe clock and wavefront paradigm is arguably the most widely accepted model for explaining the embryonic process of somitogenesis. According to this model, somitogenesis is based upon the interaction between a genetic oscillator, known as segmentation clock, and a differentiation wavefront, which provides the positional information indicating where each pair of somites is formed. Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation. Cotterell et al. called it a progressive oscillatory reaction-diffusion (PORD) model. In the Meinhardt-PORD model, somitogenesis is driven by short-range interactions and the posterior movement of the front is a local, emergent phenomenon, which is not controlled by global positional information. With this model, it is possible to explain some experimental observations that are incompatible with the clock and wavefront model. However the Meinhardt-PORD model has some important disadvantages of its own. Namely, it is quite sensitive to fluctuations and depends on very specific initial conditions (which are not biologically realistic). In this work, we propose an equivalent Meinhardt-PORD model, and then amend it to couple it with a wavefront consisting of a receding morphogen gradient. By doing so, we get a hybrid model between the Meinhardt-PORD and the clock-and-wavefront ones, which overcomes most of the deficiencies of the two originating models.

Somitogenesis, the process by which somites are formed, is an essential developmental stage in many vertebrates. This process occurs with a strikingly regular periodicity, that is preserved among embryos of a single species. The clock and wavefront paradigm is arguably the most widely accepted model for explaining somitogenesis. However, it is incapable of explaining some experimental facts, like the appearance of somites in the absence of an external wavefront (i.e. a receding morphogen gradient). Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation, and called it a progressive oscillatory reaction-diffusion (PORD) model. The Meinhardt-PORD model tackles some of the deficiencies of the clock and wave-front models, but it has some serious issues of its own. In the present work, we introduce an equivalent Meinhardt-PORD model, and then amend it to couple it with a receding morphogen gradient. By doing so, we get a hybrid model that incorporates characteristics of the Meinhardt-PORD and clock-and-wavefront models. We show that this hybrid model under-goes a bifurcation, from a stable to an unstable limit cycle, as the value of the parameter accounting for a background regulatory input (associated to the receding morphogen gradient) decreases. This bifurcation allows the model to explain why somites can form in the absence of an external wavefront, reassesses the role of the receding morphogen gradient as a conductor for somitogenesis, and makes the model behavior robust to random fluctuations, as well as independent from specific initial conditions (the latter, are two of the weak points of the Meinhardt-PORD model). We argue that this findings provide convincing evidence that reaction-diffusion and positional information (receding morphogen gradient) mechanisms could work together in somitogenesis. FcγR responses to soluble immune complexes are governed by solubility and size Fc-gamma receptor (Fc{gamma}R) activation by soluble IgG immune complexes (sICs) represents a major mechanism of inflammation in certain autoimmune diseases such as systemic lupus erythematosus (SLE). A robust and scalable test system allowing for the detection and quantification of sIC bioactivity is missing. Previously described Fc{gamma}R interaction assays are limited to certain Fc{gamma}Rs, lack scalability and flexibility, are not indicative of receptor activation or lack sensitivity towards sIC size. We developed a comprehensive reporter cell panel detecting individual activation of Fc{gamma}Rs from humans and the mouse. The reporter cell lines were integrated into an assay format that provides flexible read-outs enabling the quantification of sIC reactivity via ELISA or a fast detection using flow cytometry. This identified Fc{gamma}RIIA(H) and Fc{gamma}RIIIA as the most sIC-sensitive Fc{gamma}Rs in our test system. Applying the assay we demonstrate that sICs versus immobilized ICs are fundamentally different Fc{gamma}R-ligands with regard to Fc{gamma}R preference and signal strength. Reaching a detection limit in the very low nanomolar range, the assay proved also to be sensitive to sIC stoichiometry and size enabling for the first time a complete reproduction of the Heidelberger-Kendall precipitation curve in terms of immune receptor activation. Analyzing sera from SLE patients and mouse models of lupus and arthritis proved that sIC-dependent Fc{gamma}R activation has predictive capabilities regarding severity of SLE disease. The new methodology provides a sensitive, scalable and comprehensive tool to evaluate the size, amount and bioactivity of sICs in all settings.

One Sentence SummaryIn this study we established a comprehensive Fc{gamma}R reporter cell assay enabling the detection and quantification of soluble immune complexes generated in experimental and clinical settings. P2X2 receptor subunit interfaces are missense variant hotspots where mutations tend to increase apparent ATP affinity Background and PurposeP2X receptors (P2XRs) are trimeric ligand-gated ion channels (LGICs) that open a cation-selective pore in response to ATP binding to their large extracellular domain (ECD). The seven known P2XR subtypes can assemble as homo- or heterotrimeric complexes and contribute to numerous physiological functions, including nociception, inflammation and hearing. The overall structure of P2XRs is well established, but little is known about the spectrum and prevalence of human genetic variations and the functional implications in specific domains.

Experimental ApproachHere we examine the impact of P2X2 receptor (P2X2R) inter-subunit interface missense variants identified in the human population or through structural predictions. We test both single and double mutants through electrophysiological and biochemical approaches.

Key resultsWe demonstrate that predicted ECD inter-subunit interfaces display a higher-than-expected density of missense variations and that the majority of mutations that disrupt putative inter-subunit interactions result in channels with higher apparent ATP affinity. Lastly, we show that double mutants at the subunit interface show significant energetic coupling, especially if located in close proximity.

Conclusions and ImplicationsWe provide the first structural mapping of the mutational burden across the human population in a LGIC and show that the density of missense mutations is constrained between protein domains, indicating evolutionary selection at the domain level. Our data may indicate that, unlike other LGICs, P2X2Rs have evolved an intrinsically high threshold for activation, possibly to allow for additional modulation or as a cellular protection mechanism against overstimulation.

Bullet point summary What is already known: O_LIP2X2 receptors are ATP-activated ion channels implicated in hearing and nociceptice pathways
C_LI
 What this study adds: O_LIA structural mapping of missense variants observed in the human population
C_LIO_LIWe identify the intersubunit-interface as a variant hotspot and decipher functional impact of mutations
C_LI
 Clinical significance: O_LIThe development of both inhibitors and activators of P2X2 receptor function may be required
C_LI Cross-Frequency Multilayer Network Analysis with Bispectrum-based Functional Connectivity: A Study of Alzheimer's Disease Alzheimers disease (AD) is a neurodegenerative disorder known to affect functional connectivity (FC) across many brain regions. Linear FC measures have been applied to study the differences in AD by splitting neurophysiological signals, such as electroencephalography (EEG) recordings, into discrete frequency bands and analysing them in isolation from each other. We address this limitation by quantifying cross-frequency FC in addition to the traditional within-band approach. Cross-bispectrum, a higher-order spectral analysis approach, is used to measure the nonlinear FC and is compared with the cross-spectrum, which only measures the linear FC within bands. This work reports the first use of cross-bispectrum to reconstruct a cross-frequency FC network where each frequency band is treated as a layer in a multilayer network with both inter- and intra-layer edges. An increase of within-band FC in AD cases is observed in low-frequency bands using both methods. Cross-bispectrum also detects multiple cross-frequency differences, mainly increased FC in AD cases in{delta} -{theta} coupling. Overall, increased strength of low-frequency coupling and decreased level of high-frequency coupling is observed in AD cases in comparison to healthy controls (HC). Vulnerability analysis reveals that the integration and segregation properties of networks are enabled by different frequency couplings in AD networks compared to HCs. Moreover, the segregation property depends less on the{gamma} interactions in AD cases, suggesting a shift from high to low-frequency connectivity. The results highlight the importance of studying nonlinearity and including cross-frequency FC in characterising AD. Moreover, the results demonstrate the advantages and limitations of using cross-bispectrum to reconstruct FC networks.