JASPAR, the open access database of transcription factor-binding profiles: new content and tools in the 2008 update

MD-SVM: a novel SVM-based algorithm for the motif discovery of transcription factor binding sites

BACKGROUND

Transcription factors (TFs) play important roles in the regulation of gene expression. They can activate or block transcription of downstream genes in a manner of binding to specific genomic sequences. Therefore, motif discovery of these binding preference patterns is of central significance in the understanding of molecular regulation mechanism. Many algorithms have been proposed for the identification of transcription factor binding sites. However, it remains a challengeable problem.

RESULTS

Here, we proposed a novel motif discovery algorithm based on support vector machine (MD-SVM) to learn a discriminative model for TF binding sites. MD-SVM firstly obtains position weight matrix (PWM) from a set of training datasets. Then it translates the MD problem into a computational framework of multiple instance learning (MIL). It was applied to several real biological datasets. Results show that our algorithm outperforms MI-SVM in terms of both accuracy and specificity.

CONCLUSIONS

In this paper, we modeled the TF motif discovery problem as a MIL optimization problem. The SVM algorithm was adapted to discriminate positive and negative bags of instances. Compared to other svm-based algorithms, MD-SVM show its superiority over its competitors in term of ROC AUC. Hopefully, it could be of benefit to the research community in the understanding of molecular functions of DNA functional elements and transcription factors.

The miR206-JunD Circuit Mediates the Neurotoxic Effect of Methylmercury in Cortical Neurons

Methylmercury (MeHg) causes neuronal death through different pathways. Particularly, we found that in cortical neurons it increased the expression of Repressor Element-1 Silencing Transcription Factor (REST), histone deacetylase (HDAC)4, Specificity Protein (Sp)1, Sp4, and reduced the levels of brain-derived neurotrophic factor (BDNF). Herein, in rat cortical neurons we investigated whether microRNA (miR)206 can modulate MeHg-induced cell death by regulating REST/HDAC4/Sp1/Sp4/BDNF axis. MeHg (1 µM) reduced miR206 expression after both 12 and 24 h and miR206 transfection prevented MeHg-induced neuronal death. Furthermore, miR206 reverted MeHg-induced REST and Sp4 increase and BDNF reduction at gene and protein level, and reverted HDAC4 protein increase, but not HDAC4 mRNA upregulation. Moreover, since no miR206 seed sequences were identified in the 3'-untranslated regions (3'-UTRs) of REST and SP4, we investigated the role of JunD, that presents a consensus motif on REST, Sp4, and BDNF promoters. Indeed, MeHg increased JunD mRNA and protein levels, and JunD knockdown counteracted MeHg-induced REST, Sp4 increase, but not BDNF reduction. Furthermore, we identified a miR206 binding site in the 3'-UTR of JunD mRNA (miR206/JunD) and mutagenesis of miR206/JunD site reverted JunD luciferase activity reduction induced by miR206. Finally, miR206 prevented MeHg-increased JunD binding to REST and Sp4 promoters, and MeHg-reduced BDNF expression was determined by the increase of HDAC4 binding on BDNF promoter IV. Collectively, these results suggest that miR206 downregulation induced by MeHg exposure determines an upregulation of HDAC4, that in turn downregulated BDNF, and the activation of JunD that, by binding REST and Sp4 gene promoters, increased their expression.

HCP5 is a SMAD3-responsive long non-coding RNA that promotes lung adenocarcinoma metastasis via miR-203/SNAI axis

Introduction: Transforming growth factor-beta (TGFβ) signaling plays a vital role in lung adenocarcinoma (LUAD) progression. However, the involvement of TGFβ-regulated long non-coding RNAs (lncRNAs) in metastasis of LUAD remains poorly understood.

Methods: We performed bioinformatic analyses to identify putative lncRNAs regulated by TGF-β/SMAD3 and validated the results by quantitative PCR in LUAD cells. We performed luciferase reporter and chromatin immunoprecipitation assays to demonstrate the transcriptional regulation of the lncRNA histocompatibility leukocyte antigen complex P5 (HCP5) we decided to focus on. Stable HCP5 knockdown and HCP5-overexpressing A549 cell variants were generated respectively, to study HCP5 function and understand its mechanism of action. We also confirmed our findings in mouse xenografts and metastasis models. We analyzed the correlation between the level of lncRNA expression with EGFR, KRAS mutations, smoke state and prognostic of LUAD patients.

Results: We found that the lncRNA HCP5 is induced by TGFβ and transcriptionally regulated by SMAD3, which promotes LUAD tumor growth and metastasis. Moreover, HCP5 is overexpressed in tumor tissues of patients with LUAD, specifically in patients with EGFR and KRAS mutations and current smoker. HCP5 high expression level is positively correlated with poor prognosis of patients with LUAD. Finally, we demonstrated that upregulation of HCP5 increases the expression of Snail and Slug by sponging the microRNA-203 (miR-203) and promoting epithelial-mesenchymal transition (EMT) in LUAD cells.

Conclusions: Our work demonstrates that the lncRNA HCP5 is transcriptionally regulated by SMAD3 and acts as a new regulator in the TGFβ/SMAD signaling pathway. Therefore, HCP5 can serve as a potential therapeutic target in LUAD.

DNA-protein interaction: identification, prediction and data analysis

Life in living organisms is dependent on specific and purposeful interaction between other molecules. Such purposeful interactions make the various processes inside the cells and the bodies of living organisms possible. DNA-protein interactions, among all the types of interactions between different molecules, are of considerable importance. Currently, with the development of numerous experimental techniques, diverse methods are convenient for recognition and investigating such interactions. While the traditional experimental techniques to identify DNA-protein complexes are time-consuming and are unsuitable for genome-scale studies, the current high throughput approaches are more efficient in determining such interaction at a large-scale, but they are clearly too costly to be practice for daily applications. Hence, according to the availability of much information related to different biological sequences and clearing different dimensions of conditions in which such interactions are formed, with the developments related to the computer, mathematics, and statistics motivate scientists to develop bioinformatics tools for prediction the interaction site(s). Until now, there has been much progress in this field. In this review, the factors and conditions governing the interaction and the laboratory techniques for examining such interactions are addressed. In addition, developed bioinformatics tools are introduced and compared for this reason and, in the end, several suggestions are offered for the promotion of such tools in prediction with much more precision.

Integrated epigenomic profiling reveals endogenous retrovirus reactivation in renal cell carcinoma

BACKGROUND

Transcriptional dysregulation drives cancer formation but the underlying mechanisms are still poorly understood. Renal cell carcinoma (RCC) is the most common malignant kidney tumor which canonically activates the hypoxia-inducible transcription factor (HIF) pathway. Despite intensive study, novel therapeutic strategies to target RCC have been difficult to develop. Since the RCC epigenome is relatively understudied, we sought to elucidate key mechanisms underpinning the tumor phenotype and its clinical behavior.

METHODS

We performed genome-wide chromatin accessibility (DNase-seq) and transcriptome profiling (RNA-seq) on paired tumor/normal samples from 3 patients undergoing nephrectomy for removal of RCC. We incorporated publicly available data on HIF binding (ChIP-seq) in a RCC cell line. We performed integrated analyses of these high-resolution, genome-scale datasets together with larger transcriptomic data available through The Cancer Genome Atlas (TCGA).

FINDINGS

Though HIF transcription factors play a cardinal role in RCC oncogenesis, we found that numerous transcription factors with a RCC-selective expression pattern also demonstrated evidence of HIF binding near their gene body. Examination of chromatin accessibility profiles revealed that some of these transcription factors influenced the tumor's regulatory landscape, notably the stem cell transcription factor POU5F1 (OCT4). Elevated POU5F1 transcript levels were correlated with advanced tumor stage and poorer overall survival in RCC patients. Unexpectedly, we discovered a HIF-pathway-responsive promoter embedded within a endogenous retroviral long terminal repeat (LTR) element at the transcriptional start site of the PSOR1C3 long non-coding RNA gene upstream of POU5F1. RNA transcripts are induced from this promoter and read through PSOR1C3 into POU5F1 producing a novel POU5F1 transcript isoform. Rather than being unique to the POU5F1 locus, we found that HIF binds to several other transcriptionally active LTR elements genome-wide correlating with broad gene expression changes in RCC.

INTERPRETATION

