ZEB2 haploinsufficient Mowat-Wilson syndrome induced pluripotent stem cells show disrupted GABAergic transcriptional regulation and function.

Mowat-Wilson syndrome (MWS) is a severe neurodevelopmental disorder caused by heterozygous variants in the gene encoding transcription factor ZEB2. Affected individuals present with structural brain abnormalities, speech delay and epilepsy. In mice, conditional loss of Zeb2 causes hippocampal degeneration, altered migration and differentiation of GABAergic interneurons, a heterogeneous population of mainly inhibitory neurons of importance for maintaining normal excitability. To get insights into GABAergic development and function in MWS we investigated ZEB2 haploinsufficient induced pluripotent stem cells (iPSC) of MWS subjects together with iPSC of healthy donors. Analysis of RNA-sequencing data at two time points of GABAergic development revealed an attenuated interneuronal identity in MWS subject derived iPSC with enrichment of differentially expressed genes required for transcriptional regulation, cell fate transition and forebrain patterning. The ZEB2 haploinsufficient neural stem cells (NSCs) showed downregulation of genes required for ventral telencephalon specification, such as FOXG1, accompanied by an impaired migratory capacity. Further differentiation into GABAergic interneuronal cells uncovered upregulation of transcription factors promoting pallial and excitatory neurons whereas cortical markers were downregulated. The differentially expressed genes formed a neural protein-protein network with extensive connections to well-established epilepsy genes. Analysis of electrophysiological properties in ZEB2 haploinsufficient GABAergic cells revealed overt perturbations manifested as impaired firing of repeated action potentials. Our iPSC model of ZEB2 haploinsufficient GABAergic development thus uncovers a dysregulated gene network leading to immature interneurons with mixed identity and altered electrophysiological properties, suggesting mechanisms contributing to the neuropathogenesis and seizures in MWS.

Integration of molecular profiles in a longitudinal wellness profiling cohort.

An important aspect of precision medicine is to probe the stability in molecular profiles among healthy individuals over time. Here, we sample a longitudinal wellness cohort with 100 healthy individuals and analyze blood molecular profiles including proteomics, transcriptomics, lipidomics, metabolomics, autoantibodies and immune cell profiling, complemented with gut microbiota composition and routine clinical chemistry. Overall, our results show high variation between individuals across different molecular readouts, while the intra-individual baseline variation is low. The analyses show that each individual has a unique and stable plasma protein profile throughout the study period and that many individuals also show distinct profiles with regards to the other omics datasets, with strong underlying connections between the blood proteome and the clinical chemistry parameters. In conclusion, the results support an individual-based definition of health and show that comprehensive omics profiling in a longitudinal manner is a path forward for precision medicine.

Transcriptomes of Dravet syndrome iPSC derived GABAergic cells reveal dysregulated pathways for chromatin remodeling and neurodevelopment.

Dravet syndrome (DS) is an early onset refractory epilepsy typically caused by de novo heterozygous variants in SCN1A encoding the α-subunit of the neuronal sodium channel Nav1.1. The syndrome is characterized by age-related progression of seizures, cognitive decline and movement disorders. We hypothesized that the distinct neurodevelopmental features in DS are caused by the disruption of molecular pathways in Nav1.1 haploinsufficient cells resulting in perturbed neural differentiation and maturation. Here, we established DS-patient and control induced pluripotent stem cell derived neural progenitor cells (iPSC NPC) and GABAergic inter-neuronal (iPSC GABA) cells. The DS-patient iPSC GABA cells showed a shift in sodium current activation and a perturbed response to induced oxidative stress. Transcriptome analysis revealed specific dysregulations of genes for chromatin structure, mitotic progression, neural plasticity and excitability in DS-patient iPSC NPCs and DS-patient iPSC GABA cells versus controls. The transcription factors FOXM1 and E2F1, positive regulators of the disrupted pathways for histone modification and cell cycle regulation, were markedly up-regulated in DS-iPSC GABA lines. Our study highlights transcriptional changes and disrupted pathways of chromatin remodeling in Nav1.1 haploinsufficient GABAergic cells, providing a molecular framework that overlaps with that of neurodevelopmental disorders and other epilepsies.

