Genomic inference of a human super bottleneck in Mid-Pleistocene transition The demographic history is a foundation of human evolutionary studies. However, the ancient demographic history during the Mid-Pleistocene is poorly investigated while it is essential for understanding the early origin of humankind. Here we present the fast infinitesimal time coalescent (FitCoal) process, which allows the analytical calculation of the composite likelihood of a site frequency spectrum and provides the precise inference of demographic history. We apply it to analyze 3,154 present-day human genomic sequences. We find that African populations have passed through a population super bottleneck, a small effective size of approximately 1,280 breeding individuals between 930 and 813 thousand years ago. Further analyses confirm the existence of the super bottleneck on non-African populations although it cannot be directly inferred. This observation, together with simulation results, indicates that confounding factors, such as population structure and selection, are unlikely to affect the inference of the super bottleneck. The time interval of the super bottleneck coincides with a gap in the human fossil record in Africa and possibly marks the origin of Homo heidelbergensis. Our results provide new insights into human evolution during the Mid-Pleistocene.
Extensive transgressive gene expression in testis but not ovary in the homoploid hybrid Italian sparrow Hybridization can result in novel allelic combinations which can impact the hybrid phenotype through changes in gene expression. While mis-expression in F1 hybrids is well documented, how gene expression evolves in stabilized hybrid taxa remains an open question. As gene expression evolves in a stabilizing manner, break-up of co-evolved cis- and trans-regulatory elements could lead to transgressive patterns of gene expression in hybrids. Here, we address to what extent gonad gene expression has evolved in an established and stable homoploid hybrid, the Italian sparrow (Passer italiae). Through comparison of gene expression in gonads from individuals of the two parental species (i.e., house and Spanish sparrow) to that of Italian sparrows, we find evidence for strongly transgressive expression in male Italian sparrows - 22% of testis genes exhibit expression patterns outside the range of both parent species. In contrast, Italian sparrow ovary expression was similar to that of one of the parent species, the house sparrow (Passer domesticus). Moreover, the Italian sparrow testis transcriptome is 26 times as diverged from those of the parent species as the parental transcriptomes are from each other, despite being genetically intermediate. This highlights the potential for regulation of gene expression to produce novel variation following hybridization. Genes involved in mitochondrial respiratory chain complexes and protein synthesis are enriched in the subset that is over-dominantly expressed in Italian sparrow testis, suggesting that selection on key functions have molded the hybrid Italian sparrow transcriptome.
Dynamics of SARS-CoV-2 mutations reveals regional-specificity and similar trends of N501 and High-Frequency mutation N501Y in different levels of control measures Coronavirus disease 2019 (COVID-19) is a contagious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This disease has spread globally, causing more than 161.5 million cases and 3.3 million deaths to date. Surveillance and monitoring of new mutations in the virus genome are crucial to our understanding of the adaptation of SARS-CoV-2. Moreover, how the temporal dynamics of these mutations is influenced by control measures and non-pharmaceutical interventions (NPIs) is poorly understood. Using 1 058 020 SARS-CoV-2 from sequenced COVID-19 cases from 98 countries (totaling 714 country-month combinations), we perform a normalization by COVID-19 cases to calculate the relative frequency of SARS-CoV-2 mutations and explore their dynamics over time. We found 115 mutations estimated to be present in more than 3 % of global COVID-19 cases and determined three types of mutation dynamics: High-Frequency, Medium-Frequency, and Low-Frequency. Classification of mutations based on temporal dynamics enable us to examine viral adaptation and evaluate the effects of implemented control measures in virus evolution during the pandemic. We showed that Medium-Frequency mutations are characterized by high prevalence in specific regions and/or in constant competition with other mutations in several regions. Finally, taking N501Y mutation as representative of High-Frequency mutations, we showed that level of control measure stringency negatively correlates with the effective reproduction number of SARS-CoV-2 with High-Frequency or not-High-Frequency and both follows similar trends in different levels of stringency. Region-specific microRNA alterations in marmosets carrying SLC6A4 polymorphisms are associated to anxiety-like behavior. Psychiatric diseases such as depression and anxiety are multifactorial conditions, highly prevalent in western societies. Human studies have identified a number of high-risk genetic variants for these diseases. Among them, polymorphisms in the promoter region of the serotonin transporter gene (SLC6A4) have attracted much attention. However, due to the paucity of experimental models, molecular alterations induced by these genetic variants and how they correlate to behavioral deficits have not been examined. Marmosets have emerged as a powerful model in translational neuroscience to investigate molecular underpinnings of complex behaviors. Here, we took advantage of naturally occurring genetic polymorphisms in marmoset SLC6A4 gene that have been linked to anxiety-like behaviors. Using FACS-sorted cells from different brain regions, we revealed that marmosets bearing different SLC6A4 variants exhibit distinct microRNAs signatures in a region of the prefrontal cortex whose activity has been consistently altered in patients with depression/anxiety. We also identified DCC, a gene previously linked to these diseases, as a downstream target of the dysregulated microRNAs. Significantly, we showed that levels of both microRNAs and DCC in this region were highly correlated to anxiety-like behaviors as well as to the response to citalopram, a selective serotonin re-uptake inhibitor and widely prescribed anti-depressant. Our findings establish links between genetic variants, molecular modifications in specific cortical regions and complex behavioral/pharmacological responses, providing new insights into gene-behavior relationships underlying human psychopathology. Epitopedia: identifying molecular mimicry between pathogens and known immune epitopes Upon infection, foreign antigenic proteins stimulate the hosts immune system to produce antibodies targeting the pathogen. These antibodies bind to regions on the antibody called epitopes. Structural similarity (molecular mimicry) of epitopes between an infecting pathogen and host proteins or other pathogenic proteins the host has previously encountered can impact the host immune response to the pathogen and may lead to cross-reactive antibodies. The ability to identify potential regions of molecular mimicry in a pathogen can illuminate immune effects which are especially important to pathogen treatment and vaccine design. Here we present Epitopedia, a software pipeline that facilitates the identification of regions that may exhibit potential three-dimensional molecular mimicry between an antigenic pathogen protein and known immune epitopes as catalogued by the immune epitope database (IEDB). Epitopedia is open-source software released under the MIT license and is freely available on GitHub, including a Docker container with all other software dependencies preinstalled. We performed an analysis describing how various secondary structure states, identity between pentapeptide pairs, and identity between the parent sequences of pentapeptide pairs affects RMSD. We found that pentapeptides pairs in a helical conformation had considerably lower RMSD values than those in Extended or Coil conformations. We also found that RMSD is significantly increased when pentapeptide pairs are from non-homologous sequences. Brain endothelial STING1 activation by Plasmodium-sequestered heme promotes cerebral malaria via type I IFN response. Cerebral malaria (CM) is a life-threatening form of Plasmodium falciparum infection caused by brain inflammation. Brain endothelium dysfunction is a hallmark of CM pathology, which is also associated with the activation of the type I interferon (IFN) inflammatory pathway. The molecular triggers and sensors eliciting brain type I IFN cellular responses during CM remain largely unknown. We herein identified the stimulator of interferon response cGAMP interactor 1 (STING1) as the key innate immune sensor that induces Ifn{beta}1 transcription in the brain of mice infected with Plasmodium berghei (Pba). This STING1/IFN{beta}-mediated response increases brain CXCL10 governing the extent of brain leucocyte infiltration and blood-brain barrier (BBB) breakdown, and determining CM lethality. The critical role of brain endothelial cells (BECs) in fueling type I IFN-driven brain inflammation was demonstrated in brain endothelial-specific IFN{beta}-reporter and STING1-deficient Pba-infected mice, which are significantly protected from CM lethality. Moreover, extracellular particles (EPs) released from Pba-infected erythrocytes activated STING1-dependent type I IFN response in BECs, a response requiring intracellular acidification. Fractionation of the EPs enabled us to identify a defined fraction carrying hemoglobin degradation remnants that activates STING1/IFN{beta} in the brain endothelium, a process correlated with heme content. Notably, stimulation of STING1-deficient BECs with heme, docking experiments and in vitro binding assays unveiled that heme is a putative STING1 ligand. This work shows that heme resultant from the parasite heterotrophic activity operates as an alarmin triggering brain endothelial inflammatory responses via STING1/IFN{beta}/CXCL10 axis crucial to CM pathogenesis and lethality.

