Sequence determinants of human gene regulatory elements.

DNA can determine where and when genes are expressed, but the full set of sequence determinants that control gene expression is unknown. Here, we measured the transcriptional activity of DNA sequences that represent an ~100 times larger sequence space than the human genome using massively parallel reporter assays (MPRAs). Machine learning models revealed that transcription factors (TFs) generally act in an additive manner with weak grammar and that most enhancers increase expression from a promoter by a mechanism that does not appear to involve specific TF-TF interactions. The enhancers themselves can be classified into three types: classical, closed chromatin and chromatin dependent. We also show that few TFs are strongly active in a cell, with most activities being similar between cell types. Individual TFs can have multiple gene regulatory activities, including chromatin opening and enhancing, promoting and determining transcription start site (TSS) activity, consistent with the view that the TF binding motif is the key atomic unit of gene expression.

Human White Adipose Tissue Displays Selective Insulin Resistance in the Obese State.

Selective hepatic insulin resistance is a feature of obesity and type 2 diabetes. Whether similar mechanisms operate in white adipose tissue (WAT) of those with obesity and to what extent these are normalized by weight loss are unknown. We determined insulin sensitivity by hyperinsulinemic euglycemic clamp and insulin response in subcutaneous WAT by RNA sequencing in 23 women with obesity before and 2 years after bariatric surgery. To control for effects of surgery, women postsurgery were matched to never-obese women. Multidimensional analyses of 138 samples allowed us to classify the effects of insulin into three distinct expression responses: a common set was present in all three groups and included genes encoding several lipid/cholesterol biosynthesis enzymes; a set of obesity-attenuated genes linked to tissue remodeling and protein translation was selectively regulated in the two nonobese states; and several postobesity-enriched genes encoding proteins involved in, for example, one-carbon metabolism were only responsive to insulin in the women who had lost weight. Altogether, human WAT displays a selective insulin response in the obese state, where most genes are normalized by weight loss. This comprehensive atlas provides insights into the transcriptional effects of insulin in WAT and may identify targets to improve insulin action.

Micro-biopsy for detection of gene expression changes in ischemic swine myocardium: A pilot study.

Micro-endomyocardial biopsy (micro-EMB) is a novel catheter-based biopsy technique, aiming to increase flexibility and safety compared to conventional EMB. The technique was developed and evaluated in healthy swine. Therefore, the ability to detect disease related tissue changes could not be evaluated. The aim of the present pilot study was to investigate the ability to detect disease related gene expression changes using micro-EMB. Myocardial infarction was induced in three swine by coronary artery balloon occlusion. Micro-EMB samples (n = 164) were collected before, during, and after occlusion. RNA-sequencing was performed on 85 samples, and 53 of these were selected for bioinformatic analysis. A large number of responding genes was detected from the infarcted area (n = 1911). The early responding genes (n = 1268) were mostly related to apoptosis and inflammation. There were fewer responding genes two days after infarction (n = 6), which were related to extra-cellular matrix changes, and none after 14 days. In contrast to the infarcted area, samples harvested from a non-infarcted myocardial region showed considerably fewer regulated genes (n = 33). Deconvolution analysis, to estimate the proportion of different cell types, revealed a higher proportion of fibroblasts and a reduced proportion of cardiomyocytes two days after occlusion compared to baseline (p < 0.02 and p < 0.01, respectively. S5 File). In conclusion, this pilot study demonstrates the capabilities of micro-EMB to detect local gene expression responses at an early stage after ischemia, but not at later timepoints.

Myocardial micro-biopsy procedure for molecular characterization with increased precision and reduced trauma.