Identification of Salivary Microbiota and Its Association With Host Inflammatory Mediators in Periodontitis.

Periodontitis is a microbial-induced chronic inflammatory disease, which may not only result in tooth loss, but can also contribute to the development of various systemic diseases. The transition from healthy to diseased periodontium depends on microbial dysbiosis and impaired host immune response. Although periodontitis is a common disease as well as associated with various systemic inflammatory conditions, the taxonomic profiling of the salivary microbiota in periodontitis and its association with host immune and inflammatory mediators has not been reported. Therefore, the aim of this study was to identify key pathogens and their potential interaction with the host's inflammatory mediators in saliva samples for periodontitis risk assessment. The microbial 16S rRNA gene sequencing and the levels of inflammatory mediators were performed in saliva samples from patients with chronic periodontitis and periodontally healthy control subjects. The salivary microbial community composition differed significantly between patients with chronic periodontitis and healthy controls. Our analyses identified a number of microbes, including bacteria assigned to Eubacterium saphenum, Tannerella forsythia, Filifactor alocis, Streptococcus mitis/parasanguinis, Parvimonas micra, Prevotella sp., Phocaeicola sp., and Fretibacterium sp. as more abundant in periodontitis, compared to healthy controls. In samples from healthy individuals, we identified Campylobacter concisus, and Veillonella sp. as more abundant. Integrative analysis of the microbiota and inflammatory mediators/cytokines revealed associations that included positive correlations between the pathogens Treponema sp. and Selenomas sp. and the cytokines chitinase 3-like 1, sIL-6Rα, sTNF-R1, and gp130/sIL-6Rß. In addition, a negative correlation was identified between IL-10 and Filifactor alocis. Our results reveal distinct and disease-specific patterns of salivary microbial composition between patients with periodontitis and healthy controls, as well as significant correlations between microbiota and host-mediated inflammatory cytokines. The positive correlations between the pathogens Treponema sp. and Selenomas sp. and the cytokines chitinase 3-like 1, sIL-6Rα, sTNF-R1, and gp130/sIL-6Rß might have the future potential to serve as a combined bacteria-host salivary biomarker panel for diagnosis of the chronic infectious disease periodontitis. However, further studies are required to determine the capacity of these microbes and inflammatory mediators as a salivary biomarker panel for periodontitis.

Breast cancer quantitative proteome and proteogenomic landscape.

In the preceding decades, molecular characterization has revolutionized breast cancer (BC) research and therapeutic approaches. Presented herein, an unbiased analysis of breast tumor proteomes, inclusive of 9995 proteins quantified across all tumors, for the first time recapitulates BC subtypes. Additionally, poor-prognosis basal-like and luminal B tumors are further subdivided by immune component infiltration, suggesting the current classification is incomplete. Proteome-based networks distinguish functional protein modules for breast tumor groups, with co-expression of EGFR and MET marking ductal carcinoma in situ regions of normal-like tumors and lending to a more accurate classification of this poorly defined subtype. Genes included within prognostic mRNA panels have significantly higher than average mRNA-protein correlations, and gene copy number alterations are dampened at the protein-level; underscoring the value of proteome quantification for prognostication and phenotypic classification. Furthermore, protein products mapping to non-coding genomic regions are identified; highlighting a potential new class of tumor-specific immunotherapeutic targets.

Transcriptome and Proteome Profiling of Neural Induced Pluripotent Stem Cells from Individuals with Down Syndrome Disclose Dynamic Dysregulations of Key Pathways and Cellular Functions.

