The eukaryotic promoter database (EPD)

DeepCORE: An interpretable multi-view deep neural network model to detect co-operative regulatory elements

Gene transcription is an essential process involved in all aspects of cellular functions with significant impact on biological traits and diseases. This process is tightly regulated by multiple elements that co-operate to jointly modulate the transcription levels of target genes. To decipher the complicated regulatory network, we present a novel multi-view attention-based deep neural network that models the relationship between genetic, epigenetic, and transcriptional patterns and identifies co-operative regulatory elements (COREs). We applied this new method, named DeepCORE, to predict transcriptomes in various tissues and cell lines, which outperformed the state-of-the-art algorithms. Furthermore, DeepCORE contains an interpreter that extracts the attention values embedded in the deep neural network, maps the attended regions to putative regulatory elements, and infers COREs based on correlated attentions. The identified COREs are significantly enriched with known promoters and enhancers. Novel regulatory elements discovered by DeepCORE showed epigenetic signatures consistent with the status of histone modification marks.

GUANinE v1.0: Benchmark Datasets for Genomic AI Sequence-to-Function Models

Computational genomics increasingly relies on machine learning methods for genome interpretation, and the recent adoption of neural sequence-to-function models highlights the need for rigorous model specification and controlled evaluation, problems familiar to other fields of AI. Research strategies that have greatly benefited other fields - including benchmarking, auditing, and algorithmic fairness - are also needed to advance the field of genomic AI and to facilitate model development. Here we propose a genomic AI benchmark, GUANinE, for evaluating model generalization across a number of distinct genomic tasks. Compared to existing task formulations in computational genomics, GUANinE is large-scale, de-noised, and suitable for evaluating pretrained models. GUANinE v1.0 primarily focuses on functional genomics tasks such as functional element annotation and gene expression prediction, and it also draws upon connections to evolutionary biology through sequence conservation tasks. The current GUANinE tasks provide insight into the performance of existing genomic AI models and non-neural baselines, with opportunities to be refined, revisited, and broadened as the field matures. Finally, the GUANinE benchmark allows us to evaluate new self-supervised T5 models and explore the tradeoffs between tokenization and model performance, while showcasing the potential for self-supervision to complement existing pretraining procedures.

Coagulation factor VIIa enhances programmed death-ligand 1 expression and its stability in breast cancer cells to promote breast cancer immune evasion

BACKGROUND

Immunotherapy for breast cancer has not gained significant success. Coagulation factor VIIa (FVIIa)-tissue factor (TF) mediated activation of protease-activated receptor 2 (PAR2) is shown to promote metastasis and secretion of the immune-modulatory cytokines but the role of FVIIa in cancer immunology is still not well understood.

OBJECTIVES

Here, we aim to investigate whether FVIIa protects breast cancer cells from CD8 T-cell-mediated killing.

METHODS

Peripheral blood mononuclear cell-derived CD8 T cells were cocultured with vehicle or FVIIa pretreated MDAMB468 cells. The proliferation and activity of CD8 T cells were measured by flow cytometry and ELISA. An allograft model, using wild-type or TF/PAR2-deleted 4T1 cells, was employed to determine the effect of FVIIa on breast cancer immune evasion in vivo.

RESULTS

Here, we demonstrate that TF-FVIIa induces programmed death-ligand 1 (PD-L1) in breast cancer cells by activating PAR2. PAR2 activation triggers large tumor suppressor kinase 1 (LATS1) inactivation leading to loss of yes-associated protein (YAP)/transcriptional coactivator with PDZ-binding motif (TAZ) phosphorylation and subsequent nuclear localization of YAP/TAZ. YAP/TAZ inhibition reduces PD-L1 expression and increases CD8 T-cell activity. We further demonstrate that, apart from transcriptional induction of PD-L1, PAR2 activation also increases PD-L1 stability by enhancing its glycosylation through N-glycosyltransferases STT3A and STT3B.

CONCLUSION

In a mouse model of breast cancer, tumor cell-specific PAR2 depletion leads to PD-L1 downregulation and increases anti-PD-1 immunotherapy efficacy. In conclusion, we showed that FVIIa-mediated signaling cascade in cancer cells serves as a tumor intrinsic mechanism of immunosuppression to promote cancer immune evasion.

Genomic Variations Explorer (GenVarX): a toolset for annotating promoter and CNV regions using genotypic and phenotypic differences

The rapid growth of sequencing technology and its increasing popularity in biology-related research over the years has made whole genome re-sequencing (WGRS) data become widely available. A large amount of WGRS data can unlock the knowledge gap between genomics and phenomics through gaining an understanding of the genomic variations that can lead to phenotype changes. These genomic variations are usually comprised of allele and structural changes in DNA, and these changes can affect the regulatory mechanisms causing changes in gene expression and altering the phenotypes of organisms. In this research work, we created the GenVarX toolset, that is backed by transcription factor binding sequence data in promoter regions, the copy number variations data, SNPs and Indels data, and phenotypes data which can potentially provide insights about phenotypic differences and solve compelling questions in plant research. Analytics-wise, we have developed strategies to better utilize the WGRS data and mine the data using efficient data processing scripts, libraries, tools, and frameworks to create the interactive and visualization-enhanced GenVarX toolset that encompasses both promoter regions and copy number variation analysis components. The main capabilities of the GenVarX toolset are to provide easy-to-use interfaces for users to perform queries, visualize data, and interact with the data. Based on different input windows on the user interface, users can provide inputs corresponding to each field and submit the information as a query. The data returned on the results page is usually displayed in a tabular fashion. In addition, interactive figures are also included in the toolset to facilitate the visualization of statistical results or tool outputs. Currently, the GenVarX toolset supports soybean, rice, and Arabidopsis. The researchers can access the soybean GenVarX toolset from SoyKB via https://soykb.org/SoybeanGenVarX/, rice GenVarX toolset, and Arabidopsis GenVarX toolset from KBCommons web portal with links https://kbcommons.org/system/tools/GenVarX/Osativa and https://kbcommons.org/system/tools/GenVarX/Athaliana, respectively.

Evaluating Gene Expression and Methylation Profiles of TCF4, MBP, and EGR1 in Peripheral Blood of Drug-Free Patients with Schizophrenia: Correlations with Psychopathology, Intelligence, and Cognitive Impairment

Discovery and validation of new, reliable diagnostic and predictive biomarkers for schizophrenia (SCZ) are an ongoing effort. Here, we assessed the mRNA expression and DNA methylation of the TCF4, MBP, and EGR1 genes in the blood of patients with SCZ and evaluated their relationships to psychopathology and cognitive impairments. Quantitative real-time PCR and quantitative methylation-specific PCR methods were used to assess the expression level and promoter DNA methylation status of these genes in 70 drug-free SCZ patients and 72 healthy controls. The correlation of molecular changes with psychopathology and cognitive performance of participants was evaluated. We observed downregulation of TCF4 and upregulation of MBP mRNA levels in SCZ cases, relative to controls in our study. DNA methylation status at the promoter region of TCF4 demonstrated an altered pattern in SCZ as well. Additionally, TCF4 mRNA levels were inversely correlated with PANSS and Stroop total errors and positively correlated with WAIS total score and working memory, consistent with previous studies by our group. In contrast, MBP mRNA level was significantly positively correlated with PANSS and Stroop total errors and inversely correlated with WAIS total score and working memory. These epigenetic and expression signatures can help to assemble a peripheral biomarker-based diagnostic panel for SCZ.

The Influence of KE and EW Dipeptides in the Composition of the Thymalin Drug on Gene Expression and Protein Synthesis Involved in the Pathogenesis of COVID-19

