Single-nucleotide variant calling in single-cell sequencing data with Monopogen.

Single-cell omics technologies enable molecular characterization of diverse cell types and states, but how the resulting transcriptional and epigenetic profiles depend on the cell's genetic background remains understudied. We describe Monopogen, a computational tool to detect single-nucleotide variants (SNVs) from single-cell sequencing data. Monopogen leverages linkage disequilibrium from external reference panels to identify germline SNVs and detects putative somatic SNVs using allele cosegregating patterns at the cell population level. It can identify 100 K to 3 M germline SNVs achieving a genotyping accuracy of 95%, together with hundreds of putative somatic SNVs. Monopogen-derived genotypes enable global and local ancestry inference and identification of admixed samples. It identifies variants associated with cardiomyocyte metabolic levels and epigenomic programs. It also improves putative somatic SNV detection that enables clonal lineage tracing in primary human clonal hematopoiesis. Monopogen brings together population genetics, cell lineage tracing and single-cell omics to uncover genetic determinants of cellular processes.

A benchmark for RNA-seq deconvolution analysis under dynamic testing environments.

BACKGROUND: Deconvolution analyses have been widely used to track compositional alterations of cell types in gene expression data. Although a large number of novel methods have been developed, due to a lack of understanding of the effects of modeling assumptions and tuning parameters, it is challenging for researchers to select an optimal deconvolution method suitable for the targeted biological conditions. RESULTS: To systematically reveal the pitfalls and challenges of deconvolution analyses, we investigate the impact of several technical and biological factors including simulation model, quantification unit, component number, weight matrix, and unknown content by constructing three benchmarking frameworks. These frameworks cover comparative analysis of 11 popular deconvolution methods under 1766 conditions. CONCLUSIONS: We provide new insights to researchers for future application, standardization, and development of deconvolution tools on RNA-seq data.

Losing Dnmt3a dependent methylation in inhibitory neurons impairs neural function by a mechanism impacting Rett syndrome.

Methylated cytosine is an effector of epigenetic gene regulation. In the brain, Dnmt3a is the sole 'writer' of atypical non-CpG methylation (mCH), and MeCP2 is the only known 'reader' for mCH. We asked if MeCP2 is the sole reader for Dnmt3a dependent methylation by comparing mice lacking either protein in GABAergic inhibitory neurons. Loss of either protein causes overlapping and distinct features from the behavioral to molecular level. Loss of Dnmt3a causes global loss of mCH and a subset of mCG sites resulting in more widespread transcriptional alterations and severe neurological dysfunction than MeCP2 loss. These data suggest that MeCP2 is responsible for reading only part of the Dnmt3a dependent methylation in the brain. Importantly, the impact of MeCP2 on genes differentially expressed in both models shows a strong dependence on mCH, but not Dnmt3a dependent mCG, consistent with mCH playing a central role in the pathogenesis of Rett Syndrome.

Forniceal deep brain stimulation induces gene expression and splicing changes that promote neurogenesis and plasticity.

Clinical trials are currently underway to assess the efficacy of forniceal deep brain stimulation (DBS) for improvement of memory in Alzheimer's patients, and forniceal DBS has been shown to improve learning and memory in a mouse model of Rett syndrome (RTT), an intellectual disability disorder caused by loss-of-function mutations in MECP2. The mechanism of DBS benefits has been elusive, however, so we assessed changes in gene expression, splice isoforms, DNA methylation, and proteome following acute forniceal DBS in wild-type mice and mice lacking Mecp2. We found that DBS upregulates genes involved in synaptic function, cell survival, and neurogenesis and normalized expression of ~25% of the genes altered in Mecp2-null mice. Moreover, DBS induced expression of 17-24% of the genes downregulated in other intellectual disability mouse models and in post-mortem human brain tissue from patients with Major Depressive Disorder, suggesting forniceal DBS could benefit individuals with a variety of neuropsychiatric disorders.

Comprehensive evaluation of RNA-seq quantification methods for linearity.

