A multiscale lipid and cellular atlas of the human kidney The human kidney is composed of over 26 cell types that actively coordinate with each other to form higher-order structures, such as the nephron. It is not yet understood how these structures vary throughout a single organ or amongst the same organs within the human population. We have developed an extensive lipid and cellular atlas of the human kidney consisting of over 3 million cells comprising 75,000 functional tissue units (i.e., glomeruli, proximal tubules, distal tubules, and collecting ducts) from 13 human subjects. This atlas was developed using spatially registered and integrated technologies consisting of imaging mass spectrometry, multiplexed immunofluorescence, stained microscopy, and autofluorescence microscopy, to comprehensively probe large (i.e., centimeter-sized) areas of tissue. The cellular organization and lipid profiles of glomeruli, proximal tubules, distal tubules, and collecting ducts were discovered through these multimodal imaging data as well as their intra- and inter-subject variance. Relating the lipid profiles obtained from imaging mass spectrometry to distinct cell types obtained from immunofluorescence allowed us to hypothesize the functional role of specific phospholipids that have not previously been described. These hypotheses include subject characteristics, such as BMI and sex. The integrated data from the aforementioned datasets provide a valuable reference for kidney researchers, are publicly available through the NIH Human Biomolecular Atlas Program (https://portal.hubmapconsortium.org/), and discussed below.
Pore dynamics and asymmetric cargo loading in an encapsulin nanocompartment Encapsulins are protein nanocompartments that house various cargo enzymes, including a family of decameric ferritin-like proteins. Here, we study a recombinant Haliangium ochraceum encapsulin:encapsulated ferritin complex using electron cryo-microscopy and hydrogen/deuterium exchange mass spectrometry to gain insight into the structural relationship between the encapsulin shell and its protein cargo. An asymmetric single particle reconstruction reveals four encapsulated ferritin decamers in a tetrahedral arrangement within the encapsulin nanocompartment. This leads to a symmetry mismatch between the protein cargo and the icosahedral encapsulin shell. The encapsulated ferritin decamers are offset from the interior face of the encapsulin shell. Using HDX-MS, we observed dynamic behavior of the major five-fold pore in the encapsulin shell and show the pore opening via the movement of the encapsulin A-domain. These data will accelerate efforts to engineer the encapsulation of heterologous cargo proteins and to alter the permeability of the encapsulin shell via pore modifications.

TeaserCryo-EM and HDX-MS analysis of an encapsulin nanocompartment shows that the pores at the five-fold icosahedral vertex of the shell are flexible.
Heritable changes in division speed accompany the diversification of singleT cell fate Rapid clonal expansion of antigen specific T cells is a fundamental feature of adaptive immune responses. It enables the outgrowth of an individual T cell into thousands of clonal descendants that diversify into short-lived effectors and long-lived memory cells. Clonal expansion is thought to be programmed upon priming of a single naive T cell and then executed by homogenously fast divisions of all of its descendants. However, the actual speed of cell divisions in such an emerging  T cell family has never been measured with single-cell resolution. Here, we utilize continuous live-cell imaging in vitro to track the division speed and genealogical connections of all descendants derived from a single naive CD8+ T cell throughout up to ten divisions of activation-induced proliferation. This comprehensive mapping of T cell family trees identifies a short burst phase, in which division speed is homogenously fast and maintained independent of external cytokine availability or continued T cell receptor stimulation. Thereafter, however, division speed diversifies and model-based computational analysis using a novel Bayesian inference framework for tree-structured data reveals a segregation into heritably fast and slow dividing branches. This diversification of division speed is preceded already during the burst phase by variable expression of the interleukin-2 receptor alpha chain. Later it is accompanied by selective expression of memory marker CD62L in slower dividing branches. Taken together, these data demonstrate that T cell clonal expansion is structured into subsequent burst and diversification phases the latter of which coincides with specification of memory vs. effector fate.

SignificanceRapid clonal expansion of antigen-specific T cells is a fundamental feature of adaptive immune responses. Here, we utilize continuous live-cell imaging in vitro to track the division speed and genealogical connections of all descendants derived from a single naive CD8+ T cell throughout up to ten divisions of activation-induced proliferation. Bayesian inference of tree-structured data reveals that clonal expansion is divided into a homogenously fast burst phase encompassing two to three divisions and a subsequent diversification phase during which T cells segregate into quickly dividing effector T cells and more slowly cycling memory precursors. Our work highlights cell cycle speed as a major heritable property that is regulated in parallel to key lineage decisions of activated T cells. Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. Using total RNA quality metrics for time since deposition estimates in degrading bloodstains Determining the time since deposition (TSD) of bloodstains would provide forensic scientists with critical information regarding the timeline of the events involving bloodshed. The physicochemical changes occurring to biomolecules as a bloodstain degrades can be used to approximate the TSD of bloodstains. Our study aims to quantify the timewise degradation trends and temperature dependence found in total RNA from bloodstains without the use of amplification, expanding the scope of the RNA TSD research which has previously targeted mRNA molecules. Whole bovine blood was stored in plastic microcentrifuge tubes at 21{degrees}C or 4{degrees}C and tested over different timepoints spanning one week. Nine RNA metrics were visually assessed and quantified using linear and mixed models; the RNA Integrity Number equivalent (RINe) and the DV200 demonstrated strong negative trends over time and statistical independence. The RINe model fit was high (R2 = 0.60), and while including the biological replicate as a random effect increased the fit for all RNA metrics, no significant differences were found between biological replicates stored at the same temperature for the RINe and DV200 metrics. Importantly, this suggests that these standardized metrics can likely be directly compared between scenarios and individuals, with DV200 having an inflection point at [~]28 hrs. This study provides a novel approach for blood TSD estimates, producing metrics that are not affected by inter-individual variation and improving our understanding of the rapid degradation occurring in bloodstains.

HIGHLIGHTSO_LIAmplification-free analysis of total RNA in degrading bloodstains.
C_LIO_LIShort-term RNA degradation assessment using high-resolution size measurements.
C_LIO_LITotal RNA quality and quantity metrics were assessed across a one-week.
C_LIO_LITotal RNA quality metrics demonstrated the strongest timewise trends.
C_LIO_LIBiological replicates produced similar results for RNA quality metrics.
C_LI A distinct brain expression profile for genes within large introgression deserts and under positive selection in Homo sapiens Analyses of ancient DNA from extinct hominins have provided unique insights into the complex evolutionary history of Homo sapiens, intricately related to that of the Neanderthals and the Denisovans as revealed by several instances of admixture events. These analyses have also allowed the identification of introgression deserts: genomic regions in our species that are depleted of  archaic haplotypes. The presence of genes like FOXP2 in these deserts has been taken to be suggestive of brain-related functional differences between Homo species. Here, we seek a deeper characterization of these regions, taking into account signals of positive selection in our lineage. Analyzing publicly available transcriptomic data from the human brain at different developmental stages, we found that structures outside the cerebral neocortex, in particular the cerebellum, the striatum and the mediodorsal nucleus of the thalamus show the most divergent transcriptomic profiles when considering genes within large introgression deserts and under positive selection. TXN, a Xanthohumol Derivative, Attenuates High-Fat Diet Induced Hepatic Steatosis by Antagonizing PPARγ We previously reported xanthohumol (XN), and its synthetic derivative tetrahydro-XN (TXN) attenuates high-fat diet (HFD) induced obesity and metabolic syndrome in C57BL/6J mice. The objective of the current study was to determine the effect of XN and TXN on lipid accumulation in the liver. Non-supplemented mice were unable to adapt their caloric intake to 60% HFD, resulting in obesity and hepatic steatosis; however, TXN reduced weight gain and decreased hepatic steatosis. Liver transcriptomics indicated TXN might antagonize lipogenic PPAR{gamma} actions in vivo. XN and TXN inhibited rosiglitazone-induced 3T3-L1 cell differentiation concomitant with decreased expression of lipogenesis-related genes. A PPAR{gamma} competitive binding assay showed XN and TXN bind to PPAR{gamma} with an IC50 similar to pioglitazone and 8-10 times stronger than oleate. Molecular docking simulations demonstrated XN and TXN bind in the PPAR{gamma} ligand-binding domain pocket. Our findings are consistent with XN and TXN acting as antagonists of PPAR{gamma}. An adaptive teosinte mexicana introgression modulates phosphatidylcholine levels and is associated with maize flowering time Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp.parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability, and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern US, Canada and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time. Chromosome breaks in breast cancers occur near herpes tumor virus sequences and explain why cancer drugs stop working One explanation for why breast cancer has a relentless tendency to come back after treatment is that breast cancer phenomes are fluid, so it is difficult for therapy to target new phenotypes that arise. In breast cancer, viral interactions cause chromosome breakage and fusion, producing anaphase bridges containing two centromeres. Cell division then leads to more chromosome breaks, more mutations, and new cancer phenotypes. Analyses of public data from about 2100 breast cancers implicate Epstein-Barr virus (EBV or Human Herpes Virus 4) in the cycles that break chromosomes and produce anaphase bridges. EBV generates breakage cycles with highly fragmented chromosomes (chromothripsis) as an inherent feature. These episodes continually produce new cancer driver mutations. DNA at chromothripsis boundaries in breast cancers consistently matches EBV tumor variants. Breakpoints are often near known EBV binding sites, and breast cancer shares an EBV methylation signature with other known EBV cancers. Breast cancer breakpoints on all chromosomes cluster around the same positions as in cancers definitively associated with EBV infection (nasopharyngeal cancer and Burkitts lymphoma). The status of EBV in the microbiome and the presence of anaphase bridges are phenotypic markers of breast cancer that are not explicitly targeted by chemotherapy. The breast cancer phenome continually supplies new cancer driver genes that foster breast cancer recurrence.

SummaryBreaks in breast cancer chromosomes and in fragments of anaphase bridges occur mainly near regions that resemble Epstein-Barr virus (EBV) tumor variants

Extensive breaks are an inherent outcome of anaphase bridges and the reassembly of highly fragmented chromosomes continually generates new cancer driver mutations and new cancer phenotypes.

Treatment does not explicitly target the underlying breakage and reassembly cycles so breast cancer often comes back.

The breast cancer phenome is not fixed and can continuously generate new phenotypes

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=154 HEIGHT=200 SRC="FIGDIR/small/467751v4_ufig1.gif" ALT="Figure 1">
View larger version (38K):
org.highwire.dtl.DTLVardef@6537cforg.highwire.dtl.DTLVardef@132fa50org.highwire.dtl.DTLVardef@1cfcb04org.highwire.dtl.DTLVardef@11979aa_HPS_FORMAT_FIGEXP  M_FIG C_FIG Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. 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. Transcriptomic analysis of temporal shifts in berry development between two grapevine cultivars of the Pinot family reveals potential genes controlling ripening time BackgroundGrapevine cultivars of the Pinot family represent clonally propagated mutants with major phenotypic and physiological differences, such as different colour or shifted ripening time, as well as changes in important viticultural traits. Specifically, the cultivars  Pinot Noir (PN) and  Pinot Noir Precoce (PNP, early ripening) flower at the same time, but vary in the beginning of berry ripening (veraison) and, consequently, harvest time. In addition to genotype, seasonal climatic conditions (i.e. high temperatures) also affect ripening times. To reveal possible regulatory genes that affect the timing of veraison onset, we investigated differences in gene expression profiles between PN and PNP throughout berry development with a closely meshed time series and over two separate years.

ResultsThe difference in the duration of berry formation between PN and PNP was quantified to be approximately two weeks under the growth conditions applied, using plant material with a proven PN and PNP clonal relationship. Clusters of co-expressed genes and differentially expressed genes (DEGs) were detected which reflect the shift in the timing of veraison onset. Functional annotation of these DEGs fit to observed phenotypic and physiological changes during berry development. In total, we observed 3,342 DEGs in 2014 and 2,745 DEGs in 2017 between PN and PNP, with 1,923 DEGs across both years. Among these, 388 DEGs were identified as veraison-specific and 12 were considered as berry ripening time regulatory candidates. The expression profiles revealed two candidate genes for ripening time control which we designated VviRTIC1 and VviRTIC2 (VIT_210s0071g01145 and VIT_200s0366g00020, respectively). These genes likely contribute the phenotypic differences observed between PN and PNP.

ConclusionsMany of the 1,923 DEGs show highly similar expression profiles in both cultivars if the patterns are aligned according to developmental stage. In our work, putative genes differentially expressed between PNP and PN which could control ripening time as well as veraison-specific genes were identified. We point out connections of these genes to molecular events during berry development and discuss potential candidate genes which may control ripening time. Two of these candidates were observed to be differentially expressed in the early berry development phase. Several down-regulated genes during berry ripening are annotated as auxin response factors / ARFs. Conceivably, general changes in auxin signaling may cause the earlier ripening phenotype of PNP. Substrate specialization in microbes is driven by biochemical constraints of dynamic flux sensing Central carbon metabolism is highly conserved across microbial species, but can catalyze very different pathways depending on the organism and their ecological niche. Here, we study the dynamic re-organization of central metabolism after switches between the two major opposing pathway configurations of central carbon metabolism, glycolysis and gluconeogenesis in Escherichia coli, Pseudomonas aeruginosa and Pseudomonas putida. We combined growth dynamics and dynamic changes of intracellular metabolite levels with a coarse-grained model that integrates fluxes, regulation, protein synthesis and growth and uncovered fundamental limitations of the regulatory network: after nutrient shifts, metabolite concentrations collapse to their equilibrium, rendering the cell unable to sense which direction the flux is supposed to flow through the metabolic network. The cell can partially alleviate this by picking a preferred direction of regulation at the expense of increasing lag times in the opposite direction. Moreover, decreasing both lag times simultaneously comes at the cost of reduced growth rate or higher futile cycling between metabolic enzymes. These three trade-offs can explain why microorganisms specialize for either glycolytic or gluconeogenic substrates and can help elucidate the complex growth patterns exhibited by different microbial species.

Graphical synopsis O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

Standfirst textMicrobes face a series of fundamental trade-offs that limit their ability to optimize simultaneously for both glycolytic and gluconeogenic growth.

Bullet pointsO_LILag times between glycolysis and gluconeogenesis show asymmetry in many microbes: A long lag in one direction, but a short lag in the other.
C_LIO_LILong lag times are caused by an inability to sense fluxes after nutrient shifts.
C_LIO_LIWith existing regulation, lag time asymmetry can only be overcome by reducing either growth rate or increasing futile cycling in metabolism.
C_LI Single-Cell Multiomics Reveals Distinct Cell States at the Top of the Human Hematopoietic Hierarchy The advent of single cell (Sc) genomics has challenged the dogma of haematopoiesis as a tree-like structure of stepwise lineage commitment through distinct and increasingly restricted progenitor populations. Instead, analysis of ScRNA-seq has proposed that the earliest events in human hematopoietic stem cell (HSC) differentiation are characterized by only subtle molecular changes, with hematopoietic stem and progenitor cells (HSPCs) existing as a continuum of low-primed cell-states that gradually transition into a specific lineage (CLOUD-HSPCs). Here, we combine ScRNA-seq, ScATAC-seq and cell surface proteomics to dissect the heterogeneity of CLOUD-HSPCs at different stages of human life. Within CLOUD-HSPCs, pseudotime ordering of both mRNA and chromatin data revealed a bifurcation of megakaryocyte/erythroid and lympho/myeloid trajectories immediately downstream a subpopulation with an HSC-specific enhancer signature. Importantly, both HSCs and lineage-restricted progenitor populations could be prospectively isolated based on correlation of their molecular signatures with CD35 and CD11A expression, respectively. Moreover, we describe the changes that occur in this heterogeneity as hematopoiesis develops from neonatal to aged bone marrow, including an increase of HSCs and depletion of lympho-myeloid biased MPPs. Thus, this study dissects the heterogeneity of human CLOUD-HSPCs revealing distinct HSPC-states of relevance in homeostatic settings such as ageing. 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. A multi-state dynamic process confers mechano-adaptation to the bacterial flagellar motor Adaptation is a defining feature of living systems. The bacterial flagellar motor adapts to changes in the external mechanical environment by adding or removing torque-generating stator units. However, the molecular mechanism for mechanosensitive motor remodeling remains unclear. Here, we induced stator disassembly using electrorotation, followed by the time-dependent assembly of the individual stator units into the motor after electrorotation was terminated. From these experiments, we extracted detailed statistics of the dwell times that comprise the stochastic dynamics of the binding and unbinding of stator units. Dwell times reveal multiple timescales, indicating the existence of multiple binding states of the stator units. Based on these results, we propose a minimal model in which the stator unit can occupy four different states - two bound states with very different rates of unbinding, a diffusive unbound state, and a transiently detached state. Our minimal model quantitatively explains multiple features of the experimental data and allows us to determine the transition rates among all four states. Our experiments and modeling suggest a mechanism of mechano-adaptive remodeling of the bacterial flagellar motor in which torque generated by bound stator units controls their effective unbinding rate by modulating the transition between the two bound states. Furthermore, the binding rate of stator units with the motor has a non-monotonic dependence on the number of bound units, likely because of two counteracting effects of motor rotation on the binding process. Improving human coronavirus OC43 (HCoV-OC43) research results comparability in studies using HCoV-OC43 as a surrogate for SARS-CoV-2 The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has renewed interest in human coronaviruses that cause the common cold, particularly as research with them at biosafety level (BSL)-2 avoids the added costs and biosafety concerns that accompany work with SARS-COV-2, BSL-3 research. One of these, human coronavirus OC43 (HCoV-OC43), is a well-matched surrogate for SARS-CoV-2 because it is also a Betacoronavirus, targets the human respiratory system, is transmitted via respiratory aerosols and droplets and is relatively resistant to disinfectants. Unfortunately, growth of HCoV-OC43 in the recommended human colon cancer (HRT-18) cells does not produce obvious cytopathic effect (CPE) and its titration in these cells requires expensive antibody-based detection. Consequently, multiple quantification approaches for HCoV-OC43 using alternative cell lines exist, which complicates comparison of research results. Hence, we investigated the basic growth parameters of HCoV-OC43 infection in three of these cell lines (HRT-18, human lung fibroblasts (MRC-5) and African green monkey kidney (Vero E6) cells) including the differential development of cytopathic effect (CPE) and explored reducing the cost, time and complexity of antibody-based detection assay. Multi-step growth curves were conducted in each cell type in triplicate at a multiplicity of infection of 0.1 with daily sampling for seven days. Samples were quantified by tissue culture infectious dose50(TCID50)/ml or plaque assay (cell line dependent) and additionally analyzed on the Sartorius Virus Counter 3100 (VC), which uses flow virometry to count the total number of intact virus particles in a sample. We improved the reproducibility of a previously described antibody-based detection based TCID50 assay by identifying commercial sources for antibodies, decreasing antibody concentrations and simplifying the detection process. The growth curves demonstrated that HCoV-O43 grown in MRC-5 cells reached a peak titer of [~]107 plaque forming units/ml at two days post infection (dpi). In contrast, HCoV-OC43 grown on HRT-18 cells required six days to reach a peak titer of [~]106.5 TCID50/ml. HCoV-OC43 produced CPE in Vero E6 cells but these growth curve samples failed to produce CPE in a plaque assay after four days. Analysis of the VC data in combination with plaque and TCID50 assays together revealed that the defective:infectious virion ratio of MRC-5 propagated HCoV-OC43 was less than 3:1 for 1-6 dpi while HCoV-OC43 propagated in HRT-18 cells varied from 41:1 at 1 dpi, to 329:4 at 4 dpi to 94:1 at 7 dpi. These results should enable better comparison of extant HCoV-OC43 study results and prompt further standardization efforts. Structural-functional connectivity mapping of the insular cortex: A combined data-driven and meta-analytic topic mapping In this study, we examined structural and functional profiles of the insular cortex and mapped associations with well-described functional networks (FNs) throughout the brain using diffusion tensor imaging (DTI) and resting-state functional connectivity (RSFC) data. We used a data-driven method to independently estimate the structural-functional connectivity of the insular cortex. Data were obtained from the Human Connectome Project comprising 108 adult participants. Overall, we observed moderate to high associations between the structural and functional mapping scores of three different insular subregions: the posterior insula (associated with the sensorimotor network: RSFC, DTI = 50% and 72%, respectively), dorsal anterior insula (associated with ventral attention: RSFC, DTI = 83% and 83%, respectively), and ventral anterior insula (associated with the frontoparietal: RSFC, DTI = 42% and 89%, respectively). Further analyses utilized meta-analytic decoding maps to demonstrate specific cognitive and affective as well as gene expression profiles of the three subregions reflecting the core properties of the insular cortex. In summary, given the central role of the insular in the human brain, our results revealing correspondence between DTI and RSFC mappings provide a complementary approach and insight for clinical researchers to identify dysfunctional brain organization in various neurological disorders associated with insular pathology. Outsourced hearing in an orb-weaving spider that uses its web as an auditory sensor Hearing is a fundamental sense of many animals, including all mammals, birds, some reptiles, amphibians, fish, and arthropods. The auditory organs of these animals are extremely diverse in anatomy after hundreds of millions of years of evolution, yet all are made up of cellular tissues and are morphologically part of bodies of animals. Here we show hearing in the orb-weaving spider, Larinioides sclopetarius is not constrained by the organisms body but is extended through outsourcing hearing to its extended phenotype, the proteinaceous, self-manufactured orb-web. We find the wispy, wheel-shaped orb-web acts as a hyperacute acoustic "antenna" to capture the sound-induced air particle movements that approach the maximum physical efficiency, better than the acoustic responsivity of all previously known eardrums. By sensing the motion of web threads, the spider remotely detects and localizes the source of an incoming airborne acoustic wave such as those emitted by approaching prey or predators. By outsourcing its acoustic sensors to its web, the spider is released from body size constraints and permits the araneid spider to increase its sound-sensitive surface area enormously, up to 10,000 times greater than the spider itself. The spider also enables the flexibility to functionally adjust and regularly regenerate its "external eardrum" according to its needs. The "outsourcing" and "supersizing" of auditory function in spiders provides unique features for studying extended and regenerative sensing, and designing novel acoustic flow detectors for precise fluid dynamic measurement and manipulation. Land use change through the lens of macroecology: insights from Azorean arthropods and the Maximum Entropy Theory of Ecology Human activity and land management practices, in particular land use change, have resulted in the global loss of biodiversity. These types of disturbance affect the shape of macroecological patterns, and therefore analyzing these patterns can provide insights into how ecosystems are affected by land use change. We here use arthropod census data from 96 sites at Terceira Island in the Azores archipelago across four different land uses of increasing management intensity: native forest, exotic forest, semi-natural pasture, and intensive pasture, to examine the effects of land use type on three macroecological patterns: the species abundance distribution, the metabolic rate distribution of individuals, and the species-area relationship. The Maximum Entropy Theory of Ecology (METE) has successfully predicted these patterns across habitats and taxa in undisturbed ecosystems, and thus provides a null expectation for their shapes. Across these patterns, we find that the forest habitats are the best fit by METE predictions, while the semi-natural pasture is consistently the worst fit, and the intensive pasture is intermediately well fit. We show that the direction of failure of the METE predictions at the pasture sites is likely due to the hyper-dominance of introduced spider species present there. We hypothesize that the particularly poor fit for the semi-natural pasture is due to the mix of arthropod communities out of equilibrium, leading to greater heterogeneity in composition and complex dynamics that violate METEs assumption of static state variables. The comparative better fit for the intensive pasture plausibly results from more homogeneous arthropod communities that are well adapted to intensive management, and thus whose state variables are less in flux. Analyzing deviations from theoretical predictions across land use type provides useful information about how land use and disturbance affect ecosystems, and such comparisons could be useful across other habitats and taxa. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. Embedded 3D printing in self-healing annealable composites for precise patterning of functionally mature human neural constructs Human in vitro models of neural tissue with controllable cellular identity, tunable microenvironment, and defined spatial arrangement are needed to facilitate studies of brain development and disease. Towards this end, embedded printing in jammed microgel supports (i.e., granular gels) holds great promise as it allows precise and programmable patterning of extremely soft and compliant tissue constructs. However, in contrast to the vast material landscape available for bulk hydrogels, granular printing support formulations are restricted to a handful of materials without the ability for facile adjustment of biofunctional properties of the cellular microenvironment. Therefore, there has been a need for novel materials that take advantage of versatile biomimicry of bulk hydrogels while providing high-fidelity support for embedded printing akin to granular gels. To address this need, we present a modular platform for bioengineering of neuronal networks via direct embedded 3D printing of human stem cells inside Self-Healing Annealable Particle-Extracellular matrix (SHAPE) composites. SHAPE composites consist of soft microgels immersed in viscous extracellular-matrix solution to enable precise freeform patterning of human stem cells and consequent generation and long-term maintenance of mature subtype-specific neurons that extend projections within the volume of the annealed support. The developed approach further allows multi-ink deposition, live spatial and temporal monitoring of oxygen levels, as well as creation of vascular channels. Due to its modularity, SHAPE biomanufacturing toolbox not only offers a solution for functional modeling of mechanically sensitive neural constructs, but also has potential to be applied to a wide range of biomaterials with different crosslinking mechanisms to model tissues and diseases where recapitulation of complex architectural features and topological cues is essential. An integrative approach leads to the discovery of a novel anti-leukemic peptide from human milk Chemotherapy in childhood leukemia is associated with late morbidity in leukemic survivors, while certain patient subsets are relatively resistant to standard chemotherapy. It is therefore important to identify new agents with sensitivity and selectivity towards leukemic cells, while having less systemic toxicity. Peptide-based therapeutics has gained much attention during the last few years. Here, we used an integrative workflow combining mass spectrometric peptide library construction, in silico anticancer peptide screening, and in vitro leukemic cell studies to discover a novel anti-leukemic peptide having 3+charges and alpha-helical structure, namely HMP-S7, from human breast milk. HMP-S7 showed cytotoxic activity against four distinct leukemic cell lines in a dose-dependent manner but had no effect on solid malignancies or representative normal cells. HMP-S7 induced leukemic cell death by penetrating the plasma membrane to enter the cytoplasm and cause leakage of lactate dehydrogenase, thus acting in a membranolytic manner. Importantly, HMP-S7 exhibited anti-leukemic effect against patient-derived leukemic cells ex vivo. In conclusion, HMP-S7 is a selective anti-leukemic peptide with promise which requires further validation in preclinical and clinical studies.

TeaserIn silico screening of naturally occurring human milk peptides discovers a new anticancer peptide that kills leukemic cells in vitro and ex vivo. Designing viral diagnostics with model-based optimization Diagnostics, particularly for rapidly evolving viruses, stand to benefit from a principled, measurement-driven design that harnesses machine learning and vast genomic data--yet the capability for such design has not been previously built. Here, we develop and extensively validate an approach to designing viral diagnostics that applies a learned model within a combinatorial optimization framework. Concentrating on CRISPR-based diagnostics, we screen a library of 19,209 diagnostic-target pairs and train a deep neural network that predicts, from RNA sequence alone, diagnostic signal better than contemporary techniques. Our model then makes it possible to design assays that are maximally sensitive over the spectrum of a viruss genomic variation. We introduce ADAPT (https://adapt.guide), a system for fully-automated design, and use ADAPT to design optimal diagnostics for the 1,933 vertebrate-infecting viral species within 2 hours for most species and 24 hours for all but 3. We experimentally show ADAPTs designs are sensitive and specific down to the lineage level, including against viruses that pose challenges involving genomic variation and specificity. ADAPTs designs exhibit significantly higher fluorescence and permit lower limits of detection, across a viruss entire variation, than the outputs of standard design techniques. Our model-based optimization strategy has applications broadly to viral nucleic acid diagnostics and other sequence-based technologies, and, paired with clinical validation, could enable a critically-needed, proactive resource of assays for surveilling and responding to pathogens. oCEM: Automatic detection and analysis of overlapping co-expressed gene modules BackgroundWhen it comes to the co-expressed gene module detection, its typical challenges consist of overlap between identified modules and local co-expression in a subset of biological samples. The nature of module detection is the use of unsupervised clustering approaches and algorithms. Those methods are advanced undoubtedly, but the selection of a certain clustering method for sample- and gene-clustering tasks is separate, in which the latter task is often more complicated.

ResultsThis study presented an R-package, Overlapping CoExpressed gene Module (oCEM), armed with the decomposition methods to solve the challenges above. We also developed a novel auxiliary statistical approach to select the optimal number of principal components using a permutation procedure. We showed that oCEM outperformed state-of-the-art techniques in the ability to detect biologically relevant modules additionally.

ConclusionsoCEM helped non-technical users easily perform complicated statistical analyses and then gain robust results. oCEM and its applications, along with example data, were freely provided at https://github.com/huynguyen250896/oCEM.