SignificanceCM results from loss of blood-brain endothelial barrier function caused by unrestrained inflammatory response in the natural course of infection by Plasmodium parasites. However, the role of brain endothelium in triggering inflammatory mechanisms is still undetermined. We found that the innate immune sensor STING1 is crucial for production of IFN{beta} in brain endothelial cells in Plasmodium-infected mice. This in turn stimulates CXCL10-mediated recruitment of leukocytes and subsequent brain inflammation and tissue damage. We identified within extracellular particles released from Plasmodium-infected erythrocytes, a fraction containing products of hemoglobin degradation, namely heme, which we show can bind STING1. Our results unravel a new angle of CM pathogenesis: heme contained in particles triggers the STING/IFN{beta}/CXCL10 axis in brain endothelial cells. Electrical Stimulated GLUT4 Signaling Attenuates Critical Illness-Associated Muscle Wasting BackgroundCritical illness myopathy (CIM) is a debilitating condition characterized by the preferential loss of the motor protein myosin. CIM is a byproduct of critical care, attributed to impaired recovery, long-term complications, and mortality. CIM pathophysiology is complex, heterogeneous and remains incompletely understood, however loss of mechanical stimuli contributes to critical illness associated muscle atrophy and weakness. Passive mechanical loading (ML) and electrical stimulation (ES) therapies augment muscle mass and function. While having beneficial outcomes, the mechanistic underpinning of these therapies is less known. Therefore, here we aimed to assess the mechanism by which chronic supramaximal ES ameliorates CIM in a unique experimental rat model of critical care.

