A generic hidden Markov model for multi-parent populations A common step in the analysis of multi-parent populations is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of multi-parent populations. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly-available data on Diversity Outbred and Collaborative Cross mice
Birth Order Differences in Education Originate in Post-Natal Environments Siblings share many environments and much of their genetics. Yet, siblings turn out different. Intelligence and education are influenced by birth order, with earlier-born siblings outperforming later-borns. We investigate whether birth order differences in education are caused by biological differences present at birth, that is, genetic differences or in-utero differences. Using family data that spans two generations, combining registry, survey, and genotype information, this study is based on the Norwegian Mother, Father and Child Cohort Study (MoBa). We show that there are no genetic differences by birth order as captured by polygenic scores (PGSs) for educational attainment. Furthermore, we show that earlier-born have lower birth weight than later-born, indicating worse uterine environments. Educational outcomes are still higher for earlier-born children when we adjust for PGSs and in utero variables, indicating that birth order differences arise post-natally. Finally, we consider potential environmental influences, such as differences according to maternal age, parental educational attainment, and sibling genetic nurture. We show that birth order differences are not biological in origin, but pinning down their specific causes remains elusive.
HOW CONCERNING IS A SARS-COV-2 VARIANT OF CONCERN? COMPUTATIONAL PREDICTIONS AND THE VARIANTS LABELING SYSTEM We herein report a study to evaluate the use of computational prediction of SARS-CoV-2 genetic variations in improving the current variants labeling system. First, we reviewed the basis of the system developed by the World Health Organization (WHO) for the labeling of SARS-CoV-2 genetic variants and the adaptations made to it by the United States Center of Diseases Control (CDC). We observed that the labeling system is based upon the virus major attributes. However, we found that the labeling criteria of the SARS-CoV-2 variants derived from these attributes are not accurately defined and are used differently by the two health management agencies. Consequently, discrepancies exist between the labels given by WHO and CDC to same variants. Our observations suggest that giving the VOC label to a new variant is premature and might not be appropriate. Therefore, we carried out a comparative computational study to predict the effects of the mutations on the virus structure and functions of five VOCs. By linking these data to the criteria used by WHO and the CDC for variant labeling, we ascertained that comparative computational predictions of the impact of genetic variations are a better ground for rapid and more accurate labelling of SARS-CoV-2 variants. We propose to label all emergent variants VUM or VBM and to carry out computational predictive studies and thorough variants comparison, upon which more appropriate and informative labels can be attributed. Furthermore, harmonization of the variants labeling system would be globally beneficial to communicate about and fight COVID-19 pandemic. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. A unified mechanism for innate and learned visual landmark guidance in the insect central complex Insects can navigate efficiently in both novel and familiar environments, and this requires flexiblity in how they are guided by sensory cues. A prominent landmark, for example, can ellicit strong innate behaviours (attraction or menotaxis) but can also be used, after learning, as a specific directional cue as part of a navigation memory. However, the mechanisms that allow both pathways to co-exist, interact or override each other are largely unknown. Here we propose a model for the behavioural integration of innate and learned guidance based on the neuroanatomy of the central complex (CX), adapted to control landmark guided behaviours. We consider a reward signal provided either by an innate attraction to landmarks or a long-term visual memory in the mushroom bodies (MB) that modulates the formation of a local vector memory in the CX. Using an operant strategy for a simulated agent exploring a simple world containing a single visual cue, we show how the generated short-term memory can support both innate and learned steering behaviour. In addition, we show how this architecture is consistent with the observed effects of unilateral MB lesions in ants that cause a reversion to innate behaviour. We suggest the formation of a directional memory in the CX can be interpreted as transforming rewarding (positive or negative) sensory signals into a mapping of the environment that describes the geometrical attractiveness (or repulsion). We discuss how this scheme might represent an ideal way to combine multisensory information gathered during the exploration of an environment and support optimized cue integration.