Endomyocardial biopsy is a valuable tool in cardiac diagnostics but is limited by low diagnostic yield and significant complication risks. Meanwhile, recent developments in transcriptomic and proteomic technologies promise a wealth of biological data from minimal tissue samples. To take advantage of the minimal tissue amount needed for molecular analyses, we have developed a sub-millimeter endovascular biopsy device, considerably smaller than current clinical equipment, and devised a low-input RNA-sequencing protocol for analyzing small tissue samples. In in vivo evaluation in swine, 81% of biopsy attempts (n = 157) were successful. High quality RNA-sequencing data was generated from 91% of the sequenced cardiac micro-biopsy samples (n = 32). Gene expression signatures of samples taken with the novel device were comparable with a conventional device. No major complications were detected either during procedures or during 7 days' follow-up, despite acquiring a relatively large number of biopsies (median 30) in each animal. In conclusion, the novel device coupled with RNA-sequencing provides a feasible method to obtain molecular data from the myocardium. The method is less traumatic and has a higher flexibility compared to conventional methods, enabling safer and more targeted sampling from different parts of the myocardium.

A putative silencer variant in a spontaneous canine model of retinitis pigmentosa.

Retinitis pigmentosa (RP) is the leading cause of blindness with nearly two million people affected worldwide. Many genes have been implicated in RP, yet in 30-80% of the RP patients the genetic cause remains unknown. A similar phenotype, progressive retinal atrophy (PRA), affects many dog breeds including the Miniature Schnauzer. We performed clinical, genetic and functional experiments to identify the genetic cause of PRA in the breed. The age of onset and pattern of disease progression suggested that at least two forms of PRA, types 1 and 2 respectively, affect the breed, which was confirmed by genome-wide association study that implicated two distinct genomic loci in chromosomes 15 and X, respectively. Whole-genome sequencing revealed a fully segregating recessive regulatory variant in type 1 PRA. The associated variant has a very recent origin based on haplotype analysis and lies within a regulatory site with the predicted binding site of HAND1::TCF3 transcription factor complex. Luciferase assays suggested that mutated regulatory sequence increases expression. Case-control retinal expression comparison of six best HAND1::TCF3 target genes were analyzed with quantitative reverse-transcriptase PCR assay and indicated overexpression of EDN2 and COL9A2 in the affected retina. Defects in both EDN2 and COL9A2 have been previously associated with retinal degeneration. In summary, our study describes two genetically different forms of PRA and identifies a fully penetrant variant in type 1 form with a possible regulatory effect. This would be among the first reports of a regulatory variant in retinal degeneration in any species, and establishes a new spontaneous dog model to improve our understanding of retinal biology and gene regulation while the affected breed will benefit from a reliable genetic testing.

*-DCC: A platform to collect, annotate, and explore a large variety of sequencing experiments.

BACKGROUND: Over the past few years the variety of experimental designs and protocols for sequencing experiments increased greatly. To ensure the wide usability of the produced data beyond an individual project, rich and systematic annotation of the underlying experiments is crucial. FINDINGS: We first developed an annotation structure that captures the overall experimental design as well as the relevant details of the steps from the biological sample to the library preparation, the sequencing procedure, and the sequencing and processed files. Through various design features, such as controlled vocabularies and different field requirements, we ensured a high annotation quality, comparability, and ease of annotation. The structure can be easily adapted to a large variety of species. We then implemented the annotation strategy in a user-hosted web platform with data import, query, and export functionality. CONCLUSIONS: We present here an annotation structure and user-hosted platform for sequencing experiment data, suitable for lab-internal documentation, collaborations, and large-scale annotation efforts.

Identification and transfer of spatial transcriptomics signatures for cancer diagnosis.

