Comparative analysis of molecular representations in prediction of drug combination effects AO_SCPLOWBSTRACTC_SCPLOWApplication of machine and deep learning methods in drug discovery and cancer research has gained a considerable amount of attention in the past years. As the field grows, it becomes crucial to systematically evaluate the performance of novel computational solutions in relation to established techniques. To this end we compare rule-based and data-driven molecular representations in prediction of drug combination sensitivity and drug synergy scores using standardized results of 14 throughput screening studies, comprising 64 200 unique combinations of 4 153 molecules tested in 112 cancer cell lines. We evaluate the clustering performance of molecular representations and quantify their similarity by adapting the Centered Kernel Alignment metric. Our work demonstrates that to identify an optimal molecular representation type it is necessary to supplement quantitative benchmark results with qualitative considerations, such as model interpretability and robustness, which may vary between and throughout preclinical drug development projects.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=168 SRC="FIGDIR/small/439299v2_ufig1.gif" ALT="Figure 1">
View larger version (41K):
org.highwire.dtl.DTLVardef@1bc4509org.highwire.dtl.DTLVardef@1586640org.highwire.dtl.DTLVardef@a11129org.highwire.dtl.DTLVardef@6dc93b_HPS_FORMAT_FIGEXP  M_FIG C_FIG
Onlinemeta: A Web Server For Meta-Analysis Based On R-shiny Meta-analysis is a common statistical method used to summarize multiple studies that cover the same topic. It can provide less biased results and explain heterogeneity between studies. Although there exists a variety of meta-analysis software, they are rarely both convenient to use and capable of comprehensive analytical functions. As a result, we established a meta-analysis web tool called Meta-Analysis Online (Onlinemeta). Onlinemeta includes both risk bias analysis and meta-analysis. The risk bias analysis tool can produce heatmaps and histograms, whereas the meta-analysis tool can be used to analyze dichotomous variables, single-armed dichotomous variables, continuous variables, single-armed continuous variables, survival data, and diagnostic experiments. In addition, it can be used to generate forest plots, funnel plots, SROC curve, and crosshair plots. Onlinemeta can be found at the following address: https://smuonco.Shinyapps.io/Onlinemeta/.
A chemical screen based on an interruption of zebrafish gastrulation identifies the HTR2C inhibitor Pizotifen as a suppressor of EMT-mediated metastasis Metastasis is responsible for approximately 90% of cancer-associated mortality but few models exist that allow for rapid and effective screening of anti-metastasis drugs. Current mouse models of metastasis are too expensive and time consuming to use for rapid and high-throughput screening. Therefore, we created a unique screening concept utilizing conserved mechanisms between zebrafish gastrulation and cancer metastasis for identification of potential anti-metastatic drugs. We hypothesized that small chemicals that interrupt zebrafish gastrulation might also suppress metastatic progression of cancer cells and developed a phenotype-based chemical screen to test the hypothesis. The screen used epiboly, the first morphogenetic movement in gastrulation, as a marker and enabled 100 chemicals to be tested in five hours. The screen tested 1280 FDA-approved drugs and identified Pizotifen, an antagonist for serotonin receptor 2C (HTR2C) as an epiboly-interrupting drug. Pharmacologic and genetic inhibition of HTR2C suppressed metastatic progression in a mouse model. Blocking HTR2C with Pizotifen restored epithelial properties to metastatic cells through inhibition of Wnt-signaling. In contrast, HTR2C induced epithelial to mesenchymal transition (EMT) through activation of Wnt-signaling and promoted metastatic dissemination of human cancer cells in a zebrafish xenotransplantation model. Taken together, our concept offers a novel platform for discovery of anti-metastasis drugs. Genomic databanks and targeted assays help characterise domestic mosquito incursions Biosecurity strategies aiming to restrict the spread of invasive pests can benefit from knowing where new incursions have come from. This knowledge can be acquired using genomic databanks, where genetic variation in incursion samples is compared against variation in reference samples. Here, we use genomic databanks to investigate domestic incursions of two mosquito species in Australia. We used a deep learning method to trace a 2021 invasion of Aedes aegypti in Tennant Creek, Northern Territory, to Townsville, Queensland, and to trace Ae. albopictus incursions between several invaded and uninvaded islands in the Torres Strait. We observed high precision of tracing despite 30-70 generations separating incursions and source samples. Targeted assays also improved precision, in this case by comparing Wolbachia infection data to mitochondrial DNA variation in Ae. aegypti. For Ae. albopictus, incursions detected in 2021 were traced to different islands than incursions detected in 2019. Patterns of relatedness and inbreeding indicated that Tennant Creek was likely invaded by Ae. aegypti from one family, whereas Torres Strait incursions involved distinct kinship groups. Our results highlight the value of repeated tracing with genomic databanks that remain informative over years, and demonstrate how additional targeted assays (e.g. Wolbachia) can improve inferences. Blebbistatin protects iPSC-CMs from hypercontraction and facilitates automated patch-clamp based electrophysiological study Recently, there have been great advances in cardiovascular channelopathy modeling and drug safety pharmacology using human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). The automated patch-clamp (APC) technique overcomes the disadvantages of the manual patch-clamp (MPC) technique, which is labor intensive and gives low output. However, the application of the APC platform is still limited in iPSC-CM based research, due to the difficulty in maintaining the high quality of single iPSC-CMs during dissociation and recording. In this study, we improved the method for single iPSC-CM preparation by applying 2.5 M blebbistatin (BB, an excitation-contraction coupling uncoupler) throughout APC procedures (dissociation, filtration, storage, and recording). Under non-BB buffered condition, iPSC-CMs in suspension showed a severe bleb-like morphology. However, BB-supplement led to significant improvements in morphology and INa recording, and we even obtained several CMs that showed spontaneous action potentials with typical morphology. Furthermore, APC faithfully recapitulated the single-cell electrophysiological phenotypes of iPSC-CMs derived from Brugada syndrome patients, as detected with MPC. Our study indicates that APC is capable of replacing MPC in the modeling of cardiac channelopathies using human iPSC-CMs by providing high-quality data with higher throughput. 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. Echolocation reverses information flow in a cortical vocalization network The mammalian frontal and auditory cortices are important for vocal behaviour. Here, using local field potential recordings, we demonstrate for the first time that the timing and spatial pattern of oscillations in the fronto-auditory cortical network of vocalizing bats (Carollia perspicillata) predict the purpose of vocalization: echolocation or communication. Transfer entropy analyses revealed predominantly top-down (frontal-to-auditory cortex) information flow during spontaneous activity and pre-vocal periods. The dynamics of information flow depended on the behavioural role of the vocalization and on the timing relative to vocal onset. Remarkably, we observed the emergence of predominantly bottom-up (auditory-to-frontal cortex) information transfer patterns specific echolocation production, leading to self-directed acoustic feedback. Electrical stimulation of frontal areas selectively enhanced responses to echolocation sounds in auditory cortex. These results reveal unique changes in information flow across sensory and frontal cortices, potentially driven by the purpose of the vocalization in a highly vocal mammalian model. Inhibitory proteins block substrate access by occupying the active site cleft of Bacillus subtilis intramembrane metalloprotease SpoIVFB Intramembrane proteases function in numerous signaling pathways that impact health, but knowledge about regulation of intramembrane proteolysis is limited. We examined inhibition of intramembrane metalloprotease SpoIVFB by proteins BofA and SpoIVFA. We found that BofA residues in and near a predicted transmembrane segment are required for SpoIVFB inhibition, and cross-linking experiments indicated that this transmembrane segment occupies the SpoIVFB active site cleft. SpoIVFA is also required for SpoIVFB inhibition. The inhibitory proteins block access of the substrate N-terminal Proregion to the membrane-embedded SpoIVFB active site, based on additional cross-linking experiments; however, the inhibitory proteins did not prevent interaction between the substrate C-terminal region and the SpoIVFB soluble domain. A structural model was built of SpoIVFB in complex with BofA and parts of SpoIVFA and substrate, using partial homology and constraints from cross-linking and co-evolutionary analyses. The model predicts that conserved BofA residues interact to stabilize a transmembrane segment and a membrane-embedded C-terminal region. SpoIVFA is predicted to bridge the BofA C-terminal region and SpoIVFB, forming a membrane-embedded inhibition complex. Our results reveal a novel mechanism of intramembrane protease inhibition with clear implications for relief from inhibition in vivo and design of inhibitors as potential therapeutics. If Research is Not the Evidence, What is it? The Lack of Research Citations in Health Vlogs Doctor-influencers have the knowledge needed to understand scientific research, and they have the social media popularity to share it with large audiences. In this article, I explore the treatment of evidence in vlogs by Egyptian doctor-influencers: how they present it and what they use instead of it. I answer the following questions: 1) How do doctor-influencers present evidence on social media? 2) If evidence is not the research, then what is it? The data was collected through in-depth interviews with 12 Egyptian doctor-influencers from different specializations; two focus groups; and a critical discourse analysis of 48 of the influencers most popular and most engaging videos. I found that the doctor-influencers cited academic research only in the videos about medical controversies, new treatments, misinformation, or common medical mistakes. In most of their videos, the Egyptian doctor-influencers scarcely referred to medical research. This is because they believe that their audiences would not understand evidence due to their low educational attainment, weak research skills, and the language barrier (as most medical research is published in English rather than Arabic). These findings highlight the challenges of communicating science in societies with a low level of education. In these settings, communicators may need to move beyond mediating research studies on social media and instead focus on simplifying, translating, and connecting evidence-based information to the cultural context. Single-cell phenotypic characteristics of tolerance under recurring antibiotic exposure in Escherichia coli. Non-heritable drug resistance, such as tolerance and persistence towards antibiotics, is little characterized compared to genetic resistance. Tolerance and persistence allow cells to survive application of antibiotics that are bactericidal to non-tolerant cells. Non-heritable drug resistance challenges antibiotic treatment, particularly of recurrent infections, and have implications towards heritable resistance evolution. Tolerant cells have commonly been characterized as growth arrested cells prior and during antibiotic application that quickly resume growth post-application. Here, we explore characteristic of tolerant and susceptible E. coli single bacteria cells to different levels of recurrent antibiotic exposure and quantify their occurrence. In using a high throughput single-cell microfluidic device, we find that tolerant cells reduce their growth rate by about 50%, but contrary to previous findings do not go into growth arrest or near growth arrest. The growth reduction is induced by antibiotic exposure and not caused by a stochastic switch or predetermined state as previously described. Cells exhibiting constant intermediate growth survived best under antibiotic exposure and selection did not primarily act on fast growing cells, as expected for a {beta}-lactam antibiotic. Control experiments on population cultures confirmed and challenged scaling of single cell findings to population level processes. Our findings suggest a prevalent type of tolerance that differs from previously described tolerance and persister characteristics. Our described characteristics and its high frequency of occurrence supports acclaims of an underappreciated role of tolerant cells towards resistance evolution.