1 Author summaryIn this paper, we modeled the neural pathway allowing insects to perform landmark guided behaviours using their internal compass. First, we observed the intrinsic property of the connectome, extracted from drosophila online database, between the internal compass neurons and the steering neurons to support an oriented behaviour towards a single landmark. Then, we proposed and evaluated an adaptation of the bees path integration neural circuit, to sustain flexible landmark guidance behaviours such as attraction or menotaxis. We showed the model ability to form a memory during the exploration of the local environment to support both innate or learned navigation behaviour using a single landmark in the environment. In addtion, we demonstrated the transformation of a simple goodness/badness signal, from innate or long-term memory pathways, into an oriented steering signal that could be applied to other sensory pathways. Furthermore, by reproducing lesion experiments in the mushroom bodies of wood ants we highlight the consistency of the model with biological observation. We then discuss the different emergent properties and the potential outcome that this local, and operant, memory supports. A novel Boolean network inference strategy to model early hematopoiesis aging Hematopoietic stem cell (HSC) aging is a multifactorial event that leads to changes in HSC properties and function. These changes are intrinsically coordinated and affect the early hematopoiesis, involving hematopoietic stem and progenitor cells (HSPCs). The objective of this work is to better understand the mechanisms and factors controlling these changes. We have therefore developed an original strategy to construct a Boolean network of genes explaining the priming and homeostasis of HSCs (graphical abstract). Based on our previous scRNA-seq data, we performed an exhaustive analysis of the transcriptional network and identified active transcription modules or regulons along the differentiation trajectory of selected HSPC states. This global view of transcriptional regulation led us to focus on 15 components, 13 selected TFs (Tal1, Fli1, Gata2, Gata1, Zfpm1, Egr1, Junb, Ikzf1, Myc, Cebpa, Bclaf1, Klf1, Spi1) and 2 complexes regulating the ability of HSC to cycle (CDK4/6 - Cyclin D and CIP/KIP). We then defined the connections controlling the differentiation dynamics of HSC states and constructed an influence graph between the TFs involved in the dynamics by mixing observations from our scRNA-seq data and knowledge from the literature. Then, using answer set programming (ASP) and in silico perturbation analysis, we obtained a Boolean model which is the solution of a Boolean satisfiability problem. Finally, perturbation of the model based on age-related changes revealed important regulations, such as the overactivation of Egr1 and Junb or the loss of Cebpa activation by Gata2, which were found to be relevant for the myeloid bias of aged HSC. Our work shows the efficiency of the combination of manual and systematic methods to elaborate a Boolean model. The developed strategy led to the proposal of new regulatory mechanisms underlying the differentiation bias of aged HSCs, explaining the decreased transcriptional priming of HSCs to all mature cell types except megakaryocytes. Co-infection of Candidatus Piscichlamydia trichopodus (order Chlamydiales) and Henneguya sp (Myxosporea, Myxobolidae) in snakeskin gourami Trichopodus pectoralis (Regan 1910) The present study describes a simultaneous infection of a novel Chlamydia-like organism (CLO) with a Myxozoa parasite, Henneguya sp. in snakeskin gourami Trichopodus pectoralis in Thailand. A new CLO is proposed "Candidatus Piscichlamydia trichopodus" (CPT) based on 16S rRNA phylogenetic analysis. Systemic intracellular CPT infection was confirmed by histological examination, in situ hybridization, PCR assay, and sequencing of 16S rRNA. This novel pathogen belongs to the order Chlamydiales but differs in certain aspects from other species. The histopathological changes associated with CPT infection were different from the typical pathological lesions of epitheliocystis caused by previously known CLO. Unlike other CLO, CPT localized in the connective tissue rather than in the epithelial cells and formed smaller clumps of intracellular bacteria that stained dark blue with hematoxylin. On the other hand, typical myxospores of the genus Henneguya with tails were observed in the gill sections. Infection with Henneguya sp. resulted in extensive destruction of the gill filaments, most likely leading to respiratory distress. Due to the frequency of co-infections and the unavailability of culture methods for CLO and Henneguya sp, it was difficult to determine which pathogens were directly responsible for the associated mortality. However, co-infections may increase the negative impact on the host and the severity of the disease. Given the commercial importance of the snakeskin gourami and its significant aquaculture potential, the findings of this study are important for further studies on disease prevention. Genomic regions associate with major axes of variation driven by gas exchange and leaf construction traits in cultivated sunflower (Helianthus annuus L.) Stomata and leaf veins play an essential role in transpiration and the movement of water throughout leaves. These traits are thus thought to play a key role in the adaptation of plants to drought and a better understanding of the genetic basis of their variation and coordination could inform efforts to improve drought tolerance. Here, we explore patterns of variation and covariation in leaf anatomical traits and analyze their genetic architecture via genome-wide association (GWA) analyses in cultivated sunflower (Helianthus annuus L.). Traits related to stomatal density and morphology as well as lower order veins were manually measured from digital images while the density of minor veins was estimated using a novel deep learning approach. Leaf, stomatal, and vein traits exhibited numerous significant correlations that generally followed expectations based on functional relationships. Correlated suites of traits could further be separated along three major principal component (PC) axes that were heavily influenced by variation in traits related to gas exchange, leaf hydraulics, and leaf construction. While there was limited evidence of colocalization when individual traits were subjected to GWA analyses, major multivariate PC axes that were most strongly influenced by several traits related to gas exchange or leaf construction did exhibit significant genomic associations. These results provide insight into the genetic basis of leaf trait covariation and showcase potential targets for future efforts aimed at modifying leaf anatomical traits in sunflower.

