Multivariate Genomewide Association Analysis by Iterative Hard Thresholding 1In genome-wide association studies (GWAS), analyzing multiple correlated traits is potentially superior to conducting multiple univariate analyses. Standard methods for multivariate GWAS operate marker-by-marker and are computationally intensive. We present a penalized regression algorithm for multivariate GWAS based on iterative hard thresholding (IHT) and implement it in a convenient Julia package MendelIHT.jl (https://github.com/OpenMendel/MendelIHT.jl). In simulation studies with up to 100 traits, IHT exhibits similar true positive rates, smaller false positive rates, and faster execution times than GEMMAs linear mixed models and mv-PLINKs canonical correlation analysis. On UK Biobank data, our IHT software completed a 3-trait joint analysis in 20 hours and an 18-trait joint analysis in 53 hours, requiring up to 80GB of computer memory. In short, our software enables geneticists to fit a single regression model that simultaneously considers the effect of all SNPs and dozens of traits.
Atg11 promotes Kar9-mediated alignment of the mitotic spindle Asymmetric spindle pole body (SPB) inheritance requires a cascade of events that involve kinases, phosphatases and structural scaffold proteins including molecular motors and microtubule-associated proteins. Higher levels of an SPB component Spc72 and the spindle positioning factor Kar9 at the old SPB, which migrates to the daughter cell, ensure asymmetric SPB inheritance. In this study, by analyzing the protein interaction network of autophagy (Atg)-related proteins with those that constitute the chromosome segregation machinery, and growth inhibition kinetics in the presence of a microtubule poison, we identified Atg11 as a potential regulator of chromosome transmission from a cohort of 35 autophagy mutants in Saccharomyces cerevisiae. Cells lacking Atg11 did not show kinetochore defects but displayed a high rate of chromosome loss and delayed anaphase onset. Atg11 positively interacted with Kar9 and Kip2 and negatively with Dyn1 and Kar3 in mediating proper positioning and alignment of the mitotic spindle that was found to be under the surveillance of the spindle positioning checkpoint. Indeed, atg11{Delta} cells displayed an inverted SPB inheritance. Our experiments further revealed that via the Kar9 pathway, Atg11 promotes asymmetric localization of Spc72 and Kar9 on the old SPB. Atg11 physically interacted with Spc72 and transiently localized close to the old SPB during metaphase-to-anaphase progression. Thus, we uncover an autophagy-independent role of Atg11 in spindle alignment. This study emphasizes the importance of screens to identify factors mediating the complex and intricate crosstalk between biological processes fundamental to genomic integrity.
Phylogenomic analyses of the East Asian endemic Abelia (Caprifoliaceae) shed insights into the temporal and spatial diversification history with widespread hybridizations Background and AimsAbelia (Caprifoliaceae) is a small genus with five species, including one man-made hybrid and several natural hybrids. The genus has a discontinuous distribution in mainland China, Taiwan Island, and the Ryukyu islands, providing a model system to explore mechanisms of species dispersal in the East Asian flora. However, the current phylogenetic relationships within Abelia remain uncertain.

MethodsIn this study, we reconstructed phylogenetic relationships within Abelia using nuclear loci generated by target enrichment and plastomes from genome skimming. Divergence time estimation, ancestral area reconstruction, and ecological niche modelling (ENM) were used to examine the diversification history of Abelia.

Key ResultsWe found extensive cytonuclear discordance across the genus. By integrating lines of evidence from molecular phylogenies, divergence times, and morphology, we propose to merge A. macrotera var. zabelioides into A. uniflora. Network analyses suggested that there have been widespread and multiple hybridization events among Abelia species. These hybridization events may have contributed to the speciation mechanism and resulted in a high observed morphological diversity. The diversification of Abelia began in the early Eocene, followed by A. chinensis var. ionandra colonizing the island of Taiwan in the Middle Miocene. The ENM results suggested an expansion of climatically suitable areas during the Last Glacial Maximum and range contraction during the Last Interglacial. Disjunction between the Himalayan-Hengduan Mountain region (HHM) and the island of Taiwan is most likely the consequence of topographic isolation and postglacial contraction.

ConclusionsWe used genomic data to reconstruct the phylogeny of Abelia and found a clear pattern of reticulate evolution in the group. In addition, our results support shrinkage of postglacial range and the heterogeneity of the terrain have led to the disjunction of the mainland China-Taiwan island. This study provides important new insights into the speciation process and taxonomy of Abelia. Identification of CD4+ Sub-population of Resident Cardiac Fibroblasts Linked to Myocardial Fibrosis Infiltration with inflammatory T-cells and accumulation of cardiac myofibroblasts are hallmarks of cardiac fibrosis and maladaptive remodeling. The origin, identity, and functions of the resident cardiac cells involved in this process are, however, unclear. To determine the identity of cells contained in regions exhibiting fibrosis, mass cytometry profiling was performed using resident human ventricular cardiac fibroblasts and right ventricle autopsy tissues from individuals diagnosed with pulmonary hypertension and SUGEN/hypoxia rats. Results showed that a subpopulation of resident myocardial fibroblasts expresses increased levels of CD4+, a helper T-cell surface marker, in addition to mesenchymal markers in humans and rats. Characterization of the resident cardiac fibroblast subpopulation, both structurally and functionally, using transcriptome and secretome analysis of the secreted cytokines, chemokines, proteins, and metabolites, evidenced that IL-1{beta} induces a phenotypic switch of human cardiac fibroblasts from mesenchymal to CD4+ lymphoidal lineage in vitro. RNA sequencing (RNA-seq) analysis of FACS-sorted CD4-expressing cardiac fibroblasts further revealed that the transcriptome of such IL-1{beta}-induced CD4+ fibroblast population exhibited classical lymphoidal and stem cell-like signatures. Lastly, reversal of cell clustering, phosphorylation of MAPK p38 and NF-{kappa}B p65, and phenotypic switching was achieved with the administration of an IL-1R antagonist. In conclusion, we have identified a subpopulation of cardiac fibroblasts which exhibits structural and functional attributes of both mesenchymal and lymphoid cells which is induced by IL-1{beta}-IL-1R-NFkB pathway for differentiation of cardiac fibroblast cells. These data suggest that cardiac fibroblast transdifferentiation during inflammation may form the basis for maladaptive remodeling during myocardial fibrosis. LCK Regulates Homologous Recombination DNA Repair Identifying a New Therapeutic Target for Sensitizing PARP Inhibitors in HR Proficient Ovarian Cancer Poly-ADP Ribose Polymerase (PARP) targeted therapy is clinically approved for the treatment of homologous recombination (HR) repair deficient tumors. The remarkable success in treatment of HR repair deficient cancers has not translated to HR-proficient cancer. Our studies identify the mechanism of non-receptor lymphocyte-specific protein tyrosine kinase (LCK) in HR repair in endometrioid epithelial ovarian cancer (eEOC). LCK expression is induced and activation in the nucleus in response to DNA damage insult. LCK inhibition attenuates expression of RAD51, BRCA1, and BRCA2 proteins necessary for HR-mediated DNA repair, sufficient to suppress RAD51 foci formation, and augments {gamma}H2AX foci formation. Mechanistically, DNA damage leads to direct interaction of LCK with RAD51 and BRCA1 in a kinase dependent manner. Attenuation of LCK sensitized HR-proficient eEOC cells to PARP inhibitor in cell culture and pre-clinical mouse studies. These findings identify a new mechanism for expanding utility of PARP inhibitors in HR proficient ovarian cancer.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=120 SRC="FIGDIR/small/433791v2_ufig1.gif" ALT="Figure 1">
View larger version (33K):
org.highwire.dtl.DTLVardef@1715c60org.highwire.dtl.DTLVardef@167cb5forg.highwire.dtl.DTLVardef@158d807org.highwire.dtl.DTLVardef@1ced201_HPS_FORMAT_FIGEXP  M_FIG C_FIG Human cortical areas are sparsely connected: Combining histology with diffusion MRI to estimate the absolute number of axons The tracts between cortical areas are conceived as playing a central role in cortical information processing, but their actual numbers have never been determined in humans. Here we estimate the absolute number of axons linking cortical areas from a whole-cortex diffusion-MRI (dMRI) connectome, calibrated using the histologically-measured callosal fiber density. Median connectivity is estimated as ~6200 axons between cortical areas within-hemisphere and ~1300 axons inter-hemispherically, with axons connecting functionally-related areas surprisingly sparse. For example, we estimate that <5% of the axons in the trunk of the arcuate and superior longitudinal fasciculi connect Wernickes and Brocas areas. These results suggest that detailed information is transmitted between cortical areas either via linkage of the dense local connections or via rare, extraordinarily privileged long-range connections. Multiomic analisys reveals that viroid infection induces a temporal reprograming of plant-defence mechanisms at multiple regulatory levels Viroids are circular RNAs of minimal complexity compelled to subvert plant-regulatory networks to accomplish their infectious process. Studies focused on the response to viroid infection have mostly addressed specific regulatory levels and considered a unique infection time. Thus, much remains to be done to understand the temporal evolution and complex nature of viroid-host interactions. Here we present an integrative analysis of the temporal evolution and intensity of the genome-wide alterations in cucumber plants infected with hop stunt viroid (HSVd) by integrating differential host transcriptome, sRNAnome and methylome. Our results support that HSVd promotes the redesign of the cucumber regulatory-pathways predominantly affecting specific regulatory layers at different infection-phases. The initial response was characterized by a reconfiguration of the host-transcriptome by differential exon usage, followed by a progressive transcriptional down-regulation modulated by epigenetic changes. Regarding endogenous small RNAs, the alterations were limited and mainly occur at the late stage. The most significant host-alterations were predominantly related to the down-regulation of transcripts involved in plant-defence mechanisms, the restriction of pathogen-movement and the systemic spreading of defence signals. Altogether our data evidence the existence of a dynamic and yet poorly known arms race between the host and the viroid. We expect that these data constituting the first comprehensive map of the plant responses to a viroid infection contribute to elucidate the molecular basis of this multifaceted defence and counter-defence layout. 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. The influence of saccades on insect flight control In a way analogous to human vision, the fruit fly Drosophila melanogaster and many other flying insects generate smooth and saccadic movements to stabilize and shift their gaze in flight, respectively. It has been hypothesized that this combination of continuous and discrete movements benefits both flight stability and performance, particularly at high frequencies or speeds. Here we develop a hybrid control system model to explore the effects of saccades on the yaw stabilization reflex of D. melanogaster. Inspired from experimental data, the model includes a first order plant, a Proportional-Integral (PI) continuous controller, and a saccadic reset system that fires based on the integrated error of the continuous controller. We explore the gain, delay and switching threshold parameter space to quantify the optimum regions for stability and performance. We show that the addition of saccades to a continuous controller provides benefits to both stability and performance across a range of frequencies. Our model suggests that Drosophila operate near its optimal switching threshold for its experimental gain set. We also show that based on experimental data, D. melanogaster operates in a region that trades off performance and stability. This trade-off increases flight robustness to compensate for environmental uncertainties such as wing damage. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics. Meiotic DSB-independent role of protein phosphatase 4 in Hop1 assembly to promote meiotic chromosome axis formation in budding yeast Dynamic changes in chromosomal structure that occur during meiotic prophase play an important role in the progression of meiosis. Among these structures, meiosis-specific chromosomal axis-loop structures are important as a scaffold for integrated control between the meiotic recombination reaction and the associated checkpoint system to ensure accurate chromosome segregation. However, the molecular mechanism of the initial step of chromosome axis-loop construction is not well understood. Here, we showed that in budding yeast, protein phosphatase 4 (PP4) that primarily counteracts Mec1/Tel1 phosphorylation is required to promote the assembly of chromosomal axis components Hop1 and Red1 onto meiotic chromatin via interaction with Hop1. PP4, in contrast, did not affect Rec8 assembly. Notably, unlike the previously known function of PP4, this novel function of PP4 was independent of meiotic DSB-dependent Tel1/Mec1 kinase activities. The defect in Hop1/Red1 assembly in the absence of PP4 function was not suppressed by dysfunction of Pch2, which removes Hop1 protein from the chromosome axis. This suggests that PP4 is required for the initial step of chromatin loading of Hop1 rather than stabilization of Hop1 on axes. These results indicate that phosphorylation/dephosphorylation-mediated regulation of Hop1 recruitment onto chromatin occurs during chromosome axis construction before meiotic double-strand break formation. Flexible control of Pavlovian-instrumental transfer based on expected reward value The Pavlovian-instrumental transfer (PIT) paradigm is widely used to assay the motivational influence of reward-predictive cues, reflected by their ability to invigorate instrumental behavior. Leading theories of incentive learning assume that cues are assigned motivational properties based on the value of the rewards they predict. We outline an alternative account of PIT which recognizes that reward-predictive cues may suppress rather than motivate instrumental behavior under certain conditions. We posit that cues capable of triggering a strong reward expectancy will tend to inhibit instrumental behavior, which is exploratory by nature, in order to facilitate retrieval of the expected reward before it is lost or stolen. According to this view, the motivation to engage in instrumental behavior during a cue should be inversely related to the value of the predicted reward, since there is more to lose by failing to secure a high-value reward than a low-value reward. We tested this hypothesis in rats using a high-conflict PIT protocol which elicits competing instrumental and food-port entry responses. In Experiment 1, cues signaling different reward magnitudes (1, 3, or 9 pellets) differentially influenced these behaviors. Whereas the 1-pellet cue increased instrumental behavior, cues signaling 3 or 9 pellets suppressed instrumental behavior and elicited high levels of food-port activity. Experiment 2 found that reward-predictive cues suppressed instrumental behavior and increased food-port activity in a flexible manner that was disrupted by post-training reward devaluation. We discuss how the PIT task provides a useful tool for studying cognitive control over cue-motivated behavior in rodents. Synergistic block of SARS-CoV-2 infection by combined drug inhibition of the host entry factors PIKfyve kinase and TMPRSS2 protease Repurposing FDA-approved inhibitors able to prevent infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) could provide a rapid path to establish new therapeutic options to mitigate the effects of coronavirus disease 2019 (COVID-19). Proteolytic cleavages of the spike S protein of SARS-CoV-2, mediated by the host cell proteases cathepsin and TMPRSS2, alone or in combination, are key early activation steps required for efficient infection. The PIKfyve kinase inhibitor apilimod interferes with late endosomal viral traffic, and through an ill-defined mechanism prevents in vitro infection through late endosomes mediated by cathepsin. Similarly, inhibition of TMPRSS2 protease activity by camostat mesylate or nafamostat mesylate prevents infection mediated by the TMPRSS2-dependent and cathepsin-independent pathway. Here, we combined the use of apilimod with camostat mesylate or nafamostat mesylate and found an unexpected ~5-10-fold increase in their effectiveness to prevent SARS-CoV-2 infection in different cell types. Comparable synergism was observed using both, a chimeric vesicular stomatitis virus (VSV) containing S of SARS-CoV-2 (VSV-SARS-CoV-2) and SARS-CoV-2 virus. The substantial ~5-fold or more decrease of half maximal effective concentrations (EC50 values) suggests a plausible treatment strategy based on the combined use of these inhibitors.