Thymalin is an immunomodulatory drug containing a polypeptide extract of thymus that has demonstrated efficacy in the therapy of acute respiratory distress syndrome and chronic obstructive pulmonary disease, as well as in complex therapy related to severe COVID-19 in middle-aged and elderly patients.. KE and EW dipeptides are active substances of Thymalin. There is evidence that KE stimulates cellular immunity and nonspecific resistance in organisms, exerting an activating effect on macrophages, blood lymphocytes, thymocytes, and neutrophils, while EW reduces angiotensin-induced vasoconstriction and preserves endothelium-dependent vascular relaxation by inhibiting ACE2, the target protein of SARS-CoV-2. However, the mechanism of the immunomodulatory action of Thymalin, KE, and EW during COVID-19 remains unclear. To identify the potential mechanism of action underlying the immunomodulatory activity of Thymalin and its active components, EW and KE dipeptides, we assessed inflammatory response in the context of COVID-19. Interactions between EW and KE dipeptides and double-stranded DNA (dsDNA) were investigated by molecular modeling and docking using ICM-Pro. Analysis of the possible effect of EW and KE dipeptides on gene expression and protein synthesis involved in the pathogenesis of COVID-19 was conducted through the use of bioinformatics methods, including a search for promoter sequences in the Eukaryotic Promoter Database, the determination of genes associated with the development of COVID-19 using the PathCards database of human biological pathways (pathway unification database), identification of the relationship between proteins through cluster analysis in the STRING database ('Search Tool for Retrieval of Interacting Genes/Proteins'), and assessment of the functional enrichment of protein-protein interaction (PPI) using the terms of gene ontology (GO) and the Markov cluster algorithm (MCL). After that, in vitro studying of a lipopolysaccharide (LPS)-induced model of inflammation using human peripheral blood mononuclear cells was performed. ELISA was applied to assess the level of cytokines (IL-1β, IL-6, TNFα) in the supernatant of cells with or without the impact of EW and KE peptides. Blood samples were obtained from four donors; for each cytokine, ELISA was performed 2-4 times, with two parallel experimental or control samples for each experiment (experiments to assess the effects of peptides on LPS-stimulated cells were repeated four times, while additional experiments with unstimulated cells were performed two times). Using molecular docking, GGAG was found to be the best dsDNA sequence in the classical B-form for binding the EW dipeptide, while GCGC is the preferred dsDNA sequence in the curved nucleosomal form for the KE dipeptide. Cluster analysis revealed that potential target genes for the EW and KE peptides encode the AKT1 and AKT2 proteins involved in the development of the cytokine storm. The specific targets for the EW peptide are the ACE2 and CYSLTR1 genes, and specific target for the KE peptide is the CHUK gene. Protein products of the ACE2, CYSLTR1, and CHUK genes are functionally associated with IL-1β, IL-6, TNF-α, IL-4, and IL-10 cytokines. An in vitro model of an inflammatory reaction demonstrated that Thymalin and EW and KE dipeptides reduced the synthesis of IL-1β, IL-6, and TNF-α cytokines in human peripheral blood mononuclear cells by 1.4-6.0 times. The immunomodulatory effect of Thymalin under the inflammatory response conditions in COVID-19 is based on the potential ability of its active components, EW and KE dipeptides, to regulate protein synthesis involved in the development of the cytokine storm.

Investigating the sequence landscape in the Drosophila initiator core promoter element using an enhanced MARZ algorithm

The core promoter elements are important DNA sequences for the regulation of RNA polymerase II transcription in eukaryotic cells. Despite the broad evolutionary conservation of these elements, there is extensive variation in the nucleotide composition of the actual sequences. In this study, we aim to improve our understanding of the complexity of this sequence variation in the TATA box and initiator core promoter elements in Drosophila melanogaster. Using computational approaches, including an enhanced version of our previously developed MARZ algorithm that utilizes gapped nucleotide matrices, several sequence landscape features are uncovered, including an interdependency between the nucleotides in position 2 and 5 in the initiator. Incorporating this information in an expanded MARZ algorithm improves predictive performance for the identification of the initiator element. Overall our results demonstrate the need to carefully consider detailed sequence composition features in core promoter elements in order to make more robust and accurate bioinformatic predictions.

The paracaspase MALT1 is a downstream target of Smad3 and potentiates the crosstalk between TGF-β and NF-kB signaling pathways in cancer cells

TGF-β signaling mediates its biological effects by engaging canonical Smad proteins and crosstalking extensively with other signaling networks, including the NF-kB pathway. The paracaspase MALT1 is an intracellular signaling molecule essential for NF-kB activation downstream of several key cell surface receptors. Despite intensive research on TGF-β and NF-kB interactions, the significance of MALT1 in this context remains undecoded. Here we provide experimental evidence supporting that MALT1 functions to converge these pathways. Using A549 and Huh7 cancer cell line models, we report that TGF-β stimulation enhances MALT1 protein and transcript levels in a time- and dose-dependent manner. Systematic and selective perturbation of TGF-β signaling components identifies MALT1 as a downstream target of Smad3. Rescue experiments in SMAD3 knockout cells confirm that C-terminal phosphorylation of Smad3 is central to MALT1 induction. Corroborating these data, we document that the expression of SMAD3 and MALT1 genes are positively correlated in TCGA cohorts, and we trace the molecular basis of MALT1 elevation to promoter activation. Functional studies in parental as well as NF-kB p65 signaling reporter engineered cells conclusively reveal that MALT1 is paramount for TGF-β-stimulated nuclear translocation and transcriptional activation of NF-kB p65. Furthermore, we find that BCL10 is also implicated in TGF-β activation of NF-kB target genes, potentially coupling the TGF-β-MALT1-NF-kB signaling axis to the CARMA-BCL10-MALT1 (CBM) signalosome. The novel findings of this study indicate that MALT1 is a downstream target of the canonical TGF-β/Smad3 pathway and plays a critical role in modulating TGF-β and NF-kB crosstalk in cancer.

DeepCORE: An interpretable multi-view deep neural network model to detect co-operative regulatory elements

Gene transcription is an essential process involved in all aspects of cellular functions with significant impact on biological traits and diseases. This process is tightly regulated by multiple elements that co-operate to jointly modulate the transcription levels of target genes. To decipher the complicated regulatory network, we present a novel multi-view attention-based deep neural network that models the relationship between genetic, epigenetic, and transcriptional patterns and identifies co-operative regulatory elements (COREs). We applied this new method, named DeepCORE, to predict transcriptomes in 25 different cell lines, which outperformed the state-of-the-art algorithms. Furthermore, DeepCORE translates the attention values embedded in the neural network into interpretable information, including locations of putative regulatory elements and their correlations, which collectively implies COREs. These COREs are significantly enriched with known promoters and enhancers. Novel regulatory elements discovered by DeepCORE showed epigenetic signatures consistent with the status of histone modification marks.

Fish species lifespan prediction from promoter cytosine-phosphate-guanine density

