Reticular Dysgenesis-associated Adenylate Kinase 2 deficiency causes failure of myelopoiesis through disordered purine metabolism Reticular Dysgenesis is a particularly grave form of severe combined immunodeficiency that affects the adaptive and innate immune system. Patients suffer from congenital neutropenia, lymphopenia, and deafness. The disease is caused by biallelic loss of function in mitochondrial Adenylate Kinase 2 (AK2). AK2 mediates the phosphorylation of AMP to ADP, as substrate for ATP synthesis. Accordingly, declining oxidative phosphorylation has been postulated as the driver of disease pathology. The mechanistic basis, however, remains incompletely understood. Single cell RNA-sequencing of patient bone marrow cells implicated altered RNA catabolism and ribonucleoprotein synthesis in the pathogenesis of Reticular Dysgenesis. To investigate these findings, we developed a disease model based on CRISPR-mediated disruption of the AK2 gene in primary human hematopoietic stem cells. We found that AK2-deficient myeloid progenitor cells not only have compromised mitochondrial energy metabolism and increased AMP levels, but also NAD+ and aspartate depletion, metabolites that rely on TCA-cycle activity for regeneration and synthesis. Furthermore, AK2-deficient cells exhibited strikingly increased levels of the purine nucleotide precursor IMP, decreased cellular RNA content, ribosome subunit expression, protein synthesis and a profoundly hypo-proliferative phenotype. The rise in IMP levels stemmed from increased AMP deamination. Pharmacologic inhibition of AMP deaminase normalized IMP levels in AK2-deficient cells, but further aggravated the disease phenotype, pointing to AMP catabolism as a possible metabolic adaptation to mitigate AMP-mediated toxicity. Inducing an adenosine disequilibrium in control cells produced a similar myeloid maturation defect.

This study shows that AK2 deficiency globally curtailed mitochondrial metabolism resulting in NAD+ and aspartate deficiency and disordered purine metabolism. AMP accumulation and its detrimental effects on ribonucleotide synthesis capacity may contribute to the failure of myelopoiesis in Reticular Dysgenesis.
The adenoviral E1B-55k protein present in HEK293 cells mediates abnormal accumulation of key WNT signaling proteins in large cytoplasmic aggregates. HEK293 cells are one of the most widely used cell lines in research and HEK293 cells are frequently used as an in vitro model for studying the WNT signaling pathway. The HEK293 cell line was originally established by transfection of human embryonic kidney cells with sheared adenovirus 5 DNA and it is known that that HEK293 cells stably express the adenoviral E1A and E1B-55k proteins. Here we show that HEK293 cells display an unexpected distribution of key components of the WNT/{beta}-catenin signaling pathway where AXIN1, APC, DVL2 and tankyrase are all co-localized in large spherical cytoplasmic aggregates. The cytoplasmic aggregates are enclosed by a narrow layer of the adenoviral E1B-55k protein. Reduction of E1B-55k protein levels leads to disappearance of the cytoplasmic aggregates thus corroborating an essential role of the E1B-55k protein in mediating the formation of the aggregates. Furthermore, HEK293 cells with reduced E1B-55k protein levels display reduced levels of transcriptional activation of WNT/{beta}-catenin signaling upon stimulation by the Wnt3A agonist. The demonstrated influence of the E1B-55k protein on the cellular localization of WNT/{beta}-catenin signaling components and on transcriptional regulation of WNT/{beta}-catenin signaling asks for caution in the interpretation of data derived from the HEK293 cell line.
Differences in gait stability and acceleration characteristics between healthy young and older females Our aim was to evaluate differences in gait acceleration intensity, variability and stability of feet and trunk between older females and young females using inertial sensors. Twenty older females (OF; mean age 68.4, SD 4.1 years) and eighteen young females (YF; mean age 22.3, SD 1.7 years) were asked to walk straight for 100 meters at their preferred speed, while wearing inertial sensors on heels and lower back. We calculated spatiotemporal measures, foot and trunk acceleration characteristics and their variability, as well as trunk stability using the local divergence exponent (LDE). Two-way analysis of variance (including the factors foot and age), Students t-test, and Mann-Whitney U test were used to compare statistical differences of measures between groups. Cohens d effects were calculated for each variable. Foot maximum vertical acceleration and amplitude, trunk-foot vertical acceleration attenuation, as well as their variability were significantly smaller in OF than in YF. In contrast, trunk mediolateral acceleration amplitude, maximum vertical acceleration, and amplitude, as well as their variability were significantly larger in OF than in YF. Moreover, OF showed lower stability (i.e., higher LDE values) in mediolateral acceleration, mediolateral and vertical angular velocity of the trunk. Even though we measured healthy older females, these participants showed lower vertical foot accelerations with higher vertical trunk acceleration, lower trunk-foot vertical acceleration attenuation, less gait stability, and more variability of the trunk, and hence, were more likely to fall. These findings suggest that instrumented gait measurements may help for early detection of changes or impairments in gait performance, even before this can be observed by clinical eye or gait speed. Fatal neuroinvasion and SARS-CoV-2 tropism in K18-hACE2 mice is partially independent on hACE2 expression Animal models recapitulating distinctive features of severe COVID-19 are critical to enhance our understanding of SARS-CoV-2 pathogenesis. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. The precise mechanisms of lethality in this mouse model remain unclear. Here, we evaluated the spatiotemporal dynamics of SARS-CoV-2 infection for up to 14 days post-infection. Despite infection and moderate pneumonia, rapid clinical decline or death of mice was invariably associated with viral neuroinvasion and direct neuronal injury (including brain and spinal neurons). Neuroinvasion was observed as early as 4 dpi, with virus initially restricted to the olfactory bulb supporting axonal transport via the olfactory neuroepithelium as the earliest portal of entry. No evidence of viremia was detected suggesting neuroinvasion occurs independently of entry across the blood brain barrier. SARS-CoV-2 tropism was not restricted to ACE2-expressing cells (e.g., AT1 pneumocytes), and some ACE2-positive lineages were not associated with the presence of viral antigen (e.g., bronchiolar epithelium and brain capillaries). Detectable ACE2 expression was not observed in neurons, supporting overexpression of ACE2 in the nasal passages and neuroepithelium as more likely determinants of neuroinvasion in the K18-hACE2 model. Although our work incites caution in the utility of the K18-hACE2 model to study global aspects of SARS-CoV-2 pathogenesis, it underscores this model as a unique platform for exploring the mechanisms of SARS-CoV-2 neuropathogenesis that may have clinical relevance acknowledging the growing body of evidence that suggests COVID-19 may result in long-standing neurologic consequences.

