Putative mobilized colistin resistance genes in the human gut microbiome BackgroundThe high incidence of bacterial genes that confer resistance to last-resort antibiotics, such as colistin caused by MCR genes, poses an unprecedented threat to our civilizations health. To understand the spread, evolution, and distribution of such genes among human populations, with the final goal of diminishing their occurrence in human environments should be a priority. To tackle this problem, we investigated the distribution and prevalence of potential mcr genes in the human gut microbiome we used a set of bioinformatics tools to screen the Unified Human Gastrointestinal Genome (UHGG) collection for the presence, synteny and phylogeny of putative mcr genes, and co-located antibiotic resistance genes.

ResultsA total of 2,079 ARGs were classified as different MCR in 2,046 Metagenome assembled genomes (MAGs), present in 1,596 individuals from 41 countries, of which 215 MCRs were identified in plasmidial contigs. The genera that presented the largest number of MCR-like genes were Suterella and Parasuterella, prevalent human gut bacteria of which Suterella wadsworthensis is associated with autism. Other potential pathogens carrying MCR genes belonged to the genus Vibrio, Escherichia and Campylobacter. Finally, we identified a total of 22,746 ARGs belonging to 21 different classes in the same 2,046 MAGs, suggesting multi-resistance potential in the corresponding bacterial strains, increasing the concern of ARGs impact in the clinical settings.

