Genome-wide Gene Expression Profiling Research Articles

Systems biology comprehensive analysis on breast cancer for identification of key gene modules and genes associated with TNM-based clinical stages

Breast cancer (BC), as one of the leading causes of death among women, comprises several subtypes with controversial and poor prognosis. Considering the TNM (tumor, lymph node, metastasis) based classification for staging of breast cancer, it is essential to diagnose the disease at early stages. The present study aims to take advantage of the systems biology approach on genome wide gene expression profiling datasets to identify the potential biomarkers involved at stage I, stage II, stage III, and stage IV as well as in the integrated group. Three HER2-negative breast cancer microarray datasets were retrieved from the GEO database, including normal, stage I, stage II, stage III, and stage IV samples. Additionally, one dataset was also extracted to test the developed predictive models trained on the three datasets. The analysis of gene expression profiles to identify differentially expressed genes (DEGs) was performed after preprocessing and normalization of data. Then, statistically significant prioritized DEGs were used to construct protein–protein interaction networks for the stages for module analysis and biomarker identification. Furthermore, the prioritized DEGs were used to determine the involved GO enrichment and KEGG signaling pathways at various stages of the breast cancer. The recurrence survival rate analysis of the identified gene biomarkers was conducted based on Kaplan–Meier methodology. Furthermore, the identified genes were validated not only by using several classification models but also through screening the experimental literature reports on the target genes. Fourteen (21 genes), nine (17 genes), eight (10 genes), four (7 genes), and six (8 genes) gene modules (total of 53 unique genes out of 63 genes with involving those with the same connectivity degree) were identified for stage I, stage II, stage III, stage IV, and the integrated group. Moreover, SMC4, FN1, FOS, JUN, and KIF11 and RACGAP1 genes with the highest connectivity degrees were in module 1 for abovementioned stages, respectively. The biological processes, cellular components, and molecular functions were demonstrated for outcomes of GO analysis and KEGG pathway assessment. Additionally, the Kaplan–Meier analysis revealed that 33 genes were found to be significant while considering the recurrence-free survival rate as an alternative to overall survival rate. Furthermore, the machine learning calcification models show good performance on the determined biomarkers. Moreover, the literature reports have confirmed all of the identified gene biomarkers for breast cancer. According to the literature evidence, the identified hub genes are highly correlated with HER2-negative breast cancer. The 53-mRNA signature might be a potential gene set for TNM based stages as well as possible therapeutics with potentially good performance in predicting and managing recurrence-free survival rates at stages I, II, III, and IV as well as in the integrated group. Moreover, the identified genes for the TNM-based stages can also be used as mRNA profile signatures to determine the current stage of the breast cancer.

Integrative analysis of DNA methylation and gene expression profiles identified potential breast cancer-specific diagnostic markers.

Breast cancer is a common malignant tumor among women whose prognosis is largely determined by the period and accuracy of diagnosis. We here propose to identify a robust DNA methylation-based breast cancer-specific diagnostic signature. Genome-wide DNA methylation and gene expression profiles of breast cancer patients along with their adjacent normal tissues from the Cancer Genome Atlas (TCGA) were obtained as the training set. CpGs that with significantly elevated methylation level in breast cancer than not only their adjacent normal tissues and the other ten common cancers from TCGA but also the healthy breast tissues from the Gene Expression Omnibus (GEO) were finally remained for logistic regression analysis. Another independent breast cancer DNA methylation dataset from GEO was used as the testing set. Lots of CpGs were hyper-methylated in breast cancer samples compared with adjacent normal tissues, which tend to be negatively correlated with gene expressions. Eight CpGs located at RIIAD1, ENPP2, ESPN, and ETS1, were finally retained. The diagnostic model was reliable in separating BRCA from normal samples. Besides, chromatin accessibility status of RIIAD1, ENPP2, ESPN and ETS1 showed great differences between MCF-7 and MDA-MB-231 cell lines. In conclusion, the present study should be helpful for breast cancer early and accurate diagnosis.

Glucocorticoid regulates mesenchymal cell differentiation required for perinatal lung morphogenesis and function.

