Genome-wide strategies identify downstream target genes of chick connective tissue-associated transcription factors.

Connective tissues support organs and play crucial roles in development, homeostasis and fibrosis, yet our understanding of their formation is still limited. To gain insight into the molecular mechanisms of connective tissue specification, we selected five zinc-finger transcription factors - OSR1, OSR2, EGR1, KLF2 and KLF4 - based on their expression patterns and/or known involvement in connective tissue subtype differentiation. RNA-seq and ChIP-seq profiling of chick limb micromass cultures revealed a set of common genes regulated by all five transcription factors, which we describe as a connective tissue core expression set. This common core was enriched with genes associated with axon guidance and myofibroblast signature, including fibrosis-related genes. In addition, each transcription factor regulated a specific set of signalling molecules and extracellular matrix components. This suggests a concept whereby local molecular niches can be created by the expression of specific transcription factors impinging on the specification of local microenvironments. The regulatory network established here identifies common and distinct molecular signatures of limb connective tissue subtypes, provides novel insight into the signalling pathways governing connective tissue specification, and serves as a resource for connective tissue development.

A dual transcript-discovery approach to improve the delimitation of gene features from RNA-seq data in the chicken model.

The sequence of the chicken genome, like several other draft genome sequences, is presently not fully covered. Gaps, contigs assigned with low confidence and uncharacterized chromosomes result in gene fragmentation and imprecise gene annotation. Transcript abundance estimation from RNA sequencing (RNA-seq) data relies on read quality, library complexity and expression normalization. In addition, the quality of the genome sequence used to map sequencing reads, and the gene annotation that defines gene features, must also be taken into account. A partially covered genome sequence causes the loss of sequencing reads from the mapping step, while an inaccurate definition of gene features induces imprecise read counts from the assignment step. Both steps can significantly bias interpretation of RNA-seq data. Here, we describe a dual transcript-discovery approach combining a genome-guided gene prediction and a de novo transcriptome assembly. This dual approach enabled us to increase the assignment rate of RNA-seq data by nearly 20% as compared to when using only the chicken reference annotation, contributing therefore to a more accurate estimation of transcript abundance. More generally, this strategy could be applied to any organism with partial genome sequence and/or lacking a manually-curated reference annotation in order to improve the accuracy of gene expression studies.

Odd skipped-related 1 identifies a population of embryonic fibro-adipogenic progenitors regulating myogenesis during limb development.

Fibro-adipogenic progenitors (FAPs) are an interstitial cell population in adult skeletal muscle that support muscle regeneration. During development, interstitial muscle connective tissue (MCT) cells support proper muscle patterning, however the underlying molecular mechanisms are not well understood and it remains unclear whether adult FAPs and embryonic MCT cells share a common lineage. We show here that mouse embryonic limb MCT cells expressing the transcription factor Osr1, differentiate into fibrogenic and adipogenic cells in vivo and in vitro defining an embryonic FAP-like population. Genetic lineage tracing shows that developmental Osr1+ cells give rise to a subset of adult FAPs. Loss of Osr1 function leads to a reduction of myogenic progenitor proliferation and survival resulting in limb muscle patterning defects. Transcriptome and functional analyses reveal that Osr1+ cells provide a critical pro-myogenic niche via the production of MCT specific extracellular matrix components and secreted signaling factors.

miR-22 and miR-29a Are Members of the Androgen Receptor Cistrome Modulating LAMC1 and Mcl-1 in Prostate Cancer.

The normal prostate as well as early stages and advanced prostate cancer (PCa) require a functional androgen receptor (AR) for growth and survival. The recent discovery of microRNAs (miRNAs) as novel effector molecules of AR disclosed the existence of an intricate network between AR, miRNAs and downstream target genes. In this study DUCaP cells, characterized by high content of wild-type AR and robust AR transcriptional activity, were chosen as the main experimental model. By integrative analysis of chromatin immunoprecipitation-sequencing (ChIP-seq) and microarray expression profiling data, miRNAs putatively bound and significantly regulated by AR were identified. A direct AR regulation of miR-22, miR-29a, and miR-17-92 cluster along with their host genes was confirmed. Interestingly, endogenous levels of miR-22 and miR-29a were found to be reduced in PCa cells expressing AR. In primary tumor samples, miR-22 and miR-29a were less abundant in the cancerous tissue compared with the benign counterpart. This specific expression pattern was associated with a differential DNA methylation of the genomic AR binding sites. The identification of laminin gamma 1 (LAMC1) and myeloid cell leukemia 1 (MCL1) as direct targets of miR-22 and miR-29a, respectively, suggested a tumor-suppressive role of these miRNAs. Indeed, transfection of miRNA mimics in PCa cells induced apoptosis and diminished cell migration and viability. Collectively, these data provide additional information regarding the complex regulatory machinery that guides miRNAs activity in PCa, highlighting an important contribution of miRNAs in the AR signaling.