ConclusionThis study uncovers the diversity of MCR-like genes in the human gut microbiome. We showed the cosmopolitan distribution of these genes in individuals worldwide and the co-presence of other antibiotic resistance genes, including Extended-spectrum beta-lactamases (ESBL). Also, we described mcr-like genes fused to a PAP2-like domain in S. wadsworthensis. Although these novel sequences increase our knowledge about the diversity and evolution of mcr-like genes, their activity and a potential colistin resistance in the corresponding strains has to be experimentally validated.
SUFI: An automated approach to spectral unmixing of fluorescent multiplex images captured in mouse and postmortem human brain tissues Multispectral fluorescence imaging coupled with linear unmixing is a form of image data collection and analysis that uses multiple fluorescent dyes - each measuring a specific biological signal - that are simultaneously measured and subsequently "unmixed" to provide a read-out for each individual signal. This strategy allows for measuring multiple signals in a single data capture session - for example, multiple proteins or RNAs in tissue slices or cultured cells, but can often result in mixed signals and bleed-through problems across dyes. Existing spectral unmixing algorithms are not optimized for challenging biological specimens such as postmortem human brain tissue, and often require manual intervention to extract spectral signatures. We therefore developed an intuitive, automated, and flexible package called SUFI: spectral unmixing of fluorescent images (https://github.com/LieberInstitute/SUFI). This package unmixes multispectral fluorescence images by automating the extraction of spectral signatures using Vertex Component Analysis, and then performs one of three unmixing algorithms derived from remote sensing. We demonstrate these remote sensing algorithms performance on four unique biological datasets and compare the results to unmixing results obtained using ZEN Black software (Zeiss). We lastly integrate our unmixing pipeline into the computational tool dotdotdot that is used to quantify individual RNA transcripts at single cell resolution in intact tissues and perform differential expression analysis of smFISH data, and thereby provide a one-stop solution for multispectral fluorescence image analysis and quantification. In summary, we provide a robust, automated pipeline to assist biologists with improved spectral unmixing of multispectral fluorescence images.
Can fMRI reveal the representation of syntactic structure in the brain? While studying semantics in the brain, neuroscientists use two approaches. One is to identify areas that are correlated with semantic processing load. Another is to find areas that are predicted by the semantic representation of the stimulus words. However, in the domain of syntax, most studies have focused only on identifying areas correlated with syntactic processing load. One possible reason for this discrepancy is that representing syntactic structure in an embedding space such that it can be used to model brain activity is a non-trivial computational problem. Another possible reason is that it is unclear if the low signal-to-noise ratio of neuroimaging tools such as functional Magnetic Resonance Imaging (fMRI) can allow us to reveal correlates of complex (and perhaps subtle) syntactic representations. In this study, we propose novel multi-dimensional features that encode information about the syntactic structure of sentences. Using these features and fMRI recordings of participants reading a natural text, we model the brain representation of syntax. First, we find that our syntactic structure-based features explain additional variance in the brain activity of various parts of the language system, even after controlling for complexity metrics that capture processing load. At the same time, we see that regions well-predicted by syntactic features are distributed in the language system and are not distinguishable from those processing semantics. Spontaneous Chiral Symmetry Breaking in a Random Driven Chemical System Living systems have evolved to efficiently consume available energy sources using an elaborate circuitry of chemical reactions which, puzzlingly, bear a strict restriction to asymmetric chiral configurations. While autocatalysis is known to promote such chiral symmetry breaking, whether a similar phenomenon may also be induced in a more general class of configurable chemical systems--via energy exploitation--is a sensible yet underappreciated possibility. This work examines this question within a model of randomly generated complex chemical networks. We show that chiral symmetry breaking may occur spontaneously and generically by harnessing energy sources from external environmental drives. Key to this transition are intrinsic fluctuations of achiral-to-chiral reactions and tight matching of system configurations to the environmental drives, which together amplify and sustain diverged enantiomer distributions. These asymmetric states emerge through steep energetic transitions from the corresponding symmetric states and sharply cluster as highly-dissipating states. The results thus demonstrate a generic mechanism in which energetic drives may give rise to homochirality in an otherwise totally symmetrical environment, and from an early-life perspective, might emerge as a competitive, energy-harvesting advantage. Retraction ATPase motors from three orthologous type IVa pilus systems support promiscuous retraction of the Vibrio cholerae competence pilus Bacterial surface appendages called type IVa pili (T4aP) promote diverse activities including DNA uptake, twitching motility, and virulence. These activities rely on the ability of T4aP to dynamically extend and retract from the cell surface. Dynamic extension relies on a motor ATPase commonly called PilB. Most T4aP also rely on specific motor ATPases, commonly called PilT and PilU, to dynamically and forcefully retract. Here, we systematically assess whether motor ATPases from three orthologous T4aP can functionally complement Vibrio cholerae mutants that lack their endogenous motors. We found that the PilT and PilU retraction ATPases from the three T4aP systems tested are promiscuous and promote retraction of the V. cholerae competence T4aP despite a high degree of sequence divergence. In contrast, the orthologous extension ATPases from the same T4aP systems were not able to mediate extension of the V. cholerae competence T4aP despite exhibiting a similar degree of sequence divergence. Also, we show that one of the PilT orthologs characterized does not support PilU-dependent retraction and provide some data to indicate that the C-terminus of PilT is important for PilU-dependent retraction. Together, our data suggest that retraction ATPases may have maintained a high degree of promiscuity for promoting retraction of T4aP, while extension ATPases may have evolved to become specific for their cognate systems.

IMPORTANCEOne way that bacteria interact with their environments is via hair-like appendages called type IVa pili (T4aP). These appendages dynamically extend and retract from the cell surface via the action of distinct ATPase motors. T4aP are present in diverse bacterial species. Here, we demonstrate that retraction motors from three T4aP are promiscuous, and capable of promoting retraction of a heterologous T4aP system. By contrast, the extension ATPase motors from these same T4aP systems are specific and cannot promote extension of a heterologous T4aP. Thus, these results suggest that T4aP extension may be more tightly regulated compared to T4aP retraction. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. Superinfection exclusion: a viral strategy with short-term benefits and long-term drawbacks Viral superinfection occurs when multiple viral particles subsequently infect the same host. In nature, several viral species are found to have evolved diverse mechanisms to prevent superinfection (superinfection exclusion) but how this strategic choice impacts the fate of mutations in the viral population remains unclear. Using stochastic simulations, we find that genetic drift is suppressed when superinfection occurs, thus facilitating the fixation of beneficial mutations and the removal of deleterious ones. Interestingly, we also find that the competitive (dis)advantage associated with variations in life history parameters is not necessarily captured by the viral growth rate for either infection strategy. Putting these together, we then show that a mutant with superinfection exclusion will easily overtake a superinfecting population even if the latter has a much higher growth rate. Our findings suggest that while superinfection exclusion can negatively impact the long-term adaptation of a viral population, in the short-term it is ultimately a winning strategy.

Author summaryViral social behaviour has recently been receiving increasing attention in the context of ecological and evolutionary dynamics of viral populations. One fascinating and still relatively poorly understood example is superinfection or co-infection, which occur when multiple viruses infect the same host. Among bacteriophages, a wide range of mechanisms have been discovered that enable phage to prevent superinfection (superinfection exclusion) even at the cost of using precious resources for this purpose. What is the evolutionary impact of this strategic choice and why do so many phages exhibit this behaviour? Here, we conduct an extensive simulation study of a phage population to address this question. In particular, we investigate the fate of viral mutations arising in an environment with a constant supply of bacterial hosts designed to mimic a "turbidostat," as these are increasingly being used in laboratory evolution experiments. Our results show that allowing superinfection in the long-term yields a population which is more capable of adapting to changes in the environment. However, when in direct competition, mutants capable of preventing superinfection experience a very large advantage over their superinfecting counterparts, even if this ability comes at a significant cost to their growth rate. This indicates that while preventing superinfection can negatively impact the long-term prospects of a viral population, in the short-term it is ultimately a winning strategy. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. Sarcoplasmic reticulum calcium handling in young DBA/2J vs C57BL/10 mdx models of Duchenne muscular dystrophy The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. 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. An alternative UPF1 isoform drives conditional remodeling of nonsense-mediated mRNA decay The nonsense-mediated mRNA decay (NMD) pathway monitors translation termination to degrade transcripts with premature stop codons and regulate thousands of human genes. Here we show that an alternative mammalian-specific isoform of the core NMD factor UPF1, termed UPF1LL, enables condition-dependent remodeling of NMD specificity. Previous studies indicate that the extension of a conserved regulatory loop in the UPF1LL helicase core confers a decreased propensity to dissociate from RNA upon ATP hydrolysis relative to the major UPF1 isoform, designated UPF1 SL. Using biochemical and transcriptome-wide approaches, we find that UPF1LL overcomes the protective RNA binding proteins PTBP1 and hnRNP L to preferentially bind and down-regulate long 3UTRs normally shielded from NMD. Unexpectedly, UPF1LL supports induction of NMD on new populations of substrate mRNAs in response to activation of the integrated stress response and impaired translation efficiency. Thus, while canonical NMD is abolished by moderate translational repression, UPF1LL activity is enhanced, providing a mechanism to rapidly rewire NMD specificity in response to cellular stress. Theory and Simulations of condensin mediated loop extrusion in DNA Condensation of hundreds of mega-base-pair-long human chromosomes in a small nuclear volume is a spectacular biological phenomenon. This process is driven by the formation of chromosome loops. The ATP consuming motor, condensin, interacts with chromatin segments to actively extrude loops. Motivated by real-time imaging of loop extrusion (LE), we created an analytically solvable model, predicting the LE velocity and step size distribution as a function of external load. The theory fits the available experimental data quantitatively, and suggests that condensin must undergo a large conformational change, induced by ATP binding, bringing distant parts of the motor to proximity. Simulations using a simple model confirm that the motor transitions between an open and a closed state in order to extrude loops by a scrunching mechanism, similar to that proposed in DNA bubble formation during bacterial transcription. Changes in the orientation of the motor domains are transmitted over ~ 50 nm, connecting the motor head and the hinge, thus providing an allosteric basis for LE. The relative speed of the glucorticoid stress response varies independently of scope and is predicted by environmental variability and longevity across birds. The acute glucocorticoid response is a key mediator of the coordinated vertebrate response to unpredictable challenges. Rapid glucocorticoid increases initiate changes that allow animals to cope with stressors. The scope of the glucocorticoid response -- defined here as the absolute increase in glucocorticoids -- is associated with individual differences in performance and varies across species with environment and life history. In addition to varying in scope, responses can differ enormously in speed; however, relatively little is known about whether speed and absolute glucocorticoid levels covary, how selection shapes speed, or what aspects of speed are important. We used corticosterone samples collected at 5 time points from 1,750 individuals of 60 species of birds to ask i) how the speed and scope of the glucocorticoid response covary and ii) whether variation in absolute or relative speed is predicted by environmental context or life history. As predicted by a recent model, faster absolute glucocorticoid responses were strongly associated with a larger scope. Despite this covariation, the relative speed of the glucocorticoid response (standardized within species) varied independently of absolute scope, suggesting that selection could operate on both features independently. Species with faster relative glucocorticoid responses lived in locations with more variable temperature and had shorter lifespans. Our results suggest that rapid changes associated with the speed of the glucocorticoid response, such as those occurring through non-genomic receptors, might be an important determinant of coping ability and we emphasize the need for studies designed to measure speed independently of absolute glucocorticoid levels. Rad51 filaments assembled in the absence of the complex formed by the Rad51 paralogs Rad55 and Rad57 are outcompeted by translesion DNA polymerases on UV-induced ssDNA gaps The bypass of DNA lesions that block replicative polymerases during DNA replication relies on DNA damage tolerance pathways. The error-prone translesion synthesis (TLS) pathway relies on specialized DNA polymerases that incorporate nucleotides in front of base lesions, potentially inducing mutagenesis. Two error-free pathways can bypass the lesions: the template switching pathway, which uses the sister chromatid as a template, and the homologous recombination pathway (HR), which also can use the homologous chromosome as template. The balance between error-prone and error-free pathways controls the mutagenesis level. Therefore, it is crucial to precisely characterize factors that influence the pathway choice to better understand genetic stability at replication forks. In yeast, the complex formed by the Rad51 paralogs Rad55 and Rad57 promotes HR and template-switching at stalled replication forks. At DNA double-strand breaks (DSBs), this complex promotes Rad51 filament formation and stability, notably by counteracting the Srs2 anti-recombinase. To explore the role of the Rad55-Rad57 complex in error-free pathways, we monitored the genetic interactions between Rad55-Rad57, the translesion polymerases Pol{zeta} or Pol{eta}, and Srs2 following UV radiation that induces mostly single-strand DNA gaps. We found that the Rad55-Rad57 complex was involved in three ways. First, it protects Rad51 filaments from Srs2, as it does at DSBs. Second, it promotes Rad51 filament stability independently of Srs2. Finally, we observed that UV-induced HR is almost abolished in Rad55-Rad57 deficient cells, and is partially restored upon Pol{zeta} or Pol{eta} depletion. Hence, we propose that the Rad55-Rad57 complex is essential to promote Rad51 filament stability on single-strand DNA gaps, notably to counteract the error-prone TLS polymerases and mutagenesis.