Author SummaryResistance to antibiotics is a major cause of treatment failure in modern medicine. In addition to genetic resistance, bacterial cells can sustain periods of bactericidal application in a so-called tolerant state. Previously tolerant cells were characterized by growth arrest prior and during exposure to antibiotics. Here, we show a distinctive type of tolerance that arises out of normally growing cells, where during antibiotic exposure, cell growth is slowed but does not go into growth arrest. Cells with an intermediate growth have the best survival, which is contrary to the assumptions of exposure to a {beta}-lactam antibiotic where selection is assumed to act against the fastest growing cells and the least against non-growing cells. The novel characteristics of tolerance that we have uncovered provides bacteria with another route to evade antibiotics and that could be a considerable threat to the efficiency of antibiotic treatment. The OPAQUE1/DISCORDIA2 myosin XI is required for phragmoplast guidance during asymmetric cell division in maize Asymmetric divisions produce cells with different fates and are critical for development. Here we show that a maize myosin XI protein, OPAQUE1 (O1), is necessary for asymmetric divisions during maize stomatal development. We analyzed stomatal precursor cells prior to and during asymmetric division to determine why o1 mutants have abnormal division planes. Cell polarization and nuclear positioning occur normally in the o1 mutant, and the future site of division is correctly specified. The defect in o1 occurs during late cytokinesis, when the plant-specific phragmoplast -- made of microtubules, actin and other proteins - forms the nascent cell plate. The phragmoplast becomes misguided and does not meet the previously established division site. Initial phragmoplast guidance is correct in o1. However, as phragmoplast expansion continues, phragmoplasts in o1 stomatal precursor cells become misguided and do not meet the cortex at the established division site. To understand how this myosin protein contributes to phragmoplast guidance, we identified O1-interacting proteins. Other myosins, specific actin-binding proteins, and maize kinesins related to the Arabidopsis thaliana division site markers PHRAGMOPLAST ORIENTING KINESINs (POKs) interact with O1. We propose that different myosins are important at multiple steps of phragmoplast expansion, and the O1 actin motor and POK-like microtubule motors work together to ensure correct late-stage phragmoplast guidance. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. A humanized nanobody phage display library yields potent binders of SARS CoV-2 spike Neutralizing antibodies targeting the SARS-CoV-2 spike protein have shown a great preventative/therapeutic potential. Here, we report a rapid and efficient strategy for the development and design of SARS-CoV-2 neutralizing humanized nanobody constructs with sub-nanomolar affinities and nanomolar potencies. CryoEM-based structural analysis of the nanobodies in complex with spike revealed two distinct binding modes. The most potent nanobody, RBD-1-2G(NCATS-BL8125), tolerates the N501Y RBD mutation and remains capable of neutralizing the B.1.1.7 (Alpha) variant. Molecular dynamics simulations provide a structural basis for understanding the neutralization process of nanobodies exclusively focused on the spike-ACE2 interface with and without the N501Y mutation on RBD. A primary human airway air-lung interface (ALI) ex vivo model showed that RBD-1-2G-Fc antibody treatment was effective at reducing viral burden following WA1 and B.1.1.7 SARS-CoV-2 infections. Therefore, this presented strategy will serve as a tool to mitigate the threat of emerging SARS-CoV-2 variants. Visualization of lake succession as indicated by macrophyte species composition and trophic levels in lakes subdivided by volcanic eruption at Akan Caldera, Hokkaido, Japan It is not feasible to continuously observe ecological succession in lakes because of the long time-scales generally involved. Thus, the process has been inductively deduced by comparing many lakes with different succession states, or indirectly simulated by tracking studies of smaller water bodies, experiments using microcosms or mesocosms, and reconstruction of lake history by sediment analysis. However, the reality of succession processes in large lakes with slow succession is not well understood, and new approaches are needed. Theoretically, in a group of large and small lakes of similar ages and with similar initial and watershed environments, the rate of nutrient accumulation in each lake depends on the ratio of watershed area to lake size, and the lakes are predicted to evolve to different trophic levels over time. Here, we tested this hypothesis on the 10 lakes of varying sizes in Akan Caldera, Japan, which were formed thousands of years ago by fragmentation due to volcanic eruptions within the caldera. Topographic and water quality studies showed that the ratio of accumulated watershed area to lake size (area and volume), expressed logarithmically, had a positive linear regression with the total phosphorus concentration, an indicator of trophic level. The trophic levels of the lakes were diverse, including oligotrophic, mesotrophic, and eutrophic types in the traditional "lake type" classification based on total phosphorus and chlorophyll-a concentrations. Furthermore, 21 species of aquatic macrophytes were observed by a diving survey, and the plant species composition was classified into five groups corresponding to the trophic status of the lakes, indicating a conventional "hydrarch succession". The diversity of water quality and aquatic vegetation in a group of lakes with similar origins paves the way for new comparative studies of lakes, including large lakes. 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. Rational Exploration of Fold Atlas for Human Solute Carrier Proteins Solute carrier (SLC) superfamily is the largest group responsible for transmembrane transport of substances in human cells. It includes more than 400 members which are organized into 65 families according to their physiological function and sequence similarity. Different families of solute transporters can adopt the same or different folds that determine the working mechanism and reflect the evolutionary relationship between SLC members. Analysis of structural data in literatures shows that there are 13 different folds in the solute carrier superfamily covering 40 families and 343 members. To further study their working mechanism, we systematically explored the SLC superfamily to look for more folds. Based on our results, at least three new folds were found for SLC superfamily, and one was experimentally verified in the SLC44 family. Our work has laid a foundation and provided important insights for the systematic and comprehensive study on the structure and function of solute carriers. Human papillomavirus integration transforms chromatin to drive oncogenesis Human papillomavirus (HPV) drives almost all cervical cancers and up to [~]70% of head and neck cancers. Frequent integration into the host genome occurs only for tumourigenic strains of HPV. We hypothesized that changes in the epigenome and transcriptome contribute to the tumourigenicity of HPV. We found that viral integration events often occurred along with changes in chromatin state and expression of genes near the integration site. We investigated whether introduction of new transcription factor binding sites due to HPV integration could invoke these changes. Some regions within the HPV genome, particularly the position of a conserved CTCF binding site, showed enriched chromatin accessibility signal. ChIP-seq revealed that the conserved CTCF binding site within the HPV genome bound CTCF in 4 HPV+ cancer cell lines. Significant changes in CTCF binding pattern and increases in chromatin accessibility occurred exclusively within 100 kbp of HPV integration sites. The chromatin changes co-occurred with out-sized changes in transcription and alternative splicing of local genes. We analyzed the essentiality of genes upregulated around HPV integration sites of The Cancer Genome Atlas (TCGA) HPV+ tumours. HPV integration upregulated genes which had significantly higher essentiality scores compared to randomly selected upregulated genes from the same tumours. Our results suggest that introduction of a new CTCF binding site due to HPV integration reorganizes chromatin and upregulates genes essential for tumour viability in some HPV+ tumours. These findings emphasize a newly recognized role of HPV integration in oncogenesis. Myoparr-associated and -independent multiple roles of heterogeneous nuclear ribonucleoprotein K during skeletal muscle cell differentiation RNA-binding proteins (RBPs) regulate cell physiology via the formation of ribonucleic-protein complexes with coding and non-coding RNAs. RBPs have multiple functions in the same cells; however, the precise mechanism through which their pleiotropic functions are determined remains unknown. In this study, we revealed the multiple inhibitory functions of heterogeneous nuclear ribonucleoprotein K (hnRNPK) for myogenic differentiation. We first identified hnRNPK as a lncRNA Myoparr binding protein. Gain- and loss-of-function experiments showed that hnRNPK repressed the expression of myogenin at the transcriptional level. The hnRNPK-binding region of Myoparr was required to repress myogenin expression. Moreover, hnRNPK repressed the expression of a set of genes coding for aminoacyl-tRNA synthetases in a Myoparr-independent manner. Mechanistically, hnRNPK regulated the eIF2/Atf4 pathway, one branch of the intrinsic pathways of the endoplasmic reticulum sensors, in differentiating myoblasts. Thus, our findings demonstrate that hnRNPK plays lncRNA-associated and - independent multiple roles during myogenic differentiation, indicating that the analysis of lncRNA-binding proteins will be useful for elucidating both the physiological functions of lncRNAs and the multiple functions of RBPs. Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. Heterogeneous recruitment abilities to RNA polymerases generate nonlinear scaling of gene expression with cell volume While most genes expression levels are proportional to cell volumes, some genes exhibit nonlinear scaling between their expression levels and cell volume. Therefore, their mRNA and protein concentrations change as the cell volume increases, which often have crucial biological functions such as cell-cycle regulation. However, the biophysical mechanism underlying the nonlinear scaling between gene expression and cell volume is still unclear. In this work, we show that the nonlinear scaling is a direct consequence of the heterogeneous recruitment abilities of promoters to RNA polymerases based on a gene expression model at the whole-cell level. Those genes with weaker (stronger) recruitment abilities than the average ability spontaneously exhibit superlinear (sublinear) scaling with cell volume. Analysis of the promoter sequences and the nonlinear scaling of Saccharomyces cerevisiaes mRNA levels shows that motifs associated with transcription regulation are indeed enriched in genes exhibiting nonlinear scaling, in concert with our model. 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. Characterising ontogeny of numerosity discrimination in zebrafish reveals the use of multiple, numerical and non- numerical mechanisms. A sense of non-symbolic numerical magnitudes is widespread in the animal kingdom and has been documented in adult zebrafish. Here we investigated the ontogeny of this ability using a group size preference task in juvenile zebrafish. Fish showed group size preference from 21 days post fertilization (dpf) and reliably chose the larger group when presented with discriminations of between 1 vs. 3, 2 vs. 5 and 2 vs. 3 conspecifics but not 2 vs. 4 conspecifics. When the ratio between the number of conspecifics in each group was maintained at 1:2, fish could discriminate between 1 vs. 2 individuals and 3 vs. 6, but again, not when given a choice between 2 vs. 4 individuals. These findings are in agreement with studies in other species suggesting the systems involved in quantity representation do not operate separately from other cognitive mechanisms. Rather they suggest quantity processing in fish may be the result of an interplay between attentional, cognitive and memory-related mechanisms as in humans and other animals. Our results emphasise the potential of the use of zebrafish to explore the genetic and neural processes underlying the ontogeny and function of number cognition. Divergent Acyl Carrier Protein Decouples Mitochondrial Fe-S Cluster Biogenesis from Fatty Acid Synthesis in Malaria Parasites Plasmodium falciparum malaria parasites are early-diverging eukaryotes with many unusual metabolic adaptations. Understanding these adaptations will give insight into parasite evolution and unveil new parasite-specific drug targets. Most eukaryotic cells retain a mitochondrial fatty acid synthesis (FASII) pathway whose acyl carrier protein (mACP) and 4-phosphopantetheine (Ppant) prosthetic group provide a soluble scaffold for acyl chain synthesis. In yeast and humans, mACP also functions to biochemically couple FASII activity to electron transport chain (ETC) assembly and Fe-S cluster biogenesis. In contrast to most eukaryotes, the Plasmodium mitochondrion lacks FASII enzymes yet curiously retains a divergent mACP lacking a Ppant group. We report that ligand-dependent knockdown of mACP is lethal to parasites, indicating an essential FASII-independent function. Decyl-ubiquinone rescues parasites temporarily from death, suggesting a dominant dysfunction of the mitochondrial ETC followed by broader cellular defects. Biochemical studies reveal that Plasmodium mACP binds and stabilizes the Isd11-Nfs1 complex required for Fe-S cluster biosynthesis, despite lacking the Ppant group required for this association in other eukaryotes, and knockdown of parasite mACP causes loss of both Nfs1 and the Rieske Fe-S protein in ETC Complex III. This work reveals that Plasmodium parasites have evolved to decouple mitochondrial Fe-S cluster biogenesis from FASII activity, and this adaptation is a shared metabolic feature of other Apicomplexan pathogens, including Toxoplasma and Babesia. This discovery also highlights the ancient, fundamental role of ACP in mitochondrial Fe-S cluster biogenesis and unveils an evolutionary driving force to retain this interaction with ACP independent of its eponymous function in FASII.

