Human Genes Encoding Transcription Factors and Chromatin-Modifying Proteins Have Low Levels of Promoter Polymorphism: A Study of 1000 Genomes Project Data

Transcription factors ERα and Sox2 have differing multiphasic DNA- and RNA-binding mechanisms

Many transcription factors (TFs) have been shown to bind RNA, leading to open questions regarding the mechanism(s) of this RNA binding and its role in regulating TF activities. Here, we use biophysical assays to interrogate the k on, k off, and K d for DNA and RNA binding of two model human TFs, ERα and Sox2. Unexpectedly, we found that both proteins exhibit multiphasic nucleic acid-binding kinetics. We propose that Sox2 RNA and DNA multiphasic binding kinetics can be explained by a conventional model for sequential Sox2 monomer association and dissociation. In contrast, ERα nucleic acid binding exhibited biphasic dissociation paired with novel triphasic association behavior, in which two apparent binding transitions are separated by a 10-20 min "lag" phase depending on protein concentration. We considered several conventional models for the observed kinetic behavior, none of which adequately explained all the ERα nucleic acid-binding data. Instead, simulations with a model incorporating sequential ERα monomer association, ERα nucleic acid complex isomerization, and product "feedback" on isomerization rate recapitulated the general kinetic trends for both ERα DNA and RNA binding. Collectively, our findings reveal that Sox2 and ERα bind RNA and DNA with previously unappreciated multiphasic binding kinetics, and that their reaction mechanisms differ with ERα binding nucleic acids via a novel reaction mechanism.

TF-Prioritizer: a Java pipeline to prioritize condition-specific transcription factors

BACKGROUND

Eukaryotic gene expression is controlled by cis-regulatory elements (CREs), including promoters and enhancers, which are bound by transcription factors (TFs). Differential expression of TFs and their binding affinity at putative CREs determine tissue- and developmental-specific transcriptional activity. Consolidating genomic datasets can offer further insights into the accessibility of CREs, TF activity, and, thus, gene regulation. However, the integration and analysis of multimodal datasets are hampered by considerable technical challenges. While methods for highlighting differential TF activity from combined chromatin state data (e.g., chromatin immunoprecipitation [ChIP], ATAC, or DNase sequencing) and RNA sequencing data exist, they do not offer convenient usability, have limited support for large-scale data processing, and provide only minimal functionality for visually interpreting results.

RESULTS

We developed TF-Prioritizer, an automated pipeline that prioritizes condition-specific TFs from multimodal data and generates an interactive web report. We demonstrated its potential by identifying known TFs along with their target genes, as well as previously unreported TFs active in lactating mouse mammary glands. Additionally, we studied a variety of ENCODE datasets for cell lines K562 and MCF-7, including 12 histone modification ChIP sequencing as well as ATAC and DNase sequencing datasets, where we observe and discuss assay-specific differences.

CONCLUSION

TF-Prioritizer accepts ATAC, DNase, or ChIP sequencing and RNA sequencing data as input and identifies TFs with differential activity, thus offering an understanding of genome-wide gene regulation, potential pathogenesis, and therapeutic targets in biomedical research.

Haem oxygenase-1 and haemopexin gene polymorphisms and the risk of anti-tuberculosis drug-induced hepatotoxicity in China

Objective: To assess whether the risk of anti-tuberculosis drug-induced hepatotoxicity (ATDH) might be influenced by haem oxygenase-1 (HMOX1) and haemopexin (HPX) gene polymorphisms. Methods: A dynamic anti-tuberculosis treatment cohort was constructed, and the 1:4 matched nested case-control study was analysed. Eight single nucleotide polymorphisms (SNPs) of the two genes were selected for genotyping and Bonferroni correction was performed to correct for multiple comparison. Results: Overall, 7.8% of patients developed ATDH. SNP rs1807714 in the HMOX1 gene had decreased effects on the risk of moderate and severe hepatotoxicity under the dominant and additive models, and hepatocellular injury under the additive model. SNP rs2682099 in the HPX gene had increased effects on the risk of moderate and severe hepatotoxicity under the recessive model. However, these associations disappeared after Bonferroni correction. Conclusion: HMOX1 and HPX gene polymorphisms might not be associated with susceptibility to ATDH in the Chinese population.