Integrated transcriptomic and epigenomic analysis of matched tumor and normal tissues from even a small number of primary patient samples revealed remarkably convergent shared regulatory landscapes. Several transcription factors appear to act downstream of HIF including the potent stem cell transcription factor POU5F1. Dysregulated expression of POU5F1 is part of a larger pattern of gene expression changes in RCC that may be induced by HIF-dependent reactivation of dormant promoters embedded within endogenous retroviral LTRs.

Integrated Functional Genomic Analysis Enables Annotation of Kidney Genome-Wide Association Study Loci

BACKGROUND

Linking genetic risk loci identified by genome-wide association studies (GWAS) to their causal genes remains a major challenge. Disease-associated genetic variants are concentrated in regions containing regulatory DNA elements, such as promoters and enhancers. Although researchers have previously published DNA maps of these regulatory regions for kidney tubule cells and glomerular endothelial cells, maps for podocytes and mesangial cells have not been available.

METHODS

We generated regulatory DNA maps (DNase-seq) and paired gene expression profiles (RNA-seq) from primary outgrowth cultures of human glomeruli that were composed mainly of podocytes and mesangial cells. We generated similar datasets from renal cortex cultures, to compare with those of the glomerular cultures. Because regulatory DNA elements can act on target genes across large genomic distances, we also generated a chromatin conformation map from freshly isolated human glomeruli.

RESULTS

We identified thousands of unique regulatory DNA elements, many located close to transcription factor genes, which the glomerular and cortex samples expressed at different levels. We found that genetic variants associated with kidney diseases (GWAS) and kidney expression quantitative trait loci were enriched in regulatory DNA regions. By combining GWAS, epigenomic, and chromatin conformation data, we functionally annotated 46 kidney disease genes.

CONCLUSIONS

We demonstrate a powerful approach to functionally connect kidney disease-/trait-associated loci to their target genes by leveraging unique regulatory DNA maps and integrated epigenomic and genetic analysis. This process can be applied to other kidney cell types and will enhance our understanding of genome regulation and its effects on gene expression in kidney disease.

Genome-wide CRISPR Screens in T Helper Cells Reveal Pervasive Crosstalk between Activation and Differentiation

T helper type 2 (Th2) cells are important regulators of mammalian adaptive immunity and have relevance for infection, autoimmunity, and tumor immunology. Using a newly developed, genome-wide retroviral CRISPR knockout (KO) library, combined with RNA-seq, ATAC-seq, and ChIP-seq, we have dissected the regulatory circuitry governing activation and differentiation of these cells. Our experiments distinguish cell activation versus differentiation in a quantitative framework. We demonstrate that these two processes are tightly coupled and are jointly controlled by many transcription factors, metabolic genes, and cytokine/receptor pairs. There are only a small number of genes regulating differentiation without any role in activation. By combining biochemical and genetic data, we provide an atlas for Th2 differentiation, validating known regulators and identifying factors, such as Pparg and Bhlhe40, as part of the core regulatory network governing Th2 helper cell fates.

Identification of Novel Gata3 Distal Enhancers Active in Mouse Embryonic Lens

BACKGROUND

The tissue-specific transcriptional programs during normal development require tight control by distal cis-regulatory elements, such as enhancers, with specific DNA sequences recognized by transcription factors, coactivators, and chromatin remodeling enzymes. Gata3 is a sequence-specific DNA-binding transcription factor that regulates formation of multiple tissues and organs, including inner ear, lens, mammary gland, T-cells, urogenital system, and thyroid gland. In the eye, Gata3 has a highly restricted expression domain in the posterior part of the lens vesicle; however, the underlying regulatory mechanisms are unknown.

RESULTS

Here we describe the identification of a novel bipartite Gata3 lens-specific enhancer located ∼18 kb upstream from its transcriptional start site. We also found that a 5-kb Gata3 promoter possesses low activity in the lens. The bipartite enhancer contains arrays of AP-1, Ets-, and Smad1/5-binding sites as well as binding sites for lens-associated DNA-binding factors. Transient transfection studies of the promoter with the bipartite enhancer showed enhanced activation by BMP4 and FGF2.

CONCLUSIONS

These studies identify a novel distal enhancer of Gata3 with high activity in lens and indicate that BMP and FGF signaling can up-regulate expression of Gata3 in differentiating lens fiber cells through the identified Gata3 enhancer and promoter elements. Developmental Dynamics 247:1186-1198, 2018. © 2018 The Authors. Developmental Dynamics published by Wiley Periodicals, Inc. on behalf of American Association of Anatomists.

NFAT5 up-regulates expression of the kidney-specific ubiquitin ligase gene Rnf183 under hypertonic conditions in inner-medullary collecting duct cells

We previously reported that among the 37 RING finger protein (RNF) family members, RNF183 mRNA is specifically expressed in the kidney under normal conditions. However, the mechanism supporting its kidney-specific expression pattern remains unclear. In this study, we elucidated the mechanism of the transcriptional activation of murine Rnf183 in inner-medullary collecting duct cells. Experiments with anti-RNF183 antibody revealed that RNF183 is predominantly expressed in the renal medulla. Among the 37 RNF family members, Rnf183 mRNA expression was specifically increased in hypertonic conditions, a hallmark of the renal medulla. RNF183 up-regulation was consistent with the activation of nuclear factor of activated T cells 5 (NFAT5), a transcription factor essential for adaptation to hypertonic conditions. Accordingly, siRNA-mediated knockdown of NFAT5 down-regulated RNF183 expression. Furthermore, the -3,466 to -3,136-bp region upstream of the mouse Rnf183 promoter containing the NFAT5-binding motif is conserved among mammals. A luciferase-based reporter vector containing the NFAT5-binding site was activated in response to hypertonic stress, but was inhibited by a mutation at the NFAT5-binding site. ChIP assays revealed that the binding of NFAT5 to this DNA site is enhanced by hypertonic stress. Of note, siRNA-mediated RNF183 knockdown increased hypertonicity-induced caspase-3 activation and decreased viability of mIMCD-3 cells. These results indicate that (i) RNF183 is predominantly expressed in the normal renal medulla, (ii) NFAT5 stimulates transcriptional activation of Rnf183 by binding to its cognate binding motif in the Rnf183 promoter, and (iii) RNF183 protects renal medullary cells from hypertonicity-induced apoptosis.

Correlated Evolution of Two Copulatory Organs via a Single cis-Regulatory Nucleotide Change

Diverse traits often covary between species [1-3]. The possibility that a single mutation could contribute to the evolution of several characters between species [3] is rarely investigated as relatively few cases are dissected at the nucleotide level. Drosophila santomea has evolved additional sex comb sensory teeth on its legs and has lost two sensory bristles on its genitalia. We present evidence that a single nucleotide substitution in an enhancer of the scute gene contributes to both changes. The mutation alters a binding site for the Hox protein Abdominal-B in the developing genitalia, leading to bristle loss, and for another factor in the developing leg, leading to bristle gain. Our study suggests that morphological evolution between species can occur through a single nucleotide change affecting several sexually dimorphic traits. VIDEO ABSTRACT.

Epigenetic Study in Parkinson's Disease: A Pilot Analysis of DNA Methylation in Candidate Genes in Brain

Efforts have been made to understand the pathophysiology of Parkinson’s disease (PD). A significant number of studies have focused on genetics, despite the fact that the described pathogenic mutations have been observed only in around 10% of patients; this observation supports the fact that PD is a multifactorial disorder. Lately, differences in miRNA expression, histone modification, and DNA methylation levels have been described, highlighting the importance of epigenetic factors in PD etiology. Taking all this into consideration, we hypothesized that an alteration in the level of methylation in PD-related genes could be related to disease pathogenesis, possibly due to alterations in gene expression. After analysing promoter regions of five PD-related genes in three brain regions by pyrosequencing, we observed some differences in DNA methylation levels (hypo and hypermethylation) in substantia nigra in some CpG dinucleotides that, possibly through an alteration in Sp1 binding, could alter their expression.

Genetic determinants of childhood and adult height associated with osteosarcoma risk

BACKGROUND

Although increased height has been associated with osteosarcoma risk in previous epidemiologic studies, to the authors' knowledge the relative contribution of stature during different developmental timepoints remains unclear. Furthermore, the question of how genetic determinants of height impact osteosarcoma etiology remains unexplored. Genetic variants associated with stature in previous genome-wide association studies may be biomarkers of osteosarcoma risk.