Significance StatementPlasmodium malaria parasites are single-celled eukaryotes that evolved unusual metabolic adaptations. Parasites require a mitochondrion for blood-stage viability, but essential functions beyond the electron transport chain are sparsely understood. Unlike yeast and human cells, the Plasmodium mitochondrion lacks fatty acid synthesis enzymes but retains a divergent acyl carrier protein (mACP) incapable of tethering acyl groups. Nevertheless, mACP is essential for parasite viability by binding and stabilizing the core mitochondrial Fe-S cluster biogenesis complex via a divergent molecular interface lacking an acyl-pantetheine group that contrasts with other eukaryotes. This discovery unveils an essential metabolic adaptation in Plasmodium and other human parasites that decouples mitochondrial Fe-S cluster biogenesis from fatty acid synthesis and evolved at or near the emergence of Apicomplexan parasitism. A new method based on the von Mises-Fisher distribution shows that a minority of liver-localized CD8 T cells display hard-to-detect attraction to Plasmodium-infected hepatocytes Malaria is a disease caused by Plasmodium parasites, resulting in over 200 million infections and 400,000 deaths every year. A critical step of malaria infection is when sporozoites, injected by mosquitoes, travel to the liver and form liver stages. Malaria vaccine candidates which induce large numbers of malaria-specific CD8 T cells in mice are able to eliminate all liver stages, preventing fulminant malaria. However, how CD8 T cells find all parasites in 48 hours of the liver stage lifespan is not well understood. Using intravital microscopy of murine livers, we generated unique data on T cell search for malaria liver stages within a few hours after infection. To detect attraction of T cells to an infection site, we used the von Mises-Fisher distribution in 3D, similar to the 2D von Mises distribution previously used in ecology. Our results suggest that the vast majority (70-95%) of malaria-specific and non-specific liver-localized CD8 T cells did not display attraction towards the infection site, suggesting that the search for malaria liver stages occurs randomly. However, a small fraction (15-20%) displayed weak but detectable attraction towards parasites which already had been surrounded by several T cells. We found that speeds and turning angles correlated with attraction, suggesting that understanding mechanisms that determine the speed of T cell movement in the liver may improve the efficacy of future T cell-based vaccines. Stochastic simulations suggest that a small movement bias towards the parasite dramatically reduces the number of CD8 T cells needed to eliminate all malaria liver stages, but to detect such attraction by individual cells requires data from long imaging experiments which are not currently feasible. Importantly, as far as we know this is the first demonstration of how activated/memory CD8 T cells might search for the pathogen in nonlymphoid tissues a few hours after infection. We have also established a framework for how attraction of individual T cells towards a location in 3D can be rigorously evaluated. T606-phosphorylation deprives the function of Kaiso as a transcription factor and cancer repressor It is well known that Kaiso protein encoded by ZBTB33 gene is a transcription repressor and that Kaiso-P120ctn interaction increases the shift of Kaiso from the nucleus into the cytoplasm. However, the regulatory mechanisms of Kaiso compartmentalization are far from clear. Here, we reported that AKT1 could phosphorylate 606-threonine residue (T606) within the RSSTIP motif of Kaiso in the cytoplasm. The T606-phosphorylated Kaiso (pT606-Kaiso) could directly bind to 14-3-3 family proteins and the depletion of T606 phosphorylation by T606A mutation abolished most of the Kaiso-14-3-3 binding. In addition, the Kaiso-P120ctn interaction was essential for the pT606-Kaiso accumulation in the cytoplasm. Notably, enforced 14-3-3{sigma} (SFN) overexpression could increase the pT606-Kaiso accumulation in the cytoplasm and de-repress the transcription of Kaiso target gene CDH1. Decreased amounts of both pT606-Kaiso and CDH1 proteins were frequently observed in human gastric cancer tissues relative to paired normal controls. The mRNA levels of 14-3-3{sigma} and Kaiso target gene CDH1 were positively and significantly correlated with each other in bioinformatics analyses using publicly available RNA-seq datasets for human normal tissues (n=11688, r=0.60, p<0.001) in the GTEx project and for cancer cell lines (n=1156, r=0.41, p<0.001) in the CCLE project. Furthermore, Kaiso T606A mutant (unable to be phosphorylated) significantly increased the migration and invasion of cancer cells in vitro as well as boosted the growth of these cells in vivo. In conclusion, Kaiso could be phosphorylated by AKT1 at the T606 and the pT606-Kaiso accumulates in the cytoplasm through binding to 14-3-3/P120ctn that de-represses the expression of Kaiso target gene CDH1 in normal tissues. Decreased Kaiso phosphorylation may contribute to the development of gastrointestinal cancer. The status of Kaiso phosphorylation is a determinant factor for the role of Kaiso in the development of cancer. The Inflammatory Role of Milk Fat Globule Epidermal Growth Factor VIII in Angiotensin II Induced Arterial Remodeling BackgroundAngiotensin II (Ang II) and milk fat globule-epidermal growth factor VIII (MFG-E8) are involved in age-associated arterial remodeling; however, the inflammatory role of MFG-E8 in Ang II associated arterial remodeling with aging remains unknown.