While antenatal glucocorticoids are widely used to enhance lung function in preterm infants, cellular and molecular mechanisms by which glucocorticoid receptor (GR) signaling influences lung maturation remain poorly understood. Deletion of the glucocorticoid receptor gene (Nr3c1) from fetal pulmonary mesenchymal cells phenocopied defects caused by global Nr3c1 deletion, while lung epithelial- or endothelial-specific Nr3c1 deletion did not impair lung function at birth. We integrated genome-wide gene expression profiling, ATAC-seq, and single cell RNA-seq data in mice in which GR was deleted or activated to identify the cellular and molecular mechanisms by which glucocorticoids control prenatal lung maturation. GR enhanced differentiation of a newly defined proliferative mesenchymal progenitor cell (PMP) into matrix fibroblasts (MFBs), in part by directly activating extracellular matrix-associated target genes, including Fn1, Col16a4, and Eln and by modulating VEGF, JAK-STAT, and WNT signaling. Loss of mesenchymal GR signaling blocked fibroblast progenitor differentiation into mature MFBs, which in turn increased proliferation of SOX9+ alveolar epithelial progenitor cells and inhibited differentiation of mature alveolar type II (AT2) and AT1 cells. GR signaling controls genes required for differentiation of a subset of proliferative mesenchymal progenitors into matrix fibroblasts, in turn, regulating signals controlling AT2/AT1 progenitor cell proliferation and differentiation and identifying cells and processes by which glucocorticoid signaling regulates fetal lung maturation.

Harnessing Gene Expression Profiles for the Identification of Ex Vivo Drug Response Genes in Pediatric Acute Myeloid Leukemia.

Novel treatment strategies are of paramount importance to improve clinical outcomes in pediatric AML. Since chemotherapy is likely to remain the cornerstone of curative treatment of AML, insights in the molecular mechanisms that determine its cytotoxic effects could aid further treatment optimization. To assess which genes and pathways are implicated in tumor drug resistance, we correlated ex vivo drug response data to genome-wide gene expression profiles of 73 primary pediatric AML samples obtained at initial diagnosis. Ex vivo response of primary AML blasts towards cytarabine (Ara C), daunorubicin (DNR), etoposide (VP16), and cladribine (2-CdA) was associated with the expression of 101, 345, 206, and 599 genes, respectively (p < 0.001, FDR 0.004–0.416). Microarray based expression of multiple genes was technically validated using qRT-PCR for a selection of genes. Moreover, expression levels of BRE, HIF1A, and CLEC7A were confirmed to be significantly (p < 0.05) associated with ex vivo drug response in an independent set of 48 primary pediatric AML patients. We present unique data that addresses transcriptomic analyses of the mechanisms underlying ex vivo drug response of primary tumor samples. Our data suggest that distinct gene expression profiles are associated with ex vivo drug response, and may confer a priori drug resistance in leukemic cells. The described associations represent a fundament for the development of interventions to overcome drug resistance in AML, and maximize the benefits of current chemotherapy for sensitive patients.

SUN-LB106 The Transcriptomic Evidences on Role of Abdominal Visceral vs. Subcutaneous Adipose Tissue in the Pathophysiology of Diabetes in Asian Indian Indicates the Involvement of Both

