Mingyan Lin

RNA-Seq of human neurons derived from iPS cells reveals candidate long non-coding RNAs involved in neurogenesis and neuropsychiatric disorders.

Genome-wide expression analysis using next generation sequencing (RNA-Seq) provides an opportunity for in-depth molecular profiling of fundamental biological processes, such as cellular differentiation and malignant transformation. Differentiating human neurons derived from induced pluripotent stem cells (iPSCs) provide an ideal system for RNA-Seq since defective neurogenesis caused by abnormalities in transcription factors, DNA methylation, and chromatin modifiers lie at the heart of some neuropsychiatric disorders. As a preliminary step towards applying next generation sequencing using neurons derived from patient-specific iPSCs, we have carried out an RNA-Seq analysis on control human neurons. Dramatic changes in the expression of coding genes, long non-coding RNAs (lncRNAs), pseudogenes, and splice isoforms were seen during the transition from pluripotent stem cells to early differentiating neurons. A number of genes that undergo radical changes in expression during this transition include candidates for schizophrenia (SZ), bipolar disorder (BD) and autism spectrum disorders (ASD) that function as transcription factors and chromatin modifiers, such as POU3F2 and ZNF804A, and genes coding for cell adhesion proteins implicated in these conditions including NRXN1 and NLGN1. In addition, a number of novel lncRNAs were found to undergo dramatic changes in expression, one of which is HOTAIRM1, a regulator of several HOXA genes during myelopoiesis. The increase we observed in differentiating neurons suggests a role in neurogenesis as well. Finally, several lncRNAs that map near SNPs associated with SZ in genome wide association studies also increase during neuronal differentiation, suggesting that these novel transcripts may be abnormally regulated in a subgroup of patients.

PRC2 is recurrently inactivated through EED or SUZ12 loss in malignant peripheral nerve sheath tumors

Malignant peripheral nerve sheath tumors (MPNSTs) represent a group of highly aggressive soft-tissue sarcomas that may occur sporadically, in association with neurofibromatosis type I (NF1 associated) or after radiotherapy. Using comprehensive genomic approaches, we identified loss-of-function somatic alterations of the Polycomb repressive complex 2 (PRC2) components (EED or SUZ12) in 92% of sporadic, 70% of NF1-associated and 90% of radiotherapy-associated MPNSTs. MPNSTs with PRC2 loss showed complete loss of trimethylation at lysine 27 of histone H3 (H3K27me3) and aberrant transcriptional activation of multiple PRC2-repressed homeobox master regulators and their regulated developmental pathways. Introduction of the lost PRC2 component in a PRC2-deficient MPNST cell line restored H3K27me3 levels and decreased cell growth. Additionally, we identified frequent somatic alterations of CDKN2A (81% of all MPNSTs) and NF1 (72% of non-NF1-associated MPNSTs), both of which significantly co-occur with PRC2 alterations. The highly recurrent and specific inactivation of PRC2 components, NF1 and CDKN2A highlights their critical and potentially cooperative roles in MPNST pathogenesis.

Identification and Initial Functional Characterization of a Human Vascular Cell–Enriched Long Noncoding RNA

Objective—Long noncoding RNAs (lncRNAs) represent a rapidly growing class of RNA genes with functions related primarily to transcriptional and post-transcriptional control of gene expression. There is a paucity of information about lncRNA expression and function in human vascular cells. Thus, we set out to identify novel lncRNA genes in human vascular smooth muscle cells and to gain insight into their role in the control of smooth muscle cell phenotypes. Approach and Results—RNA sequencing (RNA-seq) of human coronary artery smooth muscle cells revealed 31 unannotated lncRNAs, including a vascular cell–enriched lncRNA (Smooth muscle and Endothelial cell–enriched migration/differentiation-associated long NonCoding RNA [SENCR]). Strand-specific reverse transcription polymerase chain reaction (PCR) and rapid amplification of cDNA ends indicate that SENCR is transcribed antisense from the 5′ end of the FLI1 gene and exists as 2 splice variants. RNA fluorescence in situ hybridization and biochemical fractionation studies demonstrate SENCR is a cytoplasmic lncRNA. Consistent with this observation, knockdown studies reveal little to no cis-acting effect of SENCR on FLI1 or neighboring gene expression. RNA-seq experiments in smooth muscle cells after SENCR knockdown disclose decreased expression of Myocardin and numerous smooth muscle contractile genes, whereas several promigratory genes are increased. Reverse transcription PCR and Western blotting experiments validate several differentially expressed genes after SENCR knockdown. Loss-of-function studies in scratch wound and Boyden chamber assays support SENCR as an inhibitor of smooth muscle cell migration. Conclusions—SENCR is a new vascular cell–enriched, cytoplasmic lncRNA that seems to stabilize the smooth muscle cell contractile phenotype.