Lifespan is a key attribute of a species' life cycle and varies extensively among major lineages of animals. In fish, lifespan varies by several orders of magnitude, with reported values ranging from less than 1 year to approximately 400 years. Lifespan information is particularly useful for species management, as it can be used to estimate invasion potential, extinction risk and sustainable harvest rates. Despite its utility, lifespan is unknown for most fish species. This is due to the difficulties associated with accurately identifying the oldest individual(s) of a given species, and/or deriving lifespan estimates that are representative for an entire species. Recently it has been shown that CpG density in gene promoter regions can be used to predict lifespan in mammals and other vertebrates, with variable accuracy across taxa. To improve accuracy of lifespan prediction in a non-mammalian vertebrate group, here we develop a fish-specific genomic lifespan predictor. Our new model includes more than eight times the number of fish species included in the previous vertebrate model (n = 442) and uses fish-specific gene promoters as reference sequences. The model predicts fish lifespan from genomic CpG density alone (measured as CpG observed/expected ratio), explaining 64% of the variance between known and predicted lifespans. The predictions are highly robust to variation in genome quality and are applicable to all classes of fish; a taxonomically diverse and speciose group. The results demonstrate the value of promoter CpG density as a universal predictor of fish lifespan that can applied where empirical data are unavailable, or impracticable to obtain.

iPromoter-Seqvec: identifying promoters using bidirectional long short-term memory and sequence-embedded features

Background

Promoters, non-coding DNA sequences located at upstream regions of the transcription start site of genes/gene clusters, are essential regulatory elements for the initiation and regulation of transcriptional processes. Furthermore, identifying promoters in DNA sequences and genomes significantly contributes to discovering entire structures of genes of interest. Therefore, exploration of promoter regions is one of the most imperative topics in molecular genetics and biology. Besides experimental techniques, computational methods have been developed to predict promoters. In this study, we propose iPromoter-Seqvec – an efficient computational model to predict TATA and non-TATA promoters in human and mouse genomes using bidirectional long short-term memory neural networks in combination with sequence-embedded features extracted from input sequences. The promoter and non-promoter sequences were retrieved from the Eukaryotic Promoter database and then were refined to create four benchmark datasets.

Results

The area under the receiver operating characteristic curve (AUCROC) and the area under the precision-recall curve (AUCPR) were used as two key metrics to evaluate model performance. Results on independent test sets showed that iPromoter-Seqvec outperformed other state-of-the-art methods with AUCROC values ranging from 0.85 to 0.99 and AUCPR values ranging from 0.86 to 0.99. Models predicting TATA promoters in both species had slightly higher predictive power compared to those predicting non-TATA promoters. With a novel idea of constructing artificial non-promoter sequences based on promoter sequences, our models were able to learn highly specific characteristics discriminating promoters from non-promoters to improve predictive efficiency.

Conclusions

iPromoter-Seqvec is a stable and robust model for predicting both TATA and non-TATA promoters in human and mouse genomes. Our proposed method was also deployed as an online web server with a user-friendly interface to support research communities. Links to our source codes and web server are available at https://github.com/mldlproject/2022-iPromoter-Seqvec.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08829-6.

Case report: Neonatal diabetes mellitus with congenital hypothyroidism as a result of biallelic heterozygous mutations in GLIS3 gene

