Phenotypic plasticity is aligned with phenological adaptation on both micro- and macroevolutionary timescales In seasonally-variable environments, phenotypic plasticity in phenology may be critical for adaptation to fluctuating environmental (temperature) conditions. Using an 18-generation longitudinal dataset from natural insect (damselfly) populations, we show that phenology has strongly advanced. Individual fitness data suggest this is likely an adaptive response towards a thermally-dependent fitness optimum. A laboratory experiment revealed that developmental plasticity qualitatively matches environmental-dependence of selection, partially explaining observed phenological advance. Expanding our analysis to the macroevolutionary level, we use a database of over 1-million occurrence records and spatiotemporally-matched temperature data from 49 Swedish Odonate species to infer macroevolutionary dynamics of phenology. Phenological plasticity was more closely aligned with adaptation for species that have recently colonized northern latitudes, with more phenological mismatch at lower latitudes. Our results show that phenological plasticity plays a key role in microevolutionary dynamics within a single species, and such plasticity may have facilitated post-Pleistocene range expansion in this insect clade.
Gene duplications and phylogenomic conflict underlie major pulses of phenotypic evolution in gymnosperms Inferring the intrinsic and extrinsic drivers of species diversification and phenotypic disparity across the Tree of Life is a major challenge in evolutionary biology. In green plants, polyploidy (or whole-genome duplication, WGD) is known to play a major role in microevolution and speciation1, but the extent to which WGD has shaped macroevolutionary patterns of diversification and phenotypic innovation across plant phylogeny remains an open question. Here we examine the relationship of various facets of genomic evolution--including gene and genome duplication, genome size, and chromosome number--with macroevolutionary patterns of phenotypic innovation, species diversification, and climatic occupancy in gymnosperms. We show that genomic changes, such as WGD and genome-size shifts, underlie the origins of most major extant gymnosperm clades, and notably our results support an ancestral WGD in the gymnosperm lineage. Spikes of gene duplication typically coincide with major spikes of phenotypic innovation, while increased rates of phenotypic evolution are typically found at nodes with high gene-tree conflict, representing historic population-level dynamics during speciation. Most shifts in gymnosperm diversification since the rise of angiosperms are decoupled from putative WGDs and instead are associated with increased rates of climatic occupancy evolution, particularly in cooler and/or more arid climatic conditions, suggesting that ecological opportunity, especially in the later Cenozoic, and environmental heterogeneity have driven a resurgence of gymnosperm diversification. Our study provides critical insight on the processes underlying diversification and phenotypic evolution in gymnosperms, with important broader implications for the major drivers of both micro- and macroevolution in plants.
Vapor exposure to Δ9-tetrahydrocannabinol (THC) slows locomotion of the Maine Lobster (Homarus americanus) RationaleDespite a long history of use in synaptic physiology, the lobster has been a neglected model for behavioral pharmacology. A restaurateur proposed that exposing lobster to cannabis smoke reduces anxiety and pain during the cooking process. It is unknown if lobster gill respiration in air would result in significant {Delta}9-tetrahydrocannabinol (THC) uptake and whether this would have any detectable behavioral effects.

ObjectiveThe primary goal was to determine tissue THC levels in the lobster after exposure to THC vapor. Secondary goals were to determine if THC vapor altered locomotor behavior or nociception.

MethodsTissue samples were collected (including muscle, brain and hemolymph) from Homarus americanus (N=3 per group) following 30 or 60 minutes of exposure to vapor generated by an e-cigarette device using THC (100 mg/mL in a propylene glycol vehicle). Separate experiments assessed locomotor behavior and hot water nociceptive responses following THC vapor exposure.

ResultsTHC vapor produced duration-related THC levels in all tissues examined. Locomotor activity was decreased (distance, speed, time-mobile) by 30 min inhalation of THC. Lobsters exhibit a temperature-dependent withdrawal response to immersion of tail, antennae or claws in warm water; this is novel evidence of thermal nociception for this species. THC exposure for 60 minutes had only marginal effect on nociception under the conditions assessed.