IMPORTANCECOVID-19 is predominantly a respiratory disease caused by SARS-CoV-2 that has infected more than 191 million people with over 4 million fatalities (2021-07-20). The development of animal models recapitulating distinctive features of severe COVID-19 is critical to enhancing our understanding of SARS-CoV-2 pathogenesis and in the evaluation of vaccine and therapeutic efficacy. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. Here, we show lethality of this model is invariably associated with viral neuroinvasion linked with viral replication and assembly. Importantly, pneumonia albeit invariably present was generally moderate with the absence of culturable infectious virus at peak neuroinvasion. The dynamics of viral neuroinvasion and pneumonia were only partially dependent on hACE2. Overall, this study provides an in-depth sequential characterization of the K18-hACE2 model following SARS-CoV-2 infection, highlighting its significance to further study the mechanisms of SARS-CoV-2 neuropathogenesis. E156G and Arg158, Phe-157/del mutation in NTD of spike protein in B.1.617.2 lineage of SARS-CoV-2 leads to immune evasion through antibody escape Emerging variants of SARS-CoV-2 with better immune escape mechanisms and higher transmissibility remains a persistent threat across the globe. B.1.617.2 (Delta) variant was first emerged from Maharashtra, India in December, 2020. This variant is classified to be a major cause and concern of the second wave of COVID-19 in India. In the present study, we explored the genomic and structural basis of this variant through computational analysis, protein modelling and molecular dynamics (MD) simulations approach. B.1.617.2 variant carried E156G and Arg158, Phe-157/del mutations in NTD of spike protein. These mutations in N-terminal domain (NTD) of spike protein of B.1.617.2 variant revealed more rigidity and reduced flexibility compared to spike protein of Wuhan isolate. Further, docking and MD simulation study with 4A8 monoclonal antibody which was reported to bind NTD of spike protein suggested reduced binding of B.1.617.2 spike protein compared to that of spike protein of Wuhan isolate. The results of the present study demonstrate the possible case of immune escape and thereby fitness advantage of the new variant and further warrants demonstration through experimental evidence. Our study identified the probable mechanism through which B.1.617.2 variant is more pathogenically evolved with higher transmissibility as compared to the wild-type.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=93 SRC="FIGDIR/small/447321v3_ufig1.gif" ALT="Figure 1">
View larger version (23K):
org.highwire.dtl.DTLVardef@b92cborg.highwire.dtl.DTLVardef@1d261f7org.highwire.dtl.DTLVardef@11da73eorg.highwire.dtl.DTLVardef@1cef6ca_HPS_FORMAT_FIGEXP  M_FIG C_FIG Community Evaluation of Glycoproteomics Informatics Solutions Reveals High-Performance Search Strategies of Serum N- and O-Glycopeptide Data Glycoproteome profiling (glycoproteomics) is a powerful yet analytically challenging research tool. The complex tandem mass spectra generated from glycopeptide mixtures require sophisticated analysis pipelines for structural determination. Diverse software aiding the process have appeared, but their relative performance remains untested. Conducted through the HUPO Human Proteome Project - Human Glycoproteomics Initiative, this community study, comprising both developers and users of glycoproteomics software, evaluates the performance of informatics solutions for system-wide glycopeptide analysis. Mass spectrometry-based glycoproteomics datasets from human serum were shared with all teams. The relative team performance for N- and O-glycopeptide data analysis was comprehensively established and validated through orthogonal performance tests. Excitingly, several high-performance glycoproteomics informatics solutions were identified. While the study illustrated that significant informatics challenges remain, as indicated by a high discordance between annotated glycopeptides, lists of high-confidence (consensus) glycopeptides were compiled from the standardised team reports. Deep analysis of the performance data revealed key performance-associated search variables and led to recommendations for improved "high coverage" and "high accuracy" glycoproteomics search strategies. This study concludes that diverse software for comprehensive glycopeptide data analysis exist, points to several high-performance search strategies, and specifies key variables that may guide future software developments and assist informatics decision-making in glycoproteomics. A New Paradigm for Investigating Real-World Social Behavior and its Neural Underpinnings Eye tracking and other behavioral measurements collected from patient-participants in their hospital rooms afford a unique opportunity to study immersive natural behavior for basic and clinical translational research. We describe an immersive social and behavioral paradigm implemented in patients undergoing evaluation for surgical treatment of epilepsy, with electrodes implanted in the brain to determine the source of their seizures. Our studies entail collecting eye tracking with other behavioral and psychophysiological measurements from patient-participants during unscripted behavior, including social interactions with clinical staff, friends and family, in the hospital room. This approach affords a unique opportunity to study the neurobiology of natural social behavior, despite the logistical, technical, and ethical challenges. Collecting neurophysiological data synchronized to behavioral and psychophysiological measures helps us to study the relationship between behavior and physiology. Combining across these rich data sources while participants eat, read, converse with friends and family, etc., enables clinical-translational research aimed at understanding the participants disorders and clinician-patient interactions, as well as basic research into natural, real-world behavior. We discuss data acquisition, quality control, annotation, and analysis pipelines that are required for our studies. We also discuss the clinical, logistical, and ethical and privacy considerations critical to working in the hospital setting. 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. Differential regulation of degradation and immune pathways underlies adaptation of the ectosymbiotic nematode Laxus oneistus to oxic-anoxic interfaces Eukaryotes may experience oxygen deprivation under both physiological and pathological conditions. Because oxygen shortage leads to a reduction in cellular energy production, all eukaryotes studied so far conserve energy by suppressing their metabolism. However, the molecular physiology of animals that naturally and repeatedly experience anoxia is underexplored. One such animal is the marine nematode Laxus oneistus. It thrives, invariably coated by its sulfur-oxidizing symbiont Candidatus Thiosymbion oneisti, in anoxic sulfidic or hypoxic sand. Here, transcriptomics and proteomics showed that, whether in anoxia or not, L. oneistus mostly expressed genes involved in ubiquitination, energy generation, oxidative stress response, immune response, development, and translation. Importantly, ubiquitination genes were also upregulated when the nematode was subjected to anoxic sulfidic conditions, together with genes involved in autophagy, detoxification and ribosome biogenesis. We hypothesize that these degradation pathways were induced to recycle damaged cellular components (mitochondria) and misfolded proteins into nutrients. Remarkably, when L. oneistus was subjected to anoxic sulfidic conditions, lectin and mucin genes were also upregulated, potentially to promote the attachment of its thiotrophic symbiont. Furthermore, the nematode appeared to survive oxygen deprivation by using an alternative electron carrier (rhodoquinone) and acceptor (fumarate), to rewire the electron transfer chain. On the other hand, under hypoxia, genes involved in costly processes (e.g., amino acid biosynthesis, development, feeding, mating) were upregulated, together with the worms Toll- like innate immunity pathway and several immune effectors (e.g., Bacterial Permeability Increasing proteins, fungicides).

