Classification of daily crop phenology in PhenoCams using deep learning and hidden markov models Near surface cameras, such as those in the PhenoCam network, are a common source of ground truth data in modelling and remote sensing studies. Despite having locations across numerous agricultural sites, few studies have used near surface cameras to track the unique phenology of croplands. Due to management activities, crops do not have a natural vegetation cycle which many phenological extraction methods are based on. For example, a field may experience abrupt changes due to harvesting and tillage throughout the year. A single camera can also record several different plants due to crop rotations, fallow fields, and cover crops. Current methods to estimate phenology metrics from image time series compress all image information into a relative greenness metric, which discards a large amount of contextual information. This can include the type of crop present, whether snow or water is present on the field, the crop phenology, or whether a field lacking green plants consists of bare soil, fully senesced plants, or plant residue. Here we developed a modelling workflow to create a daily time series of crop type and phenology, while also accounting for other factors such as obstructed images and snow covered fields. We used a mainstream deep learning image classification model, VGG16. Deep learning classification models do not have a temporal component, so to account for temporal correlation among images our workflow incorporates a hidden markov model in the post-processing. The initial image classification model had out of sample F1 scores of 0.83-0.85, which improved to 0.86-0.91 after all post-processing steps. The resulting time series show the progression of crops from emergence to harvest, and can serve as a daily, local scale dataset of field states and phenological stages for agricultural research.
Embryo development in Carica papaya Linn. Papaya (Carica papaya Linn.) is a tropical plant whose draft genome has been sequenced. Papaya produces large fruits rich in vitamins A and C and is an important cash crop in developing countries. Nonetheless, little is known about how the female gametophyte develops, how it is fertilized and how it develops into a mature seed containing an embryo and an endosperm. The Papaya female gametophyte displays a Polygonum-type architecture consisting of two synergid cells, an egg cell, a central cell, and three antipodal cells. Reports are available of the presumed existence of varieties in which cross fertilization is bypassed and autonomous development of embryos occurs (e.g., apomixis). In this study, we analyzed the development of female gametophytes in a commercial Hawaiian parental line and in the presumed apomictic Costa Rican line L1. Samples were collected before and after anthesis to compare the overall structure, size and transcriptional patterns of several genes that may be involved in egg and endosperm cell fate and proliferation. These genes were the putative papaya homologs of ARGONAUTE9 (AGO9), MEDEA (MEA), RETINOBLASTOMA RELATED-1 (RBR1), and SLOW WALKER-1 (SWA1). Our results suggest that its feasibly to identify structural features of Polygonum-type development, and that in bagged female flowers of line L1 we might have observed autonomous development of embryo-like structures. Possible downregulation of papaya homologs for AGO9, MEA, RBR1 and SWA1 was observed in embryo sacs from line L1 before and after anthesis, which may suggest a link between suspected apomixis and transcriptional downregulation of genes for RNA-directed DNA methylation, histone remodelers, and rRNA processing. Most notably, the large size of the papaya embryo sac suggests that it could be a cytological alternative to Arabidopsis thaliana for study. Significant variation in embryo sac size was observed between the varieties under study, suggesting wide differences in the genetic regulation of anatomical features.
Soluble ECM promotes organotypic formation in lung alveolar model Micropatterned suspension culture creates consistently sized and shaped cell aggregates but has not produced organotypic structures from stable cells, thus restricting its use in accurate disease modeling. Here, we show that organotypic structure is achieved in hybrid suspension culture via supplementation of soluble extracellular matrix (ECM). We created a viable lung organoid from epithelial, endothelial, and fibroblast human stable cell lines in suspension culture. We demonstrate the importance of soluble ECM in organotypic patterning with the emergence of lumen-like structures with airspace showing feasible gas exchange units, formation of branching, perfusable vasculature, and long-term 70-day maintenance of lumen structure. Our results show a dependent relationship between enhanced fibronectin fibril assembly and the incorporation of ECM in the organoid. We successfully applied this technology in modeling lung fibrosis via bleomycin induction and test a potential antifibrotic drug in vitro while maintaining fundamental cell-cell interactions in lung tissue. Our human fluorescent lung organoid (hFLO) model represents features of pulmonary fibrosis which were ameliorated by fasudil treatment. We also demonstrate a 3D culture method with potential of creating organoids from mature cells, thus opening avenues for disease modeling and regenerative medicine, enhancing understanding of lung cell biology in health and lung disease. 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. Implicating Gene and Cell Networks Responsible for Differential COVID-19 Host Responses via an Interactive Single Cell Web Portal Numerous studies have provided single-cell transcriptome profiles of host responses to SARS-CoV-2 infection. Critically lacking however is a datamine that allows users to compare and explore cell profiles to gain insights and develop new hypotheses. To accomplish this, we harmonized datasets from COVID-19 and other control condition blood, bronchoalveolar lavage, and tissue samples, and derived a compendium of gene signature modules per cell type, subtype, clinical condition, and compartment. We demonstrate approaches to probe these via a new interactive web portal (http://toppcell.cchmc.org/ COVID-19). As examples, we develop three hypotheses: (1) a multicellular signaling cascade among alternatively differentiated monocyte-derived macrophages whose tasks include T cell recruitment and activation; (2) novel platelet subtypes with drastically modulated expression of genes responsible for adhesion, coagulation and thrombosis; and (3) a multilineage cell activator network able to drive extrafollicular B maturation via an ensemble of genes strongly associated with risk for developing post-viral autoimmunity. 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. Time-lapse imaging of bacterial morphology facilitates classification of antimicrobial agents Antimicrobial resistance is a major threat to human health. Basic knowledge of antimicrobial mechanism of action (MoA) is imperative for patient care and for identification of novel antimicrobials. However, the process of antimicrobial MoA identification is relatively laborious. Here, we developed a simple, quantitative time-lapse fluorescence imaging method, Dynamic Bacterial Morphology Imaging (DBMI), to facilitate this process. It uses a membrane dye and a nucleoid dye to track the morphological changes of single Bacillus subtilis cells in response to antimicrobials for up to 60 min. DBMI of bacterial cells facilitated assignment of the MoAs of 14 distinct, known antimicrobial compounds to the five main classes. Using this method, we found that the poorly studied antimicrobial, harzianic acid, a secondary metabolite that we purified from the fungal culture of Oidiodendron flavum, targets the cell envelope. We conclude that DBMI is a simple method, which facilitates rapid classification of the MoA of antimicrobials in functionally distinct classes. Tracheal tube fusion in Drosophila involves release of extracellular vesicles from multivesicular bodies Extracellular vesicles (EVs) comprise diverse types of cell-released membranous structures that are thought to play important roles in intercellular communication. Despite extensive work on the formation and functions of EVs in cultured cells, studies of EVs in vivo have remained scarce. We report here that EVs are present in the developing lumen of tracheal tubes in Drosophila embryos. We defined two distinct EV subpopulations, one of which contains the Munc13-4 homologue Staccato (Stac) and is spatially and temporally associated with tracheal tube fusion events. The formation of Stac-positive luminal EVs depends on the tip-cell-specific GTPase Arl3, which is also required for the formation of Stac-positive multivesicular bodies, suggesting that Stac-EVs derive from fusion of Stac-MVBs with the luminal membrane in tip cells during anastomosis formation. The GTPases Rab27 and Rab35 cooperate downstream of Arl3 to promote Stac-MVB formation and tube fusion. We propose that Stac-MVBs act as membrane reservoirs that facilitate tracheal lumen fusion in a process regulated by Arl3, Rab27, Rab35, and Stac/Munc13-4. Energy efficient network activity from disparate circuit parameters Neural circuits can produce similar activity patterns from vastly different combinations of channel and synaptic conductances. These conductances are tuned for specific activity patterns but might also reflect additional constraints, such as metabolic cost or robustness to perturbations. How do such constraints influence the range of permissible conductances? Here, we investigate how metabolic cost affects the parameters of neural circuits with similar activity in a model of the pyloric network of the crab Cancer borealis. We use a novel machine learning method to identify a range of network models that can generate activity patterns matching experimental data, and find that neural circuits can consume largely different amounts of energy despite similar circuit activity. Furthermore, a reduced but still significant range of circuit parameters gives rise to energy-efficient circuits. We then examine the space of parameters of energy-efficient circuits and identify potential tuning strategies for low metabolic cost. Finally, we investigate the interaction between metabolic cost and temperature robustness. We show that metabolic cost can vary across temperatures, but that robustness to temperature changes does not necessarily incur an increased metabolic cost. Our analyses show that, despite metabolic efficiency and temperature robustness constraining circuit parameters, neural systems can generate functional, efficient, and robust network activity with widely disparate sets of conductances. Multiparametric profiling of cellular lipid metabolism in hypercholesterolemia Systematic insight into cellular dysfunctions can improve understanding of disease etiology, risk assessment and patient stratification. We present a multiparametric high-content imaging platform enabling quantification of low-density lipoprotein (LDL) uptake and lipid storage in cytoplasmic droplets of primary leukocyte subpopulations. We validated this platform with samples from 65 individuals with variable blood LDL-cholesterol (LDL-c) levels, including familial hypercholesterolemia (FH) and non-FH subjects. We integrated lipid storage data into a novel readout, lipid mobilization, measuring the efficiency with which cells deplete lipid reservoirs. Lipid mobilization correlated positively with LDL uptake and negatively with hypercholesterolemia and age, improving differentiation of individuals with normal and elevated LDL-c. Moreover, combination of cell-based readouts with a polygenic risk score for LDL-c explained hypercholesterolemia better than the genetic risk score alone. This platform provides functional insights into cellular lipid trafficking from a few mls of blood and is applicable to dissect metabolic disorders, such as hypercholesterolemia.

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

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/440471v4_ufig1.gif" ALT="Figure 1">
View larger version (84K):
org.highwire.dtl.DTLVardef@9222c1org.highwire.dtl.DTLVardef@27e289org.highwire.dtl.DTLVardef@89ba82org.highwire.dtl.DTLVardef@33c9cd_HPS_FORMAT_FIGEXP  M_FIG C_FIG 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. The maize gene maternal derepression of r1 (mdr1) encodes a DNA glycosylase that demethylates DNA and reduces siRNA expression in endosperm Demethylation of transposons can activate expression of nearby genes and cause imprinted gene expression in endosperm, and it is hypothesized to lead to expression of transposon siRNAs that reinforce silencing in the next generation through transfer either into egg or embryo. Here we describe maternal derepression of r1 (mdr1), which encodes a DNA glycosylase with homology to Arabidopsis DEMETER and which is partially responsible for demethylation of thousands of regions in endosperm. Instead of promoting siRNA expression in endosperm, MDR1 activity inhibits it. Methylation of most repetitive DNA elements in endosperm is not significantly affected by MDR1, with an exception of Helitrons. While maternally-expressed imprinted genes preferentially overlap with MDR1 demethylated regions, the majority of genes that overlap demethylated regions are not imprinted. Double mutant megagametophytes lacking both MDR1 and its close homolog DNG102 result in early seed failure, and double mutant microgametophytes fail pre-fertilization. These data establish DNA demethylation by glycosylases as essential in maize endosperm and pollen and suggest that neither transposon repression nor genomic imprinting are its main function in endosperm. Beetle bioluminescence outshines aerial predators Bioluminescence is found across the tree of life and has many functions. Yet we understand very little about its timing and origins, particularly as a predator avoidance strategy. Understanding the timing between bioluminescence and predator origins has yet to be examined and can help elucidate the evolution of this ecologically important aposematic signal. Using the most prevalent bioluminescent group, fireflies, where bioluminescence primarily functions as aposematic and sexual signals, the timing for the origins of both potential predators of fireflies and bioluminescence is explored. Divergence time estimations were performed using genomic-scale phylogenetic reconstructions, and multiple fossil calibration points, allowing for a robust estimate for the origin of lampyrid bioluminescence as a terrestrial and as an aerial signal. Our results recover the origins of terrestrial beetle bioluminescence at 141.17 (122.63-161.17) mya and firefly aerial bioluminescence at 133.18 (117.86-152.47) mya with a large dataset focused on Lampyridae; and terrestrial bioluminescence as 148.03 (130.12-166.80) mya, with the age of aerial bioluminescence at 104.97 (99.00-120.90) mya with a complementary broad Elateroidea dataset. These ages predate the origins of all known extant aerial predators (i.e., bats and birds) and support the much older terrestrial predators (assassin bugs, frogs, ground beetles, lizards, snakes, hunting spiders, and harvestmen) as the most likely drivers of terrestrial bioluminescence in beetles, and sexual signaling likely being the original function in aerial fireflies. Low heritability and high phenotypic plasticity of salivary cortisol in response to environmental heterogeneity in a wild pinniped Individuals are unique in how they interact with and respond to their environment. Correspondingly, unpredictable challenges or environmental stressors often produce an individualized response of the hypothalamic-pituitary-adrenal axis and its downstream effector cortisol. We used a fully crossed, repeated measures design to investigate the factors shaping individual variation in baseline cortisol in Antarctic fur seal pups and their mothers. Saliva samples were collected from focal individuals at two breeding colonies, one with low and the other with high density, during two consecutive years of contrasting food availability. Mothers and pups were sampled concurrently at birth and shortly before weaning, while pups were additionally sampled every 20 days. We found that heritability was low for baseline cortisol, while within-individual repeatability and among-individual variability were high. A substantial proportion of the variation in baseline cortisol could be explained in pups and mothers by a combination of intrinsic and extrinsic factors including sex, weight, day, season, and colony of birth. Our findings provide detailed insights into the individualization of endocrine phenotypes and their genetic and environmental drivers in a wild pinniped. Furthermore, the strong associations between cortisol and life history traits that we report in fur seals could have important implications for understanding the population dynamics of species impacted by environmental change. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. Ecological interactions in breast cancer: Cell facilitation promotes growth and survival under drug pressure The interplay of positive and negative interactions between drug sensitive and resistant cells determines the efficacy of treatment in heterogeneous cancer cell populations. We developed isogenic estrogen receptor positive (ER+) breast cancer cell lineages sensitive and resistant to ribociclib-induced CDK4/6 inhibition to investigate cell interactions under selective pressure. During treatment with ribociclib alone or in combination with endocrine therapies, sensitive cells survive and proliferate better when grown together with resistant cells than when grown in monoculture. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) assays detected estrone, a less active metabolite, in medium taken from each lineage but detected highly active estradiol only from resistant cells. Because sensitive cells are more dependent on estrogen signaling than resistant cells, they benefit from the increase in local estradiol from resistant cells, leading to facilitated growth when grown under treatment in coculture. We developed a suite of mathematical models to quantify the strength of competition and facilitation between cell lines and evaluate the effects of treatments that modify these interactions. These analyses support the conclusion that resistant cell lineages facilitate growth of sensitive lineages during therapy through the increased production and sharing of local estradiol and themselves suffer from increased competition and suppression as a result. These models show that blocking facilitation has the potential to control both resistant and sensitive cell populations and inhibit the emergence of a refractory population. This combination of experiments and models shows that interactions within heterogeneous cancer cell populations can combine positive and negative interactions to shape the dynamics of resistance. CRISPR-based kinome-screening revealed MINK1 as a druggable player to rewire 5FU-resistance in OSCC through AKT/MDM2/p53 axis Cisplatin, 5FU and docetaxel (TPF) are the most common chemotherapy regimen used for advanced OSCC. However, many cancer patients experience relapse, continued tumor growth, and spread due to drug resistance, which leads to treatment failure and metastatic disease. Here, using a CRISPR/Cas9 based kinome knockout screening, Misshapen-like kinase 1 (MINK1) is identified as an important mediator of 5FU resistance in OSCC. Analysis of clinical samples demonstrated significantly higher MINK1 expression in the tumor tissues of chemotherapy non-responder as compared to chemotherapy responders. The nude mice and zebrafish xenograft experiments indicate that knocking out MINK1 restores 5FU mediated cell death in chemoresistant OSCC. An antibody based phosphorylation array screen revealed MINK1 as a negative regulator of p53. Mechanistically, MINK1 modulates AKT phosphorylation at Ser473, which enables p-MDM2 (Ser 166) mediated degradation of p53. We also identified lestaurtinib as a potent inhibitor of MINK1 kinase activity. The patient derived TPF resistant cell based xenograft data suggest that lestaurtinib restores 5FU sensitivity and facilitates a significant reduction of tumor burden. Overall, our study suggests that MINK1 is a major driver of 5FU resistance in OSCC. The novel combination of MINK1 inhibitor lestaurtinib and 5FU needs further clinical investigation in advanced OSCC. Inhibition of microbial deconjugation of micellar bile acids protects against intestinal permeability and liver injury Altered host-microbe interactions and increased intestinal permeability have been implicated in the pathogenesis of a range of diseases. However, the mechanisms by which gut microbes affect epithelial barrier integrity remain unclear. Few host-produced metabolites that protect against epithelial damage have been identified, and whether microbial metabolism of host factors alters intestinal barrier function is unknown. Here, we investigate the effects of bacterial metabolism of host-produced bile acid (BA) metabolites on epithelial barrier integrity. We observe that rats fed a choline-deficient, high-fat diet (CDAHFD) exhibit reduced abundance of host-produced conjugated BAs in the intestine at early timepoints coinciding with increased permeability. We show that in vitro, conjugated BAs protect gut epithelial monolayers from damage caused by bacterially produced unconjugated BAs through micelle formation. We then demonstrate that inhibition of BA deconjugation using a small molecule inhibitor of gut bacterial bile salt hydrolase (BSH) enzymes prevents development of pathologic intestinal permeability and hepatic inflammation in CDAHFD-fed rats. Finally, we show that the predominant conjugated BAs in humans protect against epithelial barrier disruption in vitro. Our study identifies a protective role for conjugated BAs in intestinal epithelial barrier function and suggests that rational manipulation of microbial BA metabolism could be leveraged to regulate gut barrier integrity. Bone mineral density modeling via random field: normality, stationarity, sex and age dependence Background and ObjectiveCapturing the population variability of bone properties is of paramount importance to biomedical engineering. The aim of the present paper is to describe variability and correlations in bone mineral density with a spatial random field inferred from routine computed tomography data.