BACKGROUND: Distinguishing ductal carcinoma in situ (DCIS) from invasive ductal carcinoma (IDC) regions in clinical biopsies constitutes a diagnostic challenge. Spatial transcriptomics (ST) is an in situ capturing method, which allows quantification and visualization of transcriptomes in individual tissue sections. In the past, studies have shown that breast cancer samples can be used to study their transcriptomes with spatial resolution in individual tissue sections. Previously, supervised machine learning methods were used in clinical studies to predict the clinical outcomes for cancer types. METHODS: We used four publicly available ST breast cancer datasets from breast tissue sections annotated by pathologists as non-malignant, DCIS, or IDC. We trained and tested a machine learning method (support vector machine) based on the expert annotation as well as based on automatic selection of cell types by their transcriptome profiles. RESULTS: We identified expression signatures for expert annotated regions (non-malignant, DCIS, and IDC) and build machine learning models. Classification results for 798 expression signature transcripts showed high coincidence with the expert pathologist annotation for DCIS (100%) and IDC (96%). Extending our analysis to include all 25,179 expressed transcripts resulted in an accuracy of 99% for DCIS and 98% for IDC. Further, classification based on an automatically identified expression signature covering all ST spots of tissue sections resulted in prediction accuracy of 95% for DCIS and 91% for IDC. CONCLUSIONS: This concept study suggest that the ST signatures learned from expert selected breast cancer tissue sections can be used to identify breast cancer regions in whole tissue sections including regions not trained on. Furthermore, the identified expression signatures can classify cancer regions in tissue sections not used for training with high accuracy. Expert-generated but even automatically generated cancer signatures from ST data might be able to classify breast cancer regions and provide clinical decision support for pathologists in the future.

Identification of novel cerebellar developmental transcriptional regulators with motif activity analysis.

BACKGROUND: The work of the FANTOM5 Consortium has brought forth a new level of understanding of the regulation of gene transcription and the cellular processes involved in creating diversity of cell types. In this study, we extended the analysis of the FANTOM5 Cap Analysis of Gene Expression (CAGE) transcriptome data to focus on understanding the genetic regulators involved in mouse cerebellar development. RESULTS: We used the HeliScopeCAGE library sequencing on cerebellar samples over 8 embryonic and 4 early postnatal times. This study showcases temporal expression pattern changes during cerebellar development. Through a bioinformatics analysis that focused on transcription factors, their promoters and binding sites, we identified genes that appear as strong candidates for involvement in cerebellar development. We selected several candidate transcriptional regulators for validation experiments including qRT-PCR and shRNA transcript knockdown. We observed marked and reproducible developmental defects in Atf4, Rfx3, and Scrt2 knockdown embryos, which support the role of these genes in cerebellar development. CONCLUSIONS: The successful identification of these novel gene regulators in cerebellar development demonstrates that the FANTOM5 cerebellum time series is a high-quality transcriptome database for functional investigation of gene regulatory networks in cerebellar development.

Promoter Usage and Dynamics in Vascular Smooth Muscle Cells Exposed to Fibroblast Growth Factor-2 or Interleukin-1ß.

Smooth muscle cells (SMC) in blood vessels are normally growth quiescent and transcriptionally inactive. Our objective was to understand promoter usage and dynamics in SMC acutely exposed to a prototypic growth factor or pro-inflammatory cytokine. Using cap analysis gene expression (FANTOM5 project) we report differences in promoter dynamics for immediate-early genes (IEG) and other genes when SMC are exposed to fibroblast growth factor-2 or interleukin-1ß. Of the 1871 promoters responding to FGF2 or IL-1ß considerably more responded to FGF2 (68.4%) than IL-1ß (18.5%) and 13.2% responded to both. Expression clustering reveals sets of genes induced, repressed or unchanged. Among IEG responding rapidly to FGF2 or IL-1ß were FOS, FOSB and EGR-1, which mediates human SMC migration. Motif activity response analysis (MARA) indicates most transcription factor binding motifs in response to FGF2 were associated with a sharp induction at 1 h, whereas in response to IL-1ß, most motifs were associated with a biphasic change peaking generally later. MARA revealed motifs for FOS_FOS{B,L1}_JUN{B,D} and EGR-1..3 in the cluster peaking 1 h after FGF2 exposure whereas these motifs were in clusters peaking 1 h or later in response to IL-1ß. Our findings interrogating CAGE data demonstrate important differences in promoter usage and dynamics in SMC exposed to FGF2 or IL-1ß.