Introduction: Asian Indians show “thin fat phenotype” characterized by higher visceral adipose tissue(VAT) and lower subcutaneous adipose tissue(SAT) mass and their higher cardio-metabolic risk has been attributed to this fat distribution. However, the underlying molecular pathology and role of these adipose depots in the pathogenesis of T2D in them remains unknown.Hypothesis: The comparison of transcription profiles of abdominal VAT and SAT and their correlation with diabetes related intermediate phenotypic traits could shed some light on their role in the pathophysiology of diabetes.MethodologySubjects: 19 diabetics (M: F ratio, 8:11) and 16 age and BMI matched controls (M: F ratio 5:11) undergoing abdominal surgery (non-malignant and non-infective conditions).Clinical Parameters: Anthropometry, Serum glucose, insulin, HOMA-R, HbA1c, lipid profile, FFA, adipocytokines. Abdominal VAT, SAT and liver fat were estimated by MRI.Adipose tissue biopsy: SAT and VAT samples were taken during surgery. Genome-wide gene expression profiling of these biopsies was performed using Affymetrix GeneChipPrimeView® arrays. The data was submitted to NCBI-GEO (Accession # GSE78721). Selected genes were validated by qPCR. Gene set enrichment analysis (GSEA) for functional and Weighted Gene Correlation Analysis (WGCNA) for statistical comparison was done.Results:Diabetics had higher waist circumference (p=0.05), HOMA-R (p=0.0002), Visceral fat content (p=0.02) and adipocyte size (p=0.02)GSEA: diabetics vs. controls: In VAT 16 gene sets were upregulated (FDR < 25%) enriching various immune system and inflammation-related pathways. In SAT too, various inflammatory genes were upregulated however they were statistically non-significant (FDR > 25%). Moreover, 12 out of 16 significantly enriched pathways in VAT were among the top 20 pathways in SAT. GSEA in diabetics: VAT vs SAT: None of the gene sets were found significant at FDR < 25% which substantiated our hypothesis that overall pathophysioloigcal alteration in both depots are similar. WGCNA for statistical comparison of VAT and SAT depots The correlation between measures of average gene expression and overall connectivity between both depots was significantly positive. Several modules of co-expressed genes in both VAT and SAT showed positive as well as negative correlation with various intermediate phenotypic traits of diabetes. In both depots they enriched several pathways otherwise known to be associated with pathological adipose tissue like inflammation, adipogenesis etc.ConclusionsIn Asian Indians, diabetes pathology inflicts similar molecular alternations in VAT and SAT, which are more intense in the former. The role of both adipose depots in the pathophysiology of diabetes is along similar lines and they enrich several molecular pathways which are otherwise known to be implicated in pathological adipose tissue.

Examination of methionine stimulation of gene expression in dairy cow mammary epithelial cells using RNA-sequencing.

In this research communication, a cell model with elevated β-CASEIN synthesis was established by stimulating bovine mammary epithelial cells with 0.6 mM methionine, and the genome-wide gene expression profiles of methionine-stimulated cells and untreated cells were investigated by RNA sequencing. A total of 458 differentially expressed genes (DEGs; 219 upregulated and 239 downregulated) were identified between the two groups. Gene Ontology (GO) analysis showed that the two highest-ranked GO terms in 'molecular function' category were 'binding' and 'catalytic activity', suggesting that milk protein synthesis in methionine-stimulated cells requires induction of gene expression to increase metabolic activity. Kyoto Encyclopedia of Genes and Genomes analysis revealed that within the 'environmental information processing' category, the subcategory that is most highly enriched for DEGs was 'signal transduction'. cGMP-PKG, Rap1, calcium, cAMP, PI3K-AKT, MAPK, and JAK-STAT are the pathways with the highest number of DEGs, suggesting that these signaling pathways have potential roles in mediating methionine-induced milk protein synthesis in bovine mammary epithelial cells. This study provides valuable insights into the physiological and metabolic adaptations in cells stimulated with methionine. Understanding the regulation of this transition is essential for effective intervention in the lactation process.

Integrating Genome-Wide Association Studies and Gene Expression Profiles With Chemical-Genes Interaction Networks to Identify Chemicals Associated With Colorectal Cancer.

Colorectal cancer (CRC) is the third most common cancer and has the second highest mortality rate in global cancer. Exploring the associations between chemicals and CRC has great significance in prophylaxis and therapy of tumor diseases. This study aims to explore the relationships between CRC and environmental chemicals on genetic basis by bioinformatics analysis. The genome-wide association study (GWAS) datasets for CRC were obtained from the UK Biobank. The GWAS data for colon cancer (category C18) includes 2,581 individuals and 449,683 controls, while that of rectal cancer (category C20) includes 1,244 individuals and 451,020 controls. In addition, we derived CRC gene expression datasets from the NCBI-GEO (GSE106582). The chemicals related gene sets were acquired from the comparative toxicogenomics database (CTD). Transcriptome-wide association study (TWAS) analysis was applied to CRC GWAS summary data and calculated the expression association testing statistics by FUSION software. We performed chemicals related gene set enrichment analysis (GSEA) by integrating GWAS summary data, mRNA expression profiles of CRC and the CTD chemical-gene interaction networks to identify relationships between chemicals and genes of CRC. We observed several significant correlations between chemicals and CRC. Meanwhile, we also detected 5 common chemicals between colon and rectal cancer, including methylnitronitrosoguanidine, isoniazid, PD 0325901, sulindac sulfide, and importazole. Our study performed TWAS and GSEA analysis, linked prior knowledge to newly generated data and thereby helped identifying chemicals related to tumor genes, which provides new clues for revealing the associations between environmental chemicals and cancer.

Glutathione dynamics determine the therapeutic efficacy of mesenchymal stem cells for graft-versus-host disease via CREB1-NRF2 pathway.

Glutathione (GSH), the most abundant nonprotein thiol functioning as an antioxidant, plays critical roles in maintaining the core functions of mesenchymal stem cells (MSCs), which are used as a cellular immunotherapy for graft-versus-host disease (GVHD). However, the role of GSH dynamics in MSCs remains elusive. Genome-wide gene expression profiling and high-throughput live-cell imaging assays revealed that CREB1 enforced the GSH-recovering capacity (GRC) of MSCs through NRF2 by directly up-regulating NRF2 target genes responsible for GSH synthesis and redox cycling. MSCs with enhanced GSH levels and GRC mediated by CREB1-NRF2 have improved self-renewal, migratory, anti-inflammatory, and T cell suppression capacities. Administration of MSCs overexpressing CREB1-NRF2 target genes alleviated GVHD in a humanized mouse model, resulting in improved survival, decreased weight loss, and reduced histopathologic damages in GVHD target organs. Collectively, these findings demonstrate the molecular and functional importance of the CREB1-NRF2 pathway in maintaining MSC GSH dynamics, determining therapeutic outcomes for GVHD treatment.

Novel Insights into Rheumatoid Arthritis Through Characterisation of Concordant Changes in DNA Methylation and Gene Expression in Synovial Biopsies of Patients with Differing Numbers of Swollen Joints

In this study, we sought to characterize synovial tissue obtained from individuals with arthralgia and disease-specific auto-antibodies, and patients with established rheumatoid arthritis (RA) by applying an integrative multi-omics approach where we investigated changes at the level of DNA methylation and gene expression and its relation to disease pathogenesis. We performed genome-wide DNA methylation and gene expression profiling of synovial tissue obtained from the knee and ankle from 4 auto-antibody positive arthralgia patients and thirteen RA patients. Through multi-omics factor analysis we observed that the latent factor explaining the variance in gene expression and DNA methylation was associated with Swollen Joint Count 66 (SJC66) where patients with SJC66 of 9 or more displayed reasonable separation from the rest. Specifically, combined DNA methylation and gene expression patterns of patients with SJC66 ≥ 9 displayed activation of the immune response as well as dysregulation of cell adhesion pathways. Through an integrative analysis framework, we reveal differences in both methylation and expression at 59 genes associated with the difference between SJC66 ≥ 9 and 0 < SJC66 ≤ 8, for which we identify specific transcripts correlated with differential methylation. Our results highlight the utility of genome-wide multi-omics profiling of synovial samples for improved understanding of the changes associated with an increasing number of swollen joints in RA patients, and point to novel candidate targets for the treatment of the disease. Funding Statement: AYFLY was funded by the European Union’s Horizon 2020 research and innovation program under Grant Agreement No. ITN-2014-EID-641665. Declaration of Interests: AYFLY, KM, GR, EC, HDL, DFT, CL, DG, and RB are employees of GSK. EF and PPT are employees of Novartis and Kintai Therapeutics, respectively. WJJ is financially supported by GSK and Mead Johnson. All declarations are correct for the time of writing and/or conducting this study. Ethics Approval Statement: All subjects provided written informed consent and the collection and use of the samples received Institutional Review Board review and approval. Characteristics of patients included in this study are listed in Table S1.

MiR-30a and miR-379 modulate retinoic acid pathway by targeting DNA methyltransferase 3B in oral cancer

BackgroundEpigenetic silencing of retinoic acid (RA) signaling-related genes have been linked with the pathogenesis and clinical outcome in oral squamous cell carcinoma (OSCC) carcinogenesis. However, the precise mechanisms underlying the abnormal silencing of RA signaling-related genes in OSCC have not been well investigated.MethodsUsing combined analysis of genome-wide gene expression and methylation profile from 40 matched normal-tumor pairs of OSCC specimens, we found a set of retinoid signaling related genes are frequently hypermethylated and downregulated in OSCC patient samples, including alcohol dehydrogenase, iron containing 1 (ADHFE1) and aldehyde dehydrogenase 1 family, member A2 (ALDH1A2), which are the important rate-limiting enzymes in synthesis of RA. The expression of ADHFE1 and ALDH1A2 in OSCC patients was determine by quantitative real-time PCR (qRT-PCR) and immunohistochemistry. The binding sites of miR-30a and miR-379 with DNA methyltransferase 3B (DNMT3B) were predicted using a series of bioinformatic tools, and validated using dual luciferase assay and Western blot analyses. The functions of miR-30a, miR-379, and DNMT3B were accessed by growth and colony formation analyses using gain- and loss-of-function approaches. Chromatin immunoprecipitation (ChIP) was performed to explore the molecular mechanisms by arecoline and 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) treatment.ResultsWe demonstrated that deregulated miR-30a and miR-379 could represent a mechanism for the silencing of ADHFE1 and ALDH1A2 in OSCC through targeting DNMT3B. Ectopic expression of miR-30a and miR-379 could induce re-expression of methylation-silenced ADHFE1 and ALDH1A2, and lead to growth inhibition in oral cancer cells. Furthermore, the dysregulation of the miRNAs and DNMT-3B may result from exposure to tobacco smoking and betel quid chewing.ConclusionsOur results demonstrate that tobacco smoking and betel quid chewing could repress miR-30a and miR-379, which upregulate the DNMT3B expression, in turn, lead to the hypermethylation of ADHFE1 and ALDH1A genes, consequently, promote the oncogenic activity. These findings highlight the potential use of retinoids in combination with epigenetic modifiers for the prevention or treatment of oral cancer.

