Discriminative stimulus properties of chlordiazepoxide and zolpidem. Agonist and antagonist effects of CGS 9896 and ZK 91296

Astrocyte-Like Cells Transcriptome Changes After Exposure to a Low and Non-cytotoxic MeHg Concentration

The central nervous system is the main target of MeHg toxicity and glial cells are the first line of defense; however, their true role remains unclear. This study aimed to identify the global map of human glial-like (U87) cells transcriptome after exposure to a non-toxic and non-lethal MeHg concentration and to investigate the related molecular changes. U87 cells were exposed upon 0.1, 0.5, and 1 µM MeHg for 4 and 24 h. Although no changes were observed in the percentage of viable cells, the metabolic viability was significantly decreased after exposure to 1 µM MeHg for 24 h; thus, the non-toxic concentration of 0.1 µM MeHg was chosen to perform microarray analysis. Significant changes in U87 cells transcriptome were observed only after 24 h. The expression of 392 genes was down regulated while 431 genes were up-regulated. Gene ontology showed alterations in biological processes (75%), cellular components (21%), and molecular functions (4%). The main pathways showed by KEGG and Reactome were cell cycle regulation and Rho GTPase signaling. The complex mechanism of U87 cells response against MeHg exposure indicates that even a low and non-toxic concentration is able to alter the gene expression profile.

Metabolic labeling of secreted matrix to investigate cell-material interactions in tissue engineering and mechanobiology

Re-creating features of the native extracellular matrix (ECM) with engineered biomaterials has become a valuable tool to probe the influence of ECM properties on cellular functions (e.g., differentiation) and toward the engineering of tissues. However, characterization of newly secreted (nascent) matrix and turnover, which are important in the context of cells interacting with these biomaterials, has been limited by a lack of tools. We developed a protocol to visualize and quantify the spatiotemporal evolution of newly synthesized and deposited matrix by cells that are either cultured atop (2D) or embedded within (3D) biomaterial systems (e.g., hydrogels, fibrous matrices). This technique relies on the incorporation of a noncanonical amino acid (azidohomoalanine) into proteins as they are synthesized. Deposited nascent ECM components are then visualized with fluorescent cyclooctynes via copper-free cycloaddition for spatiotemporal analysis or modified with cleavable biotin probes for identification. Here we describe the preparation of hyaluronic acid hydrogels through ultraviolet or visible light induced cross-linking for 2D and 3D cell culture, as well as the fluorescent labeling of nascent ECM deposited by cells during culture. We also provide protocols for secondary immunofluorescence of specific ECM components and ImageJ-based ECM quantification methods. Hyaluronic acid polymer synthesis takes 2 weeks to complete, and hydrogel formation for 2D or 3D cell culture is performed in 2-3 h. Lastly, we detail the identification of nascent proteins, including enrichment, preparation and analysis with mass spectrometry, which can be completed in 10 d.

Quantitative multiplexed proteomics analysis reveals reshaping of the lysine 2-hydroxyisobutyrylome in Fusarium graminearum by tebuconazole

Backgrounds

Lysine 2-hydroxyisobutyrylation (Khib) is a newly discovered posttranslational modification (PTM) and has been identified in several prokaryotic and eukaryotic organisms. Fusarium graminearum, a major pathogen of Fusarium head blight (FHB) in cereal crops, can cause considerable yield loss and produce various mycotoxins that threaten human health. The application of chemical fungicides such as tebuconazole (TEC) remains the major method to control this pathogen. However, the distribution of Khib in F. graminearum and whether Khib is remodified in response to fungicide stress remain unknown.

Results

Here, we carried out a proteome-wide analysis of Khib in F. graminearum, identifying the reshaping of the lysine 2-hydroxyisobutyrylome by tebuconazole, using the most recently developed high-resolution LC–MS/MS technique in combination with high-specific affinity enrichment. Specifically, 3501 Khib sites on 1049 proteins were identified, and 1083 Khib sites on 556 modified proteins normalized to the total protein content were changed significantly after TEC treatment. Bioinformatics analysis showed that Khib proteins are involved in a wide range of biological processes and may be involved in virulence and deoxynivalenol (DON) production, as well as sterol biosynthesis, in F. graminearum.

Conclusions

Here, we provided a wealth of resources for further study of the roles of Khib in the fungicide resistance of F. graminearum. The results enhanced our understanding of this PTM in filamentous ascomycete fungi and provided insight into the remodification of Khib sites during azole fungicide challenge in F. graminearum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08372-4.

Endothelial Cell RNA-Seq Data: Differential Expression and Functional Enrichment Analyses to Study Phenotypic Switching

RNA-seq is a common approach used to explore gene expression data between experimental conditions or cell types and ultimately leads to information that can shed light on the biological processes involved and inform further hypotheses. While the protocols required to generate samples for sequencing can be performed in most research facilities, the resulting computational analysis is often an area in which researchers have little experience. Here we present a user-friendly bioinformatics workflow which describes the methods required to take raw data produced by RNA sequencing to interpretable results. Widely used and well documented tools are applied. Data quality assessment and read trimming were performed by FastQC and Cutadapt, respectively. Following this, STAR was utilized to map the trimmed reads to a reference genome and the alignment was analyzed by Qualimap. The subsequent mapped reads were quantified by featureCounts. DESeq2 was used to normalize and perform differential expression analysis on the quantified reads, identifying differentially expressed genes and preparing the data for functional enrichment analysis. Gene set enrichment analysis identified enriched gene sets from the normalized count data and clusterProfiler was used to perform functional enrichment against the GO, KEGG, and Reactome databases. Example figures of the functional enrichment analysis results were also generated. The example data used in the workflow are derived from HUVECs, an in vitro model used in the study of endothelial cells, published and publicly available for download from the European Nucleotide Archive.

Pathogen resistance in Sphagneticola trilobata (Singapore daisy): molecular associations and differentially expressed genes in response to disease from a widespread fungus

Understanding the molecular associations underlying pathogen resistance in invasive plant species is likely to provide useful insights into the effective control of alien plants, thereby facilitating the conservation of native biodiversity. In the current study, we investigated pathogen resistance in an invasive clonal plant, Sphagneticola trilobata, at the molecular level. Sphagneticola trilobata (i.e., Singapore daisy) is a noxious weed that affects both terrestrial and aquatic ecosystems, and is less affected by pathogens in the wild than co-occurring native species. We used Illumina sequencing to investigate the transcriptome of S. trilobata following infection by a globally distributed generalist pathogen (Rhizoctonia solani). RNA was extracted from leaves of inoculated and un-inoculated control plants, and a draft transcriptome of S. trilobata was generated to examine the molecular response of this species following infection. We obtained a total of 49,961,014 (94.3%) clean reads for control (un-inoculated plants) and 54,182,844 (94.5%) for the infected treatment (inoculated with R. solani). Our analyses facilitated the discovery of 117,768 de novo assembled contigs and 78,916 unigenes. Of these, we identified 3506 differentially expressed genes and 60 hormones associated with pathogen resistance. Numerous genes, including candidate genes, were associated with plant-pathogen interactions and stress response in S. trilobata. Many recognitions, signaling, and defense genes were differentially regulated between treatments, which were confirmed by qRT-PCR. Overall, our findings improve our understanding of the genes and molecular associations involved in plant defense of a rapidly spreading invasive clonal weed, and serve as a valuable resource for further work on mechanism of disease resistance and managing invasive plants.

Integration of spatial and single-cell transcriptomic data elucidates mouse organogenesis

Molecular profiling of single cells has advanced our knowledge of the molecular basis of development. However, current approaches mostly rely on dissociating cells from tissues, thereby losing the crucial spatial context of regulatory processes. Here, we apply an image-based single-cell transcriptomics method, sequential fluorescence in situ hybridization (seqFISH), to detect mRNAs for 387 target genes in tissue sections of mouse embryos at the 8-12 somite stage. By integrating spatial context and multiplexed transcriptional measurements with two single-cell transcriptome atlases, we characterize cell types across the embryo and demonstrate that spatially resolved expression of genes not profiled by seqFISH can be imputed. We use this high-resolution spatial map to characterize fundamental steps in the patterning of the midbrain-hindbrain boundary (MHB) and the developing gut tube. We uncover axes of cell differentiation that are not apparent from single-cell RNA-sequencing (scRNA-seq) data, such as early dorsal-ventral separation of esophageal and tracheal progenitor populations in the gut tube. Our method provides an approach for studying cell fate decisions in complex tissues and development.