MethodsRats were subjected to 8 days critical care conditions entailing deep sedation, controlled mechanical ventilation, and immobilization with and without direct soleus ES. Muscle size and function were assessed at the single cell level. RNAseq and Western blotting were employed to understand the mechanisms driving ES muscle outcomes in CIM.

ResultsFollowing 8 days of controlled mechanical ventilation and immobilization, soleus muscle mass, Myosin:Actin ratio and single muscle fiber maximum force normalized to cross-sectional area (specific force) were reduced by 40-50% (p< 0.0001). ES significantly reduced the loss of soleus muscle fiber cross-sectional area (CSA) and Myosin:Actin ratio by approximately 30% (p< 0.05) yet failed to effect specific force. RNAseq pathway analysis revealed downregulation of insulin signaling in the soleus muscle following critical care and GLUT4 trafficking was reduced by 55% leading to an 85% reduction of muscle glycogen content (p< 0.01). ES promoted phosphofructokinase and insulin signaling pathways to control levels (p< 0.05), consistent with the maintenance of GLUT4 translocation and glycogen levels. AMPK, but not AKT, signaling pathway was stimulated following ES, where the downstream target TBC1D4 increased 3 logFC (p= 0.029) and AMPK-specific P-TBC1D4 levels were increased approximately 2-fold (p= 0.06). Reduction of muscle protein degradation rather than protein synthesis promoted soleus CSA, as ES reduced E3 ubiquitin proteins, Atrogin-1 (p= 0.006) and MuRF1 (p= 0.08) by approximately 50%, downstream of AMPK-FoxO3.