The ankyrin repeat gene family in Capsicum spp: Genome-wide survey, characterization and gene expression profile

The ankyrin (ANK) repeat protein family is largely distributed across plants and has been found to participate in multiple processes such as plant growth and development, hormone response, response to biotic and abiotic stresses. It is considered as one of the major markers of capsaicin content in pepper fruits. In this study, we performed a genome-wide identification and expression analysis of genes encoding ANK proteins in three Capsicum species: Capsicum baccatum, Capsicum annuum and Capsicum chinense. We identified a total of 87, 85 and 96 ANK genes in C. baccatum, C. annuum and C. chinense genomes, respectively. Next, we performed a comprehensive bioinformatics analysis of the Capsicum ANK gene family including gene chromosomal localization, Cis-elements, conserved motif identification, intron/exon structural patterns and gene ontology classification as well as profile expression. Phylogenetic and domain organization analysis grouped the Capsicum ANK gene family into ten subfamilies distributed across all 12 pepper chromosomes at different densities. Analysis of the expression of ANK genes in leaf and pepper fruits suggested that the ANKs have specific expression patterns at various developmental stages in placenta tissue. Our results provide valuable information for further studies of the evolution, classification and putative functions of ANK genes in pepper.

CIP2A Constrains Th17 Differentiation by Modulating STAT3 Signaling.