Down syndrome (DS) or trisomy 21 (T21) is a leading genetic cause of intellectual disability. To gain insights into dynamics of molecular perturbations during neurogenesis in DS, we established a model using induced pluripotent stem cells (iPSC) with transcriptome profiles comparable to that of normal fetal brain development. When applied on iPSCs with T21, transcriptome and proteome signatures at two stages of differentiation revealed strong temporal dynamics of dysregulated genes, proteins and pathways belonging to 11 major functional clusters. DNA replication, synaptic maturation and neuroactive clusters were disturbed at the early differentiation time point accompanied by a skewed transition from the neural progenitor cell stage and reduced cellular growth. With differentiation, growth factor and extracellular matrix, oxidative phosphorylation and glycolysis emerged as major perturbed clusters. Furthermore, we identified a marked dysregulation of a set of genes encoded by chromosome 21 including an early upregulation of the hub gene APP, supporting its role for disturbed neurogenesis, and the transcription factors OLIG1, OLIG2 and RUNX1, consistent with deficient myelination and neuronal differentiation. Taken together, our findings highlight novel sequential and differentiation-dependent dynamics of disturbed functions, pathways and elements in T21 neurogenesis, providing further insights into developmental abnormalities of the DS brain.

A primer on deep learning in genomics.

Deep learning methods are a class of machine learning techniques capable of identifying highly complex patterns in large datasets. Here, we provide a perspective and primer on deep learning applications for genome analysis. We discuss successful applications in the fields of regulatory genomics, variant calling and pathogenicity scores. We include general guidance for how to effectively use deep learning methods as well as a practical guide to tools and resources. This primer is accompanied by an interactive online tutorial.

Stationary and portable sequencing-based approaches for tracing wastewater contamination in urban stormwater systems.

Urban sewer systems consist of wastewater and stormwater sewers, of which only wastewater is processed before being discharged. Occasionally, misconnections or damages in the network occur, resulting in untreated wastewater entering natural water bodies via the stormwater system. Cultivation of faecal indicator bacteria (e.g. Escherichia coli; E. coli) is the current standard for tracing wastewater contamination. This method is cheap but has limited specificity and mobility. Here, we compared the E. coli culturing approach with two sequencing-based methodologies (Illumina MiSeq 16S rRNA gene amplicon sequencing and Oxford Nanopore MinION shotgun metagenomic sequencing), analysing 73 stormwater samples collected in Stockholm. High correlations were obtained between E. coli culturing counts and frequencies of human gut microbiome amplicon sequences, indicating E. coli is indeed a good indicator of faecal contamination. However, the amplicon data further holds information on contamination source or alternatively how much time has elapsed since the faecal matter has entered the system. Shotgun metagenomic sequencing on a subset of the samples using a portable real-time sequencer, MinION, correlated well with the amplicon sequencing data. This study demonstrates the use of DNA sequencing to detect human faecal contamination in stormwater systems and the potential of tracing faecal contamination directly in the field.

Publisher Correction: Discovery of coding regions in the human genome by integrated proteogenomics analysis workflow.

In the original version of this Article, extraneous text not belonging to the Article was accidentally appended to the results section. This error has now been corrected in both the PDF and HTML versions of the Article.

Transcriptome profiling of the interconnection of pathways involved in malignant transformation and response to hypoxia.

In tumor tissues, hypoxia is a commonly observed feature resulting from rapidly proliferating cancer cells outgrowing their surrounding vasculature network. Transformed cancer cells are known to exhibit phenotypic alterations, enabling continuous proliferation despite a limited oxygen supply. The four-step isogenic BJ cell model enables studies of defined steps of tumorigenesis: the normal, immortalized, transformed, and metastasizing stages. By transcriptome profiling under atmospheric and moderate hypoxic (3% O2) conditions, we observed that despite being highly similar, the four cell lines of the BJ model responded strikingly different to hypoxia. Besides corroborating many of the known responses to hypoxia, we demonstrate that the transcriptome adaptation to moderate hypoxia resembles the process of malignant transformation. The transformed cells displayed a distinct capability of metabolic switching, reflected in reversed gene expression patterns for several genes involved in oxidative phosphorylation and glycolytic pathways. By profiling the stage-specific responses to hypoxia, we identified ASS1 as a potential prognostic marker in hypoxic tumors. This study demonstrates the usefulness of the BJ cell model for highlighting the interconnection of pathways involved in malignant transformation and hypoxic response.