METHODS

The authors tested the associations between osteosarcoma risk and polygenic scores for adult height (416 variants), childhood height (6 variants), and birth length (5 variants) in 864 osteosarcoma cases and 1879 controls of European ancestry.

RESULTS

Each standard deviation increase in the polygenic score for adult height, corresponding to a 1.7-cm increase in stature, was found to be associated with a 1.10-fold increase in the risk of osteosarcoma (95% confidence interval [95% CI], 1.01-1.19; P =.027). Each standard deviation increase in the polygenic score for childhood height, corresponding to a 0.5-cm increase in stature, was associated with a 1.10-fold increase in the risk of osteosarcoma (95% CI, 1.01-1.20; P =.023). The polygenic score for birth length was not found to be associated with osteosarcoma risk (P =.11). When adult and childhood height scores were modeled together, they were found to be independently associated with osteosarcoma risk (P =.037 and P = .043, respectively). An expression quantitative trait locus for cartilage intermediate layer protein 2 (CILP2), rs8103992, was significantly associated with osteosarcoma risk after adjustment for multiple comparisons (odds ratio, 1.35; 95% CI, 1.16-1.56 [P = 7.93×10-5 and Padjusted =.034]).

CONCLUSIONS

A genetic propensity for taller adult and childhood height attainments contributed independently to osteosarcoma risk in the current study data. These results suggest that the biological pathways affecting normal bone growth may be involved in osteosarcoma etiology.

A coding and non-coding transcriptomic perspective on the genomics of human metabolic disease

Genome-wide association studies (GWAS), relying on hundreds of thousands of individuals, have revealed >200 genomic loci linked to metabolic disease (MD). Loss of insulin sensitivity (IS) is a key component of MD and we hypothesized that discovery of a robust IS transcriptome would help reveal the underlying genomic structure of MD. Using 1,012 human skeletal muscle samples, detailed physiology and a tissue-optimized approach for the quantification of coding (>18,000) and non-coding (>15,000) RNA (ncRNA), we identified 332 fasting IS-related genes (CORE-IS). Over 200 had a proven role in the biochemistry of insulin and/or metabolism or were located at GWAS MD loci. Over 50% of the CORE-IS genes responded to clinical treatment; 16 quantitatively tracking changes in IS across four independent studies (P = 0.0000053: negatively: AGL, G0S2, KPNA2, PGM2, RND3 and TSPAN9 and positively: ALDH6A1, DHTKD1, ECHDC3, MCCC1, OARD1, PCYT2, PRRX1, SGCG, SLC43A1 and SMIM8). A network of ncRNA positively related to IS and interacted with RNA coding for viral response proteins (P < 1 × 10-48), while reduced amino acid catabolic gene expression occurred without a change in expression of oxidative-phosphorylation genes. We illustrate that combining in-depth physiological phenotyping with robust RNA profiling methods, identifies molecular networks which are highly consistent with the genetics and biochemistry of human metabolic disease.

YY1 binds to the E3' enhancer and inhibits the expression of the immunoglobulin κ gene via epigenetic modifications

The rearrangement and expression of immunoglobulin genes are regulated by enhancers and their binding transcriptional factors that activate or suppress the activities of the enhancers. The immunoglobulin κ (Igκ) gene locus has three important enhancers: the intrinsic enhancer (Ei), 3' enhancer (E3'), and distal enhancer (Ed). Ei and E3' are both required for Igκ gene rearrangement during early stages of B-cell development, whereas optimal expression of the rearranged Igκ gene relies on both E3' and Ed. The transcription factor YY1 affects the expression of many genes involved in B-cell development, probably by mediating interactions between their enhancers and promoters. Herein, we found that YY1 binds to the E3' enhancer and suppresses Igκ expression in B lymphoma cells by epigenetically modifying the enhancer. Knocking down YY1 enhanced Igκ expression, which was associated with increased levels of E2A (encoded by the TCF3 gene) and its binding to the E3' enhancer. Moreover, in germinal centre B cells and plasma cells, YY1 expression was reversely associated with Igκ levels, implying that YY1 might facilitate antibody affinity maturation in germinal centre B cells through the transient attenuation of Igκ expression.

DIVERSITY in binding, regulation, and evolution revealed from high-throughput ChIP

Genome-wide in vivo protein-DNA interactions are routinely mapped using high-throughput chromatin immunoprecipitation (ChIP). ChIP-reported regions are typically investigated for enriched sequence-motifs, which are likely to model the DNA-binding specificity of the profiled protein and/or of co-occurring proteins. However, simple enrichment analyses can miss insights into the binding-activity of the protein. Note that ChIP reports regions making direct contact with the protein as well as those binding through intermediaries. For example, consider a ChIP experiment targeting protein X, which binds DNA at its cognate sites, but simultaneously interacts with four other proteins. Each of these proteins also binds to its own specific cognate sites along distant parts of the genome, a scenario consistent with the current view of transcriptional hubs and chromatin loops. Since ChIP will pull down all X-associated regions, the final reported data will be a union of five distinct sets of regions, each containing binding sites of one of the five proteins, respectively. Characterizing all five different motifs and the corresponding sets is important to interpret the ChIP experiment and ultimately, the role of X in regulation. We present diversity which attempts exactly this: it partitions the data so that each partition can be characterized with its own de novo motif. Diversity uses a Bayesian approach to identify the optimal number of motifs and the associated partitions, which together explain the entire dataset. This is in contrast to standard motif finders, which report motifs individually enriched in the data, but do not necessarily explain all reported regions. We show that the different motifs and associated regions identified by diversity give insights into the various complexes that may be forming along the chromatin, something that has so far not been attempted from ChIP data. Webserver at http://diversity.ncl.res.in/; standalone (Mac OS X/Linux) from https://github.com/NarlikarLab/DIVERSITY/releases/tag/v1.0.0.

A high-throughput chromatin immunoprecipitation (ChIP) experiment identifies genomic regions bound by a protein in vivo. Current motif-discovery approaches seek an enriched motif signature in the reported regions, which they can attribute to the protein’s binding preferences. However, Diversity models the fact that since a ChIP experiment pulls down regions participating in all complexes involving the profiled protein, the reported regions are in all likelihood, a collection of different types of protein-DNA contacts. Diversity asks a different question: what sequence component caused a specific region to be reported in a ChIP experiment? The answer, in combination with additional data such as sequence conservation, SNPs, chromatin structure, downstream gene-expression, etc. can yield insights into the diverse regulatory mechanisms at play. The added benefits of a webserver and a standalone parallel version make diversity a practical tool for discovering new biology from ChIP experiments.

Selection on the regulation of sympathetic nervous activity in humans and chimpanzees

Adrenergic α2C receptor (ADRA2C) is an inhibitory modulator of the sympathetic nervous system. Knockout mice for this gene show physiological and behavioural alterations that are associated with the fight-or-flight response. There is evidence of positive selection on the regulation of this gene during chicken domestication. Here, we find that the neuronal expression of ADRA2C is lower in human and chimpanzee than in other primates. On the basis of three-dimensional chromatin structure, we identified a cis-regulatory region whose DNA sequences have been significantly accelerated in human and chimpanzee. Active histone modification marks this region in rhesus macaque but not in human and chimpanzee; instead, repressive marks are enriched in various human brain samples. This region contains two neuron-restrictive silencer factor (NRSF) binding motifs, each of which harbours a polymorphism. Our genotyping and analysis of population genome data indicate that at both polymorphic sites, the derived allele has reached fixation in humans and chimpanzees but not in bonobos, whereas only the ancestral allele is present among macaques. Our CRISPR/Cas9 genome editing and reporter assays show that both derived nucleotides repress ADRA2C, most likely by increasing NRSF binding. In addition, we detected signatures of recent positive selection for lower neuronal ADRA2C expression in humans. Our findings indicate that there has been selective pressure for enhanced sympathetic nervous activity in the evolution of humans and chimpanzees.

The HNF1α-regulated lncRNA HNF1A-AS1 reverses the malignancy of hepatocellular carcinoma by enhancing the phosphatase activity of SHP-1

BACKGROUND