ConclusionsVapor exposure of lobsters, using an e-cigarette based model, produces dose-dependent THC levels in all tissues and reduces locomotor activity. Hot water nociception was temperature dependent, but only minimal anti-nociceptive effect of THC exposure was confirmed. The use of principal component and factor analysis to measure fundamental cognitive processes in neuropsychological data For years, dissociation studies on neurological single cases with brain lesions were the dominant method to infer fundamental cognitive functions in neuropsychology. In contrast, the association between deficits was considered to be of less epistemological value and even misleading. Still, principal component analysis (PCA) - an associational method for dimensionality reduction - recently became popular for the identification of fundamental functions. The current study evaluated the ability of PCA and factor analysis (FA) to overcome the association problem in behavioural data of neurological patients and to identify the fundamental variables underlying a battery of measures. Synthetic data were simulated to resemble neuropsychological data with typical dissociation patterns and, thereby, typical patterns of dependence between variables. In most experiments, PCA and FA succeeded to measure the underlying target variables with high up to almost perfect precision. However, this success was fragile and relied on the success of factor rotation, which failed its intended purpose when no test scores existed that primarily measured each underlying target variable. Further, commonly used strategies to estimate the number of meaningful factors appear to underfactor neuropsychological data, thereby consistently underestimating the dimensionality of the data. Finally, simulations suggested a high potential of PCA to denoise data, with factor rotation providing an additional filter function. This can be invaluable in neuropsychology, where measures are often inherently noisy, and PCA can be superior to common compound measures, such as the arithmetic mean, in the measurement of variables with high reliability. In summary, PCA and FA appear to be powerful tools in neuropsychology that are well capable to infer fundamental cognitive functions with high precision, but the typical structure of neuropsychological data and limited informative value of associations in neuropsychology place clear limitations and a risk of a complete failure on the methods. Strong immobilization of C. elegans en masse on cultivation plates by thermoelectric cooling for high-throughput submicron-resolution microscopy Despite its profound impact on biology, most high-resolution in vivo microscopy approaches remain low throughput because current immobilization techniques require significant manual effort. We greatly accelerate imaging of the nematode Caenorhabditis elegans by implementing a simple cooling approach to easily immobilize entire populations directly on their cultivation plates. We optimize and characterize cooling immobilization. Counterintuitively, relatively warmer temperatures immobilize animals significantly more effectively than colder temperatures utilized in prior studies. This enhanced immobilization enables clear submicron-resolution fluorescence imaging, which is challenging to achieve with most current immobilization techniques. We demonstrate 64x magnification 3D imaging and timelapse recording of neurons in adults and embryos without motion blur. Compared to standard azide immobilization, cooling immobilization reduces the animal preparation and recovery time by >98%, significantly increasing experimental speed. By obviating individual animal manipulation, our approach could also empower automated imaging of large C. elegans populations within standard experimental setups and workflows. Naturally occurring combinations of receptors from single cell transcriptomics in endothelial cells VEGF inhibitor drugs have been successful, especially in ophthalmology, but not all patients respond to them. Combinations of drugs are likely to be needed for more effective therapies of angiogenesis-related diseases. In this paper we describe naturally occurring combinations of receptors in endothelial cells that might help to understand how cells communicate and to identify targets for drug combinations. We also develop and share a new software tool called DECNEO to identify them.

Single-cell gene expression data are used to identify a set of co-expressed endothelial cell receptors, conserved among species (mice and humans) and enriched, within a network, of connections to up-regulated genes. This set includes several receptors previously shown to play a role in angiogenesis. Multiple statistical tests from large datasets, including an independent validation set, support the reproducibility, evolutionary conservation and role in angiogenesis of these naturally occurring combinations of receptors. We also show tissue-specific combinations and, in the case of choroid endothelial cells, consistency with both well-established and recent experimental findings, presented in a separate paper.