Discovery of coding regions in the human genome by integrated proteogenomics analysis workflow.

Proteogenomics enable the discovery of novel peptides (from unannotated genomic protein-coding loci) and single amino acid variant peptides (derived from single-nucleotide polymorphisms and mutations). Increasing the reliability of these identifications is crucial to ensure their usefulness for genome annotation and potential application as neoantigens in cancer immunotherapy. We here present integrated proteogenomics analysis workflow (IPAW), which combines peptide discovery, curation, and validation. IPAW includes the SpectrumAI tool for automated inspection of MS/MS spectra, eliminating false identifications of single-residue substitution peptides. We employ IPAW to analyze two proteomics data sets acquired from A431 cells and five normal human tissues using extended (pH range, 3-10) high-resolution isoelectric focusing (HiRIEF) pre-fractionation and TMT-based peptide quantitation. The IPAW results provide evidence for the translation of pseudogenes, lncRNAs, short ORFs, alternative ORFs, N-terminal extensions, and intronic sequences. Moreover, our quantitative analysis indicates that protein production from certain pseudogenes and lncRNAs is tissue specific.

Impact of genetic risk loci for multiple sclerosis on expression of proximal genes in patients.

Despite advancements in genetic studies, it is difficult to understand and characterize the functional relevance of disease-associated genetic variants, especially in the context of a complex multifactorial disease such as multiple sclerosis (MS). As a large proportion of expression quantitative trait loci (eQTLs) are context-specific, we performed RNA-Seq in peripheral blood mononuclear cells from MS patients (n = 145) to identify eQTLs in regions centered on 109 MS risk single nucleotide polymorphisms and 7 associated human leukocyte antigen variants. We identified 77 statistically significant eQTL associations, including pseudogenes and non-coding RNAs. Thirty-eight out of 40 testable eQTL effects were colocalized with the disease association signal. As many eQTLs are tissue specific, we aimed to detail their significance in different cell types. Approximately 70% of the eQTLs were replicated and characterized in at least one major peripheral blood mononuclear cell-derived cell type. Furthermore, 40% of eQTLs were found to be more pronounced in MS patients compared with non-inflammatory neurological diseases patients. In addition, we found two single nucleotide polymorphisms to be significantly associated with the proportions of three different cell types. Mapping to enhancer histone marks and predicted transcription factor binding sites added additional functional evidence for eight eQTL regions. As an example, we found that rs71624119, shared with three other autoimmune diseases and located in a primed enhancer (H3K4me1) with potential binding for STAT transcription factors, significantly associates with ANKRD55 expression. This study provides many novel and validated targets for future functional characterization of MS and other diseases.

The role of Cdx2 as a lineage specific transcriptional repressor for pluripotent network during the first developmental cell lineage segregation.

The first cellular differentiation event in mouse development leads to the formation of the blastocyst consisting of the inner cell mass (ICM) and trophectoderm (TE). The transcription factor CDX2 is required for proper TE specification, where it promotes expression of TE genes, and represses expression of Pou5f1 (OCT4). However its downstream network in the developing embryo is not fully characterized. Here, we performed high-throughput single embryo qPCR analysis in Cdx2 null embryos to identify CDX2-regulated targets in vivo. To identify genes likely to be regulated by CDX2 directly, we performed CDX2 ChIP-Seq on trophoblast stem (TS) cells. In addition, we examined the dynamics of gene expression changes using inducible CDX2 embryonic stem (ES) cells, so that we could predict which CDX2-bound genes are activated or repressed by CDX2 binding. By integrating these data with observations of chromatin modifications, we identify putative novel regulatory elements that repress gene expression in a lineage-specific manner. Interestingly, we found CDX2 binding sites within regulatory elements of key pluripotent genes such as Pou5f1 and Nanog, pointing to the existence of a novel mechanism by which CDX2 maintains repression of OCT4 in trophoblast. Our study proposes a general mechanism in regulating lineage segregation during mammalian development.