Conserved temporal ordering of promoter activation implicates common mechanisms governing the immediate early response across cell types and stimuli.

The promoters of immediate early genes (IEGs) are rapidly activated in response to an external stimulus. These genes, also known as primary response genes, have been identified in a range of cell types, under diverse extracellular signals and using varying experimental protocols. Whereas genomic dissection on a case-by-case basis has not resulted in a comprehensive catalogue of IEGs, a rigorous meta-analysis of eight genome-wide FANTOM5 CAGE (cap analysis of gene expression) time course datasets reveals successive waves of promoter activation in IEGs, recapitulating known relationships between cell types and stimuli: we obtain a set of 57 (42 protein-coding) candidate IEGs possessing promoters that consistently drive a rapid but transient increase in expression over time. These genes show significant enrichment for known IEGs reported previously, pathways associated with the immediate early response, and include a number of non-coding RNAs with roles in proliferation and differentiation. Surprisingly, we also find strong conservation of the ordering of activation for these genes, such that 77 pairwise promoter activation orderings are conserved. Using the leverage of comprehensive CAGE time series data across cell types, we also document the extensive alternative promoter usage by such genes, which is likely to have been a barrier to their discovery until now. The common activation ordering of the core set of early-responding genes we identify may indicate conserved underlying regulatory mechanisms. By contrast, the considerably larger number of transiently activated genes that are specific to each cell type and stimulus illustrates the breadth of the primary response.

Characterization of the human RFX transcription factor family by regulatory and target gene analysis.

BACKGROUND: Evolutionarily conserved RFX transcription factors (TFs) regulate their target genes through a DNA sequence motif called the X-box. Thereby they regulate cellular specialization and terminal differentiation. Here, we provide a comprehensive analysis of all the eight human RFX genes (RFX1-8), their spatial and temporal expression profiles, potential upstream regulators and target genes. RESULTS: We extracted all known human RFX1-8 gene expression profiles from the FANTOM5 database derived from transcription start site (TSS) activity as captured by Cap Analysis of Gene Expression (CAGE) technology. RFX genes are broadly (RFX1-3, RFX5, RFX7) and specifically (RFX4, RFX6) expressed in different cell types, with high expression in four organ systems: immune system, gastrointestinal tract, reproductive system and nervous system. Tissue type specific expression profiles link defined RFX family members with the target gene batteries they regulate. We experimentally confirmed novel TSS locations and characterized the previously undescribed RFX8 to be lowly expressed. RFX tissue and cell type specificity arises mainly from differences in TSS architecture. RFX transcript isoforms lacking a DNA binding domain (DBD) open up new possibilities for combinatorial target gene regulation. Our results favor a new grouping of the RFX family based on protein domain composition. We uncovered and experimentally confirmed the TFs SP2 and ESR1 as upstream regulators of specific RFX genes. Using TF binding profiles from the JASPAR database, we determined relevant patterns of X-box motif positioning with respect to gene TSS locations of human RFX target genes. CONCLUSIONS: The wealth of data we provide will serve as the basis for precisely determining the roles RFX TFs play in human development and disease.

Correction to: Relatively frequent switching of transcription start sites during cerebellar development.

CORRECTION: The authors of the original article [1] would like to recognize the critical contribution of core members of the FANTOM5 Consortium, who played the critical role of HeliScopeCAGE sequencing experiments, quality control of tag reads and processing of the raw sequencing data.

Discovery of Transcription Factors Novel to Mouse Cerebellar Granule Cell Development Through Laser-Capture Microdissection.