In conclusion, we hypothesize that, in anoxic sulfidic sand, L. oneistus upregulates degradation processes, rewires oxidative phosphorylation and by reinforces its coat of bacterial sulfur-oxidizers. In upper sand layers, instead, it appears to produce broad-range antimicrobials and to exploit oxygen for biosynthesis and development. Using networks to analyze and visualize the distribution of overlapping reading frames in virus genomes Gene overlap occurs when two or more genes are encoded by the same nucleotides. This phenomenon is found in all taxonomic domains, but is particularly common in viruses, where it may increase the information content of compact genomes or influence the creation of new genes. Here we report a global comparative study of overlapping reading frames (OvRFs) of 12,609 virus reference genomes in the NCBI database. We retrieved metadata associated with all annotated reading frames in each genome record to calculate the number, length, and frameshift of OvRFs. Our results show that while the number of OvRFs increases with genome length, they tend to be shorter in longer genomes. The majority of overlaps involve +2 frameshifts, predominantly found in ds-DNA viruses. However, the longest overlaps involve no shift in reading frame (+0), increasing the selective burden of the same nucleotide positions within codons, instead of exposing additional sites to purifying selection. Next, we develop a new graph-based representation of the distribution of OvRFs among the reading frames of genomes in a given virus family. In the absence of an unambiguous partition of reading frames by homology at this taxonomic level, we used an alignment-free k-mer based approach to cluster protein coding sequences by similarity. We connect these clusters with two types of directed edges to indicate (1) that constituent reading frames are adjacent in one or more genomes, and (2) that the reading frames overlap. These adjacency graphs not only provide a natural visualization scheme, but also a novel statistical framework for analyzing the effects of gene- and genome-level attributes on the frequencies of overlaps. Ustilago maydis telomere protein Pot1 harbors an extra N-terminal OB fold and regulates homology-directed DNA repair factors in a dichotomous and context-dependent manner The telomere G-strand binding protein Pot1 plays multifaceted roles in telomere maintenance and protection. We examined the structure and activities of Pot1 in Ustilago maydis, a fungal model that recapitulates key features of mammalian telomere regulation. Compared to the well-characterized primate and fission yeast Pot1 orthologs, UmPot1 harbors an extra N-terminal OB-fold domain (OB-N), which was recently shown to be present in most metazoans. UmPot1 binds directly to Rad51 and regulates the latters strand exchange activity. Deleting the OB-N domain, which is implicated in Rad51-binding, caused telomere shortening, suggesting that Pot1-Rad51 interaction facilitates telomere replication. Depleting Pot1 through transcriptional repression triggered growth arrest as well as rampant recombination, leading to multiple telomere aberrations. In addition, telomere repeat RNAs transcribed from both the G- and C-strand were dramatically up-regulated, and this was accompanied by elevated levels of telomere RNA-DNA hybrids. Telomere abnormalities of pot1-deficient cells were suppressed, and cell viability was restored by the deletion of genes encoding Rad51 or Brh2 (the BRCA2 ortholog), indicating that homology-directed repair (HDR) proteins are key mediators of telomere aberrations and cellular toxicity. Together, these observations underscore the complex physical and functional interactions between Pot1 and DNA repair factors, leading to context-dependent and dichotomous effects of HDR proteins on telomere maintenance and protection. 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. Highly efficient synthetic CRISPR RNA/Cas9-based mutagenesis for rapid cardiovascular phenotypic screening in F0 zebrafish The zebrafish is a valuable vertebrate model to study cardiovascular formation and function due to the facile visualization and rapid development of the circulatory system in its externally growing embryos. Despite having distinct advantages, zebrafish have paralogs of many important genes, making reverse genetics approaches inefficient since generating animals bearing multiple gene mutations requires substantial efforts. Here, we present a simple and robust synthetic CRISPR RNA/Cas9-based mutagenesis approach for generating biallelic F0 zebrafish knockouts. Using a dual-guide synthetic CRISPR RNA/Cas9 ribonucleoprotein (dgRNP) system, we compared the efficiency of biallelic gene disruptions following the injections of one, two, and three dgRNPs per gene into the cytoplasm or yolk. We show that simultaneous cytoplasmic injections of three distinct dgRNPs per gene into one-cell stage embryos resulted in the most efficient and consistent biallelic gene disruptions. Importantly, this triple dgRNP approach enables efficient inactivation of cell autonomous and cell non-autonomous gene function, likely due to the low mosaicism of biallelic disruptions. In support of this finding, we provide evidence that the F0 animals generated by this method fully phenocopied the endothelial and peri-vascular defects observed in corresponding stable mutant homozygotes. Moreover, this approach faithfully recapitulated the trunk vessel phenotypes resulting from the genetic interaction between two vegfr2 zebrafish paralogs. Mechanistically, investigation of genome editing and mRNA decay indicates that the combined mutagenic actions of three dgRNPs per gene lead to an increased probability of frameshift mutations, enabling efficient biallelic gene disruptions. Therefore, our approach offers a highly robust genetic platform to quickly assess novel and redundant gene function in F0 zebrafish. Individual-based multiple-unit dissimilarity: novel indices and null model for assessing temporal variability in community composition Beta-diversity was originally defined spatially, i.e., as variation in community composition among sites in a region. However, the concept of beta-diversity has since been expanded to temporal contexts. This is referred to as "temporal beta-diversity", and most approaches are simply an extension of spatial beta-diversity. The persistence and turnover of individuals over time is a unique feature of temporal beta-diversity. Nakadai (2020) introduced the "individual-based beta-diversity" concept, and provided novel indices to evaluate individual turnover and compositional shift by comparing individual turnover between two periods at a given site. However, the proposed individual-based indices are applicable only to pairwise dissimilarity, not to multiple-temporal (or more generally, multiple-unit) dissimilarity. Here, individual-based beta-diversity indices are extended to multiple-unit cases. In addition, a novel type of random permutation criterion related to these multiple-unit indices for detecting patterns of individual persistence is introduced in the present study. To demonstrate the usage the properties of these indices compared to average pairwise measures, I applied them to a dataset for a permanent 50-ha forest dynamics plot on Barro Colorado Island in Panama. Information regarding "individuals" is generally missing from community ecology and biodiversity studies of temporal dynamics. In this context, the methods proposed here are expected to be useful for addressing a wide range of research questions regarding temporal changes in biodiversity, especially studies using traditional individual-tracked forest monitoring data. Volatile social environments can favour investments in quality over quantity of social relationships Cooperation does not occur in a vacuum: interactions develop over time in social groups that undergo demographic changes. Intuition suggests that stable social environments favour developing few but strong reciprocal relationships (a  focused strategy), while volatile social environments favour the opposite: more but weaker social relationships (a  diversifying strategy). We model reciprocal investments under a quality-quantity tradeoff for social relationships. We find that volatility, counterintuitively, can favour a focused strategy. This result becomes explicable through applying the theory of antagonistic pleiotropy, originally developed for senescence, to social life. Diversifying strategies show superior performance later in life, but with costs paid at young ages while the social network is slowly being built. Under volatile environments, many individuals die before reaching sufficiently old ages to reap the benefits. Social strategies that do well early in life are then favoured: a focused strategy leads individuals to form their first few social bonds quickly and to make strong use of existing bonds. Our model highlights the importance of pleiotropy and population age structure for the evolution of cooperative strategies and other social traits, and shows that it is not sufficient to reflect on the fate of survivors only, when evaluating the benefits of social strategies. Bulk-assembly of monodisperse coacervates and giant unilamellar vesicles with programmable hierarchical complexity Compartmentalization is an essential step for the emergence of protocellular life on the early earth and the endeavor to bottom-up construct synthetic cells in the lab. Among manifold natural strategies for compartmentalization, assembly of lipid vesicles or coacervate droplets with resemblance to cell envelope or cytoplasm have emerged as two dominant paradigms. However, the spontaneous assembly of lipid vesicles or coacervates generally results in polydisperse microcompartments and lacks the ability to efficiently integrate functional building blocks towards defined complexity, which give rise to the difficulty to achieve the order and complexity towards the scope of that of life via spontaneous assembly. Herein, we show that the interplay of coacervates and colloidal particles provides a bulk-assembly approach to form monodisperse coacervates and giant unilamellar vesicles (GUVs) with programmable hierarchical complexity. The hierarchical complexity of the synthetic cell microcompartments can be delicately engineered from the control of the spatial or temporal organization of components in entity unit to manifold forms of interconnected synthetic cell consortia, achieving a compartment hierarchy analogous to the organelle-cell-tissue structure. The increase in complexity gives rise to emergent properties, for instance, collective morphology evolution of coacervate assemblages, gated permeability of GUVs without sophisticated protein machinery, and remarkable structural and functional stability of GUVs at extreme conditions. This work paves an unprecedented step to form monodisperse, hierarchical coacervates and GUVs via bulk-assembly, which provides a prebiotically plausible approach for the proper integration of inanimate matter towards the order and complexity of life and monodisperse, hierarchical microdroplets for follow-up scientific and technological applications. Innovative GenExpA software for selecting suitable reference genes for reliable normalization of gene expression in melanoma The algorithms commonly used to select the best stable reference gene in RT-qPCR data analysis have their limitations. We showed that simple selection of the reference gene or pair of genes with the lowest stability value from the pool of potential reference genes - a commonly used approach - is not sufficient to accurately and reliably normalize the target gene transcript and can lead to biologically incorrect conclusions. For reliable assessment of changes in a target gene expression level, we propose our innovative GenExpA software, which works in a manner independent of the experimental model and the normalizer used. GenExpA software selects the best reference by combining the NormFinder algorithm with progressive removal of the least stable gene from the candidate genes in a given experimental model and in the set of daughter models assigned to it. The reliability of references is validated based on the consistency of the statistical analyses of normalized target gene expression levels through all models, described by the coherence score (CS). The use of the CS value imparts a new quality to qPCR analysis because it clarifies how low the stability value of reference must be in order for biologically correct conclusions to be drawn. We tested our method on qPCR data for the B4GALT genes family in melanoma, which is characterized by a high mutation rate, and in melanocytes. GenExpA is available at https://github.com/DorotaHojaLukowicz/GenExpA or https://www.sciencemarket.pl/baza-programow-open-source#oferty.