Our previous study has demonstrated that hepatocyte nuclear factor 1α (HNF1α) exerts potent therapeutic effects on hepatocellular carcinoma (HCC). However, the molecular mechanisms by which HNF1α reverses HCC malignancy need to be further elucidated.

METHODS

lncRNA microarray was performed to identify the long noncoding RNAs (lncRNAs) regulated by HNF1α. Chromatin immunoprecipitation and luciferase reporter assays were applied to clarify the mechanism of the transcriptional regulation of HNF1α to HNF1A antisense RNA 1 (HNF1A-AS1). The effect of HNF1A-AS1 on HCC malignancy was evaluated in vitro and in vivo. RNA pulldown, RNA-binding protein immunoprecipitation and the Bio-Layer Interferometry assay were used to validate the interaction of HNF1A-AS1 and Src homology region 2 domain-containing phosphatase 1 (SHP-1).

RESULTS

HNF1α regulated the expression of a subset of lncRNAs in HCC cells. Among these lncRNAs, the expression levels of HNF1A-AS1 were notably correlated with HNF1α levels in HCC cells and human HCC tissues. HNF1α activated the transcription of HNF1A-AS1 by directly binding to its promoter region. HNF1A-AS1 inhibited the growth and the metastasis of HCC cells in vitro and in vivo. Moreover, knockdown of HNF1A-AS1 reversed the suppressive effects of HNF1α on the migration and invasion of HCC cells. Importantly, HNF1A-AS1 directly bound to the C-terminal of SHP-1 with a high binding affinity (KD = 59.57 ± 14.29 nM) and increased the phosphatase activity of SHP-1. Inhibition of SHP-1 enzymatic activity substantially reversed the HNF1α- or HNF1A-AS1-induced reduction on the metastatic property of HCC cells.

CONCLUSIONS

Our data revealed that HNF1A-AS1 is a direct transactivation target of HNF1α in HCC cells and involved in the anti-HCC effect of HNF1α. HNF1A-AS1 functions as phosphatase activator through the direct interaction with SHP-1. These findings suggest that regulation of the HNF1α/HNF1A-AS1/SHP-1 axis may have beneficial effects in the treatment of HCC.

Local sequence features that influence AP-1 cis-regulatory activity

In the genome, most occurrences of transcription factor binding sites (TFBS) have no cis-regulatory activity, which suggests that flanking sequences contain information that distinguishes functional from nonfunctional TFBS. We interrogated the role of flanking sequences near Activator Protein 1 (AP-1) binding sites that reside in DNase I Hypersensitive Sites (DHS) and regions annotated as Enhancers. In these regions, we found that sequence features directly adjacent to the core motif distinguish high from low activity AP-1 sites. Some nearby features are motifs for other TFs that genetically interact with the AP-1 site. Other features are extensions of the AP-1 core motif, which cause the extended sites to match motifs of multiple AP-1 binding proteins. Computational models trained on these data distinguish between sequences with high and low activity AP-1 sites and also predict changes in cis-regulatory activity due to mutations in AP-1 core sites and their flanking sequences. Our results suggest that extended AP-1 binding sites, together with adjacent binding sites for additional TFs, encode part of the information that governs TFBS activity in the genome.

The Ensembl Genome Browser: Strategies for Accessing Eukaryotic Genome Data

The Ensembl Genome Browser provides a wealth of freely available genomic data that can be accessed for many purposes by genetics, genomics, and molecular biology researchers. Herein we present two protocols for exploring different aspects of these data: a phenotype and its associated variants and genes, and a promoter and the epigenetic marks and protein-binding activity associated with it. These workflows illustrate a subset of the data types available through the Ensembl Browser, and can be considered a springboard for further exploration.

Apoptosis signal-regulating kinase 1 mediates the inhibitory effect of hepatocyte nuclear factor-4α on hepatocellular carcinoma

Previous studies provided substantial evidence of a striking suppressive effect of hepatocyte nuclear factor 4α (HNF4α) on hepatocellular carcinoma (HCC). Apoptosis signal-regulating kinase 1 (ASK1) is involved in death receptor-mediated apoptosis and may acts as a tumor suppressor in hepatocarcinogenesis. However, the status and function of ASK1 during HCC progression are unclear. In this study, we found that HNF4α increased ASK1 expression by directly binding to its promoter. ASK1 expression was dramatically suppressed and correlated with HNF4α levels in HCC tissues. Reduced ASK1 expression was associated with aggressive tumors and poor prognosis for human HCC. Moreover, ASK1 inhibited the malignant phenotype of HCC cells in vitro. Intratumoral ASK1 injection significantly suppressed the growth of subcutaneous HCC xenografts in nude mice. More interestingly, systemic ASK1 delivery strikingly inhibited the growth of orthotopic HCC nodules in NOD/SCID mice. In addition, inhibition of endogenous ASK1 partially reversed the suppressive effects of HNF4α on HCC. Collectively, this study highlights the suppressive effect of ASK1 on HCC and its biological significance in HCC development. These outcomes broaden the knowledge of ASK1 function in HCC progression, and provide a novel potential prognostic biomarker and therapeutic target for advanced HCC.

Eukaryotic and prokaryotic promoter databases as valuable tools in exploring the regulation of gene transcription: a comprehensive overview

The complete exploration of the regulation of gene expression remains one of the top-priority goals for researchers. As the regulation is mainly controlled at the level of transcription by promoters, study on promoters and findings are of great importance. This review summarizes forty selected databases that centralize experimental and theoretical knowledge regarding the organization of promoters, interacting transcription factors (TFs) and microRNAs (miRNAs) in many eukaryotic and prokaryotic species. The presented databases offer researchers valuable support in elucidating the regulation of gene transcription.

DnrI of Streptomyces peucetius binds to the resistance genes, drrAB and drrC but is activated by daunorubicin

The master regulator, DnrI of Streptomyces peucetius is a member of the family of transcriptional activator, Streptomyces antibiotic regulatory proteins (SARP), which controls the biosynthesis of antitumor anthracycline, daunorubicin (DNR) and doxorubicin (DXR). The binding of DnrI to the heptameric repeat sequence found within the -35 promoter region of biosynthetic gene, dpsE activates it. To combat the increased level of intracellular DNR, the cell has developed self resistance mechanism mediated by drrAB and drrC genes which are regulated by regulatory genes. We find that a drug non-producing mutant, ΔdpsA, showed sensitive phenotype in plate assay along with an increased level of dnrI transcript. Whereas the mutant grown in the presence of DNR showed a resistant phenotype with a six and eight folds increase in drrAB and drrC transcripts respectively. Computational studies followed by molecular docking showed that DnrI bound as a monomer to a slightly modified heptameric DNA motif, 5'-ACACGCA in drrA and 5'-ACAACCT in drrC which was also proved by electrophoretic mobility shift assay. These findings confirm that DnrI belongs to winged helix-turn-helix DNA-binding protein with Tetratricopeptide Repeat domain. The transcriptional regulator DnrI binds to the resistance genes at specific sites but they are activated only when an increased load of intracellular DNR is sensed.

Ubiquitin-specific peptidase 7 (USP7)-mediated deubiquitination of the histone deacetylase SIRT7 regulates gluconeogenesis

Sirtuin 7 (SIRT7), a member of the NAD⁺-dependent class III histone deacetylases, is involved in the regulation of various cellular processes and in resisting various stresses, such as hypoxia, low glucose levels, and DNA damage. Interestingly, SIRT7 is linked to the control of glycolysis, suggesting a role in glucose metabolism. Given the important roles of SIRT7, it is critical to clarify how SIRT7 activity is potentially regulated. It has been reported that some transcriptional and post-transcriptional regulatory mechanisms are involved. However, little is known how SIRT7 is regulated by the post-translational modifications. Here, we identified ubiquitin-specific peptidase 7 (USP7), a deubiquitinase, as a negative regulator of SIRT7. We showed that USP7 interacts with SIRT7 both in vitro and in vivo, and we further demonstrated that SIRT7 undergoes endogenous Lys-63-linked polyubiquitination, which is removed by USP7. Although the USP7-mediated deubiquitination of SIRT7 had no effect on its stability, the deubiquitination repressed its enzymatic activity. We also showed that USP7 coordinates with SIRT7 to regulate the expression of glucose-6-phosphatase catalytic subunit (G6PC), a gluconeogenic gene. USP7 depletion by RNA interference increased both G6PC expression and SIRT7 enzymatic activity. Moreover, SIRT7 targeted the G6PC promoter through the transcription factor ELK4 but not through forkhead box O1 (FoxO1). In summary, SIRT7 is a USP7 substrate and has a novel role as a regulator of gluconeogenesis. Our study may provide the basis for new clinical approaches to treat metabolic disorders related to glucose metabolism.