Dissecting role of founder mutation p.V727M in GNE in Indian HIBM cohort

GNE gene-specific c.2179G>A(p.V727M) is a key alteration reported in patients with hereditary inclusion body myopathy (HIBM) and represents an ethnic founder mutation in the Indian cohort. However, the underlying role of this mutation in pathogenesis remains largely unknown. Thus, in this study, we aimed to access possible mechanisms of V727M mutation that could be leading to myopathy. We evaluated various in silico tools to predict the effect of this mutation on pathogenicity, structural or possible interactions, that could induce myopathy. Our results propose that V727M mutation could induce deleterious effects or pathogenicity and affect the stability of GNE protein. Analysis of differential genes reported in the V727 mutant case suggests that it can affect GNE protein interaction with Myc-proto-oncogene (MYC) transcription factor. Our in silico analysis also suggests a possible interaction between GNE ManNac-kinase domain with MYC protein at the C-terminal DNA-binding domain. MYC targets reported in skeletal muscles via ChIP-seq suggest that it plays a key role in regulating the expression of many genes reported differentially expressed in V727M-mutated HIBMs. We conclude that V727M mutation could alter the interaction of GNE with MYC thereby altering transcription of sialyltransferase and neuromuscular genes, thus understanding these effects could pave the way for developing effective therapies against HIBM.

Evidence of widespread, independent sequence signature for transcription factor cobinding

Transcription factors (TFs) are the vocabulary that genomes use to regulate gene expression and phenotypes. The interactions among TFs enrich this vocabulary and orchestrate diverse biological processes. Although simple models identify open chromatin and the presence of TF motifs as the two major contributors to TF binding patterns, it remains elusive what contributes to the in vivo TF cobinding landscape. In this study, we developed a machine learning algorithm to explore the contributors of the cobinding patterns. The algorithm substantially outperforms the state-of-the-field models for TF cobinding prediction. Game theory-based feature importance analysis reveals that, for most of the TF pairs we studied, independent motif sequences contribute one or more of the two TFs under investigation to their cobinding patterns. Such independent motif sequences include, but are not limited to, transcription initiation-related proteins and known TF complexes. We found the motif sequence signatures and the TFs are rarely mutual, corroborating a hierarchical and directional organization of the regulatory network and refuting the possibility of artifacts caused by shared sequence similarity with the TFs under investigation. We modeled such regulatory language with directed graphs, which reveal shared, global factors that are related to many binding and cobinding patterns.

Cellular reprogramming: Mathematics meets medicine

Generating needed cell types using cellular reprogramming is a promising strategy for restoring tissue function in injury or disease. A common method for reprogramming is addition of one or more transcription factors that confer a new function or identity. Advancements in transcription factor selection and delivery have culminated in successful grafting of autologous reprogrammed cells, an early demonstration of their clinical utility. Though cellular reprogramming has been successful in a number of settings, identification of appropriate transcription factors for a particular transformation has been challenging. Computational methods enable more sophisticated prediction of relevant transcription factors for reprogramming by leveraging gene expression data of initial and target cell types, and are built on mathematical frameworks ranging from information theory to control theory. This review highlights the utility and impact of these mathematical frameworks in the field of cellular reprogramming. This article is categorized under: Reproductive System Diseases > Reproductive System Diseases>Genetics/Genomics/Epigenetics Reproductive System Diseases > Reproductive System Diseases>Stem Cells and Development Reproductive System Diseases > Reproductive System Diseases>Computational Models.

