An image reconstruction framework for characterizing initial visual encoding We developed an image-computable observer model of the initial visual encoding that operates on natural image input, based on the framework of Bayesian image reconstruction from the excitations of the retinal cone mosaic. Our model extends previous work on ideal observer analysis and evaluation of performance beyond psychophysical discrimination, takes into account the statistical regularities of the visual environment, and provides a unifying framework for answering a wide range of questions regarding the visual front end. Using the error in the reconstructions as a metric, we analyzed variations of the number of different photoreceptor types on human retina as an optimal design problem. In addition, the reconstructions allow both visualization and quantification of information loss due to physiological optics and cone mosaic sampling, and how these vary with eccentricity. Furthermore, in simulations of color deficiencies and interferometric experiments, we found that the reconstructed images provide a reasonable proxy for modeling subjects percepts. Lastly, we used the reconstruction-based observer for the analysis of psychophysical threshold, and found notable interactions between spatial frequency and chromatic direction in the resulting spatial contrast sensitivity function. Our method is widely applicable to experiments and applications in which the initial visual encoding plays an important role.
Region-specific microRNA alterations in marmosets carrying SLC6A4 polymorphisms are associated to anxiety-like behavior. Psychiatric diseases such as depression and anxiety are multifactorial conditions, highly prevalent in western societies. Human studies have identified a number of high-risk genetic variants for these diseases. Among them, polymorphisms in the promoter region of the serotonin transporter gene (SLC6A4) have attracted much attention. However, due to the paucity of experimental models, molecular alterations induced by these genetic variants and how they correlate to behavioral deficits have not been examined. Marmosets have emerged as a powerful model in translational neuroscience to investigate molecular underpinnings of complex behaviors. Here, we took advantage of naturally occurring genetic polymorphisms in marmoset SLC6A4 gene that have been linked to anxiety-like behaviors. Using FACS-sorted cells from different brain regions, we revealed that marmosets bearing different SLC6A4 variants exhibit distinct microRNAs signatures in a region of the prefrontal cortex whose activity has been consistently altered in patients with depression/anxiety. We also identified DCC, a gene previously linked to these diseases, as a downstream target of the dysregulated microRNAs. Significantly, we showed that levels of both microRNAs and DCC in this region were highly correlated to anxiety-like behaviors as well as to the response to citalopram, a selective serotonin re-uptake inhibitor and widely prescribed anti-depressant. Our findings establish links between genetic variants, molecular modifications in specific cortical regions and complex behavioral/pharmacological responses, providing new insights into gene-behavior relationships underlying human psychopathology.
Development of a Novel PROTAC using the Nucleic Acid Aptamer as a Targeting Ligand for Tumor Selective Degradation of Nucleolin PROteolysis TArgeting Chimeras (PROTACs) induce targeted protein degradation by hijacking the intracellular ubiquitin proteasome system, thus emerging as a new strategy for drug development. However, most PROTACs generated lack cell-type selectivity and have poorly soluble in water. To address this drawback, we developed a novel PROTAC ZL216 using aptamer AS1411 as a targeting ligand of nucleolin to conjugate with a small molecule ligand of E3 ligase VHL, which shows high aqueous solubility and serum stability. Based on the differential expression of nucleolin on the cell surface, ZL216 could bind to and internalize into breast cancer cells, but not normal breast cells. Furthermore, we revealed that ZL216 promoted the formation of a nucleolin-ZL216-VHL ternary complex in breast cancer cells and potently induced nucleolin degradation in vitro and in vivo. As a result, ZL216 inhibited the proliferation and migration of breast cancer cells. These studies demonstrate that in addition to peptides and small molecule compounds, nuclei acid aptamers can also be used to generate PROTACs, which broadens the toolbox constructing PROTACs and provides a promising strategy for development of tumor-selective PROTACs. A potent human monoclonal antibody with pan-neutralizing activities directly dislocates S trimer of SARS-CoV-2 through binding both up and down forms of RBD The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a global pandemic of novel corona virus disease (COVID-19). The neutralizing monoclonal antibodies (mAbs) targeting the receptor binding domain (RBD) of SARS-CoV-2 are among the most promising strategies to prevent and treat COVID-19. However, SARS-CoV-2 variants of concern (VOCs) profoundly reduced the efficacies of most of mAbs and vaccines approved for clinical use. Herein, we demonstrated mAb 35B5 efficiently neutralizes both wild-type (WT) SARS-CoV-2 and VOCs, including B.1.617.2 (delta) variant, in vitro and in vivo. Cryo-electron microscopy (cryo-EM) revealed that 35B5 neutralizes SARS-CoV-2 by targeting a unique epitope that avoids the prevailing mutation sites on RBD identified in circulating VOCs, providing the molecular basis for its pan-neutralizing efficacy. The 35B5-binding epitope could also be exploited for the rational design of a universal SARS-CoV-2 vaccine. 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. A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Identification of COVID-19 prognostic markers and therapeutic targets through meta-analysis and validation of Omics data from nasopharyngeal samples While our battle with the COVID-19 pandemic continues, a multitude of Omics data has been generated from patient samples in various studies, which remains to be translated. We conducted a meta-analysis of published transcriptome and proteome profiles of nasal swab and bronchioalveolar lavage fluid (BALF) samples of COVID-19 patients, to shortlist high confidence upregulated host factors. Subsequently, mRNA overexpression of selected genes was validated in nasal swab/BALF samples from a cohort of COVID-19 positive/negative, symptomatic/asymptomatic individuals. Analysis of these data revealed S100 family genes (S100A6, S100A8, S100A9, and S100P) as prognostic markers of COVID-19 disease. Furthermore, Thioredoxin gene (TXN) was identified as a significant upregulated host factor in our overlap analysis. An FDA-approved drug Auranofin, which inhibits Thioredoxin reduction, was found to mitigate SARS-CoV-2 replication in vitro and in vivo in the hamster challenge model. Overall, this study translates COVID-19 host response Big Data into potential clinical interventions. Sexual selection for males with beneficial mutations Sexual selection is the process by which traits providing a mating advantage are favoured. Theoretical treatments of the evolution of sex by sexual selection propose that it operates by reducing the load of deleterious mutations. Here, we postulate instead that sexual selection primarily acts through females preferentially mating with males carrying beneficial mutations. We used simulation and analytical modelling to investigate the evolutionary dynamics of beneficial mutations in the presence of sexual selection. We found that female choice for males with beneficial mutations had a much greater impact on genetic quality than choice for males with low mutational load. We also relaxed the typical assumption of a fixed mutation rate. For deleterious mutations, mutation rate should always be minimized, but when rare beneficial mutations can occur, female choice for males with those rare beneficial mutations could overcome a decline in average fitness and allow an increase in mutation rate. We propose that sexual selection for beneficial mutations could overcome the  two-fold cost of sex much more readily than choice for males with low mutational load and may therefore be a more powerful explanation for the prevalence of sexual reproduction than the existing theory. If sexual selection results in higher fitness at higher mutation rates and if the variability produced by mutation itself promotes sexual selection, then a feedback loop between these two factors could have had a decisive role in driving adaptation. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. 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. Maf/ham1-like pyrophosphatases of non-canonical nucleotides are host-specific partners of viral RNA-dependent RNA polymerases Cassava brown streak disease (CBSD), dubbed the "Ebola of plants", is a serious threat to food security in Africa caused by two viruses of the family Potyviridae: cassava brown streak virus (CBSV) and Ugandan (U)CBSV. Intriguingly, U/CBSV, along with another member of this family and one secoviridae, are the only known RNA viruses encoding a protein of the Maf/ham1-like family, a group of widespread pyrophosphatase of non-canonical nucleotides (ITPase) expressed by all living organisms. Despite the socio-economic impact of CDSD, the relevance and role of this atypical viral factor has not been yet established. Here, using an infectious cDNA clone and reverse genetics, we demonstrate that UCBSV requires the ITPase activity for infectivity in cassava, but not in the model plant Nicotiana benthamiana. HPLC MS/MS experiments showed that, quite likely, this host-specific constraint is due to an unexpected high concentration of non-canonical nucleotides in cassava. Finally, protein analyses and experimental evolution of mutant viruses indicated that keeping a fraction of the yielded UCBSV ITPase covalently bound to the viral RNA-dependent RNA polymerase (RdRP) optimizes viral fitness, and this seems to be a feature shared by the other members of the Potyviridae family expressing Maf/ham1-like proteins. All in all, our work (i) reveals that the over-accumulation of non-canonical nucleotides in the host might have a key role in antiviral defense, and (ii) provides the first example of an RdRP-ITPase partnership, reinforcing the idea that RNA viruses are incredibly versatile at adaptation to different host setups. Paradoxical Roles of Peritumoral Myofibroblasts in Intrahepatic Cholangiocarcinoma Growth and Metastasis Intrahepatic cholangiocarcinoma (iCCA) is an aggressive type of primary liver cancer characterized by its highly desmoplastic stroma. Compared to the ample knowledge on cancer-associated fibroblast residing within iCCA tumor mass, little is known about the fibrotic response from the tumor-surrounding liver and its role in iCCA development. In this study, we find a significant accumulation of peritumoral myofibroblasts (pMFs) in both patient iCCA and tumors from an iCCA orthotopic allograft mouse model. Using tumor and liver spheroid coculture we show that iCCA-liver interaction induces rapid pMF accumulation at the interface. We find pMFs placed around iCCA spheroids exert a strong suppressive effect on tumor cells growth, in contrast to the pro-proliferative effect of MFs mixed within tumor spheroids. However, prolonged iCCA-pMF interaction elicits tumor cell dissemination in vitro. We find an upregulation of Vcam1 in tumor cells in the early phase of iCCA-pMF interaction both in vitro and in vivo which is downregulated when tumor cells disseminate. Blocking Vcam1 activity in iCCA allograft mouse models slows primary tumor growth but lead to increased tumor metastasis. Our data suggest that pMFs are beyond simple pro- or anti-tumorigenic in iCCA, with the ability to suppress tumor growth but elicit tumor cell dissemination. HSP90 modulates T2R bitter taste receptor nitric oxide production and innate immune responses in human airway epithelial cells and macrophages Bitter taste receptors (T2Rs) are G protein-coupled receptors (GPCRs) expressed in various cell types including ciliated airway epithelial cells and macrophages. T2Rs in these two airway innate immune cell types are activated by bitter products, including some secreted by Pseudomonas aeruginosa, leading to Ca2+-dependent activation of endothelial nitric oxide (NO) synthase (eNOS). NO enhances mucociliary clearance and has direct antibacterial effects in ciliated epithelial cells and increases phagocytosis by macrophages. Using biochemistry and live cell imaging, we explored the role of heat shock protein 90 (HSP90) in regulating T2R-dependent NO pathways in primary sinonasal epithelial cells, primary monocyte-derived macrophages, and a human bronchiolar cell line (H441). We used immunofluorescence to show that H441 cells express eNOS and certain T2Rs and that the bitterant denatonium benzoate activates NO production in an HSP90-dependent manner in cells grown either as submerged cultures and at air liquid interface. In primary sinonasal epithelial cells, we determined that HSP-90 inhibition reduces T2R-stimulated NO production and ciliary beating which are crucial for pathogen clearance. In primary monocyte-derived macrophages, we found that HSP-90 is integral to T2R-stimulated NO production and phagocytosis of FITC-labeled Escherichia coli and pHrodo-Staphylococcus aureus. Our study demonstrates that HSP90 serves an innate immune role by regulating NO production downstream of T2R signaling by augmenting eNOS activation without impairing upstream calcium signaling. These findings suggest that HSP90 plays an important role in airway antibacterial innate immunity and may be an important target in airway diseases like chronic rhinosinusitis, asthma, or cystic fibrosis. Multiple brace root phenotypes promote anchorage and limit root lodging in maize Plant mechanical failure (lodging) causes global yield losses of 7-66% in cereal crops. We have previously shown that the above-ground nodal roots (brace roots) in maize are critical for anchorage. However, it is unknown how brace root phenotypes vary across genotypes and the functional consequence of this variation. This study quantifies the contribution of brace roots to anchorage, brace root traits, plant height, and root lodging susceptibility in 52 maize inbred lines. We show that the contribution of brace roots to anchorage and root lodging susceptibility varies among genotypes and this contribution can be explained by plant architectural variation. Additionally, supervised machine learning models were developed and show that multiple plant architectural phenotypes can predict the contribution of brace roots to anchorage and root lodging susceptibility. Together these data define the plant architectures that are important in lodging resistance and show that the contribution of brace roots to anchorage is a good proxy for root lodging susceptibility. Darwin's small and medium ground finches might have taste preferences, but not for human foods Urbanization is rapidly changing ecological niches. On the inhabited Galapagos Islands, Darwins finches consume human-introduced foods preferentially; however, it remains unclear why. Here we presented pastry with flavour profiles typical of human foods (oily, salty, sweet) to small ground finches (Geospiza fuliginosa) and medium ground finches (Geospiza fortis) to test if latent taste preferences might drive selection of human foods. If human-food flavours were consumed more than a neutral or bitter control only at sites with human foods, then we predicted tastes were acquired after urbanization; however, if no site-differences were found then this would indicate latent taste preferences. Contrary to both predictions, we found little evidence that human-food flavours were preferred compared to control flavours at any site. Instead, finches showed a weak aversion to oily foods, but only at remote (no human foods present) sites. This was further supported by behavioural responses, with beak wiping occurring more often at remote sites after finches tasted flavours associated with human foods. Our results suggest, therefore, that while Darwins finches may have acquired a tolerance to the flavours of human food, latent taste preferences are unlikely to have played a major role in their dietary response to increased urbanization. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. Modelling porcine reproductive and respiratory syndrome virus dissemination dynamics to quantifythe contribution of multiple modes of transmission: between-farm pig and vehicle movements, farmto-farm proximity, feed ingredients, and re-break Accounting for multiple modes of livestock disease dissemination in epidemiological models remains a challenge. We developed and calibrated a mathematical model for transmission of porcine reproductive and respiratory syndrome virus (PRRSV), tailored to fit nine modes of between-farm transmission pathways including: farm-to-farm proximity (local transmission), contact network of batches of pigs transferred between farms (pig movements), re-break probabilities for farms with previous PRRSV outbreaks, with the addition of four different contact networks of transportation vehicles (vehicles to transport pigs to farms, pigs to markets, feed and crew) and the amount of animal by-products within feed ingredients (e.g. animal fat or meat and bone meal). The model was calibrated on weekly PRRSV outbreaks data. We assessed the role of each transmission pathway considering the dynamics of specific types of production (i.e., sow, nursery). Although our results estimated that the networks formed by transportation vehicles were more densely connected than the network of pigs transported between-farms, pig movements and farm proximity were the main PRRSV transmission routes regardless of farm types. Among the four vehicle networks, vehicles transporting pigs to farms explained a large proportion of infections, sow = 20.9%; nursery = 15%; and finisher = 20.6%. The animal by-products showed a limited association with PRRSV outbreaks through descriptive analysis, and our model results showed that the contribution of animal fat contributed only 2.5% and meat and bone meal only 0.03% of the infected sow farms. Our work demonstrated the contribution of multiple routes of PRRSV dissemination, which has not been deeply explored before. It also provides strong evidence to support the need for cautious, measured PRRSV control strategies for transportation vehicles and further research for feed by-products modeling. Finally, this study provides valuable information and opportunities for the swine industry to focus effort on the most relevant modes of PRRSV between-farm transmission. Structural variant-based pangenome construction has low sensitivity to variability of haplotype-resolved bovine assemblies Advantages of pangenomes over linear reference assemblies for genome research have recently been established. However, potential effects of sequence platform and assembly approach, or of combining assemblies created by different approaches, on pangenome construction have not been investigated. We generated haplotype-resolved assemblies from the offspring of three bovine trios representing increasing levels of heterozygosity that each demonstrate a substantial improvement in contiguity, completeness, and accuracy over the current Bos taurus reference genome. Diploid coverage as low as 20x for HiFi or 60x for ONT was sufficient to produce two haplotype-resolved assemblies meeting the standards set by the Vertebrate Genome Project. Structural variant-based pangenomes created from the haplotype-resolved assemblies demonstrated significant consensus regardless of sequence platform, assembler algorithm, or coverage. Inspecting pangenome topologies identified 90 thousand structural variants including 931 overlapping with coding sequences; this approach revealed variants affecting QRICH2, PRDM9, HSPA1A, TAS2R46, and GC that have potential to affect phenotype. 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. RapidoPGS: A rapid polygenic score calculator for summary GWAS data without a test dataset MotivationPolygenic scores (PGS) aim to genetically predict complex traits at an individual level. PGS are typically trained on genome-wide association summary statistics and require an independent test dataset to tune parameters. More recent methods allow parameters to be tuned on the training data, removing the need for independent test data, but approaches are computationally intensive. Based on fine-mapping principles, we present RapidoPGS, a flexible and fast method to compute PGS requiring summary-level GWAS datasets only, with little computational requirements and no test data required for parameter tuning.