HighlightsGenExpA - next-generation software for normalizer selection and validation

The GenExpA tool defines how low the stability value of the reference should be

The GenExpA tool determines the level of gene expression analysis robustness

The GenExpA tool increases the speed of analysis and reduces cost Synthetic Cooling Agents in US-marketed E-cigarette Refill Liquids and Popular Disposable E-cigarettes: Chemical Analysis and Risk Assessment. BackgroundMenthol, through its cooling and pleasant sensory effects, facilitates smoking and tobacco product initiation, resulting in the high popularity of mint/menthol-flavored E-cigarettes. More recently, E-cigarette vendors started marketing synthetic cooling agents as additives that impart a cooling effect but lack a characteristic minty odor. Knowledge about content of synthetic coolants in US-marketed E-cigarette products and associated health risks is limited.

MethodsE-liquid vendor sites were searched with the terms "koolada", "kool/cool", "ice", or WS-3/WS-23, denoting individual cooling agents, and relevant refill E-liquids were purchased. "Ice" flavor varieties of Puffbar, the most popular disposable E-cigarette brand, were compared with non-"Ice" varieties. E-liquids were characterized, and synthetic coolants quantified using GC/MS. Margin of exposure (MOE), a risk assessment parameter, was calculated to assess the risk associated with synthetic coolant exposure from E-cigarette use.