The results and methods presented here advance the understanding of signaling to endothelial cells. The methods are generally applicable to the decoding of intercellular combinations of signals. Elucidation of SARS-CoV-2 budding mechanisms through molecular dynamics simulations of M and E protein complexes SARS-CoV-2 and other coronaviruses pose major threats to global health, yet computational efforts to understand them have largely overlooked the process of budding, a key part of the coronavirus life cycle. When expressed together, coronavirus M and E proteins are sufficient to facilitate budding into the ER-Golgi intermediate compartment (ERGIC). To help elucidate budding, we ran atomistic molecular dynamics (MD) simulations using the Feig laboratorys refined structural models of the SARS-CoV-2 M protein dimer and E protein pentamer. Our MD simulations consisted of M protein dimers and E protein pentamers in patches of membrane. By examining where these proteins induced membrane curvature in silico, we obtained insights around how the budding process may occur. Multiple M protein dimers acted together to induce global membrane curvature through protein-lipid interactions while E protein pentamers kept the membrane planar. These results could eventually help guide development of antiviral therapeutics which inhibit coronavirus budding. Hnf1b-CreER causes efficient recombination of a Rosa26-RFP reporter in duct and islet δ cells The Hnf1b-CreERT2 BAC transgenic (Tg(Hnf1b-cre/ERT2)1Jfer) has been used extensively to trace the progeny of pancreatic ducts in developmental, regeneration, or cancer models. Hnf1b-CreERT2 transgenics have been used to show that cells from that form a duct-like plexus in the embryonic pancreas are bipotent duct-endocrine progenitors, whereas adult mouse duct cells are not a common source of {beta} cells in various regenerative settings. The interpretation of such genetic lineage tracing studies is critically dependent on a correct understanding of the cell type specificity of recombinase activity with each reporter system. We have re-examined the performance of Hnf1b-CreERT2 with a Rosa26-RFP reporter transgene. This showed inducible recombination of up to 96% adult duct cells, a much higher efficiency than previously used reporter transgenes. Despite this high duct-cell excision, recombination in  and {beta} cells remained very low, similar to previously used reporters. However, nearly half of somatostatin-expressing {delta} cells showed reporter activation, which was due to Cre expression in {delta} cells rather than to duct to {delta} cell conversions. The high recombination efficiency in duct cells indicates that the Hnf1b-CreERT2 model can be useful for both ductal fate mapping and genetic inactivation studies. The recombination in {delta} cells does not modify the interpretation of studies that failed to show duct conversions to other cell types, but needs to be considered if this model is used in studies that aim to modify the plasticity of pancreatic duct cells. Protein reservoirs of seeds are composites of amyloid and amyloid-like structures facilitating germination The functional roles of amyloids are known across diverse species, but their natural presence in plants remains underexplored. Based on the structure and function of protein bodies found in other organisms, we hypothesize that the seed storage protein bodies (SSPB) may have amyloid structures for controlling physiological functions. Here we show that monocot and dicot SSPB exhibit a speckled-pattern of amyloids interspersed in an amyloid-like matrix, suggesting their composite nature. This is confirmed by amyloid-specific probes, amyloidogenic region scans, infrared signatures, reconstitution studies, and differential degradation during germination. The role of amyloidic composites in seed germination is proved by amyloid content, degradation inhibition and relevant germination parameters, using in-vitro SSPB, ex-vivo protoplasts and in-situ seed sections. The results open new research questions in understanding molecular details of SSPB assembly-disassembly and utilization, in addition to evolutionary roles of amyloid structure and function. Contrasting resistance and resilience to light variation of the coupled oxic and anoxic componentsof an experimental microbial ecosystem Understanding how microbial communities of aquatic ecosystems respond to environmental change remains a critical challenge in microbial ecology. In this study, we used phototrophic oxic-anoxic micro-ecosystems to understand how the functioning and diversity of aerobic and anaerobic lake analog communities is affected by light deprivation. Continuous measurements were performed to describe oxygen dynamics (mean/min/max/amplitude) and time-series of full-length 16S rRNA sequencing were used to quantify changes in alpha- and beta-diversity. In the top oxic layer, oxygen concentration decreased significantly under light deprivation, but showed resilience in mean, minimum and maximum after light conditions were restored. Only the amplitude of diurnal fluctuations in oxygen concentrations did not recover fully, and instead tended to remain lower in treated ecosystems. Alpha-diversity of the top oxic layer communities showed a delayed increase after light conditions were restored, and was not resilient. In contrast, alpha-diversity of the anoxic bottom layer communities increased due to the stressor, but was resilient. Community composition changed significantly during light deprivation, showed resilience in the oxic layer and lack of resilience in the anoxic layer. Alpha-diversity and the amplitude of oxygen within and among treatments were strongly correlated, suggesting that higher diversity could lead to less variable oxygen concentrations, or vice versa. Our experiment showed that light deprivation induces multifaceted responses of community function and structure, hence focusing on a single stability component could potentially be misleading. 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. Fatal COVID-19 and non-COVID-19 Acute Respiratory Distress Syndrome is Associated with Incomplete Alveolar Type 1 Epithelial Cell Differentiation from the Transitional State Without Fibrosis ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. Human papillomavirus integration transforms chromatin to drive oncogenesis Human papillomavirus (HPV) drives almost all cervical cancers and up to [~]70% of head and neck cancers. Frequent integration into the host genome occurs only for tumourigenic strains of HPV. We hypothesized that changes in the epigenome and transcriptome contribute to the tumourigenicity of HPV. We found that viral integration events often occurred along with changes in chromatin state and expression of genes near the integration site. We investigated whether introduction of new transcription factor binding sites due to HPV integration could invoke these changes. Some regions within the HPV genome, particularly the position of a conserved CTCF binding site, showed enriched chromatin accessibility signal. ChIP-seq revealed that the conserved CTCF binding site within the HPV genome bound CTCF in 4 HPV+ cancer cell lines. Significant changes in CTCF binding pattern and increases in chromatin accessibility occurred exclusively within 100 kbp of HPV integration sites. The chromatin changes co-occurred with out-sized changes in transcription and alternative splicing of local genes. We analyzed the essentiality of genes upregulated around HPV integration sites of The Cancer Genome Atlas (TCGA) HPV+ tumours. HPV integration upregulated genes which had significantly higher essentiality scores compared to randomly selected upregulated genes from the same tumours. Our results suggest that introduction of a new CTCF binding site due to HPV integration reorganizes chromatin and upregulates genes essential for tumour viability in some HPV+ tumours. These findings emphasize a newly recognized role of HPV integration in oncogenesis. Creation of a single cell RNA seq meta-atlas to define human liver immune homeostasis The liver is unique in both its ability to maintain immune homeostasis and in its potential for immune tolerance following solid organ transplantation. Single-cell RNA sequencing (scRNA seq) is a powerful approach to generate highly dimensional transcriptome data to understand cellular phenotypes. However, when scRNA data is produced by different groups, with different data models, different standards, and samples processed in different ways, it can be challenging to draw meaningful conclusions from the aggregated data. The goal of this study was to establish a method to combine  human liver scRNA seq datasets by 1) characterizing the heterogeneity between studies and 2) using the meta-atlas to define the dominant phenotypes across immune cell subpopulations in healthy human liver. Publicly available scRNA seq data generated from liver samples obtained from a combined total of 17 patients and [~]32,000 cells were analyzed. Liver-specific immune cells (CD45+) were extracted from each dataset, and immune cell subpopulations (myeloid cells, NK and T cells, plasma cells, and B cells) were examined using dimensionality reduction (UMAP), differential gene expression, and ingenuity pathway analysis. All datasets co-clustered, but cell proportions differed between studies. Gene expression correlation demonstrated similarity across all studies, and canonical pathways that differed between datasets were related to cell stress and oxidative phosphorylation rather than immune-related function. Next, a meta-atlas was generated via data integration and compared against PBMC data to define gene signatures for each hepatic immune subpopulation. This analysis defined key features of hepatic immune homeostasis, with decreased expression across immunologic pathways and enhancement of pathways involved with cell death. This method for meta-analysis of scRNA seq data provides a novel approach to broadly define the features of human liver immune homeostasis. Specific pathways and cellular phenotypes described in this human liver immune meta-atlas provide a critical reference point for further study of immune mediated disease processes within the liver. 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. First evidence of in vitro cytotoxic effects of marine microlitter on Merluccius merluccius and Mullus barbatus, two Mediterranean commercial fish species Marine litter is composed mainly of plastics and is recognized as a serious threats to marine ecosystems. Ecotoxicological approaches have started elucidating the potential severity of microplastics (MPs) in controlled laboratory studies with pristine materials but no information exist on marine environmental microlitter as a whole. Here, we characterized the litter in the coastal Northern Tyrrhenian sea and in the stomach of two fish species of socio-economic importance, and exposed primary cell cultures of mucosal and lymphoid organs to marine microlitter for evaluating possible cytotoxic effects. An average of 0.30 {+/-} 0.02 microlitter items m-3 was found in water samples. FT-IR analysis revealed that plastic particles, namely HDPE, polyamide and polypropylene were present in 100% and 83.3% of Merluccius merluccius and Mullus barbatus analyzed, which overall ingested 14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual, respectively. Moreover, microlitter was confirmed as a vector of microorganisms. Lastly, the apical end-point of viability was found to be significantly reduced in splenic cells exposed in vitro to two microlitter conditions. Considering the role of the spleen in the mounting of adaptive immune responses, our results warrant more in-depth investigations for clarifying the actual susceptibility of these two species to anthropogenic microlitter.