Significance StatementUsing traditional and automated/high-throughput (using a novel deep learning approach) phenotyping methods we studied leaf anatomical variation in sunflower. Genome-wide association (GWA) analyses identified numerous genomic regions underlying individual trait variation and regions underlying major multivariate axes of phenotypic variation. These results illustrate the value of employing a multivariate approach to GWA analyses and shed light on the extent to which leaf trait (co-)variation can be genetically decoupled to explore novel phenotypic space. Source-sink connectivity: A novel interictal EEG marker for seizure localization Over 15 million epilepsy patients worldwide have drug-resistant epilepsy (DRE). Successful surgery is a standard of care treatment for DRE but can only be achieved through complete resection or disconnection of the epileptogenic zone (EZ), the brain region(s) where seizures originate. Surgical success rates vary between 20-80% because no clinically validated biological markers of the EZ exist. Localizing the EZ is a costly and time-consuming process beginning with non-invasive neuroimaging and often followed by days to weeks of intracranial EEG (iEEG) monitoring. Clinicians visually inspect iEEG data to identify abnormal activity (e.g., low-voltage high frequency activity) on individual channels occurring immediately before seizures or spikes that occur on interictal iEEG (i.e., between seizures). In the end, the clinical standard mainly relies on a small proportion of the iEEG data captured to assist in EZ localization (minutes of seizure data versus days of recordings), missing opportunities to leverage these largely ignored interictal data to better diagnose and treat patients.