ResultsWS-3 was detected in 24/25 refill E-liquids analyzed. All Puffbar flavor varieties contained either WS-23 (13/14) or WS-3 (5/14), in both "Ice"- and non-"Ice" flavors. Modeling consumption of WS-3 from vaped E-liquids, resulted in MOEs below the safe margin of 100 for most daily use scenarios. MOEs for WS-23 were <100 for 10/13 Puffbar flavors in all use scenarios. Puffbar power specifications are identical to Juul devices.

ConclusionsSynthetic cooling agents (WS-3/WS-23) were present in US-marketed E-cigarettes, at levels that may result in consumer exposures exceeding safety thresholds set by regulatory agencies. Synthetic coolants are not only found in mint-or menthol-flavored products, but also in fruit- and candy-flavored products, including popular disposable E-cigarette products such as Puffbar.

ImplicationsSynthetic cooling agents are widely used in "kool/cool"- and "ice"-flavored E-liquids and in E-liquids without these labels, both as a potential replacement for menthol or to add cooling  notes to non-menthol flavors. These agents may be used to bypass current and future regulatory limits on menthol content in tobacco products, and not just E-cigarettes. Since synthetic cooling agents are odorless, they may not fall under the category of "characterizing flavor", potentially circumventing regulatory measures based on this concept. Regulators need to consider the additional health risks associated with exposure to synthetic cooling agents. Royal knifefish generate powerful suction feeding through large neurocranial elevation and high epaxial muscle power Suction feeding in ray-finned fishes involves powerful buccal cavity expansion to accelerate water and food into the mouth. Previous XROMM studies in largemouth bass (Micropterus salmoides), bluegill sunfish (Lepomis macrochirus), and channel catfish (Ictalurus punctatus) have shown that more than 90% of suction power in high performance strikes comes from the axial musculature. Thus, the shape of the axial muscles and skeleton may impact suction feeding mechanics. Royal knifefish (Chitala blanci) have an unusual postcranial morphology, with a ventrally flexed vertebral column and relatively large mass of epaxial muscle. Based on their body shape, we hypothesized that royal knifefish would generate high power strikes by utilizing large neurocranial elevation, vertebral column extension, and epaxial shortening. As predicted, C. blanci generated high suction expansion power compared to the other three species studied to date (up to 160 W), which was achieved by increasing both the rate of volume change and the intraoral subambient pressure. The large epaxial muscle (25% of body mass) shortened at high velocities to produce large neurocranial elevation and vertebral extension (up to 41 deg, combined), as well as high muscle mass-specific power (up to 800 W kg-1). For the highest power strikes, axial muscles generated 95% of the power, and 64% of the axial muscle mass consisted of the epaxial muscles. The epaxial-dominated suction expansion of royal knifefish supports our hypothesis that postcranial morphology may be a strong predictor of suction feeding biomechanics.