HighlightsO_LI0.30 {+/-} 0.02 microlitter items m-3 were found at the surface of coastal Northern Tyrrhenian sea
C_LIO_LI14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual were retrieved from the stomach of hakes and mullets
C_LIO_LIThe ingested microlitter contained plastic items
C_LIO_LIMicrolitter was validated as a carrier of bacteria, fungi and flagellates
C_LIO_LISplenic cells exposed to two microlitter conditions for 72 hours suffered cytotoxicity
C_LI Intermittent compression increases flow through a gel: a mechanical model perhaps relevant to flow of cerebrospinal fluid Exchange of molecules between cerebrospinal fluid (CSF) and brain cells contributes to brain function and protection from dementia. Despite widespread acceptance of the  glymphatic theory, the route by which CSF is brought close enough to the neural tissue for solutes to be exchanged by extracellular diffusion is not entirely clear. Exogenous molecules injected into CSF are observed to reach the basement lamina that surrounds the dense capillary network. Transport of solutes by diffusion along the basement lamina, a gel of macromolcules about 100 nm thick, would be too slow; bulk flow in a static geometry would require unphysiologically high pressures. However, it is known that the pulsation of blood aids transport of CSF, and we hypothesized that this is because the pulsation intermittently squeezes the pericapillary lamina. In a primitive mimicry, we have tested whether intermittent squeezing increases flow through an agar gel. In all but one of 216 tests, pulsation caused a reversible increase, sometimes by a factor of 100 or more. The enhancement was greatest for frequencies 5-11 Hz and, over the tested range of pressure heads (20 - 50 cmH2O), was greatest for the lowest pressure. The results suggest a reason why exercise slows the aging of the brain. The TPLATE subunit is essential for structural assembly of the endocytic TSET complex All eukaryotic cells rely on endocytosis to regulate the plasma membrane proteome and lipidome. Most eukaryotic groups, with the exception of fungi and animals, have retained the evolutionary ancient TSET complex as a regulator of endocytosis. Despite the presence of similar building blocks in TSET, compared to other coatomer complexes, structural insight into this adaptor complex is lacking. Here, we elucidate the molecular architecture of the octameric plant TSET complex (TPLATE complex/TPC) using an integrative structural approach. This allowed us to describe a plant-specific connection between the TML subunit and the AtEH/Pan1 proteins and show a direct interaction between the complex and the plasma membrane without the need for any additional protein factors. Furthermore, we identify the appendage of TPLATE as crucial for complex assembly. Structural elucidation of this ancient adaptor complex vastly advances our functional as well as evolutionary insight into the process of endocytosis.

