Detailed characterization of the mouse embryonic stem cell transcriptome reveals novel genes and intergenic splicing associated with pluripotency

<p>Abstract</p> <p>Background</p> <p>Transcriptional control of embryonic stem (ES) cell pluripotency has been a subject of intense study. Transcriptional regulators including Oct4 (Oct3/4 index), Sox2 and Nanog are fundamental for maintaining the undifferentiated state. However, the ES cell transcriptome is not limited to their targets, and exhibits considerable complexity when assayed with microarray, MPSS, cDNA/EST sequencing, and SAGE technologies. To identify novel genes associated with pluripotency, we globally searched for ES transcripts not corresponding to known genes, validated their sequences, determined their expression profiles, and employed RNAi to test their function.</p> <p>Results</p> <p>Gene Identification Signature (GIS) analysis, a SAGE derivative distinguished by paired 5' and 3' transcript end tags, identified 153 candidate novel transcriptional units (TUs) distinct from known genes in a mouse E14 ES mRNA library. We focused on 16 TUs free of artefacts and mapping discrepancies, five of which were validated by RTPCR product sequencing. Two of the TUs were revealed by annotation to represent novel protein-coding genes: a PRY-domain cluster member and a KRAB-domain zinc finger. The other three TUs represented intergenic splicing events involving adjacent, functionally unrelated protein-coding genes transcribed in the same orientation, with one event potentially encoding a fusion protein containing domains from both component genes (Clk2 and Scamp3). Expression profiling using embryonic samples and adult tissue panels confirmed that three of the TUs were unique to or most highly expressed in ES cells. Expression levels of all five TUs dropped dramatically during three distinct chemically induced differentiation treatments of ES cells in culture. However, siRNA knockdowns of the TUs did not alter mRNA levels of pluripotency or differentiation markers, and did not affect cell morphology.</p> <p>Conclusion</p> <p>Transcriptome libraries retain considerable potential for novel gene discovery despite massive recent cDNA and EST sequencing efforts; cDNA and EST evidence for these ES cell TUs had been limited or absent. RTPCR and full-length sequencing remain essential in resolving the bottleneck between numerous candidate novel transcripts inferred from high-throughput sequencing and the small fraction that can be validated. RNAi results indicate that, despite their strong association with pluripotency, these five transcriptomic novelties may not be required for maintaining it.</p

Non-coding regions of nuclear-DNA-encoded mitochondrial genes and intergenic sequences are targeted by autoantibodies in breast cancer

Autoantibodies against mitochondrial-derived antigens play a key role in chronic tissue inflammation in autoimmune disorders and cancers. Here, we identify autoreactive nuclear genomic DNA (nDNA)-encoded mitochondrial gene products (GAPDH, PKM2, GSTP1, SPATA5, MFF, TSPOAP1, PHB2, COA4, and HAGH) recognized by breast cancer (BC) patients\u27 sera as nonself, supporting a direct relationship of mitochondrial autoimmunity to breast carcinogenesis. Autoreactivity of multiple nDNA-encoded mitochondrial gene products was mapped to protein-coding regions, 3\u27 untranslated regions (UTRs), as well as introns. In addition, autoantibodies in BC sera targeted intergenic sequences that may be parts of long non-coding RNA (lncRNA) genes, including LINC02381 and other putative lncRNA neighbors of the protein-coding genes ERCC4, CXCL13, SOX3, PCDH1, EDDM3B, and GRB2. Increasing evidence indicates that lncRNAs play a key role in carcinogenesis. Consistent with this, our findings suggest that lncRNAs, as well as mRNAs of nDNA-encoded mitochondrial genes, mechanistically contribute to BC progression. This work supports a new paradigm of breast carcinogenesis based on a globally dysfunctional genome with altered function of multiple mitochondrial and non-mitochondrial oncogenic pathways caused by the effects of autoreactivity-induced dysregulation of multiple genes and their products. This autoimmunity-based model of carcinogenesis will open novel avenues for BC treatment

Mice and Men: Their Promoter Properties