Intracranial EEG offers a unique opportunity to observe epileptic cortical network dynamics but waiting for seizures increases patient risks associated with invasive monitoring. In this study, we aim to leverage interictal iEEG data by developing a new network-based interictal iEEG marker of the EZ. We hypothesize that when a patient is not clinically seizing, it is because the EZ is inhibited by other regions. We developed an algorithm that identifies two groups of nodes from the interictal iEEG network: those that are continuously inhibiting a set of neighboring nodes ("sources") and the inhibited nodes themselves ("sinks"). Specifically, patient-specific dynamical network models (DNMs) were estimated from minutes of iEEG and their connectivity properties revealed top sources and sinks in the network, with each node being quantified by source-sink metrics (SSMs). We validated the SSMs in a retrospective analysis of 65 patients by using the SSMs of the annotated EZ to predict surgical outcomes. The SSMs predicted outcomes with an accuracy of 79% compared to an accuracy of 43% for clinicians predictions (surgical success rate of this dataset). In failed outcomes, we identified regions of the brain with high SSMs that were untreated. When compared to high frequency oscillations, the most commonly proposed interictal iEEG feature for EZ localization, SSMs outperformed in predictive power (by a factor of 1.2) suggesting SSMs may be an interictal iEEG fingerprint of the EZ. 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. Stitching and registering highly multiplexed whole slide images of tissues and tumors using ASHLAR MotivationStitching microscope images into a mosaic is an essential step in the analysis and visualization of large biological specimens, particularly human and animal tissues. Recent approaches to highly-multiplexed imaging generate high-plex data from sequential rounds of lower-plex imaging. These multiplexed imaging methods promise to yield precise molecular single-cell data and information on cellular neighborhoods and tissue architecture. However, attaining mosaic images with single-cell accuracy requires robust image stitching and image registration capabilities that are not met by existing methods.

ResultsWe describe the development and testing of ASHLAR, a Python tool for coordinated stitching and registration of 103 or more individual multiplexed images to generate accurate whole-slide mosaics. ASHLAR reads image formats from most commercial microscopes and slide scanners, and we show that it performs better than existing open source and commercial software. ASHLAR outputs standard OME-TIFF images that are ready for analysis by other open-source tools and recently developed image analysis pipelines.