Neonatal diabetes mellitus with congenital hypothyroidism (NDH) syndrome (MIM# 610199) is a rare disease caused by autosomal recessive mutations in the GLIS3 gene. GLIS3 is an important transcription factor that might acts as both a repressor and activator of transcription. To date, 22 cases of NDH syndrome from 16 families and 11 countries have been described. Herein, we report a child who developed diabetes during the first week of age. Additionally, she suffered from congenital hypothyroidism, cardiac abnormalities, and polycystic kidney disease. Genetic analysis revealed that patient is a carrier of two novel heterozygous mutations, p.Pro444fsdelG (c.1330delC) and p.His647Arg (c.1940A > G) in the GLIS3 gene. Each was inherited from clinically healthy father and mother, respectively. Bioinformatic tools (SIFT, PolyPhen2, PROVEAN and SWISS-MODEL) declared that the p.His647Arg (c.1940A > G) variant has strong detrimental effect and disturbs Kruppel-like zinc finger domain. All but one so far described cases of NDH syndrome have been caused by homozygous of GLIS3, making the described case the second case of pathogenic, compound heterozygosity of GLIS3 worldwide posing substantial clinical novelty and detailing an interesting interplay between the observed variants and GLIS3 expression, which seems to be autoregulated. Hence, the damaging missense mutation may further reduce the expression of any remaining functional alleles. This case report expands our understanding of the clinical phenotype, treatment approaches, and outcome of infants with GLIS3 mutations and indicates the need for further research to deepen our understanding of the role of GLIS3.

Genome-wide DNA methylation profiling of stomach cancer in the ethnic population of Mizoram, North East India

Stomach cancer is the fifth most common cancer in terms of prevalence and incidence and the fourth leading cause of mortality in men and women worldwide. It is well-established that aberrant DNA methylation in cells can lead to carcinogenesis. The primary objective of our study was to investigate the aberrant DNA methylation status of genes associated with stomach cancer with a particular reference to the ethnic population of Mizoram, North East India. The site-level analysis identified 2883 CpG sites differentially methylated, representing ~922 genes. Out of which 476 Differentially Methylated Positions (DMPs) were promoter-associated, 452 DMPs were hypermethylated, and 24 were hypomethylated. The region-level analysis identified 462 Differentially Methylated Regions (DMRs) corresponding to ~320 genes, of which ~281 genes were hypermethylated and ~ 40 genes were hypomethylated. TCGA analysis showed that some of the genes had been previously implicated in other cancers including stomach cancer. Five hypermethylated genes were selected as candidate genes for further investigations and they have shown to be novel and could serve as candidate hypermethylation biomarkers for stomach cancer in this particular ethnic group.

Therapeutic application of PPE2 protein of Mycobacterium tuberculosis in inhibiting tissue inflammation

There is an increasing need to develop biological anti-inflammatory agents that are more targeted, effective, and with lesser side effects as compared to conventional chemical drugs. In the present study, we found that Mycobacterium tuberculosis protein PPE2 and a synthetic derivative peptide can suppress the mast cell population and inhibit several vasoactive and fibrogenic mediators and pro-inflammatory cytokines induced by mast cells in formalin-induced tissue injury. PPE2 was found to inhibit transcription from the promoter of stem cell factor, important for mast cell maintenance and migration. Thus, PPE2/peptide can be used as a potent nonsteroidal therapeutic agent for the treatment of inflammation and tissue injury.

Critical assessment of computational tools for prokaryotic and eukaryotic promoter prediction

Promoters are crucial regulatory DNA regions for gene transcriptional activation. Rapid advances in next-generation sequencing technologies have accelerated the accumulation of genome sequences, providing increased training data to inform computational approaches for both prokaryotic and eukaryotic promoter prediction. However, it remains a significant challenge to accurately identify species-specific promoter sequences using computational approaches. To advance computational support for promoter prediction, in this study, we curated 58 comprehensive, up-to-date, benchmark datasets for 7 different species (i.e. Escherichia coli, Bacillus subtilis, Homo sapiens, Mus musculus, Arabidopsis thaliana, Zea mays and Drosophila melanogaster) to assist the research community to assess the relative functionality of alternative approaches and support future research on both prokaryotic and eukaryotic promoters. We revisited 106 predictors published since 2000 for promoter identification (40 for prokaryotic promoter, 61 for eukaryotic promoter, and 5 for both). We systematically evaluated their training datasets, computational methodologies, calculated features, performance and software usability. On the basis of these benchmark datasets, we benchmarked 19 predictors with functioning webservers/local tools and assessed their prediction performance. We found that deep learning and traditional machine learning-based approaches generally outperformed scoring function-based approaches. Taken together, the curated benchmark dataset repository and the benchmarking analysis in this study serve to inform the design and implementation of computational approaches for promoter prediction and facilitate more rigorous comparison of new techniques in the future.

lncRNAfunc: a knowledgebase of lncRNA function in human cancer

The long non-coding RNAs associating with other molecules can coordinate several physiological processes and their dysfunction can impact diverse human diseases. To date, systematic and intensive annotations on diverse interaction regulations of lncRNAs in human cancer were not available. Here, we built lncRNAfunc, a knowledgebase of lncRNA function in human cancer at https://ccsm.uth.edu/lncRNAfunc, aiming to provide a resource and reference for providing therapeutically targetable lncRNAs and intensive interaction regulations. To do this, we collected 15 900 lncRNAs across 33 cancer types from TCGA. For individual lncRNAs, we performed multiple interaction analyses of different biomolecules including DNA, RNA, and protein levels. Our intensive studies of lncRNAs provide diverse potential mechanisms of lncRNAs that regulate gene expression through binding enhancers and 3'-UTRs of genes, competing for miRNA binding sites with mRNAs, recruiting the transcription factors to gene promoters. Furthermore, we investigated lncRNAs that potentially affect the alternative splicing events through interacting with RNA binding Proteins. We also performed multiple functional annotations including cancer stage-associated lncRNAs, RNA A-to-I editing event-associated lncRNAs, and lncRNA expression quantitative trait loci. lncRNAfunc is a unique resource for cancer research communities to help better understand potential lncRNA regulations and therapeutic lncRNA targets.

The -592C>A Variation of IL-10 Gene and Susceptibility to Chronic Periodontitis: A Genetic Association Study and In-Silico Analysis

OBJECTIVES

Chronic periodontitis (CP) is a common inflammatory disorder with a considerable impact of genetic variations in the interleukin family on predisposition to this disease. This study aimed to investigate the association between the -592C>A polymorphism of the interleukin 10 (IL-10) gene with CP risk in an Iranian population. This experimental study was followed by a meta-analysis and in silico examination.

METHODS

In a case-control study, 270 subjects, including 135 patients with CP and 135 healthy controls, were enrolled. The -592C>A genotyping was performed using the PCR-RFLP method. In the meta-analysis, valid databases were systematically searched to identify relevant studies. Pooled odds ratios (ORs) with 95% confidence intervals (CIs) were computed to examine the association between -592C>A and CP. In silico analysis was conducted using bioinformatics tools to evaluate the impact of the -592C>A polymorphism on IL-10 gene function.

RESULTS

Our case-control study revealed a significant association between polymorphism and CP risk. Overall, we found significant associations between -592C>A genetic variation and CP and stratified meta-analysis. In silico analysis revealed that this polymorphism could change the pattern of the transcription binding site upstream of the IL-10 gene. It may also alter the hsa-miR-101-3p miRNA-targeted sequence upstream of IL-10.

CONCLUSIONS

Based on our results, the -592C>A variation in IL-10 may be a genetic risk factor for susceptibility to chronic periodontitis. However, further studies in different ethnicities and results adjusted for clinical and demographic characteristics are needed to obtain more accurate deductions.

Relationship between Functional miR-143/145 Cluster Variants and Susceptibility to Type 2 Diabetes Mellitus: A Preliminary Case-Control Study and Bioinformatics Analyses

Purpose: To investigate the link between two variants (rs4705342 and rs4705343) in the promoter of the miR-143/145 cluster with Type 2 diabetes mellitus (T2DM) risk. Methods:A total of 1200 subjects were genotyped using the ARMS-PCR method. Results: The rs4705342 variant enhanced the risk of T2DM under codominant CC (OR = 3.24; 95% CI: 1.89-5.60), recessive TT+TC (OR = 3.02; 95% CI: 1.77-5.17), and dominant TC+CC (OR = 1.35; 95% CI: 1.08-1.71) genetic models. Individuals carrying the C allele of rs4705342 conferred a 1.43 fold increased risk of T2DM. As regards rs4705343, decreased risk of T2DM was observed under codominant TC (OR = 0.53; 95% CI: 0.42-0.67), over-dominant TT+CC (OR = 0.51; 95% CI: 0.40-0.64), and dominant TC+CC (OR = 0.59; 95% CI: 0.48-0.75) models. Haplotype analysis of the variants showed a 1.941-fold increased risk of T2DM regarding the C T combination. Significant associations were noticed between different haplotypes and lipid indices of T2DM patients. There were no notable changes in p-values after adjustment for BMI. Computational analysis revealed that miR143 and/or miR145 target important genes involved in glucose and lipid metabolism. Conclusions: Functional miR-143/145 variants might influence the risk of T2DM. Hence, clarifying the precise regulatory mechanisms of gene expression in the development of T2DM will significantly guide researchers to find a novel target for therapeutic intervention.

BioMaster: An Integrated Database and Analytic Platform to Provide Comprehensive Information About BioBrick Parts

Synthetic biology seeks to create new biological parts, devices, and systems, and to reconfigure existing natural biological systems for custom-designed purposes. The standardized BioBrick parts are the foundation of synthetic biology. The incomplete and flawed metadata of BioBrick parts, however, are a major obstacle for designing genetic circuit easily, quickly, and accurately. Here, a database termed BioMaster http://www.biomaster-uestc.cn was developed to extensively complement information about BioBrick parts, which includes 47,934 items of BioBrick parts from the international Genetically Engineered Machine (iGEM) Registry with more comprehensive information integrated from 10 databases, providing corresponding information about functions, activities, interactions, and related literature. Moreover, BioMaster is also a user-friendly platform for retrieval and analyses of relevant information on BioBrick parts.

Differential Methylation in Promoter Regions of the Genes NR3C1 and HSP90AA1, Involved in the Regulation, and Bioavailability of Cortisol in Leukocytes of Women With Preeclampsia

Introduction: Hypertensive disorders are of interest in obstetrics and gynecology because they are the second place among causes of maternal mortality and a source of complications in the short, mid, and long term. Even if the pathophysiological process behind preeclampsia (PE) is still unknown, stress factors have been revealed to play an important role in the genesis of this pathologic process.

Methods: A case-control study was designed with the purpose of determining if there is a differential methylation in NR3C1, HSD11B2, CYP11A1, CRHBP, TEAD3, and HSP90AA1 genes, related to signaling of the hypothalamic–pituitary–adrenal axis, and its regulation on early-onset PE (EOPE).

Results: A total of 20 cases and 20 controls were studied by DNA methylation analysis, demonstrating differences among groups in the percentage of methylation of the NR3C1 gene. After a contingency analysis, an odds ratio (OR) for PE of 12.25 was identified for NR3C1 and 9.9 for HSP90AA1 genes. NR3C1, TEAD3, and HSP90AA1 genes showed a positive correlation with the systolic and diastolic blood pressure levels with a p ≤ 0.05.

Conclusion: This study found a differential methylation in the glucocorticoid receptor (GR) NR3C1 and its co-chaperone HSP90AA1 in women with PE, with a possible regulatory role in the response to stress in pregnancy and is a likely physiopathological mechanism in PE.

Antibiotic drug discovery: Challenges and perspectives in the light of emerging antibiotic resistance

Amid a rising threat of antimicrobial resistance in a global scenario, our huge investments and high-throughput technologies injected for rejuvenating the key therapeutic scaffolds to suppress these rising superbugs has been diminishing severely. This has grasped world-wide attention, with increased consideration being given to the discovery of new chemical entities. Research has now proven that the relatively tiny and simpler microbes possess enhanced capability of generating novel and diverse chemical constituents with huge therapeutic leads. The usage of these beneficial organisms could help in producing new chemical scaffolds that govern the power to suppress the spread of obnoxious superbugs. Here in this review, we have explicitly focused on several appealing strategies employed for the generation of new chemical scaffolds. Also, efforts on providing novel insights on some of the unresolved questions in the production of metabolites, metabolic profiling and also the serendipity of getting "hit molecules" have been rigorously discussed. However, we are highly aware that biosynthetic pathway of different classes of secondary metabolites and their biosynthetic route is a vast topic, thus we have avoided discussion on this topic.

Genetic and epigenetic analysis of the BAX and BCL2 in the placenta of pregnant women complicated by preeclampsia

The current study examined the effects of BAX and BCL2 polymorphisms and methylation as well as mRNA expression on susceptibility to PE. After delivery, the placentas were collected from 92 women with PE, as well as 106 normotensive pregnant women. The BAX rs4645878 and BCL2 rs2279115 polymorphisms were genotyped by the PCR-RFLP method. Methylation-specific PCR (MSP) was used for analysis of promoter methylation. mRNA expression was assayed by Quantitative RT-PCR. In addition, in silico analysis was performed by bioinformatics tools. There was no relationship between PE and placental BAX rs4645878 and BCL2 rs2279115 polymorphisms. The groups were not significantly different regarding the promoter methylation of BAX gene. Nonetheless, the MM status of BCL2 promoter had a significantly higher frequency in the PE group and was associated with 2.7-fold higher risk of PE (OR = 2.7, 95% CI = 1.3-5.6; P = 0.01). The relative mRNA expression of BCL2 was decreased in the placentas of PE women (P < 0.0001). The expression of BAX gene was not significantly different between the two groups. There was no association between placental BAX rs4645878 and BCL2 rs2279115 polymorphisms and mRNA expression levels. In silico analysis indicated that BAX rs4645878 and BCL2 rs2279115 polymorphisms were located in the core recognition site of different transcription factors and these substitutions of wild allele resulted in the loss and/ or change of these binding sites and subsequently may alter BCL2 and BAX expression. This study showed that the BAX and BCL2 polymorphisms and BAX promoter methylation were not associated with PE risk. The BCL2 promoter methylation was associated with lower BCL2 expression and higher PE susceptibility.

Characterization of the CCAAT-binding transcription factor complex in the plant pathogenic fungus Fusarium graminearum

The CCAAT sequence is a ubiquitous cis-element of eukaryotic promoters, and genes containing CCAAT sequences have been shown to be activated by the CCAAT-binding transcription factor complex in several eukaryotic model organisms. In general, CCAAT-binding transcription factors form heterodimers or heterotrimeric complexes that bind to CCAAT sequences within the promoters of target genes and regulate various cellular processes. To date, except Hap complex, CCAAT-binding complex has been rarely reported in fungi. In this study, we characterized two CCAAT-binding transcription factors (Fct1 and Fct2) in the plant pathogenic fungus Fusarium graminearum. Previously, FCT1 and FCT2 were shown to be related to DNA damage response among eight CCAAT-binding transcription factors in F. graminearum. We demonstrate that the nuclear CCAAT-binding complex of F. graminearum has important functions in various fungal developmental processes, not just DNA damage response but virulence and mycotoxin production. Moreover, the results of biochemical and genetic analyses revealed that Fct1 and Fct2 may form a complex and play distinct roles among the eight CCAAT-binding transcription factors encoded by F. graminearum. To the best of our knowledge, the results of this study represent a substantial advancement in our understanding of the molecular mechanisms underlying the functions of CCAAT-binding factors in eukaryotes.

Association of A-197G polymorphism in interleukin-17 gene with chronic periodontitis: Evidence from six case-control studies with a computational biology approach

AIM

The aim of the present study was to evaluate the association of interleukin-17 (IL-17) A-197G gene polymorphism with chronic periodontitis (CP) in a case-control study, a meta-analysis, and an in silico approach.

METHODS

In the case-control study, 122 cases with CP and 126 healthy controls were recruited; IL-17 A-197G genotyping was performed by polymerase chain reaction-restriction fragment length polymorphism. In the meta-analysis, comprehensive literature retrieval was performed on valid databases to identify relevant studies. Bioinformatics tools were employed to investigate the effects of A-197G transition on the promoter region of IL-17.

RESULTS

Our case-control study revealed a significant association between IL-17 A-197G transition and CP. The overall meta-analysis revealed significant associations between the IL-17 A-197G polymorphism and CP risk in homozygote co-dominant and recessive models. The stratified analysis also showed a statistically significant association between the mentioned transition and CP risk in the Caucasian population. The in silico analysis revealed that the A-197G polymorphism could make changes in protein-binding sites of the IL-17 promoter region.

CONCLUSIONS

Our study supports that IL-17 A-197G transition could be a genetic risk factor for CP. However, further studies with a larger sample size among different ethnicities are required to obtain a more accurate conclusion.

DeePromoter: Robust Promoter Predictor Using Deep Learning

The promoter region is located near the transcription start sites and regulates transcription initiation of the gene by controlling the binding of RNA polymerase. Thus, promoter region recognition is an important area of interest in the field of bioinformatics. Numerous tools for promoter prediction were proposed. However, the reliability of these tools still needs to be improved. In this work, we propose a robust deep learning model, called DeePromoter, to analyze the characteristics of the short eukaryotic promoter sequences, and accurately recognize the human and mouse promoter sequences. DeePromoter combines a convolutional neural network (CNN) and a long short-term memory (LSTM). Additionally, instead of using non-promoter regions of the genome as a negative set, we derive a more challenging negative set from the promoter sequences. The proposed negative set reconstruction method improves the discrimination ability and significantly reduces the number of false positive predictions. Consequently, DeePromoter outperforms the previously proposed promoter prediction tools. In addition, a web-server for promoter prediction is developed based on the proposed methods and made available at https://home.jbnu.ac.kr/NSCL/deepromoter.htm.

Aberrant expression of the microtubule-associated protein tau is an independent prognostic feature in prostate cancer

Background

Microtubule-associated protein Tau (MAPT) overexpression has been linked to poor prognosis and decreased response to taxane-based therapies in several cancer types, but its relevance in prostate cancer is unknown.

Methods

In this study, MAPT expression was analyzed by immunohistochemistry on a tissue microarray containing 17,747 prostate cancers.

Results

MAPT was absent in normal prostate epithelial cells but detectable in 1004 (8.2%) of 12,313 interpretable cancers. Its expression was associated with advanced tumor stage, high Gleason grade, positive lymph nodes, and early biochemical recurrence (p < 0.0001 each). For example, MAPT was found in 3.6% of 2072 Gleason ≤3 + 3 cancers but in 14.4% of 704 Gleason ≥4 + 4 cancers. High-level MAPT staining was also linked to TMPRSS2:ERG fusions (p < 0.0001). MAPT staining was seen in 15.2 and 16% of cancers with TMPRSS2:ERG fusion detected by immunohistochemistry and fluorescence in-situ hybridization, but in only 3.5 and 3.9% of cancers without ERG staining or ERG rearrangements. Moreover, an association was found between MAPT expression and PTEN deletions, with 19% MAPT positivity in 948 PTEN deleted cancers but only 7% MAPT positivity in 3895 tumors with normal PTEN copy numbers (p < 0.0001). Multivariate analysis revealed that the prognostic value of MAPT was independent from established parameters. Conventional large section analyses showed intratumoral MAPT heterogeneity in all three analyzed cancers.

Conclusions

The results of our study identify MAPT, as a moderate prognostic marker in prostate cancer, whose clinical impact, however, may be limited due to the rarity and heterogeneity of its expression.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5390-1) contains supplementary material, which is available to authorized users.