Laser-capture microdissection was used to isolate external germinal layer tissue from three developmental periods of mouse cerebellar development: embryonic days 13, 15, and 18. The cerebellar granule cell-enriched mRNA library was generated with next-generation sequencing using the Helicos technology. Our objective was to discover transcriptional regulators that could be important for the development of cerebellar granule cells-the most numerous neuron in the central nervous system. Through differential expression analysis, we have identified 82 differentially expressed transcription factors (TFs) from a total of 1311 differentially expressed genes. In addition, with TF-binding sequence analysis, we have identified 46 TF candidates that could be key regulators responsible for the variation in the granule cell transcriptome between developmental stages. Altogether, we identified 125 potential TFs (82 from differential expression analysis, 46 from motif analysis with 3 overlaps in the two sets). From this gene set, 37 TFs are considered novel due to the lack of previous knowledge about their roles in cerebellar development. The results from transcriptome-wide analyses were validated with existing online databases, qRT-PCR, and in situ hybridization. This study provides an initial insight into the TFs of cerebellar granule cells that might be important for development and provide valuable information for further functional studies on these transcriptional regulators.

Long-Term Improvement in Aortic Pulse Wave Velocity After Weight Loss Can Be Predicted by White Adipose Tissue Factors.

BACKGROUND: Arterial stiffness, measured by pulse wave velocity (PWV), is linked to obesity, cardiovascular disease, and all-cause mortality. Short-term weight loss improves PWV, but the long-term effects are unknown. We investigated the effect of pronounced long-term weight loss on PWV and whether anthropometric/metabolic parameters and/or white adipose tissue (WAT) phenotype could predict this change in PWV. METHODS: Eighty-two obese subjects were examined before and 2 years after Roux-en-Y gastric bypass. Analyses included anthropometrics, routine clinical chemistry, and hyperinsulinemic-euglycemic clamp. Arterial stiffness was measured as aortic PWV (aPWV) using the Arteriograph device. WAT mass and distribution were assessed by dual-X-ray absorptiometry. Baseline visceral and subcutaneous WAT samples were obtained to measure adipocyte cell size. Transcriptomic profiling of subcutaneous WAT was performed in a subset of subjects (n = 30). RESULTS: At the 2-year follow-up, there were significant decreases in body mass index (39.4 ± 3.5 kg/m2 vs. 26.6 ± 3.4 kg/m2; P < 0.0001) and aPWV (7.8 ± 1.5 m/s vs. 7.2 ± 1.4 m/s; P = 0.006). Multiple regression analyses showed that baseline subcutaneous adipocyte volume was associated with a reduction in aPWV (P = 0.014), after adjusting for confounders. Expression analyses of 52 genes implicated in arterial stiffness showed that only one, COL4A1, independently predicted improvements in aPWV after adjusting for confounders (P = 0.006). CONCLUSIONS: Bariatric surgery leads to long-term reduction in aPWV. This improvement can be independently predicted by subcutaneous adipocyte volume and WAT COL4A1 expression, which suggests that subcutaneous WAT has a role in regulating aPWV. CLINICAL TRIALS REGISTRATION: Trial Number NCT01727245 (clinicaltrials.gov).

Acute doses of caffeine shift nervous system cell expression profiles toward promotion of neuronal projection growth.

Caffeine is a widely consumed psychoactive substance, but little is known about the effects of caffeine stimulation on global gene expression changes in neurons. Here, we conducted gene expression profiling of human neuroepithelial stem cell-derived neurons, stimulated with normal consumption levels of caffeine (3 µM and 10 µM), over a period of 9 h. We found dosage-dependent activation of immediate early genes after 1 h. Neuronal projection development processes were up-regulated and negative regulation of axon extension processes were down-regulated at 3 h. In addition, genes involved in extracellular matrix organization, response for wound healing, and regulation of immune system processes were down-regulated by caffeine at 3 h. This study identified novel genes within the neuronal projection guidance pathways that respond to acute caffeine stimulation and suggests potential mechanisms for the effects of caffeine on neuronal cells.

An integrated expression atlas of miRNAs and their promoters in human and mouse.

