Functional analysis and identification of cis-regulatory elements of human chromosome 21 gene promoters

Differentially Expressed Genes Related to Flowering Transition between Once- and Continuous-Flowering Roses

Roses are the most important cut flower crops and widely used woody ornamental plants in gardens throughout the world, and they are model plants for studying the continuous-flowering trait of woody plants. To analyze the molecular regulation mechanism of continuous flowering, comparative transcriptome data of once- and continuous-flowering roses in our previous study were used to conduct weighted gene co-expression network analysis (WGCNA) to obtain the candidate genes related to flowering transitions. The expression patterns of candidate genes at different developmental stages between Rosa chinensis “Old Blush” (continuous-flowering cultivar) and R. “Huan Die” (once-flowering cultivar) were investigated, and the relationship of the key gene with the endogenous hormone was analyzed. The results showed that the expression trends of VIN3-LIKE 1 (VIL1), FRIGIDA- LIKE 3 (FRI3), APETALA 2- LIKE (AP2-like) and CONSTANS-LIKE 2 (CO-like 2) genes were significantly different between “Old Blush” and “Huan Die”, and the expression trends of SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) and CO-like 2 were consistent in the flowering transition of “Old Blush” under different environments. The changes in cytokinin and gibberellic acid (GA₃) content were different in the two rose cultivars. The overall change trend of the abscisic acid and GA₃ in the flowering transition of “Old Blush” under different environments was consistent. The promoter sequence of CO-like 2 contained a P-box element associated with gibberellin response, as well as binding sites for transcription factors. In a word, we found CO-like 2 associated with continuous flowering and some factors that may synergistically regulate continuous flowering. The results provided a reference for elucidating the molecular regulatory mechanisms of continuous-flowering traits in roses.

Computational annotation of miRNA transcription start sites

Motivation

MicroRNAs (miRNAs) are small noncoding RNAs that play important roles in gene regulation and phenotype development. The identification of miRNA transcription start sites (TSSs) is critical to understand the functional roles of miRNA genes and their transcriptional regulation. Unlike protein-coding genes, miRNA TSSs are not directly detectable from conventional RNA-Seq experiments due to miRNA-specific process of biogenesis. In the past decade, large-scale genome-wide TSS-Seq and transcription activation marker profiling data have become available, based on which, many computational methods have been developed. These methods have greatly advanced genome-wide miRNA TSS annotation.

Results

In this study, we summarized recent computational methods and their results on miRNA TSS annotation. We collected and performed a comparative analysis of miRNA TSS annotations from 14 representative studies. We further compiled a robust set of miRNA TSSs (RSmirT) that are supported by multiple studies. Integrative genomic and epigenomic data analysis on RSmirT revealed the genomic and epigenomic features of miRNA TSSs as well as their relations to protein-coding and long non-coding genes.

Contact

xiaoman@mail.ucf.edu, haihu@cs.ucf.edu

Identification of novel long non-coding RNA biomarkers for prognosis prediction of papillary thyroid cancer

Papillary thyroid carcinoma (PTC) is the most frequent type of malignant thyroid tumor. Several lncRNA signatures have been established for prognosis prediction in some cancers. However, the prognostic value of lncRNAs has not been investigated in PTC yet. In this study, we performed genome-wide analysis of lncRNA expression profiles in a large cohort of PTC patients from The Cancer Genome Atlas and identified 111 differentially expressed lncRNAs between tumor and non-tumor samples and between recurrent and recurrence-free samples. From the 111 differentially expressed lncRNAs, four independent lncRNA biomarkers associated with prognosis were identified and were integrated into a four-lncRNA signature which classified the patients of training dataset into the high-risk group and low-risk group with significantly different overall survival (p=0.016, log-rank test). The prognostic value of the four-lncRNA signature was validated in the independent testing dataset. Multivariate analysis indicated that the four-lncRNA signature was an independent prognostic predictor. Moreover, identified four lncRNA biomarkers demonstrates good performance in predicting disease recurrence with AUC of 0.833 using leave one out cross-validation. Our study not only highlighted the potential role for lncRNAs to improve clinical prognosis prediction in patients with PTC and but also provided alternative biomarkers and therapeutic targets for PTC patients.