Prediction of Human Immunodeficiency Virus Type 1 Subtype-Specific Off-Target Effects Arising from CRISPR-Cas9 Gene Editing Therapy

Chronic human immunodeficiency virus type 1 (HIV-1) disease is characterized by the retention of provirus within latently infected cells. Anti-HIV-1 CRISPR-Cas9 gene editing is an attractive strategy to excise or inactivate the HIV-1 genome. Recent strategies have focused on designing gRNAs that target the long terminal repeat (LTR) because 5' and 3' LTR symmetry can facilitate proviral excision. However, the promiscuity of CRISPR-Cas9 gene editing system necessitates the investigation of potential off-target effects. Here, potential gRNAs designed from HIV-1 phylogenetic subtypes using the CRISPRseek tool were investigated. Across the LTR, it was found that certain regions show higher human homology than others. When using recommended cutoffs, 96.40% of gRNAs were predicted to have no high probability off-target effects. Given this observation, while high-probability off-target effects are a potential danger, they can be avoided with proper gRNA design.

Genetic polymorphisms and haplotypes of the DJ-1 gene promoter associated with the susceptibility to male infertility

PURPOSE

In this study, we evaluate the relationship between genetic polymorphisms of the DJ-1 gene, g.-6_+10del, and g.168_185del with male infertility susceptibility.

