ABSTRACT
Emerging evidence supports that mitochondrial dysfunction contributes to systemic lupus erythematosus (SLE) pathogenesis. Here we show that programmed mitochondrial removal, a hallmark of mammalian erythropoiesis, is defective in SLE. Specifically, we demonstrate that during human erythroid cell maturation, a hypoxia-inducible factor (HIF)-mediated metabolic switch is responsible for the activation of the ubiquitin-proteasome system (UPS), which precedes and is necessary for the autophagic removal of mitochondria. A defect in this pathway leads to accumulation of red blood cells (RBCs) carrying mitochondria (Mito+ RBCs) in SLE patients and in correlation with disease activity. Antibody-mediated internalization of Mito+ RBCs induces type I interferon (IFN) production through activation of cGAS in macrophages. Accordingly, SLE patients carrying both Mito+ RBCs and opsonizing antibodies display the highest levels of blood IFN-stimulated gene (ISG) signatures, a distinctive feature of SLE.
Subject(s)
Interferon Type I/metabolism , Lupus Erythematosus, Systemic/metabolism , Mitochondria/metabolism , Myeloid Cells/metabolism , Adolescent , Basic Helix-Loop-Helix Transcription Factors/metabolism , Child , Child, Preschool , Erythroblasts/metabolism , Erythroblasts/ultrastructure , Erythrocytes/metabolism , Erythropoiesis , Humans , Mitophagy , Proteasome Endopeptidase Complex/metabolism , Ubiquitin/metabolismABSTRACT
Opsonization of red blood cells that retain mitochondria (Mito+ RBCs), a feature of systemic lupus erythematosus (SLE), triggers type I interferon (IFN) production in macrophages. We report that monocytes (Mos) co-produce IFN and mature interleukin-1ß (mIL-1ß) upon Mito+ RBC opsonization. IFN expression depended on cyclic GMP-AMP synthase (cGAS) and RIG-I-like receptors' (RLRs) sensing of Mito+ RBC-derived mitochondrial DNA (mtDNA) and mtRNA, respectively. Interleukin-1ß (IL-1ß) production was initiated by the RLR antiviral signaling adaptor (MAVS) pathway recognition of Mito+ RBC-derived mtRNA. This led to the cytosolic release of Mo mtDNA, which activated the inflammasome. Importantly, mIL-1ß secretion was independent of gasdermin D (GSDMD) and pyroptosis but relied on IFN-inducible myxovirus-resistant protein 1 (MxA), which facilitated the incorporation of mIL-1ß into a trans-Golgi network (TGN)-mediated secretory pathway. RBC internalization identified a subset of blood Mo expressing IFN-stimulated genes (ISGs) that released mIL-1ß and expanded in SLE patients with active disease.
ABSTRACT
SUMMARY: Subcluster analysis is a powerful means to improve clustering and characterization of single cell RNA-Seq data. However, there are no existing tools to systematically integrate results from multiple subclusters, which creates hurdles for accurate data quantification, visualization, and interpretation in downstream analysis. To address this issue, we developed Ragas, an R package that integrates multi-level subclustering objects for streamlined analysis and visualization. A new data structure was implemented to seamlessly connect and assemble miscellaneous single cell analyses from different levels of subclustering, along with several new or enhanced visualization functions. Moreover, a re-projection algorithm was developed to integrate nearest-neighbor graphs from multiple subclusters in order to maximize their separability on the combined cell embeddings, which significantly improved the presentation of rare and homogeneous subpopulations. AVAILABILITY AND IMPLEMENTATION: The Ragas package and its documentation can be accessed through https://github.com/jig4003/Ragas and its source code is also available at https://zenodo.org/records/11244921.
Subject(s)
Algorithms , Single-Cell Analysis , Software , Single-Cell Analysis/methods , Cluster Analysis , Humans , RNA-Seq/methods , Sequence Analysis, RNA/methodsABSTRACT
Femtosecond laser is a promising surface treatment tool for zirconia implant. In this study, the fatigue behavior of zirconia specimens with microgrooved surfaces formed by femtosecond laser is reported. One hundred sixty CAD/CAM zirconia bars (20 mm × 4 mm × 1.4 mm) were evenly divided into four groups with different surface: as sintered; sandblasted with 110 µm Al2O3; femtosecond laser produced microgrooves having 50 µm width, 30 µm depth, and 100 µm pitch; microgrooves having 30 µm width, 20 µm depth, and 60 µm pitch. The femtosecond laser formed micro/nanostructured microgrooves with precise size on zirconia surfaces. XRD analysis indicated that microgrooved surface showed no obvious tetragonal-to-monoclinic phase transformation. The fatigue strength of sandblasted specimens (728 MPa) was significantly higher than that of as sintered specimens (570 MPa). However, the fatigue strength of specimens with microgrooved surface decreased to about 360-380 MPa. The results suggest femtosecond laser is an effective technique to regulate the surface microtopography of zirconia, while further investigations are needed to improve its fatigue behavior.
Subject(s)
Lasers , Zirconium , Surface Properties , Microscopy, Electron, Scanning , Materials Testing , Ceramics , Dental MaterialsABSTRACT
The 17q21 asthma susceptibility locus includes asthma risk alleles associated with decreased sphingolipid synthesis, likely resulting from increased expression of ORMDL3. ORMDL3 inhibits serine-palmitoyl transferase (SPT), the rate-limiting enzyme of de novo sphingolipid synthesis. There is evidence that decreased sphingolipid synthesis is critical to asthma pathogenesis. Children with asthma and 17q21 asthma risk alleles display decreased sphingolipid synthesis in blood cells. Reduced SPT activity results in airway hyperreactivity, a hallmark feature of asthma. 17q21 asthma risk alleles are also linked to childhood infections with human rhinovirus (RV). This study evaluates the interaction of RV with the de novo sphingolipid synthesis pathway, and the alterative effects of concurrent SPT inhibition in SPT-deficient mice and human airway epithelial cells. In mice, RV infection shifted lung sphingolipid synthesis gene expression to a pattern that resembles genetic SPT deficiency, including decreased expression of Sptssa, a small SPT subunit. This pattern was pronounced in lung epithelial cellular adhesion molecule (EpCAM+) cells and reproduced in human bronchial epithelial cells. RV did not affect Sptssa expression in lung CD45+ immune cells. RV increased sphingolipids unique to the de novo synthesis pathway in mouse lung and human airway epithelial cells. Interestingly, these de novo sphingolipid species were reduced in the blood of RV-infected wild-type mice. RV exacerbated SPT deficiency-associated airway hyperreactivity. Airway inflammation was similar in RV-infected wild-type and SPT-deficient mice. This study reveals the effects of RV infection on the de novo sphingolipid synthesis pathway, elucidating a potential mechanistic link between 17q21 asthma risk alleles and rhinoviral infection.
Subject(s)
Membrane Proteins , Rhinovirus , Animals , Child , Humans , Lung/metabolism , Membrane Proteins/metabolism , Mice , Serine C-Palmitoyltransferase/genetics , Serine C-Palmitoyltransferase/metabolism , Sphingolipids/metabolismABSTRACT
The pathogenesis of subsquamous intestinal metaplasia (SSIM), in which glands of Barrett's esophagus (BE) are buried under esophageal squamous epithelium, is unknown. In a rat model of reflux esophagitis, we found that columnar-lined esophagus developed via a wound-healing process involving epithelial-mesenchymal plasticity (EMP) that buried glands under ulcerated squamous epithelium. To explore a role for reflux-induced EMP in BE, we established and characterized human Barrett's organoids and sought evidence of EMP after treatment with acidic bile salts (AB). We optimized media to grow human BE organoids from immortalized human Barrett's cells and from BE biopsies from seven patients, and we characterized histological, morphological, and molecular features of organoid development. Features and markers of EMP were explored following organoid exposure to AB, with and without a collagen I (COL1) matrix to simulate a wound-healing environment. All media successfully initiated organoid growth, but advanced DMEM/F12 (aDMEM) was best at sustaining organoid viability. Using aDMEM, organoids comprising nongoblet and goblet columnar cells that expressed gastric and intestinal cell markers were generated from BE biopsies of all seven patients. After AB treatment, early-stage Barrett's organoids exhibited EMP with loss of membranous E-cadherin and increased protrusive cell migration, events significantly enhanced by COL1. Using human BE biopsies, we have established Barrett's organoids that recapitulate key histological and molecular features of BE to serve as high-fidelity BE models. Our findings suggest that reflux can induce EMP in human BE, potentially enabling Barrett's cells to migrate under adjacent squamous epithelium to form SSIM.NEW & NOTEWORTHY Using Barrett's esophagus (BE) biopsies, we established organoids recapitulating key BE features. During early stages of organoid development, a GERD-like wound environment-induced features of epithelial-mesenchymal plasticity (EMP) in Barrett's progenitor cells, suggesting that reflux-induced EMP can enable Barrett's cells to migrate underneath squamous epithelium to form subsquamous intestinal metaplasia, a condition that may underlie Barrett's cancers that escape detection by endoscopic surveillance, and recurrences of Barrett's metaplasia following endoscopic eradication therapy.
Subject(s)
Barrett Esophagus , Carcinoma, Squamous Cell , Esophageal Neoplasms , Esophagitis, Peptic , Gastroesophageal Reflux , Animals , Barrett Esophagus/pathology , Bile Acids and Salts/pharmacology , Carcinoma, Squamous Cell/complications , Esophageal Neoplasms/complications , Esophageal Neoplasms/pathology , Gastroesophageal Reflux/complications , Humans , Metaplasia , Organoids/pathology , RatsABSTRACT
BACKGROUND & AIMS: In upper airway cells, T helper 2 cytokines that signal through interleukin-4 (IL-4) receptor-α have been shown to stimulate eotaxin-3 secretion via a nongastric proton pump (ngH+,K+ATPase). To seek novel targets for eosinophilic esophagitis (EoE) treatments, we evaluated ngH+,K+ATPase expression in EoE squamous cells, and explored molecular pathways involved in eotaxin-3 secretion by IL-4 receptor-α signaling. METHODS: ngH+,K+ATPase expression in EoE cells was evaluated by quantitative real-time polymerase chain reaction and Western blotting. IL-4-stimulated eotaxin-3 secretion was measured by enzyme-linked immunosorbent assay after treatment with omeprazole, SCH 28080 (potassium-competitive acid blocker), ethylene glycol-bis(ß-aminoethyl)-N,N,N',N'-tetraacetoxymethyl ester (calcium chelator), 2-aminoethoxydiphenyl borate (inhibitor of endoplasmic reticulum calcium release), verapamil, and diltiazem (L-type calcium channel inhibitors). Intracellular calcium transients were measured by Fluo-4 fluorescence. Key experiments were confirmed in EoE primary cells and in RNA sequencing datasets from mucosal biopsies of patients with EoE and controls. RESULTS: EoE cells expressed ngH+,K+ATPase messenger RNA and protein. Omeprazole and SCH 28080 decreased IL-4-stimulated eotaxin-3 secretion. IL-4 increased intracellular calcium transients, and IL-4-stimulated eotaxin-3 secretion was blocked by ethylene glycol-bis(ß-aminoethyl)-N,N,N',N'-tetraacetoxymethyl ester, 2-aminoethoxydiphenyl borate, verapamil, and diltiazem. The combination of omeprazole and verapamil suppressed IL-4-stimulated eotaxin-3 secretion more than either agent alone. EoE biopsies expressed higher ngH+,K+ATPase and exhibited more calcium signaling than controls. CONCLUSIONS: EoE cells express a nongastric proton pump that mediates T helper 2 cytokine-stimulated eotaxin-3 secretion. IL-4 induces calcium release from the endoplasmic reticulum and calcium entry via L-type calcium channels, increasing intracellular calcium that contributes to eotaxin-3 secretion by EoE cells. L-type calcium channel inhibitors block T helper 2 cytokine-stimulated eotaxin-3 secretion, suggesting a potential role for these agents in EoE treatment.
Subject(s)
Chemokine CCL26/metabolism , Eosinophilic Esophagitis/metabolism , Eosinophilic Esophagitis/pathology , Epithelial Cells/metabolism , H(+)-K(+)-Exchanging ATPase/genetics , H(+)-K(+)-Exchanging ATPase/metabolism , Biological Transport/drug effects , Boron Compounds/pharmacology , Calcium/metabolism , Calcium Channel Blockers/pharmacology , Cell Line , Diltiazem/pharmacology , Egtazic Acid/analogs & derivatives , Egtazic Acid/pharmacology , Esophageal Mucosa/metabolism , Esophageal Mucosa/pathology , Famotidine/pharmacology , Female , Histamine H2 Antagonists/pharmacology , Humans , Interleukin-4 Receptor alpha Subunit/metabolism , Male , Omeprazole/pharmacology , Primary Cell Culture , Proton Pump Inhibitors/pharmacology , Proton Pumps/drug effects , Proton Pumps/metabolism , RNA, Messenger/metabolism , Ranitidine/pharmacology , Signal Transduction/drug effects , Th2 Cells/metabolism , Verapamil/pharmacologyABSTRACT
Rationale: A human model to better understand tuberculosis immunopathogenesis and facilitate vaccine development is urgently needed.Objectives: We evaluated the feasibility, safety, and immunogenicity of live bacillus Calmette-Guérin (BCG) in a lung-oriented controlled human infection model.Methods: We recruited 106 healthy South African participants with varying degrees of tuberculosis susceptibility. Live BCG, sterile PPD, and saline were bronchoscopically instilled into separate lung segments (n = 65). A control group (n = 34) underwent a single bronchoscopy without challenge. The primary outcome was safety. Cellular and antibody immune signatures were identified in BAL before and 3 days after challenge using flow cytometry, ELISA, RNA sequencing, and mass spectrometry.Measurements and Main Results: The frequency of adverse events was low (9.4%; n = 10), similar in the challenge versus control groups (P = 0.8), and all adverse events were mild and managed conservatively in an outpatient setting. The optimal PPD and BCG dose was 0.5 TU and 104 cfu, respectively, based on changes in BAL cellular profiles (P = 0.02) and antibody responses (P = 0.01) at incremental doses before versus after challenge. At 104 versus 103 cfu BCG, there was a significant increase in number of differentially expressed genes (367 vs. 3; P < 0.001) and dysregulated proteins (64 vs. 0; P < 0.001). Immune responses were highly setting specific (in vitro vs. in vivo) and compartment specific (BAL vs. blood) and localized to the challenged lung segments.Conclusions: A lung-oriented mycobacterial controlled human infection model using live BCG and PPD is feasible and safe. These data inform the study of tuberculosis immunopathogenesis and strategies for evaluation and development of tuberculosis vaccine candidates.
Subject(s)
Adjuvants, Immunologic/administration & dosage , BCG Vaccine/administration & dosage , Bronchoscopy , Immunogenicity, Vaccine , Tuberculin/administration & dosage , Tuberculosis/prevention & control , Administration, Topical , Adult , Feasibility Studies , Female , Humans , Immunity, Mucosal , Male , Young AdultABSTRACT
MOTIVATION: Transcriptome-based computational drug repurposing has attracted considerable interest by bringing about faster and more cost-effective drug discovery. Nevertheless, key limitations of the current drug connectivity-mapping paradigm have been long overlooked, including the lack of effective means to determine optimal query gene signatures. RESULTS: The novel approach Dr Insight implements a frame-breaking statistical model for the 'hand-shake' between disease and drug data. The genome-wide screening of concordantly expressed genes (CEGs) eliminates the need for subjective selection of query signatures, added to eliciting better proxy for potential disease-specific drug targets. Extensive comparisons on simulated and real cancer datasets have validated the superior performance of Dr Insight over several popular drug-repurposing methods to detect known cancer drugs and drug-target interactions. A proof-of-concept trial using the TCGA breast cancer dataset demonstrates the application of Dr Insight for a comprehensive analysis, from redirection of drug therapies, to a systematic construction of disease-specific drug-target networks. AVAILABILITY AND IMPLEMENTATION: Dr Insight R package is available at https://cran.r-project.org/web/packages/DrInsight/index.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Drug Repositioning , Drug Discovery , Models, Statistical , Software , TranscriptomeABSTRACT
AIMS/HYPOTHESIS: Pancreatic islets produce non-coding microRNAs (miRNAs) that regulate islet cell function and survival. Our earlier investigations revealed that human islets undergo significant damage due to various types of stresses following transplantation and release miRNAs. Here, we sought to identify and validate exosomal miRNAs (exo-miRNAs) produced by human islets under conditions of cellular stress, preceding loss of cell function and death. We also aimed to identify islet stress signalling pathways targeted by exo-miRNAs to elucidate potential regulatory roles in islet cell stress. METHODS: Human islets were subjected to proinflammatory cytokine and hypoxic cell stress and miRNA from exosomes was isolated for RNA sequencing and analysis. Stress-induced exo-miRNAs were evaluated for kinetics of expression and release by intact islets for up to 48 h exposure to cytokines and hypoxia. A subset of stress-induced exo-miRNAs were assessed for recovery and detection as biomarkers of islet cell stress in a diabetic nude mouse xenotransplant model and in patients undergoing total pancreatectomy with islet auto-transplantation (TPIAT). Genes and signalling pathways targeted by stress-induced exo-miRNAs were identified by Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and direct interactions of miRNAs with downstream signalling targets were validated in human islet cells using the miRNA Tests for Read Analysis and Prediction (MirTrap) system. RESULTS: Global exo-miRNA sequencing revealed that 879 miRNA species were released from human islets and 190 islet exo-miRNAs were differentially expressed in response to proinflammatory cytokines, hypoxia or both. Release of exo-miRNAs hsa-miR-29b-3p and hsa-miR-216a-5p was detected within 6 h of exposure to cytokines and hypoxia. The remaining subset of stress-induced exo-miRNAs, including hsa-miR-148a-3p and islet cell damage marker hsa-miR-375, showed delayed release at 24-48 h, correlating with apoptosis and cell death. Stress and damage exo-miRNAs were significantly elevated in the circulation in human-to-mouse xenotransplant models and in human transplant recipients. Elevated blood exo-miRNAs negatively correlated with post-transplant islet function based on comparisons of stress and damage exo-miRNA indices with Secretory Unit of Islet Transplant Objects (SUITO) indices. KEGG analysis and further validation of exo-miRNA targets by MirTrap analysis revealed significant enrichment of islet mRNAs involved in phosphoinositide 3-kinase/Akt and mitogen-activated protein kinase signalling pathways. CONCLUSIONS/INTERPRETATION: The study identifies exo-miRNAs differentially expressed and released by islets in response to damage and stress. These exo-miRNAs could serve as potential biomarkers for assessing islet damage and predicting outcomes in islet transplantation. Notably, exo-miRNAs 29b-3p and 216a-5p could be detected in islets prior to damage-released miRNAs and indicators of cellular apoptosis and death. Thus, these stress-induced exo-miRNAs may have potential diagnostic value for detecting early islet stress prior to progressive loss of islet cell mass and function. Further investigations are warranted to investigate the utility of these exo-miRNAs as early indicators of islet cell stress during prediabetic conditions.
Subject(s)
MicroRNAs/metabolism , Animals , Exosomes/metabolism , Humans , Hypoxia/metabolism , Immunoblotting , Islets of Langerhans/metabolism , Male , Mice , Mice, Nude , MicroRNAs/genetics , Signal Transduction/genetics , Signal Transduction/physiology , Transcriptome/geneticsABSTRACT
Motivation: NGS techniques have been widely applied in genetic and epigenetic studies. Multiple ChIP-seq and RNA-seq profiles can now be jointly used to infer functional regulatory networks (FRNs). However, existing methods suffer from either oversimplified assumption on transcription factor (TF) regulation or slow convergence of sampling for FRN inference from large-scale ChIP-seq and time-course RNA-seq data. Results: We developed an efficient Bayesian integration method (CRNET) for FRN inference using a two-stage Gibbs sampler to estimate iteratively hidden TF activities and the posterior probabilities of binding events. A novel statistic measure that jointly considers regulation strength and regression error enables the sampling process of CRNET to converge quickly, thus making CRNET very efficient for large-scale FRN inference. Experiments on synthetic and benchmark data showed a significantly improved performance of CRNET when compared with existing methods. CRNET was applied to breast cancer data to identify FRNs functional at promoter or enhancer regions in breast cancer MCF-7 cells. Transcription factor MYC is predicted as a key functional factor in both promoter and enhancer FRNs. We experimentally validated the regulation effects of MYC on CRNET-predicted target genes using appropriate RNAi approaches in MCF-7 cells. Availability and implementation: R scripts of CRNET are available at http://www.cbil.ece.vt.edu/software.htm. Contact: xuan@vt.edu. Supplementary information: Supplementary data are available at Bioinformatics online.
Subject(s)
RNA/genetics , Sequence Analysis, RNA/methods , Bayes Theorem , Breast Neoplasms/genetics , Humans , Promoter Regions, Genetic , Transcription Factors/metabolismABSTRACT
The human liver's capacity to rapidly regenerate to a full-sized functional organ after resection has allowed successful outcomes for living donor liver transplantation (LDLT) procedures. However, the ability to detect and track physiological changes occurring during liver regeneration after resection and throughout the restoration process is still lacking. We performed a comprehensive whole-transcriptome RNA sequencing analysis of liver and circulating blood tissue from 12 healthy LDLT donors to define biomarker signatures for monitoring physiological activities during liver regeneration at 14 time points for up to a 1-year procedural follow-up. LDLT donor liver tissue differentially expressed 1238 coding and noncoding genes after resection, and an additional 1260 genes were selectively regulated after LDLT. A total of 15,011 RNA transcript species were identified in the blood in response to liver resection. The transcripts most highly regulated were sequentially expressed within 3 distinct peaks that correlated with sets of functional genes involved in the induction of liver resection-specific innate immune response (peak 1), activation of the complement system (peak 2), and platelet activation and erythropoiesis (peak 3). Each peak corresponded with progressive phases of extracellular matrix degradation, remodeling, and organization during liver restoration. These processes could be tracked by distinct molecular signatures of up-regulated and down-regulated gene profiles in the blood during phases of liver repair and regeneration. In conclusion, the results establish temporal and dynamic transcriptional patterns of gene expression following surgical liver resection that can be detected in the blood and potentially used as biomarker signatures for monitoring phases of liver regeneration.
Subject(s)
Hepatectomy/adverse effects , Liver Regeneration/genetics , Liver/physiology , Living Donors , RNA-Seq , Adult , Biomarkers/blood , Cohort Studies , Female , Gene Expression Regulation/physiology , Humans , Liver/surgery , Liver Transplantation , Male , Middle Aged , Tissue and Organ Procurement , Young AdultABSTRACT
Although the anticancer properties of oligomeric proanthocyanidins (OPCs) from grape seeds have been well recognized, the molecular mechanisms by which they exert anticancer effects are poorly understood. In this study, through comprehensive RNA-sequencing-based gene expression profiling in multiple colorectal cancer cell lines, we for the first time illuminate the genome-wide effects of OPCs from grape seeds in colorectal cancer. Our data revealed that OPCs affect several key cancer-associated genes. In particular, genes involved in cell cycle and DNA replication were most significantly and consistently altered by OPCs across multiple cell lines. Intriguingly, our in vivo experiments showed that OPCs were significantly more potent at decreasing xenograft tumor growth compared with the unfractionated grape seed extract (GSE) that includes the larger polymers of proanthocyanidins. These findings were further confirmed in colorectal cancer patient-derived organoids, wherein OPCs more potently inhibited the formation of organoids compared with GSE. Furthermore, we validated alteration of cell cycle and DNA replication-associated genes in cancer cell lines, mice xenografts as well as patient-derived organoids. Overall, this study provides an unbiased and comprehensive look at the mechanisms by which OPCs exert anticancer properties in colorectal cancer.
Subject(s)
Antineoplastic Agents, Phytogenic/pharmacology , Antioxidants/pharmacology , Colorectal Neoplasms/drug therapy , Grape Seed Extract/pharmacology , Proanthocyanidins/pharmacology , Seeds/chemistry , Vitis/chemistry , Animals , Caco-2 Cells , Cell Cycle/drug effects , Cell Line, Tumor , DNA Replication/drug effects , HCT116 Cells , HT29 Cells , Humans , Male , Mice , Mice, NudeABSTRACT
MOTIVATION: The advent of high-throughput DNA methylation profiling techniques has enabled the possibility of accurate identification of differentially methylated genes for cancer research. The large number of measured loci facilitates whole genome methylation study, yet posing great challenges for differential methylation detection due to the high variability in tumor samples. RESULTS: We have developed a novel probabilistic approach, D: ifferential M: ethylation detection using a hierarchical B: ayesian model exploiting L: ocal D: ependency (DM-BLD), to detect differentially methylated genes based on a Bayesian framework. The DM-BLD approach features a joint model to capture both the local dependency of measured loci and the dependency of methylation change in samples. Specifically, the local dependency is modeled by Leroux conditional autoregressive structure; the dependency of methylation changes is modeled by a discrete Markov random field. A hierarchical Bayesian model is developed to fully take into account the local dependency for differential analysis, in which differential states are embedded as hidden variables. Simulation studies demonstrate that DM-BLD outperforms existing methods for differential methylation detection, particularly when the methylation change is moderate and the variability of methylation in samples is high. DM-BLD has been applied to breast cancer data to identify important methylated genes (such as polycomb target genes and genes involved in transcription factor activity) associated with breast cancer recurrence. AVAILABILITY AND IMPLEMENTATION: A Matlab package of DM-BLD is available at http://www.cbil.ece.vt.edu/software.htm CONTACT: Xuan@vt.eduSupplementary information: Supplementary data are available at Bioinformatics online.
Subject(s)
Breast Neoplasms/genetics , DNA Methylation , Sequence Analysis, DNA/methods , Software , Bayes Theorem , Breast Neoplasms/metabolism , Breast Neoplasms/pathology , DNA, Neoplasm , Female , Genomics/methods , Humans , Neoplasm Recurrence, Local/geneticsABSTRACT
MOTIVATION: Gene set analysis is a powerful tool to study the coordinative change of time-course data. However, most existing methods only model the overall change of a gene set, yet completely overlooked heterogeneous time-dependent changes within sub-sets of genes. RESULTS: We have developed a novel statistical method, Phantom, to investigate gene set heterogeneity. Phantom employs the principle of multi-objective optimization to assess the heterogeneity inside a gene set, which also accounts for the temporal dependency in time-course data. Phantom improves the performance of gene set based methods to detect biological changes across time. AVAILABILITY AND IMPLEMENTATION: Phantom webpage can be accessed at: http://www.baylorhealth.edu/Phantom . R package of Phantom is available at https://cran.r-project.org/web/packages/phantom/index.html . CONTACT: jinghua.gu@bswhealth.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Computational Biology/methods , Gene Expression Regulation , Models, Genetic , Software , Humans , Influenza, Human/geneticsABSTRACT
OBJECTIVE: To evaluate the level of serum carbohydrate antigen-125(CA125) and its related factors in patients with bronchiectasis. METHODS: The clinical data of 504 patients with bronchiectasis in Zhejiang Putuo People's Hospital from June 2009 to June 2014 were collected in the study.The patients were divided into CA125 elevated group and CA125 normal group according to serum CA125 level,and the differences of serum CA125,age,gender, white blood cell(WBC),C-reactive protein(CRP), blood glucose and other test indicators were compared between two groups. RESULTS: There were 276 patients including 117 male and 159 female with elevated serum CA125.Their mean age was(66.3±13.1)years and the mean level of CA125 was(83.70±43.87) U/mL. There were 228 patients including 84 male and 144 female with normal CA125 levels. Their mean age was(67.5±10.5) years and the mean level of CA125 was(20.68±9.67)U/mL.The peripheral blood WBC in patients with CA125 elevated group[(10.08±5.68)×10(9)/L] was significantly higher than that in CA125 normal group[7.73±3.46)×10(9)/L], the difference was statistically significant(P<0.05).The medium of CRP level in patients with CA125 elevated group[22.98(3.18~196.88)mg/L] was significantly higher than that in CA125 normal group[6.34(0.50~97.66)mg/L](P<0.05). Correlation analysis showed that CA125 was positively correlated with WBC and CRP(P<0.05). Stepwise regression analysis showed that CRP was the only independent prognostic factors of CA125. Paired t test showed the presence of CA125 serum in patients with bronchiectasis had a significant difference between before and after anti-infection therapy(P<0.05). CONCLUSION: The serum levels of CA125 rise in patients with bronchiectasis,while it decrease after anti-infection therapy.CRP is an independent associated factor of serum CA125 level.
Subject(s)
Bronchiectasis/blood , CA-125 Antigen/blood , Aged , Blood Glucose/analysis , Bronchiectasis/therapy , C-Reactive Protein/analysis , Female , Humans , Male , Middle AgedABSTRACT
BACKGROUND: Recent advances in RNA sequencing (RNA-Seq) technology have offered unprecedented scope and resolution for transcriptome analysis. However, precise quantification of mRNA abundance and identification of differentially expressed genes are complicated due to biological and technical variations in RNA-Seq data. RESULTS: We systematically study the variation in count data and dissect the sources of variation into between-sample variation and within-sample variation. A novel Bayesian framework is developed for joint estimate of gene level mRNA abundance and differential state, which models the intrinsic variability in RNA-Seq to improve the estimation. Specifically, a Poisson-Lognormal model is incorporated into the Bayesian framework to model within-sample variation; a Gamma-Gamma model is then used to model between-sample variation, which accounts for over-dispersion of read counts among multiple samples. Simulation studies, where sequencing counts are synthesized based on parameters learned from real datasets, have demonstrated the advantage of the proposed method in both quantification of mRNA abundance and identification of differentially expressed genes. Moreover, performance comparison on data from the Sequencing Quality Control (SEQC) Project with ERCC spike-in controls has shown that the proposed method outperforms existing RNA-Seq methods in differential analysis. Application on breast cancer dataset has further illustrated that the proposed Bayesian model can 'blindly' estimate sources of variation caused by sequencing biases. CONCLUSIONS: We have developed a novel Bayesian hierarchical approach to investigate within-sample and between-sample variations in RNA-Seq data. Simulation and real data applications have validated desirable performance of the proposed method. The software package is available at http://www.cbil.ece.vt.edu/software.htm.
Subject(s)
Bayes Theorem , Gene Expression Profiling/methods , RNA, Messenger/genetics , Sequence Analysis, RNA/methods , Computer Simulation , Gene Expression Regulation , Humans , Models, Genetic , SoftwareABSTRACT
OBJECTIVE: To evaluate the quality of life in patients with chronic obstructive pulmonary disease (COPD) by COPD Assessment Test (CAT) and Body Mass Index, Airflow Obstruction, Dyspnea, Exercise Capacity Index (BODE). METHODS: One hundred patients with stable COPD admitted in Putuo People's Hospital were recruited in the study. CAT and BODE index were measured for each patient.The deaths and frequency of exacerbations were recorded during 3-year follow-up period,and the correlation between CAT and BODE in evaluating COPD prognosis was analyzed. RESULTS: There were 28, 30, 29 and 13 patients with CAT score of 1, 2, 3 and 4, respectively; while there were 31, 29, 28 and 12 cases with BODE scores of 1, 2, 3 and 4. CAT scores were well correlated with BODE evaluation in terms of overall score and scores of 4 items (r= -0.237, -0.772, 0.789, -0.767, 0.888, respectively, Ps<0.05). COPD exacerbation incidence and mortality increased with the increasing CAT levels. The rank sum test showed that there were no significant differences between CAT and BODE index in the frequency of acute exacerbation(P<0.05); and in the death toll, the difference was not significant(1 group Χ2=0.919, 2 group Χ2=0.001, 3 group Χ2=0.177,4 group Χ2=0.322, Ps>0.05). CONCLUSION: CAT is relevant to BODE in evaluating incidence of exacerbation and mortality for patients with COPD and CAT is more easily to be applied.
Subject(s)
Pulmonary Disease, Chronic Obstructive/diagnosis , Quality of Life , Aged , Aged, 80 and over , Female , Humans , Male , Middle Aged , PrognosisABSTRACT
In recent years, m6A modifications in RNA transcripts have arisen as a hot topic in cancer research. Indeed, a number of independent studies have elaborated that the m6A modification impacts the behavior of tumor cells and tumor-infiltrating immune cells, altering tumor cell metabolism along with the differentiation and functional activity of immune cells. This review elaborates on the links between RNA m6A modifications, tumor cell metabolism, and immune cell behavior, discussing this topic from the viewpoint of reciprocal regulation through "RNA m6A-tumor cell metabolism-immune cell behavior" and "RNA m6A-immune cell behavior-tumor cell metabolism" axes. In addition, we discuss the various factors affecting RNA m6A modifications in the tumor microenvironment, particularly the effects of hypoxia associated with cancer cell metabolism along with immune cell-secreted cytokines. Our analysis proposes the conclusion that RNA m6A modifications support widespread interactions between tumor metabolism and tumor immunity. With the current viewpoint that long-term cancer control must tackle cancer cell malignant behavior while strengthening anti-tumor immunity, the recognition of RNA m6A modifications as a key factor provides a new direction for the targeted therapy of tumors. This article is categorized under: RNA Processing > RNA Editing and Modification RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.
Subject(s)
Neoplasms , RNA Methylation , Humans , RNA Processing, Post-Transcriptional , Neoplasms/genetics , Biological Transport , RNA , Tumor MicroenvironmentABSTRACT
Genetic variants can alter the profile of heritable molecules such as small RNAs in sperm and oocytes, and in this manner ancestral genetic variants can have a significant effect on offspring phenotypes even if they are not themselves inherited. Here we show that wild type female mice descended from ancestors with a mutation in the mammalian germ cell gene Khdc3 have hepatic metabolic defects that persist over multiple generations. We find that genetically wild type females descended from Khdc3 mutants have transcriptional dysregulation of critical hepatic metabolic genes, which persist over multiple generations and pass through both female and male lineages. This was associated with dysregulation of hepatically-metabolized molecules in the blood of these wild type mice with mutational ancestry. The oocytes of Khdc3-null females, as well as their wild type descendants, had dysregulation of multiple small RNAs, suggesting that these epigenetic changes in the gametes transmit the phenotype between generations. Our results demonstrate that ancestral mutation in Khdc3 can produce transgenerational inherited phenotypes, potentially indefinitely.