Sox5 regulates beta-cell phenotype and is reduced in type 2 diabetes

Type 2 diabetes (T2D) is characterized by insulin resistance and impaired insulin secretion, but the mechanisms underlying insulin secretion failure are not completely understood. Here, we show that a set of co-expressed genes, which is enriched for genes with islet-selective open chromatin, is associated with T2D. These genes are perturbed in T2D and have a similar expression pattern to that of dedifferentiated islets. We identify Sox5 as a regulator of the module. Sox5 knockdown induces gene expression changes similar to those observed in T2D and diabetic animals and has profound effects on insulin secretion, including reduced depolarization-evoked Ca²⁺-influx and β-cell exocytosis. SOX5 overexpression reverses the expression perturbations observed in a mouse model of T2D, increases the expression of key β-cell genes and improves glucose-stimulated insulin secretion in human islets from donors with T2D. We suggest that human islets in T2D display changes reminiscent of dedifferentiation and highlight SOX5 as a regulator of β-cell phenotype and function.

Canonical and alternative transcript expression of PAX6 and CXCR4 in pancreatic cancer

Pancreatic cancer is a lethal disease with a propensity for invading and metastasizing into the surrounding tissues, including the liver and intestines. A number of factors are aberrantly overexpressed in this tumor type and actively promote cancer progression and metastasis. The present study demonstrates that paired box transcription factor 6 (PAX6) and C-X-C chemokine receptor 4 (CXCR4) are frequently co-expressed in primary pancreatic adenocarcinoma tumors and established cell lines. Expression analysis methods used in the present study included evaluation of protein expression by western blot analysis and immunofluorescence, transcript expression levels by reverse transcription-quantitative polymerase chain reaction (RT-qPCR), and luciferase assays utilizing regulatory elements from the CXCR4 gene locus. Canonical PAX6 and alternative splice variant PAX6(5a) proteins are expressed in pancreatic cancer and can drive gene expression through a conserved enhancer element within the first intron of the CXCR4 gene. As demonstrated by the introduction of an exogenous reporter construct with or without the intronic enhancer, loss of this element inhibited gene expression within numerous pancreatic cancer cell lines including Panc1, MIA-PaCa2 and BxPC3. All of the pancreatic cancer cell lines expressed the canonical CXCR4B transcript in addition to the alternatively spliced variant CXCR4A as determined by RT-qPCR experiments. The discovery of variant transcripts in pancreatic cancer cells may provide new candidates for future targeted therapies.

Epigenomic annotation of noncoding mutations identifies mutated pathways in primary liver cancer

Evidence that noncoding mutation can result in cancer driver events is mounting. However, it is more difficult to assign molecular biological consequences to noncoding mutations than to coding mutations, and a typical cancer genome contains many more noncoding mutations than protein-coding mutations. Accordingly, parsing functional noncoding mutation signal from noise remains an important challenge. Here we use an empirical approach to identify putatively functional noncoding somatic single nucleotide variants (SNVs) from liver cancer genomes. Annotation of candidate variants by publicly available epigenome datasets finds that 40.5% of SNVs fall in regulatory elements. When assigned to specific regulatory elements, we find that the distribution of regulatory element mutation mirrors that of nonsynonymous coding mutation, where few regulatory elements are recurrently mutated in a patient population but many are singly mutated. We find potential gain-of-binding site events among candidate SNVs, suggesting a mechanism of action for these variants. When aggregating noncoding somatic mutation in promoters, we find that genes in the ERBB signaling and MAPK signaling pathways are significantly enriched for promoter mutations. Altogether, our results suggest that functional somatic SNVs in cancer are sporadic, but occasionally occur in regulatory elements and may affect phenotype by creating binding sites for transcriptional regulators. Accordingly, we propose that noncoding mutation should be formally accounted for when determining gene- and pathway-mutation burden in cancer.

A scalable algorithm for structure identification of complex gene regulatory network from temporal expression data

Background

Gene regulatory interactions are of fundamental importance to various biological functions and processes. However, only a few previous computational studies have claimed success in revealing genome-wide regulatory landscapes from temporal gene expression data, especially for complex eukaryotes like human. Moreover, recent work suggests that these methods still suffer from the curse of dimensionality if a network size increases to 100 or higher.

Results

Here we present a novel scalable algorithm for identifying genome-wide gene regulatory network (GRN) structures, and we have verified the algorithm performances by extensive simulation studies based on the DREAM challenge benchmark data. The highlight of our method is that its superior performance does not degenerate even for a network size on the order of 10⁴, and is thus readily applicable to large-scale complex networks. Such a breakthrough is achieved by considering both prior biological knowledge and multiple topological properties (i.e., sparsity and hub gene structure) of complex networks in the regularized formulation. We also validate and illustrate the application of our algorithm in practice using the time-course gene expression data from a study on human respiratory epithelial cells in response to influenza A virus (IAV) infection, as well as the CHIP-seq data from ENCODE on transcription factor (TF) and target gene interactions. An interesting finding, owing to the proposed algorithm, is that the biggest hub structures (e.g., top ten) in the GRN all center at some transcription factors in the context of epithelial cell infection by IAV.

Conclusions

The proposed algorithm is the first scalable method for large complex network structure identification. The GRN structure identified by our algorithm could reveal possible biological links and help researchers to choose which gene functions to investigate in a biological event. The algorithm described in this article is implemented in MATLAB ^Ⓡ, and the source code is freely available from https://github.com/Hongyu-Miao/DMI.git.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1489-z) contains supplementary material, which is available to authorized users.

Regulatory elements in molecular networks

Regulatory elements determine the connectivity of molecular networks and mediate a variety of regulatory processes ranging from DNA looping to transcriptional, posttranscriptional, and posttranslational regulation. This review highlights our current understanding of the different types of regulatory elements found in molecular networks with a focus on DNA regulatory elements. We highlight technical advances and current challenges for the mapping of regulatory elements at the genome-wide scale, and describe new computational methods to uncover these elements via reconstructing regulatory networks from large genomic datasets. WIREs Syst Biol Med 2017, 9:e1374. doi: 10.1002/wsbm.1374 For further resources related to this article, please visit the WIREs website.

Genome-wide identification of regulatory elements in Sertoli cells

A current goal of molecular biology is to identify transcriptional networks that regulate cell differentiation. However, identifying functional gene regulatory elements has been challenging in the context of developing tissues where material is limited and cell types are mixed. To identify regulatory sites during sex determination, we subjected Sertoli cells from mouse fetal testes to DNaseI-seq and ChIP-seq for H3K27ac. DNaseI-seq identified putative regulatory sites around genes enriched in Sertoli and pregranulosa cells; however, active enhancers marked by H3K27ac were enriched proximal to only Sertoli-enriched genes. Sequence analysis identified putative binding sites of known and novel transcription factors likely controlling Sertoli cell differentiation. As a validation of this approach, we identified a novel Sertoli cell enhancer upstream of Wt1, and used it to drive expression of a transgenic reporter in Sertoli cells. This work furthers our understanding of the complex genetic network that underlies sex determination and identifies regions that potentially harbor non-coding mutations underlying disorders of sexual development.

Cell-lineage specificity and role of AP-1 in the prostate fibroblast androgen receptor cistrome

Androgen receptor (AR) signalling in fibroblasts is important in prostate development and carcinogenesis, and is inversely related to prostate cancer mortality. However, the molecular mechanisms of AR action in fibroblasts and other non-epithelial cell types are largely unknown. The genome-wide DNA binding profile of AR in human prostate fibroblasts was identified by chromatin immunoprecipitation sequencing (ChIP-Seq), and found to be common to other fibroblast lines but disparate from AR cistromes of prostate cancer cells and tissue. Although AR binding sites specific to fibroblasts were less well conserved evolutionarily than those shared with cancer epithelia, they were likewise correlated with androgen regulation of fibroblast gene expression. Whereas FOXA1 is the key pioneer factor of AR in cancer epithelia, our data indicated that AP-1 likely plays a more important role in the AR cistrome in fibroblasts. The specificity of AP-1 and FOXA1 to binding in these cells is demonstrated using immunoblot and immunohistochemistry. Importantly, we find the fibroblast cistrome is represented in whole tissue/in vivo ChIP-seq studies at both genomic and resulting protein levels, highlighting the importance of the stroma in whole tissue -omic studies. This is the first nuclear receptor ChIP-seq study in prostatic fibroblasts, and provides novel insight into the action of fibroblast AR in prostate cancer.

