In situ Patch-seq analysis of microglia reveals a lack of stress genes as found in FACS-isolated microglia.

We applied the patch-seq technique to harvest transcripts from individual microglial cells from cortex, hippocampus and corpus callosum of acute brain slices from adult mice. After recording membrane currents with the patch-clamp technique, the cytoplasm was collected via the pipette and underwent adapted SMART-seq2 preparation with subsequent sequencing. On average, 4138 genes were detected in 113 cells from hippocampus, corpus callosum and cortex, including microglia markers such as Tmem119, P2ry12 and Siglec-H. Comparing our dataset to previously published single cell mRNA sequencing data from FACS-isolated microglia indicated that two clusters of cells were absent in our patch-seq dataset. Pathway analysis of marker genes in FACS-specific clusters revealed association with microglial activation and stress response. This indicates that under normal conditions microglia in situ lack transcripts associated with a stress-response, and that the microglia-isolation procedure by mechanical dissociation and FACS triggers the expression of genes related to activation and stress.

High Yap and Mll1 promote a persistent regenerative cell state induced by Notch signaling and loss of p53.

Specified intestinal epithelial cells reprogram and contribute to the regeneration and renewal of the epithelium upon injury. Mutations that deregulate such renewal processes may contribute to tumorigenesis. Using intestinal organoids, we show that concomitant activation of Notch signaling and ablation of p53 induce a highly proliferative and regenerative cell state, which is associated with increased levels of Yap and the histone methyltransferase Mll1. The induced signaling system orchestrates high proliferation, self-renewal, and niche-factor-independent growth, and elevates the trimethylation of histone 3 at lysine 4 (H3K4me3). We demonstrate that Yap and Mll1 are also elevated in patient-derived colorectal cancer (CRC) organoids and control growth and viability. Our data suggest that Notch activation and p53 ablation induce a signaling circuitry involving Yap and the epigenetic regulator Mll1, which locks cells in a proliferative and regenerative state that renders them susceptible for tumorigenesis.

Dissecting Herpes Simplex Virus 1-Induced Host Shutoff at the RNA Level.

Herpes simplex virus 1 (HSV-1) induces a profound host shutoff during lytic infection. The virion host shutoff (vhs) protein plays a key role in this process by efficiently cleaving host and viral mRNAs. Furthermore, the onset of viral DNA replication is accompanied by a rapid decline in host transcriptional activity. To dissect relative contributions of both mechanisms and elucidate gene-specific host transcriptional responses throughout the first 8 h of lytic HSV-1 infection, we used transcriptome sequencing of total, newly transcribed (4sU-labeled) and chromatin-associated RNA in wild-type (WT) and Δvhs mutant infection of primary human fibroblasts. Following virus entry, vhs activity rapidly plateaued at an elimination rate of around 30% of cellular mRNAs per hour until 8 h postinfection (p.i.). In parallel, host transcriptional activity dropped to 10 to 20%. While the combined effects of both phenomena dominated infection-induced changes in total RNA, extensive gene-specific transcriptional regulation was observable in chromatin-associated RNA and was surprisingly concordant between WT and Δvhs infections. Both induced strong transcriptional upregulation of a small subset of genes that were poorly expressed prior to infection but already primed by H3K4me3 histone marks at their promoters. Most interestingly, analysis of chromatin-associated RNA revealed vhs-nuclease-activity-dependent transcriptional downregulation of at least 150 cellular genes, in particular of many integrin adhesome and extracellular matrix components. This was accompanied by a vhs-dependent reduction in protein levels by 8 h p.i. for many of these genes. In summary, our study provides a comprehensive picture of the molecular mechanisms that govern cellular RNA metabolism during the first 8 h of lytic HSV-1 infection.IMPORTANCE The HSV-1 virion host shutoff (vhs) protein efficiently cleaves both host and viral mRNAs in a translation-dependent manner. In this study, we model and quantify changes in vhs activity, as well as virus-induced global loss of host transcriptional activity, during productive HSV-1 infection. In general, HSV-1-induced alterations in total RNA levels were dominated by these two global effects. In contrast, chromatin-associated RNA depicted gene-specific transcriptional changes. This revealed highly concordant transcriptional changes in WT and Δvhs infections, confirmed DUX4 as a key transcriptional regulator in HSV-1 infection, and identified vhs-dependent transcriptional downregulation of the integrin adhesome and extracellular matrix components. The latter explained seemingly gene-specific effects previously attributed to vhs-mediated mRNA degradation and resulted in a concordant loss in protein levels by 8 h p.i. for many of the respective genes.