ConclusionsES maintained GLUT4 translocation through increased AMPK-TBC1D4 signaling leading to improved muscle glucose homeostasis. Soleus CSA and myosin content was promoted through reduced protein degradation via AMPK-FoxO3 E3 ligases, Atrogin-1 and MuRF1. These results demonstrate chronic supramaximal ES reduces critical care associated muscle wasting, preserved glucose signaling and reduced muscle protein degradation in CIM. Evaluating the Anti-aging Effects of Chinese Herbal Medicine Rhodiola rosea L. in Cultured Keratinocytes Among the many health benefits Chinese herbal medicine (CHM) presents, anti-aging is of special interest (Shen, Jiang, Yang, Wang, & Zhu, 2017; Zhao & Luo, 2017). Reported to possess anti-aging effects, the CHM Rhodiola rosea L., known colloquially as the "golden root", has been widely incorporated in various drinks, daily supplements, and even cosmetics. This study investigates the effects of commercial R. rosea extracts and natural R. rosea roots on preventing UV-induced photoaging of the skin and correlates such effects with the composition of active ingredients in the extracts. To simulate the photoaging process, drug treated HaCaT cells were exposed to UVA and UVB radiation. The pharmacological anti-aging effects of R. rosea extracts were evaluated qualitatively and quantitatively through confocal immunofluorescence images with {gamma}-H2A.X marker and telomerase activity assay (Telo TAGGG Telomerase PCR-ELISA). Preparatory thin layer chromatography and high-performance liquid chromatography were performed to isolate and quantify active ingredients. Cultured HaCaT cells showed morphological change after exposure to both UVA (>15.0 J) and UVB (>2000 mJ). The photoaging of keratinocytes was rescued by pretreating cells with R. rosea extracts as well as salidroside and rosavin active ingredients (P < 0.05). R. rosea-treated cells were characterized by increased telomerase activity and fewer {gamma}-H2A.X foci compared to that of the control. Extracts with better preventative effects contained higher salidroside and rosavin content. The findings in this study reaffirm R. roseas efficacy as an anti-aging remedy and provide a basis for CHMs integration into the mainstream of global healthcare. 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. Geometric machine learning informed by ground truth: Recovery of conformational continuum from single-particle cryo-EM data of biomolecules This work is based on the manifold-embedding approach to study biological molecules exhibiting continuous conformational changes. Previous work established a method capable of reconstructing 3D movies and accompanying energetics of atomic-level structures from single-particle cryo-EM images of macromolecules displaying multiple conformational degrees of freedom. Here, we introduce an unsupervised geometric machine learning approach that is informed by detailed heuristic analysis of manifolds formed by simulated heterogeneous cryo-EM datasets generated from an atomic structure. These simulated data were generated with increasing complexity to account for multiple conformational motions, state occupancies and typical microscope parameters in a wide range of signal-to-noise ratios. Using these datasets as ground-truth, we provide detailed exposition of our findings using several conformational motions while exploring the available parameter space. Guided by these insights, we build a framework to leverage the high-dimensional geometric information obtained towards reconstituting a quasi-continuum of conformational states in the form of a free-energy landscape and respective 3D density maps for all states therein. As shown by a direct comparison of results, this framework offers substantial improvements relative to the previous work. Synergy of chemotherapy and macrophage depletion leads to B cell-mediated T cell memory activation and durable triple negative breast cancer regression Immunosuppressive elements within the tumor microenvironment such as Tumor Associated Macrophages (TAMs) can present a barrier to successful anti-tumor responses by cytolytic T cells. We employed preclinical syngeneic p53 null mouse models of triple negative breast cancer (TNBC) to develop a treatment regimen that harnessed the immunostimulatory effects of low-dose cyclophosphamide coupled with the pharmacologic inhibition of TAMs using either a small molecule CSF1R inhibitor or an anti-CSF1R antibody. This therapeutic combination was used to successfully treat several highly aggressive TNBC murine mammary tumors and lung metastasis. Using this regimen and single cell RNA sequencing we characterized tumor infiltrating lymphocytes (TILs) including helper T cells and antigen-presenting B cells that were highly enriched in good responders to combination therapy. Using high dimensional imaging techniques, we identified the close spatial localization of B220+ CD86+ activated B cells and CD4+ T cells in tertiary lymphoid structures that were present up to 6 weeks post-treatment in one model that also exhibited long-term tumor regression post-treatment. We also characterized the transcriptional and metabolic heterogeneity of TAMS in these two closely related claudin-low/mesenchymal subtype tumor models with differential treatment responses. A murine TAM signature derived from the T12 model is highly expressed and conserved in human claudin-low breast cancers, and high expression of the T12 signature correlated with reduced overall survival. This T12 tumor TAM signature may help identify human claudin-low breast cancer patients that will benefit from the combination of cyclophosphamide and anti-CSF1R therapy. These studies illustrate the complexity of the tumor immune microenvironment and highlight different immune responses that result from rationale combinations of immunotherapy.

