Search

The ability to study the bone microenvironment of failed fracture healing may lead to biomarkers for fracture nonunion. Herein the authors describe a technique for isolating individual cells suitable for single-cell RNA sequencing analyses from intramedullary canal tissue collected by reaming during surgery. The purpose was to detail challenges and solutions inherent to the collection and processing of intramedullary canal tissue samples. The authors then examined single-cell RNA sequencing data from fresh and reanimated samples to demonstrate the feasibility of this approach for prospective studies.

Kidney intercalated cells are involved in acid-base homeostasis via vacuolar ATPase expression. Here we report six human intercalated cell subtypes, including hybrid principal-intercalated cells identified from single cell transcriptomics. Phagosome maturation is a biological process that increases in biological pathway analysis rank following exposure to uropathogenic Escherichia coli in two of the intercalated cell subtypes. Real time confocal microscopy visualization of murine renal tubules perfused with green fluorescent protein expressing Escherichia coli or pHrodo Green E. coli BioParticles demonstrates that intercalated cells actively phagocytose bacteria then acidify phagolysosomes. Additionally, intercalated cells have increased vacuolar ATPase expression following in vivo experimental UTI. Taken together, intercalated cells exhibit a transcriptional response conducive to the kidney’s defense, engulf bacteria and acidify the internalized bacteria. Intercalated cells represent an epithelial cell with characteristics of professional phagocytes like macrophages., Kidney intercalated cells are involved in acid-base homeostasis in the kidneys. Here, the authors use single cell transcriptomics and find that interalated cells exhibit a transcriptional response conducive to defense and can engulf and acidify internalized bacteria, similarly to professional phagocytes.

Diffuse large B cell lymphoma (DLBCL) comprises distinct entities due to its heterogeneity. The currently used international prognostic index for DLBCL prognosis prediction is only based on clinical factors and cannot reflect the molecular mechanisms underlying its progression. Here, we aimed to establish a long noncoding RNAs (lncRNA)-based signature for DLBCL prognosis prediction. The data were retrieved from the Gene Expression Omnibus and The Cancer Genome Atlas database. After identifying the differentially expressed lncRNAs (DELs), univariate COX regression, LASSO regression, and stepwise regression analysis were performed to construct a 6-lncRNA risk score system. Kaplan-Meier survival presented that the high-risk group had a significantly poorer overall survival. Based on the risk score and clinical characters, a nomogram was established, which had better predictive accuracy than each factor alone. Finally, weighted gene co-expression network analysis showed that these lncRNAs might regulate immune response, metabolism process, and signal transduction to influence the outcome. Conclusively, our model and nomogram could be reliable prognostic tools for DLBCL patients.

The RNA-binding protein QKI belongs to the hnRNP K-homology domain protein family, a well-known regulator of pre-mRNA alternative splicing and is associated with several neurodevelopmental disorders. Qki is found highly expressed in developing and adult hearts. By employing the human embryonic stem cell (hESC) to cardiomyocyte differentiation system and generating QKI-deficient hESCs (hESCs-QKIdel) using CRISPR/Cas9 gene editing technology, we analyze the physiological role of QKI in cardiomyocyte differentiation, maturation, and contractile function. hESCs-QKIdel largely maintain normal pluripotency and normal differentiation potential for the generation of early cardiogenic progenitors, but they fail to transition into functional cardiomyocytes. In this work, by using a series of transcriptomic, cell and biochemical analyses, and the Qki-deficient mouse model, we demonstrate that QKI is indispensable to cardiac sarcomerogenesis and cardiac function through its regulation of alternative splicing in genes involved in Z-disc formation and contractile physiology, suggesting that QKI is associated with the pathogenesis of certain forms of cardiomyopathies., RNA binding protein Quaking (QKI) is known for its broad function in pre-mRNA splicing and modification and its association with several neurodevelopmental disorders. Here the authors reveal that QKI-mediated regulation of RNA splicing is indispensable to cardiac development and contractile physiology.

Radiologic techniques remain the main method for early detection for breast cancer and are critical to achieve a favorable outcome from cancer. However, more sensitive detection methods to complement radiologic techniques are needed to enhance early detection and treatment strategies. Using our recently established culturing method that allows propagation of normal and cancerous breast epithelial cells of luminal origin, flow cytometry characterization, and genomic sequencing, we show that cancer cells can be detected in breast milk. Cells derived from milk from the breast with cancer were enriched for CD49f+/EpCAM−, CD44+/CD24−, and CD271+ cancer stem–like cells (CSC). These CSCs carried mutations within the cytoplasmic retention domain of HDAC6, stop/gain insertion in MORF4L1, and deletion mutations within SWI/SNF complex component SMARCC2. CSCs were sensitive to HDAC6 inhibitors, BET bromodomain inhibitors, and EZH2 inhibitors, as mutations in SWI/SNF complex components are known to increase sensitivity to these drugs. Among cells derived from breast milk of additional ten women not known to have breast cancer, two of them contained cells that were enriched for the CSC phenotype and carried mutations in NF1 or KMT2D, which are frequently mutated in breast cancer. Breast milk–derived cells with NF1 mutations also carried copy-number variations in CDKN2C, PTEN, and REL genes. The approach described here may enable rapid cancer cell characterization including driver mutation detection and therapeutic screening for pregnancy/postpartum breast cancers. Furthermore, this method can be developed as a surveillance or early detection tool for women at high risk for developing breast cancer. Significance: These findings describe how a simple method for characterization of cancer cells in pregnancy and postpartum breast cancer can be exploited as a surveillance tool for women at risk of developing breast cancer.

The skin is a multilayered organ, equipped with appendages (that is, follicles and glands), that is critical for regulating body temperature and the retention of bodily fluids, guarding against external stresses and mediating the sensation of touch and pain1,2. Reconstructing appendage-bearing skin in cultures and in bioengineered grafts is a biomedical challenge that has yet to be met3–9. Here we report an organoid culture system that generates complex skin from human pluripotent stem cells. We use stepwise modulation of the transforming growth factor β (TGFβ) and fibroblast growth factor (FGF) signalling pathways to co-induce cranial epithelial cells and neural crest cells within a spherical cell aggregate. During an incubation period of 4–5 months, we observe the emergence of a cyst-like skin organoid composed of stratified epidermis, fat-rich dermis and pigmented hair follicles that are equipped with sebaceous glands. A network of sensory neurons and Schwann cells form nerve-like bundles that target Merkel cells in organoid hair follicles, mimicking the neural circuitry associated with human touch. Single-cell RNA sequencing and direct comparison to fetal specimens suggest that the skin organoids are equivalent to the facial skin of human fetuses in the second trimester of development. Moreover, we show that skin organoids form planar hair-bearing skin when grafted onto nude mice. Together, our results demonstrate that nearly complete skin can self-assemble in vitro and be used to reconstitute skin in vivo. We anticipate that our skin organoids will provide a foundation for future studies of human skin development, disease modelling and reconstructive surgery. Skin organoids generated in vitro from human pluripotent stem cells form complex, multilayered skin tissue with hair follicles, sebaceous glands and neural circuitry, and integrate with endogenous skin when grafted onto immunocompromised mice.

The neutrophil-lymphocyte ratio (NLR) and the platelet-lymphocyte ratio (PLR) are markers of systemic inflammation. However, there is little evidence of the value of inflammation in the early diagnosis of gastric cancer (GC). A total of 2,606 patients diagnosed with GC in the past three years and 3,219 healthy controls over the same period were included in this study. Peripheral blood samples were obtained to analyze the NLR, PLR, carcinoembryonic antigen (CEA), and carbohydrate antigen 19-9 (CA19-9). The optimal cutoff levels for the NLR and PLR were defined by receiver operating characteristic (ROC) curve analysis (NLR=2.258, PLR=147.368). The value of different biomarkers for diagnosing GC was compared by the area under the curve (AUC). The NLR and PLR showed diagnostic sensitivity in GC (AUC=0.715, AUC=0.707). Using the Bonferroni correction, the NLR and PLR were superior to CEA and CA19-9 in the diagnosis of GC (P<0.0001). The systemic inflammatory markers were significantly higher in the early stage of GC than tumor markers. After grouping patients and healthy controls by gender, we found that the diagnostic significance of combined NLR and PLR for GC was greater in male patients than in female patients (P<0.0001). The diagnostic value of the NLR and PLR in GC is higher than that of the traditional tumor markers CEA and CA19-9. Systemic markers of inflammation are more valuable in male than female patients.