High-throughput miRNA and mRNA sequencing of paired colorectal normal, tumor and metastasis tissues and bioinformatic modeling of miRNA-1 therapeutic applications.

MiRNAs are discussed as diagnostic and therapeutic molecules. However, effective miRNA drug treatments with miRNAs are, so far, hampered by the complexity of the miRNA networks. To identify potential miRNA drugs in colorectal cancer, we profiled miRNA and mRNA expression in matching normal, tumor and metastasis tissues of eight patients by Illumina sequencing. We validated six miRNAs in a large tissue screen containing 16 additional tumor entities and identified miRNA-1, miRNA-129, miRNA-497 and miRNA-215 as constantly de-regulated within the majority of cancers. Of these, we investigated miRNA-1 as representative in a systems-biology simulation of cellular cancer models implemented in PyBioS and assessed the effects of depletion as well as overexpression in terms of miRNA-1 as a potential treatment option. In this system, miRNA-1 treatment reverted the disease phenotype with different effectiveness among the patients. Scoring the gene expression changes obtained through mRNA-Seq from the same patients we show that the combination of deep sequencing and systems biological modeling can help to identify patient-specific responses to miRNA treatments. We present this data as guideline for future pre-clinical assessments of new and personalized therapeutic options.

ERG induces epigenetic activation of Tudor domain-containing protein 1 (TDRD1) in ERG rearrangement-positive prostate cancer.

BACKGROUND: Overexpression of ERG transcription factor due to genomic ERG-rearrangements defines a separate molecular subtype of prostate tumors. One of the consequences of ERG accumulation is modulation of the cell's gene expression profile. Tudor domain-containing protein 1 gene (TDRD1) was reported to be differentially expressed between TMPRSS2:ERG-negative and TMPRSS2:ERG-positive prostate cancer. The aim of our study was to provide a mechanistic explanation for the transcriptional activation of TDRD1 in ERG rearrangement-positive prostate tumors. METHODOLOGY/PRINCIPAL FINDINGS: Gene expression measurements by real-time quantitative PCR revealed a remarkable co-expression of TDRD1 and ERG (r(2) = 0.77) but not ETV1 (r(2)<0.01) in human prostate cancer in vivo. DNA methylation analysis by MeDIP-Seq and bisulfite sequencing showed that TDRD1 expression is inversely correlated with DNA methylation at the TDRD1 promoter in vitro and in vivo (ρ = -0.57). Accordingly, demethylation of the TDRD1 promoter in TMPRSS2:ERG-negative prostate cancer cells by DNA methyltransferase inhibitors resulted in TDRD1 induction. By manipulation of ERG dosage through gene silencing and forced expression we show that ERG governs loss of DNA methylation at the TDRD1 promoter-associated CpG island, leading to TDRD1 overexpression. CONCLUSIONS/SIGNIFICANCE: We demonstrate that ERG is capable of disrupting a tissue-specific DNA methylation pattern at the TDRD1 promoter. As a result, TDRD1 becomes transcriptionally activated in TMPRSS2:ERG-positive prostate cancer. Given the prevalence of ERG fusions, TDRD1 overexpression is a common alteration in human prostate cancer which may be exploited for diagnostic or therapeutic procedures.

Genome-wide DNA methylation events in TMPRSS2-ERG fusion-negative prostate cancers implicate an EZH2-dependent mechanism with miR-26a hypermethylation.

UNLABELLED: Prostate cancer is the second most common cancer among men worldwide. Alterations in the DNA methylation pattern can be one of the leading causes for prostate cancer formation. This study is the first high-throughput sequencing study investigating genome-wide DNA methylation patterns in a large cohort of 51 tumor and 53 benign prostate samples using methylated DNA immunoprecipitation sequencing. Comparative analyses identified more than 147,000 cancer-associated epigenetic alterations. In addition, global methylation patterns show significant differences based on the TMPRSS2-ERG rearrangement status. We propose the hypermethylation of miR-26a as an alternative pathway of ERG rearrangement-independent EZH2 activation. The observed increase in differential methylation events in fusion-negative tumors can explain the tumorigenic process in the absence of genomic rearrangements. SIGNIFICANCE: In contrast to TMPRSS2-ERG -rearranged tumors, the pathomechanism for gene fusion-negative tumors is completely unclear. Using a sequencing-based approach, our work uncovers significant global epigenetic alterations in TMPRSS2-ERG gene fusion-negative tumors and provides a mechanistic explanation for the tumor formation process.