The epigenetic regulator Mll1 is required for Wnt-driven intestinal tumorigenesis and cancer stemness.

Wnt/ß-catenin signaling is crucial for intestinal carcinogenesis and the maintenance of intestinal cancer stem cells. Here we identify the histone methyltransferase Mll1 as a regulator of Wnt-driven intestinal cancer. Mll1 is highly expressed in Lgr5+ stem cells and human colon carcinomas with increased nuclear ß-catenin. High levels of MLL1 are associated with poor survival of colon cancer patients. The genetic ablation of Mll1 in mice prevents Wnt/ß-catenin-driven adenoma formation from Lgr5+ intestinal stem cells. Ablation of Mll1 decreases the self-renewal of human colon cancer spheres and halts tumor growth of xenografts. Mll1 controls the expression of stem cell genes including the Wnt/ß-catenin target gene Lgr5. Upon the loss of Mll1, histone methylation at the stem cell promoters switches from activating H3K4 tri-methylation to repressive H3K27 tri-methylation, indicating that Mll1 sustains stem cell gene expression by antagonizing gene silencing through polycomb repressive complex 2 (PRC2)-mediated H3K27 tri-methylation. Transcriptome profiling of Wnt-mutated intestinal tumor-initiating cells reveals that Mll1 regulates Gata4/6 transcription factors, known to sustain cancer stemness and to control goblet cell differentiation. Our results demonstrate that Mll1 is an essential epigenetic regulator of Wnt/ß-catenin-induced intestinal tumorigenesis and cancer stemness.

Comparative Analysis of the Secretome and Interactome of Trypanosoma cruzi and Trypanosoma rangeli Reveals Species Specific Immune Response Modulating Proteins.

Chagas disease, a zoonosis caused by the flagellate protozoan Trypanosoma cruzi, is a chronic and systemic parasitic infection that affects ~5-7 million people worldwide, mainly in Latin America. Chagas disease is an emerging public health problem due to the lack of vaccines and effective treatments. According to recent studies, several T. cruzi secreted proteins interact with the human host during cell invasion. Moreover, some comparative studies with T. rangeli, which is non-pathogenic in humans, have been performed to identify proteins directly involved in the pathogenesis of the disease. In this study, we present an integrated analysis of canonical putative secreted proteins (PSPs) from both species. Additionally, we propose an interactome with human host and gene family clusters, and a phylogenetic inference of a selected protein. In total, we identified 322 exclusively PSPs in T. cruzi and 202 in T. rangeli. Among the PSPs identified in T. cruzi, we found several trans-sialidases, mucins, MASPs, proteins with phospholipase 2 domains (PLA2-like), and proteins with Hsp70 domains (Hsp70-like) which have been previously characterized and demonstrated to be related to T. cruzi virulence. PSPs found in T. rangeli were related to protozoan metabolism, specifically carboxylases and phosphatases. Furthermore, we also identified PSPs that may interact with the human immune system, including heat shock and MASP proteins, but in a lower number compared to T. cruzi. Interestingly, we describe a hypothetical hybrid interactome of PSPs which reveals that T. cruzi secreted molecules may be down-regulating IL-17 whilst T. rangeli may enhance the production of IL-15. These results will pave the way for a better understanding of the pathophysiology of Chagas disease and may ultimately lead to the identification of molecular targets, such as key PSPs, that could be used to minimize the health outcomes of Chagas disease by modulating the immune response triggered by T. cruzi infection.

Correction to: Comprehensive gene expression meta-analysis identifies signature genes that distinguish microglia from peripheral monocytes/macrophages in health and glioma.

The original publication of this article [1] contained 3 minor errors in Figs. 1, 3 and 5. In this correction article the updated figures are published. The figure captions describe the updated information in these figures.

Comprehensive gene expression meta-analysis identifies signature genes that distinguish microglia from peripheral monocytes/macrophages in health and glioma.