METHODS

Four hundred and twenty-two male infertile patients and 285 fertile male controls were recruited. Genotyping was performed by polymerase chain reaction. In silico analysis was performed by EPD, ElemeNT, SNPnexus, and PROMO to predict the potential functions of rs901561484 and rs373653682 polymorphisms.

RESULTS

The Del (D) allele carriers of DJ-1 g.-6_+10del polymorphism were significantly associated with the risk of male infertility in total infertile, asthenozoospermia, and oligoasthenozoospermia patients. Moreover, the Del (D) allele of DJ-1 g.-6_+10del polymorphism significantly increased in total male infertile, asthenozoospermia, and oligoasthenozoospermia groups. In addition, the frequencies of different genotypes and the Del allele and Dup allele carriers of DJ-1 g.168_185del gene polymorphisms were associated with male infertility in total infertile and four different sub-group patients. Furthermore, haplotype analysis of DJ-1 g.-6_+10del and g.168_185del polymorphisms revealed that the D-Dup and I-Del haplotype frequencies significantly increased the risk of male infertility, while I-Ins haplotypes were associated with a decreased risk of male infertility in total and sub-group patients. The in silico analysis showed that the presence of Ins and/or Dup alleles of the DJ-1 g.-6_+10del and g.168_185del polymorphisms could provide additional binding sites of more nuclear factors and probably affect transcriptional activity.

CONCLUSION(S)

Our study presents evidence of a strong association between functional polymorphisms of the DJ-1 promoter, g.-6_+10del, and g.168_185del with the risk of male infertility.

XRCC5 VNTR, XRCC6 -61C>G, and XRCC7 6721G>T Gene Polymorphisms Associated with Male Infertility Risk: Evidences from Case-Control and In Silico Studies

We evaluate the association between genetic polymorphisms of XRCC5 VNTR, XRCC6 -61C>G, and XRCC7 6721G>T with male infertility susceptibility. A total of 392 men including 178 infertile males (102 idiopathic azoospermia and 76 severe oligozoospermia) and 214 healthy controls were recruited. XRCC6 -61C>G and XRCC7 6721G>T genotyping was performed by PCR-RFLP whereas XRCC5 VNTR was performed by PCR. The 2R allele and 2R allele carriers of XRCC5 VNTR polymorphism significantly decreased risk of male infertility. The mutant GG genotypes and carriers of the CG and GG genotypes of XRCC6 -61C>G showed increased risk for the male infertility. Furthermore, the G allele of the XRCC6 -61C>G was correlated with increased susceptibility to male infertility. Likewise, the T allele of the XRCC7 6721G>T polymorphism was associated with increased susceptibility to male infertility in azoospermia. In silico analysis predicted that the presence of tandem repeats in XRCC5 gene prompter can be sequence to bind to more nuclear factors. Also, rs2267437 (C>G) variant was located in a well-conserved region in XRCC6 promoter and this variation might lead to differential allelic expression. The XRCC7 6721G>T gene polymorphism occurred in an acceptor-splicing site, but this polymorphism has no severe modification on XRCC7 mRNA splicing. Our results indicate the association of XRCC5 VNTR, XRCC6 -61C>G, and XRCC7 6721G>T gene polymorphisms with male infertility in Iranian men.

Deciphering the molecular mechanisms behind cellulase production in Trichoderma reesei, the hyper-cellulolytic filamentous fungus

The filamentous fungus Trichoderma reesei is a potent cellulase producer and the best-studied cellulolytic fungus. A lot of investigations not only on glycoside hydrolases produced by T. reesei, but also on the machinery controlling gene expression of these enzyme have made this fungus a model organism for cellulolytic fungi. We have investigated the T. reesei strain including mutants developed in Japan in detail to understand the molecular mechanisms that control the cellulase gene expression, the biochemical and morphological aspects that could favor this phenotype, and have attempted to generate novel strains that may be appropriate for industrial use. Subsequently, we developed recombinant strains by combination of these insights and the heterologous-efficient saccharifing enzymes. Resulting enzyme preparations were highly effective for saccharification of various biomass. In this review, we present some of the salient findings from the recent biochemical, morphological, and molecular analyses of this remarkable cellulase hyper-producing fungus.

Recognition Code of ZNF191(243-368) and Its Interaction with DNA

ZNF191(243-368) is the C-terminal region of ZNF191 which contains a putative DNA-binding domain of four Cys2His2 zinc finger motifs. In this study, an expression vector of a fusion protein of ZNF191(243-368) with glutathione-S-transferase (GST) was constructed and transformed into Escherichia coli BL21. The fusion protein GST-ZNF191(243-368) was expressed using this vector to investigate the protein-DNA binding reaction through an affinity selection strategy on the basis of the binding quality of the zinc finger domain. Results showed that ZNF191(243-368) can selectively bind with sequences and react with genes which contain an AGGG core. However, the recognition mechanism of Cys2His2 zinc finger proteins to DNA warrants further investigation.

Genome-wide analysis of tandem repeats in plants and green algae