ResultsWe show that RapidoPGS performs slightly less well than two out of three other widely-used PGS methods (LDpred2, PRScs, and SBayesR) for case-control datasets, with median r2 difference: -0.0092, -0.0042, and 0.0064, respectively, but up to 17,000-fold faster with reduced computational requirements. RapidoPGS is implemented in R and can work with user-supplied summary statistics or download them from the GWAS catalog.

Availability and implementationOur method is available with a GPL license as an R package from GitHub. The pH Dependence of Niclosamide Solubility, Dissolution, and Morphology Motivates Potentially Universal Mucin-Penetrating Nasal and Throat Sprays for COVID19, its Contagious Variants, and Other Respiratory Viral Infections MotivationWith the coronavirus pandemic still raging, prophylactic nasal and early treatment throat sprays could help prevent infection and reduce viral load. Niclosamide has the potential to treat a broad range of viral infections if local bioavailability is optimized as mucin-penetrating solutions as opposed to more traditional microparticle-based sprays that cannot penetrate the mucin.

ExperimentalpH-dependence of supernatant concentrations and dissolution rates of niclosamide were measured in buffered solutions by Nanodrop-UV/Vis-spectroscopy for niclosamide from different suppliers, as precipitated material, and as cosolvates. Data was compared to predictions from Henderson Hasselbalch and precipitation pH models. Optimal microscopy was used to observe the morphologies of precipitated and converted niclosamide.