Using the two largest collections of Mus musculus and Homo sapiens transcription start sites (TSSs) determined based on CAGE tags, ditags, full-length cDNAs, and other transcript data, we describe the compositional landscape surrounding TSSs with the aim of gaining better insight into the properties of mammalian promoters. We classified TSSs into four types based on compositional properties of regions immediately surrounding them. These properties highlighted distinctive features in the extended core promoters that helped us delineate boundaries of the transcription initiation domain space for both species. The TSS types were analyzed for associations with initiating dinucleotides, CpG islands, TATA boxes, and an extensive collection of statistically significant cis-elements in mouse and human. We found that different TSS types show preferences for different sets of initiating dinucleotides and cis-elements. Through Gene Ontology and eVOC categories and tissue expression libraries we linked TSS characteristics to expression. Moreover, we show a link of TSS characteristics to very specific genomic organization in an example of immune-response-related genes (GO:0006955). Our results shed light on the global properties of the two transcriptomes not revealed before and therefore provide the framework for better understanding of the transcriptional mechanisms in the two species, as well as a framework for development of new and more efficient promoter- and gene-finding tools

YY1 directly interacts with myocardin to repress the triad myocardin/SRF/CArG box-mediated smooth muscle gene transcription during smooth muscle phenotypic modulation

Yin Yang 1 (YY1) regulates gene transcription in a variety of biological processes. In this study, we aim to determine the role of YY1 in vascular smooth muscle cell (VSMC) phenotypic modulation both in vivo and in vitro. Here we show that vascular injury in rodent carotid arteries induces YY1 expression along with reduced expression of smooth muscle differentiation markers in the carotids. Consistent with this finding, YY1 expression is induced in differentiated VSMCs in response to serum stimulation. To determine the underlying molecular mechanisms, we found that YY1 suppresses the transcription of CArG box-dependent SMC-specific genes including SM22α, SMα-actin and SMMHC. Interestingly, YY1 suppresses the transcriptional activity of the SM22α promoter by hindering the binding of serum response factor (SRF) to the proximal CArG box. YY1 also suppresses the transcription and the transactivation of myocardin (MYOCD), a master regulator for SMC-specific gene transcription by binding to SRF to form the MYOCD/SRF/CArG box triad (known as the ternary complex). Mechanistically, YY1 directly interacts with MYOCD to competitively displace MYOCD from SRF. This is the first evidence showing that YY1 inhibits SMC differentiation by directly targeting MYOCD. These findings provide new mechanistic insights into the regulatory mechanisms that govern SMC phenotypic modulation in the pathogenesis of vascular diseases

Complex Loci in Human and Mouse Genomes

Mammalian genomes harbor a larger than expected number of complex loci, in which multiple genes are coupled by shared transcribed regions in antisense orientation and/or by bidirectional core promoters. To determine the incidence, functional significance, and evolutionary context of mammalian complex loci, we identified and characterized 5,248 cis–antisense pairs, 1,638 bidirectional promoters, and 1,153 chains of multiple cis–antisense and/or bidirectionally promoted pairs from 36,606 mouse transcriptional units (TUs), along with 6,141 cis–antisense pairs, 2,113 bidirectional promoters, and 1,480 chains from 42,887 human TUs. In both human and mouse, 25% of TUs resided in cis–antisense pairs, only 17% of which were conserved between the two organisms, indicating frequent species specificity of antisense gene arrangements. A sampling approach indicated that over 40% of all TUs might actually be in cis–antisense pairs, and that only a minority of these arrangements are likely to be conserved between human and mouse. Bidirectional promoters were characterized by variable transcriptional start sites and an identifiable midpoint at which overall sequence composition changed strand and the direction of transcriptional initiation switched. In microarray data covering a wide range of mouse tissues, genes in cis–antisense and bidirectionally promoted arrangement showed a higher probability of being coordinately expressed than random pairs of genes. In a case study on homeotic loci, we observed extensive transcription of nonconserved sequences on the noncoding strand, implying that the presence rather than the sequence of these transcripts is of functional importance. Complex loci are ubiquitous, host numerous nonconserved gene structures and lineage-specific exonification events, and may have a cis-regulatory impact on the member genes

Whole-Genome Cartography of Estrogen Receptor α Binding Sites