In vivo transcriptional governance of hair follicle stem cells by canonical Wnt regulators

Hair follicle stem cells (HFSCs) regenerate hair in response to Wnt signalling. Here, we unfold genome-wide transcriptional and chromatin landscapes of β-catenin–TCF3/4–TLE circuitry, and genetically dissect their biological roles within the native HFSC niche. We show that during HFSC quiescence, TCF3, TCF4 and TLE (Groucho) bind coordinately and transcriptionally repress Wnt target genes. We also show that β-catenin is dispensable for HFSC viability, and if TCF3/4 levels are sufficiently reduced, it is dispensable for proliferation. However, β-catenin is essential to activate genes that launch hair follicle fate and suppress sebocyte fate determination. TCF3/4 deficiency mimics Wnt–β-catenin-dependent activation of these hair follicle fate targets; TCF3 overexpression parallels their TLE4-dependent suppression. Our studies unveil TCF3/4–TLE histone deacetylases as a repressive rheostat, whose action can be relieved by Wnt–β-catenin signalling. When TCF3/4 and TLE levels are high, HFSCs can maintain stemness, but remain quiescent. When these levels drop or when Wnt–β-catenin levels rise, this balance is shifted and hair regeneration initiates.

SOX9: a stem cell transcriptional regulator of secreted niche signaling factors

Hair follicles (HFs) undergo cyclical periods of growth, which are fueled by stem cells (SCs) at the base of the resting follicle. HF-SC formation occurs during HF development and requires transcription factor SOX9. Whether and how SOX9 functions in HF-SC maintenance remain unknown. By conditionally targeting Sox9 in adult HF-SCs, we show that SOX9 is essential for maintaining them. SOX9-deficient HF-SCs still transition from quiescence to proliferation and launch the subsequent hair cycle. However, once activated, bulge HF-SCs begin to differentiate into epidermal cells, which naturally lack SOX9. In addition, as HF-SC numbers dwindle, outer root sheath production is not sustained, and HF downgrowth arrests prematurely. Probing the mechanism, we used RNA sequencing (RNA-seq) to identify SOX9-dependent transcriptional changes and chromatin immunoprecipitation (ChIP) and deep sequencing (ChIP-seq) to identify SOX9-bound genes in HF-SCs. Intriguingly, a large cohort of SOX9-sensitive targets encode extracellular factors, most notably enhancers of Activin/pSMAD2 signaling. Moreover, compromising Activin signaling recapitulates SOX9-dependent defects, and Activin partially rescues them. Overall, our findings reveal roles for SOX9 in regulating adult HF-SC maintenance and suppressing epidermal differentiation in the niche. In addition, our studies expose a role for SCs in coordinating their own behavior in part through non-cell-autonomous signaling within the niche.

Allele-Biased Expression in Differentiating Human Neurons: Implications for Neuropsychiatric Disorders