BACKGROUND: Deconvolution is a mathematical process of resolving an observed function into its constituent elements. In the field of biomedical research, deconvolution analysis is applied to obtain single cell-type or tissue specific signatures from a mixed signal and most of them follow the linearity assumption. Although recent development of next generation sequencing technology suggests RNA-seq as a fast and accurate method for obtaining transcriptomic profiles, few studies have been conducted to investigate best RNA-seq quantification methods that yield the optimum linear space for deconvolution analysis. RESULTS: Using a benchmark RNA-seq dataset, we investigated the linearity of abundance estimated from seven most popular RNA-seq quantification methods both at the gene and isoform levels. Linearity is evaluated through parameter estimation, concordance analysis and residual analysis based on a multiple linear regression model. Results show that count data gives poor parameter estimations, large intercepts and high inter-sample variability; while TPM value from Kallisto and Salmon shows high linearity in all analyses. CONCLUSIONS: Salmon and Kallisto TPM data gives the best fit to the linear model studied. This suggests that TPM values estimated from Salmon and Kallisto are the ideal RNA-seq measurements for deconvolution studies.

Codon-level co-occurrences of germline variants and somatic mutations in cancer are rare but often lead to incorrect variant annotation and underestimated impact prediction.

Cancer cells explore a broad mutational landscape, bringing the possibility that tumor-specific somatic mutations could fall in the same codons as germline SNVs and leverage their presence to produce substitutions with a larger impact on protein function. While multiple, temporally consecutive mutations to the same codon have in the past been detected in the germline, this phenomenon has not yet been explored in the context of germline-somatic variant co-occurrences during cancer development. We examined germline context at somatic mutation sites for 1395 patients across four cancer cohorts (breast, skin, colon, and head and neck) and found 392 codon-level co-occurrences between germline and somatic variants, including over a dozen in well-known cancer genes. We found that for the majority of these co-occurrence events, traditional somatic calling led to an inaccurate representation of the protein site and a significantly lower predicted impact on protein fitness. We conclude that these events often lead to imprecise annotation of somatic variants but do not appear to be a frequent source of driver events during cancer development.

Mitochondrial electron transport chain identified as a novel molecular target of SPIO nanoparticles mediated cancer-specific cytotoxicity.

Superparamagnetic iron oxide nanoparticles (SPIONs) are highly cytotoxic and target cancer cells with high specificity; however, the mechanism by which SPIONs induce cancer cell-specific cytotoxicity remains unclear. Herein, the molecular mechanism of SPION-induced cancer cell-specific cytotoxicity to cancer cells is clarified through DNA microarray and bioinformatics analyses. SPIONs can interference with the mitochondrial electron transport chain (METC) in cancer cells, which further affects the production of ATP, mitochondrial membrane potential, and microdistribution of calcium, and induces cell apoptosis. Additionally, SPIONs induce the formation of reactive oxygen species in mitochondria; these reactive oxygen species trigger cancer-specific cytotoxicity due to the lower antioxidative capacity of cancer cells. Moreover, the DNA microarray and gene ontology analyses revealed that SPIONs elevate the expression of metallothioneins in both normal and cancer cells but decrease the expression of METC genes in cancer cells. Overall, these results suggest that SPIONs induce cancer cell death by targeting the METC, which is helpful for designing anti-cancer nanotheranostics and evaluating the safety of future nanomedicines.

The resurrection genome of Boea hygrometrica: A blueprint for survival of dehydration.