SummaryCancerous Inhibitor of Protein Phosphatase 2A (CIP2A) is an oncogene and a potential cancer therapy target protein. Accordingly, a better understanding of the physiological function of CIP2A, especially in the context of immune cells, is a prerequisite for its exploitation in cancer therapy. Here, we report that CIP2A negatively regulates interleukin (IL)-17 production by Th17 cells in human and mouse. Interestingly, concomitant with increased IL-17 production, CIP2A-deficient Th17 cells had increased strength and duration of STAT3 phosphorylation. We analyzed the interactome of phosphorylated STAT3 in CIP2A-deficient and CIP2A-sufficient Th17 cells and indicated together with genome-wide gene expression profiling, a role of Acylglycerol Kinase (AGK) in the regulation of Th17 differentiation by CIP2A. We demonstrated that CIP2A regulates the strength of the interaction between AGK and STAT3, and thereby modulates STAT3 phosphorylation and expression of IL-17 in Th17 cells.

Genome-wide gene expression profiling of the melon fly, Zeugodacus cucurbitae, during thirteen life stages

The melon fly, Zeugodacus cucurbitae (Coquillett), is an important destructive pest worldwide. Functional studies of the genes associated with development and reproduction during different life stages are limited in Z. cucurbitae. There have yet to be comprehensive transcriptomic resources for genetic and functional genomic studies to identify the molecular mechanisms related to its development and reproduction. In this study, we comprehensively sequenced the transcriptomes of four different developmental stages: egg, larva, pupa, and adults. Using the Illumina RNA-Seq technology, we constructed 52 libraries from 13 stages with four biological replicates in each and generated 435.61 Gb clean reads. We comprehensively characterized the transcriptomes with high-coverage mapping to the reference genome. A total of 13,760 genes were mapped to the reference genome, and another 4481 genes were characterized as new genes. Finally, 14,931 genes (81.85%) were functionally annotated against six annotation databases. This study provides the first comprehensive transcriptome data of all developmental stages of Z. cucurbitae, and will serve as a valuable resource for future genetic and functional studies.

Integrated Genome-Wide Analysis of an Isogenic Pair of Pseudomonas aeruginosa Clinical Isolates with Differential Antimicrobial Resistance to Ceftolozane/Tazobactam, Ceftazidime/Avibactam, and Piperacillin/Tazobactam.