Stochastic processes and imprinting, along with genetic factors, lead to monoallelic or allele-biased gene expression. Stochastic monoallelic expression fine-tunes information processing in immune cells and the olfactory system, and imprinting plays an important role in development. Recent studies suggest that both stochastic events and imprinting may be more widespread than previously considered. We are interested in allele-biased gene expression occurring in the brain because parent-of-origin effects suggestive of imprinting appear to play a role in the transmission of schizophrenia (SZ) and autism spectrum disorders (ASD) in some families. In addition, allele-biased expression could help explain monozygotic (MZ) twin discordance and reduced penetrance. The ability to study allele-biased expression in human neurons has been transformed with the advent of induced pluripotent stem cell (iPSC) technology and next generation sequencing. Using transcriptome sequencing (RNA-Seq) we identified 801 genes in differentiating neurons that were expressed in an allele-biased manner. These included a number of putative SZ and ASD candidates, such as A2BP1 (RBFOX1), ERBB4, NLGN4X, NRG1, NRG3, NRXN1, and NLGN1. Overall, there was a modest enrichment for SZ and ASD candidate genes among those that showed evidence for allele-biased expression (chi-square, p = 0.02). In addition to helping explain MZ twin discordance and reduced penetrance, the capacity to group many candidate genes affecting a variety of molecular and cellular pathways under a common regulatory process – allele-biased expression – could have therapeutic implications.

CRISPR/Cas9-mediated heterozygous knockout of the autism gene CHD8 and characterization of its transcriptional networks in neurodevelopment

BackgroundDisruptive mutation in the CHD8 gene is one of the top genetic risk factors in autism spectrum disorders (ASDs). Previous analyses of genome-wide CHD8 occupancy and reduced expression of CHD8 by shRNA knockdown in committed neural cells showed that CHD8 regulates multiple cell processes critical for neural functions, and its targets are enriched with ASD-associated genes.MethodsTo further understand the molecular links between CHD8 functions and ASD, we have applied the CRISPR/Cas9 technology to knockout one copy of CHD8 in induced pluripotent stem cells (iPSCs) to better mimic the loss-of-function status that would exist in the developing human embryo prior to neuronal differentiation. We then carried out transcriptomic and bioinformatic analyses of neural progenitors and neurons derived from the CHD8 mutant iPSCs.ResultsTranscriptome profiling revealed that CHD8 hemizygosity (CHD8+/−) affected the expression of several thousands of genes in neural progenitors and early differentiating neurons. The differentially expressed genes were enriched for functions of neural development, β-catenin/Wnt signaling, extracellular matrix, and skeletal system development. They also exhibited significant overlap with genes previously associated with autism and schizophrenia, as well as the downstream transcriptional targets of multiple genes implicated in autism. Providing important insight into how CHD8 mutations might give rise to macrocephaly, we found that seven of the twelve genes associated with human brain volume or head size by genome-wide association studies (e.g., HGMA2) were dysregulated in CHD8+/− neural progenitors or neurons.ConclusionsWe have established a renewable source of CHD8+/− iPSC lines that would be valuable for investigating the molecular and cellular functions of CHD8. Transcriptomic profiling showed that CHD8 regulates multiple genes implicated in ASD pathogenesis and genes associated with brain volume.

MicroRNA Profiling of Neurons Generated Using Induced Pluripotent Stem Cells Derived from Patients with Schizophrenia and Schizoaffective Disorder, and 22q11.2 Del.

We are using induced pluripotent stem cell (iPSC) technology to study neuropsychiatric disorders associated with 22q11.2 microdeletions (del), the most common known schizophrenia (SZ)-associated genetic factor. Several genes in the region have been implicated; a promising candidate is DGCR8, which codes for a protein involved in microRNA (miRNA) biogenesis. We carried out miRNA expression profiling (miRNA-seq) on neurons generated from iPSCs derived from controls and SZ patients with 22q11.2 del. Using thresholds of p<0.01 for nominal significance and 1.5-fold differences in expression, 45 differentially expressed miRNAs were detected (13 lower in SZ and 32 higher). Of these, 6 were significantly down-regulated in patients after correcting for genome wide significance (FDR<0.05), including 4 miRNAs that map to the 22q11.2 del region. In addition, a nominally significant increase in the expression of several miRNAs was found in the 22q11.2 neurons that were previously found to be differentially expressed in autopsy samples and peripheral blood in SZ and autism spectrum disorders (e.g., miR-34, miR-4449, miR-146b-3p, and miR-23a-5p). Pathway and function analysis of predicted mRNA targets of the differentially expressed miRNAs showed enrichment for genes involved in neurological disease and psychological disorders for both up and down regulated miRNAs. Our findings suggest that: i. neurons with 22q11.2 del recapitulate the miRNA expression patterns expected of 22q11.2 haploinsufficiency, ii. differentially expressed miRNAs previously identified using autopsy samples and peripheral cells, both of which have significant methodological problems, are indeed disrupted in neuropsychiatric disorders and likely have an underlying genetic basis.