Methods and ResultsIn this study, 30-week-old MFG-E8 knock out (KO) and age-matched wild-type (WT) mice were infused with Ang II or saline. After infusion the with Ang II, the aortic molecular, cellular, and structural remodeling were observed in mice and compared to those infused with saline, but these effects were dependent on the expression of MFG-E8: (1) In the WT mice, Ang II infusion substantially increased intimal-medial thickness, elastic lamina degradation, collagen deposition, and the proliferation of VSMCs; in contrast, in the KO mice, these effects were significantly reduced; (2) In the WT mice, Ang II treatment significantly increased the activation and expression of MMP2, TGF-{beta}1, and its downstream signaling molecule p-SMAD2, and collagen type I production, however, in the KO mice, these molecular effects were significantly reduced; (3) In the WT mice, Ang II treatment increased inflammatory p-NF-{kappa}B p65, MCP1, TNF-, ICAM1, and VCAM1 molecular expression, while conversely, in the KO mice, no significant inflammatory changes were found; (4) Importantly, compared to untreated control mice with a wide range of age from 4-96 weeks, Ang II infused "younger" mice produced an "older" arterial inflammatory phenotype, which was alleviated by MFG-E8 deficiency.

ConclusionsMFG-E8 mediates Ang II associated arterial inflammatory remodeling. Targeting MFG-E8 is a novel molecular approach to curb adverse arterial remodeling during aging and hypertension.

CLINICAL PERSPECTIVEO_ST_ABSWhat Is New?C_ST_ABS* Both Ang II and MFG-E8 increases are involved in proinflammatory arterial remodeling mediating the molecular, cellular and tissue events in aging and hypertension.
* MFG-E8 is essential for Ang II induced and age-associated adverse arterial remodeling via the increase of proinflammation, intimal medial thickening, elastin fragmentation, collagen deposition, and VSMC proliferation.


What Are the Clinical Implications?Since MFG-E8 mediates Ang II induced proinflammation in arterial wall remodeling in aging and hypertension, targeting MFG-E8 is a potential molecular approach to curb inflammatory arterial remodeling, maintaining the health of the vascular system during aging and hypertension.