MethodsRandom fields were simulated by transforming pairwise uncorrelated Gaussian random variables into correlated variables through the spectral decomposition of an age-detrended correlation matrix. The validity of the random field model was demonstrated in the spatiotemporal analysis of bone mineral density. The similarity between the computed tomography samples and those generated via random fields was analyzed with the energy distance metric.

ResultsThe random field of bone mineral density was found to be approximately Gaussian/slightly left-skewed/strongly right-skewed at various locations. However, average bone density could be simulated well with the proposed Gaussian random field for which the energy distance, i.e., a measure that quantifies discrepancies between two distribution functions, is convergent with respect to the number of correlation eigenpairs.

ConclusionsThe proposed random field model allows the enhancement of computational biomechanical models with variability in bone mineral density, which could increase the usability of the model and provides a step forward in in-silico medicine. Inhibitors of ERp44, PDIA1, and AGR2 induce disulfide-mediated oligomerization of Death Receptors 4 and 5 and cancer cell death Breast cancer mortality remains unacceptably high, indicating a need for safer and more effective therapeutic agents. Disulfide bond Disrupting Agents (DDAs) were previously identified as a novel class of anticancer compounds that selectively kill cancers that overexpress the Epidermal Growth Factor Receptor (EGFR) or its family member HER2. DDAs kill EGFR+ and HER2+ cancer cells via the parallel downregulation of EGFR, HER2, and HER3 and activation/oligomerization of Death Receptors 4 and 5 (DR4/5). However, the mechanisms by which DDAs mediate these effects are unknown. Affinity purification analyses employing biotinylated-DDAs reveal that the Protein Disulfide Isomerase (PDI) family members AGR2, PDIA1, and ERp44 are DDA target proteins. Further analyses demonstrate that shRNA-mediated knockdown of AGR2 and ERp44, or expression of ERp44 mutants, enhance basal and DDA-induced DR5 oligomerization. DDA treatment of breast cancer cells disrupts PDIA1 and ERp44 mixed disulfide bonds with their client proteins. Together, the results herein reveal DDAs as the first small molecule, active site inhibitors of AGR2 and ERp44, and demonstrate roles for AGR2 and ERp44 in regulating the activity, stability, and localization of DR4 and DR5, and activation of Caspase 8. Targeted biallelic integration of an inducible Caspase 9 suicide gene for safer cellular therapies prevents development of drug-resistant escapees in human iPSCs Drug-inducible suicide systems may help to minimize risks of cellular therapies due to the tumor forming potential of human induced pluripotent stem cells (hiPSCs). Recent research challenged the usefulness of such systems since rare drug-resistant subclones were observed that showed elimination or silencing of the transgene.