Transcriptional interference by small transcripts in proximal promoter regions

Proximal promoter regions (PPR) are heavily transcribed yielding different types of small RNAs. The act of transcription within PPRs might regulate downstream gene expression via transcriptional interference (TI). For analysis, we investigated capped and polyadenylated small RNA transcripts within PPRs of human RefSeq genes in eight different cell lines. Transcripts of our datasets overlapped with experimentally determined transcription factor binding sites (TFBS). For TFBSs intersected by these small RNA transcripts, we established negative correlation of sRNA expression levels and transcription factor (TF) DNA binding affinities; suggesting that the transcripts acted via TI. Accordingly, datasets were designated as TFbiTrs (TF-binding interfering transcripts). Expression of most TFbiTrs was restricted to certain cell lines. This facilitated the analysis of effects related to TFbiTr expression for the same RefSeq genes across cell lines. We consistently uncovered higher relative TF/DNA binding affinities and concomitantly higher expression levels for RefSeq genes in the absence of TFbiTrs. Analysis of corresponding chromatin landscapes supported these results. ChIA-PET revealed the participation of distal enhancers in TFbiTr transcription. Enhancers regulating TFbiTrs, in effect, act as repressors for corresponding downstream RefSeq genes. We demonstrate the significant impact of TI on gene expression using selected small RNA datasets.

Human endogenous retrovirus HERV-K(HML-2) activity in prostate cancer is dominated by a few loci

BACKGROUND

Increased expression of human endogenous retroviruses, especially HERV-K(HML-2) proviruses, has recently been associated with prostate carcinoma progression. In particular, a HML-2 locus in chromosome 22q11.23 (H22q) is upregulated in many cases. We therefore aimed at delineating the extent and repertoire of HML-2 transcription in prostate cancer tissues and cell lines and to define the transcription pattern and biological effects of H22q.

METHODS

Sanger and high throughput amplicon sequencing was used to define the repertoire of expressed HML-2 in a selected set of samples. qRT-PCR was used to quantify expression of selected proviruses in an extended set of prostate cancer tissues. Transcription factor binding sites (TFBS) were compared bioinformatically using the Transfac database. Expression of H22q was further characterized by siRNA-mediated knockdown, 5' RACE mapping of transcriptional start sites (TSS) and identification of splice sites. Functional effects of H22q knockdown were investigated by viability and apoptosis assays.

RESULTS

In addition to H22q, a limited number of other proviruses were found expressed by sequencing. Of these, provirus ERVK-5 and to a lesser degree ERVK-15 were frequently upregulated in prostate cancer. In contrast, expression of ERVK-24, predominant in germ cell tumors, was not detectable in prostatic tissues. While HML-2 LTRs contain binding sites for the androgen receptor and cofactors, no consistent differences in transcription factor binding sites were found between expressed and non-expressed proviruses. The H22q locus contains two 5'-LTRs of which the upstream LTR is predominantly used in prostatic cells, with an imprecise TSS. Splicing of H22q transcripts is complex, generating, among others, a transcript with an Np9-like ORF. Knockdown of H22q did not significantly affect proliferation or apoptosis of prostate cancer cells.

CONCLUSIONS

Our findings further underline that HML-2 expression is commonly highly tissue-specific. In prostate cancer, a limited number of loci become activated, especially H22q and ERVK-5. As expressed and non-expressed proviruses do not differ significantly in TFBS, tissue- and tumor-specific expression may be governed primarily by chromatin context. Overexpression of HML-2 H22q is more likely consequence than cause of prostate cancer progression.

High-fat diet caused widespread epigenomic differences on hepatic methylome in rat