Author summaryProcessive and accurate DNA polymerases replicate genomic DNA during the S phase of each cell cycle. DNA base lesions on template DNA block these polymerases and result in an accumulation of single-stranded DNA gaps behind moving replication forks. These gaps are filled-in by error-prone and error-free pathways. In this work, we show that the complex made by the Rad51 paralogs Rad55 and Rad57 is essential for the error-free homologous recombination gap-filling pathway when DNA replication is stalled by UV-induced DNA lesions, but not for DNA double strand break repair. Interestingly, we found that homologous recombination is efficiently outcompeted by error-prone translesion DNA polymerases in Rad55-Rad57-deficient cells. We propose that the Rad55-Rad57 complex is essential for Rad51 filament stability at UV-induced DNA gaps to promote efficient error-free homologous recombination. Furthermore, our study in yeast predicts that inhibitors of error-prone DNA polymerases might selectively target cancer cells in which RAD51 paralogs are mutated. Soluble ECM promotes organotypic formation in lung alveolar model Micropatterned suspension culture creates consistently sized and shaped cell aggregates but has not produced organotypic structures from stable cells, thus restricting its use in accurate disease modeling. Here, we show that organotypic structure is achieved in hybrid suspension culture via supplementation of soluble extracellular matrix (ECM). We created a viable lung organoid from epithelial, endothelial, and fibroblast human stable cell lines in suspension culture. We demonstrate the importance of soluble ECM in organotypic patterning with the emergence of lumen-like structures with airspace showing feasible gas exchange units, formation of branching, perfusable vasculature, and long-term 70-day maintenance of lumen structure. Our results show a dependent relationship between enhanced fibronectin fibril assembly and the incorporation of ECM in the organoid. We successfully applied this technology in modeling lung fibrosis via bleomycin induction and test a potential antifibrotic drug in vitro while maintaining fundamental cell-cell interactions in lung tissue. Our human fluorescent lung organoid (hFLO) model represents features of pulmonary fibrosis which were ameliorated by fasudil treatment. We also demonstrate a 3D culture method with potential of creating organoids from mature cells, thus opening avenues for disease modeling and regenerative medicine, enhancing understanding of lung cell biology in health and lung disease. PKR-mediated Stress Response Enhances Dengue and Zika Virus Replication The mechanisms by which Flaviviruses use non-canonical translation to support their replication in host cells are largely unknown. Here we investigated how the integrated stress response (ISR), which promotes translational arrest by eIF2 phosphorylation (p-eIF2), regulates Flavivirus replication. During Dengue virus (DENV) and Zika virus (ZIKV) infection, eIF2 activation peaked at 24 hours post infection and was dependent on PKR but not type I interferon. The ISR is activated downstream of p-eIF2 during infection with either virus, but translation arrest only occurred following DENV infection. Despite this difference, both DENV and ZIKV replication was impaired in cells lacking PKR, independent of IFN-I/NF-kB signaling or cell viability. By using a ZIKV 5UTR reporter system as a model, we found that this region of the genome is sufficient to promote an enhancement of viral mRNA translation in the presence of an active ISR. Together we provide evidence that Flaviviruses escape ISR translational arrest and co-opt this response to increase viral replication. 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. Genomic sequencing confirms absence of introgression despite past hybridisation between a critically endangered bird and its common congener Genetic swamping resulting from interspecific hybridisation can increase extinction risk for threatened species. The development of high-throughput and reduced-representation genomic sequencing and analyses to generate large numbers of high resolution genomic markers has the potential to reveal introgression previously undetected using small numbers of genetic markers. However, few studies to date have implemented genomic tools to assess the extent of interspecific hybridisation in threatened species. Here we investigate the utility of genome-wide single nucleotide polymorphisms (SNPs) to detect introgression resulting from past interspecific hybridisation in one of the worlds rarest birds. Anthropogenic impacts have resulted in hybridisation and subsequent backcrossing of the critically endangered Aotearoa New Zealand endemic kak[i] (black stilts; Himantopus novaezelandiae) with the non-threatened self-introduced congeneric poaka (Aotearoa New Zealand population of pied stilts, Himantopus himantopus leucocephalus), yet genetic analyses with a limited set of microsatellite markers revealed no evidence of introgression of poaka genetic material in kak[i], excluding one individual. We use genomic data for [~]63% of the wild adult kak[i] population to reassess the extent of introgression resulting from hybridisation between kak[i] and poaka. Consistent with previous genetic analyses, we detected no introgression from poaka into kak[i]. These collective results indicate that, for kak[i], existing microsatellite markers provide a robust, cost-effective approach to detect cryptic hybrids. Further, for well-differentiated species, the use of genomic markers may not be required to detect admixed individuals. Darwin's small and medium ground finches might have taste preferences, but not for human foods Urbanization is rapidly changing ecological niches. On the inhabited Galapagos Islands, Darwins finches consume human-introduced foods preferentially; however, it remains unclear why. Here we presented pastry with flavour profiles typical of human foods (oily, salty, sweet) to small ground finches (Geospiza fuliginosa) and medium ground finches (Geospiza fortis) to test if latent taste preferences might drive selection of human foods. If human-food flavours were consumed more than a neutral or bitter control only at sites with human foods, then we predicted tastes were acquired after urbanization; however, if no site-differences were found then this would indicate latent taste preferences. Contrary to both predictions, we found little evidence that human-food flavours were preferred compared to control flavours at any site. Instead, finches showed a weak aversion to oily foods, but only at remote (no human foods present) sites. This was further supported by behavioural responses, with beak wiping occurring more often at remote sites after finches tasted flavours associated with human foods. Our results suggest, therefore, that while Darwins finches may have acquired a tolerance to the flavours of human food, latent taste preferences are unlikely to have played a major role in their dietary response to increased urbanization. Multi-omic analysis of the Arabidopsis clock activator mutant rve 4 6 8 reveals connections to carbohydrate metabolism and proteasome regulation Plants are able to sense changes in their light environments, such as the onset of day and night, as well as anticipate these changes in order to adapt and survive. Central to this ability is the plant circadian clock, a molecular circuit that precisely orchestrates plant cell processes over the course of a day. REVEILLE proteins (RVEs) are recently discovered members of the plant circadian circuitry that activate the evening complex and PRR genes to maintain regular circadian oscillation. The RVE 8 protein and its two homologs, RVE 4 and 6, have been shown to limit the length of the circadian period, with rve 4 6 8 triple-knockout plants possessing an elongated period along with increased leaf surface area, biomass, cell size and delayed flowering relative to wild-type Col-0 plants. Here, using a multi-omics approach consisting of phenomics, transcriptomics, proteomics, and metabolomics we draw novel connections between RVE8-like proteins and a number of core plant cell processes. In particular, we reveal that loss of RVE8-like proteins results in altered carbohydrate, organic acid and lipid metabolism, including a starch excess phenotype at dawn. We further demonstrate that rve 4 6 8 plants have lower levels of 20S proteasome subunits and possess significantly reduced proteasome activity, potentially explaining the increase in cell-size observed in RVE8-like mutants. Overall, this robust, multi-omic dataset, provides substantial new insights into the far reaching impact RVE8-like proteins have on the diel plant cell environment. Integrated analysis of cervical squamous cell carcinoma cohorts from three continents reveals conserved subtypes of prognostic significance. Human papillomavirus (HPV)-associated cervical cancer represents one of the leading causes of cancer death worldwide. Although low-middle income countries are disproportionately affected, our knowledge of the disease predominantly originates from populations in high-income countries. Using the largest multi-omic analysis of cervical squamous cell carcinoma (CSCC) to date, totalling 643 tumours and representing patient populations from the USA, Europe and Sub-Saharan Africa, we identify two CSCC subtypes (C1 and C2) with differing prognosis. C1 tumours are largely HPV16-driven, display increased cytotoxic T-lymphocyte infiltration and frequently harbour PIK3CA and EP300 mutations. C2 tumours are associated with shorter overall survival, are frequently driven by HPVs from the HPV18-containing alpha-7 clade, harbour alterations in the Hippo signalling pathway and increased expression of immune checkpoint genes, B7-H3 (also known as CD276) and NT5E (also known as CD73) and PD-L2 (also known as PDCD1LG2). In conclusion, we identify two novel, therapy-relevant CSCC subtypes that share the same defining characteristics across three geographically diverse cohorts. Single-cell Transcriptomic Profiling Unveils Cardiac Cell-type Specific Response to Maternal Hyperglycemia Underlying the Risk of Congenital Heart Defects Congenital heart disease (CHD) is the most prevalent structural malformations of the heart affecting [~]1% of live births. To date, both damaging genetic variations and adverse environmental exposure such as maternal diabetes have been found to cause CHD. Clinical studies show [~]fivefold higher risk of CHD in the offspring of mothers with pregestational diabetes. Maternal pregestational diabetes affects the gene regulatory networks key to proper cardiac development in the fetus. However, the cell-type specificity of these gene regulatory responses to maternal diabetes and their association with the observed cardiac defects in the fetuses remains unknown. To uncover the transcriptional responses to maternal diabetes in the early embryonic heart, we used an established murine model of pregestational diabetes. In this model, we have previously demonstrated an increased incidence of CHD. Here, we show maternal hyperglycemia (matHG) elicits diverse cellular responses during heart development by single-cell RNA-sequencing in embryonic hearts exposed to control and matHG environment. Through differential gene-expression and pseudotime trajectory analyses of this data, we identified changes in lineage specifying transcription factors, predominantly affecting Isl1+ second heart field progenitors and Tnnt2+ cardiomyocytes with matHG. Using in vivo cell-lineage tracing studies, we confirmed that matHG exposure leads to impaired second heart field-derived cardiomyocyte differentiation. Finally, this work identifies matHG-mediated transcriptional determinants in cardiac cell lineages elevate CHD risk and show perturbations in Isl1-dependent gene-regulatory network (Isl1-GRN) affect cardiomyocyte differentiation. Functional analysis of this GRN in cardiac progenitor cells will provide further mechanistic insights into matHG-induced severity of CHD associated with diabetic pregnancies. Sablefish (Anoplopoma fimbria Pallas, 1814) plasma biochemistry and hematology reference intervals including blood cell morphology Plasma biochemistry and hematology reference intervals are integral health assessment tools in all medical fields, including aquatic animal health. As sablefish (Anoplopoma fimbria) are becoming aquaculturally and economically more important, this manuscript provides essential reference intervals (RI) for their plasma biochemistry and hematology along with reference photomicrographs of blood cells in healthy, fasted sablefish. Blood cell morphology can differ between fish species. In addition, blood cell counts and blood chemistry can vary between fish species, demographics, water conditions, seasons, diets, and culture systems, which precludes the use of RIs from other fish species. For this study, blood was collected for plasma biochemistry and hematology analysis between June 20 and July 18, 2019, from healthy, yearling sablefish, hatched and reared in captivity on a commercial diet. Overnight fast of 16-18 hours did not sufficiently reduce lipids in the blood, which led to visible lipemia and frequent rupture of blood cells during analysis. Therefore, sablefish should be fasted for 24 to 36 hours before blood is collected to reduce hematology artifacts or possible reagent interference in plasma biochemistry analysis. Lymphocytes were the most dominant leukocytes (98%), while eosinophils were rare, and basophils were not detected in sablefish. Neutrophils were very large cells with Dohle bodies. In mammals and avian species, Dohle bodies are usually signs of toxic change from inflammation, but no such association was found in these fish. In conclusion, lipemia can interfere with sablefish blood analysis, and available removal methods should be evaluated as fasting for up to 36 h might not always be feasible. Also, more studies are required to establish RI for different developmental stages and rearing conditions. Communities of tropical soil fungi differ between burned and unburned forest, with corresponding changes in plant community composition, litter and soil chemistry Fire is predicted to be more severe and frequent in forests of the Australian Monsoon Tropics over the coming decades. The way in which groups of ecologically important soil fungi respond to disturbance caused by fire has not been studied in Australian tropical forest ecosystems. Ectomycorrhizal (EM) fungi are important tree symbionts and saprotrophic fungi drive soil nutrient cycles. We analysed both publicly-available environmental DNA sequence data as well as soil chemistry data to test a hypothesis that burned areas in a contiguous tropical forest have different community composition and diversity of EM and saprotrophic soil fungi relative to nearby unburned sites. We tested this hypothesis by measuring community-level taxonomic composition, fungal diversity, species richness and evenness. We determined whether changes in fungal communities were associated with fire-altered soil chemical/physical properties, vegetation types, or the direct effect of fire. Soil fungi differed in abundance and community phylogenetic structure between forest sites that had experienced fire, and sites dominated by unburned forest. EM community composition differed between unburned and burned sites, which had more mycorrhizal hosts including Corymbia intermedia, Acacia flavescens and Acacia midgleyi. Lower diversity of saprotrophic fungi was correlated with lower soil nutrient levels and different litter composition at burned sites. Pyrophilic, truffle-like EM fungi that rely on mycophagous mammals for dispersal were abundant at recently burned sites. We conclude that EM fungi show different patterns of diversity in burned tropical forest, likely driven by changing plant communities, whereas differences in saprotrophic fungal communities of burned sites may be driven by by reduced soil nutrient levels and altered litter composition.