Monocytes/macrophages have begun to emerge as key cellular modulators of brain homeostasis and central nervous system (CNS) disease. In the healthy brain, resident microglia are the predominant macrophage cell population; however, under conditions of blood-brain barrier leakage, peripheral monocytes/macrophages can infiltrate the brain and participate in CNS disease pathogenesis. Distinguishing these two populations is often challenging, owing to a paucity of universally accepted and reliable markers. To identify discriminatory marker sets for microglia and peripheral monocytes/macrophages, we employed a large meta-analytic approach using five published murine transcriptional datasets. Following hierarchical clustering, we filtered the top differentially expressed genes (DEGs) through a brain cell type-specific sequencing database, which led to the identification of eight microglia and eight peripheral monocyte/macrophage markers. We then validated their differential expression, leveraging a published single cell RNA sequencing dataset and quantitative RT-PCR using freshly isolated microglia and peripheral monocytes/macrophages from two different mouse strains. We further verified the translation of these DEGs at the protein level. As top microglia DEGs, we identified P2ry12, Tmem119, Slc2a5 and Fcrls, whereas Emilin2, Gda, Hp and Sell emerged as the best DEGs for identifying peripheral monocytes/macrophages. Lastly, we evaluated their utility in discriminating monocyte/macrophage populations in the setting of brain pathology (glioma), and found that these DEG sets distinguished glioma-associated microglia from macrophages in both RCAS and GL261 mouse models of glioblastoma. Taken together, this unbiased bioinformatic approach facilitated the discovery of a robust set of microglia and peripheral monocyte/macrophage expression markers to discriminate these monocyte populations in both health and disease.

The Wnt-Driven Mll1 Epigenome Regulates Salivary Gland and Head and Neck Cancer.

We identified a regulatory system that acts downstream of Wnt/ß-catenin signaling in salivary gland and head and neck carcinomas. We show in a mouse tumor model of K14-Cre-induced Wnt/ß-catenin gain-of-function and Bmpr1a loss-of-function mutations that tumor-propagating cells exhibit increased Mll1 activity and genome-wide increased H3K4 tri-methylation at promoters. Null mutations of Mll1 in tumor mice and in xenotransplanted human head and neck tumors resulted in loss of self-renewal of tumor-propagating cells and in block of tumor formation but did not alter normal tissue homeostasis. CRISPR/Cas9 mutagenesis and pharmacological interference of Mll1 at sequences that inhibit essential protein-protein interactions or the SET enzyme active site also blocked the self-renewal of mouse and human tumor-propagating cells. Our work provides strong genetic evidence for a crucial role of Mll1 in solid tumors. Moreover, inhibitors targeting specific Mll1 interactions might offer additional directions for therapies to treat these aggressive tumors.

Fructooligosaccharide intake promotes epigenetic changes in the intestinal mucosa in growing and ageing rats.

PURPOSE: The aim of this study was to investigate the relationship between fructooligosaccharide (FOS) intake at different life stages of Wistar rats and its stimulatory effects on intestinal parameters. METHODS: Recently weaned and ageing female rats were divided into growing and ageing treatments, which were fed diets that partially replaced sucrose with FOS for 12 weeks. RESULTS: Dietary FOS intake induced a significant increase in the numbers of Bifidobacterium and Lactobacillus in growing rats. FOS intake was associated with increased butyric acid levels and a reduced pH of the caecal contents at both ages. Differential gene expression patterns were observed by microarray analysis of growing and ageing animals fed the FOS diet. A total of 133 genes showed detectable changes in expression in the growing rats, while there were only 19 gene expression changes in ageing rats fed with FOS. CONCLUSION: These results suggest that dietary FOS intake may be beneficial for some parameters of intestinal health in growing rats.

Severe DCM phenotype of patient harboring RBM20 mutation S635A can be modeled by patient-specific induced pluripotent stem cell-derived cardiomyocytes.