Tandem repeats (TRs) extensively exist in the genomes of prokaryotes and eukaryotes. Based on the sequenced genomes and gene annotations of 31 plant and algal species in Phytozome version 8.0 (http://www.phytozome.net/), we examined TRs in a genome-wide scale, characterized their distributions and motif features, and explored their putative biological functions. Among the 31 species, no significant correlation was detected between the TR density and genome size. Interestingly, green alga Chlamydomonas reinhardtii (42,059 bp/Mbp) and castor bean Ricinus communis (55,454 bp/Mbp) showed much higher TR densities than all other species (13,209 bp/Mbp on average). In the 29 land plants, including 22 dicots, 5 monocots, and 2 bryophytes, 5′-UTR and upstream intergenic 200-nt (UI200) regions had the first and second highest TR densities, whereas in the two green algae (C. reinhardtii and Volvox carteri) the first and second highest densities were found in intron and coding sequence (CDS) regions, respectively. In CDS regions, trinucleotide and hexanucleotide motifs were those most frequently represented in all species. In intron regions, especially in the two green algae, significantly more TRs were detected near the intron–exon junctions. Within intergenic regions in dicots and monocots, more TRs were found near both the 5′ and 3′ ends of genes. GO annotation in two green algae revealed that the genes with TRs in introns are significantly involved in transcriptional and translational processing. As the first systematic examination of TRs in plant and green algal genomes, our study showed that TRs displayed nonrandom distribution for both intragenic and intergenic regions, suggesting that they have potential roles in transcriptional or translational regulation in plants and green algae.

Transcriptional regulation of immediate-early gene response by THOC5, a member of mRNA export complex, contributes to the M-CSF-induced macrophage differentiation

Hematopoiesis and commitment to a restricted lineage are guided by a timely expressed set of cytokine receptors and their downstream transcription factors. A member of the mRNA export complex, THOC5 (suppressors of the transcriptional defects of hpr1 delta by overexpression complex 5) is a substrate for several tyrosine kinases such as macrophage colony-stimulating factor (M-CSF) receptor and various leukemogenic tyrosine kinases, such as Bcr-Abl, or NPM-ALK. THOC5 tyrosine phosphorylation is elevated in stem cells from patients with chronic myeloid leukemia, suggesting that THOC5 may be involved in leukemia development. THOC5 is also an essential element in the maintenance of hematopoiesis in adult mice. In this report, we show that THOC5 is located in the nuclear speckles, and that it is translocated from the nucleus to cytoplasm during M-CSF-induced bone marrow-derived macrophage differentiation. Furthermore, we have identified THOC5 target genes by trancriptome analysis, using tamoxifen-inducible THOC5 knockout macrophages. Although only 99 genes were downregulated in THOC5-depleted macrophages, half of the genes are involved in differentiation and/or migration. These include well-known regulators of myeloid differentiation inhibitor of DNA binding (Id)1, Id3, Smad family member 6 (Smad6) and Homeobox (Hox)A1. In addition, a subset of M-CSF-inducible genes, such as Ets family mRNAs are THOC5 target mRNAs. Upon depletion of THOC5, unspliced v-ets erythroblastosis virus E26 oncogene homolog (Ets1) mRNA was accumulated in the nucleus. Furthermore, THOC5 was recruited to chromatin where Ets1 was transcribed and bound to unspliced and spliced Ets1 transcripts, indicating that THOC5 has a role in processing/export of M-CSF-inducible genes. In conclusion, regulation of immediate-early gene response by THOC5, a member of mRNA export complex contributes to the M-CSF-induced macrophage differentiation.

Computational tools and resources for prediction and analysis of gene regulatory regions in the chick genome

The discovery of cis-regulatory elements is a challenging problem in bioinformatics, owing to distal locations and context-specific roles of these elements in controlling gene regulation. Here we review the current bioinformatics methodologies and resources available for systematic discovery of cis-acting regulatory elements and conserved transcription factor binding sites in the chick genome. In addition, we propose and make available, a novel workflow using computational tools that integrate CTCF analysis to predict putative insulator elements, enhancer prediction, and TFBS analysis. To demonstrate the usefulness of this computational workflow, we then use it to analyze the locus of the gene Sox2 whose developmental expression is known to be controlled by a complex array of cis-acting regulatory elements. The workflow accurately predicts most of the experimentally verified elements along with some that have not yet been discovered. A web version of the CTCF tool, together with instructions for using the workflow can be accessed from http://toolshed.g2.bx.psu.edu/view/mkhan1980/ctcf_analysis. For local installation of the tool, relevant Perl scripts and instructions are provided in the directory named "code" in the supplementary materials.

Gene promoters show chromosome-specificity and reveal chromosome territories in humans

Background

Gene promoters have guided evolution processes for millions of years. It seems that they were the main engine responsible for the integration of different mutations favorable for the environmental conditions. In cooperation with different transcription factors and other biochemical components, these regulatory regions dictate the synthesis frequency of RNA molecules. Predominantly in the last decade, it has become clear that nuclear organization impacts upon gene regulation. To fully understand the connections between Homo sapiens chromosomes and their gene promoters, we analyzed 1200 promoter sequences using our Kappa Index of Coincidence method.

Results

In order to measure the structural similarity of gene promoters, we used two-dimensional image-based patterns obtained through Kappa Index of Coincidence (Kappa IC) and (C+G)% values. The center of weight of each promoter pattern indicated a structure similarity between promoters of each chromosome. Furthermore, the proximity of chromosomes seems to be in accordance to the structural similarity of their gene promoters. The arrangement of chromosomes according to Kappa IC values of promoters, shows a striking symmetry between the chromosome length and the structure of promoters located on them. High Kappa IC and (C+G)% values of gene promoters were also directly associated with the most frequent genetic diseases. Taking into consideration these observations, a general hypothesis for the evolutionary dynamics of the genome has been proposed. In this hypothesis, heterochromatin and euchromatin domains exchange DNA sequences according to a difference in the rate of Slipped Strand Mispairing and point mutations.

Conclusions

In this paper we showed that gene promoters appear to be specific to each chromosome. Furthermore, the proximity between chromosomes seems to be in accordance to the structural similarity of their gene promoters. Our findings are based on comprehensive data from Transcriptional Regulatory Element Database and a new computer model whose core is using Kappa index of coincidence.

Eukaryotic genomes may exhibit up to 10 generic classes of gene promoters

BACKGROUND

The main function of gene promoters appears to be the integration of different gene products in their biological pathways in order to maintain homeostasis. Generally, promoters have been classified in two major classes, namely TATA and CpG. Nevertheless, many genes using the same combinatorial formation of transcription factors have different gene expression patterns. Accordingly, we tried to ask ourselves some fundamental questions: Why certain genes have an overall predisposition for higher gene expression levels than others? What causes such a predisposition? Is there a structural relationship of these sequences in different tissues? Is there a strong phylogenetic relationship between promoters of closely related species?

RESULTS

In order to gain valuable insights into different promoter regions, we obtained a series of image-based patterns which allowed us to identify 10 generic classes of promoters. A comprehensive analysis was undertaken for promoter sequences from Arabidopsis thaliana, Drosophila melanogaster, Homo sapiens and Oryza sativa, and a more extensive analysis of tissue-specific promoters in humans. We observed a clear preference for these species to use certain classes of promoters for specific biological processes. Moreover, in humans, we found that different tissues use distinct classes of promoters, reflecting an emerging promoter network. Depending on the tissue type, comparisons made between these classes of promoters reveal a complementarity between their patterns whereas some other classes of promoters have been observed to occur in competition. Furthermore, we also noticed the existence of some transitional states between these classes of promoters that may explain certain evolutionary mechanisms, which suggest a possible predisposition for specific levels of gene expression and perhaps for a different number of factors responsible for triggering gene expression. Our conclusions are based on comprehensive data from three different databases and a new computer model whose core is using Kappa index of coincidence.

CONCLUSIONS

To fully understand the connections between gene promoters and gene expression, we analyzed thousands of promoter sequences using our Kappa Index of Coincidence method and a specialized Optical Character Recognition (OCR) neural network. Under our criteria, 10 classes of promoters were detected. In addition, the existence of "transitional" promoters suggests that there is an evolutionary weighted continuum between classes, depending perhaps upon changes in their gene products.

Cryptic transcripts from a ubiquitous plasmid origin of replication confound tests for cis-regulatory function

A vast amount of research on the regulation of gene expression has relied on plasmid reporter assays. In this study, we show that plasmids widely used for this purpose constitutively produce substantial amounts of RNA from a TATA-containing cryptic promoter within the origin of replication. Readthrough of these RNAs into the intended transcriptional unit potently stimulated reporter activity when the inserted test sequence contained a 3′ splice site (ss). We show that two human sequences, originally reported to be internal ribosome entry sites and later to instead be promoters, mimic both types of element in dicistronic reporter assays by causing these cryptic readthrough transcripts to splice in patterns that allow efficient translation of the downstream cistron. Introduction of test sequences containing 3′ ss into monocistronic luciferase reporter vectors widely used in the study of transcriptional regulation also created the false appearance of promoter function via the same mechanism. Across a large number of variants of these plasmids, we found a very highly significant correlation between reporter activity and levels of such spliced readthrough transcripts. Computational estimation of the frequency of cryptic 3′ ss in genomic sequences suggests that misattribution of cis-regulatory function may be a common occurrence.

The development of RAD51C, Cystatin A, p53 and Nrf2 luciferase-reporter assays in metabolically competent HepG2 cells for the assessment of mechanism-based genotoxicity and of oxidative stress in the early research phase of drug development

Four different mechanism-based high-throughput luciferase-reporter assays were developed in human HepG2 cells, which contain phase I and II metabolic activity and a functionally active p53 protein. The promoter regions of RAD51C and Cystatin A, as well as the responsive element of the p53 protein, were selected for the generation of the genotoxicity reporter assays. Moreover, a luciferase-based reporter assay was generated that measures the activation of the Nrf2 oxidative stress pathway. Validation with respect to the ECVAM compound list [D. Kirkland, P. Kasper, L. Muller, R. Corvi, G. Speit, Recommended lists of genotoxic and non-genotoxic chemicals for assessment of the performance of new or improved genotoxicity tests: a follow-up to an ECVAM workshop, Mutat. Res. 653 (2008) 99-108] resulted in an overall sensitivity of the HepG2 genotoxicity reporter assays for genotoxicity of 85% (17/20). The specificity and predictivity were high with 81% (34/42) and 82% (51/62), respectively. Various compounds had a positive score although metabolic activation was needed. The HepG2 reporter data were also compared with the available data on bacterial mutagenicity (Ames test), in vitro clastogenicity and in vivo clastogenicity for an additional set of 192 compounds. The predictivity for mutagenicity results was 74% (sensitivity, 61%, 30/49; specificity, 80%, 77/96) and for in vitro clastogenicity 59% (sensitivity, 45%, 35/78; specificity 83%, 38/46). The correlation between results from the HepG2 genotoxicity reporter assays and in vivo clastogenicity was much higher with 77% (sensitivity, 74%, 28/38; specificity 81%, 26/32). Results from the Nrf2 reporter assay showed that a large number of genotoxic compounds activated the Nrf2 oxidative stress pathway. In conclusion, four high-throughput mechanism-based reporter assays in the HepG2 cell line were developed, which can be applied for screening in the early research phase of drug development. The use of these assays in combination with the previously validated Vitotox and RadarScreen assays will certainly reduce the attrition rate due to genotoxicity in the developmental phase of drug development.

Epithelial cells in fetal intestine produce chemerin to recruit macrophages

Macrophages are first seen in the fetal intestine at 11-12 wk and rapidly increase in number during the 12- to 22-wk period of gestation. The development of macrophage populations in the fetal intestine precedes the appearance of lymphocytes and neutrophils and does not require the presence of dietary or microbial antigens. In this study, we investigated the role of chemerin, a recently discovered, relatively selective chemoattractant for macrophages, in the recruitment of macrophage precursors to the fetal intestine. Chemerin mRNA/protein expression was measured in jejunoileal tissue from 10- to 24-wk human fetuses, neonates operated for intestinal obstruction, and adults undergoing bariatric surgery. The expression of chemerin in intestinal epithelial cells (IECs) was confirmed by using cultured primary IECs and IEC-like cell lines in vitro. The regulatory mechanisms involved in chemerin expression were investigated by in silico and immunolocalization techniques. IECs in the fetal, but not mature, intestine express chemerin. Chemerin expression peaked in the fetal intestine at 20-24 wk and then decreased to original low levels by full term. During the 10- to 24-wk period, chemerin accounted for most of the macrophage chemotactic activity of cultured fetal IECs. The maturational changes in chemerin expression correlated with the expression of retinoic acid receptor-beta in the intestine. Chemerin is an important mediator of epithelial-macrophage cross talk in the fetal/premature, but not in the mature, intestine. Understanding the regulation of the gut macrophage pool is an important step in development of novel strategies to boost mucosal immunity in premature infants and other patient populations at risk of microbial translocation.

TATA box polymorphisms in human gene promoters and associated hereditary pathologies

TATA-binding protein (TBP) is the first basal factor that recognizes and binds a TATA box on TATA-containing gene promoters transcribed by RNA polymerase II. Data available in the literature are indicative of admissible variability of the TATA box. The TATA box flanking sequences can influence TBP affinity as well as the level of basal and activated transcription. The possibility of mediated involvement in in vivo gene expression regulation of the TBP interactions with variant TATA boxes is supported by data on TATA box polymorphisms and associated human hereditary pathologies. A table containing data on TATA element polymorphisms in human gene promoters (about 40 mutations have been described), associated with particular pathologies, their short functional characteristics, and manifestation mechanisms of TATA-box SNPs is presented. Four classes of polymorphisms are considered: TATA box polymorphisms that weaken and enhance promoter, polymorphisms causing TATA box emergence and disappearance, and human virus TATA box polymorphisms. The described examples are indicative of the polymorphism-associated severe pathologies like thalassemia, the increased risk of hepatocellular carcinoma, sensitivity to H. pylori infection, oral cavity and lung cancers, arterial hypertension, etc.

Structure, expression, and biological function of INSM1 transcription factor in neuroendocrine differentiation

Zinc-finger transcription factors are DNA-binding proteins that are implicated in many diverse biological functions. INSM1 (formerly IA-1) contains five zinc-finger motifs and functions as a transcription factor. INSM1 protein structure is highly conserved in homologues of different species. It is predominantly expressed in developing neuroendocrine tissues and the nervous system in mammals. INSM1 represents an important player in early embryonic neurogenesis. In pancreatic endocrine cell differentiation, Ngn3 first activates INSM1 and subsequently NeuroD/beta2. Conversely, INSM1 exerts a feedback mechanism to suppress NeuroD/beta2 and its own gene expression. INSM1 gene ablation in the mouse results in the impairment of pancreatic endocrine cell maturation. Further, deletion of INSM1 severely impairs catecholamine biosynthesis and secretion from the adrenal gland that results in early embryonic lethality. Genetically, INSM1 acts as a downstream factor of Mash 1 and Phox2b in the differentiation of the sympatho-adrenal lineage. In the developing neocortex, mouse embryos lacking INSM1 expression contain half the number of basal progenitors and show a reduction in cortical plate radial thickness. Cell signaling studies reveal that INSM1 contributes to the induction of cell cycle arrest/exit necessary to facilitate cellular differentiation. INSM1 is highly expressed in tumors of neuroendocrine origin. Hence, its promoter could serve as a tumor-specific promoter that drives a specific targeted cancer gene therapy for the treatment of neuroendocrine tumors. Taken together, all of these features of INSM1 strongly support its role as an important regulator during neuroendocrine differentiation.

Interferon regulatory factors are transcriptional regulators of adipogenesis

We have sought to identify transcriptional pathways in adipogenesis using an integrated experimental and computational approach. Here, we employ high-throughput DNase hypersensitivity analysis to find regions of altered chromatin structure surrounding key adipocyte genes. Regions that display differentiation-dependent changes in hypersensitivity were used to predict binding sites for proteins involved in adipogenesis. A high-scoring example was a binding motif for interferon regulatory factor (IRF) family members. Expression of all nine mammalian IRF mRNAs is regulated during adipogenesis, and several bind to the identified motifs in a differentiation-dependent manner. Furthermore, several IRF proteins repress differentiation. This analysis suggests an important role for IRF proteins in adipocyte biology and demonstrates the utility of this approach in identifying cis- and trans-acting factors not previously suspected to participate in adipogenesis.