SignificanceA treatment regimen that harnessed the immunostimulatory effects of cyclophosphamide coupled with the inhibition of CSF1R was used to successfully treat several highly aggressive claudin-low TNBC murine mammary tumors and lung metastasis. A novel method for producing functionalized vesicles that efficiently deliver oligonucleotides in vitro in cancer cells and in vivo in mice Nano-based delivery systems have enhanced our ability to administer and target drugs and macromolecules to their targets. Oligonucleotide drugs have great therapeutic potential but often have off-target effects and stability issues. Therefore, they are often encapsulated in vesicles with targeting ligands such as antibodies (Ab) to deliver their cargo. Herein, we describe a novel, scalable and straightforward approach to producing functionalized vesicles called the "Functionalized Lipid Insertion Method," which differs from the older "Detergent-Dialysis Method." The latter method required excess detergent and extensive dialysis over many hours to produce the functionalized vesicles. With our method, only the functionalized lipid is detergent-solubilized during the engineering of the vesicle. The approach reduces the dialysis time, keeps the vesicles intact while orienting the targeting moieties of the functionalized lipid toward the outside of the vesicle. Pilot in vitro and in vivo experiments was performed to show the feasibility of our method. Dynamic light scattering (DLS) experiments suggested that the original vesicular structure was relatively unperturbed, and the functionalized lipid was inserted externally. Our approach efficiently delivered oligonucleotides and affected the function of liver cancer HepG2 cells. Furthermore, functionalized vesicles achieved targeted delivery of oligonucleotides in mice without inducing a significant innate immune response. The industrial and therapeutic significance and implications of functionalized vesicles produced by our method are also discussed. Additional experiments and analyses are recommended to bring out the full potential of this molecular delivery technology. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. Protective host-dependent antagonism among 	Pseudomonas in the Arabidopsis phyllosphere The plant microbiome is a rich biotic environment, comprising numerous taxa. The community structure of these colonizers is constrained by multiple factors, including host-microbe and microbe-microbe interactions, as well as the interplay between the two. While much can be learned from pairwise relationships between individual hosts and microbes, or individual microbes with themselves, the ensemble of interrelations between the host and microbial consortia may lead to different outcomes that are not easily predicted from the individual interactions. Their study can thus provide new insights into the complex relationship between plants and microbes. Of particular importance is how strain-specific such plant-microbe-microbe interactions are, and how they eventually affect plant health. Here, we test strain-level interactions in the phyllosphere between groups of co-existing commensal and pathogenic Pseudomonas among each other and with A. thaliana, by employing synthetic communities of genome-barcoded isolates. We found that commensal Pseudomonas prompted a host response leading to a selective inhibition of a specific pathogenic lineage, resulting in plant protection. The extent of plant protection, however, was dependent on plant genotype, indicating that these effects are host-mediated. There were similar genotype-specific effects on the microbe side, as we could pinpoint an individual Pseudomonas isolate as the predominant cause for this differential interaction. Collectively, our work highlights how within-species genetic differences on both the host and microbe side can have profound effects on host-microbe-microbe dynamics. The paradigm that we have established provides a platform for the study of host-dependent microbe-microbe competition and cooperation in the A. thaliana-Pseudomonas system. 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. Diet effects on mouse meiotic recombination: a warning for recombination studies Meiotic recombination is a critical process for sexually reproducing organisms. This exchange of genetic information between homologous chromosomes during meiosis is important not only because it generates genetic diversity, but also because it is often required for proper chromosome segregation. Consequently, the frequency and distribution of crossovers are tightly controlled to ensure fertility and offspring viability. However, in many systems it has been shown that environmental factors can alter the frequency of crossover events. Two studies in flies and yeast point to nutritional status affecting the frequency of crossing over. However, this question remains unexplored in mammals. Here we test how crossover frequency varies in response to diet in Mus musculus males. We use immunohistochemistry to estimate crossover frequency in multiple genotypes under two diet treatments. Our results indicate that while crossover frequency was unaffected by diet in some strains, other strains were sensitive even to small composition changes between two common laboratory chows. Therefore, recombination is both resistant and sensitive to certain dietary changes in a strain-dependent manner and, hence, this response is genetically determined. Our study is the first to report a nutrition effect on genome-wide levels of recombination. Moreover, our work highlights the importance of controlling diet in recombination studies and may point to diet as a potential source of variability among studies, which is relevant for reproducibility. 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. Towards a mathematical understanding of colonization resistance in multispecies microbial communities Microbial community composition and dynamics are key to health and disease. Explaining the forces generating and shaping diversity in the microbial consortia making up our bodys defenses is a major aim of current research in microbiology. For this, tractable models are needed, that bridge the gap between observations of patterns and underlying mechanisms. While most microbial dynamics models are based on the Lotka-Volterra framework, we still do not have an analytic quantity for colonization resistance, by which a microbial systems fitness as a whole can be understood. Here, inspired by an epidemiological perspective, we propose a rather general modeling framework whereby colonization resistance can be clearly mathematically defined and studied. In our model, N similar species interact with each other through a co-colonization interaction network encompassing pairwise competition and cooperation, abstractly mirroring how organisms effectively modify their micro-scale environment in relation to others. This formulation relies on a generic notion of shared resources between members of a consortium, yielding explicit frequency-dependent dynamics among N species, in the form of a replicator equation, and offering a precise definition of colonization resistance. We demonstrate that colonization resistance arises and evolves naturally in a multispecies system as a collective quadratic term in a replicator equation, describing dynamic mean invasion fitness. Each pairwise invasion growth rate between two ecological partners, [Formula], is derived explicitly from species asymmetries and mean traits. This makes the systemic colonization resistance [Formula] also an emergent function of global mean-field parameters and trait variation architecture, weighted by the evolving relative abundances among species. In particular, if the underlying invasion fitness matrix {Lambda} displays species-specific  invasiveness or  invasibility structure, colonization resistance will be insensitive to mean micro-scale cooperation or competition. However, in general, colonization resistance depends on and may undergo critical transitions with changes in mean  environment, e.g. cooperation and growth level in a community. We illustrate several key links between our proposed measure of colonization resistance and invader success, including sensitivity to timing, and to the intrinsic pairwise invasion architecture of the resident community. Our simulations reveal that symmetric and invader-driven mutual invasion among resident species tend to maximize systemic colonization resistance to outsiders, when compared to resident-driven, anti-symmetric, almost anti-symmetric and random {Lambda} structures. We contend this modeling approach is a powerful new avenue to study, test and validate interaction networks and invasion topologies in diverse microbial consortia, and quantify analytically their role in colonization resistance, system function, and invasibility. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Coupled myovascular expansion directs cardiac growth and regeneration BackgroundHeart regeneration requires multiple cell types to enable cardiomyocyte (CM) proliferation. How these cells interact to create growth niches is unclear. Here we profile proliferation kinetics of cardiac endothelial cells (CECs) and CMs in the neonatal mouse heart and find that they are spatiotemporally coupled.