Availability and implementationASHLAR is written in Python and available under an MIT license at https://github.com/labsyspharm/ashlar. An informational website with user guides and test data is available at https://labsyspharm.github.io/ashlar/. Negative DNA supercoiling makes protein-mediated looping deterministic and ergodic within the bacterial doubling time Protein-mediated DNA looping is fundamental to gene regulation and such loops occur stochastically in purified systems. Additional proteins increase the probability of looping, but these probabilities maintain a broad distribution. For example, the probability of lac repressor-mediated looping in individual molecules ranged 0-100%, and individual molecules exhibited representative behavior only in observations lasting an hour or more. Titrating with HU protein progressively compacted the DNA without narrowing the 0-100% distribution. Increased negative supercoiling produced an ensemble of molecules in which all individual molecules more closely resembled the average. Furthermore, in only twelve minutes of observation, well within the doubling time of the bacterium, most molecules exhibited the looping probability of the ensemble. DNA supercoiling, an inherent feature of all genomes, appears to impose time-constrained, emergent behavior on otherwise random molecular activity. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. The compound of baicalein, wogonin, and oroxylin-A inhibits EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. Macroscopic control of synchronous electrical signaling with chemically-excited gene expression Excitable cells can convert electrical signals into chemical outputs to facilitate the active transport of information across larger distances. This electrical-to-chemical conversion requires a tightly regulated expression of ion channels. Alterations of ion channel expression provide landmarks of numerous pathological diseases, such as cardiac arrhythmia, epilepsy, or cancer. Although the activity of ion channels can be locally regulated by external light or chemical stimulus, it remains challenging to coordinate the expression of ion channels on extended spatial-temporal scales in a non-invasive manner. Here, we have engineered yeast S. cerevisiae to read and convert local chemical concentrations into a dynamic electrical field distributed across cell populations. The core mechanism encodes a chemically-excitable dual-feedback gene circuit that precisely tunes the expression domain of potassium channels, globally coordinating cyclic firing of the plasma membrane potential (PMP). We demonstrate that this mechanism leverages an engineered constitutively open bacterial potassium channel KcsA to directly couple chemical stimuli with ion flux through gene expression and it can interface with the host ion channels through the pulsatile production of toxins. Our study provides a robust synthetic transcriptional toolbox underlying the conversion of local chemical environments into spatiotemporally organized electrical impulses for various cellular engineering, synthetic biology, and potential therapeutic applications. The human brain networks mediating the vestibular sensation of self-motion Vestibular Agnosia - where peripheral vestibular activation triggers the usual reflex nystagmus response but with attenuated or no self-motion perception - is found in brain disease with disrupted cortical network functioning, e.g. traumatic brain injury (TBI) or neurodegeneration (Parkinsons Disease). Patients with acute focal hemispheric lesions (e.g. stroke) do not manifest vestibular agnosia. Thus brain network mapping techniques, e.g. resting state functional MRI (rsfMRI), are needed to interrogate functional brain networks mediating vestibular agnosia. Whole-brain rsfMRI was acquired from 39 prospectively recruited acute TBI patients with preserved peripheral vestibular function, along with self-motion perceptual thresholds during passive yaw rotations in the dark. Following quality-control checks, 25 patient scans were analyzed. TBI patients were classified as having vestibular agnosia (n = 11) or not (n = 14) via laboratory testing of self-motion perception. Using independent component analysis, we found altered functional connectivity in the right superior longitudinal fasciculus and left rostral prefrontal cortex in vestibular agnosia. Moreover, regions of interest analyses showed both inter-hemispheric and intra-hemispheric network disruption in vestibular agnosia. In conclusion, our results show that vestibular agnosia is mediated by bilateral anterior and posterior network dysfunction and reveal the distributed brain mechanisms mediating vestibular self-motion perception. Loss of SUMO-specific protease 2 causes isolated glucocorticoid deficiency by blocking adrenal cortex zonal transdifferentiation SUMOylation is a dynamic posttranslational modification, that provides fine-tuning of protein function involved in the cellular response to stress, differentiation, and tissue development. In the adrenal cortex, an emblematic endocrine organ that mediates adaptation to physiological demands, the SUMOylation gradient is inversely correlated with the gradient of cellular differentiation raising important questions about its role in functional zonation and the response to stress. Considering that SUMO-specific protease 2 (SENP2), a deSUMOylating enzyme, is upregulated by ACTH/PKA signalling within the zona Fasciculata (zF), we generated mice with adrenal-specific Senp2 loss to address these questions. Disruption of SENP2 activity in steroidogenic cells leads to specific hypoplasia of the zF, a blunted reponses to ACTH and isolated glucocorticoid deficiency. Mechanistically, overSUMOylation resulting from SENP2 loss shifts the balance between ACTH/PKA and WNT/{beta}-catenin signalling leading to repression of PKA activity and ectopic activation of {beta}-catenin. At the cellular level, this blocks transdifferentiation of {beta}-catenin-positive zona Glomerulosa cells into zF cells and sensitises them to premature apoptosis. Our findings indicate that the SUMO pathway is critical for adrenal homeostasis and stress responsiveness. An endosome-associated actin network involved in directed apical plasma membrane growth The subcellular membrane-trafficking system plays many roles in cell growth and morphogenesis, from bulk membrane provision to targeted delivery of membrane proteins and other cargos. In tracheal terminal cells of the Drosophila respiratory system, endocytosis and trafficking through late endosomes balances membrane delivery between the basal plasma membrane and the apical membrane that forms a subcellular tube, but it has been unclear how the direction of growth of the subcellular tube with the overall growth of the cell is coordinated. We show here that endosomes are involved in this process as coordinators of F-actin. Actin assembles around late endocytic vesicles in the growth cone of the cell, reaching from the tip of the subcellular tube to the leading filopodia of the basal membrane. Preventing actin nucleation at late endosomes disturbs the directionality of tube growth, uncoupling it from the direction of cell elongation. Severing actin in this area affects tube integrity. Our findings show a new role for late endosomes in directing cell morphogenesis by organizing actin, in addition to their known role in membrane and protein trafficking. Transmembrane helices are an overlooked and evolutionarily conserved source of major histocompatibility complex class I and II epitopes Cytolytic T cell responses are predicted to be biased towards membrane proteins. The peptide-binding grooves of most haplotypes of histocompatibility complex class I (MHC-I) are relatively hydrophobic, therefore peptide fragments derived from human transmembrane helices (TMHs) are predicted to be presented more often as would be expected based on their abundance in the proteome. However, the physiological reason of why membrane proteins might be over-presented is unclear. In this study, we show that the over-presentation of TMH-derived peptides is general, as it is predicted for bacteria and viruses and for both MHCI and MHC-II. Moreover, we show that TMHs are evolutionarily more conserved, because single nucleotide polymorphisms (SNPs) are present relatively less frequently in TMH-coding chromosomal regions compared to regions coding for extracellular and cytoplasmic protein regions. Thus, our findings suggest that both cytolytic and helper T cells respond more to membrane proteins, because these are evolutionary more conserved. We speculate that TMHs therefore are less prone to escape mutations that enable pathogens to evade T cell responses. Stability of the genetic structure of a wild wheat population across microhabitats over 36 years Long-term genetic studies of wild populations are very scarce but are essential for connecting ecological and population genetics models, and are also necessary for understanding how climate change can affect biodiversity. We present a study of the wild ancestor of modern wheat, carried out over a 36-year period during which atmospheric CO2 increased by 70 ppm and temperature by 1.5{degrees}C. We genotyped 832 individuals that had been collected along the same transects during the course of the experiment. To our surprise, very little changed. Genotypes were clustered into ecological microhabitats over scales of tens of metres, and this clustering was stable over the 36 generations of the study. Simulations indicate that neutral processes alone are unlikely to fully explain this stability, suggesting a role for fine-scale local adaptation. Our results provide a rare insight into the dynamics of a wild population at the scale of individuals. The neurogenic fate of the hindbrain boundaries relies on Notch-dependent asymmetric cell divisions The generation of cell diversity in the central nervous system occurs during embryogenesis and requires a precise balance between cell proliferation, commitment to specific fates, and further neuronal differentiation. The cellular and molecular mechanisms regulating this balance in the embryonic brain are still poorly understood. Here we study how the neurogenic capacity in the embryonic hindbrain is spatiotemporally allocated, and how the fate and the growth of the hindbrain boundary cells are regulated. By generating a CRISPR-based knock-in zebrafish transgenic line to specifically label the hindbrain boundaries, we unveiled that boundary cells undergo a functional transition to become neurogenic during hindbrain segmentation concurrently as they maintain the progenitor cell pool. Boundary cells engaged in neurogenesis coinciding with the onset of Notch signaling, which triggered their asymmetrical cell division. Our findings reveal that distinct neurogenic phases take place during hindbrain growth and suggest that boundary cells contribute to refine the final number, identity, and proportion of neurons in the brain.