"Drying without dying" is an essential trait in land plant evolution. Unraveling how a unique group of angiosperms, the Resurrection Plants, survive desiccation of their leaves and roots has been hampered by the lack of a foundational genome perspective. Here we report the ∼1,691-Mb sequenced genome of Boea hygrometrica, an important resurrection plant model. The sequence revealed evidence for two historical genome-wide duplication events, a compliment of 49,374 protein-coding genes, 29.15% of which are unique (orphan) to Boea and 20% of which (9,888) significantly respond to desiccation at the transcript level. Expansion of early light-inducible protein (ELIP) and 5S rRNA genes highlights the importance of the protection of the photosynthetic apparatus during drying and the rapid resumption of protein synthesis in the resurrection capability of Boea. Transcriptome analysis reveals extensive alternative splicing of transcripts and a focus on cellular protection strategies. The lack of desiccation tolerance-specific genome organizational features suggests the resurrection phenotype evolved mainly by an alteration in the control of dehydration response genes.

3-Anhydro-6-hydroxy-ophiobolin A, a fungal sesterterpene from Bipolaris oryzae induced autophagy and promoted the degradation of α-synuclein in PC12 cells.

Autophagy is defined as an evolutionarily conserved process responsible for degradation of the cytoplasmic components including protein aggregates via the lysosomal machinery. Increasing evidence has linked defective autophagic degradation of protein aggregates with the pathogenesis of neurodegenerative disorders, and it is suggested that promotion of autophagy is regarded as a potential therapeutic for these diseases including Parkinson's disease (PD). Here we identified, 3-anhydro-6-hydroxy-ophiobolin A (X15-2), an ophiobolin derivative from Bipolaris oryzae that can strongly induce autophagic degradation of α-synuclein, the major constituent of Lewy bodies. We showed that X15-2 induced autophagy is dependent on both Beclin1 and Beclin2. Knockout of ATG5 by CRISPER/Cas9 prevented X15-2 induced autophagy and degradation of α-synuclein. Mechanistically, we showed that X15-2 induces ROS and the activation of JNK signaling for the autophagic degradation of α-synuclein in PC12 cells.

ADReCS: an ontology database for aiding standardization and hierarchical classification of adverse drug reaction terms.