The ability to generate patient-specific induced pluripotent stem cells (iPSCs) provides a unique opportunity for modeling heart disease in vitro. In this study, we generated iPSCs from a patient with dilated cardiomyopathy (DCM) caused by a missense mutation S635A in RNA-binding motif protein 20 (RBM20) and investigated the functionality and cell biology of cardiomyocytes (CMs) derived from patient-specific iPSCs (RBM20-iPSCs). The RBM20-iPSC-CMs showed abnormal distribution of sarcomeric α-actinin and defective calcium handling compared to control-iPSC-CMs, suggesting disorganized myofilament structure and altered calcium machinery in CMs of the RBM20 patient. Engineered heart muscles (EHMs) from RBM20-iPSC-CMs showed that not only active force generation was impaired in RBM20-EHMs but also passive stress of the tissue was decreased, suggesting a higher visco-elasticity of RBM20-EHMs. Furthermore, we observed a reduced titin (TTN) N2B-isoform expression in RBM20-iPSC-CMs by demonstrating a reduction of exon skipping in the PEVK region of TTN and an inhibition of TTN isoform switch. In contrast, in control-iPSC-CMs both TTN isoforms N2B and N2BA were expressed, indicating that the TTN isoform switch occurs already during early cardiogenesis. Using next generation RNA sequencing, we mapped transcriptome and splicing target profiles of RBM20-iPSC-CMs and identified different cardiac gene networks in response to the analyzed RBM20 mutation in cardiac-specific processes. These findings shed the first light on molecular mechanisms of RBM20-dependent pathological cardiac remodeling leading to DCM. Our data demonstrate that iPSC-CMs coupled with EHMs provide a powerful tool for evaluating disease-relevant functional defects and for a deeper mechanistic understanding of alternative splicing-related cardiac diseases.

Identification of candidate genes for drought tolerance in coffee by high-throughput sequencing in the shoot apex of different Coffea arabica cultivars.

BACKGROUND: Drought is a widespread limiting factor in coffee plants. It affects plant development, fruit production, bean development and consequently beverage quality. Genetic diversity for drought tolerance exists within the coffee genus. However, the molecular mechanisms underlying the adaptation of coffee plants to drought are largely unknown. In this study, we compared the molecular responses to drought in two commercial cultivars (IAPAR59, drought-tolerant and Rubi, drought-susceptible) of Coffea arabica grown in the field under control (irrigation) and drought conditions using the pyrosequencing of RNA extracted from shoot apices and analysing the expression of 38 candidate genes. RESULTS: Pyrosequencing from shoot apices generated a total of 34.7 Mbp and 535,544 reads enabling the identification of 43,087 clusters (41,512 contigs and 1,575 singletons). These data included 17,719 clusters (16,238 contigs and 1,575 singletons) exclusively from 454 sequencing reads, along with 25,368 hybrid clusters assembled with 454 sequences. The comparison of DNA libraries identified new candidate genes (n = 20) presenting differential expression between IAPAR59 and Rubi and/or drought conditions. Their expression was monitored in plagiotropic buds, together with those of other (n = 18) candidates genes. Under drought conditions, up-regulated expression was observed in IAPAR59 but not in Rubi for CaSTK1 (protein kinase), CaSAMT1 (SAM-dependent methyltransferase), CaSLP1 (plant development) and CaMAS1 (ABA biosynthesis). Interestingly, the expression of lipid-transfer protein (nsLTP) genes was also highly up-regulated under drought conditions in IAPAR59. This may have been related to the thicker cuticle observed on the abaxial leaf surface in IAPAR59 compared to Rubi. CONCLUSIONS: The full transcriptome assembly of C. arabica, followed by functional annotation, enabled us to identify differentially expressed genes related to drought conditions. Using these data, candidate genes were selected and their differential expression profiles were confirmed by qPCR experiments in plagiotropic buds of IAPAR59 and Rubi under drought conditions. As regards the genes up-regulated under drought conditions, specifically in the drought-tolerant IAPAR59, several corresponded to orphan genes but also to genes coding proteins involved in signal transduction pathways, as well as ABA and lipid metabolism, for example. The identification of these genes should help advance our understanding of the genetic determinism of drought tolerance in coffee.

Transcription factors involved in prostate gland adaptation to androgen deprivation.