SUMMARYHindbrain boundary cells undergo a functional transition to become neurogenic concurrently as they are maintained as a progenitor cell pool. This involves a behavioral switch from symmetrically dividing progenitor cells to asymmetrically dividing progenitors, which depends on Notch-activity. DeGlyPHER: an ultrasensitive method for analysis of viral spike N-glycoforms Viruses can evade the host immune system by displaying numerous glycans on their surface "spike-proteins" that cover immune epitopes. We have developed an ultrasensitive "single pot" method to assess glycan occupancy and the extent of glycan processing from high-mannose to complex forms at each N-glycosylation site. Though aimed at characterizing glycosylation of viral spike-proteins as potential vaccines, this method is applicable for analysis of site-specific glycosylation of any glycoprotein. Revisiting regulatory decoherence: Accounting for temporal bias in a co-expression analysis reveals novel candidates controlling environmental response Environment can alter the degree of phenotypic variation and covariation, potentially influencing evolutionary trajectories. However, environment-driven changes in phenotypic variation remain understudied. In an effort to exploit the abundance of RNASequencing data now available, an increasing number of ecological studies rely on population-level correlation to characterize the plastic response of the entire transcriptome and to identify environmentally responsive molecular pathways. These studies are fundamentally interested in identifying groups of genes that respond in concert to environmental shifts. We show that population-level differential co-expression exhibits biases when capturing changes of regulatory activity and strength in rice plants responding to elevated temperature. One possible cause of this bias is regulatory saturation, the observation that detectable co-variance between a regulator and its target may be low as their transcript abundances are induced. This phenomenon appears to be particularly acute for rapid-onset environmental stressors. However, our results suggest that temporal correlations may be a reliable means to detect transient regulatory activity following rapid onset environmental perturbations such as temperature stress. Such temporal bias is likely to confound the studies of phenotypic integration, where high-order organismal traits are hypothesized to be more integrated with strong correlation under stressful conditions, while recent transcriptome studies exhibited weaker coexpression between genes under stressful conditions. Collectively, our results point to the need to account for the nuances of molecular interactions and the possibly confounding effects that these can introduce into conventional approaches to study transcriptome datasets. O-antigen structure and surface hydrophobicity determines aggregate assembly type in Pseudomonas aeruginosa Bacteria live in spatially organized aggregates during chronic infections, where they adapt to the host environment, evade immune responses and resist therapeutic interventions. Although it is known that environmental factors such as polymers influence bacterial aggregation, it is not clear how bacterial adaptation during chronic infection impacts the formation and spatial organization of aggregates in the presence of polymers. Here we show that in an in vitro model of cystic fibrosis (CF) containing the polymers eDNA and mucin, O-specific antigen is a major factor in determining the formation of two distinct aggregate assembly types of Pseudomonas aeruginosa due to alterations in cell surface hydrophobicity. Our findings suggest that during chronic infection, interplay between cell surface properties and polymers in the environment may influence the formation and structure of bacterial aggregates, which would shed new light on the fitness costs and benefits of O-antigen production in environments such as CF lungs.