Phytophthora sojae is an oomycete that causes stem and root rot disease in soybean. P. sojae delivers many RxLR effector proteins, including Avr1b, into host cells to promote infection. We show here that Avr1b interacts with the soybean U-box protein, GmPUB1-1, in yeast two-hybrid, pull down, and bimolecular fluorescence complementation (BIFC) assays. GmPUB1-1, and a homeologous copy GmPUB1-2, are induced by infection and encode 403 amino acid proteins with U-Box domains at their N-termini. Non-synonymous mutations in the Avr1b C-terminus that abolish suppression of cell death also abolished the interaction of Avr1b with GmPUB1-1, while deletion of the GmPUB1-1 C-terminus, but not the U box, abolished the interaction. BIFC experiments suggested that the GmPUB1-1-Avr1b complex is targeted to the nucleus. In vitro ubiquitination assays demonstrated that GmPUB1-1 possesses E3 ligase activity. Silencing of the GmPUB1 genes in soybean cotyledons resulted in loss of recognition of Avr1b by gene products encoded by Rps1-b and Rps1-k. The recognition of Avr1k (which did not interact with GmPUB1-1) by Rps1-k plants was not, however, affected following GmPUB1-1 silencing. Furthermore, over-expression of GmPUB1-1 in particle bombardment experiments triggered cell death suggesting that GmPUB1 may be a positive regulator of effector-triggered immunity. In a yeast two-hybrid system, GmPUB1-1 also interacted with a number of other RxLR effectors including Avr1d, while Avr1b and Avr1d interacted with a number of other infection-induced GmPUB proteins, suggesting that the pathogen uses a multiplex of interactions of RxLR effectors with GmPUB proteins to modulate host immunity.

Myasthenia gravis (MG) is an autoimmune disease associated with autoantibody production that leads to skeletal muscle weakness. The molecular mechanisms underlying MG are not fully understood. We analyzed the gene expression profile (GSE85452) and methylation profile (GSE85647) of MG samples from the GEO database to identify aberrantly methylated-differentially expressed genes. By integrating the datasets, we identified 143 hypermethylation-low expression genes and 91 hypomethylation-high expression genes. Then we constructed PPI network and ceRNA networks by these genes. Phosphatase and tensin homolog (PTEN) and Abelson tyrosine-protein kinase (ABL)1 were critical genes in both PPI networks and ceRNA networks. And potential MG associated lncRNAs were selected by comprehensive analysis of the critical genes and ceRNA networks. In the hypermethylation-low expression genes associated ceRNA network, sirtuin (SIRT)1 was the most important gene and the lncRNA HLA complex (HC) P5 had the highest connection degree. Meanwhile, PTEN was the most important gene and the lncRNA LINC00173 had the highest connection degree in the hypomethylation-high expression genes associated ceRNA network. LINC00173 was validated to be upregulated in MG patients by qRT-PCR (P = 0.005), which indicated LINC00173 might be a potential biomarker for MG. These results provide a basis for future studies on the molecular pathogenesis of MG.

Background & Aims: Liver fibrosis is a pathological healing process resulting from hepatic stellate cell (HSC) activation and the generation of myofibroblasts from activated HSCs. The precise underlying mechanisms of liver fibrogenesis are still largely vague due to lack of understanding the functional heterogeneity of activated HSCs during liver injury. Approach and Results: In this study, to define the mechanism of HSC activation, we performed the transcriptomic analysis at single-cell resolution (scRNA-seq) on HSCs in mice treated with carbon tetrachloride (CCl4). By employing LRAT-Cre:Rosa26mT/mG mice, we were able to isolate an activated GFP-positive HSC lineage derived cell population by fluorescence-activated cell sorter (FACS). A total of 8 HSC subpopulations were identified based on an unsupervised analysis. Each HSC cluster displayed a unique transcriptomic profile, despite all clusters expressing common mouse HSC marker genes. We demonstrated that one of the HSC subpopulations expressed high levels of mitosis regulatory genes, velocity, and monocle analysis indicated that these HSCs are at transitioning and proliferating phases at the beginning of HSCs activation and will eventually give rise to several other HSC subtypes. We also demonstrated cell clusters representing HSC-derived mature myofibroblast populations that express myofibroblasts hallmark genes with unique contractile properties. Most importantly, we found a novel HSC cluster that is likely to be critical in liver regeneration, immune reaction, and vascular remodeling, in which the unique profiles of genes such as Rgs5, Angptl6, and Meg3 are highly expressed. Lastly, we demonstrated that the heterogeneity of HSCs in the injured mouse livers is closely similar to that of cirrhotic human livers. Conclusions: Collectively, our scRNA-seq data provided insight into the landscape of activated HSC populations and the dynamic transitional pathway from HSC to myofibroblasts in response to liver injury.

Genome-wide association studies (GWAS) of complex traits, such as alcohol use disorders (AUD), usually identify variants in non-coding regions and cannot by themselves distinguish whether the associated variants are functional or in linkage disequilibrium with the functional variants. Transcriptome studies can identify genes whose expression differs between alcoholics and controls. To test which variants associated with AUD may cause expression differences, we integrated data from deep RNA-seq and GWAS of four postmortem brain regions from 30 subjects with AUD and 30 controls to analyze allele-specific expression (ASE). We identified 88 genes with differential ASE in subjects with AUD compared to controls. Next, to test one potential mechanism contributing to the differential ASE, we analyzed single nucleotide polymorphisms (SNPs) in the 3′ untranslated regions (3′UTR) of these genes. Of the 88 genes with differential ASE, 61 genes contained 437 SNPs in the 3′UTR with at least one heterozygote among the subjects studied. Using a modified PASSPORT-seq (parallel assessment of polymorphisms in miRNA target-sites by sequencing) assay, we identified 25 SNPs that affected RNA levels in a consistent manner in two neuroblastoma cell lines, SH-SY5Y and SK-N-BE(2). Many of these SNPs are in binding sites of miRNAs and RNA-binding proteins, indicating that these SNPs are likely causal variants of AUD-associated differential ASE. In sum, we demonstrate that a combination of computational and experimental approaches provides a powerful strategy to uncover functionally relevant variants associated with the risk for AUD.

AimsTo examine the role of the basic Helix-loop-Helix (bHLH) transcription factor HAND1 in embryonic and adult myocardium.Methods and resultsHand1 is expressed within the cardiomyocytes of the left ventricle (LV) and myocardial cuff between embryonic days (E) 9.5–13.5. Hand gene dosage plays an important role in ventricular morphology and the contribution of Hand1 to congenital heart defects requires further interrogation. Conditional ablation of Hand1 was carried out using either Nkx2.5 knockin Cre (Nkx2.5Cre) or α-myosin heavy chain Cre (αMhc-Cre) driver. Interrogation of transcriptome data via ingenuity pathway analysis reveals several gene regulatory pathways disrupted including translation and cardiac hypertrophy-related pathways. Embryo and adult hearts were subjected to histological, functional, and molecular analyses. Myocardial deletion of Hand1 results in morphological defects that include cardiac conduction system defects, survivable interventricular septal defects, and abnormal LV papillary muscles (PMs). Resulting Hand1 conditional mutants are born at Mendelian frequencies; but the morphological alterations acquired during cardiac development result in, the mice developing diastolic heart failure.ConclusionCollectively, these data reveal that HAND1 contributes to the morphogenic patterning and maturation of cardiomyocytes during embryogenesis and although survivable, indicates a role for Hand1 within the developing conduction system and PM development.

Bone marrow-derived mesenchymal stem cells (MSCs) improve cardiac function after ischemia/reperfusion injury, in part, due to the release of cytoprotective paracrine factors. Toll-like receptor 4 (TLR4) is expressed in MSCs and regulates the expression of cytoprotective factors, cytokines, and chemokines. Lipopolysaccharide (LPS) stimulation of TLR4 activates two distinct signaling pathways that are either MyD88 dependent or MyD88 independent/TIR-domain-containing adapter-inducing interferon-β (TRIF) dependent. While it was reported previously that LPS treatment improved MSC-mediated cardioprotection, the mechanism underlying such improved effect remains unknown. To study the role of MyD88 signaling in MSC cardioprotective activity, wild type (WT) and MyD88(−/−) MSCs were treated with LPS (200 ng/mL) for 24 h. WT and MyD88(−/−) MSCs with or without LPS pretreatment were infused into the coronary circulation of isolated mouse hearts (Langendorff model) and then subjected to ischemia (25 min) and reperfusion (50 min). Saline served as a negative control. Both untreated and LPS-pretreated WT MSCs significantly improved postischemic recovery of myocardial function of isolated mouse hearts, as evidenced by improved left ventricular developed pressure and ventricular contractility assessment (ie, the rate of left ventricle pressure change over time, ± dp/dt). LPS-pretreated WT MSCs conferred better cardiac function recovery than untreated MSCs; however, such effect of LPS was abolished when using MyD88(−/−) MSCs. In addition, LPS stimulated stat3 activity in WT MSCs, but not MyD88(−/−) MSCs. stat3 small interfering RNA abolished the effect of LPS in improving the cardioprotection of WT MSCs. In conclusion, this study demonstrates that LPS improves MSC-mediated cardioprotection by MyD88-dependent activation of stat3.