IMPORTANCEInfection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is causing the coronavirus disease 2019 (COVID-2019) global pandemic. There are ongoing efforts to uncover effective antiviral agents that could mitigate the severity of the disease by controlling the ensuing viral replication. Promising candidates include small molecules that inhibit the enzymatic activities of host proteins, thus preventing SARS-CoV-2 entry and infection. They include Apilimod, an inhibitor of PIKfyve kinase and camostat mesylate and nafamostat mesylate, inhibitors of TMPRSS2 protease. Our research is significant for having uncovered an unexpected synergism in the effective inhibitory activity of apilimod used together with camostat mesylate or with nafamostat mesylate. Feline tooth resorption findings from a questionnaire survey on 8115 Finnish cats and a case-control study for a more defined population drawn from the survey Tooth resorption (TR) is one of the most common dental diseases in cats. It is a painful disease characterised by progressive dental destruction, which eventually results in loss of teeth. The aetiology of the disease remains unclear, but associations with old age, breed, other oral and dental diseases, and certain environmental factors have been suspected. We collected health and environmental information of 8115 Finnish cats from 41 breeds through an online survey for cat owners to study TR in Finland with an observational study. The frequency of TR that the owner reported as veterinarian-diagnosed was 3.9% in the entire data (316/8115) and 15% in cats that the owner reported as having veterinarian-diagnosed oral or dental disease (316/2070). We determined the risk factors and breed variation of feline TR by logistic regression with the entire data and a case-control study in a more defined population drawn from the entire data. In the entire data, we defined those cats that the owners had reported as being diagnosed with TR by a veterinarian as TR-cats and compared them to all other cats. In the case-control study, we included only cats with owner-reported veterinarian-diagnosed oral and dental diseases. We defined cats with TR-diagnosis made under sedation as TR-cases. Other veterinarian-diagnosed TR-cats were excluded. Cats with other oral and dental diseases were used as controls. The model with the entire population included 308 TR-cats compared with 7508 cats, and the case-control study included 198 TR-cats and 1683 controls. Although the diagnoses were reportedly made by a veterinarian, the owner answered the survey and might have understood, remembered, or interpreted the diagnoses in a way that could affect the results. The results indicated an increasing risk of TR with age. Dental calculus, gingivitis, and periodontitis were potential risk factors for TR. These findings and the interaction between dental calculus and gingivitis suggest that inflammatory changes caused by dental calculus increase the risk of TR. Keeping food available constantly was a potential protective factor. We found that Cornish Rex, European, and Ragdoll breeds were at higher risk for TR. The observed differences between breeds highlight a genetic contribution. GeneCup: mine PubMed for gene relationships using custom ontology and deep learning Interpreting and integrating results from omics studies typically requires a comprehensive and time consuming survey of extant literature. Here, we introduce GeneCup, an easy to use literature mining web service that searches all PubMed abstracts for user-provided gene symbols in conjunction with a set of custom keywords organized into a customized ontology, as well as results from human genome-wide association studies (GWAS). As an example, we organized over 300 keywords related to drug addiction into seven categories. The literature search is conducted by querying the NIH PubMed server using a programming interface, which is followed by retrieving abstracts from a local copy of the PubMed archive. The main results presented to the user are individual sentences containing the gene symbol, organized by the keywords they also contain. These sentences are presented through an interactive graphical interface or as tables. GWAS results are displayed using a similar method. All results are linked to the original abstract in PubMed. In addition, a convolutional neural network is employed to distinguish sentences describing systemic stress from those describing cellular stress. The automated and comprehensive search strategy provided by GeneCup facilitates the integration of new discoveries from omic studies with existing literature. GeneCup is free and open source software. The source code of GeneCup and the link to a running instance is available at https://github.com/hakangunturkun/GeneCup 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 Quantifying the Morphology and Mechanisms of Cancer Progression in 3D in-vitro environments: Integrating Experiments and Multiscale Models Mathematical models of cancer growth have become increasingly more accurate both in the space and time domains. However, the limited amount of data typically available has resulted in a larger number of qualitative rather than quantitative studies. In the present study, we provide an integrated experimental-computational framework for the quantification of the morphological characteristics and the mechanistic modelling of cancer progression in 3D environments. The proposed framework allows for the calibration of multiscale, spatiotemporal models of cancer growth using state-of-the-art 3D cell culture data, and their validation based on the resulting experimental morphological patterns. Its implementation enabled us to pursue two goals; first, the quantitative description of the morphology of cancer progression in 3D cultures, and second, the relation of tumour morphology with underlying biophysical mechanisms that govern cancer growth and migration. We applied this framework to the study of the spatiotemporal progression of Triple Negative Breast Cancer cells cultured in Matrigel scaffolds, and validated the hypothesis of chemotactic migration using a multiscale Keller-Segel model. The results revealed transient, non-random spatial distributions of cancer cells that consist of clustered, and dispersion patterns. The proposed model was able to describe the general characteristics of the experimental observations and suggests that cancer cells exhibited chemotactic migration and accumulation, as well as random motion during the examined time period of development. To our knowledge, this is the first time that a multiscale model is used to quantify the relationship between the spatial patterns and the underlying mechanisms of cancer growth in 3D environments. A phenomics approach for in vitro antiviral drug discovery BackgroundThe current COVID-19 pandemic has highlighted the need for new and fast methods to identify novel or repurposed therapeutic drugs. Here we present a method for untargeted phenotypic drug screening of virus-infected cells, combining Cell Painting with antibody-based detection of viral infection in a single assay. We designed an image analysis pipeline for segmentation and classification of virus-infected and non-infected cells, followed by extraction of morphological properties.