Multidrug-resistant (MDR) Pseudomonas aeruginosa is one of the main causes of morbidity and mortality in hospitalized patients and the leading cause of nosocomial infections. We investigated, here, two MDR P. aeruginosa clinical isolates from a hospitalized patient with differential antimicrobial resistance to ceftazidime/avibactam (CZA), ceftolozane/tazobactam (C/T), and piperacillin/tazobactam (P/T). Their assembled complete genomes revealed they belonged to ST235, a widespread MDR clone; and were isogenic with only a single nucleotide variant, causing G183D mutation in AmpC β-lactamase, responsible for a phenotypic change from susceptible to resistant to CZA and C/T. Further epigenomic profiling uncovered two conserved DNA methylation motifs targeted by two distinct putative methyltransferase-containing restriction-modification systems, respectively; more intriguingly, there was a significant difference between the paired isolates in the pattern of genomic DNA methylation and modifications. Moreover, genome-wide gene expression profiling demonstrated the inheritable genomic methylation and modification induced 14 genes being differentially regulated, of which only toxR (downregulated), a regulatory transcription factor, had its promoter region differentially methylate and modified. Since highly expressed opdQ encodes an OprD porin family protein, therefore, we proposed an epigenetic regulation of opdQ expression pertinent to the phenotypic change of P. aeruginosa from resistant to susceptible to P/T. The disclosed epigenetic mechanism controlling phenotypic antimicrobial resistance deserves further experimental investigation.

Jiaji (EX-B2)-Based Electroacupuncture Preconditioning Attenuates Early Ischaemia Reperfusion Injury in the Rat Myocardium.

Background Acupuncture preconditioning was able to reduce the extent of ischaemia reperfusion (I/R) injury. Previous studies have shown that electroacupuncture (EA) pretreatment at T4-T5 Jiaji (EX-B2) acupoints had cardioprotective effects against myocardial I/R injury. However, the molecular mechanism remains inconclusive. Methods Wistar rats were pretreated with electroacupuncture for 7 days at the Neiguan (PC6), T4-T5 Jiaji (EX-B2), Yanglingquan (GB34), and Quchi (LI11) acupoints, which belong to different meridians. Then, we investigated the genome-wide gene expression profiles of rats prestimulated at these acupoints after I/R injury. Results Our study revealed previously unknown cardioprotective roles of T4-T5 Jiaji (EX-B2) acupoints in the I/R progression. The extent of myocardial injury was significantly decreased in the Jiaji group compared with the I/R group. In addition, our data are among the first to link the EA preconditioning at Neiguan (PC6) acupoints and circadian rhythm in the I/R model. Also, for the first time, we explored the meridian and acupoint specificity involved in EA pretreatment at the heart meridian, in which Yanglingquan and Quchi acupoints were selected as the control group for heart-divergent-meridian and nonheart-meridian acupoints. Conclusions The present study suggested that EA pretreatment at Jiaji alters genome-wide gene expression and protects the rat myocardium against I/R injury, which are most likely through neurohumoral regulation.

Genetic Signatures of Acute Asthma Exacerbation Related With Ineffective Response to Corticosteroid.

PurposeAcute exacerbation (AE) is an important domain of asthma management and may be related with ineffective response to corticosteroid. This study aimed to find mechanisms of AE using genome-wide gene expression profiles of blood cells from asthmatics and its perturbation by in vitro dexamethasone (Dex)-treatment.MethodsWe utilized lymphoblastoid B cells from 107 childhood asthmatics and peripheral blood mononuclear cells from 29 adult asthmatics who were treated with inhaled corticosteroids. We searched for a preserved co-expression gene module significantly associated with the AE rate in both cohorts and measured expression changes of genes belong to this module after Dex-treatment.ResultsWe identified a preserved module composed of 77 genes. Among them, expressions of 2 genes (EIF2AK2 and NOL11) decreased significantly after Dex-treatment in both cohorts. EIF2AK2, a key gene acting antiviral defense mechanism, showed significantly higher expressions in asthmatics with AE. The protein repair pathway was enriched significantly in 64 genes which belong to the preserved module but showed no expression differences after Dex-treatment in both cohorts. Among them, MSRA and MSRB2 may play key roles by controlling oxidative stress.ConclusionsMany genes belong to the AE rate-associated and preserved module identified in blood cells from childhood and adults asthmatics showed no expression changes after in vitro Dex-treatment. These findings suggest that we may need alternative treatment options to corticosteroids to prevent AE. EIF2AK2, MSRA and MSRB2 expressions on blood cells may help us select AE-susceptible asthmatics and adjust treatments to prevent AE.

Optimized nitrogen management enhances lodging resistance of rice and its morpho-anatomical, mechanical, and molecular mechanisms