Summary It is generally considered that mouse embryonic stem cell (ESC) differentiation into blood cells in vitro recapitulates yolk sac (YS) hematopoiesis. As such, similar to YS-derived B-progenitors, we demonstrate here that ESC-derived B-progenitors differentiate into B-1 and marginal zone B cells, but not B-2 cells in immunodeficient mice after transplantation. ESC-derived B-1 cells were maintained in the recipients for more than 6 months, secreting natural IgM antibodies in vivo. Gene expression profiling displayed a close relationship between ESC- and YS-derived B-1 progenitors. Because there are no hematopoietic stem cells (HSCs) detectable in our ESC differentiation culture, successful long-term engraftment of ESC-derived functional B-1 cells supports the presence of HSC-independent B-1 cell development., Graphical Abstract, Highlights • ESC-derived B-progenitors mature into B-1 cells and MZ B cells in vivo • ESC-derived B-1 cells engrafted in vivo long-term and secrete natural antibodies • ESC-derived B-progenitors are molecularly similar to YS-derived B-progenitors • Long-term B-1 cell engraftment represents HSC-independent lymphopoiesis, Yoshimoto and colleagues demonstrate that mouse ESC-derived B-progenitors engraft in the immunodeficient mice long-term without gene manipulation as IgM-secreting B-1 and MZ B cells, but not conventional B-2 cells. These ESC-derived B-progenitors are molecularly similar to YS-derived B-progenitors. Since no HSCs are produced in their ESC differentiation system, this result supports the emerging paradigm-shifting concept of HSC-independent lymphopoiesis.

Single-cell RNA sequencing reveals gene expression differences between individual cells and also identifies different cell populations that are present in the bulk starting material. To obtain an accurate assessment of patient samples, single-cell suspensions need to be generated as soon as possible once the tissue or sample has been collected. However, this requirement poses logistical challenges for experimental designs involving multiple samples from the same subject since these samples would ideally be processed at the same time to minimize technical variation in data analysis. Although cryopreservation has been shown to largely preserve the transcriptome, it is unclear whether the freeze-thaw process might alter gene expression profiles in a cell-type specific manner or whether changes in cell-type proportions might also occur. To address these questions in the context of multiple myeloma clinical studies, we performed single-cell RNA sequencing (scRNA-seq) to compare fresh and frozen cells isolated from bone marrow aspirates of six multiple myeloma patients, analyzing both myeloma cells (CD138+) and cells constituting the microenvironment (CD138−). We found that cryopreservation using 90% fetal calf serum and 10% dimethyl sulfoxide resulted in highly consistent gene expression profiles when comparing fresh and frozen samples from the same patient for both CD138+ myeloma cells (R ≥ 0.96) and for CD138– cells (R ≥ 0.9). We also demonstrate that CD138– cell-type proportions showed minimal alterations, which were mainly related to small differences in immune cell subtype sensitivity to the freeze-thaw procedures. Therefore, when processing fresh multiple myeloma samples is not feasible, cryopreservation is a useful option in single-cell profiling studies.

Improved understanding of local breast biology that favors the development of estrogen receptor negative (ER-) breast cancer (BC) would foster better prevention strategies. We have previously shown that overexpression of specific lipid metabolism genes is associated with the development of ER- BC. We now report results of exposure of MCF-10A cells and mammary organoids to representative medium- and long-chain polyunsaturated fatty acids. This exposure caused a dynamic and profound change in gene expression, accompanied by changes in chromatin packing density, chromatin accessibility and histone posttranslational modifications (PTMs). We identified 38 metabolic reactions that showed significantly increased activity, including reactions related to one-carbon metabolism. Among these reactions are those that produce S-adenosyl-L-methionine for histone PTMs. Utilizing both an in-vitro model and samples from women at high risk for ER- BC, we show that lipid exposure engenders gene expression, signaling pathway activation, and histone marks associated with the development of ER- BC.

T helper cell differentiation requires lineage-defining transcription factors and factors that have shared expression among multiple subsets. BATF is required for development of multiple Th subsets but functions in a lineage-specific manner. BATF is required for IL-9 production in Th9 cells but in contrast to its function as a pioneer factor in Th17 cells, BATF is neither sufficient nor required for accessibility at the Il9 locus. Here we show that STAT5 is the earliest factor binding and remodeling the Il9 locus to allow BATF binding in both mouse and human Th9 cultures. The ability of STAT5 to mediate accessibility for BATF is observed in other Th lineages and allows acquisition of the IL-9-secreting phenotype. STAT5 and BATF convert Th17 cells into cells that mediate IL-9-dependent effects in allergic airway inflammation and anti-tumor immunity. Thus, BATF requires the STAT5 signal to mediate plasticity at the Il9 locus., BATF is a transcription factor that is needed for IL-9 production by T helper 9 cells. Here the authors show that STAT5 is needed at the Il9 locus to enable BATF to function in this manner and that this interaction can reprogram other T helper subsets into IL-9 producing cells, thus regulating the immune response to disease.

Background: Early human heart and brain development simultaneously occur during embryogenesis. Notably, in human newborns, congenital heart defects strongly associate with neurodevelopmental abnormalities, suggesting a common gene/complex underlying both cardiogenesis and neurogenesis. However, due to lack of in vivo studies, the molecular mechanisms that govern both early human heart and brain development remain elusive. The evolutionarily conserved ATP-dependent SWI/SNF complex is one of the largest chromatin remodeling complexes, consisting of ~15 subunits, including SMARCA2 (also known as BRM) or SMARCA4 (also known as BRG1) as the ATPase catalytic subunit. Several BRG1-associated factors (BAFs), such as ARID1A (Baf250a), have DNA binding capacity and assemble with either BRM or BRG1 to form a functional chromatin-remodeling complex. A single amino acid mutation (Arid1aV 1068G/V1068G ), impaired Arid1a-DNA interactions and resulted in both cardiac neural defects. Mutations in 4 different SWI/SNF subunits including ARID1A/B were identified in human congenital syndromes that include both neural and cardiac defects. Results: Here, we report ARID1A, which is a DNA-binding-subunit of the SWI/SNF epigenetic complex, controls both neurogenesis and cardiogenesis from human embryonic stem cells (hESCs) via employing distinct mechanisms. CRISPR/Cas-9 knockout of ARID1A (ARID1A -/- ) led to spontaneous differentiation of neural cells together with globally enhanced expression of neurogenic genes in undifferentiated hESCs. Additionally, when compared with WT hESCs, cardiac differentiation from ARID1A -/- hESCs was prominently suppressed, whereas neural differentiation was significantly promoted. Whole genome-wide ChIP-seq and ATAC-seq analyses revealed that ARID1A was required to open chromatin accessibility on promoters of essential cardiogenic genes, and temporally associated with key cardiogenic transcriptional factors T and MEF2C during early cardiac development. However, during neural development, transcription of most essential neurogenic genes was dependent on ARID1A and ARID1A could interact with REST, which is a known transcriptional repressor. Conclusions: We uncovered the key and opposite roles by ARID1A to govern both early human cardiac and neural development and characterized the mechanisms. We found global chromatin accessibility on cardiogenic genes was dependent on ARID1A, whereas transcriptional activity of neurogenic genes was regulated by ARID1A, possibly through ARID1A-REST interaction.

Background Early human heart and brain development simultaneously occur during embryogenesis. Notably, in human newborns, congenital heart defects strongly associate with neurodevelopmental abnormalities, suggesting a common gene or complex underlying both cardiogenesis and neurogenesis. However, due to lack of in vivo studies, the molecular mechanisms that govern both early human heart and brain development remain elusive. Results Here, we report ARID1A, a DNA-binding subunit of the SWI/SNF epigenetic complex, controls both neurogenesis and cardiogenesis from human embryonic stem cells (hESCs) through distinct mechanisms. Knockout-of-ARID1A (ARID1A−/−) leads to spontaneous differentiation of neural cells together with globally enhanced expression of neurogenic genes in undifferentiated hESCs. Additionally, when compared with WT hESCs, cardiac differentiation from ARID1A −/− hESCs is prominently suppressed, whereas neural differentiation is significantly promoted. Whole genome-wide scRNA-seq, ATAC-seq, and ChIP-seq analyses reveal that ARID1A is required to open chromatin accessibility on promoters of essential cardiogenic genes, and temporally associated with key cardiogenic transcriptional factors T and MEF2C during early cardiac development. However, during early neural development, transcription of most essential neurogenic genes is dependent on ARID1A, which can interact with a known neural restrictive silencer factor REST/NRSF. Conclusions We uncover the opposite roles by ARID1A to govern both early cardiac and neural development from pluripotent stem cells. Global chromatin accessibility on cardiogenic genes is dependent on ARID1A, whereas transcriptional activity of neurogenic genes is under control by ARID1A, possibly through ARID1A-REST/NRSF interaction.

