The race for the discovery of enhancers at a genome-wide scale has been on since the commencement of next generation sequencing decades after the discovery of the first enhancer, SV40. A few enhancer-predicting features such as chromatin feature, histone modifications and sequence feature had been implemented with varying success rates. However, to date, there is no consensus yet on the single enhancer marker that can be employed to ultimately distinguish and uncover enhancers from the enormous genomic regions. Many supervised, unsupervised and semi-supervised computational approaches had emerged to complement and facilitate experimental approaches in enhancer discovery. In this review, we placed our focus on the recently emerged enhancer predictor tools that work on general enhancer features such as sequences, chromatin states and histone modifications, eRNA and of multiple feature approach. Comparisons of their prediction methods and outcomes were done across their functionally similar counterparts. We provide some recommendations and insights for future development of more comprehensive and robust tools.
Using Genetic Algorithm, this paper presents a modelling method to generate novel logical-based features from DNA sequences enriched with H3K4mel histone signatures. Current histone signature is mostly represented using k-mers content features incapable of representing all the possible complex interactions of various DNA segments. The main contributions are, among others: (a) demonstrating that there are complex interactions among sequence segments in the histone regions; (b) developing a parse tree representation of the logical complex features. The proposed novel feature is compared to the k-mers content features using datasets from the mouse (mm9) genome. Evaluation results show that the new feature improves the prediction performance as shown by f-measure for all datasets tested. Also, it is discovered that tree-based features generated from a single chromosome can be generalized to predict histone marks in other chromosomes not used in the training. These findings have a great impact on feature design considerations for histone signatures as well as other classifier design features.
Progress in the functional studies of human olfactory receptors has been largely hampered by the lack of a reliable experimental model system. Although transgenic approaches in mice could characterize the function of individual olfactory receptors, the presence of over 300 functional genes in the human genome becomes a daunting task. Thus, the characterization of individuals with a genetic susceptibility to altered olfaction coupled with the absence of particular olfactory receptor genes will allow phenotype/genotype correlations and vindicate the function of specific olfactory receptors with their cognate ligands. We characterized a 118 kb β-globin deletion and found that its 3' end breakpoint extends to the neighboring olfactory receptor region downstream of the β-globin gene cluster. This deletion encompasses six contiguous olfactory receptor genes (OR51V1, OR52Z1, OR51A1P, OR52A1, OR52A5, and OR52A4) all of which are expressed in the brain. Topology analysis of the encoded proteins from these olfactory receptor genes revealed that OR52Z1, OR52A1, OR52A5, and OR52A4 are predicted to be functional receptors as they display integral characteristics of G-proteins coupled receptors. Individuals homozygous for the 118 kb β-globin deletion are afflicted with β-thalassemia due to a homozygous deletion of the β-globin gene and have no alleles for the above mentioned olfactory receptors genes. This is the first example of a homozygous deletion of olfactory receptor genes in human. Although altered olfaction remains to be ascertained in these individuals, such a study can be carried out in β-thalassemia patients from Malaysia, Indonesia and the Philippines where this mutation is common. Furthermore, OR52A1 contains a γ-globin enhancer, which was previously shown to confer continuous expression of the fetal γ-globin genes. Thus, the hypothesis that β-thalassemia individuals, who are homozygous for the 118 kb deletion, may also have an exacerbation of their anemia due to the deletion of two copies of the γ-globin enhancer element is worthy of consideration.
Transcriptional networks are critical for the establishment of tissue-specific cellular states in health and disease, including cancer. Yet, the transcriptional circuits that control carcinogenesis remain poorly understood. Here we report that Kruppel like factor 6 (KLF6), a transcription factor of the zinc finger family, regulates lipid homeostasis in clear cell renal cell carcinoma (ccRCC). We show that KLF6 supports the expression of lipid metabolism genes and promotes the expression of PDGFB, which activates mTOR signalling and the downstream lipid metabolism regulators SREBF1 and SREBF2. KLF6 expression is driven by a robust super enhancer that integrates signals from multiple pathways, including the ccRCC-initiating VHL-HIF2A pathway. These results suggest an underlying mechanism for high mTOR activity in ccRCC cells. More generally, the link between super enhancer-driven transcriptional networks and essential metabolic pathways may provide clues to the mechanisms that maintain the stability of cell identity-defining transcriptional programmes in cancer.
Genome-wide association studies have identified SNPs near ZNF365 at 10q21.2 that are associated with both breast cancer risk and mammographic density. To identify the most likely causal SNPs, we fine mapped the association signal by genotyping 428 SNPs across the region in 89,050 European and 12,893 Asian case and control subjects from the Breast Cancer Association Consortium. We identified four independent sets of correlated, highly trait-associated variants (iCHAVs), three of which were located within ZNF365. The most strongly risk-associated SNP, rs10995201 in iCHAV1, showed clear evidence of association with both estrogen receptor (ER)-positive (OR = 0.85 [0.82-0.88]) and ER-negative (OR = 0.87 [0.82-0.91]) disease, and was also the SNP most strongly associated with percent mammographic density. iCHAV2 (lead SNP, chr10: 64,258,684:D) and iCHAV3 (lead SNP, rs7922449) were also associated with ER-positive (OR = 0.93 [0.91-0.95] and OR = 1.06 [1.03-1.09]) and ER-negative (OR = 0.95 [0.91-0.98] and OR = 1.08 [1.04-1.13]) disease. There was weaker evidence for iCHAV4, located 5' of ADO, associated only with ER-positive breast cancer (OR = 0.93 [0.90-0.96]). We found 12, 17, 18, and 2 candidate causal SNPs for breast cancer in iCHAVs 1-4, respectively. Chromosome conformation capture analysis showed that iCHAV2 interacts with the ZNF365 and NRBF2 (more than 600 kb away) promoters in normal and cancerous breast epithelial cells. Luciferase assays did not identify SNPs that affect transactivation of ZNF365, but identified a protective haplotype in iCHAV2, associated with silencing of the NRBF2 promoter, implicating this gene in the etiology of breast cancer.
Genome-wide association studies (GWASs) have revealed increased breast cancer risk associated with multiple genetic variants at 5p12. Here, we report the fine mapping of this locus using data from 104,660 subjects from 50 case-control studies in the Breast Cancer Association Consortium (BCAC). With data for 3,365 genotyped and imputed SNPs across a 1 Mb region (positions 44,394,495-45,364,167; NCBI build 37), we found evidence for at least three independent signals: the strongest signal, consisting of a single SNP rs10941679, was associated with risk of estrogen-receptor-positive (ER+) breast cancer (per-g allele OR ER+ = 1.15; 95% CI 1.13-1.18; p = 8.35 × 10-30). After adjustment for rs10941679, we detected signal 2, consisting of 38 SNPs more strongly associated with ER-negative (ER-) breast cancer (lead SNP rs6864776: per-a allele OR ER- = 1.10; 95% CI 1.05-1.14; p conditional = 1.44 × 10-12), and a single signal 3 SNP (rs200229088: per-t allele OR ER+ = 1.12; 95% CI 1.09-1.15; p conditional = 1.12 × 10-05). Expression quantitative trait locus analysis in normal breast tissues and breast tumors showed that the g (risk) allele of rs10941679 was associated with increased expression of FGF10 and MRPS30. Functional assays demonstrated that SNP rs10941679 maps to an enhancer element that physically interacts with the FGF10 and MRPS30 promoter regions in breast cancer cell lines. FGF10 is an oncogene that binds to FGFR2 and is overexpressed in ∼10% of human breast cancers, whereas MRPS30 plays a key role in apoptosis. These data suggest that the strongest signal of association at 5p12 is mediated through coordinated activation of FGF10 and MRPS30, two candidate genes for breast cancer pathogenesis.