The pan-cancer pathological regulatory landscape

Dysregulation of the normal gene expression program is the cause of a broad range of diseases, including cancer. Detecting the specific perturbed regulators that have an effect on the generation and the development of the disease is crucial for understanding the disease mechanism and for taking decisions on efficient preventive and curative therapies. Moreover, detecting such perturbations at the patient level is even more important from the perspective of personalized medicine. We applied the Transcription Factor Target Enrichment Analysis, a method that detects the activity of transcription factors based on the quantification of the collective transcriptional activation of their targets, to a large collection of 5607 cancer samples covering eleven cancer types. We produced for the first time a comprehensive catalogue of altered transcription factor activities in cancer, a considerable number of them significantly associated to patient’s survival. Moreover, we described several interesting TFs whose activity do not change substantially in the cancer with respect to the normal tissue but ultimately play an important role in patient prognostic determination, which suggest they might be promising therapeutic targets. An additional advantage of this method is that it allows obtaining personalized TF activity estimations for individual patients.

Most of the tight positional conservation of transcription factor binding sites near the transcription start site reflects their co-localization within regulatory modules

Background

Transcription factors (TFs) form complexes that bind regulatory modules (RMs) within DNA, to control specific sets of genes. Some transcription factor binding sites (TFBSs) near the transcription start site (TSS) display tight positional preferences relative to the TSS. Furthermore, near the TSS, RMs can co-localize TFBSs with each other and the TSS. The proportion of TFBS positional preferences due to TFBS co-localization within RMs is unknown, however. ChIP experiments confirm co-localization of some TFBSs genome-wide, including near the TSS, but they typically examine only a few TFs at a time, using non-physiological conditions that can vary from lab to lab. In contrast, sequence analysis can examine many TFs uniformly and methodically, broadly surveying the co-localization of TFBSs with tight positional preferences relative to the TSS.

Results

Our statistics found 43 significant sets of human motifs in the JASPAR TF Database with positional preferences relative to the TSS, with 38 preferences tight (±5 bp). Each set of motifs corresponded to a gene group of 135 to 3304 genes, with 42/43 (98%) gene groups independently validated by DAVID, a gene ontology database, with FDR < 0.05. Motifs corresponding to two TFBSs in a RM should co-occur more than by chance alone, enriching the intersection of the gene groups corresponding to the two TFs. Thus, a gene-group intersection systematically enriched beyond chance alone provides evidence that the two TFs participate in an RM. Of the 903 = 43*42/2 intersections of the 43 significant gene groups, we found 768/903 (85%) pairs of gene groups with significantly enriched intersections, with 564/768 (73%) intersections independently validated by DAVID with FDR < 0.05. A user-friendly web site at http://go.usa.gov/3kjsH permits biologists to explore the interaction network of our TFBSs to identify candidate subunit RMs.

Conclusions

Gene duplication and convergent evolution within a genome provide obvious biological mechanisms for replicating an RM near the TSS that binds a particular TF subunit. Of all intersections of our 43 significant gene groups, 85% were significantly enriched, with 73% of the significant enrichments independently validated by gene ontology. The co-localization of TFBSs within RMs therefore likely explains much of the tight TFBS positional preferences near the TSS.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1354-5) contains supplementary material, which is available to authorized users.

Identifying changes in punitive transcriptional factor binding sites from regulatory single nucleotide polymorphisms that are significantly associated with disease or sickness

AIM

To identify punitive transcriptional factor binding sites (TFBS) from regulatory single nucleotide polymorphisms (rSNPs) that are significantly associated with disease.

METHODS

The genome-wide association studies have provided us with nearly 6500 disease or trait-predisposing SNPs where 93% are located within non-coding regions such as gene regulatory or intergenic areas of the genome. In the regulatory region of a gene, a SNP can change the DNA sequence of a transcriptional factor (TF) motif and in turn may affect the process of gene regulation. SNP changes that affect gene expression and impact gene regulatory sequences such as promoters, enhancers, and silencers are known as rSNPs. Computational tools can be used to identify unique punitive TFBS created by rSNPs that are associated with disease or sickness. Computational analysis was used to identify punitive TFBS generated by the alleles of these rSNPs.

RESULTS

rSNPs within nine genes that have been significantly associated with disease or sickness were used to illustrate the tremendous diversity of punitive unique TFBS that can be generated by their alleles. The genes studied are the adrenergic, beta, receptor kinase 1, the v-akt murine thymoma viral oncogene homolog 3, the activating transcription factor 3, the type 2 demodkinase gene, the endothetal Per-Arnt-Sim domain protein 1, the lysosomal acid lipase A, the signal Transducer and Activator of Transcription 4, the thromboxane A2 receptor and the vascular endothelial growth factor A. From this sampling of SNPs among the nine genes, there are 73 potential unique TFBS generated by the common alleles compared to 124 generated by the minor alleles indicating the tremendous diversity of potential TFs that are capable of regulating these genes.

CONCLUSION

From the diversity of unique punitive binding sites for TFs, it was found that some TFs play a role in the disease or sickness being studied.

Understanding the combinatorial action of transcription factors and microRNA regulation from regions of open chromatin

Transcriptional regulatory cascades are always triggered through the combinatorial interplay between transcription factors (TFs) and microRNAs (miRNAs) in eukaryotes. However, it is still a very substantial undertaking to dynamically profile their coordinated actions. In this work, we compared the differences in TFBS numbers between miRNA targets and non-targets, and found that miRNA targets tend to have more TFBS numbers. With the attempt to comprehensively understand the combinatorial action of TF and miRNA regulation from regions of open chromatin, we retrieved recently published DNase I hypersensitive sites (DHSs) across different human cell lines. The result showed that the differences are more statistically significant in DHS regions than non-DHS regions. Next, we trained classifiers for miRNA targets and non-targets. The result showed that TFBSs located in DHS regions achieved a competitive performance when discriminating miRNA targets and non-targets, whereas the performance of classifiers using TFBSs located in non-DHS regions is close to that of a random classifier. After the DHSs were divided into intergenic, transcription start sites (TSSs) and gene body DHS regions based on their genomic locations, only TFBSs located in TSS DHS regions provided a competitive performance. Our results provide us a clue that the coordinated activity of miRNAs and TFs describing the mechanism of gene expression control should be examined in a dynamic perspective.

Homeodomain proteins in action: similar DNA binding preferences, highly variable connectivity

Homeodomain proteins are evolutionary conserved proteins present in the entire eukaryote kingdom. They execute functions that are essential for life, both in developing and adult organisms. Most homeodomain proteins act as transcription factors and bind DNA to control the activity of other genes. In contrast to their similar DNA binding specificity, homeodomain proteins execute highly diverse and context-dependent functions. Several factors, including genome accessibility, DNA shape, combinatorial binding and the ability to interact with many transcriptional partners, diversify the activity of homeodomain proteins and culminate in the activation of highly dynamic, context-specific transcriptional programs. Clarifying how homeodomain transcription factors work is central to our understanding of development, disease and evolution.

Methodology for single nucleotide polymorphism selection in promoter regions for clinical use. An example of its applicability