Expression quantitative trait loci (eQTL) analysis is useful for identifying genetic variants correlated with gene expression, however, it cannot distinguish between causal and nearby non-functional variants. Because the majority of disease-associated SNPs are located in regulatory regions, they can impact allele-specific binding (ASB) of transcription factors and result in differential expression of the target gene alleles. In this study, our aim was to identify functional single-nucleotide polymorphisms (SNPs) that alter transcriptional regulation and thus, potentially impact cellular function. Here, we present regSNPs-ASB, a generalized linear model-based approach to identify regulatory SNPs that are located in transcription factor binding sites. The input for this model includes ATAC-seq (assay for transposase-accessible chromatin with high-throughput sequencing) raw read counts from heterozygous loci, where differential transposase-cleavage patterns between two alleles indicate preferential transcription factor binding to one of the alleles. Using regSNPs-ASB, we identified 53 regulatory SNPs in human MCF-7 breast cancer cells and 125 regulatory SNPs in human mesenchymal stem cells (MSC). By integrating the regSNPs-ASB output with RNA-seq experimental data and publicly available chromatin interaction data from MCF-7 cells, we found that these 53 regulatory SNPs were associated with 74 potential target genes and that 32 (43%) of these genes showed significant allele-specific expression. By comparing all of the MCF-7 and MSC regulatory SNPs to the eQTLs in the Genome-Tissue Expression (GTEx) Project database, we found that 30% (16/53) of the regulatory SNPs in MCF-7 and 43% (52/122) of the regulatory SNPs in MSC were also in eQTL regions. The enrichment of regulatory SNPs in eQTLs indicated that many of them are likely responsible for allelic differences in gene expression (chi-square test, p-value < 0.01). In summary, we conclude that regSNPs-ASB is a useful tool for identifying causal variants from ATAC-seq data. This new computational tool will enable efficient prioritization of genetic variants identified as eQTL for further studies to validate their causal regulatory function. Ultimately, identifying causal genetic variants will further our understanding of the underlying molecular mechanisms of disease and the eventual development of potential therapeutic targets.

SummarySingle cell RNA sequencing is an evolving field to elucidate cellular architecture of adult organs. Using normal breast tissues from healthy volunteers and a rapid procurement/processing/sequencing protocol, 13 breast epithelial cell clusters were identified. Approximately 90% of breast cancers were enriched for cell-of-origin signatures derived from differentiated luminal clusters and two minor luminal progenitor clusters. Expression of cell cycle and chromosome segregation-related genes were higher in one of the minor clusters and breast tumors with this cluster signature displayed the highest mutation rate and poor outcome. We identified TBX3 and PDK4 as genes co-expressed with estrogen receptor (ER) in the normal breasts and their expression analyses in >550 breast cancers enabled prognostically relevant cell-of-origin based subclassification of ER+ breast cancers.SignificanceThis study elucidates different epithelial cell types of the normal breasts and identifies a minor subpopulation of cells from which the majority of breast cancers may originate. This observation should help to develop methods to characterize breast tumors based on cell-of-origin. Although it was suggested that intrinsic subtypes of breast cancers have distinct cells of origins, this study suggests multiple cell-of-origin for an intrinsic subtype of breast cancer, including for hormone responsive breast cancers. Cell-of-origin signatures allowed survival-associated subclassification of intrinsic subtypes. Critically, this normal breast cell atlas would allow for the classification of genes differentially expressed in a breast tumor compared to normal breast due to the cell-of-origin of tumor and those that are acquired due to genomic aberrations.

Background While many studies have described the impact of prenatal opioid exposure on development, possible mechanisms for how opioids exert developmental impairments remain elusive. Emerging evidence indicates disruptions in the maternal gut microbiome can alter offspring development; however, no studies to date have examined the impact of maternal opioid treatment on maternal-offspring microbiome dysbiosis. Methods A mouse model of prenatal methadone exposure (PME) was employed to assess the impact of maternal opioid treatment on the microbiome of methadone-treated dams (MD) and their offspring. Fecal samples were collected from dams (n = 8 per treatment), one male and one female offspring per dam (n = 8 offspring per sex per treatment) for 16S rRNA sequencing. Results Methadone treatment significantly increased the microbial diversity and led to an expansion in family level bacterial abundance. Correlational analysis revealed significant positive associations between dam and offspring measures of diversity indicating methadone-induced shifts in the microbial communities are shared between dam and offspring. Sixteen features in dams and 10 features in offspring were significantly differentially abundant between treatment groups with many features corresponding to the Lachnospiraceae NK4A136 genus. Of the six features identified as differentially abundant in both MD and PME offspring, all were assigned to the Lachnospiraceae NK4A136 group, and the abundances demonstrated strong positive correlations between dam and offspring. Conclusions These preliminary findings indicate that maternal opioid treatment during pregnancy alters the composition of the maternal microbiome, and this opioid-induced shift is similarly observed in offspring which could contribute to the impaired developmental phenotypes previously described.

The aim of the study was to investigate the effect of astaxanthin and its possible mechanisms on acute cerebral infarction (ACI) in rat model. Male Sprague Dawley rats were randomly divided into sham group, model group, and astaxanthin-treated groups (20, 40, and 80 mg/kg). Neurological examination, the ratio of cerebral edema, and histopathology changes were assessed. Moreover, some oxidative stress markers were obtained for biochemical analysis, and the expression of neurotrophic factors gene was detected by real-time polymerase chain reaction (RT-PCR) method. The results showed that treatment with astaxanthin notably reduced neurological deficit scores and the ratio of cerebral edema compared with the model group. Meanwhile, astaxanthin increased the activity of catalase, superoxide dismutase, and glutathioneperoxidase as well as decreased the content of malondialdehyde in brain tissue. RT-PCR results showed that the expression of brain-derived neurotrophic factor (BDNF) and nerve growth factor (NGF) mRNA were increased with astaxanthin treatment. The results indicated that astaxanthin could ameliorate ACI followed by suppressing oxidative stress and upregulating the expression of BDNF and NGF mRNA.