Methods and ResultsWe show that coupled myovascular expansion during cardiac growth or regeneration is dependent upon VEGF-VEGFR2 signaling, as genetic deletion of Vegfr2 from CECs or inhibition of VEGFA abrogates both CEC and CM proliferation. Repair of cryoinjury displays poor spatial coupling of CEC and CM proliferation. Boosting CEC density after cryoinjury with virus encoding Vegfa enhances regeneration. Using Mendelian randomization, we demonstrate that circulating VEGFA levels are positively linked with human myocardial mass, suggesting that Vegfa can stimulate human cardiac growth.

ConclusionsOur work demonstrates the importance of coupled CEC and CM expansion and reveals a myovascular niche that may be therapeutically targeted for heart regeneration. Col6a1+/CD201+ mesenchymal cells regulate intestinal morphogenesis and homeostasis Intestinal mesenchymal cells encompass multiple subsets, whose origins, functions, and pathophysiological importance are still not clear. Here, we used the Col6a1Cre mouse, which targets telocytes and perivascular cells that can be further distinguished by the combination of the CD201, PDGFR and SMA markers. Developmental studies revealed that the Col6a1Cre mouse also targets mesenchymal aggregates that are crucial for intestinal morphogenesis and patterning, suggesting an ontogenic relationship between them and homeostatic telocytes. Cell depletion experiments in adulthood showed that Col6a1+/CD201+ mesenchymal cells regulate homeostatic enteroendocrine cell differentiation and epithelial proliferation. During acute colitis, they expressed an inflammatory and extracellular matrix remodeling gene signature, but they also retained their properties and topology. Notably, both in homeostasis and tissue regeneration, they were dispensable for normal organ architecture, while CD34+ mesenchymal cells expanded, localised at the top of the crypts and showed increased expression of villous-associated morphogenetic factors, providing thus evidence for the plasticity potential of distinct mesenchymal populations in the intestine. Our results provide a comprehensive analysis of the identities, origin, and functional significance of distinct mesenchymal populations in the intestine. N-VEGF, the autoregulatory arm of VEGF-A Vascular endothelial growth factor A (VEGF-A) is a secreted protein that stimulates angiogenesis in response to hypoxia. Under hypoxic conditions, a non-canonical long isoform called L-VEGF is concomitantly expressed with VEGF-A. Once translated, L-VEGF it is proteolytically cleaved to generate N-VEGF and VEGF-A. Interestingly, while VEGF-A is secreted and affects the surrounding cells, N-VEGF is mobilized to the nucleus. This suggests that N-VEGF participates in transcriptional response to hypoxia. In this study, we performed a series of complementary experiments to examine the functional role of N-VEGF. Strikingly, we found that the mere expression of N-VEGF followed by its hypoxia-independent mobilization to the nucleus was sufficient to induce key genes associated with angiogenesis, such as Hif1, VEGF-A isoforms, as well as genes associated with cell survival under hypoxia. Complementarily, when N-VEGF was genetically depleted, key hypoxia-induced genes were downregulated and cells were significantly susceptible to hypoxia-mediated apoptosis. This is the first reports of N-VEGF serving as an autoregulatory arm of VEGF-A. Further experiments will be needed to determine the role of N-VEGF in cancer and embryogenesis. Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable and highly reproducible in terms of microbial composition, function and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation. Variation in female songbird state determines signal strength needed to evoke copulation It is the female response to male signals that determines courtship success. In most songbirds, females control reproduction via the copulation solicitation display (CSD), an innate, stereotyped posture produced in direct response to male displays. Because CSD can be elicited in the absence of males by the presentation of recorded song, CSD production enables investigations into the effects of underlying signal features and behavioral state on female mating preferences. Using computer vision to quantify CSD trajectory in female brown-headed cowbirds (Molothrus ater), we show that both song quality and a females internal state predict CSD production, as well as the onset latency and duration of the display. We also show that CSD can be produced in a graded fashion based on both signal strength and internal state. These results emphasize the importance of underlying receiver state in determining behavioral responses and suggest that female responsiveness acts in conjunction with male signal strength to determine the efficacy of male courtship. A biochemically-realisable relational model of the self-manufacturing cell As shown by Hofmeyr, the processes in the living cell can be divided into three classes of efficient causes that produce each other, so making the cell closed to efficient causation, the hallmark of an organism. They are the enzyme catalysts of covalent metabolic chemistry, the intracellular milieu that drives the supramolecular processes of chaperone-assisted folding and self-assembly of polypeptides and nucleic acids into functional catalysts and transporters, and the membrane transporters that maintain the intracellular milieu, in particular its electrolyte composition. Each class of efficient cause can be modelled as a relational diagram in the form of a mapping in graph-theoretic form, and a minimal model of a self-manufacturing system that is closed to efficient causation can be constructed from these three mappings using the formalism of relational biology. This Fabrication-Assembly or (F,A)-system serves as an alternative to Robert Rosens replicative Metabolism-Repair or (M,R)-system, which has been notoriously problematic to realise in terms of real biochemical processes. A key feature of the model is the explicit incorporation of formal cause, which arrests the infinite regress that plagues all relational models of the cell. The (F,A)-system is extended into a detailed relational model of the self-manufacturing cell that has a clear biochemical realisation. This (F,A) cell model allows the interpretation and visualisation of concepts such as the metabolism and repair components of Rosens (M,R)-system, John von Neumanns universal constructor, Howard Pattees symbol-function split via the symbol-folding transformation, Marcello Barbieris genotype-ribotype-phenotype ontology, and Tibor Gantis chemoton. Response of soil and vegetation in a warm-temperate Pine forest to intensive biomass harvests, phosphorus fertilisation, and wood ash application. Background and AimsConcerns about climate change and carbon economy have prompted the promotion of alternative energy sources, including forest-based bioenergy. An evaluation of the environmental consequences of intensive harvests (stumps and roots, and also branches and foliage) for energy wood supply, and use of wood-ash recycling as a compensatory practice, helps in the evaluation of the use of forest biomass for energy production.