Transcriptome Comparison of Human Neurons Generated Using Induced Pluripotent Stem Cells Derived from Dental Pulp and Skin Fibroblasts

Induced pluripotent stem cell (iPSC) technology is providing an opportunity to study neuropsychiatric disorders through the capacity to grow patient-specific neurons in vitro. Skin fibroblasts obtained by biopsy have been the most reliable source of cells for reprogramming. However, using other somatic cells obtained by less invasive means would be ideal, especially in children with autism spectrum disorders (ASD) and other neurodevelopmental conditions. In addition to fibroblasts, iPSCs have been developed from cord blood, lymphocytes, hair keratinocytes, and dental pulp from deciduous teeth. Of these, dental pulp would be a good source for neurodevelopmental disorders in children because obtaining material is non-invasive. We investigated its suitability for disease modeling by carrying out gene expression profiling, using RNA-seq, on differentiated neurons derived from iPSCs made from dental pulp extracted from deciduous teeth (T-iPSCs) and fibroblasts (F-iPSCs). This is the first RNA-seq analysis comparing gene expression profiles in neurons derived from iPSCs made from different somatic cells. For the most part, gene expression profiles were quite similar with only 329 genes showing differential expression at a nominally significant p-value (p<0.05), of which 63 remained significant after correcting for genome-wide analysis (FDR <0.05). The most striking difference was the lower level of expression detected for numerous members of the all four HOX gene families in neurons derived from T-iPSCs. In addition, an increased level of expression was seen for several transcription factors expressed in the developing forebrain (FOXP2, OTX1, and LHX2, for example). Overall, pathway analysis revealed that differentially expressed genes that showed higher levels of expression in neurons derived from T-iPSCs were enriched for genes implicated in schizophrenia (SZ). The findings suggest that neurons derived from T-iPSCs are suitable for disease-modeling neuropsychiatric disorder and may have some advantages over those derived from F-iPSCs.

Characterization of human pseudogene-derived non-coding RNAs for functional potential.

Thousands of pseudogenes exist in the human genome and many are transcribed, but their functional potential remains elusive and understudied. To explore these issues systematically, we first developed a computational pipeline to identify transcribed pseudogenes from RNA-Seq data. Applying the pipeline to datasets from 16 distinct normal human tissues identified ∼3,000 pseudogenes that could produce non-coding RNAs in a manner of low abundance but high tissue specificity under normal physiological conditions. Cross-tissue comparison revealed that the transcriptional profiles of pseudogenes and their parent genes showed mostly positive correlations, suggesting that pseudogene transcription could have a positive effect on the expression of their parent genes, perhaps by functioning as competing endogenous RNAs (ceRNAs), as previously suggested and demonstrated with the PTEN pseudogene, PTENP1. Our analysis of the ENCODE project data also found many transcriptionally active pseudogenes in the GM12878 and K562 cell lines; moreover, it showed that many human pseudogenes produced small RNAs (sRNAs) and some pseudogene-derived sRNAs, especially those from antisense strands, exhibited evidence of interfering with gene expression. Further integrated analysis of transcriptomics and epigenomics data, however, demonstrated that trimethylation of histone 3 at lysine 9 (H3K9me3), a posttranslational modification typically associated with gene repression and heterochromatin, was enriched at many transcribed pseudogenes in a transcription-level dependent manner in the two cell lines. The H3K9me3 enrichment was more prominent in pseudogenes that produced sRNAs at pseudogene loci and their adjacent regions, an observation further supported by the co-enrichment of SETDB1 (a H3K9 methyltransferase), suggesting that pseudogene sRNAs may have a role in regional chromatin repression. Taken together, our comprehensive and systematic characterization of pseudogene transcription uncovers a complex picture of how pseudogene ncRNAs could influence gene and pseudogene expression, at both epigenetic and post-transcriptional levels.