Long non-coding RNAs (lncRNAs) are continuously transcribed and are involved in various cellular activities. However, their contributions to the occurrence and development of germinal center B-cell (GCB)-like diffuse large B-cell lymphoma (DLBCL) remain largely unknown. We applied microarray technology to profile the expression of lncRNAs in two different GCB-DLBCL cell lines (OCI-ly1 and OCI-ly19) and normal B lymphocytes. We demonstrated that 21,539 lncRNAs were expressed in all of the samples analyzed. This included 1,648 lncRNAs that showed a ≥2-fold upregulation and 2,671 lncRNAs that displayed a ≥2-fold downregulation in tumor cell lines (P

Multiple myeloma (MM) is a common hematological malignancy. To identify key genes and microRNAs in MM, we downloaded two gene expression profiles (GSE16558 and GSE47552) and two microRNA expression profiles (GSE17498 and GSE16558) from the Gene Expression Omnibus (GEO) database. A total of 596 differentially expressed genes (DEGs) and 39 differentially expressed microRNAs (DEMs) were screened out. Pathway analysis showed that upregulated genes were mainly enriched in the "B cell receptor signaling pathway", "Cell cycle" and "NF-kappa B signaling pathway", whereas downregulated genes were mainly enriched in the "Ribosome", "FoxO signaling pathway" and "p53 signaling pathway". We subsequently constructed a protein-protein interaction network of DEGs consisting of 277 genes and 563 interactions. In addition, 32 genes with high degrees in the network were identified as hub genes in MM, e.g. HDAC2, RBBP4, CREB1, and RB1. Additionally, we constructed a microRNA-mRNA regulatory network depicting interactions between DEMs and their targets, including the miR-135b-GADD45A and miR-148a-USPL1 pairs. In conclusion, the results of this data mining and integration help reveal the molecular basis of MM pathogenesis as well as potential biomarkers and therapeutic targets for MM diagnosis and treatment.

Myasthenia gravis (MG) is an autoimmune disorder resulting from antibodies against the proteins at the neuromuscular junction. Emerging evidence indicates that long non-coding RNAs (lncRNAs), acting as competing endogenous RNAs (ceRNAs), are involved in various diseases. However, the regulatory mechanisms of ceRNAs underlying MG remain largely unknown. In this study, we constructed a lncRNA-mediated ceRNA network involved in MG using a multi-step computational strategy. Functional annotation analysis suggests that these lncRNAs may play crucial roles in the immunological mechanism underlying MG. Importantly, through manual literature mining, we found that lncRNA SNHG16 (small nucleolar RNA host gene 16), acting as a ceRNA, plays important roles in the immune processes. Further experiments showed that SNHG16 expression was upregulated in peripheral blood mononuclear cells (PBMCs) from MG patients compared to healthy controls. Luciferase reporter assays confirmed that SNHG16 is a target of the microRNA (miRNA) let-7c-5p. Subsequent experiments indicated that SNHG16 regulates the expression of the key MG gene interleukin (IL)-10 by sponging let-7c-5p in a ceRNA manner. Furthermore, functional assays showed that SNHG16 inhibits Jurkat cell apoptosis and promotes cell proliferation by sponging let-7c-5p. Our study will contribute to a deeper understanding of the regulatory mechanism of MG and will potentially provide new therapeutic targets for MG patients.

RNA Polymerase II (RNAPII) transcription termination is regulated by the phosphorylation status of the C-terminal domain (CTD). Using disruption-compensation (DisCo) protein-protein interaction network analysis, interaction changes were observed within the termination machinery as a consequence of deletion of the serine 5 RNAPII CTD phosphatase Rtr1. Interactions between RNAPII and the cleavage factor IA (CF1A) subunit Pcf11 were reduced in rtr1Δ, whereas interactions with the CTD and RNA-binding termination factor Nrd1 were increased. These changes could be the result of altered interactions between the termination machinery and/or increased levels of premature termination of RNAPII. Transcriptome analysis in rtr1Δ cells found decreased pervasive transcription and a shift in balance of expression of sense and antisense transcripts. Globally, rtr1Δ leads to decreases in noncoding RNAs that are linked to the Nrd1, Nab3 and Sen1 (NNS)-dependent RNAPII termination pathway. Genome-wide analysis of RNAPII and Nrd1 occupancy suggests that loss of RTR1 leads to increased termination at noncoding genes and increased efficiency of snRNA termination. Additionally, premature termination increases globally at protein-coding genes where NNS is recruited during early elongation. The effects of rtr1Δ on RNA expression levels were erased following deletion of the exosome subunit Rrp6, which works with the NNS complex to rapidly degrade terminated noncoding RNAs. Overall, these data suggest that Rtr1 restricts the NNS-dependent termination pathway in WT cells to prevent premature RNAPII termination of mRNAs and ncRNAs. Additionally, Rtr1 phosphatase activity facilitates low-level elongation of noncoding transcripts that impact the transcriptome through RNAPII interference.AUTHOR SUMMARYMany cellular RNAs including those that encode for proteins are produced by the enzyme RNA Polymerase II. In this work, we have defined a new role for the phosphatase Rtr1 in the regulation of RNA Polymerase II progression from the start of transcription to the 3’ end of the gene where the nascent RNA from protein-coding genes is typically cleaved and polyadenylated. Deletion of the gene that encodes RTR1 leads to changes in the interactions between RNA polymerase II and the termination machinery. Rtr1 loss also causes early termination of RNA Polymerase II at many of its target gene types including protein coding genes and noncoding RNAs. Evidence suggests that the premature termination observed in RTR1 knockout cells occurs through the termination factor and RNA binding protein Nrd1 and its binding partner Nab3. Additionally, many of the prematurely terminated noncoding RNA transcripts are degraded by the Rrp6-containing nuclear exosome, a known component of the Nrd1-Nab3 termination coupled RNA degradation pathway. These findings suggest that Rtr1 normally promotes elongation of RNA Polymerase II transcripts through preventation of Nrd1-directed termination.

Dual RNA sequencing (RNA-seq) offers the promise of determining an interactome at a transcriptional level between a bacterium and the host but has yet to be done on any bacterial infection in human tissue. We performed dual RNA-seq and metabolomics analyses on wounded and infected sites following experimental infection of the arm with H. ducreyi. Our results suggest that H. ducreyi survives in an abscess by utilizing l-ascorbate as an alternative carbon source, possibly taking advantage of host ascorbic acid recycling, and that H. ducreyi also adapts by upregulating genes involved in anaerobic metabolism and inorganic ion and nutrient transport. To our knowledge, this is the first description of an interaction network between a bacterium and the human host at a site of infection., A major gap in understanding infectious diseases is the lack of information about molecular interaction networks between pathogens and the human host. Haemophilus ducreyi causes the genital ulcer disease chancroid in adults and is a leading cause of cutaneous ulcers in children in the tropics. We developed a model in which human volunteers are infected on the upper arm with H. ducreyi until they develop pustules. To define the H. ducreyi and human interactome, we determined bacterial and host transcriptomic and host metabolomic changes in pustules. We found that in vivo H. ducreyi transcripts were distinct from those in the inocula, as were host transcripts in pustule and wounded control sites. Many of the upregulated H. ducreyi genes were found to be involved in ascorbic acid and anaerobic metabolism and inorganic ion/nutrient transport. The top 20 significantly expressed human pathways showed that all were involved in immune responses. We generated a bipartite network for interactions between host and bacterial gene transcription; multiple positively correlated networks contained H. ducreyi genes involved in anaerobic metabolism and host genes involved with the immune response. Metabolomic studies showed that pustule and wounded samples had different metabolite compositions; the top ion pathway involved ascorbate and aldarate metabolism, which correlated with the H. ducreyi transcriptional response and upregulation of host genes involved in ascorbic acid recycling. These data show that an interactome exists between H. ducreyi and the human host and suggest that H. ducreyi exploits the metabolic niche created by the host immune response.

Clonal hematopoiesis of indeterminate potential (CHIP) increases with age and is associated with increased risks of hematological malignancies. WhileTP53mutations have been identified in CHIP, the molecular mechanisms by which mutant p53 promotes hematopoietic stem and progenitor cell (HSPC) expansion are largely unknown. We discovered that mutant p53 confers a competitive advantage to HSPCs following transplantation and promotes HSPC expansion after radiation-induced stress. Mechanistically, mutant p53 interacts with EZH2 and enhances its association with the chromatin, thereby increasing the levels of H3K27me3 in genes regulating HSPC self-renewal and differentiation. Further, genetic and pharmacological inhibition of EZH2 decrease the repopulating potential of p53 mutant HSPCs. Thus, we have uncovered an epigenetic mechanism by which mutant p53 drives clonal hematopoiesis. Our work will likely establish epigenetic regulator EZH2 as a novel therapeutic target for preventing CHIP progression and treating hematological malignancies withTP53mutations.

Clonal hematopoiesis of indeterminate potential (CHIP) increases with age and is associated with increased risks of hematological malignancies. While TP53 mutations have been identified in CHIP, the molecular mechanisms by which mutant p53 promotes hematopoietic stem and progenitor cell (HSPC) expansion are largely unknown. Here we discover that mutant p53 confers a competitive advantage to HSPCs following transplantation and promotes HSPC expansion after radiation-induced stress. Mechanistically, mutant p53 interacts with EZH2 and enhances its association with the chromatin, thereby increasing the levels of H3K27me3 in genes regulating HSPC self-renewal and differentiation. Furthermore, genetic and pharmacological inhibition of EZH2 decreases the repopulating potential of p53 mutant HSPCs. Thus, we uncover an epigenetic mechanism by which mutant p53 drives clonal hematopoiesis. Our work will likely establish epigenetic regulator EZH2 as a novel therapeutic target for preventing CHIP progression and treating hematological malignancies with TP53 mutations.

Transcriptome studies can identify genes whose expression differs between alcoholics and controls. To test which variants associated with alcohol use disorder (AUDs) may cause expression differences, we integrated deep RNA-seq and genome-wide association studies (GWAS) data from four postmortem brain regions of 30 AUDs subjects and 30 controls (social/non-drinkers) and analyzed allele-specific expression (ASE). We identified 90 genes with differential ASE in subjects with AUDs compared to controls. Of these, 61 genes contained 437 single nucleotide polymorphisms (SNPs) in the 3’ untranslated regions (3’UTR) with at least one heterozygote among the subjects studied. Using a modified PASSPORT-seq (parallel assessment of polymorphisms in miRNA target-sites by sequencing) assay, we identified 25 SNPs that showed affected RNA levels in a consistent manner in two neuroblastoma cell lines, SH-SY5Y and SK-N-BE(2). Many of these are in binding sites of miRNAs and RNA binding proteins, indicating that these SNPs are likely causal variants of AUD-associated differential ASE.

Purpose Epithelial to mesenchymal transition (EMT) can contribute to gastric cancer (GC) progression and recurrence following therapy. Tumor-associated neutrophils (TANs) are associated with poor outcomes in a variety of cancers. However, it is not clear whether TANs interact with the EMT process during GC development. Methods Immunohistochemistry was performed to examine the distribution and levels of CD66 + neutrophils in samples from 327 patients with GC. CD66b + TANs were isolated either directly from GC cell suspensions or were conditioned from healthy donor peripheral blood polymorphonuclear neutrophils (PMNs) stimulated with tumor tissue culture supernatants (TTCS) and placed into co-culture with MKN45 or MKN74 cells, after which migration, invasion and EMT were measured. Interleukin-17a (IL-17a) was blocked with a polyclonal antibody, and the STAT3 pathway was blocked with the specific inhibitor AG490. Results Neutrophils were widely distributed in gastric tissues of patients with GC and were enriched predominantly at the invasion margin. Neutrophil levels at the invasion margin were an independent predictor of poor disease-free survival (DFS) and disease-specific survival (DSS). IL-17a + neutrophils constituted a large portion of IL-17a-producing cells in GC, and IL-17a was produced at the highest levels in co-culture compared with that in TANs not undergoing co-culture. TANs enhanced the migration, invasion and EMT of GC cells through the secretion of IL-17a, which activated the Janus kinase 2/signal transducers and activators of transcription (JAK2/STAT3) pathway in GC cells, while deprivation of IL-17a using a neutralizing antibody or inhibition of the JAK2/STAT3 pathway with AG490 markedly reversed these TAN-induced phenotypes in GC cells induced by TANs. Conclusions Neutrophils correlate with tumor stage and predict poor prognosis in GC. TANs produce IL-17a, which promotes EMT of GC cells through JAK2/STAT3 signalling. Blockade of IL-17a signalling with a neutralizing antibody inhibits TAN-stimulated activity in GC cells. Therefore, IL-17a-targeted therapy might be used to treat patients with GC. Electronic supplementary material The online version of this article (10.1186/s13046-018-1003-0) contains supplementary material, which is available to authorized users.

Summary Single-cell RNA sequencing (scRNA-seq) is an evolving technology used to elucidate the cellular architecture of adult organs. Previous scRNA-seq on breast tissue utilized reduction mammoplasty samples, which are often histologically abnormal. We report a rapid tissue collection/processing protocol to perform scRNA-seq of breast biopsies of healthy women and identify 23 breast epithelial cell clusters. Putative cell-of-origin signatures derived from these clusters are applied to analyze transcriptomes of ~3,000 breast cancers. Gene signatures derived from mature luminal cell clusters are enriched in ~68% of breast cancers, whereas a signature from a luminal progenitor cluster is enriched in ~20% of breast cancers. Overexpression of luminal progenitor cluster-derived signatures in HER2+, but not in other subtypes, is associated with unfavorable outcome. We identify TBX3 and PDK4 as genes co-expressed with estrogen receptor (ER) in the normal breasts, and their expression analyses in >550 breast cancers enable prognostically relevant subclassification of ER+ breast cancers., Graphical abstract, Highlights Healthy breast contains 23 subclusters of epithelial cells Breast cancers may originate from 3 luminal mature and 1 progenitor subclusters TBX3 and PDK4, co-expressed with estrogen receptor (ER), subclassify ER+ breast cancers, Bhat-Nakshatri et al. describe the single-cell atlas and document 23 epithelial cell subclusters in the healthy breast. Although experimentally validating cell of origin of a tumor is technically challenging, overlap analysis of subcluster-enriched gene signatures with breast tumor transcriptomes revealed dominant representation of differentiated luminal subclusters-derived signatures in breast cancers.

IMPACT: Screening the effect of thousands of non-coding genetic variants will help identify variants important in the etiology of diseases OBJECTIVES/GOALS: Massively parallel reporter assays (MPRAs) can experimentally evaluate the impact of genetic variants on gene expression. In this study, our objective was to systematically evaluate the functional activity of 3’-UTR SNPs associated with neurological disorders and use those results to help understand their contributions to disease etiology. METHODS/STUDY POPULATION: To choose variants to evaluate with the MPRA, we first gathered SNPs from the GWAS Catalog that were associated with any neurological disorder trait with p-value < 10-5. For each SNP, we identified the region that was in linkage disequilibrium (r2 > 0.8) and retrieved all the common 3’-UTR SNPs (allele-frequency > 0.05) within that region. We used an MPRA to measure the impact of these 3’-UTR variants in SH-SY5Y neuroblastoma cells and a microglial cell line. These results were then used to train a deep-learning model to predict the impact of variants and identify features that contribute to the predictions. RESULTS/ANTICIPATED RESULTS: Of the 13,515 3’-UTR SNPs tested, 400 and 657 significantly impacted gene expression in SH-SY5Y and microglia, respectively. Of the 84 SNPs significantly impacted in both cells, the direction of impact was the same in 81. The direction of eQTL in GTEx tissues agreed with the assay SNP effect in SH-SY5Y cells but not microglial cells. The deep-learning model predicted sequence activity level correlated with the experimental activity level (Spearman’s corr = 0.45). The deep-learning model identified several predictive motifs similar to motifs of RNA-binding proteins. DISCUSSION/SIGNIFICANCE OF FINDINGS: This study demonstrates that MPRAs can be used to evaluate the effect of non-coding variants, and the results can be used to train a machine learning model and interpret its predictions. Together, these can help identify causal variants and further understand the etiology of diseases.

The aim of the present study was to analyze the selenoprotein expression levels in gastric cancer patients. We enrolled 40 patients (29 males, 11 females) who were recently diagnosed with gastric cancer and 50 healthy people (30 males, 20 females) as controls. The expression of 25 selenoprotein genes (Dio1, Dio2, Dio3, Gpx1, Gpx2, Gpx3, Gpx4, Gpx6, SelH, SelI, SelK, SelM, SelN, SelO, SelP, SelS, SelT, SelV, SelW, SelX, Sel15, Sps2, TR1, TR2, and TR3) in human gastric cancer tissues, para-carcinoma tissues, adjacent normal gastric tissues, erythrocytes, and lymphocytes in the gastric cancer group and healthy control group was analyzed by qRT-PCR. Here, we showed that among the 25 selenoproteins, 13 selenoproteins in erythrocytes (Gpx1, Gpx4, Sel15, TR1, TR2, SelH, SelK, SelM, SelO, SelS, SelV, SelW, and Sps2), 15 selenoproteins in lymphocytes (Gpx1, Gpx4, Sel15, TR1, TR2, SelH, SelK, SelN, SelO, SelS, SelT, SelV, SelX, SelW, and Sps2) and 13 selenoproteins in gastric cancer and para-carcinoma tissues (Dio1, Dio2, Dio3, Gpx1, Gpx4, Sel15, SelH, SelK, SelM, SelS, SelT, SelW, and Sps2) were significantly decreased (P

Cachexia affects the majority of cancer patients, with currently no effective treatments. Cachexia is defined by increased fatigue and loss of muscle function resulting from muscle and fat depletion. Previous studies suggest that chemotherapy may contribute to cachexia, although the causes responsible for this association are not clear. The purpose of this study was to investigate the mechanism(s) associated with chemotherapy-related effects on body composition and muscle function. Normal mice were administered chemotherapy regimens used for the treatment of colorectal cancer, such as Folfox (5-FU, leucovorin, oxaliplatin) or Folfiri (5-FU, leucovorin, irinotecan) for 5 weeks. The animals that received chemotherapy exhibited concurrent loss of muscle mass and muscle weakness. Consistently with previous findings, muscle wasting was associated with up-regulation of ERK1/2 and p38 MAPKs. No changes in ubiquitin-dependent proteolysis or in the expression of TGFβ-family members were detected. Further, marked decreases in mitochondrial content, associated with abnormalities at the sarcomeric level and with increase in the number of glycolytic fibers were observed in the muscle of mice receiving chemotherapy. Finally, ACVR2B/Fc or PD98059 prevented Folfiri-associated ERK1/2 activation and myofiber atrophy in C2C12 cultures. Our findings demonstrate that chemotherapy promotes MAPK-dependent muscle atrophy as well as mitochondrial depletion and alterations of the sarcomeric units. Therefore, these findings suggest that chemotherapy potentially plays a causative role in the occurrence of muscle loss and weakness. Moreover, the present observations provide a strong rationale for testing ACVR2B/Fc or MEK1 inhibitors in combination with anticancer drugs as novel strategies aimed at preventing chemotherapy-associated muscle atrophy.

Haemophilus ducreyi causes the sexually transmitted disease chancroid in adults and cutaneous ulcers in children. In humans, H. ducreyi resides in an abscess and must adapt to a variety of stresses. Previous studies (D. Gangaiah, M. Labandeira-Rey, X. Zhang, K. R. Fortney, S. Ellinger, B. Zwickl, B. Baker, Y. Liu, D. M. Janowicz, B. P. Katz, C. A. Brautigam, R. S. Munson, Jr., E. J. Hansen, and S. M. Spinola, mBio 5:e01081-13, 2014, http://dx.doi.org/10.1128/mBio.01081-13 ) suggested that H. ducreyi encounters growth conditions in human lesions resembling those found in stationary phase. However, how H. ducreyi transcriptionally responds to stress during human infection is unknown. Here, we determined the H. ducreyi transcriptome in biopsy specimens of human lesions and compared it to the transcriptomes of bacteria grown to mid-log, transition, and stationary phases. Multidimensional scaling showed that the in vivo transcriptome is distinct from those of in vitro growth. Compared to the inoculum (mid-log-phase bacteria), H. ducreyi harvested from pustules differentially expressed ∼93 genes, of which 62 were upregulated. The upregulated genes encode homologs of proteins involved in nutrient transport, alternative carbon pathways ( l -ascorbate utilization and metabolism), growth arrest response, heat shock response, DNA recombination, and anaerobiosis. H. ducreyi upregulated few genes ( hgbA , flp-tad , and lspB-lspA2 ) encoding virulence determinants required for human infection. Most genes regulated by CpxRA, RpoE, Hfq, (p)ppGpp, and DksA, which control the expression of virulence determinants and adaptation to a variety of stresses, were not differentially expressed in vivo , suggesting that these systems are cycling on and off during infection. Taken together, these data suggest that the in vivo transcriptome is distinct from those of in vitro growth and that adaptation to nutrient stress and anaerobiosis is crucial for H. ducreyi survival in humans.

Lysine acetylation is a critical posttranslational modification that influences protein activity, stability, and binding properties. The acetylation of histone proteins in particular is a well-characterized feature of gene expression regulation. In the protozoan parasite Toxoplasma gondii , a number of lysine acetyltransferases (KATs) contribute to gene expression and are essential for parasite viability. The natural product garcinol was recently reported to inhibit enzymatic activities of GCN5 and p300 family KATs in other species. Here we show that garcinol inhibits TgGCN5b, the only nuclear GCN5 family KAT known to be required for Toxoplasma tachyzoite replication. Treatment of tachyzoites with garcinol led to a reduction of global lysine acetylation, particularly on histone H3 and TgGCN5b itself. We also performed transcriptome sequencing (RNA-seq), which revealed increasing aberrant gene expression coincident with increasing concentrations of garcinol. The majority of the genes that were most significantly affected by garcinol were also associated with TgGCN5b in a previously reported chromatin immunoprecipitation assay with microarray technology (ChIP-chip) analysis. The dysregulated gene expression induced by garcinol significantly inhibits Toxoplasma tachyzoite replication, and the concentrations used exhibit no overt toxicity on human host cells. Garcinol also inhibits Plasmodium falciparum asexual replication with a 50% inhibitory concentration (IC 50 ) similar to that for Toxoplasma . Together, these data support that pharmacological inhibition of TgGCN5b leads to a catastrophic failure in gene expression control that prevents parasite replication.

The thiazolidinedione pioglitazone, which is also a PPAR-γ agonist, now is widely used in patients with hypercholesterolemia and hypertriglyceridemia. NF-κB is a ubiquitously expressed transcription factor controlling the expression of numerous genes involved in inflammation. The aim of the present study was to evaluate whether the activation of PPAR-γ attenuates the cisplatin-induced NF-κB activation in cisplatin nephrotoxicity. The results showed that the PPAR-γ agonist pioglitazone decreased the expression of NF-κB p65 transcription target genes (e.g., IL-6, IL-1β, and TNF-α) and inhibited histological injury and inflammatory cells infiltration in cisplatin nephrotoxicity. The suppression of NF-κB activity following pioglitazone treatment inhibited the cisplatin-induced IκB-α degredation and NF-κB p65 subunit translocation. Translocation of the NF-κB p65 subunit depends on p65 acetylation, which primarily regulated by SIRT1 or p300. Notably, AMP kinase (AMPK) activation not only decreased the phosphorylation, activation and p65 interaction of p300 but also increased SIRT1 expression, activation and p65 binding, thus leading to a significant reduction in p65 acetylation. Interestingly, the reduction of IL-6, TNF-α and IL-1β, the inhibition of histological injury and the inflammatory cells infiltration following pioglitazone treatment in cisplatin nephrotoxicity were attenuated after treatment with the PPAR-γ antagonist GW9662. These results suggest that the PPAR-γ agonist pioglitazone prevents NF-κB activation in cisplatin nephrotoxicity through a reduction in p65 acetylation via the AMPK-SIRT1/p300 pathway.

Obesity is closely associated with adipose tissue inflammation and insulin resistance. Dysglycemia and type 2 diabetes results when islet &beta, cells fail to maintain appropriate insulin secretion in the face of insulin resistance. To clarify the early transcriptional events leading to &beta, cell failure in the setting of obesity, we fed male C57BL/6J mice an obesogenic, high-fat diet (60% kcal from fat) or a control diet (10% kcal from fat) for one week, and islets from these mice (from four high-fat- and three control-fed mice) were subjected to single-cell RNA sequencing (sc-RNAseq) analysis. Islet endocrine cell types (&alpha, cells, &beta, cells, &delta, cells, PP cells) and other resident cell types (macrophages, T cells) were annotated by transcript profiles and visualized using Uniform Manifold Approximation and Projection for Dimension Reduction (UMAP) plots. UMAP analysis revealed distinct cell clusters (11 for &beta, cells, 5 for &alpha, cells, 3 for &delta, cells, PP cells, ductal cells, endothelial cells), emphasizing the heterogeneity of cell populations in the islet. Collectively, the clusters containing the majority of &beta, cells showed the fewest gene expression changes, whereas clusters harboring the minority of &beta, cells showed the most changes. We identified that distinct &beta, cell clusters downregulate genes associated with the endoplasmic reticulum stress response and upregulate genes associated with insulin secretion, whereas others upregulate genes that impair insulin secretion, cell proliferation, and cell survival. Moreover, all &beta, cell clusters negatively regulate genes associated with immune response activation. Glucagon-producing &alpha, cells exhibited patterns similar to &beta, cells but, again, in clusters containing the minority of &alpha, cells. Our data indicate that an early transcriptional response in islets to an obesogenic diet reflects an attempt by distinct populations of &beta, cells to augment or impair cellular function and/or reduce inflammatory responses as possible harbingers of ensuing insulin resistance.

Summary Brain metastasis (br-met) develops in an immunologically unique br-met niche. Central nervous system-native myeloid cells (CNS-myeloids) and bone-marrow-derived myeloid cells (BMDMs) cooperatively regulate brain immunity. The phenotypic heterogeneity and specific roles of these myeloid subsets in shaping the br-met niche to regulate br-met outgrowth have not been fully revealed. Applying multimodal single-cell analyses, we elucidated a heterogeneous but spatially defined CNS-myeloid response during br-met outgrowth. We found Ccr2+ BMDMs minimally influenced br-met while CNS-myeloid promoted br-met outgrowth. Additionally, br-met-associated CNS-myeloid exhibited downregulation of Cx3cr1. Cx3cr1 knockout in CNS-myeloid increased br-met incidence, leading to an enriched interferon response signature and Cxcl10 upregulation. Significantly, neutralization of Cxcl10 reduced br-met, while rCxcl10 increased br-met and recruited VISTAHi PD-L1+ CNS-myeloid to br-met lesions. Inhibiting VISTA- and PD-L1-signaling relieved immune suppression and reduced br-met burden. Our results demonstrate that loss of Cx3cr1 in CNS-myeloid triggers a Cxcl10-mediated vicious cycle, cultivating a br-met-promoting, immune-suppressive niche.

Clonal hematopoiesis of indeterminate potential (CHIP) increases with age and is associated with increased risks of hematological malignancies. While TP53 mutations have been identified in CHIP, the molecular mechanisms by which mutant p53 promotes hematopoietic stem and progenitor cell (HSPC) expansion are largely unknown. Here we discover that mutant p53 confers a competitive advantage to HSPCs following transplantation and promotes HSPC expansion after radiation-induced stress. Mechanistically, mutant p53 interacts with EZH2 and enhances its association with the chromatin, thereby increasing the levels of H3K27me3 in genes regulating HSPC self-renewal and differentiation. Furthermore, genetic and pharmacological inhibition of EZH2 decreases the repopulating potential of p53 mutant HSPCs. Thus, we uncover an epigenetic mechanism by which mutant p53 drives clonal hematopoiesis. Our work will likely establish epigenetic regulator EZH2 as a novel therapeutic target for preventing CHIP progression and treating hematological malignancies with TP53 mutations., Ageing is associated with clonal hematopoiesis of indeterminate potential (CHIP), which is linked to increased risks of hematological malignancies. Here the authors uncover an epigenetic mechanism through which mutant p53 drives clonal hematopoiesis through interaction with EZH2.

In response to vascular injury, vascular smooth muscle cells (VSMCs) may switch from a contractile to a proliferative phenotype thereby contributing to neointima formation. Previous studies showed that the long noncoding RNA (lncRNA) NEAT1 is critical for paraspeckle formation and tumorigenesis by promoting cell proliferation and migration. However, the role of NEAT1 in VSMC phenotypic modulation is unknown. Herein we showed that NEAT1 expression was induced in VSMCs during phenotypic switching in vivo and in vitro. Silencing NEAT1 in VSMCs resulted in enhanced expression of SM-specific genes while attenuating VSMC proliferation and migration. Conversely, overexpression of NEAT1 in VSMCs had opposite effects. These in vitro findings were further supported by in vivo studies in which NEAT1 knockout mice exhibited significantly decreased neointima formation following vascular injury, due to attenuated VSMC proliferation. Mechanistic studies demonstrated that NEAT1 sequesters the key chromatin modifier WDR5 (WD Repeat Domain 5) from SM-specific gene loci, thereby initiating an epigenetic “off” state, resulting in down-regulation of SM-specific gene expression. Taken together, we demonstrated an unexpected role of the lncRNA NEAT1 in regulating phenotypic switching by repressing SM-contractile gene expression through an epigenetic regulatory mechanism. Our data suggest that NEAT1 is a therapeutic target for treating occlusive vascular diseases.

Next-generation sequencing (NGS) studies have identified large numbers of genetic variants that are predicted to alter miRNA–mRNA interactions. We developed a novel high-throughput bioassay, PASSPORT-seq, that can functionally test in parallel 100s of these variants in miRNA binding sites (mirSNPs). The results are highly reproducible across both technical and biological replicates. The utility of the bioassay was demonstrated by testing 100 mirSNPs in HEK293, HepG2, and HeLa cells. The results of several of the variants were validated in all three cell lines using traditional individual luciferase assays. Fifty-five mirSNPs were functional in at least one of three cell lines (FDR ≤ 0.05); 11, 36, and 27 of them were functional in HEK293, HepG2, and HeLa cells, respectively. Only four of the variants were functional in all three cell lines, which demonstrates the cell-type specific effects of mirSNPs and the importance of testing the mirSNPs in multiple cell lines. Using PASSPORT-seq, we functionally tested 111 variants in the 3′ UTR of 17 pharmacogenes that are predicted to alter miRNA regulation. Thirty-three of the variants tested were functional in at least one cell line.

B-cell translocation gene 1 (BTG1) is a member of the BTG/transducer of Erb family. The present study evaluated the impact of BTG1 gene expression on the clinical outcome of diffuse large B-cell lymphoma (DLBCL) and investigated potential mechanisms using the Gene Expression Omnibus (GEO) database. The gene expression profile datasets GSE31312, GSE10846, GSE65420 and GSE87371 were downloaded from the GEO database. BTG1 expression and clinicopathological data were obtained from the GSE31312 dataset. In 498 cases, the expression of BTG1 in DLBCL was associated with treatment response (χ2=19.020; P

CpxRA is an envelope stress response system found in all members of the family Enterobacteriaceae ; CpxA has kinase activity for CpxR and phosphatase activity for phospho-CpxR, a transcription factor. CpxR also accepts phosphate groups from acetyl phosphate, a glucose metabolite. Activation of CpxR increases the transcription of genes encoding membrane repair and downregulates virulence determinants. We hypothesized that activation of CpxR could serve as an antimicrobial/antivirulence strategy and discovered compounds that activate CpxR in Escherichia coli by inhibiting CpxA phosphatase activity. As a prelude to testing such compounds in vivo , here we constructed cpxA (in the presence of glucose, CpxR is activated because of a lack of CpxA phosphatase) and cpxR (system absent) deletion mutants of uropathogenic E. coli (UPEC) CFT073. By RNA sequencing, few transcriptional differences were noted between the cpxR mutant and its parent, but in the cpxA mutant, several UPEC virulence determinants were downregulated, including the fim and pap operons, and it exhibited reduced mannose-sensitive hemagglutination of guinea pig red blood cells in vitro . In competition experiments with mice, both mutants were less fit than the parent in the urine, bladder, and kidney; these fitness defects were complemented in trans . Unexpectedly, in single-strain challenges, only the cpxA mutant was attenuated for virulence in the kidney but not in the bladder or urine. For the cpxA mutant, this may be due to the preferential use of amino acids over glucose as a carbon source in the bladder and urine by UPEC. These studies suggest that CpxA phosphatase inhibitors may have some utility for treating complex urinary tract infections.

During infection, Neisseria gonorrhoeae senses and responds to stress; such responses may be modulated by MisRS (NGO0177 and NGO0176), a two-component system that is a homolog of CpxRA. In Escherichia coli , CpxRA senses and responds to envelope stress; CpxA is a sensor kinase/phosphatase for CpxR, a response regulator. When a cpxA mutant is grown in medium containing glucose, CpxR is phosphorylated by acetyl phosphate but cannot be dephosphorylated, resulting in constitutive activation. Kandler and coworkers (J. L. Kandler, C. L. Holley, J. L. Reimche, V. Dhulipala, J. T. Balthazar, A. Muszyński, R. W. Carlson, and W. M. Shafer, Antimicrob Agents Chemother 60:4690–4700, 2016, https://doi.org/10.1128/AAC.00823-16 ) showed that MisR (CpxR) is required for the maintenance of membrane integrity and resistance to antimicrobial peptides, suggesting a role in gonococcal survival in vivo . Here, we evaluated the contributions of MisR and MisS (CpxA) to gonococcal infection in a murine model of cervicovaginal colonization and identified MisR-regulated genes using RNA sequencing (RNA-Seq). The deletion of misR or misS severely reduced the capacity of N. gonorrhoeae to colonize mice or maintain infection over a 7-day period and reduced microbial fitness after exposure to heat shock. Compared to the wild type (WT), the inactivation of misR identified 157 differentially regulated genes, most of which encoded putative envelope proteins. The inactivation of misS identified 17 differentially regulated genes compared to the WT and 139 differentially regulated genes compared to the misR mutant, 111 of which overlapped those differentially expressed in the comparison of the WT versus the misR mutant. These data indicate that an intact MisRS system is required for gonococcal infection of mice. Provided the MisR is constitutively phosphorylated in the misS mutant, the data suggest that controlled but not constitutive activation is required for gonococcal infection in mice.

Integration of viral vectors into a host genome is associated with insertional mutagenesis and subjects in clinical gene therapy trials must be monitored for this adverse event. Several PCR based methods such as ligase-mediated (LM) PCR, linear-amplification-mediated (LAM) PCR and non-restrictive (nr) LAM PCR were developed to identify sites of vector integration. Coupling the power of next-generation sequencing technologies with various PCR approaches will provide a comprehensive and genome-wide profiling of insertion sites and increase throughput. In this bioinformatics study, we aimed to develop and apply quality metrics to viral insertion data obtained using next-generation sequencing. We developed five simple metrics for assessing next-generation sequencing data from different PCR products and showed how the metrics can be used to objectively compare runs performed with the same methodology as well as data generated using different PCR techniques. The results will help researchers troubleshoot complex methodologies, understand the quality of sequencing data, and provide a starting point for developing standardization of vector insertion site data analysis.

A variety of microbial communities and their genes (the microbiome) exist throughout the human body, with fundamental roles in human health and disease. The National Institutes of Health (NIH)-funded Human Microbiome Project Consortium has established a population-scale framework to develop metagenomic protocols, resulting in a broad range of quality-controlled resources and data including standardized methods for creating, processing and interpreting distinct types of high-throughput metagenomic data available to the scientific community. Here we present resources from a population of 242 healthy adults sampled at 15 or 18 body sites up to three times, which have generated 5,177 microbial taxonomic profiles from 16S ribosomal RNA genes and over 3.5 terabases of metagenomic sequence so far. In parallel, approximately 800 reference strains isolated from the human body have been sequenced. Collectively, these data represent the largest resource describing the abundance and variety of the human microbiome, while providing a framework for current and future studies.