Complete genome sequence of Bacillus velezensis YYC, a bacterium isolated from the tomato rhizosphere

The Bacillus velezensis YYC strain was isolated from the tomato rhizosphere. In a previous experiment, it increased tomato growth and induced systemic resistance against Ralstonia solanacearum. However, information on its genomic content is lacking. The complete genome sequence of the bacterium was described in this study. The genome size was 3,973,236 bp and consisted of 4034 genes in total, with a mean G + C content of 46.52%. In addition, 86 tRNAs and 27 ribosomal RNAs were identified. Fourteen clusters of secondary metabolites were identified. The KEGG database analysis showed that 69 genes were related to quorum sensing, which were important for microbe-plant interaction. In addition, genes involved in promoting plant growth and triggering plant immunity were identified from the genome. Based on digital DNA-DNA hybridizations (dDDH), B. velezensis YYC was most closely related with B. velezensis FZB42. The complete genome data of B. velezensis YYC will provide a basis for explanation of its growth-promoting mechanism and biocontrol mechanism.

The Rat Genome Database (RGD) facilitates genomic and phenotypic data integration across multiple species for biomedical research

Model organism research is essential for discovering the mechanisms of human diseases by defining biologically meaningful gene to disease relationships. The Rat Genome Database (RGD, ( https://rgd.mcw.edu )) is a cross-species knowledgebase and the premier online resource for rat genetic and physiologic data. This rich resource is enhanced by the inclusion and integration of comparative data for human and mouse, as well as other human disease models including chinchilla, dog, bonobo, pig, 13-lined ground squirrel, green monkey, and naked mole-rat. Functional information has been added to records via the assignment of annotations based on sequence similarity to human, rat, and mouse genes. RGD has also imported well-supported cross-species data from external resources. To enable use of these data, RGD has developed a robust infrastructure of standardized ontologies, data formats, and disease- and species-centric portals, complemented with a suite of innovative tools for discovery and analysis. Using examples of single-gene and polygenic human diseases, we illustrate how data from multiple species can help to identify or confirm a gene as involved in a disease and to identify model organisms that can be studied to understand the pathophysiology of a gene or pathway. The ultimate aim of this report is to demonstrate the utility of RGD not only as the core resource for the rat research community but also as a source of bioinformatic tools to support a wider audience, empowering the search for appropriate models for human afflictions.

Description and genome analysis of Microvirga antarctica sp. nov., a novel pink-pigmented psychrotolerant bacterium isolated from Antarctic soil

A novel pink-pigmented bacterium, designated strain 3D7^T, was isolated during an investigation of potential psychrotolerant species from Antarctic soil. Cells of the isolate were observed to be rod-shaped (0.7-0.9 × 1.0-2.2 µm), Gram-stain negative and non-motile. It was able to grow at 4-32  °C, pH 7.0-10.0 and in the presence of 0-3% (w/v) NaCl. Phylogenetic analysis based on 16S rRNA gene sequences showed that strain 3D7^T belongs to the genus Microvirga and was most closely related to 'Microvirga brassicacearum' CDVBN77^T (98.3%), Microvirga subterranea DSM 14364 ^T (96.8%), Microvirga guangxiensis 25B^T (96.5%) and Microvirga aerophila DSM 21344 ^T (96.5%). The predominant quinone was ubiquinone 10 (Q-10), and the major fatty acids were summed feature 8 (C_18:1ω7c and/or C_18:1ω6c) and C_19:0 cyclo ω8c. The predominant polar lipids were phosphatidylcholine and phosphatidylethanolamine. The genomic DNA G + C content of strain 3D7^T was 63.5 mol%. Its genome sequence showed genes encoding phosphatases and lipases. Genetic machinery related to carbohydrate-active enzymes and secondary metabolites were also observed. The average nucleotide identity and digital DNA-DNA hybridization values based on whole genome sequences of strain 3D7^T and its closely related species were below the threshold range for species determination. Phenotypic, chemotaxonomic, phylogenetic and genomic analyses suggested that strain 3D7^T represents a novel species of the genus Microvirga, for which the name Microvirga antarctica sp. nov. is proposed. The type strain is 3D7^T (= CGMCC 1.13821^T = KCTC 72465^T).

Transcriptome repository of North-Western Himalayan endangered medicinal herbs: a paramount approach illuminating molecular perspective of phytoactive molecules and secondary metabolism

Medicinal plants of the North-Western Himalayan region are known for their unprecedented biodiversity and valuable secondary metabolites that are unique to this dynamic geo-climatic region. From ancient times these medicinal herbs have been used traditionally for their therapeutic potentials. But from the last 2 decades increasing pharmaceutical demand, illegal and unorganized trade of these medicinal plants have accelerated the rate of over-exploitation in a non-scientific manner. In addition, climate change and anthropogenic activities also affected their natural habitat and driving most of these endemic plant species to critically endangered that foresee peril of mass extinction from this eco-region. Hence there is an urgent need for developing alternative sustainable approaches and policies to utilize this natural bioresource ensuring simultaneous conservation. Hither, arise the advent of sequencing-based transcriptomic studies significantly contributes to better understand the background of important metabolic pathways and related genes/enzymes of high-value medicinal herbs, in the absence of genomic information. The use of comparative transcriptomics in conjunction with biochemical techniques in North-Western Himalayan medicinal plants has resulted in significant advances in the identification of the molecular players involved in the production of secondary metabolic pathways over the last decade. This information could be used to further engineer metabolic pathways and breeding programs, ultimately leading to the development of in vitro systems dedicated to the production of pharmaceutically important secondary metabolites at the industrial level. Collectively, successful adoption of these approaches can certainly ensure the sustainable utilization of Himalayan bioresource by reducing the pressure on the wild population of these critically endangered medicinal herbs. This review provides novel insight as a transcriptome-based bioresource repository for the understanding of important secondary metabolic pathways genes/enzymes and metabolism of endangered high-value North-Western Himalayan medicinal herbs, so that researchers across the globe can effectively utilize this information for devising effective strategies for the production of pharmaceutically important compounds and their scale-up for sustainable usage and take a step forward in omics-based conservation genetics.

Nesterenkonia ebinurensis sp. nov., a Novel Actinobacterium Isolated From Populus euphratica

During our studies focused on the microorganism diversity and community structure of Populus euphratica at Ebinur lake wetland nature reserve in Xinjiang Uyghur Autonomous Region, PR China, a Gram-positive, aerobic, short rod-shaped bacterium without flagellum, designated as MD2^T, was isolated from a piece of resina on Populus euphratica. The isolate grew at temperature of 10-45 °C (optimum 37 °C), pH of 7.0-12.0 (optimum pH 8.0) and NaCl concentration of 1-18% (optimum 3%, w/v). Based on the 16S rRNA gene sequences and the phylogenetic analysis, the strain shared the highest sequence similarities to Nesterenkonia alkaliphila JCM 19766^T (96.3%), Nesterenkonia populi KCTC 29119^T (95.9%), Nesterenkonia alba CCTCC AB 207011^T (95.5%), and was placed within the radiation of Nesterenkonia species in the phylogenetic trees. The draft genome of the isolate was sequenced, which comprised 3,739,891 bp with G + C content of 63 mol%, and was annotated to contain 3614 protein-coding genes, 44 tRNA genes and 5 rRNA genes. Chemotaxonomic analysis indicated that the main respiratory quinones were MK-8 and MK-9, the predominant cellular fatty acids were anteiso-C_17:0, anteiso-C_15:0 and iso-C_16:0, the major polar lipids consisted of diphosphatidylglycerol, phosphatidylglycerol andphosphatidylinositol. According to the phenotypic, chemotaxonomic and phylogenetic features, strain MD2^T is considered to represent a novel species, for which the name Nesterenkonia ebinurensis sp. nov. is proposed. The type strain is MD2^T (= KCTC 52999^T = MCCC 1K03343^T).

Early-life trauma endophenotypes and brain circuit-gene expression relationships in functional neurological (conversion) disorder

Functional neurological (conversion) disorder (FND) is a neuropsychiatric condition whereby individuals present with sensorimotor symptoms incompatible with other neurological disorders. Early-life maltreatment (ELM) is a risk factor for developing FND, yet few studies have investigated brain network-trauma relationships in this population. In this neuroimaging-gene expression study, we used two graph theory approaches to elucidate ELM subtype effects on resting-state functional connectivity architecture in 30 patients with motor FND. Twenty-one individuals with comparable depression, anxiety, and ELM scores were used as psychiatric controls. Thereafter, we compared trauma endophenotypes in FND with regional differences in transcriptional gene expression as measured by the Allen Human Brain Atlas (AHBA). In FND patients only, we found that early-life physical abuse severity, and to a lesser extent physical neglect, correlated with corticolimbic weighted-degree functional connectivity. Connectivity profiles influenced by physical abuse occurred in limbic (amygdalar-hippocampal), paralimbic (cingulo-insular and ventromedial prefrontal), and cognitive control (ventrolateral prefrontal) areas, as well as in sensorimotor and visual cortices. These findings held adjusting for individual differences in depression/anxiety, PTSD, and motor phenotypes. In FND, physical abuse also correlated with amygdala and insula coupling to motor cortices. In exploratory analyses, physical abuse correlated connectivity maps overlapped with the AHBA spatial expression of three gene clusters: (i) neuronal morphogenesis and synaptic transmission genes in limbic/paralimbic areas; (ii) locomotory behavior and neuronal generation genes in left-lateralized structures; and (iii) nervous system development and cell motility genes in right-lateralized structures. These circuit-specific architectural profiles related to individual differences in childhood physical abuse burden advance our understanding of the pathophysiology of FND.

Exploring tissue architecture using spatial transcriptomics

Deciphering the principles and mechanisms by which gene activity orchestrates complex cellular arrangements in multicellular organisms has far-reaching implications for research in the life sciences. Recent technological advances in next-generation sequencing- and imaging-based approaches have established the power of spatial transcriptomics to measure expression levels of all or most genes systematically throughout tissue space, and have been adopted to generate biological insights in neuroscience, development and plant biology as well as to investigate a range of disease contexts, including cancer. Similar to datasets made possible by genomic sequencing and population health surveys, the large-scale atlases generated by this technology lend themselves to exploratory data analysis for hypothesis generation. Here we review spatial transcriptomic technologies and describe the repertoire of operations available for paths of analysis of the resulting data. Spatial transcriptomics can also be deployed for hypothesis testing using experimental designs that compare time points or conditions-including genetic or environmental perturbations. Finally, spatial transcriptomic data are naturally amenable to integration with other data modalities, providing an expandable framework for insight into tissue organization.

A proximity-dependent biotinylation map of a human cell

Compartmentalization is a defining characteristic of eukaryotic cells, and partitions distinct biochemical processes into discrete subcellular locations. Microscopy¹ and biochemical fractionation coupled with mass spectrometry^2-4 have defined the proteomes of a variety of different organelles, but many intracellular compartments have remained refractory to such approaches. Proximity-dependent biotinylation techniques such as BioID provide an alternative approach to define the composition of cellular compartments in living cells^5-7. Here we present a BioID-based map of a human cell on the basis of 192 subcellular markers, and define the intracellular locations of 4,145 unique proteins in HEK293 cells. Our localization predictions exceed the specificity of previous approaches, and enabled the discovery of proteins at the interface between the mitochondrial outer membrane and the endoplasmic reticulum that are crucial for mitochondrial homeostasis. On the basis of this dataset, we created humancellmap.org as a community resource that provides online tools for localization analysis of user BioID data, and demonstrate how this resource can be used to understand BioID results better.

A genome-wide atlas of co-essential modules assigns function to uncharacterized genes

A central question in the post-genomic era is how genes interact to form biological pathways. Measurements of gene dependency across hundreds of cell lines have been used to cluster genes into 'co-essential' pathways, but this approach has been limited by ubiquitous false positives. In the present study, we develop a statistical method that enables robust identification of gene co-essentiality and yields a genome-wide set of functional modules. This atlas recapitulates diverse pathways and protein complexes, and predicts the functions of 108 uncharacterized genes. Validating top predictions, we show that TMEM189 encodes plasmanylethanolamine desaturase, a key enzyme for plasmalogen synthesis. We also show that C15orf57 encodes a protein that binds the AP2 complex, localizes to clathrin-coated pits and enables efficient transferrin uptake. Finally, we provide an interactive webtool for the community to explore our results, which establish co-essentiality profiling as a powerful resource for biological pathway identification and discovery of new gene functions.

PCprophet: a framework for protein complex prediction and differential analysis using proteomic data

Despite the availability of methods for analyzing protein complexes, systematic analysis of complexes under multiple conditions remains challenging. Approaches based on biochemical fractionation of intact, native complexes and correlation of protein profiles have shown promise. However, most approaches for interpreting cofractionation datasets to yield complex composition and rearrangements between samples depend considerably on protein-protein interaction inference. We introduce PCprophet, a toolkit built on size exclusion chromatography-sequential window acquisition of all theoretical mass spectrometry (SEC-SWATH-MS) data to predict protein complexes and characterize their changes across experimental conditions. We demonstrate improved performance of PCprophet over state-of-the-art approaches and introduce a Bayesian approach to analyze altered protein-protein interactions across conditions. We provide both command-line and graphical interfaces to support the application of PCprophet to any cofractionation MS dataset, independent of separation or quantitative liquid chromatography-MS workflow, for the detection and quantitative tracking of protein complexes and their physiological dynamics.

27-Hydroxycholesterol is a specific factor in the neoplastic microenvironment of HCC that causes MDR via GRP75 regulation of the redox balance and metabolic reprogramming

OBJECTIVE

Due to the tissue specificity of the liver, long-term exposure to a high concentration of 27-hydroxycholesterol (27HC) is a special characteristic of the tumour microenvironment in hepatocellular carcinoma (HCC). However, what occurs after HCC cells are long-term exposure to 27HC and the molecular mechanisms involved remain largely unexamined.

METHODS

A long-term 27HC-treated HepG2 cell line and the xenografts in nude mice were used as experimental models. Molecular mechanisms were investigated using bioinformatics analysis and molecular biological experiments.

RESULTS

Here, we found that by inducing an increase in oxidative stress signalling, 27HC activated glucose-regulated protein 75 (GRP75). On the one hand, GRP75 resulted in a change in the redox balance by regulating ROS generation and antioxidant system activity via affecting MMP, NRF2, HO-1, and NQO1 levels. On the other hand, GRP75 modified the metabolic reprogramming process by regulating key factors (HIF-1α, p-Akt, and c-myc) and glucose uptake, facilitating HCC cell growth in the inhospitable microenvironment. These two factors caused HCC cells to resist 27HC-induced cytotoxicity and attain multidrug resistance (MDR).

CONCLUSIONS

Our present study not only identified 27HC, a characteristic component of the neoplastic microenvironment of HCC that causes MDR via GRP75 to regulate the redox balance and metabolic reprogramming, but also revealed that targeted intervention by the "switch"-like molecule GRP75 could reverse the effect of 27HC from cancer promotion to cytotoxicity in HCC, suggesting a new strategy for specific intervention of HCC.

Genomic insights into population history and biological adaptation in Oceania

The Pacific region is of major importance for addressing questions regarding human dispersals, interactions with archaic hominins and natural selection processes¹. However, the demographic and adaptive history of Oceanian populations remains largely uncharacterized. Here we report high-coverage genomes of 317 individuals from 20 populations from the Pacific region. We find that the ancestors of Papuan-related ('Near Oceanian') groups underwent a strong bottleneck before the settlement of the region, and separated around 20,000-40,000 years ago. We infer that the East Asian ancestors of Pacific populations may have diverged from Taiwanese Indigenous peoples before the Neolithic expansion, which is thought to have started from Taiwan around 5,000 years ago^2-4. Additionally, this dispersal was not followed by an immediate, single admixture event with Near Oceanian populations, but involved recurrent episodes of genetic interactions. Our analyses reveal marked differences in the proportion and nature of Denisovan heritage among Pacific groups, suggesting that independent interbreeding with highly structured archaic populations occurred. Furthermore, whereas introgression of Neanderthal genetic information facilitated the adaptation of modern humans related to multiple phenotypes (for example, metabolism, pigmentation and neuronal development), Denisovan introgression was primarily beneficial for immune-related functions. Finally, we report evidence of selective sweeps and polygenic adaptation associated with pathogen exposure and lipid metabolism in the Pacific region, increasing our understanding of the mechanisms of biological adaptation to island environments.

Identification of four key prognostic genes and three potential drugs in human papillomavirus negative head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is a common tumor worldwide with poor prognosis. The pathogenesis of human papillomavirus (HPV)-positive and HPV-negative HNSCCs differs. However, few studies have considered the HPV status when identifying biomarkers for HNSCC. Thus, the identification of biomarkers for HPV-positive and HPV-negative HNSCCs is urgently needed.

METHODS

Three microarray datasets from Gene Expression Omnibus (GEO) were analyzed, and the differentially expressed genes (DEGs) were obtained. Then, functional enrichment pathway analysis was performed and protein-protein interaction (PPI) networks were constructed. The expression of hub genes at both the mRNA and protein level was determined in Oncomine, The Cancer Genome Atlas (TCGA) and the Human Protein Atlas (HPA). In addition, survival analysis of the patient stratified by HPV status and the expression levels of key genes were performed based on TCGA data. The role of AREG, STAG3, CAV1 and C19orf57 in cancer were analyzed through Gene set enrichment analysis (GSEA). The top ten small molecule drugs were identified and the therapeutic value of zonisamide, NVP-AUY922, PP-2 and fostamatinib was further evaluated in six HPV-negative HNSCC cell lines. Finally, the therapeutic value of NVP-AUY922 was tested in vivo based on three HPV-negative HNSCC models, and statistical analysis was performed.

RESULTS

In total, 47 DEGs were obtained, 11 of which were identified as hub genes. Biological process analysis indicated that the hub genes were associated with the G1/S transition of the mitotic cell cycle. Survival analysis uncovered that the prognostic value of AREG, STAG3, C19orf57 and CAV1 differed between HPV-positive and HPV-negative patients. Gene set enrichment analysis (GSEA) showed the role of AREG, STAG3 and CAV1 in dysregulated pathways of tumor. Ten small molecules were identified as potential drugs specifically for HPV-positive or HPV-negative patients; three-NVP-AUY922, fostamatinib and PP-2-greatly inhibited the proliferation of six HPV-negative HNSCC cell lines in vitro, and NVP-AUY922 inhibited three HPV-negative HNSCC xenografts in vivo.

CONCLUSIONS

In conclusion, AREG, STAG3, C19orf57 and CAV1 are key prognostic factors and potential therapeutic targets in HPV-negative HNSCC. NVP-AUY922, fostamatinib and PP-2 may be effective drugs for HPV-negative HNSCC.

Human Chorionic Gonadotropin Modulates the Transcriptome of the Myometrium and Cervix in Late Gestation

Human chorionic gonadotropin (hCG) is a critical hormone for the establishment and maintenance of pregnancy. hCG administration prevents the onset of preterm labor in mice; yet, the transcriptomic changes associated with this tocolytic effect that take place in the myometrium and cervix have not been elucidated. Herein, we implemented both discovery and targeted approaches to investigate the transcriptome of the myometrium and cervix after hCG administration. Pregnant mice were intraperitoneally injected with 10 IU of hCG on 13.0, 15.0, and 17.0 days post coitum, and the myometrium and cervix were collected. RNA sequencing was performed to determine differentially expressed genes, enriched biological processes, and impacted KEGG pathways. Multiplex qRT-PCR was performed to investigate the expression of targeted contractility- and inflammation-associated transcripts. hCG administration caused the differential expression of 720 genes in the myometrium. Among the downregulated genes, enriched biological processes were primarily associated with regulation of transcription. hCG administration downregulated key contractility genes, Gja1 and Oxtr, but upregulated the prostaglandin-related genes Ptgfr and Ptgs2 and altered the expression of inflammation-related genes in the myometrium. In the cervix, hCG administration caused differential expression of 3348 genes that were related to inflammation and host defense, among others. The downregulation of key contractility genes and upregulation of prostaglandin-related genes were also observed in the cervix. Thus, hCG exerts tocolytic and immunomodulatory effects in late gestation by altering biological processes in the myometrium and cervix, which should be taken into account when considering hCG as a potential treatment to prevent the premature onset of labor.

Regional Differences in S-Nitrosylation in the Cortex, Striatum, and Hippocampus of Juvenile Male Mice

Nitric oxide (NO) is a multifunctional neurotransmitter that plays a major role in neuronal and synaptic functions. S-nitrosylation (SNO), the NO-mediated protein posttransitional modification (PTM), is known to regulate physiological and pathological processes in the brain. However, the physiological role in different neuroanatomical brain regions has not been well investigated. To understand the role of SNO in the brain of juvenile WT mice, we used SNOTRAP technology. We mapped the SNO-proteome in three different neuroanatomical regions: cortex, striatum, and hippocampus. By conducting systems biology analysis, we found that the three brain regions share similar biological processes (BP) including biogenesis and developmental processes. Exclusive and different BP and molecular functions were found for each of the regions. Unraveling the BP and signaling mechanisms of SNO in the cortex, striatum, and hippocampus may help to understand the functional differences between the three regions under physiological conditions.

TDP-43 aggregation induced by oxidative stress causes global mitochondrial imbalance in ALS

Amyotrophic lateral sclerosis (ALS) was initially thought to be associated with oxidative stress when it was first linked to mutant superoxide dismutase 1 (SOD1). The subsequent discovery of ALS-linked genes functioning in RNA processing and proteostasis raised the question of how different biological pathways converge to cause the disease. Both familial and sporadic ALS are characterized by the aggregation of the essential DNA- and RNA-binding protein TDP-43, suggesting a central role in ALS etiology. Here we report that TDP-43 aggregation in neuronal cells of mouse and human origin causes sensitivity to oxidative stress. Aggregated TDP-43 sequesters specific microRNAs (miRNAs) and proteins, leading to increased levels of some proteins while functionally depleting others. Many of those functionally perturbed gene products are nuclear-genome-encoded mitochondrial proteins, and their dysregulation causes a global mitochondrial imbalance that augments oxidative stress. We propose that this stress-aggregation cycle may underlie ALS onset and progression.

Transcriptome of nasopharyngeal samples from COVID-19 patients and a comparative analysis with other SARS-CoV-2 infection models reveal disparate host responses against SARS-CoV-2

BACKGROUND

Although it is becoming evident that individual's immune system has a decisive influence on SARS-CoV-2 disease progression, pathogenesis is largely unknown. In this study, we aimed to profile the host transcriptome of COVID-19 patients from nasopharyngeal samples along with virus genomic features isolated from respective host, and a comparative analyses of differential host responses in various SARS-CoV-2 infection systems.

RESULTS

Unique and rare missense mutations in 3C-like protease observed in all of our reported isolates. Functional enrichment analyses exhibited that the host induced responses are mediated by innate immunity, interferon, and cytokine stimulation. Surprisingly, induction of apoptosis, phagosome, antigen presentation, hypoxia response was lacking within these patients. Upregulation of immune and cytokine signaling genes such as CCL4, TNFA, IL6, IL1A, CCL2, CXCL2, IFN, and CCR1 were observed in lungs. Lungs lacked the overexpression of ACE2 as suspected, however, high ACE2 but low DPP4 expression was observed in nasopharyngeal cells. Interestingly, directly or indirectly, viral proteins specially non-structural protein mediated overexpression of integrins such as ITGAV, ITGA6, ITGB7, ITGB3, ITGA2B, ITGA5, ITGA6, ITGA9, ITGA4, ITGAE, and ITGA8 in lungs compared to nasopharyngeal samples suggesting the possible way of enhanced invasion. Furthermore, we found comparatively highly expressed transcription factors such as CBP, CEBP, NFAT, ATF3, GATA6, HDAC2, TCF12 which have pivotal roles in lung injury.

CONCLUSIONS

Even though this study incorporates a limited number of cases, our data will provide valuable insights in developing potential studies to elucidate the differential host responses on the viral pathogenesis in COVID-19, and incorporation of further data will enrich the search of an effective therapeutics.

The Rcs stress response inversely controls surface and CRISPR-Cas adaptive immunity to discriminate plasmids and phages

Bacteria harbour multiple innate defences and adaptive CRISPR-Cas systems that provide immunity against bacteriophages and mobile genetic elements. Although some bacteria modulate defences in response to population density, stress and metabolic state, a lack of high-throughput methods to systematically reveal regulators has hampered efforts to understand when and how immune strategies are deployed. We developed a robust approach called SorTn-seq, which combines saturation transposon mutagenesis, fluorescence-activated cell sorting and deep sequencing to characterize regulatory networks controlling CRISPR-Cas immunity in Serratia sp. ATCC 39006. We applied our technology to assess csm gene expression for ~300,000 mutants and uncovered multiple pathways regulating type III-A CRISPR-Cas expression. Mutation of igaA or mdoG activated the Rcs outer-membrane stress response, eliciting cell-surface-based innate immunity against diverse phages via the transcriptional regulators RcsB and RcsA. Activation of this Rcs phosphorelay concomitantly attenuated adaptive immunity by three distinct type I and III CRISPR-Cas systems. Rcs-mediated repression of CRISPR-Cas defence enabled increased acquisition and retention of plasmids. Dual downregulation of cell-surface receptors and adaptive immunity in response to stress by the Rcs pathway enables protection from phage infection without preventing the uptake of plasmids that may harbour beneficial traits.

Gene signature of children with severe respiratory syncytial virus infection

BACKGROUND

The limited treatment options for children with severe respiratory syncytial virus (RSV) infection highlights the need for a comprehensive understanding of the host cellular response during infection. We aimed to identify host genes that are associated with severe RSV disease and to identify drugs that can be repurposed for the treatment of severe RSV infection.

METHODS

We examined clinical data and blood samples from 37 hospitalized children (29 mild and 8 severe) with RSV infection. We tested RNA from blood samples using next-generation sequencing to profile global mRNA expression and identify cellular processes.

RESULTS

Retractions, decreased breath sounds, and tachypnea were associated with disease severity. We observed upregulation of genes related to neutrophil, inflammatory response, blood coagulation, and downregulation of genes related to T cell response in children with severe RSV. Using network-based approach, 43 drugs were identified that are predicted to interact with the gene products of these differentially expressed genes.

CONCLUSIONS

These results suggest that the changes in the expression pattern in the innate and adaptive immune responses may be associated with RSV clinical severity. Compounds that target these cellular processes can be repositioned as candidate drugs in the treatment of severe RSV.

IMPACT

Neutrophil, inflammation, and blood coagulation genes are upregulated in children with severe RSV infection. Expression of T cell response genes are suppressed in cases of severe RSV. Genes identified in this study can contribute in understanding the pathogenesis of RSV disease severity. Drugs that target cellular processes associated with severe RSV can be repositioned as potential therapeutic options.

Donkey genomes provide new insights into domestication and selection for coat color

Current knowledge about the evolutionary history of donkeys is still incomplete due to the lack of archeological and whole-genome diversity data. To fill this gap, we have de novo assembled a chromosome-level reference genome of one male Dezhou donkey and analyzed the genomes of 126 domestic donkeys and seven wild asses. Population genomics analyses indicate that donkeys were domesticated in Africa and conclusively show reduced levels of Y chromosome variability and discordant paternal and maternal histories, possibly reflecting the consequences of reproductive management. We also investigate the genetic basis of coat color. While wild asses show diluted gray pigmentation (Dun phenotype), domestic donkeys display non-diluted black or chestnut coat colors (non-Dun) that were probably established during domestication. Here, we show that the non-Dun phenotype is caused by a 1 bp deletion downstream of the TBX3 gene, which decreases the expression of this gene and its inhibitory effect on pigment deposition.

A new donkey reference genome and comparisons with wild asses yields insights into the evolutionary history of donkey domestication and identifies a genetic variant that results in the non-Dun coat colours of domestic donkeys.

Genome-wide identification of low phosphorus responsive microRNAs in two soybean genotypes by high-throughput sequencing

MicroRNAs (miRNAs) have been reported to be correlated with various stress responses in soybean, but only a few miRNAs have been demonstrated to respond to low phosphorus (LP) stress. To unravel the response mechanisms of miRNAs to low-P stress, the roots of two representative soybean genotypes with different P efficiency, Nannong94-156 (a LP-tolerant genotype) and Bogao (a LP-sensitive genotype), were used for the construction of RNA sequencing (RNA-seq) libraries under low/normal-P treatment by high-throughput sequencing. In total, 603 existing miRNAs and 1699 novel miRNAs belonging to 248 and 1582 families in all samples were identified, respectively. Among these miRNAs, 777 miRNAs were differentially expressed (DE) across different P levels and genotypes. Furthermore, putative targets of DE miRNAs were predicted, and these miRNAs mainly targeted ERF (ethylene responsive factor), auxin response factors (ARF), zinc finger protein, MYB, and NAC domain transcription factors. Gene ontology (GO) analysis showed that targets of DE miRNAs were significantly enriched in binding, metabolic processes, biological regulation, response to stress, and phosphorus metabolic processes. In addition, the expression profiles of chosen P-responsive miRNAs and target genes were validated by quantitative real-time PCR (qRT-PCR). Our study focused on genome-wide miRNA identification in two representative soybean genotypes under low-P stress. Overall, the DE miRNAs across different P levels and genotypes and their putative target genes will provide useful information for further study of miRNAs mediating low-P response and facilitate improvements in soybean breeding.

RUNX3 methylation drives hypoxia-induced cell proliferation and antiapoptosis in early tumorigenesis

Inactivation of tumor suppressor Runt-related transcription factor 3 (RUNX3) plays an important role during early tumorigenesis. However, posttranslational modifications (PTM)-based mechanism for the inactivation of RUNX3 under hypoxia is still not fully understood. Here, we demonstrate a mechanism that G9a, lysine-specific methyltransferase (KMT), modulates RUNX3 through PTM under hypoxia. Hypoxia significantly increased G9a protein level and G9a interacted with RUNX3 Runt domain, which led to increased methylation of RUNX3 at K129 and K171. This methylation inactivated transactivation activity of RUNX3 by reducing interactions with CBFβ and p300 cofactors, as well as reducing acetylation of RUNX3 by p300, which is involved in nucleocytoplasmic transport by importin-α1. G9a-mediated methylation of RUNX3 under hypoxia promotes cancer cell proliferation by increasing cell cycle or cell division, while suppresses immune response and apoptosis, thereby promoting tumor growth during early tumorigenesis. Our results demonstrate the molecular mechanism of RUNX3 inactivation by G9a-mediated methylation for cell proliferation and antiapoptosis under hypoxia, which can be a therapeutic or preventive target to control tumor growth during early tumorigenesis.

Maternal choline supplementation ameliorates Alzheimer's disease pathology by reducing brain homocysteine levels across multiple generations

The lack of effective treatments for Alzheimer's disease (AD) is alarming, considering the number of people currently affected by this disorder and the projected increase over the next few decades. Elevated homocysteine (Hcy) levels double the risk of developing AD. Choline, a primary dietary source of methyl groups, converts Hcy to methionine and reduces age-dependent cognitive decline. Here, we tested the transgenerational benefits of maternal choline supplementation (ChS; 5.0 g/kg choline chloride) in two generations (Gen) of APP/PS1 mice. We first exposed 2.5-month-old mice to the ChS diet and allowed them to breed with each other to generate Gen-1 mice. Gen-1 mice were exposed to the ChS diet only during gestation and lactation; once weaned at postnatal day 21, Gen-1 mice were then kept on the control diet for the remainder of their life. We also bred a subset of Gen-1 mice to each other and obtained Gen-2 mice; these mice were never exposed to ChS. We found that ChS reduced Aβ load and microglia activation, and improved cognitive deficits in old Gen-1 and Gen-2 APP/PS1 mice. Mechanistically, these changes were linked to a reduction in brain Hcy levels in both generations. Further, RNA-Seq data from APP/PS1 hippocampal tissue revealed that ChS significantly changed the expression of 27 genes. These genes were enriched for inflammation, histone modifications, and neuronal death functional classes. Our results are the first to demonstrate a transgenerational benefit of ChS and suggest that modifying the maternal diet with additional choline reduces AD pathology across multiple generations.

Systems biology reveals reprogramming of the S-nitroso-proteome in the cortical and striatal regions of mice during aging process

Cell aging depends on the rate of cumulative oxidative and nitrosative damage to DNA and proteins. Accumulated data indicate the involvement of protein S-nitrosylation (SNO), the nitric oxide (NO)-mediated posttranslational modification (PTM) of cysteine thiols, in different brain disorders. However, the changes and involvement of SNO in aging including the development of the organism from juvenile to adult state is still unknown. In this study, using the state-of-the-art mass spectrometry technology to identify S-nitrosylated proteins combined with large-scale computational biology, we tested the S-nitroso-proteome in juvenile and adult mice in both cortical and striatal regions. We found reprogramming of the S-nitroso-proteome in adult mice of both cortex and striatum regions. Significant biological processes and protein–protein clusters associated with synaptic and neuronal terms were enriched in adult mice. Extensive quantitative analysis revealed a large set of potentially pathological proteins that were significantly upregulated in adult mice. Our approach, combined with large scale computational biology allowed us to perform a system-level characterization and identification of the key proteins and biological processes that can serve as drug targets for aging and brain disorders in future studies.

Differential proteomics study of postharvest Volvariella volvacea during storage at 4 °C

The postharvest storage of Volvariella volvacea is an important factor limiting the industry development. Low-temperature storage is the traditional storage method used for most edible fungi, but V. volvacea undergoes autolysis at low temperature. To understand the molecular mechanism underlying the low-temperature autolysis of V. volvacea after harvesting, fruiting bodies of V. volvacea strain V23 were stored at 4 °C. Based on our previous study, in which the changes of morphological and physiological indexes during storage for 0, 6, 12, 24, 30, 36, 48 and 60 h were measured; four time points, namely, 0, 12, 24 and 60 h, were selected for this differential proteomics study. The proteomic changes in the postharvest storage samples were studied by isobaric tags for relative and absolute quantification-coupled two-dimensional liquid chromatography-tandem mass spectrometry (2D LC–MS/MS). A total of 2,063 proteins were identified, and 192 differentially expressed proteins (DEPs), including 24 up-regulated proteins and 168 down-regulated proteins, were detected after 12 h of storage. After 24 h of storage, 234 DEPs, including 48 up-regulated and 186 down-regulated proteins, were observed, and after 60 h, 415 DEPs, including 65 up-regulated proteins and 350 down-regulated proteins, were observed. An in-depth data analysis showed that the DEPs participated in various cellular processes, particularly metabolic processes. In this study, we combined Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses, and the results focused on oxidative phosphorylation and ubiquitin mediated proteolysis pathways. In addition, sdh2, uba1 and ubc1 was confirmed by quantitative real-time polymerase chain reaction, and the results showed that the expression of these genes were consistent with their protein level. Based on the literature and our results, it is speculated that the identified DEPs, such as ATP1, SDH2, COR1, UBA1, COX4, UBC1 and SKP1 play a key role in the low-temperature autolysis of V. volvacea.

Investigating higher-order interactions in single-cell data with scHOT

Single-cell genomics has transformed our ability to examine cell fate choice. Examining cells along a computationally ordered 'pseudotime' offers the potential to unpick subtle changes in variability and covariation among key genes. We describe an approach, scHOT-single-cell higher-order testing-which provides a flexible and statistically robust framework for identifying changes in higher-order interactions among genes. scHOT can be applied for cells along a continuous trajectory or across space and accommodates various higher-order measurements including variability or correlation. We demonstrate the use of scHOT by studying coordinated changes in higher-order interactions during embryonic development of the mouse liver. Additionally, scHOT identifies subtle changes in gene-gene correlations across space using spatially resolved transcriptomics data from the mouse olfactory bulb. scHOT meaningfully adds to first-order differential expression testing and provides a framework for interrogating higher-order interactions using single-cell data.

Transcriptomic differences between male and female Trachycarpus fortunei

Trachycarpus fortunei (Hook.) is a typical dioecious plant, which has important economic value. There is currently no sex identification method for the early stages of T. fortunei growth. The aim of this study was to obtain expression and site differences between male and female T. fortunei transcriptomes. Using the Illumina sequencing platform, the transcriptomes of T. fortunei male and female plants were sequenced. By analyzing transcriptomic differences, the chromosomal helical binding protein (CHD1), serine/threonine protein kinase (STPK), cytochrome P450 716B1, and UPF0136 were found to be specifically expressed in T. fortunei males. After single nucleotide polymorphism (SNP) detection, a total of 12 male specific sites were found and the THUMP domain protein homologs were found to be male-biased expressed. Cytokinin dehydrogenase 6 (CKX6) was upregulated in male flowers and the lower concentrations of cytokinin (CTK) may be more conducive to male flower development. During new leaf growth, flavonoid and flavonol biosynthesis were initiated. Additionally, the flavonoids, 3',5'-hydroxylase (F3'5'H), flavonoids 3'-hydroxylase, were upregulated, which may cause the pale yellow phenotype. Based on these data, it can be concluded that inter-sex differentially expressed genes (DEGs) and specific SNP loci may be associated with sex determination in T. fortunei.

A bioinformatics investigation into molecular mechanism of Yinzhihuang granules for treating hepatitis B by network pharmacology and molecular docking verification

Yinzhihuang granules (YZHG) is a patented Chinese medicine for the treatment of hepatitis B. This study aimed to investigate the intrinsic mechanisms of YZHG in the treatment of hepatitis B and to provide new evidence and insights for its clinical application. The chemical compounds of YZHG were searched in the CNKI and PUBMED databases, and their putative targets were then predicted through a search of the SuperPred and Swiss Target Prediction databases. In addition, the targets of hepatitis B were obtained from TTD, PharmGKB and DisGeNET. The abovementioned data were visualized using Cytoscape 3.7.1, and network construction identified a total of 13 potential targets of YZHG in the treatment of hepatitis B. Molecular docking verification showed that CDK6, CDK2, TP53 and BRCA1 might be strongly correlated with hepatitis B treatment. Furthermore, GO and KEGG analyses indicated that the treatment of hepatitis B by YZHG might be related to positive regulation of transcription, positive regulation of gene expression, the hepatitis B pathway and the viral carcinogenesis pathway. Network pharmacology intuitively shows the multicomponent, multitarget and multichannel pharmacological effects of YZHG in the treatment of hepatitis B and provides a scientific basis for its mechanism of action.

Gold/alpha-lactalbumin nanoprobes for the imaging and treatment of breast cancer

Theranostic agents should ideally be renally cleared and biodegradable. Here, we report the synthesis, characterization and theranostic applications of fluorescent ultrasmall gold quantum clusters that are stabilized by the milk metalloprotein alpha-lactalbumin. We synthesized three types of these nanoprobes that together display fluorescence across the visible and near-infrared spectra when excited at a single wavelength through optical colour coding. In live tumour-bearing mice, the near-infrared nanoprobe generates contrast for fluorescence, X-ray computed tomography and magnetic resonance imaging, and exhibits long circulation times, low accumulation in the reticuloendothelial system, sustained tumour retention, insignificant toxicity and renal clearance. An intravenously administrated near-infrared nanoprobe with a large Stokes shift facilitated the detection and image-guided resection of breast tumours in vivo using a smartphone with modified optics. Moreover, the partially unfolded structure of alpha-lactalbumin in the nanoprobe helps with the formation of an anti-cancer lipoprotein complex with oleic acid that triggers the inhibition of the MAPK and PI3K-AKT pathways, immunogenic cell death and the recruitment of infiltrating macrophages. The biodegradability and safety profile of the nanoprobes make them suitable for the systemic detection and localized treatment of cancer.

Manumycin polyketides act as molecular glues between UBR7 and P53

Molecular glues are an intriguing therapeutic modality that harness small-molecules to induce interactions between proteins that typically do not interact. However, such molecules are rare and have been discovered fortuitously, thus limiting their potential as a general strategy for therapeutic intervention. We postulated that natural products bearing one or more electrophilic sites may be an unexplored source of new molecular glues, potentially acting through multi-covalent attachment. Using chemoproteomic platforms, we show that members of the manumycin family of polyketides, which bear multiple potentially reactive sites, target C374 of the putative E3 ligase UBR7 in breast cancer cells and engage in molecular glue interactions with the neo-substrate tumor-suppressor TP53, leading to p53 transcriptional activation and cell death. Our results reveal a novel anti-cancer mechanism of this natural product family and highlight the potential for combining chemoproteomics and multi-covalent natural products for the discovery of new molecular glues.

Transcriptional regulation of strigolactone signalling in Arabidopsis

Plant hormones known as strigolactones control plant development and interactions between host plants and symbiotic fungi or parasitic weeds^1-4. In Arabidopsis thaliana and rice, the proteins DWARF14 (D14), MORE AXILLARY GROWTH 2 (MAX2), SUPPRESSOR OF MAX2-LIKE 6, 7 and 8 (SMXL6, SMXL7 and SMXL8) and their orthologues form a complex upon strigolactone perception and play a central part in strigolactone signalling^5-10. However, whether and how strigolactones activate downstream transcription remains largely unknown. Here we use a synthetic strigolactone to identify 401 strigolactone-responsive genes in Arabidopsis, and show that these plant hormones regulate shoot branching, leaf shape and anthocyanin accumulation mainly through transcriptional activation of the BRANCHED 1, TCP DOMAIN PROTEIN 1 and PRODUCTION OF ANTHOCYANIN PIGMENT 1 genes. We find that SMXL6 targets 729 genes in the Arabidopsis genome and represses the transcription of SMXL6, SMXL7 and SMXL8 by binding directly to their promoters, showing that SMXL6 serves as an autoregulated transcription factor to maintain the homeostasis of strigolactone signalling. These findings reveal an unanticipated mechanism through which a transcriptional repressor of hormone signalling can directly recognize DNA and regulate transcription in higher plants.

Large-scale network analysis captures biological features of bacterial plasmids

Many bacteria can exchange genetic material through horizontal gene transfer (HGT) mediated by plasmids and plasmid-borne transposable elements. Here, we study the population structure and dynamics of over 10,000 bacterial plasmids, by quantifying their genetic similarities and reconstructing a network based on their shared k-mer content. We use a community detection algorithm to assign plasmids into cliques, which correlate with plasmid gene content, bacterial host range, GC content, and existing classifications based on replicon and mobility (MOB) types. Further analysis of plasmid population structure allows us to uncover candidates for yet undescribed replicon genes, and to identify transposable elements as the main drivers of HGT at broad phylogenetic scales. Our work illustrates the potential of network-based analyses of the bacterial 'mobilome' and opens up the prospect of a natural, exhaustive classification framework for bacterial plasmids.

Transcript expression-aware annotation improves rare variant interpretation

The acceleration of DNA sequencing in samples from patients and population studies has resulted in extensive catalogues of human genetic variation, but the interpretation of rare genetic variants remains problematic. A notable example of this challenge is the existence of disruptive variants in dosage-sensitive disease genes, even in apparently healthy individuals. Here, by manual curation of putative loss-of-function (pLoF) variants in haploinsufficient disease genes in the Genome Aggregation Database (gnomAD)1, we show that one explanation for this paradox involves alternative splicing of mRNA, which allows exons of a gene to be expressed at varying levels across different cell types. Currently, no existing annotation tool systematically incorporates information about exon expression into the interpretation of variants. We develop a transcript-level annotation metric known as the 'proportion expressed across transcripts', which quantifies isoform expression for variants. We calculate this metric using 11,706 tissue samples from the Genotype Tissue Expression (GTEx) project2 and show that it can differentiate between weakly and highly evolutionarily conserved exons, a proxy for functional importance. We demonstrate that expression-based annotation selectively filters 22.8% of falsely annotated pLoF variants found in haploinsufficient disease genes in gnomAD, while removing less than 4% of high-confidence pathogenic variants in the same genes. Finally, we apply our expression filter to the analysis of de novo variants in patients with autism spectrum disorder and intellectual disability or developmental disorders to show that pLoF variants in weakly expressed regions have similar effect sizes to those of synonymous variants, whereas pLoF variants in highly expressed exons are most strongly enriched among cases. Our annotation is fast, flexible and generalizable, making it possible for any variant file to be annotated with any isoform expression dataset, and will be valuable for the genetic diagnosis of rare diseases, the analysis of rare variant burden in complex disorders, and the curation and prioritization of variants in recall-by-genotype studies.

CRISPR-mediated gene targeting of CK1δ/ε leads to enhanced understanding of their role in endocytosis via phosphoregulation of GAPVD1

Human casein kinase 1 delta (CK1δ) and epsilon (CK1ε) are members of a conserved family of abundant, ubiquitously expressed serine/threonine kinases that regulate multiple cellular processes including circadian rhythm and endocytosis. Here, we have investigated the localization and interactomes of endogenously tagged CK1δ and CK1ε during interphase and mitosis. CK1δ and CK1ε localize to centrosomes throughout the cell cycle, and in interphase cells to the nucleus, and in both a diffuse and punctate pattern in the cytoplasm. Also, for the first time, they were detected at the midbody during cell division. Mass spectrometry analysis identified a total of 181 proteins co-purifying with a Venus multifunctional (VM)-tagged CK1δ and/or CK1ε. GTPase-activating protein and VPS9 domain-containing protein 1 (GAPVD1), a protein required for efficient endocytosis, was consistently one of the most abundant interacting partners. We demonstrate that GAPVD1 is a substrate of CK1δ/ε with up to 38 phosphorylated residues in vitro and in vivo. Wildtype and a phosphomimetic mutant of GAPVD1, but not a phospho-ablating mutant, were able to rescue defects in transferrin and EGF internalization caused by loss of endogenous GAPVD1. Our results indicate that GAPVD1 is an important interacting partner and substrate of CK1δ/ε for endocytosis.

Identification, quantification and bioinformatic analysis of RNA-dependent proteins by RNase treatment and density gradient ultracentrifugation using R-DeeP

Analysis of RNA-protein complexes is central to understanding the molecular circuitry governing cellular processes. In recent years, several proteome-wide studies have been dedicated to the identification of RNA-binding proteins. Here, we describe in detail R-DeeP, an approach built on RNA dependence, defined as the ability of a protein to engage in protein complexes only in the presence of RNA, involving direct or indirect interaction with RNA. This approach provides-for the first time, to our knowledge-quantitative information on the fraction of a protein associated with RNA-protein complexes. R-DeeP is independent of any potentially biased purification procedures. It is based on cellular lysate fractionation by density gradient ultracentrifugation and subsequent analysis by proteome-wide mass spectrometry (MS) or individual western blotting. The comparison of lysates with and without previous RNase treatment enables the identification of differences in the apparent molecular weight and, hence, the size of the complexes. In combination with information from databases of protein-protein complexes, R-DeeP facilitates the computational reconstruction of protein complexes from proteins migrating in the same fraction. In addition, we developed a pipeline for the statistical analysis of the MS dataset to automatically identify RNA-dependent proteins (proteins whose interactome depends on RNA). With this protocol, the individual analysis of proteins of interest by western blotting can be completed within 1-2 weeks. For proteome-wide studies, additional time is needed for the integration of the proteomic and statistical analyses. In the future, R-DeeP can be extended to other fractionation techniques, such as chromatography.

A census of pathway maps in cancer systems biology

A key goal of cancer systems biology is to use big data to elucidate the molecular networks by which cancer develops. However, to date there has been no systematic evaluation of how far these efforts have progressed. In this Analysis, we survey six major systems biology approaches for mapping and modelling cancer pathways with attention to how well their resulting network maps cover and enhance current knowledge. Our sample of 2,070 systems biology maps captures all literature-curated cancer pathways with significant enrichment, although the strong tendency is for these maps to recover isolated mechanisms rather than entire integrated processes. Systems biology maps also identify previously underappreciated functions, such as a potential role for human papillomavirus-induced chromosomal alterations in ovarian tumorigenesis, and they add new genes to known cancer pathways, such as those related to metabolism, Hippo signalling and immunity. Notably, we find that many cancer networks have been provided only in journal figures and not for programmatic access, underscoring the need to deposit network maps in community databases to ensure they can be readily accessed. Finally, few of these findings have yet been clinically translated, leaving ample opportunity for future translational studies. Periodic surveys of cancer pathway maps, such as the one reported here, are critical to assess progress in the field and identify underserved areas of methodology and cancer biology.

Integrative differential expression and gene set enrichment analysis using summary statistics for scRNA-seq studies

Differential expression (DE) analysis and gene set enrichment (GSE) analysis are commonly applied in single cell RNA sequencing (scRNA-seq) studies. Here, we develop an integrative and scalable computational method, iDEA, to perform joint DE and GSE analysis through a hierarchical Bayesian framework. By integrating DE and GSE analyses, iDEA can improve the power and consistency of DE analysis and the accuracy of GSE analysis. Importantly, iDEA uses only DE summary statistics as input, enabling effective data modeling through complementing and pairing with various existing DE methods. We illustrate the benefits of iDEA with extensive simulations. We also apply iDEA to analyze three scRNA-seq data sets, where iDEA achieves up to five-fold power gain over existing GSE methods and up to 64% power gain over existing DE methods. The power gain brought by iDEA allows us to identify many pathways that would not be identified by existing approaches in these data.

Differential expression (DE) and gene set enrichment (GSE) analysis tend to be carried out separately. Here, the authors present iDEA (integrative Differential expression and gene set Enrichment Analysis) for the analysis of scRNAseq data which uses a Baysian approach to jointly model DE and GSE for improved power in both tasks.

Sjögren syndrome/scleroderma autoantigen 1 is a direct Tankyrase binding partner in cancer cells

Sjögren syndrome/scleroderma autoantigen 1 (SSSCA1) was first described as an auto-antigen over-expressed in Sjögren’s syndrome and in scleroderma patients. SSSCA1 has been linked to mitosis and centromere association and as a potential marker candidate in diverse solid cancers. Here we characterize SSSCA1 for the first time, to our knowledge, at the molecular, structural and subcellular level. We have determined the crystal structure of a zinc finger fold, a zinc ribbon domain type 2 (ZNRD2), at 2.3 Å resolution. We show that the C-terminal domain serves a dual function as it both behaves as the interaction site to Tankyrase 1 (TNKS1) and as a nuclear export signal. We identify TNKS1 as a direct binding partner of SSSCA1, map the binding site to TNKS1 ankyrin repeat cluster 2 (ARC2) and thus define a new binding sequence. We experimentally verify and map a new nuclear export signal sequence in SSSCA1.

Perdreau-Dahl et al. systematically characterise Sjögren syndrome/scleroderma autoantigen 1 at the molecular, structural and subcellular level. They show that the C-terminal domain serves a dual function as it both behaves as the interaction site to Tankyrase 1 and as a nuclear export signal.

Acid ceramidase of macrophages traps herpes simplex virus in multivesicular bodies and protects from severe disease

Macrophages have important protective functions during infection with herpes simplex virus type 1 (HSV-1). However, molecular mechanisms that restrict viral propagation and protect from severe disease are unclear. Here we show that macrophages take up HSV-1 via endocytosis and transport the virions into multivesicular bodies (MVBs). In MVBs, acid ceramidase (aCDase) converts ceramide into sphingosine and increases the formation of sphingosine-rich intraluminal vesicles (ILVs). Once HSV-1 particles reach MVBs, sphingosine-rich ILVs bind to HSV-1 particles, which restricts fusion with the limiting endosomal membrane and prevents cellular infection. Lack of aCDase in macrophage cultures or in Asah1^-/- mice results in replication of HSV-1 and Asah1^-/- mice die soon after systemic or intravaginal inoculation. The treatment of macrophages with sphingosine enhancing compounds blocks HSV-1 propagation, suggesting a therapeutic potential of this pathway. In conclusion, aCDase loads ILVs with sphingosine, which prevents HSV-1 capsids from penetrating into the cytosol.

Transcriptional diversity and bioenergetic shift in human breast cancer metastasis revealed by single-cell RNA sequencing

Although metastasis remains the cause of most cancer-related mortality, mechanisms governing seeding in distal tissues are poorly understood. Here, we establish a robust method for the identification of global transcriptomic changes in rare metastatic cells during seeding using single-cell RNA sequencing and patient-derived-xenograft models of breast cancer. We find that both primary tumours and micrometastases display transcriptional heterogeneity but micrometastases harbour a distinct transcriptome program conserved across patient-derived-xenograft models that is highly predictive of poor survival of patients. Pathway analysis revealed mitochondrial oxidative phosphorylation as the top pathway upregulated in micrometastases, in contrast to higher levels of glycolytic enzymes in primary tumour cells, which we corroborated by flow cytometric and metabolomic analyses. Pharmacological inhibition of oxidative phosphorylation dramatically attenuated metastatic seeding in the lungs, which demonstrates the functional importance of oxidative phosphorylation in metastasis and highlights its potential as a therapeutic target to prevent metastatic spread in patients with breast cancer.

Gene regulatory network analysis with drug sensitivity reveals synergistic effects of combinatory chemotherapy in gastric cancer

The combination of docetaxel, cisplatin, and fluorouracil (DCF) is highly synergistic in advanced gastric cancer. We aimed to explain these synergistic effects at the molecular level. Thus, we constructed a weighted correlation network using the differentially expressed genes between Stage I and IV gastric cancer based on The Cancer Genome Atlas (TCGA), and three modules were derived. Next, we investigated the correlation between the eigengene of the expression of the gene network modules and the chemotherapeutic drug response to DCF from the Genomics of Drug Sensitivity in Cancer (GDSC) database. The three modules were associated with functions related to cell migration, angiogenesis, and the immune response. The eigengenes of the three modules had a high correlation with DCF (−0.41, −0.40, and −0.15). The eigengenes of the three modules tended to increase as the stage increased. Advanced gastric cancer was affected by the interaction the among modules with three functions, namely cell migration, angiogenesis, and the immune response, all of which are related to metastasis. The weighted correlation network analysis model proved the complementary effects of DCF at the molecular level and thus, could be used as a unique methodology to determine the optimal combination of chemotherapy drugs for patients with gastric cancer.

Toxicological response of the model fungus Saccharomyces cerevisiae to different concentrations of commercial graphene nanoplatelets

Graphene nanomaterials have attracted a great interest during the last years for different applications, but their possible impact on different biological systems remains unclear. Here, an assessment to understand the toxicity of commercial polycarboxylate functionalized graphene nanoplatelets (GN) on the unicellular fungal model Saccharomyces cerevisiae was performed. While cell proliferation was not negatively affected even in the presence of 800 mg L-1 of the nanomaterial for 24 hours, oxidative stress was induced at a lower concentration (160 mg L-1), after short exposure periods (2 and 4 hours). No DNA damage was observed under a comet assay analysis under the studied conditions. In addition, to pinpoint the molecular mechanisms behind the early oxidative damage induced by GN and to identify possible toxicity pathways, the transcriptome of S. cerevisiae exposed to 160 and 800 mg L-1 of GN was studied. Both GN concentrations induced expression changes in a common group of genes (337), many of them related to the fungal response to reduce the nanoparticles toxicity and to maintain cell homeostasis. Also, a high number of genes were only differentially expressed in the GN800 condition (3254), indicating that high GN concentrations can induce severe changes in the physiological state of the yeast.

Restricted immunological and cellular pathways are shared by murine models of chronic alcohol consumption

Murine models of chronic alcohol consumption are frequently used to investigate alcoholic liver injury and define new therapeutic targets. Lieber-DeCarli diet (LD) and Meadows-Cook diet (MC) are the most accepted models of chronic alcohol consumption. It is unclear how similar these models are at the cellular, immunologic, and transcriptome levels. We investigated the common and specific pathways of LD and MC models. Livers from LD and MC mice were subjected to histologic changes, hepatic leukocyte population, hepatic transcripts level related to leukocyte recruitment, and hepatic RNA-seq analysis. Cross-species comparison was performed using the alcoholic liver disease (ALD) transcriptomic public dataset. Despite LD mice have increased liver injury and steatosis by alcohol exposure, the number of CD45⁺ cells were reduced. Opposite, MC mice have an increased number of monocytes/liver by alcohol. The pattern of chemokine gradient, adhesion molecules, and cytokine transcripts is highly specific for each model, not shared with advanced human alcoholic liver disease. Moreover, hepatic RNA-seq revealed a limited and restricted number of shared genes differentially changed by alcohol exposure in these 2 models. Thus, mechanisms involved in alcohol tissue injury are model-dependent at multiple levels and raise the consideration of significant pathophysiological diversity of human alcoholic liver injury.

Arginine π-stacking drives binding to fibrils of the Alzheimer protein Tau

Aggregation of the Tau protein into fibrils defines progression of neurodegenerative diseases, including Alzheimer’s Disease. The molecular basis for potentially toxic reactions of Tau aggregates is poorly understood. Here we show that π-stacking by Arginine side-chains drives protein binding to Tau fibrils. We mapped an aggregation-dependent interaction pattern of Tau. Fibrils recruit specifically aberrant interactors characterised by intrinsically disordered regions of atypical sequence features. Arginine residues are key to initiate these aberrant interactions. Crucial for scavenging is the guanidinium group of its side chain, not its charge, indicating a key role of π-stacking chemistry for driving aberrant fibril interactions. Remarkably, despite the non-hydrophobic interaction mode, the molecular chaperone Hsp90 can modulate aberrant fibril binding. Together, our data present a molecular mode of action for derailment of protein-protein interaction by neurotoxic fibrils.

Tau fibril formation is a hallmark of Alzheimer’s disease. Here the authors reveal an aggregation-dependent protein interaction pattern of Tau and further show that π-stacking of the arginine side-chains drives aberrant protein binding to Tau fibrils.