Erratum to: Identification of mutations, gene expression changes and fusion transcripts by whole transcriptome RNAseq in docetaxel resistant prostate cancer cells.

[This corrects the article DOI: 10.1186/s40064-016-3543-0.].

Identification of mutations, gene expression changes and fusion transcripts by whole transcriptome RNAseq in docetaxel resistant prostate cancer cells.

Docetaxel has been the standard first-line therapy in metastatic castration resistant prostate cancer. The survival benefit is, however, limited by either primary or acquired resistance. In this study, Du145 prostate cancer cells were converted to docetaxel-resistant cells Du145-R and Du145-RB by in vitro culturing. Next generation RNAseq was employed to analyze these cell lines. Forty-two genes were identified to have acquired mutations after the resistance development, of which thirty-four were found to have mutations in published sequencing studies using prostate cancer samples from patients. Fourteen novel and 2 previously known fusion genes were inferred from the RNA-seq data, and 13 of these were validated by RT-PCR and/or re-sequencing. Four in-frame fusion transcripts could be transcribed into fusion proteins in stably transfected HEK293 cells, including MYH9-EIF3D and LDLR-RPL31P11, which were specific identified or up-regulated in the docetaxel resistant DU145 cells. A panel of 615 gene transcripts was identified to have significantly changed expression profile in the docetaxel resistant cells. These transcriptional changes have potential for further study as predictive biomarkers and as targets of docetaxel treatment.

The Impact of Endurance Training on Human Skeletal Muscle Memory, Global Isoform Expression and Novel Transcripts.

Regularly performed endurance training has many beneficial effects on health and skeletal muscle function, and can be used to prevent and treat common diseases e.g. cardiovascular disease, type II diabetes and obesity. The molecular adaptation mechanisms regulating these effects are incompletely understood. To date, global transcriptome changes in skeletal muscles have been studied at the gene level only. Therefore, global isoform expression changes following exercise training in humans are unknown. Also, the effects of repeated interventions on transcriptional memory or training response have not been studied before. In this study, 23 individuals trained one leg for three months. Nine months later, 12 of the same subjects trained both legs in a second training period. Skeletal muscle biopsies were obtained from both legs before and after both training periods. RNA sequencing analysis of all 119 skeletal muscle biopsies showed that training altered the expression of 3,404 gene isoforms, mainly associated with oxidative ATP production. Fifty-four genes had isoforms that changed in opposite directions. Training altered expression of 34 novel transcripts, all with protein-coding potential. After nine months of detraining, no training-induced transcriptome differences were detected between the previously trained and untrained legs. Although there were several differences in the physiological and transcriptional responses to repeated training, no coherent evidence of an endurance training induced transcriptional skeletal muscle memory was found. This human lifestyle intervention induced differential expression of thousands of isoforms and several transcripts from unannotated regions of the genome. It is likely that the observed isoform expression changes reflect adaptational mechanisms and processes that provide the functional and health benefits of regular physical activity.

Visualization and analysis of gene expression in tissue sections by spatial transcriptomics.

Analysis of the pattern of proteins or messengerRNAs (mRNAs) in histological tissue sections is a cornerstone in biomedical research and diagnostics. This typically involves the visualization of a few proteins or expressed genes at a time. We have devised a strategy, which we call "spatial transcriptomics," that allows visualization and quantitative analysis of the transcriptome with spatial resolution in individual tissue sections. By positioning histological sections on arrayed reverse transcription primers with unique positional barcodes, we demonstrate high-quality RNA-sequencing data with maintained two-dimensional positional information from the mouse brain and human breast cancer. Spatial transcriptomics provides quantitative gene expression data and visualization of the distribution of mRNAs within tissue sections and enables novel types of bioinformatics analyses, valuable in research and diagnostics.