Increasing evidence shows that improved nitrogen management can enhance lodging resistance and lower internodes play a key role in the lodging resistance of rice. However, little is known about the cellular and molecular mechanisms underlying the enhanced lodging resistance under improved nitrogen management. In the present study, two rice varieties, with contrasting lodging resistance, were grown under optimized N management (OPT) and farmers’ fertilizer practices. Under OPT, the lower internodes of both cultivars were shorter but the upper internodes were longer, while both culm diameter and wall thickness of lower internodes were dramatically increased. Microscopic examination showed that the culm wall of lower internodes under OPT contained more sclerenchyma cells beneath epidermis and vascular bundle sheath. The genome-wide gene expression profiling revealed that transcription of genes encoding cell wall loosening factors was down-regulated while transcription of genes participating in lignin and starch synthesis was up-regulated under OPT, resulting in inhibition of longitudinal growth, promotion in transverse growth of lower internodes and enhancement of lodging resistance. This is the first comprehensive report on the morpho-anatomical, mechanical, and molecular mechanisms of lodging resistance of rice under optimized N management.

Transcriptome resilience predicts thermotolerance in Caenorhabditis elegans

BackgroundThe detrimental effects of a short bout of stress can persist and potentially turn lethal, long after the return to normal conditions. Thermotolerance, which is the capacity of an organism to withstand relatively extreme temperatures, is influenced by the response during stress exposure, as well as the recovery process afterwards. While heat-shock response mechanisms have been studied intensively, predicting thermal tolerance remains a challenge.ResultsHere, we use the nematode Caenorhabditis elegans to measure transcriptional resilience to heat stress and predict thermotolerance. Using principal component analysis in combination with genome-wide gene expression profiles collected in three high-resolution time series during control, heat stress, and recovery conditions, we infer a quantitative scale capturing the extent of stress-induced transcriptome dynamics in a single value. This scale provides a basis for evaluating transcriptome resilience, defined here as the ability to depart from stress-expression dynamics during recovery. Independent replication across multiple highly divergent genotypes reveals that the transcriptional resilience parameter measured after a spike in temperature is quantitatively linked to long-term survival after heat stress.ConclusionOur findings imply that thermotolerance is an intrinsic property that pre-determines long-term outcome of stress and can be predicted by the transcriptional resilience parameter. Inferring the transcriptional resilience parameters of higher organisms could aid in evaluating rehabilitation strategies after stresses such as disease and trauma.

H19 promote calcium oxalate nephrocalcinosis-induced renal tubular epithelial cell injury via a ceRNA pathway.

BackgroundIntrarenal calcium oxalate (CaOx) crystals induce inflammation and kidney tubular cell injury, which are processes that involve TLR4/NF-κB signalling. A recent genome-wide gene expression profile analysis of Randall's plaques in CaOx stone patients revealed that the expression of the long noncoding RNA H19 was significantly upregulated. However, to date, its role in kidney CaOx stones has not been reported.MethodA Gene Expression Omnibus (GEO) dataset was utilized to analyse gene expression profiles. Luciferase reporter, Western blotting, qRT-PCR, immunofluorescence staining and reactive oxygen species (ROS) assays were employed to study the molecular mechanism of HMGB1/TLR4/NF-κB regulation by H19 and miR-216b. In vitro and in vivo assays were performed to further confirm the proinflammatory and prooxidative stress effects.FindingH19 expression was significantly increased and positively correlated with the expression levels of HMGB1, TLR4 and NF-κB in Randall's plaques and glyoxylate-induced CaOx nephrocalcinosis mouse models. H19 interacted with miR-216b and suppressed its expression. Additionally, miR-216b inhibited HMGB1 expression by directly binding its 3′-untranslated region. Moreover, H19 downregulation inhibited HMGB1, TLR4 and NF-κB expression and suppressed CaOx nephrocalcinosis-induced renal tubular epithelial cell injury, NADPH oxidase, and oxidative stress in vivo and in vitro. Interestingly, miR-216b inhibition partially reversed the inhibitory effect of H19 knockdown on HMGB1 expression.InterpretationWe determined that H19 might serve as a facilitator in the process of CaOx nephrocalcinosis-induced oxidative stress and renal tubular epithelial cell injury, and we revealed that the interaction between H19 and miR-216b could exert its effect via the HMGB1/TLR4/NF-κB pathway.FundingThis work was supported by the National Nature Science Foundation of China (Nos. 8196030190, 8190033175, 81370805, 81470935, 81900645, 81500534, and 81602236).

Integrative Analysis of Pediatric Acute Leukemia Identifies Immature Subtypes That Span a T Lineage and Myeloid Continuum with Distinct Prognoses