Deletion of the transcription factor Prox-1 specifically in the renal distal convoluted tubule causes hypomagnesemia via reduced expression of TRPM6 and NCC

The renal distal convoluted tubule (DCT) is critical for the fine-tuning of urinary ion excretion and the control of blood pressure. Ion transport along the DCT is tightly controlled by posttranscriptional mechanisms including a complex interplay of kinases, phosphatases, and ubiquitin ligases. Previous work identified the transcription factor Prox-1 as a gene significantly enriched in the DCT of adult mice. To test if Prox-1 contributes to the transcriptional regulation of DCT function and structure, we developed a novel mouse model (NCC^cre:Prox-1^flox/flox) for an inducible deletion of Prox-1 specifically in the DCT. The deletion of Prox-1 had no obvious impact on DCT structure and growth independent whether the deletion was achieved in newborn or adult mice. Furthermore, DCT-specific Prox-1 deficiency did not alter DCT-proliferation in response to loop diuretic treatment. Likewise, the DCT-specific deletion of Prox-1 did not cause other gross phenotypic abnormalities. Body weight, urinary volume, Na⁺ and K⁺ excretion as well as plasma Na⁺, K⁺, and aldosterone levels were similar in Prox-1_DCT^KO and Prox-1_DCT^Ctrl mice. However, Prox-1_DCT^KO mice exhibited a significant hypomagnesemia with a profound downregulation of the DCT-specific apical Mg²⁺ channel TRPM6 and the NaCl cotransporter (NCC) at both mRNA and protein levels. The expression of other proteins involved in distal tubule Mg²⁺ and Na⁺ handling was not affected. Thus, Prox-1 is a DCT-enriched transcription factor that does not control DCT growth but contributes to the molecular control of DCT-dependent Mg²⁺ homeostasis in the adult kidney.

Exome-wide search and functional annotation of genes associated in patients with severe tick-borne encephalitis in a Russian population

Background

Tick-borne encephalitis (TBE) is a viral infectious disease caused by tick-borne encephalitis virus (TBEV). TBEV infection is responsible for a variety of clinical manifestations ranging from mild fever to severe neurological illness. Genetic factors involved in the host response to TBEV that may potentially play a role in the severity of the disease are still poorly understood. In this study, using whole-exome sequencing, we aimed to identify genetic variants and genes associated with severe forms of TBE as well as biological pathways through which the identified variants may influence the severity of the disease.

Results

Whole-exome sequencing data analysis was performed on 22 Russian patients with severe forms of TBE and 17 Russian individuals from the control group. We identified 2407 candidate genes harboring rare, potentially pathogenic variants in exomes of patients with TBE and not containing any rare, potentially pathogenic variants in exomes of individuals from the control group. According to DAVID tool, this set of 2407 genes was enriched with genes involved in extracellular matrix proteoglycans pathway and genes encoding proteins located at the cell periphery. A total of 154 genes/proteins from these functional groups have been shown to be involved in protein-protein interactions (PPIs) with the known candidate genes/proteins extracted from TBEVHostDB database. By ranking these genes according to the number of rare harmful minor alleles, we identified two genes (MSR1 and LMO7), harboring five minor alleles, and three genes (FLNA, PALLD, PKD1) harboring four minor alleles.

When considering genes harboring genetic variants associated with severe forms of TBE at the suggestive P-value < 0.01, 46 genes containing harmful variants were identified. Out of these 46 genes, eight (MAP4, WDFY4, ACTRT2, KLHL25, MAP2K3, MBD1, OR10J1, and OR2T34) were additionally found among genes containing rare pathogenic variants identified in patients with TBE; and five genes (WDFY4,ALK, MAP4, BNIPL, EPPK1) were found to encode proteins that are involved in PPIs with proteins encoded by genes from TBEVHostDB. Three genes out of five (MAP4, EPPK1, ALK) were found to encode proteins located at cell periphery.

Conclusions