Cryptosporidium as a testbed for single cell genome characterization of unicellular eukaryotes.

BACKGROUND: Infectious disease involving multiple genetically distinct populations of pathogens is frequently concurrent, but difficult to detect or describe with current routine methodology. Cryptosporidium sp. is a widespread gastrointestinal protozoan of global significance in both animals and humans. It cannot be easily maintained in culture and infections of multiple strains have been reported. To explore the potential use of single cell genomics methodology for revealing genome-level variation in clinical samples from Cryptosporidium-infected hosts, we sorted individual oocysts for subsequent genome amplification and full-genome sequencing. RESULTS: Cells were identified with fluorescent antibodies with an 80 % success rate for the entire single cell genomics workflow, demonstrating that the methodology can be applied directly to purified fecal samples. Ten amplified genomes from sorted single cells were selected for genome sequencing and compared both to the original population and a reference genome in order to evaluate the accuracy and performance of the method. Single cell genome coverage was on average 81 % even with a moderate sequencing effort and by combining the 10 single cell genomes, the full genome was accounted for. By a comparison to the original sample, biological variation could be distinguished and separated from noise introduced in the amplification. CONCLUSIONS: As a proof of principle, we have demonstrated the power of applying single cell genomics to dissect infectious disease caused by closely related parasite species or subtypes. The workflow can easily be expanded and adapted to target other protozoans, and potential applications include mapping genome-encoded traits, virulence, pathogenicity, host specificity and resistance at the level of cells as truly meaningful biological units.

Distinct transcription factor complexes act on a permissive chromatin landscape to establish regionalized gene expression in CNS stem cells.

Spatially distinct gene expression profiles in neural stem cells (NSCs) are a prerequisite to the formation of neuronal diversity, but how these arise from the regulatory interactions between chromatin accessibility and transcription factor activity has remained unclear. Here, we demonstrate that, despite their distinct gene expression profiles, NSCs of the mouse cortex and spinal cord share the majority of their DNase I hypersensitive sites (DHSs). Regardless of this similarity, domain-specific gene expression is highly correlated with the relative accessibility of associated DHSs, as determined by sequence read density. Notably, the binding pattern of the general NSC transcription factor SOX2 is also largely cell type specific and coincides with an enrichment of LHX2 motifs in the cortex and HOXA9 motifs in the spinal cord. Interestingly, in a zebrafish reporter gene system, these motifs were critical determinants of patterned gene expression along the rostral-caudal axis. Our findings establish a predictive model for patterned NSC gene expression, whereby domain-specific expression of LHX2 and HOX proteins act on their target motifs within commonly accessible cis-regulatory regions to specify SOX2 binding. In turn, this binding correlates strongly with these DHSs relative accessibility-a robust predictor of neighboring gene expression.

Evolutionary history of host use, rather than plant phylogeny, determines gene expression in a generalist butterfly.

BACKGROUND: Although most insect species are specialized on one or few groups of plants, there are phytophagous insects that seem to use virtually any kind of plant as food. Understanding the nature of this ability to feed on a wide repertoire of plants is crucial for the control of pest species and for the elucidation of the macroevolutionary mechanisms of speciation and diversification of insect herbivores. Here we studied Vanessa cardui, the species with the widest diet breadth among butterflies and a potential insect pest, by comparing tissue-specific transcriptomes from caterpillars that were reared on different host plants. We tested whether the similarities of gene-expression response reflect the evolutionary history of adaptation to these plants in the Vanessa and related genera, against the null hypothesis of transcriptional profiles reflecting plant phylogenetic relatedness. RESULT: Using both unsupervised and supervised methods of data analysis, we found that the tissue-specific patterns of caterpillar gene expression are better explained by the evolutionary history of adaptation of the insects to the plants than by plant phylogeny. CONCLUSION: Our findings suggest that V. cardui may use two sets of expressed genes to achieve polyphagy, one associated with the ancestral capability to consume Rosids and Asterids, and another allowing the caterpillar to incorporate a wide range of novel host-plants.