Genetic variability in humans can explain many differences in disease risk factors. Polymorphism-related studies focus mainly on the single nucleotide polymorphisms (SNPs) of coding regions of the genes. SNPs on DNA binding motifs of the promoter region have been less explored. On a recent study of SNPs in patients with non-Hodgkin lymphomas we faced the problem of SNP selection from promoter regions and developed a practical methodology for clinical studies. The process consists in identifying SNPs in the coding and promoter regions of the antigen-processing system using the 'dbSNP' database. With the 'HapMap' program, we select SNPs with frequencies >20% in Caucasian populations. For coding regions, we sought biologically and clinically relevant SNPs described in the literature. For the promoter regions, we determined their chromosomal location on 'QiagenSABioscience' site database. The nucleotide sequence of ancestral and variant alleles is available in the 'dbSNP'. These sequences were used in 'Promoter TESS' to determine binding differences of transcription factors. Each sequence may have affinity to different TFs. Thus, SNP selection on the promoter regions was based in the differences on TF binding pattern between the old and the new allele. The potential clinical relevance of the new TFs was also evaluated before the final selection. With this approach, we found that almost half of the relevant SNP fall within the promoter region. In conclusion, we were able to develop a methodology of oriented selection of promoter regions of human genes, comparing the TF with affinity to the ancestral allele with the TF to a variant allele. We selected those SNPs that change the TF's affinity to a pattern with functional significance.

Full-Length Isoform Sequencing Reveals Novel Transcripts and Substantial Transcriptional Overlaps in a Herpesvirus

Whole transcriptome studies have become essential for understanding the complexity of genetic regulation. However, the conventionally applied short-read sequencing platforms cannot be used to reliably distinguish between many transcript isoforms. The Pacific Biosciences (PacBio) RS II platform is capable of reading long nucleic acid stretches in a single sequencing run. The pseudorabies virus (PRV) is an excellent system to study herpesvirus gene expression and potential interactions between the transcriptional units. In this work, non-amplified and amplified isoform sequencing protocols were used to characterize the poly(A⁺) fraction of the lytic transcriptome of PRV, with the aim of a complete transcriptional annotation of the viral genes. The analyses revealed a previously unrecognized complexity of the PRV transcriptome including the discovery of novel protein-coding and non-coding genes, novel mono- and polycistronic transcription units, as well as extensive transcriptional overlaps between neighboring and distal genes. This study identified non-coding transcripts overlapping all three replication origins of the PRV, which might play a role in the control of DNA synthesis. We additionally established the relative expression levels of gene products. Our investigations revealed that the whole PRV genome is utilized for transcription, including both DNA strands in all coding and intergenic regions. The genome-wide occurrence of transcript overlaps suggests a crosstalk between genes through a network formed by interacting transcriptional machineries with a potential function in the control of gene expression.

Enhancer Variants Synergistically Drive Dysfunction of a Gene Regulatory Network In Hirschsprung Disease

Common sequence variants in cis-regulatory elements (CREs) are suspected etiological causes of complex disorders. We previously identified an intronic enhancer variant in the RET gene disrupting SOX10 binding and increasing Hirschsprung disease (HSCR) risk 4-fold. We now show that two other functionally independent CRE variants, one binding Gata2 and the other binding Rarb, also reduce Ret expression and increase risk 2- and 1.7-fold. By studying human and mouse fetal gut tissues and cell lines, we demonstrate that reduced RET expression propagates throughout its gene regulatory network, exerting effects on both its positive and negative feedback components. We also provide evidence that the presence of a combination of CRE variants synergistically reduces RET expression and its effects throughout the GRN. These studies show how the effects of functionally independent non-coding variants in a coordinated gene regulatory network amplify their individually small effects, providing a model for complex disorders.

Selected heterozygosity at cis-regulatory sequences increases the expression homogeneity of a cell population in humans

Background

Examples of heterozygote advantage in humans are scarce and limited to protein-coding sequences. Here, we attempt a genome-wide functional inference of advantageous heterozygosity at cis-regulatory regions.

Results

The single-nucleotide polymorphisms bearing the signatures of balancing selection are enriched in active cis-regulatory regions of immune cells and epithelial cells, the latter of which provide barrier function and innate immunity. Examples associated with ancient trans-specific balancing selection are also discovered. Allelic imbalance in chromatin accessibility and divergence in transcription factor motif sequences indicate that these balanced polymorphisms cause distinct regulatory variation. However, a majority of these variants show no association with the expression level of the target gene. Instead, single-cell experimental data for gene expression and chromatin accessibility demonstrate that heterozygous sequences can lower cell-to-cell variability in proportion to selection strengths. This negative correlation is more pronounced for highly expressed genes and consistently observed when using different data and methods. Based on mathematical modeling, we hypothesize that extrinsic noise from fluctuations in transcription factor activity may be amplified in homozygotes, whereas it is buffered in heterozygotes. While high expression levels are coupled with intrinsic noise reduction, regulatory heterozygosity can contribute to the suppression of extrinsic noise.

Conclusions

This mechanism may confer a selective advantage by increasing cell population homogeneity and thereby enhancing the collective action of the cells, especially of those involved in the defense systems in humans.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1027-8) contains supplementary material, which is available to authorized users.

Transcriptome Analysis on Monocytes from Patients with Neovascular Age-Related Macular Degeneration

Mononuclear phagocytes (MPs), including monocytes/macrophages, play complex roles in age-related macular degeneration (AMD) pathogenesis. We reported altered gene-expression signature in peripheral blood mononuclear cells from AMD patients, and a chemokine receptor signature on AMD monocytes. To obtain comprehensive understanding of MP involvement, particularly in peripheral circulation in AMD, we performed global gene expression analysis in monocytes. We separated monocytes from treatment-naïve neovascular AMD (nvAMD) patients (n = 14) and age-matched controls (n = 15), and performed microarray and bioinformatics analysis. Quantitative real-time PCR was performed on other sets of nvAMD (n = 25), atrophic AMD (n = 21), and controls (n = 28) for validation. This validated microarray genes (like TMEM176A/B and FOSB) tested, including differences between nvAMD and atrophic AMD. We identified 2,165 differentially-expressed genes (P < 0.05), including 79 genes with log2 fold change ≥1.5 between nvAMD and controls. Functional annotation using DAVID and TANGO demonstrated immune response alterations in AMD monocytes (FDR-P <0.05), validated by randomized data comparison (P < 0.0001). GSEA, ISMARA, and MEME analysis found immune enrichment and specific involved microRNAs. Enrichment of differentially-expressed genes in monocytes was found in retina via SAGE data-mining. These genes were enriched in non-classical vs. classical monocyte subsets (P < 0.05). Therefore, global gene expression analysis in AMD monocytes reveals an altered immune-related signature, further implicating systemic MP activation in AMD.

PU.1 is involved in the immune response to Aspergillus fumigatus through upregulating Dectin-1 expression

Background

Invasive aspergillosis is a life-threatening disease, and its incidence has increased in the recent past. Dectin-1 recognizes β-glucans and mediates innate immune responses to Aspergillus fumigatus. Transcription factor PU.1 has been the focus of recent research due to its role in inflammation and infection. However, its role in Dectin-1 regulation during A. fumigatus infection remains to be elucidated.

Methods

THP-1 cells were stimulated with A. fumigatus conidia. We then used real-time RT-PCR, Western blot, and immunofluorescence assays to analyze the mRNA and protein levels and cellular distribution, respectively, of Dectin-1 and PU.1 in stimulated THP-1 cells. Additionally, we used the luciferase reporter assays, chromatin immunoprecipitation (ChIP) assays, electrophoretic mobility shift assays (EMSA), and RNA interference experiments to investigate the role of PU.1 in Dectin-1 regulation.

Results

Our results revealed that Dectin-1 mRNA and protein levels as well as the PU.1 protein level were increased in THP-1 cells stimulated with A. fumigatus conidia, while the mRNA expression level did not significantly change between the stimulated and control groups. We also observed that PU.1 translocated into the nucleus in stimulated THP-1 cells. The results of the luciferase reporter assay showed that PU.1 promoted human Dectin-1 (hDectin-1) gene activity. ChIP and EMSA indicated that PU.1 could bind with hDectin-1 gene promoter at three potential transcription factor-binding sites (TFBSs). In addition, knockdown of PU.1 significantly decreased Dectin-1 expression.

Conclusions

This study demonstrated the novel role of PU.1 in the immune response to A. fumigatus through upregulation of Dectin-1 expression and its translocation to the nucleus in A. fumigatus-stimulated THP-1 cells.

Next generation sequencing technology and genomewide data analysis: Perspectives for retinal research