A high-fat (HF) diet is associated with progression of liver diseases. To illustrate genome-wide landscape of DNA methylation in liver of rats fed either a control or HF diet, two enrichment-based methods, namely methyl-DNA immunoprecipitation assay with high-throughput sequencing (MeDIP-seq) and methylation-sensitive restriction enzyme sequencing (MRE-seq), were performed in our study. Rats fed with the HF diet exhibited an increased body weight and liver fat accumulation compared with that of the control group when they were 12 wk of age. Genome-wide analysis of differentially methylated regions (DMRs) showed that 12,494 DMRs induced by HF diet were hypomethylated and 6,404 were hypermethylated. DMRs with gene annotations [differentially methylated genes (DMGs)] were further analyzed to show gene-specific methylation profile. There were 88, 2,680, and 95 hypomethylated DMGs identified with changes in DNA methylation in the promoter, intragenic and downstream regions, respectively, compared with fewer hypermethylated DMGs (45, 1,623, and 50 in the respective regions). Some of these genes also contained an ACGT cis-acting motif whose DNA methylation status may affect gene expression. Pathway analysis showed that these DMGs were involved in critical hepatic signaling networks related to hepatic development. Therefore, HF diet had global impacts on DNA methylation profile in the liver of rats, leading to differential expression of genes in hepatic pathways that may involve in functional changes in liver development.

Resveratrol compounds inhibit human holocarboxylase synthetase and cause a lean phenotype in Drosophila melanogaster

Holocarboxylase synthetase (HLCS) is the sole protein-biotin ligase in the human proteome. HLCS has key regulatory functions in intermediary metabolism, including fatty acid metabolism, and in gene repression through epigenetic mechanisms. The objective of this study was to identify food-borne inhibitors of HLCS that alter HLCS-dependent pathways in metabolism and gene regulation. When libraries of extracts from natural products and chemically pure compounds were screened for HLCS inhibitor activity, resveratrol compounds in grape materials caused an HLCS inhibition of >98% in vitro. The potency of these compounds was piceatannol>resveratrol>piceid. Grape-borne compounds other than resveratrol metabolites also contributed toward HLCS inhibition, e.g., p-coumaric acid and cyanidin chloride. HLCS inhibitors had meaningful effects on body fat mass. When Drosophila melanogaster brummer mutants, which are genetically predisposed to storing excess amounts of lipids, were fed diets enriched with grape leaf extracts and piceid, body fat mass decreased by more than 30% in males and females. However, Drosophila responded to inhibitor treatment with an increase in the expression of HLCS, which elicited an increase in the abundance of biotinylated carboxylases in vivo. We conclude that mechanisms other than inhibition of HLCS cause body fat loss in flies. We propose that the primary candidate is the inhibition of the insulin receptor/Akt signaling pathway.

Novel roles of holocarboxylase synthetase in gene regulation and intermediary metabolism

The role of holocarboxylase synthetase (HLCS) in catalyzing the covalent binding of biotin to the five biotin-dependent carboxylases in humans is well established, as are the essential roles of these carboxylases in the metabolism of fatty acids, the catabolism of leucine, and gluconeogenesis. This review examines recent discoveries regarding the roles of HLCS in assembling a multiprotein gene repression complex in chromatin. In addition, emerging evidence suggests that the number of biotinylated proteins is far larger than previously assumed and includes members of the heat-shock superfamily of proteins and proteins coded by the ENO1 gene. Evidence is presented linking biotinylation of heat-shock proteins HSP60 and HSP72 with redox biology and immune function, respectively, and biotinylation of the two ENO1 gene products MBP-1 and ENO1 with tumor suppression and glycolysis, respectively.

Three promoters regulate the transcriptional activity of the human holocarboxylase synthetase gene

Holocarboxylase synthetase (HLCS) is the only protein biotin ligase in the human proteome. HLCS-dependent biotinylation of carboxylases plays crucial roles in macronutrient metabolism. HLCS appears to be an essential part of multiprotein complexes in the chromatin that cause gene repression and contribute toward genome stability. Consistent with these essential functions, HLCS knockdown causes strong phenotypes including shortened life span and low stress resistance in Drosophila melanogaster, and de-repression of long-terminal repeats in humans, other mammalian cell lines and Drosophila. Despite previous observations that the expression of HLCS depends on biotin status in rats and in human cell lines, little is known about the regulation of HLCS expression. The goal of this study was to identify promoters that regulate the expression of the human HLCS gene. Initially, the human HLCS locus was interrogated in silico using predictors of promoters including sequences of HLCS mRNA and expressed sequence tags, CpG islands, histone marks denoting transcriptionally poised chromatin, transcription factor binding sites and DNaseI hypersensitive regions. Our predictions revealed three putative HLCS promoters, denoted P1, P2 and P3. Promoters lacked a TATA box, which is typical for housekeeping genes. When the three promoters were cloned into a luciferase reporter plasmid, reporter gene activity was at least three times background noise in human breast, colon and kidney cell lines; activities consistently followed the pattern P1>>P3>P2. Promoter activity depended on the concentration of biotin in culture media, but the effect was moderate. We conclude that we have identified promoters in the human HLCS gene.