SUMMARY STATEMENTRoyal knifefish rely on their distinct postcranial morphology--with a curved vertebral column and large dorsal body muscles--to produce large neurocranial elevation and powerful suction feeding. Dramatic consequences of reducing erythrocyte membrane cholesterol on Plasmodium falciparum Cholesterol is the most abundant lipid in the erythrocyte. During its blood stage development, the malaria parasite establishes an active cholesterol gradient across the various membrane systems within the infected erythrocyte. Interestingly, some antimalarial compounds have recently been shown to disrupt cholesterol homeostasis in intraerythrocytic stages of Plasmodium falciparum. These studies point to the importance of cholesterol for parasite growth. Previously, reduction of cholesterol from the erythrocyte membrane by treatment with methyl-{beta}-cyclodextrin (M{beta}CD) was shown to inhibit parasite invasion and growth. In addition, M{beta}CD treatment of trophozoite stage P. falciparum was shown to result in parasite expulsion from the host cell. We have revisited these phenomena by using live video microscopy, ultrastructural analysis, and response to antimalarial compounds. By using time-lapse video microscopy of fluorescently tagged parasites, we show that M{beta}CD treatment for just 30 min causes dramatic expulsion of the trophozoite stage parasites. This forceful expulsion occurs within 10 sec. Remarkably, the plasma membrane of the host cell from which the parasite has been expelled does not appear to be compromised. The parasitophorous vacuolar membrane (PVM) continued to surround the extruded parasite, but the PVM appeared damaged. Treatment with antimalarial compounds targeting PfATP4 or PfNCR1 prevented M{beta}CD-mediated extrusion of the parasites, pointing to a potential role of cholesterol dynamics underlying the expulsion phenomena. We also confirmed the essential role of erythrocyte plasma membrane cholesterol for invasion and growth of P. falciparum. This defect can be partially complemented by cholesterol and desmosterol but not with epicholesterol, revealing stereospecificity underlying cholesterol function. Overall, our studies advance previous observations and reveal unusual cell biological features underlying cholesterol depletion of infected erythrocyte plasma membrane.