Graphical abstract O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY Multi-omics investigation of tacrolimus nephrotoxicity BackgroundTacrolimus, an immunosuppressive drug prescribed to a majority of transplanted patients is nephrotoxic, through still unclear mechanisms. This study aims to evaluate the impact of tacrolimus on a lineage of proximal tubular cells using a multi-omics approach.

MethodsLLC-PK1 cells were exposed to 5 M of tacrolimus for 24h. Intracellular proteins and metabolites, and extracellular metabolites were extracted and analysed by LC-MS/MS. The transcriptional expression of the dysregulated proteins PCK-1, FBP1 and FBP2 was measured using RT-qPCR.

ResultsIn our cell model, tacrolimus impacted different metabolic pathways including those of arginine (e.g., citrulline, ornithine) (p < 0.0001), amino acids (e.g., valine, isoleucine, aspartic acid) (p < 0.0001) and pyrimidine (p<0.01). In addition, it induced oxidative stress (p < 0.01) as shown by a decrease in total cell glutathione quantity. It impacted cell energy through an increase in Krebs cycle intermediates (e.g., citrate, aconitate, fumarate) (p < 0.01) and down-regulation of PCK-1 (p < 0.05) and FPB1 (p < 0.01), which are key enzymes in gluconeogenesis. Apart from glucose synthesis, gluconeogenesis is an important process in kidney mediated acid-base balance control.