Acute myeloid leukemia (AML) comprises a heterogeneous group of malignancies that are linked by the presence of blasts displaying morphologic and immunophenotypic features of myeloid cell differentiation. With the development of genome-wide gene expression profiling (GEP), array-base comparative genomic hybridization methodologies, and next generation sequencing technologies, the field has gained a greater understanding of the molecular features of this malignancy. Several pathologic lesions have been found to have prognostic implications contributing to a continuous refinement of risk stratification in the context of modern therapy. Recently, the Children's Oncology Group (COG)-National Cancer Institute (NCI) TARGET AML initiative molecularly characterized 993 pediatric AML cases including 197 specimens that underwent comprehensive whole genome sequencing. Of these, 94 carried one of three oncogenic fusions known to be strong drivers of leukemogenesis: RUNX1-RUNX1T1, CBFB-MYH11 and KMT2A rearrangements (KMT2Ar). Among all the alterations detected only ten occurred in more than 5% of subjects, all of which had been previously described. This suggested that low-frequency molecular subsets may exist that require larger cohorts to fully elucidate. To address this limitation, we selected 122 pediatric AML specimens that lacked RUNX1-RUNX1T1, CBFB-MYH11 and KMT2Ar by clinical testing for whole genome (WGS), exome (WES) and RNA (RNAseq) sequencing to enrich for cases that carry low-frequency events. GEP coupled with somatic mutation calls and outcome data were utilized to identify distinct molecular subtypes with prognostic implications. Structural variations, copy number alterations, single nucleotide variations and indels were determined by our established pipelines, as well as an evaluation for regulatory rearrangements driving oncogene overexpression through enhancer hijacking. In addition to known AML somatic mutations and rearrangements in genes such as CEBPA, GATA2, NPM1, WT1, FLT3, NRAS, KRAS, ETV6, Cohesin, NUP98 and KAT6A, we identified rare novel events in known oncogenic drivers. These include a GATA2-ITD as well as the repositioning of a distal MYC enhancer to ectopically activate either the MECOM or BCL11B loci. Interestingly, several AML cases carrying loss of function mutations in polycomb repressive complex 2 (PRC2) genes were found to resemble an early T-cell precursor acute lymphoblastic leukemia (ETP-ALL) GEP by gene set enrichment analysis. ETP-ALL exhibits aberrant expression of stem cell and myeloid markers and has been shown to have a GEP consistent with transformation of a stem cell progenitor. Further, mixed phenotype acute leukemias (MPAL) with T and myeloid lineage characteristics have been previously suggested to be in this spectrum of immature leukemias. We therefore hypothesized that these PRC2-mutated AML cases represented the myeloid end of this continuum. To provide global transcriptional context to these ETP-like AMLs and evaluate a comprehensive cohort encompassing a range of pediatric myeloid malignancies, we integrated results from previously published AML, MPAL, acute megakaryoblastic Leukemia (AMKL), and ETP-ALL datasets that had RNAseq and either WES or WGS available for a total of 436 cases. t-SNE visualization using a 381 gene list derived from the top 100 most variably expressed transcripts within each cohort revealed a clear molecular classifier identifying groups that had consistent mutational compositions and disease outcomes but were agnostic of immunophenotype. This approach allowed the distinction of 63 ETP-like cases comprising a mixture of AML, MPAL, and ETP-ALL leukemias which fell into two subgroups distinguished by FLT3-ITD and PRC2 alterations. Irrespective of treatment approach, FLT3-ITD positive ETP-like leukemias enjoyed a favorable outcome whereas those with PRC2 mutations had a poor prognosis. Our data support a refined classification of pediatric myeloid malignancies based on molecular determinants that can be used for risk stratification in therapeutic trials. Disclosures Gruber: Bristol-Myers Squibb: Consultancy. Rubnitz:AbbVie: Research Funding. Reinhardt:Jazz: Other: Participation in Advisory Boards, Research Funding; CSL Behring: Research Funding; Novartis: Other: Participation in Advisory Boards; Roche: Research Funding. Locatelli:bluebird bio: Consultancy; Bellicum: Consultancy, Membership on an entity's Board of Directors or advisory committees; Amgen: Honoraria, Membership on an entity's Board of Directors or advisory committees; Miltenyi: Honoraria; Novartis: Consultancy, Membership on an entity's Board of Directors or advisory committees. Zwaan:Roche: Consultancy; Servier: Consultancy; Daiichi Sankyo: Consultancy; Novartis: Consultancy; Sanofi: Consultancy; Pfizer: Consultancy, Research Funding; BMS: Research Funding; Celgene: Consultancy, Research Funding; Jazz pharmaceuticals: Other: Travel support; Janssen: Consultancy; Incyte: Consultancy.