ImportanceMalaria remains a major challenge in much of the world. Symptoms of malaria are caused by the growth of parasites belonging to Plasmodium spp. inside the red blood cells (RBC), leading to their destruction. The parasite depends upon its host for much of its nutritional need. Cholesterol is a major lipid in the RBC plasma membrane, which is the only source of this lipid for malarial parasites. We have previously shown that certain new antimalarial compounds disrupt cholesterol homeostasis in P. falciparum. Here we use live time-lapse video microscopy to show dramatic expulsion of the parasite from the host RBC when the cholesterol content of the RBC is reduced. Remarkably, this expulsion is inhibited by the antimalarials that disrupt lipid homeostasis. We also show stereospecificity of cholesterol in supporting parasite growth inside RBC. Overall, these results point to a critical role of cholesterol in physiology of malaria parasites. Choice of pre-processing pipeline influences clustering quality of scRNA-seq datasets BackgroundSingle-cell RNA sequencing (scRNA-seq) has quickly become one of the most dominant techniques in modern transcriptome assessment. In particular, 10X Genomics Chromium system, with its high throughput approach, turn key and thorough user guide made this cutting-edge technique accessible to many laboratories using diverse animal models. However, standard pre-processing, including the alignment and cell filtering pipelines might not be ideal for every organism or tissue. Here we applied an alternative strategy, based on the pseudoaligner kallisto, on twenty-two publicly available single cell sequencing datasets from a wide range of tissues of eight organisms and compared the results with the standard 10X Genomics Cell Ranger pipeline.

ResultsIn most of the tested samples, kallisto produced higher sequencing read alignment rates and total gene detection rates in comparison to Cell Ranger. Although datasets processed with Cell Ranger had higher cell counts, outside of human and mouse datasets, these additional cells were routinely of low quality, containing low gene detection rates. Thorough downstream analysis of one kallisto processed dataset, obtained from the zebrafish pineal gland, revealed clearer clustering, allowing the identification of an additional photoreceptor cell type that previously went undetected. The finding of the new cluster suggests that the photoreceptive pineal gland is essentially a bi-chromatic tissue containing both green and red cone-like photoreceptors and implies that the alignment and pre-processing pipeline can affect the discovery of biologically-relevant cell types.