MicroSCALE screening reveals genetic modifiers of therapeutic response in melanoma

Cell microarrays are a promising tool for performing large-scale functional genomic screening in mammalian cells at reasonable cost, but owing to technical limitations they have been restricted for use with a narrow range of cell lines and short-term assays. Here, we describe MicroSCALE (Microarrays of Spatially Confined Adhesive Lentiviral Features), a cell microarray-based platform that enables application of this technology to a wide range of cell types and longer-term assays. We used MicroSCALE to uncover kinases that when overexpressed partially desensitized B-RAFV600E-mutant melanoma cells to inhibitors of the mitogen-activated protein kinase kinase kinase (MAPKKK) RAF, the MAPKKs MEK1 and 2 (MEK1/2, mitogen-activated protein kinase kinase 1 and 2), mTOR (mammalian target of rapamycin), or PI3K (phosphatidylinositol 3-kinase). These screens indicated that cells treated with inhibitors acting through common mechanisms were affected by a similar profile of overexpressed proteins. In contrast, screens involving inhibitors acting through distinct mechanisms yielded unique profiles, a finding that has potential relevance for small-molecule target identification and combination drugging studies. Further, by integrating large-scale functional screening results with cancer cell line gene expression and pharmacological sensitivity data, we validated the nuclear factor κB pathway as a potential mediator of resistance to MAPK pathway inhibitors. The MicroSCALE platform described here may enable new classes of large-scale, resource-efficient screens that were not previously feasible, including those involving combinations of cell lines, perturbations, and assay outputs or those involving limited numbers of cells and limited or expensive reagents.

Transcription factor binding predictions using TRAP for the analysis of ChIP-seq data and regulatory SNPs

The transcription factor affinity prediction (TRAP) method calculates the affinity of transcription factors for DNA sequences on the basis of a biophysical model. This method has proven to be useful for several applications, including for determining the putative target genes of a given factor. This protocol covers two other applications: (i) determining which transcription factors have the highest affinity in a set of sequences (illustrated with chromatin immunoprecipitation-sequencing (ChIP-seq) peaks), and (ii) finding which factor is the most affected by a regulatory single-nucleotide polymorphism. The protocol describes how to use the TRAP web tools to address these questions, and it also presents a way to run TRAP on random control sequences to better estimate the significance of the results. All of the tools are fully available online and do not need any additional installation. The complete protocol takes about 45 min, but each individual tool runs in a few minutes.

Effects of single-nucleotide polymorphisms in the human holocarboxylase synthetase gene on enzyme catalysis

Holocarboxylase synthetase (HLCS) is a biotin protein ligase, which has a pivotal role in biotin-dependent metabolic pathways and epigenetic phenomena in humans. Knockdown of HLCS produces phenotypes such as heat susceptibility and decreased life span in Drosophila melanogaster, whereas knockout of HLCS appears to be embryonic lethal. HLCS comprises 726 amino acids in four domains. More than 2500 single-nucleotide polymorphisms (SNPs) have been identified in human HLCS. Here, we tested the hypotheses that HLCS SNPs impair enzyme activity, and that biotin supplementation restores the activities of HLCS variants to wild-type levels. We used an in silico approach to identify five SNPs that alter the amino acid sequence in the N-terminal, central, and C-terminal domains in human HLCS. Recombinant HLCS was used for enzyme kinetics analyses of HLCS variants, wild-type HLCS, and the L216R mutant, which has a biotin ligase activity near zero. The biotin affinity of variant Q699R is lower than that of the wild-type control, but the maximal activity was restored to that of wild-type HLCS when assay mixtures were supplemented with biotin. In contrast, the biotin affinities of HLCS variants V96F and G510R are not significantly different from the wild-type control, but their maximal activities remained moderately lower than that of wild-type HLCS even when assay mixtures were supplemented with biotin. The V96 L SNP did not alter enzyme kinetics. Our findings suggest that individuals with HLCS SNPs may benefit from supplemental biotin, yet to different extents depending on the genotype.