Whole-exome sequencing followed by systems biology approach enabled to identify eight candidate genes (MAP4, WDFY4, ACTRT2, KLHL25, MAP2K3, MBD1, OR10J1, and OR2T34) that can potentially determine predisposition to severe forms of TBE. Analyses of the genetic risk factors for severe forms of TBE revealed a significant enrichment with genes controlling extracellular matrix proteoglycans pathway as well as genes encoding components of cell periphery.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0503-x) contains supplementary material, which is available to authorized users.

Applications of the 1000 Genomes Project resources

The 1000 Genomes Project created a valuable, worldwide reference for human genetic variation. Common uses of the 1000 Genomes dataset include genotype imputation supporting Genome-wide Association Studies, mapping expression Quantitative Trait Loci, filtering non-pathogenic variants from exome, whole genome and cancer genome sequencing projects, and genetic analysis of population structure and molecular evolution. In this article, we will highlight some of the multiple ways that the 1000 Genomes data can be and has been utilized for genetic studies.

A database of human genes and a gene network involved in response to tick-borne encephalitis virus infection

BACKGROUND

Tick-borne encephalitis is caused by the neurotropic, positive-sense RNA virus, tick-borne encephalitis virus (TBEV). TBEV infection can lead to a variety of clinical manifestations ranging from slight fever to severe neurological illness. Very little is known about genetic factors predisposing to severe forms of disease caused by TBEV. The aims of the study were to compile a catalog of human genes involved in response to TBEV infection and to rank genes from the catalog based on the number of neighbors in the network of pairwise interactions involving these genes and TBEV RNA or proteins.

RESULTS

Based on manual review and curation of scientific publications a catalog comprising 140 human genes involved in response to TBEV infection was developed. To provide access to data on all genes, the TBEVhostDB web resource ( http://icg.nsc.ru/TBEVHostDB/ ) was created. We reconstructed a network formed by pairwise interactions between TBEV virion itself, viral RNA and viral proteins and 140 genes/proteins from TBEVHostDB. Genes were ranked according to the number of interactions in the network. Two genes/proteins (CCR5 and IFNAR1) that had maximal number of interactions were revealed. It was found that the subnetworks formed by CCR5 and IFNAR1 and their neighbors were a fragments of two key pathways functioning during the course of tick-borne encephalitis: (1) the attenuation of interferon-I signaling pathway by the TBEV NS5 protein that targeted peptidase D; (2) proinflammation and tissue damage pathway triggered by chemokine receptor CCR5 interacting with CD4, CCL3, CCL4, CCL2. Among nine genes associated with severe forms of TBEV infection, three genes/proteins (CCR5, IL10, ARID1B) were found to have protein-protein interactions within the network, and two genes/proteins (IFNL3 and the IL10, that was just mentioned) were up- or down-regulated in response to TBEV infection. Based on this finding, potential mechanisms for participation of CCR5, IL10, ARID1B, and IFNL3 in the host response to TBEV infection were suggested.

CONCLUSIONS

A database comprising 140 human genes involved in response to TBEV infection was compiled and the TBEVHostDB web resource, providing access to all genes was created. This is the first effort of integrating and unifying data on genetic factors that may predispose to severe forms of diseases caused by TBEV. The TBEVHostDB could potentially be used for assessment of risk factors for severe forms of tick-borne encephalitis and for the design of personalized pharmacological strategies for the treatment of TBEV infection.

A compendium and functional characterization of mammalian genes involved in adaptation to Arctic or Antarctic environments

Background

Many mammals are well adapted to surviving in extremely cold environments. These species have likely accumulated genetic changes that help them efficiently cope with low temperatures. It is not known whether the same genes related to cold adaptation in one species would be under selection in another species. The aims of this study therefore were: to create a compendium of mammalian genes related to adaptations to a low temperature environment; to identify genes related to cold tolerance that have been subjected to independent positive selection in several species; to determine promising candidate genes/pathways/organs for further empirical research on cold adaptation in mammals.

Results

After a search for publications containing keywords: “whole genome”, “transcriptome or exome sequencing data”, and “genome-wide genotyping array data” authors looked for information related to genetic signatures ascribable to positive selection in Arctic or Antarctic mammalian species. Publications related to Human, Arctic fox, Yakut horse, Mammoth, Polar bear, and Minke whale were chosen. The compendium of genes that potentially underwent positive selection in >1 of these six species consisted of 416 genes. Twelve of them showed traces of positive selection in three species. Gene ontology term enrichment analysis of 416 genes from the compendium has revealed 13 terms relevant to the scope of this study. We found that enriched terms were relevant to three major groups: terms associated with collagen proteins and the extracellular matrix; terms associated with the anatomy and physiology of cilium; terms associated with docking. We further revealed that genes from compendium were over-represented in the lists of genes expressed in the lung and liver.

Conclusions

A compendium combining mammalian genes involved in adaptation to cold environment was designed, based on the intersection of positively selected genes from six Arctic and Antarctic species. The compendium contained 416 genes that have been positively selected in at least two species. However, we did not reveal any positively selected genes that would be related to cold adaptation in all species from our list. But, our work points to several strong candidate genes involved in mechanisms and biochemical pathways related to cold adaptation response in different species.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0580-9) contains supplementary material, which is available to authorized users.