We have introduced a drug-inducible Caspase9 suicide system (iCASP9) into the AAVS1 safe harbor locus of hiPSCs. In these cells, apoptosis could be efficiently induced in vitro. In mice, drug treatment generally led to rapid elimination of teratomas, but individual animals subsequently formed tumor tissue from monoallelic iCASP9 hiPSCs. Very rare drug-resistant subclones of monoallelic iCASP9 hiPSCs appeared in vitro with frequencies of ~ 3x10-8. Transgene elimination, presumably via Loss of Heterozygosity (LoH), or methylation of the CAG promoter but not methylation of the O_SCPLOWPPPC_SCPLOWO_SCPCAP1C_SCPCAPO_SCPLOWRC_SCPLOWO_SCPCAP12C_SCPCAPO_SCPLOWCC_SCPLOW locus were identified as underlying mechanisms. In contrast, we never observed any escapees from biallelic iCASP9 hiPSCs, even after treatment of up to 0.8 billion hiPSCs.

In conclusion, biallelic integration of an iCASP9 system in the AAVS1 locus may substantially contribute to the safety level of iPSC-based therapies, which should be calculated by relating clonal escapee frequencies to the cell number in tumors of a size that is readily detectable during routine screening procedures. Evaluation of methods for estimating coalescence times using ancestral recombination graphs The ancestral recombination graph (ARG) is a structure that describes the joint genealogies of sampled DNA sequences along the genome. Recent computational methods have made impressive progress towards scalably estimating whole-genome genealogies. In addition to inferring the ARG, some of these methods can also provide ARGs sampled from a defined posterior distribution. Obtaining good samples of ARGs is crucial for quantifying statistical uncertainty and for estimating population genetic parameters such as effective population size, mutation rate, and allele age. Here, we use standard neutral coalescent simulations to benchmark the estimates of pairwise coalescence times from three popular ARG inference programs: ARGweaver, Relate, and tsinfer+tsdate. We compare 1) the true coalescence times to the inferred times at each locus; 2) the distribution of coalescence times across all loci to the expected exponential distribution; 3) whether the sampled coalescence times have the properties expected of a valid posterior distribution. We find that inferred coalescence times at each locus are most accurate in ARGweaver, and often more accurate in Relate than in tsinfer+tsdate. However, all three methods tend to overestimate small coalescence times and underestimate large ones. Lastly, the posterior distribution of ARGweaver is closer to the expected posterior distribution than Relates, but this higher accuracy comes at a substantial trade-off in scalability. The best choice of method will depend on the number and length of input sequences and on the goal of downstream analyses, and we provide guidelines for the best practices. The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. D155Y Substitution of SARS-CoV-2 ORF3a Weakens Binding with Caveolin-1 The clinical manifestation of the recent pandemic COVID-19, caused by the novel SARS-CoV-2 virus, varies from mild to severe respiratory illness. Although environmental, demographic and co-morbidity factors have an impact on the severity of the disease, contribution of the mutations in each of the viral genes towards the degree of severity needs a deeper understanding for designing a better therapeutic approach against COVID-19. Open Reading Frame-3a (ORF3a) protein has been found to be mutated at several positions. In this work, we have studied the effect of one of the most frequently occurring mutants, D155Y of ORF3a protein, found in Indian COVID-19 patients. Using computational simulations we demonstrated that the substitution at 155th changed the amino acids involved in salt bridge formation, hydrogen-bond occupancy, interactome clusters, and the stability of the protein compared with the other substitutions found in Indian patients. Protein-protein docking using HADDOCK analysis revealed that substitution D155Y weakened the binding affinity of ORF3a with caveolin-1 compared with the other substitutions, suggesting its importance in the overall stability of ORF3a-caveolin-1 complex, which may modulate the virulence property of SARS-CoV-2. Conformist social learning inoculates against adverse risk aversion Given the ubiquity of potentially adverse behavioural bias owing to myopic trial-and-error learning, it seems paradoxical that improvements in decision-making performance through conformist social learning, a process widely considered to be bias amplification, still prevail in animal collective behaviour. Here we show, through model analyses and large-scale interactive behavioural experiments with 585 human subjects, that conformist influence can indeed promote favourable risk taking in repeated experience-based decision making, even though many individuals are systematically biased towards adverse risk aversion. Although strong positive feedback conferred by copying the majoritys behaviour could result in unfavourable informational cascades, our differential equation model of collective behavioural dynamics identified a key role for increasing exploration by negative feedback arising when a weak minority influence undermines the inherent behavioural bias. This  collective behavioural rescue, emerging through coordination of positive and negative feedback, highlights a benefit of collective learning in a broader range of environmental conditions than previously assumed and resolves the ostensible paradox of adaptive collective behavioural flexibility under conformist influences. Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less-functional Misfolded protein conformations with decreased functionality can bypass the proteostasis machinery and remain soluble in vivo. This is an unexpected phenomenon as several cellular quality control mechanisms have evolved to rid cells of misfolded proteins. Three questions, then, are: how is it structurally possible for long-lived, soluble, misfolded proteins to bypass the proteostasis machinery and processes? How widespread are these soluble, misfolded states across the proteome? And how long do they persist for? Here, we address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of all proteins exhibit subpopulations of misfolded conformations that are likely to bypass molecular chaperones, avoid aggregation, and not be rapidly degraded. These misfolded states may persist for months or longer for some proteins. Structurally characterizing these misfolded states, we observe they have a large amount of native structure, but also contain localized misfolded regions from non-native changes in entanglement, in which a protein segment threads through a loop formed by another portion of the protein that is not found in the native state. The surface properties of these misfolded states are native like, suggesting they may bypass the proteostasis machinery and its regulatory processes to remain soluble, while their entanglements make these states long-lived kinetic traps, as disentanglement requires unfolding of already folded portions of the protein. In terms of function, we predict that one-third of proteins have subpopulations that misfold into less-functional states that have structurally perturbed functional sites yet remain soluble. Data from limited-proteolysis mass spectrometry experiments, which interrogate the misfolded conformations populated by proteins upon unfolding and refolding, are consistent with the structural changes seen in the entangled states of glycerol-3-phosphate dehydrogenase upon misfolding. These results provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass cellular quality controls, and indicate, unexpectedly, this may be a wide-spread phenomenon in proteomes. Such entanglements are observed in many native structures, suggesting the non-native entanglements we observe are plausible. More broadly, these near-native entangled structures suggest a hypothesis for how synonymous mutations can modulate downstream protein structure and function, with these mutations partitioning nascent proteins between these kinetically trapped states.