ResultsWe show that the methodology can successfully capture virus-induced phenotypic signatures of MRC-5 human lung fibroblasts infected with Human coronavirus 229E (CoV-229E). Moreover, we demonstrate that our method can be used in phenotypic drug screening using a panel of nine host- and virus-targeting antivirals. Treatment with effective antiviral compounds reversed the morphological profile of the host cells towards a non-infected state.

ConclusionsThe method can be used in drug discovery for morphological profiling of novel antiviral compounds on both infected and non-infected cells. Assigning metabolic rate measurements to torpor and euthermy in heterothermic endotherms: 'torpor', a new package for R. Torpor is a state of controlled reduction of metabolic rate (M) in endotherms. Assigning measurements of M to torpor or euthermy can be challenging, especially when the difference between euthermic M and torpid M is small, in species defending a high minimal body temperature in torpor, in thermolabile species, and slightly below the thermoneutral zone (TNZ). Here, we propose a novel method for distinguishing torpor from euthermy. We use the variation in M measured during euthermic rest and torpor at varying ambient temperatures (Ta) to objectively estimate the lower critical temperature (Tlc) of the TNZ and to assign measurements to torpor, euthermic rest or rest within TNZ. In addition, this method allows the prediction of M during euthermic rest and torpor at varying Ta, including resting M within the TNZ. The present method has shown highly satisfactory results using 28 published sets of metabolic data obtained by respirometry on 26 species of mammals. Ultimately, this novel method aims to facilitate analysis of respirometry data in heterothermic endotherms. Finally, the development of the associated R-package (torpor) will enable widespread use of the method amongst biologists.