Using a chromatin immunoprecipitation-paired end diTag cloning and sequencing strategy, we mapped estrogen receptor α (ERα) binding sites in MCF-7 breast cancer cells. We identified 1,234 high confidence binding clusters of which 94% are projected to be bona fide ERα binding regions. Only 5% of the mapped estrogen receptor binding sites are located within 5 kb upstream of the transcriptional start sites of adjacent genes, regions containing the proximal promoters, whereas vast majority of the sites are mapped to intronic or distal locations (>5 kb from 5′ and 3′ ends of adjacent transcript), suggesting transcriptional regulatory mechanisms over significant physical distances. Of all the identified sites, 71% harbored putative full estrogen response elements (EREs), 25% bore ERE half sites, and only 4% had no recognizable ERE sequences. Genes in the vicinity of ERα binding sites were enriched for regulation by estradiol in MCF-7 cells, and their expression profiles in patient samples segregate ERα-positive from ERα-negative breast tumors. The expression dynamics of the genes adjacent to ERα binding sites suggest a direct induction of gene expression through binding to ERE-like sequences, whereas transcriptional repression by ERα appears to be through indirect mechanisms. Our analysis also indicates a number of candidate transcription factor binding sites adjacent to occupied EREs at frequencies much greater than by chance, including the previously reported FOXA1 sites, and demonstrate the potential involvement of one such putative adjacent factor, Sp1, in the global regulation of ERα target genes. Unexpectedly, we found that only 22%–24% of the bona fide human ERα binding sites were overlapping conserved regions in whole genome vertebrate alignments, which suggest limited conservation of functional binding sites. Taken together, this genome-scale analysis suggests complex but definable rules governing ERα binding and gene regulation

Non-coding regions of nuclear-DNA-encoded mitochondrial genes and intergenic sequences are targeted by autoantibodies in breast cancer

Autoantibodies against mitochondrial-derived antigens play a key role in chronic tissue inflammation in autoimmune disorders and cancers. Here, we identify autoreactive nuclear genomic DNA (nDNA)-encoded mitochondrial gene products (GAPDH, PKM2, GSTP1, SPATA5, MFF, TSPOAP1, PHB2, COA4, and HAGH) recognized by breast cancer (BC) patients’ sera as nonself, supporting a direct relationship of mitochondrial autoimmunity to breast carcinogenesis. Autoreactivity of multiple nDNA-encoded mitochondrial gene products was mapped to protein-coding regions, 3’ untranslated regions (UTRs), as well as introns. In addition, autoantibodies in BC sera targeted intergenic sequences that may be parts of long non-coding RNA (lncRNA) genes, including LINC02381 and other putative lncRNA neighbors of the protein-coding genes ERCC4, CXCL13, SOX3, PCDH1, EDDM3B, and GRB2. Increasing evidence indicates that lncRNAs play a key role in carcinogenesis. Consistent with this, our findings suggest that lncRNAs, as well as mRNAs of nDNA-encoded mitochondrial genes, mechanistically contribute to BC progression. This work supports a new paradigm of breast carcinogenesis based on a globally dysfunctional genome with altered function of multiple mitochondrial and non-mitochondrial oncogenic pathways caused by the effects of autoreactivity-induced dysregulation of multiple genes and their products. This autoimmunity-based model of carcinogenesis will open novel avenues for BC treatment

De novo sequencing, assembly and analysis of eight different transcriptomes from the Malayan pangolin

Pangolins are scale-covered mammals, containing eight endangered species. Maintaining pangolins in captivity is a significant challenge, in part because little is known about their genetics. Here we provide the first large-scale sequencing of the critically endangered Manis javanica transcriptomes from eight different organs using Illumina HiSeq technology, yielding ~75 Giga bases and 89,754 unigenes. We found some unigenes involved in the insect hormone biosynthesis pathway and also 747 lipids metabolism-related unigenes that may be insightful to understand the lipid metabolism system in pangolins. Comparative analysis between M. javanica and other mammals revealed many pangolin-specific genes significantly over-represented in stress-related processes, cell proliferation and external stimulus, probably reflecting the traits and adaptations of the analyzed pregnant female M. javanica. Our study provides an invaluable resource for future functional works that may be highly relevant for the conservation of pangolins