ResultsSupernatant-concentrations of niclosamide increased with increasing pH, from 1.77uM at pH 3.66 to 30uM at pH 8, and more rapidly from 90uM at pH8.5 to 300uM at pH 9.1, reaching 641uM at pH 9.5. Logarithmic rates for dissolution increased by [~]3x for pHs 8.62 to 9.44. However, when precipitated from supersaturated solution, niclosamide equilibrated to much lower final supernatant concentrations, reflective of more stable polymorphs at each pH that were also apparent for niclosamide from other suppliers and cosolvates.

ConclusionsGiven niclosamides activity against COVID19, its more contagious variants, and other respiratory viral infections, these niclosamide solutions, that put the virus in lockdown, could represent universal prophylactic nasal and early treatment throat sprays. As solutions they would be the simplest and potentially most effective formulations from both an efficacy standpoint as well as manufacturing and distribution, with no cold chain. They now just need testing. The Digital 3D-Atlas MAKER (DAMAKER): a dynamic and expandable digital 3D-tool for monitoring the temporal changes in tissue growth during hindbrain morphogenesis Reconstruction of prototypic three-dimensional (3D) atlases at the scale of whole tissues or organs requires specific methods to be developed. We have established a digital 3D-atlas maker (DAMAKER) and built a digital 3D-temporal atlas to monitor the changes in the growth of the tissue and neuronal differentiation domain in the zebrafish hindbrain. DAMAKER integrates spatial and temporal data from cell populations, neuronal differentiation and brain morphogenesis, through in vivo imaging techniques paired with image analyses and segmentation tools. First, we generated a 3D-reference from several imaged hindbrains and segmented them using a trainable tool; these were aligned using rigid registration, revealing distribution of neuronal differentiation growth patterns along the axes. Second, we quantified the dynamic growth of the neuronal differentiation domain vs. the progenitor domain, and by in vivo neuronal birthdating experiments we generated a digital 3D-temporal map of the neuronal growth in the whole hindbrain, revealing the spatiotemporal dynamics of neuronal differentiation upon morphogenesis. Last, we applied it to glutamatergic and GABAergic neurons, as proof-of-concept that the digital 3D-temporal map could be used as a proxy to infer neuronal birthdate. As this protocol uses open-access tools and algorithms, it can be shared for standardized, accessible, tissue-wide cell population atlas construction. 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. Regulation of Myoepithelial Differentiation in 3-Dimensional Culture The salivary gland can be permanently impaired by radiation treatment for head and neck cancers. Efforts at tissue regeneration have focused on saliva-producing acinar cells. However, myoepithelial cells are also critical to gland function, but mechanisms that regulate their differentiation are poorly defined. To study myoepithelial differentiation, we employed mSG-PAC1 murine salivary gland epithelial cells. We demonstrate that mSG-PAC1 spheroids exhibit phenotypic plasticity between pro-acinar and myoepithelial cell fates. Increased expression of pro-acinar/acinar or myoepithelial RNAs was identified from spheroids cultured under different media conditions by microarray followed by gene-set enrichment analysis. Spheroids cultured with different medium components expressed proteins typical of either acinar or myoepithelial cells, as detected by immunocytochemistry. We demonstrate that the pattern of TAZ expression in the epithelial compartment of the differentiating murine salivary gland correlates with the expression of the myoepithelial marker alpha-SMA, as is the case for TAZ expression in mSG-PAC1 spheroids. Our analysis also indicates that YAP/TAZ target genes are upregulated together with myoepithelial markers. Importantly, siRNA targeting of TAZ expression in mSG-PAC1 spheroids diminished the expression of myoepithelial markers. Our results in this in vitro cell model implicate TAZ signaling in myoepithelial differentiation. A general framework to learn tertiary structure for protein sequence characterization During the past five years, deep-learning algorithms have enabled ground-breaking progress towards the prediction of tertiary structure from a protein sequence. Very recently, we developed SAdLSA, a new computational algorithm for protein sequence comparison via deep-learning of protein structural alignments. SAdLSA shows significant improvement over established sequence alignment methods. In this contribution, we show that SAdLSA provides a general machine-learning framework for structurally characterizing protein sequences. By aligning a protein sequence against itself, SAdLSA generates a fold distogram for the input sequence, including challenging cases whose structural folds were not present in the training set. About 70% of the predicted distograms are statistically significant. Although at present the accuracy of the intra-sequence distogram predicted by SAdLSA self-alignment is not as good as deep-learning algorithms specifically trained for distogram prediction, it is remarkable that the prediction of single protein structures is encoded by an algorithm that learns ensembles of pairwise structural comparisons, without being explicitly trained to recognize individual structural folds. As such, SAdLSA can not only predict protein folds for individual sequences, but also detects subtle, yet significant, structural relationships between multiple protein sequences using the same deep-learning neural network. The former reduces to a special case in this general framework for protein sequence annotation. Probing the effect of glycosaminoglycan depletion on integrin interactions with collagen I fibrils in the native ECM environment. Fibrillar collagen-integrin interactions in the extracellular matrix (ECM) regulate a multitude of cellular processes and cell signalling. Collagen I fibrils serve as the molecular scaffolding for connective tissues throughout the human body and are the most abundant protein building blocks in the ECM. The ECM environment is diverse, made up of several ECM proteins, enzymes, and proteoglycans. The contents of the ECM environment are modulated by disease and aging and may influence these critical collagen-integrin interactions. In particular, glycosaminoglycans (GAGs), anionic polysaccharides that decorate proteoglycans, become depleted in the ECM with natural aging and their mis-regulation has been linked to cancers and other diseases. The impact of GAG concentration in the ECM environment on collagen interactions is not well understood. Here, we integrate protein adhesion assays with liquid high resolution atomic force microscopy (AFM) to assess the affects of GAG depletion on the interaction of collagen I fibrils with the integrin 2I domain. Adhesion assays demonstrate that 2I preferentially binds to GAG-depleted collagen I fibrils. By amplitude modulated AFM in air and in solution, we find that GAG-depleted collagen I fibrils retain structural features of the native fibrils, including their characteristic D-banding pattern, a key structural motif. AFM fast force mapping in solution shows that GAG depletion reduces the stiffness of individual fibrils, lowering the indentation modulus by half compared to native fibrils. Together these results shed new light on how GAGs influence collagen I fibril- integrin interactions and may aid in strategies to treat diseases that result from GAG misregulation.

Statement for broader audienceAging and disease result in mis-regulation of glycosaminoglycan (GAG) levels in the extracellular matrix (ECM), which may affect fibrillar collagen interactions that are vital for cellular processes. Here, we characterize the impact of GAG depletion on collagen-integrin 2I domain interactions and collagen fibril topography and stiffness. We show that GAG depletion increases collagen-2I binding and reduces stiffness in comparison to native fibrils. These results may inform on strategies for treating GAG mis-regulation.