Graphical abstract O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible The members of the tribe Brassiceae share a whole genome triplication (WGT), and one proposed model for its formation is a "two-step" pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using POInT (the Polyploidy Orthology Inference Tool). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring prior to the arrival of the third. We also, for the 90,000 individual genes in our study, make parental "subgenome" assignments, inferring, with measured uncertainty, which of the progenitor genomes of the allohexaploidy each gene derives from. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, where subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty. The initial response of females towards congeneric males matches the propensity to hybridise in Ophthalmotilapia. AO_SCPLOWBSTRACTC_SCPLOWCichlid radiations often harbour closely related species with overlapping niches and distribution ranges. Such species sometimes hybridise in nature, which raises the question how can they coexist. This also holds for the Tanganyika mouthbrooders Ophthalmotilapia ventralis and O. nasuta. Earlier studies found indications of asymmetrical hybridisation with females of O. ventralis accepting males of O. nasuta, but not the other way around. We hypothesised that this was due to differences in the capacity for species recognition. Given the higher propensity of O. ventralis females towards hybridisation, we expect a reduced ability for species recognition in O. ventralis females, compared to O. nasuta females. We staged two experiments, one focusing on 22 female O. nasuta and one on 21 female O. ventralis. These fish were placed in one half of a tank and briefly exposed to a conspecific or a heterospecific male, a conspecific female, or nothing (control). Female response was evaluated by scoring six tracking parameters and by noting the occurrence of ten discrete behaviours before and during the encounter. Females always responded to the presence of another fish by approaching it. Remarkably, for both O. nasuta and O. ventralis, we did not find a different response between encounters with conspecific males and females. However, in agreement with our hypothesis, O. nasuta females behaved differently towards conspecific or heterospecific males, whereas O. ventralis females did not. When presented with a heterospecific male, O. nasuta females performed a lower number of  ram behaviours. Additionally, they never displayed the  flee behaviour, a component of the species mating repertoire that was seen in all but one of the presentations with a conspecific male. Our findings show that differences in species recognition at first encounter predict to a large degree the outcome of the mating process, even in the absence of mating behaviour. Deficiency of gluconeogenic enzyme PCK1 promotes murine NASH progression and fibrosis by activation PI3K/AKT/PDGF axis Background and AimsNonalcoholic steatohepatitis (NASH) is a chronic liver disease characterized by hepatic lipid accumulation, inflammation, and progressive fibrosis. However, the pathomechanisms underlying NASH are incompletely explored. Phosphoenolpyruvate carboxykinase 1 (PCK1) catalyzes the first rate-limiting step of gluconeogenesis. This study was designed to determine the role of PCK1 in regulating NASH progression.