ConclusionThe variations found using this multi-omics approach clearly point towards a dysregulation of energy production in epithelial cells of the renal tubule, and potentially of their functions, that may be implicated in tacrolimus nephrotoxicity in the clinics. Multiparametric Platform for Profiling Lipid Trafficking in Human Leukocytes: Application for Hypercholesterolemia Systematic insight into cellular dysfunctions can improve understanding of disease etiology, risk assessment and patient stratification. We present a multiparametric high-content imaging platform enabling quantification of low-density lipoprotein (LDL) uptake and lipid storage in cytoplasmic droplets of primary leukocyte subpopulations. We validated this platform with samples from 65 individuals with variable blood LDL-cholesterol (LDL-c) levels, including familial hypercholesterolemia (FH) and non-FH subjects. We integrated lipid storage data into a novel readout, lipid mobilization, measuring the efficiency with which cells deplete lipid reservoirs. Lipid mobilization correlated positively with LDL uptake and negatively with hypercholesterolemia and age, improving differentiation of individuals with normal and elevated LDL-c. Moreover, combination of cell-based readouts with a polygenic risk score for LDL-c explained hypercholesterolemia better than the genetic risk score alone. This platform provides functional insights into cellular lipid trafficking from a few mls of blood and is applicable to dissect metabolic disorders, such as hypercholesterolemia.

MotivationWe have limited information on how cellular lipid uptake and processing differ between individuals and influence the development of metabolic diseases, such as hypercholesterolemia. Available assays are labor intensive, require skilled personnel and are difficult to scale to higher throughput, making it challenging to obtain systematic functional cell-based data from individuals. To overcome this problem, we established a scalable automated analysis pipeline enabling reliable quantification of multiple cellular readouts, including lipid uptake, storage and mobilization, from different white blood cell populations. This approach provides new personalized insights into the cellular basis of hypercholesterolemia and obesity.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/440471v4_ufig1.gif" ALT="Figure 1">
View larger version (84K):
org.highwire.dtl.DTLVardef@9222c1org.highwire.dtl.DTLVardef@27e289org.highwire.dtl.DTLVardef@89ba82org.highwire.dtl.DTLVardef@33c9cd_HPS_FORMAT_FIGEXP  M_FIG C_FIG