ConclusionWhile Cell Ranger favors higher cell numbers, using kallisto results in datasets with higher median gene detection per cell. We could demonstrate that cell type identification was not hampered by the lower cell count, but in fact improved as a result of the high gene detection rate and the more stringent filtering. Depending on the acquired dataset, it can be beneficial to favor high quality cells and accept a lower cell count, leading to an improved classification of cell types. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Immunogenicity of SARS-CoV-2 trimetric spike protein associated to Poly(I:C) plus Alum The SARS-CoV-2 pandemic has had a social and economic impact worldwide, and vaccination is an efficient strategy for diminishing those damages. New adjuvant formulations are required for the high vaccine demands, especially adjuvant formulations that induce a Th1 phenotype. Herein we assess a vaccination strategy using a combination of Alum and polyinosinic:polycytidylic acid (Poly(I:C)) adjuvants plus the SARS-CoV-2 spike protein in a prefusion trimeric conformation by an intradermal (ID) route. We found high levels of IgG anti-spike antibodies in the serum by enzyme linked immunosorbent assay (ELISA) and high neutralizing titers against SARS-CoV-2 in vitro by neutralization assay, after one or two boosts. By evaluating the production of IgG subtypes, as expected, we found that formulations containing Poly(I:C) induced IgG2a whereas Alum did not. The combination of these two adjuvants induced high levels of both IgG1 and IgG2a. In addition, cellular immune responses of CD4+ and CD8+ T cells producing interferon-gamma were equivalent, demonstrating that the Alum + Poly(I:C) combination supported a Th1 profile. Based on the high neutralizing titers, we evaluated B cells in the germinal centers, which are specific for receptor-binding domain (RBD) and spike, and observed that more positive B cells were induced upon the Alum + Poly(I:C) combination. Moreover, these B cells produced antibodies against both RBD and non-RBD sites. We also studied the impact of this vaccination preparation (spike protein with Alum + Poly(I:C)) in the lungs of mice challenged with inactivated SARS-CoV-2 virus. We found a production of IgG, but not IgA, and a reduction in neutrophil recruitment in the bronchoalveolar lavage fluid (BALF) of mice, suggesting that our immunization scheme reduced lung inflammation. Altogether, our data suggest that Alum and Poly(I:C) together is a possible adjuvant combination for vaccines against SARS-CoV-2 by the intradermal route. The quantitative genetics of the prevalence of infectious diseases: hidden genetic variation due to Indirect Genetic Effects dominates heritable variation and response to selection Infectious diseases have profound effects on life, both in nature and agriculture. Despite the availability of well-established epidemiological theory, however, a quantitative genetic theory of the host population for the endemic prevalence of infectious diseases is almost entirely lacking. While several studies have demonstrated the relevance of the transmission dynamics of infectious diseases for heritable variation and response to selection of the host population, our current theoretical framework of quantitative genetics excludes these dynamics. As a consequence, we do not know which genetic effects of the host population determine the prevalence of an infection, and have no concepts of breeding value and heritable variation for endemic prevalence.

Here we integrate quantitative genetics and epidemiology, and propose a quantitative genetic theory for R0 and for the endemic prevalence of an infectious disease. We first identify the genetic factors that determine the prevalence of an infection, using an approach founded in epidemiological theory. Subsequently we investigate the population level consequences of individual genetic variation, both for R0 and for the endemic prevalence. Next, we present expressions for the breeding value and heritable variation, for both prevalence and individual binary disease status, and show that these parameters depend strongly on the level of the prevalence. Results show that heritable variation for endemic prevalence is substantially greater than currently believed, and increases when prevalence approaches zero, while heritability of individual disease status goes to zero. As a consequence, response of prevalence to selection accelerates considerably when prevalence goes down, in contrast to predictions from classical theory. Finally, we show that most of the heritable variation for the endemic prevalence of an infection is hidden due to indirect genetic effects, suggesting a key role for kin-group selection both in the evolutionary history of current populations and for genetic improvement strategies in animals and plants. Rattus Economicus: A neural model for rational regulation of water-seeking effort in rodents In the laboratory, animals motivation to work tends to be positively correlated with reward magnitude. But in nature, rewards earned by work are essential to survival (e.g., working to find water), and the payoff of that work can vary on long timescales (e.g., seasonally). Under these constraints, the strategy of working less when rewards are small could be fatal. We found that instead, rats in a closed economy did more work for water rewards when the rewards were stably smaller, a phenomenon also observed in human labor supply curves. Like human consumers, rats also showed elasticity of demand, consuming far more water per day when its price in effort was lower. The neural mechanisms underlying such "rational" market behaviors remain largely unexplored. We propose a dynamic utility maximization model that can account for the dependence of rat labor supply (trials/day) on the wage rate (ml/trial), and also predict the temporal dynamics of when rats work. Based on data from mice, we hypothesize that SFOGLUT neurons in lamina terminalis continuously compute the instantaneous marginal utility of voluntary work for water reward, and causally determine the amount and timing of work. Spatial organization of chromosomes leads to heterogeneous chromatin motion and drives the liquid- or gel-like dynamical behavior of chromatin Chromosome organization and dynamics are involved in regulating many fundamental processes such as gene transcription and DNA repair. Experiments unveiled that chromatin motion is highly heterogeneous inside cell nuclei, ranging from a liquid-like, mobile state to a gel-like, rigid regime. Using polymer modeling, we investigate how these different physical states and dynamical heterogeneities may emerge from the same structural mechanisms. We found that the formation of topologically-associating domains (TADs) is a key driver of chromatin motion heterogeneity. In particular, we demonstrated that the local degree of compaction of the TAD regulates the transition from a weakly compact, fluid state of chromatin to a more compact, gel state exhibiting anomalous diffusion and coherent motion. Our work provides a comprehensive study of chromosome dynamics and a unified view of chromatin motion enabling to interpret the wide variety of dynamical behaviors observed experimentally across different biological conditions, suggesting that the  liquid or  solid behaviour of chromatin are in fact two sides of the same coin.