MethodsLiver metabolism, hepatic steatosis, and fibrosis were evaluated at 24 weeks in liver-specific Pck1-knockout mice fed with NASH diet or chow diet. Gain- or loss-of-function approaches were used to explore the underlying mechanism in vitro. AKT and RhoA inhibitors were evaluated for NASH treatment in vivo.

ResultsHepatic PCK1 was downregulated in patients with NASH and mouse models of NASH. Mice with liver Pck1 deficiency displayed hepatic lipid disorder and liver injury fed with normal diet, while showed aggravated fibrosis and inflammation when fed NASH diet. Mechanistically, PCK1 deficiency upregulated genes involved in fatty acid transport and lipid droplet formation. Moreover, metabolomics analysis showed the accumulation of glycerol 3-phosphate, the substrate of triglyceride synthesis. Furthermore, the loss of hepatic PCK1 could activate the RhoA/PI3K/AKT pathway, which leads to increased secretion of PDGF-AA and promotes the activation of hepatic stellate cells. Accordingly, treatment with RhoA and AKT inhibitors alleviated NASH progression in the presence of Pck1 deletion in vivo.

ConclusionsPCK1 deficiency plays a key role in the development of hepatic steatosis and fibrosis by facilitating the RhoA/PI3K/AKT/PDGF-AA axis. These findings provide a novel insight into therapeutic approaches for the treatment of NASH. Evidence for a hydrogen sulfide-sensing E3 ligase in yeast In yeast, control of sulfur amino acid metabolism relies upon Met4, a transcription factor which activates the expression of a network of enzymes responsible for the biosynthesis of cysteine and methionine. In times of sulfur abundance, the activity of Met4 is repressed via ubiquitination by the SCFMet30 E3 ubiquitin ligase, but the mechanism by which the F-box protein Met30 senses sulfur status to tune its E3 ligase activity remains unresolved. Herein, we show that Met30 responds to flux through the transsulfuration pathway to regulate the MET gene transcriptional program. In particular, Met30 is responsive to the biological gas hydrogen sulfide, which is sufficient to induce ubiquitination of Met4 in vivo. Additionally, we identify important cysteine residues in Met30s WD-40 repeat region that sense the availability of sulfur in the cell. Our findings reveal how SCFMet30 dynamically senses the flow of sulfur metabolites through the transsulfuration pathway to regulate synthesis of these special amino acids. Identifying and correcting multiple sources of misspecification in GWAS summary statistics for polygenic scores Publicly available genome-wide association studies (GWAS) summary statistics exhibit uneven quality, which can impact the validity of follow-up analyses. First, we present an overview of possible misspecifications that come with GWAS summary statistics. Then, in both simulations and real data analyses, we show that additional information such as imputation INFO scores, allele frequencies, and per-variant sample sizes in GWAS summary statistics can be used to detect possible issues and correct for misspecifications in the GWAS summary statistics. One important motivation for us is to improve the predictive performance of polygenic scores built from these summary statistics. Unfortunately, due to the lack of reporting standards for GWAS summary statistics, this additional information is not systematically reported. We also show that using well-matched LD references can improve model fit and translate into more accurate prediction. Finally, we discuss how to make polygenic score methods such as lassosum and LDpred2 more robust to these misspecifications to improve their predictive power. Tourmaline: a containerized workflow for rapid and iterable amplicon sequence analysis using QIIME 2 and Snakemake BackgroundAmplicon sequencing (metabarcoding) is a common method to survey diversity of environmental communities whereby a single genetic locus is amplified and sequenced from the DNA of whole or partial organisms, organismal traces (e.g., skin, mucus, feces), or microbes in an environmental sample. Several software packages exist for analyzing amplicon data, among which QIIME 2 has emerged as a popular option because of its broad functionality, plugin architecture, provenance tracking, and interactive visualizations. However, each new analysis requires the user to keep track of input and output file names, parameters, and commands; this lack of automation and standardization is inefficient and creates barriers to meta-analysis and sharing of results.