Cytosine methylation in miR-153 gene promoters increases the expression of holocarboxylase synthetase, thereby increasing the abundance of histone H4 biotinylation marks in HEK-293 human kidney cells

Holocarboxylase synthetase (HCS) plays an essential role in catalyzing the biotinylation of carboxylases and histones. Biotinylated carboxylases are important for the metabolism of glucose, lipids and leucine; biotinylation of histones plays important roles in gene regulation and genome stability. Recently, we reported that HCS activity is partly regulated by subcellular translocation events and by miR-539. Here we tested the hypothesis that the HCS 3'-untranslated region (3'-UTR) contains binding sites for miR other than miR-539. A binding site for miR-153 was predicted to reside in the HCS 3'-UTR by using in silico analyses. When miR-153 site was overexpressed in transgenic HEK-293 human embryonic kidney cells, the abundance of HCS mRNA decreased by 77% compared with controls. In silico analyses also predicted three putative cytosine methylation sites in two miR-153 genes; the existence of these sites was confirmed by methylation-sensitive polymerase chain reaction. When cytosines were demethylated by treatment with 5-aza-2'-deoxycytidine, the abundance of miR-153 increased by more than 25 times compared with untreated controls, and this increase coincided with low levels of HCS and histone biotinylation. Together, this study provides novel insights into the mechanisms of novel epigenetic synergies among folate-dependent methylation events, miR and histone biotinylation.

The role of holocarboxylase synthetase in genome stability is mediated partly by epigenomic synergies between methylation and biotinylation events

Holocarboxylase synthetase (HLCS) catalyzes the covalent binding of biotin to histones. Biotinylated histones are gene repression marks and are particularly enriched in long terminal repeats, telomeres, and other repeat regions. The effects of HLCS in gene regulation are mediated by its physical interactions with chromatin proteins such as histone H3, DNMT1, MeCP2, and EHMT-1. It appears that histone biotinylation depends on prior methylation of cytosines. De-repression of long terminal repeats in biotin- or HLCS-deficient cell cultures and organisms is associated with genome instability.

The BTB and CNC homology 1 (BACH1) target genes are involved in the oxidative stress response and in control of the cell cycle

The regulation of gene expression in response to environmental signals and metabolic imbalances is a key step in maintaining cellular homeostasis. BTB and CNC homology 1 (BACH1) is a heme-binding transcription factor repressing the transcription from a subset of MAF recognition elements at low intracellular heme levels. Upon heme binding, BACH1 is released from the MAF recognition elements, resulting in increased expression of antioxidant response genes. To systematically address the gene regulatory networks involving BACH1, we combined chromatin immunoprecipitation sequencing analysis of BACH1 target genes in HEK 293 cells with knockdown of BACH1 using three independent types of small interfering RNAs followed by transcriptome profiling using microarrays. The 59 BACH1 target genes identified by chromatin immunoprecipitation sequencing were found highly enriched in genes showing expression changes after BACH1 knockdown, demonstrating the impact of BACH1 repression on transcription. In addition to known and new BACH1 targets involved in heme degradation (HMOX1, FTL, FTH1, ME1, and SLC48A1) and redox regulation (GCLC, GCLM, and SLC7A11), we also discovered BACH1 target genes affecting cell cycle and apoptosis pathways (ITPR2, CALM1, SQSTM1, TFE3, EWSR1, CDK6, BCL2L11, and MAFG) as well as subcellular transport processes (CLSTN1, PSAP, MAPT, and vault RNA). The newly identified impact of BACH1 on genes involved in neurodegenerative processes and proliferation provides an interesting basis for future dissection of BACH1-mediated gene repression in neurodegeneration and virus-induced cancerogenesis.