Adverse drug reactions (ADRs) are noxious and unexpected effects during normal drug therapy. They have caused significant clinical burden and been responsible for a large portion of new drug development failure. Molecular understanding and in silico evaluation of drug (or candidate) safety in laboratory is thus so desired, and unfortunately has been largely hindered by misuse of ADR terms. The growing impact of bioinformatics and systems biology in toxicological research also requires a specialized ADR term system that works beyond a simple glossary. Adverse Drug Reaction Classification System (ADReCS; http://bioinf.xmu.edu.cn/ADReCS) is a comprehensive ADR ontology database that provides not only ADR standardization but also hierarchical classification of ADR terms. The ADR terms were pre-assigned with unique digital IDs and at the same time were well organized into a four-level ADR hierarchy tree for building an ADR-ADR relation. Currently, the database covers 6544 standard ADR terms and 34,796 synonyms. It also incorporates information of 1355 single active ingredient drugs and 134,022 drug-ADR pairs. In summary, ADReCS offers an opportunity for direct computation on ADR terms and also provides clues to mining common features underlying ADRs.

The BCL2L1 and PGAM5 axis defines hypoxia-induced receptor-mediated mitophagy.

Receptor-mediated mitophagy is one of the major mechanisms of mitochondrial quality control essential for cell survival. We previously have identified FUNDC1 as a mitophagy receptor for selectively removing damaged mitochondria in mammalian systems. A critical unanswered question is how receptor-mediated mitophagy is regulated in response to cellular and environmental cues. Here, we report the striking finding that BCL2L1/Bcl-xL, but not BCL2, suppresses mitophagy mediated by FUNDC1 through its BH3 domain. Mechanistically, we demonstrate that BCL2L1, but not BCL2, interacts with and inhibits PGAM5, a mitochondrially localized phosphatase, to prevent the dephosphorylation of FUNDC1 at serine 13 (Ser13), which activates hypoxia-induced mitophagy. Our results showed that the BCL2L1-PGAM5-FUNDC1 axis is critical for receptor-mediated mitophagy in response to hypoxia and that BCL2L1 possesses unique functions distinct from BCL2.

A small natural molecule promotes mitochondrial fusion through inhibition of the deubiquitinase USP30.

Mitochondrial fusion is a highly coordinated process that mixes and unifies the mitochondrial compartment for normal mitochondrial functions and mitochondrial DNA inheritance. Dysregulated mitochondrial fusion causes mitochondrial fragmentation, abnormal mitochondrial physiology and inheritance, and has been causally linked with a number of neuronal diseases. Here, we identified a diterpenoid derivative 15-oxospiramilactone (S3) that potently induced mitochondrial fusion to restore the mitochondrial network and oxidative respiration in cells that are deficient in either Mfn1 or Mfn2. A mitochondria-localized deubiquitinase USP30 is a target of S3. The inhibition of USP30 by S3 leads to an increase of non-degradative ubiquitination of Mfn1/2, which enhances Mfn1 and Mfn2 activity and promotes mitochondrial fusion. Thus, through the use of an inhibitor of USP30, our study uncovers an unconventional function of non-degradative ubiquitination of Mfns in promoting mitochondrial fusion.

Establishment of a human colorectal cancer cell line P6C with stem cell properties and resistance to chemotherapeutic drugs.

AIM: Cancer stem cells have the capacity to initiate and sustain tumor growth. In this study, we established a CD44(+) colorectal cancer stem cell line with particular emphasis on its self-renewal capacity, enhanced tumor initiation and drug resistance. METHODS: Fresh colon cancer and paired normal colon tissues were collected from 13 patients who had not received chemotherapy or radiotherapy prior to surgery. Among the 6 single-cell derived clones, only the P6C cell line was cultured for more than 20 passages in serial culture and formed holoclones with high efficiency, and then the stemness gene expression, colony formation, tumorigenicity and drug sensitivities of the P6C cell line were examined. RESULTS: Stemness proteins, including c-Myc, Oct3/4, Nanog, Lgr5, and SOX2, were highly expressed in the P6C cell line. Oct3/4-positive P6C cells mostly generated holoclones through symmetric division, while a small number of P6C cells generated meroclones through asymmetric division. P6C cells stably expressed CD44 and possessed a high capacity to form tumor spheres. A single cell-derived sphere was capable of generating xenograft tumors in nude mice. Compared to SW480 and HCT116 colorectal cancer cells, P6C cells were highly resistant to Camptothecin and 5-fluorouracil, the commonly used chemotherapeutic agents to treat colorectal cancers. CONCLUSION: We established a colorectal cancer stem cell line P6C with a high tumorigenic capacity and the characteristics of normal stem cells. It will benefit the mechanistic studies on cancer stem cells and the development of drugs that specifically target the cancer stem cells.

Reciprocal interactions between tumor-associated macrophages and CD44-positive cancer cells via osteopontin/CD44 promote tumorigenicity in colorectal cancer.

PURPOSE: CD44 is of functional importance for tumor initiation and progression in colorectal cancer, but how this molecule benefits cancer cells from the tumor microenvironment, especially tumor-associated macrophages (TAM), remains poorly defined. EXPERIMENTAL DESIGN: In vivo tumorigenic assays were conducted to assess the role of murine TAMs in the tumorigenesis of human colorectal cancer cells. Both in vitro and in vivo osteopontin (OPN) expression levels in TAMs were examined by immunohistochemistry, quantitative PCR, and Western blotting. Soft agar colony formation assays were used to estimate the clonogenicity of colorectal cancer cells that had received different treatments. The relationships between the expression levels of OPN, CD44v6, and CD68 and clinical prognosis were evaluated by tissue microarray analysis. RESULTS: We found that macrophages, when coinjected or cocultured with CD44-positive colorectal cancer cells, were able to produce higher levels of OPN, which in turn facilitated the tumorigenicity and clonogenicity of the colorectal cancer cells. The knockdown of CD44 or treatment with blocking antibodies to CD44 attenuated OPN secretion. OPN, through binding to its receptor CD44, activated c-jun-NH(2)-kinase signaling and promoted the clonogenicity of colorectal cancer cells. Moreover, tissue microarray data have shown that OPN expression, in combination with CD44v6, has a negative correlation with colorectal cancer patient survival. CONCLUSIONS: These results suggest that the OPN-CD44 interaction is important for colorectal cancer progression and could serve as a potential therapeutic target for the treatment of colorectal cancer.

Mitochondrial outer-membrane protein FUNDC1 mediates hypoxia-induced mitophagy in mammalian cells.

Accumulating evidence has shown that dysfunctional mitochondria can be selectively removed by mitophagy. Dysregulation of mitophagy is implicated in the development of neurodegenerative disease and metabolic disorders. How individual mitochondria are recognized for removal and how this process is regulated remain poorly understood. Here we report that FUNDC1, an integral mitochondrial outer-membrane protein, is a receptor for hypoxia-induced mitophagy. FUNDC1 interacted with LC3 through its typical LC3-binding motif Y(18)xxL(21), and mutation of the LC3-interaction region impaired its interaction with LC3 and the subsequent induction of mitophagy. Knockdown of endogenous FUNDC1 significantly prevented hypoxia-induced mitophagy, which could be reversed by the expression of wild-type FUNDC1, but not LC3-interaction-deficient FUNDC1 mutants. Mechanistic studies further revealed that hypoxia induced dephosphorylation of FUNDC1 and enhanced its interaction with LC3 for selective mitophagy. Our findings thus offer insights into mitochondrial quality control in mammalian cells.

Caspase cleavage of cytochrome c1 disrupts mitochondrial function and enhances cytochrome c release.

Mitochondrial catastrophe can be the cause or consequence of apoptosis and is associated with a number of pathophysiological conditions. The exact relationship between mitochondrial catastrophe and caspase activation is not completely understood. Here we addressed the underlying mechanism, explaining how activated caspase could feedback to attack mitochondria to amplify further cytochrome c (cyto.c) release. We discovered that cytochrome c1 (cyto.c1) in the bc1 complex of the mitochondrial respiration chain was a novel substrate of caspase 3 (casp.3). We found that cyto.c1 was cleaved at the site of D106, which is critical for binding with cyto.c, following apoptotic stresses or targeted expression of casp.3 into the mitochondrial intermembrane space. We demonstrated that this cleavage was closely linked with further cyto.c release and mitochondrial catastrophe. These mitochondrial events could be effectively blocked by expressing non-cleavable cyto.c1 (D106A) or by caspase inhibitor z-VAD-fmk. Our results demonstrate that the cleavage of cyto.c1 represents a critical step for the feedback amplification of cyto.c release by caspases and subsequent mitochondrial catastrophe.

Natural diterpenoid compound elevates expression of Bim protein, which interacts with antiapoptotic protein Bcl-2, converting it to proapoptotic Bax-like molecule.

Overwhelming evidence indicates that Bax and Bak are indispensable for mediating cytochrome c release from mitochondria during apoptosis. Here we report a Bax/Bak-independent mechanism of cytochrome c release and apoptosis. We identified a natural diterpenoid compound that induced apoptosis in bax/bak double knock-out murine embryonic fibroblasts and substantially reduced the tumor growth from these cells implanted in mice. Treatment with the compound significantly increased expression of Bim, which migrated to mitochondria, altering the conformation of and forming oligomers with resident Bcl-2 to induce cytochrome c release and caspase activation. Importantly, purified Bim and Bcl-2 proteins cooperated to permeabilize a model mitochondrial outer membrane; this was accompanied by oligomerization of these proteins and deep embedding of Bcl-2 in the membrane. Therefore, the diterpenoid compound induces a structural and functional conversion of Bcl-2 through Bim to permeabilize the mitochondrial outer membrane, thereby inducing apoptosis independently of Bax and Bak. Because Bcl-2 family proteins play important roles in cancer development and relapse, this novel cell death mechanism can be explored for developing more effective anticancer therapeutics.

Morphine induces Beclin 1- and ATG5-dependent autophagy in human neuroblastoma SH-SY5Y cells and in the rat hippocampus.

Chronic exposure to morphine can induce drug addiction and neural injury, but the exact mechanism is not fully understood. Here we show that morphine induces autophagy in neuroblastoma SH-SY5Y cells and in the rat hippocampus. Pharmacological approach shows that this effect appears to be mediated by PTX-sensitive G protein-coupled receptors signaling cascade. Morphine increases Beclin 1 expression and reduces the interaction between Beclin 1 and Bcl-2, thus releasing Beclin 1 for its pro-autophagic activity. Bcl-2 overexpression inhibits morphine-induced autophagy, whereas knockdown of Beclin 1 or knockout of ATG5 prevents morphine-induced autophagy. In addition, chronic treatment with morphine induces cell death, which is increased by autophagy inhibition through Beclin 1 RNAi. Our data are the first to reveal that Beclin 1 and ATG5 play key roles in morphine-induced autophagy, which may contribute to morphine-induced neuronal injury.

Beclin 1 cleavage by caspase-3 inactivates autophagy and promotes apoptosis.

Autophagy and apoptosis are both highly regulated biological processes that play essential roles in tissue homeostasis, development and diseases. Autophagy is also described as a mechanism of death pathways, however, the precise mechanism of how autophagy links to cell death remains to be fully understood. Beclin 1 is a dual regulator for both autophagy and apoptosis. In this study we found that Beclin 1 was a substrate of caspase-3 with two cleavage sites at positions 124 and 149, respectively. Furthermore, the autophagosome formation occurred, followed by the appearance of morphological hallmarks of apoptosis after staurosporine treatment. The cleavage products of Beclin 1 reduced autophagy and promoted apoptosis in HeLa cells and the cells in which Beclin 1 was stably knocked down by specific shRNA. In addition, the cleavage of Beclin 1 resulted in abrogating the interaction between Bcl-2 with Beclin 1, which could be blocked by z-VAD-fmk. Thus, our results suggest that the cleavage of Beclin 1 by caspase-3 may contribute to inactivate autophagy leading towards augmented apoptosis.

Targeted expression of uncoupling protein 2 to mouse liver increases the susceptibility to lipopolysaccharide/galactosamine-induced acute liver injury.

UNLABELLED: Normal hepatocytes do not express endogenous uncoupling protein 2 (UCP2) in adult liver, although Kupffer cells do, and it is strikingly induced in hepatocytes in steatotic liver and obese conditions. However, the direct link of UCP2 with the pathogenic development of liver diseases and liver injury remains elusive. Here we report that targeted expression of UCP2 to mouse liver increases susceptibility to acute liver injury induced by lipopolysaccharide (LPS) and galactosamine (GalN). UCP2 appears to enhance proton leak, leading to mild uncoupling in a guanosine diphosphate-repressible manner. Indeed, mitochondria from the genetically manipulated mouse liver have increased state 4 respiration, lower respiratory control ratio, and reduced adenosine triphosphate (ATP) levels, which altered mitochondrial physiology. To address the underlying mechanism of how UCP2 and the reduced energy coupling efficiency enhance cell death in mouse liver, we show that the reduced ATP levels lead to activation of 5'AMP-activated protein kinase (AMPK) and its downstream effector, c-Jun N-terminal kinase; thus, the increased sensitivity toward LPS/GalN-induces apoptosis. Importantly, we show that inhibition of UCP2 activity by its pharmacological inhibitor genipin prevents LPS/GalN-induced ATP reduction, AMPK activation, and apoptosis. Also, inhibition of ATP production by oligomycin promotes LPS/GalN-induced cell death both in vivo and in vitro. CONCLUSION: Our results clearly show that targeted expression of UCP2 in liver may result in compromised mitochondrial physiology that contributes to enhanced cell death and suggests a potential role of UCP2 in the development of liver diseases.