FindingsWe developed Tourmaline, a Python-based workflow that implements QIIME 2 and is built using the Snakemake workflow management system. Starting from a configuration file that defines parameters and input files--a reference database, a sample metadata file, and a manifest or archive of FASTQ sequences--it uses QIIME 2 to run either the DADA2 or Deblur denoising algorithm, assigns taxonomy to the resulting representative sequences, performs analyses of taxonomic, alpha, and beta diversity, and generates an HTML report summarizing and linking to the output files. Features include support for multiple cores, automatic determination of trimming parameters using quality scores, representative sequence filtering (taxonomy, length, abundance, prevalence, or ID), support for multiple taxonomic classification and sequence alignment methods, outlier detection, and automated initialization of a new analysis using previous settings. The workflow runs natively on Linux and macOS or via a Docker container. We ran Tourmaline on a 16S rRNA amplicon dataset from Lake Erie surface water, showing its utility for parameter optimization and the ability to easily view interactive visualizations through the HTML report, QIIME 2 viewer, and R- and Python-based Jupyter notebooks.

ConclusionsAutomated workflows like Tourmaline enable rapid analysis of environmental and biomedical amplicon data, decreasing the time from data generation to actionable results. Tourmaline is available for download at github.com/aomlomics/tourmaline. Honeycomb crystallography: comb formation under geometric frustrations 1As honeybees build their nests in pre-existing tree cavities, they must deal with dealing with the presence of geometric constraints, resulting in non-regular hexagons and topological defects in the comb. In this work, we study how bees adapt to their environment in order to regulate the comb structure. Specifically, we identify the irregularities in honeycomb structure in the presence of various geometric frustrations. We 3D-print experimental frames with a variety of constraints imposed on the imprinted foundations. The combs constructed by the bees show clear evidence of reoccurring patterns built by bees in response to specific geometric frustrations on these starter frames. Furthermore, using an experimental-modeling framework, we demonstrate that these patterns can be successfully modeled and replicated through a simulated annealing process, in which the minimized potential is a variation of the Lennard-Jones potential that only considers first-neighbor interactions according to a Delaunay triangulation. Our simulation results not only confirm the connection between honeycomb structures and other crystal systems such as graphene, but also show that irregularities in the honeycomb structure can be explained as the result of local interactions between honeybees and their immediate surroundings, leading to emergent global order. Additionally, our computational model can be used to describe specific strategies that bees use to effectively solve each geometric mismatch problem while minimizing cost of comb building. 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.