MicroRNAs (miRNAs) are short non-coding RNAs with key roles in cellular regulation. As part of the fifth edition of the Functional Annotation of Mammalian Genome (FANTOM5) project, we created an integrated expression atlas of miRNAs and their promoters by deep-sequencing 492 short RNA (sRNA) libraries, with matching Cap Analysis Gene Expression (CAGE) data, from 396 human and 47 mouse RNA samples. Promoters were identified for 1,357 human and 804 mouse miRNAs and showed strong sequence conservation between species. We also found that primary and mature miRNA expression levels were correlated, allowing us to use the primary miRNA measurements as a proxy for mature miRNA levels in a total of 1,829 human and 1,029 mouse CAGE libraries. We thus provide a broad atlas of miRNA expression and promoters in primary mammalian cells, establishing a foundation for detailed analysis of miRNA expression patterns and transcriptional control regions.

An atlas of human long non-coding RNAs with accurate 5' ends.

Long non-coding RNAs (lncRNAs) are largely heterogeneous and functionally uncharacterized. Here, using FANTOM5 cap analysis of gene expression (CAGE) data, we integrate multiple transcript collections to generate a comprehensive atlas of 27,919 human lncRNA genes with high-confidence 5' ends and expression profiles across 1,829 samples from the major human primary cell types and tissues. Genomic and epigenomic classification of these lncRNAs reveals that most intergenic lncRNAs originate from enhancers rather than from promoters. Incorporating genetic and expression data, we show that lncRNAs overlapping trait-associated single nucleotide polymorphisms are specifically expressed in cell types relevant to the traits, implicating these lncRNAs in multiple diseases. We further demonstrate that lncRNAs overlapping expression quantitative trait loci (eQTL)-associated single nucleotide polymorphisms of messenger RNAs are co-expressed with the corresponding messenger RNAs, suggesting their potential roles in transcriptional regulation. Combining these findings with conservation data, we identify 19,175 potentially functional lncRNAs in the human genome.

Transcriptional Dynamics During Human Adipogenesis and Its Link to Adipose Morphology and Distribution.

White adipose tissue (WAT) can develop into several phenotypes with different pathophysiological impact on type 2 diabetes. To better understand the adipogenic process, the transcriptional events that occur during in vitro differentiation of human adipocytes were investigated and the findings linked to WAT phenotypes. Single-molecule transcriptional profiling provided a detailed map of the expressional changes of genes, enhancers, and long noncoding RNAs, where different types of transcripts share common dynamics during differentiation. Common signatures include early downregulated, transient, and late induced transcripts, all of which are linked to distinct developmental processes during adipogenesis. Enhancers expressed during adipogenesis overlap significantly with genetic variants associated with WAT distribution. Transiently expressed and late induced genes are associated with hypertrophic WAT (few but large fat cells), a phenotype closely linked to insulin resistance and type 2 diabetes. Transcription factors that are expressed early or transiently affect differentiation and adipocyte function and are controlled by several well-known upstream regulators such as glucocorticosteroids, insulin, cAMP, and thyroid hormones. Taken together, our results suggest a complex but highly coordinated regulation of adipogenesis.

The Adipose Transcriptional Response to Insulin Is Determined by Obesity, Not Insulin Sensitivity.

Metabolically healthy obese subjects display preserved insulin sensitivity and a beneficial white adipose tissue gene expression pattern. However, this observation stems from fasting studies when insulin levels are low. We investigated adipose gene expression by 5'Cap-mRNA sequencing in 17 healthy non-obese (NO), 21 insulin-sensitive severely obese (ISO), and 30 insulin-resistant severely obese (IRO) subjects, before and 2 hr into a hyperinsulinemic euglycemic clamp. ISO and IRO subjects displayed a clear but globally similar transcriptional response to insulin, which differed from the small effects observed in NO subjects. In the obese, 231 genes were altered; 71 were enriched in ISO subjects (e.g., phosphorylation processes), and 52 were enriched in IRO subjects (e.g., cellular stimuli). Common cardio-metabolic risk factors and gender do not influence these findings. This study demonstrates that differences in the acute transcriptional response to insulin are primarily driven by obesity per se, challenging the notion of healthy obese adipose tissue, at least in severe obesity.