Androgens regulate prostate physiology, and exert their effects through the androgen receptor. We hypothesized that androgen deprivation needs additional transcription factors to orchestrate the changes taking place in the gland after castration and for the adaptation of the epithelial cells to the androgen-deprived environment, ultimately contributing to the origin of castration-resistant prostate cancer. This study was undertaken to identify transcription factors that regulate gene expression after androgen deprivation by castration (Cas). For the sake of comparison, we extended the analysis to the effects of administration of a high dose of 17ß-estradiol (E2) and a combination of both (Cas+E2). We approached this by (i) identifying gene expression profiles and enrichment terms, and by searching for transcription factors in the derived regulatory pathways; and (ii) by determining the density of putative transcription factor binding sites in the proximal promoter of the 10 most up- or down-regulated genes in each experimental group in comparison to the controls Gapdh and Tbp7. Filtering and validation confirmed the expression and localized EVI1 (Mecom), NFY, ELK1, GATA2, MYBL1, MYBL2, and NFkB family members (NFkB1, NFkB2, REL, RELA and RELB) in the epithelial and/or stromal cells. These transcription factors represent major regulators of epithelial cell survival and immaturity as well as an adaptation of the gland as an immune barrier in the absence of functional stimulation by androgens. Elk1 was expressed in smooth muscle cells and was up-regulated after day 4. Evi1 and Nfy genes are expressed in both epithelium and stroma, but were apparently not affected by androgen deprivation.

IIS--Integrated Interactome System: a web-based platform for the annotation, analysis and visualization of protein-metabolite-gene-drug interactions by integrating a variety of data sources and tools.

BACKGROUND: High-throughput screening of physical, genetic and chemical-genetic interactions brings important perspectives in the Systems Biology field, as the analysis of these interactions provides new insights into protein/gene function, cellular metabolic variations and the validation of therapeutic targets and drug design. However, such analysis depends on a pipeline connecting different tools that can automatically integrate data from diverse sources and result in a more comprehensive dataset that can be properly interpreted. RESULTS: We describe here the Integrated Interactome System (IIS), an integrative platform with a web-based interface for the annotation, analysis and visualization of the interaction profiles of proteins/genes, metabolites and drugs of interest. IIS works in four connected modules: (i) Submission module, which receives raw data derived from Sanger sequencing (e.g. two-hybrid system); (ii) Search module, which enables the user to search for the processed reads to be assembled into contigs/singlets, or for lists of proteins/genes, metabolites and drugs of interest, and add them to the project; (iii) Annotation module, which assigns annotations from several databases for the contigs/singlets or lists of proteins/genes, generating tables with automatic annotation that can be manually curated; and (iv) Interactome module, which maps the contigs/singlets or the uploaded lists to entries in our integrated database, building networks that gather novel identified interactions, protein and metabolite expression/concentration levels, subcellular localization and computed topological metrics, GO biological processes and KEGG pathways enrichment. This module generates a XGMML file that can be imported into Cytoscape or be visualized directly on the web. CONCLUSIONS: We have developed IIS by the integration of diverse databases following the need of appropriate tools for a systematic analysis of physical, genetic and chemical-genetic interactions. IIS was validated with yeast two-hybrid, proteomics and metabolomics datasets, but it is also extendable to other datasets. IIS is freely available online at: http://www.lge.ibi.unicamp.br/lnbio/IIS/.

The genome sequence of Leishmania (Leishmania) amazonensis: functional annotation and extended analysis of gene models.

We present the sequencing and annotation of the Leishmania (Leishmania) amazonensis genome, an etiological agent of human cutaneous leishmaniasis in the Amazon region of Brazil. L. (L.) amazonensis shares features with Leishmania (L.) mexicana but also exhibits unique characteristics regarding geographical distribution and clinical manifestations of cutaneous lesions (e.g. borderline disseminated cutaneous leishmaniasis). Predicted genes were scored for orthologous gene families and conserved domains in comparison with other human pathogenic Leishmania spp. Carboxypeptidase, aminotransferase, and 3'-nucleotidase genes and ATPase, thioredoxin, and chaperone-related domains were represented more abundantly in L. (L.) amazonensis and L. (L.) mexicana species. Phylogenetic analysis revealed that these two species share groups of amastin surface proteins unique to the genus that could be related to specific features of disease outcomes and host cell interactions. Additionally, we describe a hypothetical hybrid interactome of potentially secreted L. (L.) amazonensis proteins and host proteins under the assumption that parasite factors mimic their mammalian counterparts. The model predicts an interaction between an L. (L.) amazonensis heat-shock protein and mammalian Toll-like receptor 9, which is implicated in important immune responses such as cytokine and nitric oxide production. The analysis presented here represents valuable information for future studies of leishmaniasis pathogenicity and treatment.