MethodsWe made use of records from a split-plot experimental site crossing four different intensities of biomass harvesting (Stem-Only Harvest [SOH], Aboveground Additional Harvest [AAH], Belowground Additional Harvest [BAH], and Whole-Tree Harvest [WTH]) and three compensation methods (control [C], wood ash application [A] and phosphorus fertilisation [P]) to evaluate, in the 11-years-old stand (maritime pine: Pinus pinaster) that followed the biomass exports of the former stand, their effects on nutrient budgets, tree growth, soil fertility, chemical properties and soil carbon. This site is located in a forest on a poor soil, under a warm temperate climate (SW France).

Key resultsDespite their low additional biomass exports (+10% for AAH to +34% for WTH), the non-conventional harvest practices exported much higher quantities of nutrients than the conventional SOH technique (for example +145% for N and K in WTH). Consequently, these treatments had negative effects on the soil nutritive status. Additional biomass harvests impacted the soil organic matter content, with negative effects on P-organic, soil cation exchange capacity, exchangeable Ca, and most extractible nutrients. However, data suggested that tree growth and foliage nutrient content had not yet been significantly impacted by harvest treatments, whereas tree nutritional status was improved by P-fertiliser or wood ash. As expected, we observed a positive effect of wood ash application on soil pH and nutrient content but, like additional harvests, wood ash application decreased the pool of soil organic carbon (~10% of the initial stock with ~7% of N-total losses).

ConclusionsOverall, this factorial experiment showed that exporting more forest biomass due to the additional harvesting of tree canopies, stumps and roots had negative consequences on the ecosystem biogeochemistry. Additional harvests have aggravated the poverty of the already oligotrophic soil, and decreased the soil organic carbon content. Importantly, applying nutrients as fertiliser or wood ash did not compensate for all the negative impacts of biomass exports and the method of wood ash recycling in forests could even decrease the soil organic carbon.