The advent of high throughput next generation sequencing (NGS) has accelerated the pace of discovery of disease-associated genetic variants and genomewide profiling of expressed sequences and epigenetic marks, thereby permitting systems-based analyses of ocular development and disease. Rapid evolution of NGS and associated methodologies presents significant challenges in acquisition, management, and analysis of large data sets and for extracting biologically or clinically relevant information. Here we illustrate the basic design of commonly used NGS-based methods, specifically whole exome sequencing, transcriptome, and epigenome profiling, and provide recommendations for data analyses. We briefly discuss systems biology approaches for integrating multiple data sets to elucidate gene regulatory or disease networks. While we provide examples from the retina, the NGS guidelines reviewed here are applicable to other tissues/cell types as well.

Gene regulatory mechanisms underpinning prostate cancer susceptibility

Molecular characterization of genome-wide association study (GWAS) loci can uncover key genes and biological mechanisms underpinning complex traits and diseases. Here we present deep, high-throughput characterization of gene regulatory mechanisms underlying prostate cancer risk loci. Our methodology integrates data from 295 prostate cancer chromatin immunoprecipitation and sequencing experiments with genotype and gene expression data from 602 prostate tumor samples. The analysis identifies new gene regulatory mechanisms affected by risk locus SNPs, including widespread disruption of ternary androgen receptor (AR)-FOXA1 and AR-HOXB13 complexes and competitive binding mechanisms. We identify 57 expression quantitative trait loci at 35 risk loci, which we validate through analysis of allele-specific expression. We further validate predicted regulatory SNPs and target genes in prostate cancer cell line models. Finally, our integrated analysis can be accessed through an interactive visualization tool. This analysis elucidates how genome sequence variation affects disease predisposition via gene regulatory mechanisms and identifies relevant genes for downstream biomarker and drug development.

Natural selection in a population of Drosophila melanogaster explained by changes in gene expression caused by sequence variation in core promoter regions

Background

Understanding the evolutionary forces that influence variation in gene regulatory regions in natural populations is an important challenge for evolutionary biology because natural selection for such variations could promote adaptive phenotypic evolution. Recently, whole-genome sequence analyses have identified regulatory regions subject to natural selection. However, these studies could not identify the relationship between sequence variation in the detected regions and change in gene expression levels. We analyzed sequence variations in core promoter regions, which are critical regions for gene regulation in higher eukaryotes, in a natural population of Drosophila melanogaster, and identified core promoter sequence variations associated with differences in gene expression levels subjected to natural selection.

Results

Among the core promoter regions whose sequence variation could change transcription factor binding sites and explain differences in expression levels, three core promoter regions were detected as candidates associated with purifying selection or selective sweep and seven as candidates associated with balancing selection, excluding the possibility of linkage between these regions and core promoter regions. CHKov1, which confers resistance to the sigma virus and related insecticides, was identified as core promoter regions that has been subject to selective sweep, although it could not be denied that selection for variation in core promoter regions was due to linked single nucleotide polymorphisms in the regulatory region outside core promoter regions. Nucleotide changes in core promoter regions of CHKov1 caused the loss of two basal transcription factor binding sites and acquisition of one transcription factor binding site, resulting in decreased gene expression levels. Of nine core promoter regions regions associated with balancing selection, brat, and CG9044 are associated with neuromuscular junction development, and Nmda1 are associated with learning, behavioral plasticity, and memory. Diversity of neural and behavioral traits may have been maintained by balancing selection.

Conclusions

Our results revealed the evolutionary process occurring by natural selection for differences in gene expression levels caused by sequence variation in core promoter regions in a natural population. The sequences of core promoter regions were diverse even within the population, possibly providing a source for natural selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0606-3) contains supplementary material, which is available to authorized users.

Relevance of Sp Binding Site Polymorphism in WWOX for Treatment Outcome in Pancreatic Cancer

Background:

A genome-wide association study (GWAS) suggested inherited genetic single-nucleotide polymorphisms (SNPs) affecting overall survival (OS) in advanced pancreatic cancer. To identify robust clinical biomarkers, we tested the strongest reported candidate loci in an independent patient cohort, assessed cellular drug sensitivity, and evaluated molecular effects.

Methods:

This study comprised 381 patients with histologically verified pancreatic ductal adenocarcinoma treated with gemcitabine-based chemotherapy. The primary outcome was the relationship between germline polymorphisms and OS. Functional assays addressed pharmacological dose-response effects in lymphoblastoid cell lines (LCLs) and pancreatic cancer cell lines (including upon RNAi), gene expression analyses, and allele-specific transcription factor binding. All statistical tests were two-sided.

Results:

The A allele (26% in Caucasians) at SNP rs11644322 in the putative tumor suppressor gene WWOX conferred worse prognosis. Median OS was 14 months (95% confidence interval [CI] = 12 to 15 months), 13 months (95% CI = 11 to 15 months), and nine months (95% CI = 7 to 12 months) for the GG, GA, and AA genotypes, respectively (P_trend < .001 for trend in univariate log-rank assuming a codominant mode of inheritance; advanced disease subgroup P_trend < .001). Mean OS was 25 months (95% CI = 21 to 29 months), 19 months (95% CI = 15 to 22 months), and 13 months (95% CI = 10 to 16 months), respectively. This effect held true after adjustment for age, performance status according to Eastern Cooperative Oncology Group classification, TNM, grading, and resection status and was comparable with the strongest established prognostic factors in multivariable analysis. Consistently, reduced responsiveness to gemcitabine, but not 5-fluorouracil, along with lower WWOX expression was demonstrated in LCLs harboring the AA genotype. Likewise, RNAi-mediated WWOX knockdown in pancreatic cancer cells confirmed differential cytostatic drug sensitivity. In electrophoretic mobility shift assays, the A allele exhibited weaker binding of Sp family members Sp1/Sp3.

Conclusions:

WWOX rs11644322 represents a major predictive factor in gemcitabine-treated pancreatic cancer. Decreased WWOX expression may interfere with gemcitabine sensitivity, and allele-specific binding at rs11644332 might be a causative molecular mechanism behind the observed clinical associations.

Erosion of Conserved Binding Sites in Personal Genomes Points to Medical Histories

Although many human diseases have a genetic component involving many loci, the majority of studies are statistically underpowered to isolate the many contributing variants, raising the question of the existence of alternate processes to identify disease mutations. To address this question, we collect ancestral transcription factor binding sites disrupted by an individual's variants and then look for their most significant congregation next to a group of functionally related genes. Strikingly, when the method is applied to five different full human genomes, the top enriched function for each is invariably reflective of their very different medical histories. For example, our method implicates "abnormal cardiac output" for a patient with a longstanding family history of heart disease, "decreased circulating sodium level" for an individual with hypertension, and other biologically appealing links for medical histories spanning narcolepsy to axonal neuropathy. Our results suggest that erosion of gene regulation by mutation load significantly contributes to observed heritable phenotypes that manifest in the medical history. The test we developed exposes a hitherto hidden layer of personal variants that promise to shed new light on human disease penetrance, expressivity and the sensitivity with which we can detect them.

Small Genetic Circuits and MicroRNAs: Big Players in Polymerase II Transcriptional Control in Plants

RNA Polymerase II (Pol II) regulatory cascades involving transcription factors (TFs) and their targets orchestrate the genetic circuitry of every eukaryotic organism. In order to understand how these cascades function, they can be dissected into small genetic networks, each containing just a few Pol II transcribed genes, that generate specific signal-processing outcomes. Small RNA regulatory circuits involve direct regulation of a small RNA by a TF and/or direct regulation of a TF by a small RNA and have been shown to play unique roles in many organisms. Here, we will focus on small RNA regulatory circuits containing Pol II transcribed microRNAs (miRNAs). While the role of miRNA-containing regulatory circuits as modular building blocks for the function of complex networks has long been on the forefront of studies in the animal kingdom, plant studies are poised to take a lead role in this area because of their advantages in probing transcriptional and posttranscriptional control of Pol II genes. The relative simplicity of tissue- and cell-type organization, miRNA targeting, and genomic structure make the Arabidopsis thaliana plant model uniquely amenable for small RNA regulatory circuit studies in a multicellular organism. In this Review, we cover analysis, tools, and validation methods for probing the component interactions in miRNA-containing regulatory circuits. We then review the important roles that plant miRNAs are playing in these circuits and summarize methods for the identification of small genetic circuits that strongly influence plant function. We conclude by noting areas of opportunity where new plant studies are imminently needed.