Summary statementThe presented method and its associated R-package (torpor) enable the assignment of metabolic rate measurements to torpor or euthermy, ultimately improving the standardization of respirometry analyses in heterotherms. 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. "Adopt a Microorganism" methodology enhances the discourse richness of high school students in both hybrid learning and remote education The Internet has changed the way teachers and students access information and build knowledge. The recent COVID-19 pandemic has created challenges for both teachers and students and a demand for new methodologies of remote learning. In the life sciences, mixing online content with practical activities represents an even greater challenge. In microbiology, the implementation of an active teaching methodology, the #Adopt project, based on the social network Facebook(R), represents an excellent option for connecting remote education with classroom activities. In 2020, the version applied in high school, "Adopt a Microorganism", was adapted to meet the demands of emergency remote education owing to the suppression of face-to-face activities caused by the pandemic. In the present study, we assessed how the change in methodology impacted the discourse richness of students from high school integrated with technical education in the Business Administration program of the Federal Institute of Sao Paulo, Sorocaba Campus. Three questionnaires related to the groups of microorganisms (Archaea, Bacteria, Virus, Fungi, and Protozoan) were applied. The students responses in the 2019 and 2020 classes were compared concerning content richness and multiplicity of concepts through the application of the Shannon diversity index, an approach that is generally used to assess biodiversity in different environments. The observed results suggest that remote learning provided students with a conceptual basis and richness of content equivalent to that achieved by students subjected to the hybrid teaching model. In conclusion, this study suggests that the #Adopt project methodology increases students discourse richness in microbiology even without face-to-face traditional classes. Inflammatory response in hematopoietic stem and progenitor cells triggered by activating SHP2 mutations evokes blood defects Gain-of-function mutations in the protein-tyrosine phosphatase SHP2 are the most frequently occurring mutations in sporadic juvenile myelomonocytic leukemia (JMML) and JMML-like myeloproliferative neoplasm (MPN) associated with Noonan syndrome (NS). Hematopoietic stem and progenitor cells (HSPCs) are the disease propagating cells of JMML. Here, we explored transcriptomes of HSPCs with Shp2 mutations derived from JMML patients and a novel NS zebrafish model. In addition to major NS traits, CRISPR/Cas9 knock-in Shp2D61G mutant zebrafish recapitulated a JMML-like MPN phenotype, including myeloid lineage hyperproliferation, ex vivo growth of myeloid colonies and in vivo transplantability of HSPCs. Single cell mRNA sequencing of HSPCs from Shp2D61G zebrafish embryos and bulk sequencing of HSPCs from JMML patients revealed an overlapping inflammatory gene expression pattern. Strikingly, an anti-inflammatory agent rescued JMML-like MPN in Shp2D61G zebrafish embryos. Our results indicate that a common inflammatory response was triggered in the HSPCs from sporadic JMML patients and syndromic NS zebrafish, which potentiated MPN and may represent a future target for JMML therapies. Identification of a conserved neutralizing epitope present on spike proteins from all highly pathogenic coronaviruses Three pathogenic human coronaviruses have emerged within the last 20 years, with SARS-CoV-2 causing a global pandemic. Although therapeutic antibodies targeting the SARS-CoV-2 spike currently focus on the poorly conserved receptor-binding domain, targeting essential neutralizing epitopes on the more conserved S2 domain may provide broader protection. We report an antibody binding an epitope conserved in the pre-fusion core of MERS-CoV, SARS-CoV and SARS-CoV-2 spike S2 domains. Antibody 3A3 binds a conformational epitope with ~2.5 nM affinity and neutralizes spike from SARS-CoV, SARS-CoV-2 and variants of concern in in vitro pseudovirus assays. Hydrogen-deuterium exchange mass spectrometry identified residues 980-1006 in the flexible hinge region at the S2 apex as the 3A3 epitope, suggesting 3A3 prevents the S2 conformational rearrangements required for conversion to the spike post-fusion state and virus-host cell fusion. This work defines a conserved vulnerable site on the SARS-CoV-2 S2 domain and guides the design of pan-protective spike immunogens. Diversity of SARS-CoV-2 genome among various strains identified in Lucknow, Uttar Pradesh Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has emerged as a significant challenge worldwide. Rapid genome sequencing of SARS-CoV-2 is going on across the globe to detect mutations and genomic modifications in SARS-CoV-2. In this study, we have sequenced twenty-three SARS-CoV-2 positive samples collected during the first pandemic from the state of Uttar Pradesh, India. We observed a range of already reported mutations (2-22), including; D614G, L452R, Q613H, Q677H, T1027I in the S gene; S194L in the N gene; Q57H, L106F, T175I in the ORF3. Few unreported mutations such as P309S in the ORF1ab gene; T379I in the N gene; and L52F, V77I in the ORF3a gene were also detected. Phylogenetic genome analysis showed similarity with other SARS-CoV-2 viruses reported from Uttar Pradesh. The observed mutations may be associated with SARS-CoV-2 virus pathogenicity or disease severity. No impact of biocontrol agent's predation cues on development time or size of surviving Aedes albopictus under optimal nutritional availability Cyclopoid copepods have been applied successfully to limit populations of highly invasive Aedes albopictus mosquitoes that can transmit diseases of public health importance. However, there is concern that changes in certain mosquito traits, induced by exposure to copepod predation, might increase the risk of disease transmission. In this study, third instar Ae. albopictus larvae, hereafter referred to as "focal individuals", were exposed to Megacyclops viridis predator cues associated with both consumption of newly-hatched mosquito larvae and attacks on focal individuals. The number of newly-hatched larvae surrounding each focal larva was held constant to control for density effects on size, and the focal individuals day of pupation and wing length were recorded for each replicate.

Exposing late instar Ae. albopictus to predation decreased their chances of surviving to adulthood, and three focal larvae that died in the predator treatment showed signs of melanization, indicative of wounding. Among surviving focal Ae. albopictus, no significant difference in either pupation day or wing length was observed due to copepod predation. The absence of a significant sublethal impact from M. viridis copepod predation on surviving later-stage larvae in this analysis supports the use of M. viridis as a biocontrol agent against Ae. albopictus.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=106 SRC="FIGDIR/small/460697v2_ufig1.gif" ALT="Figure 1">
View larger version (9K):
org.highwire.dtl.DTLVardef@1abaabborg.highwire.dtl.DTLVardef@1561651org.highwire.dtl.DTLVardef@d325a7org.highwire.dtl.DTLVardef@4d1830_HPS_FORMAT_FIGEXP  M_FIG C_FIG 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.