ImportanceDuring chronic infection, several factors contribute to the biogeography of microbial communities. Heterogeneous populations of Pseudomonas aeruginosa form aggregates in cystic fibrosis airways, however, the impact of this population heterogeneity on spatial organization and aggregate assembly is not well understood. In this study we found that changes in O-specific antigen determine the spatial organization of P. aeruginosa cells by altering the relative cell surface hydrophobicity. This finding suggests a role for O-antigen in regulating P. aeruginosa aggregate size and shape in cystic fibrosis airways. An ex vivo human bone marrow milieu identifies CDH2 as a driver of niche-mediated treatment resistance in leukaemia Being able to clinically target the interaction of leukaemic cells with their microenvironment has remained a key therapeutic obstacle. Leukaemia cells re-program their microenvironment to provide support and protection from standard chemotherapy, molecularly targeted therapies as well as immunotherapy. To address this challenge, we have developed experimentally accessible human induced pluripotent stem cell engineered (iPSC) niches ex vivo to reveal insights into druggable cancer-niche interactions. We show that mesenchymal (iMSC) and vascular niche-like (iANG) cells support ex vivo proliferation of patient-derived leukaemia cells, impact dormancy and mediate therapy resistance. Mesenchymal stem cells protected both cycling and non- cycling blasts against dexamethasone treatment while vascular niche-like cells only protected dormant cells. Growth support and protection from dexamethasone induced- apoptosis was dependent on direct cell-cell contact and was mediated by N-cadherin (CDH2). To explore the therapeutic potential of disrupting this cell-cell interaction, we tested the CDH2 antagonist ADH-1 (previously in phase I / II for solid tumours) in a very aggressive patient-derived xenograft acute lymphoblastic leukaemia mouse model. ADH-1 showed high in vivo efficacy and combination therapy with ADH-1 and dexamethasone was superior compared with dexamethasone alone. These findings provide a proof-of-concept starting point to develop novel and potentially safer anti- cancer therapeutics that target niche-mediated cancer cell dependencies in haematological malignancies.

One Sentence SummaryCDH2 drives niche-dependent treatment resistance targetable via ADH-1, a low toxic potential candidate for future leukaemia clinical trials. 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.