A compendium of human genes regulating feeding behavior and body weight, its functional characterization and identification of GWAS genes involved in brain-specific PPI network

Background

Obesity is heritable. It predisposes to many diseases. The objectives of this study were to create a compendium of genes relevant to feeding behavior (FB) and/or body weight (BW) regulation; to construct and to analyze networks formed by associations between genes/proteins; and to identify the most significant genes, biological processes/pathways, and tissues/organs involved in BW regulation.

Results

The compendium of genes controlling FB or BW includes 578 human genes. Candidate genes were identified from various sources, including previously published original research and review articles, GWAS meta-analyses, and OMIM (Online Mendelian Inheritance in Man). All genes were ranked according to knowledge about their biological role in body weight regulation and classified according to expression patterns or functional characteristics. Substantial and overrepresented numbers of genes from the compendium encoded cell surface receptors, signaling molecules (hormones, neuropeptides, cytokines), transcription factors, signal transduction proteins, cilium and BBSome components, and lipid binding proteins or were present in the brain-specific list of tissue-enriched genes identified with TSEA tool. We identified 27 pathways from KEGG, REACTOME and BIOCARTA whose genes were overrepresented in the compendium. Networks formed by physical interactions or homological relationships between proteins or interactions between proteins involved in biochemical/signaling pathways were reconstructed and analyzed. Subnetworks and clusters identified by the MCODE tool included genes/proteins associated with cilium morphogenesis, signal transduction proteins (particularly, G protein–coupled receptors, kinases or proteins involved in response to insulin stimulus) and transcription regulation (particularly nuclear receptors). We ranked GWAS genes according to the number of neighbors in three networks and revealed 22 GWAS genes involved in the brain-specific PPI network. On the base of the most reliable PPIs functioning in the brain tissue, new regulatory schemes interpreting relevance to BW regulation are proposed for three GWAS genes (ETV5, LRP1B, and NDUFS3).

Conclusions

A compendium comprising 578 human genes controlling FB or BW was designed, and the most significant functional groups of genes, biological processes/pathways, and tissues/organs involved in BW regulation were revealed. We ranked genes from the GWAS meta-analysis set according to the number and quality of associations in the networks and then according to their involvement in the brain-specific PPI network and proposed new regulatory schemes involving three GWAS genes (ETV5, LRP1B, and NDUFS3) in BW regulation. The compendium is expected to be useful for pathology risk estimation and for design of new pharmacological approaches in the treatment of human obesity.

Electronic supplementary material

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