Increasing the density of markers around a major QTL controlling resistance to angular leaf spot in common bean.

Angular leaf spot (ALS) causes major yield losses in the common bean (Phaseolus vulgaris L.), an important protein source in the human diet. This study describes the saturation around a major quantitative trait locus (QTL) region, ALS10.1, controlling resistance to ALS located on linkage group Pv10 and explores the genomic context of this region using available data from the P. vulgaris genome sequence. DArT-derived markers (STS-DArT) selected by bulk segregant analysis and SCAR and SSR markers were used to increase the resolution of the QTL, reducing the confidence interval of ALS10.1 from 13.4 to 3.0 cM. The position of the SSR ATA220 coincided with the maximum LOD score of the QTL. Moreover, a new QTL (ALS10.2(UC)) was identified at the end of the same linkage group. Sequence analysis using the P. vulgaris genome located ten SSRs and seven STS-DArT on chromosome 10 (Pv10). Coincident linkage and genome positions of five markers enabled the definition of a core region for ALS10.1 spanning 5.3 Mb. These markers are linked to putative genes related to disease resistance such as glycosyl transferase, ankyrin repeat-containing, phospholipase, and squamosa-promoter binding protein. Synteny analysis between ALS10.1 markers and the genome of soybean suggested a dynamic evolution of this locus in the common bean. The present study resulted in the identification of new candidate genes and markers closely linked to a major ALS disease resistance QTL, which can be used in marker-assisted selection, fine mapping and positional QTL cloning.

Identification of SNPs in RNA-seq data of two cultivars of Glycine max (soybean) differing in drought resistance.

The legume Glycine max (soybean) plays an important economic role in the international commodities market, with a world production of almost 260 million tons for the 2009/2010 harvest. The increase in drought events in the last decade has caused production losses in recent harvests. This fact compels us to understand the drought tolerance mechanisms in soybean, taking into account its variability among commercial and developing cultivars. In order to identify single nucleotide polymorphisms (SNPs) in genes up-regulated during drought stress, we evaluated suppression subtractive libraries (SSH) from two contrasting cultivars upon water deprivation: sensitive (BR 16) and tolerant (Embrapa 48). A total of 2,222 soybean genes were up-regulated in both cultivars. Our method identified more than 6,000 SNPs in tolerant and sensitive Brazilian cultivars in those drought stress related genes. Among these SNPs, 165 (in 127 genes) are positioned at soybean chromosome ends, including transcription factors (MYB, WRKY) related to tolerance to abiotic stress.

A high-throughput data mining of single nucleotide polymorphisms in Coffea species expressed sequence tags suggests differential homeologous gene expression in the allotetraploid Coffea arabica.

Polyploidization constitutes a common mode of evolution in flowering plants. This event provides the raw material for the divergence of function in homeologous genes, leading to phenotypic novelty that can contribute to the success of polyploids in nature or their selection for use in agriculture. Mounting evidence underlined the existence of homeologous expression biases in polyploid genomes; however, strategies to analyze such transcriptome regulation remained scarce. Important factors regarding homeologous expression biases remain to be explored, such as whether this phenomenon influences specific genes, how paralogs are affected by genome doubling, and what is the importance of the variability of homeologous expression bias to genotype differences. This study reports the expressed sequence tag assembly of the allopolyploid Coffea arabica and one of its direct ancestors, Coffea canephora. The assembly was used for the discovery of single nucleotide polymorphisms through the identification of high-quality discrepancies in overlapped expressed sequence tags and for gene expression information indirectly estimated by the transcript redundancy. Sequence diversity profiles were evaluated within C. arabica (Ca) and C. canephora (Cc) and used to deduce the transcript contribution of the Coffea eugenioides (Ce) ancestor. The assignment of the C. arabica haplotypes to the C. canephora (CaCc) or C. eugenioides (CaCe) ancestral genomes allowed us to analyze gene expression contributions of each subgenome in C. arabica. In silico data were validated by the quantitative polymerase chain reaction and allele-specific combination TaqMAMA-based method. The presence of differential expression of C. arabica homeologous genes and its implications in coffee gene expression, ontology, and physiology are discussed.