FORGEdb: a tool for identifying candidate functional variants and uncovering target genes and mechanisms for complex diseases

The majority of disease-associated variants identified through genome-wide association studies are located outside of protein-coding regions. Prioritizing candidate regulatory variants and gene targets to identify potential biological mechanisms for further functional experiments can be challenging. To address this challenge, we developed FORGEdb ( https://forgedb.cancer.gov/ ; https://forge2.altiusinstitute.org/files/forgedb.html ; and https://doi.org/10.5281/zenodo.10067458 ), a standalone and web-based tool that integrates multiple datasets, delivering information on associated regulatory elements, transcription factor binding sites, and target genes for over 37 million variants. FORGEdb scores provide researchers with a quantitative assessment of the relative importance of each variant for targeted functional experiments.

Bioinformatics Research and qRT-PCR Verify Hub Genes and a Transcription Factor-MicroRNA Feedback Network in Intervertebral Disc Degeneration

The present study explores the potentials of bioinformatics analysis to identify hub genes linked to intervertebral disc degeneration (IDD) and explored the potential molecular mechanism of transcription factor-microRNA regulatory network. Furthermore, the hub genes were identified through quantitative reverse transcriptase PCR (qRT-PCR). GEO database expression profile datasets for candidate genes (GSE124272) were downloaded. Genes that were differentially expressed (DEGs) were detected utilizing limma technique in the R programming language. Search Tool for the Retrieval of Interacting Genes/Proteins and NetworkAnalyst software identified hub genes. The Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis as well as Gene Ontology annotation of the DEGs were performed using Metascape. Using Bioinformatics data from the TRRUST, StarBase, and TransmiR databases, a TF-miRNA-hub genes network was constructed. qRT-PCR was utilized to confirm the result. As compared to healthy persons, 521 DEGs, comprising 203 down-regulated and 318 up-regulated genes, as well as 7 core genes, were found in people with IDD. Analysis revealed that all seven essential genes were under-expressed. qRT-PCR further confirmed the low expression of these seven important genes. Based on the TRRUST database, 16 TFs that could target five junction genes were then predicted. According to the StarBase database, four miRNAs were linked to crucial genes, while the TransmiR database predicted regulatory connections between four miRNAs and five TFs. The expression of the TP53-(hsa-miR-183-5p)-CCNB1 TF-miRNA-mRNA interaction network was discovered to be correlated with IDD. Throughout this investigation, a network of TF-miRNA-mRNA connections was built for investigation of the probable molecular mechanisms responsible for IDD. The identification of hub genes associated with IDD may reveal promising IDD treatment strategies.

Epigenetic regulation of human-specific gene expression in the prefrontal cortex

BACKGROUND

Changes in gene expression levels during brain development are thought to have played an important role in the evolution of human cognition. With the advent of high-throughput sequencing technologies, changes in brain developmental expression patterns, as well as human-specific brain gene expression, have been characterized. However, interpreting the origin of evolutionarily advanced cognition in human brains requires a deeper understanding of the regulation of gene expression, including the epigenomic context, along the primate genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) to measure the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 acetylation (H3K27ac), both of which are associated with transcriptional activation in the prefrontal cortex of humans, chimpanzees, and rhesus macaques.

RESULTS

We found a discrete functional association, in which ^H3K4me3HP gain was significantly associated with myelination assembly and signaling transmission, while ^H3K4me3HP loss played a vital role in synaptic activity. Moreover, ^H3K27acHP gain was enriched in interneuron and oligodendrocyte markers, and ^H3K27acHP loss was enriched in CA1 pyramidal neuron markers. Using strand-specific RNA sequencing (ssRNA-seq), we first demonstrated that approximately 7 and 2% of human-specific expressed genes were epigenetically marked by ^H3K4me3HP and ^H3K27acHP, respectively, providing robust support for causal involvement of histones in gene expression. We also revealed the co-activation role of epigenetic modification and transcription factors in human-specific transcriptome evolution. Mechanistically, histone-modifying enzymes at least partially contribute to an epigenetic disturbance among primates, especially for the H3K27ac epigenomic marker. In line with this, peaks enriched in the macaque lineage were found to be driven by upregulated acetyl enzymes.

CONCLUSIONS

Our results comprehensively elucidated a causal species-specific gene-histone-enzyme landscape in the prefrontal cortex and highlighted the regulatory interaction that drove transcriptional activation.

A translatome-transcriptome multi-omics gene regulatory network reveals the complicated functional landscape of maize

BACKGROUND

Maize (Zea mays L.) is one of the most important crops worldwide. Although sophisticated maize gene regulatory networks (GRNs) have been constructed for functional genomics and phenotypic dissection, a multi-omics GRN connecting the translatome and transcriptome is lacking, hampering our understanding and exploration of the maize regulatome.

RESULTS

We collect spatio-temporal translatome and transcriptome data and systematically explore the landscape of gene transcription and translation across 33 tissues or developmental stages of maize. Using this comprehensive transcriptome and translatome atlas, we construct a multi-omics GRN integrating mRNAs and translated mRNAs, demonstrating that translatome-related GRNs outperform GRNs solely using transcriptomic data and inter-omics GRNs outperform intra-omics GRNs in most cases. With the aid of the multi-omics GRN, we reconcile some known regulatory networks. We identify a novel transcription factor, ZmGRF6, which is associated with growth. Furthermore, we characterize a function related to drought response for the classic transcription factor ZmMYB31.

CONCLUSIONS

Our findings provide insights into spatio-temporal changes across maize development at both the transcriptome and translatome levels. Multi-omics GRNs represent a useful resource for dissection of the regulatory mechanisms underlying phenotypic variation.

Identification and Characterization of Transcription Factors Regulating Terpenoid Indole Alkaloid Biosynthesis in Catharanthus roseus

Biosynthesis of the therapeutically valuable terpenoid indole alkaloids (TIAs), in the medicinal plant Catharanthus roseus, is one of the most elaborate and complex metabolic processes. Although genomic and transcriptomic resources have significantly accelerated gene discovery in the TIA pathway, relatively few genes of transcription factors (TFs) have been identified and characterized thus far. Systematic identification of TFs and elucidation of their functions are crucial for understanding TIA pathway regulation. The successful discovery of TFs in the TIA pathway has relied mostly on three different approaches, (1) identification of cis-regulatory motifs (CRMs) present in the pathway gene promoters as they often provide clues on potential TFs that bind to the promoters, (2) co-expression analysis, based on the assumption that TFs regulating a metabolic or developmental pathway exhibit similar spatiotemporal expression as the pathway genes, and (3) isolation of homologs of TFs known to regulate structurally similar or diverse specialized metabolites in different plant species. TFs regulating TIA pathway have been isolated using either an individual or a combination of the three approaches. Here we describe transcriptome-based coexpression analysis and cis-element determination to identify TFs in C. roseus. In addition, we describe the protocols for generation of transgenic hairy roots, Agrobacterium infiltration of flowers, and electrophoretic mobility shift assay (EMSA). The methods described here are useful for the identification and characterization of potential TFs involved in the regulation of special metabolism in other medicinal plants.

DAP-Seq Identification of Transcription Factor-Binding Sites in Potato

Plant growth and adaptation to environmental fluctuations involve a tight control of cellular processes which, to a great extent, are mediated by changes at the transcriptional level. This regulation is exerted by transcription factors (TFs), a group of regulatory proteins that control gene expression by directly binding to the gene promoter regions via their cognate TF-binding sites (TFBS). The nature of TFBS defines the pattern of expression of the various plant loci, the precise combinatorial assembly of these elements being key in conferring plant's adaptation ability and in domestication. As such, TFs are main potential targets for biotechnological interventions, prompting in the last decade notable protein-DNA interaction efforts toward definition of their TFBS. Distinct methods based on in vivo or in vitro approaches defined the TFBS for many TFs, mainly in Arabidopsis, but comprehensive information on the transcriptional networks for many regulators is still lacking, especially in crops. In this chapter, detailed protocols for DAP-seq studies to unbiased identification of TFBS in potato are provided. This methodology relies on the affinity purification of genomic DNA-protein complexes in vitro, and high-throughput sequencing of the eluted DNA fragments. DAP-seq outperforms other in vitro DNA-motif definition strategies, such as protein-binding microarrays and SELEX-seq, since the protein of interest is directly bound to the genomic DNA extracted from plants yielding all the potential sites bound by the TF in the genome. Actually, data generated from DAP-seq experiments are highly similar to those out of ChIP-seq methods, but are generated much faster. We also provide a standard procedure to the analysis of the DAP-seq data, addressed to non-experienced users, that involves two consecutive steps: (1) processing of raw data (trimming, filtering, and read alignment) and (2) peak calling and identification of enriched motifs. This method allows identification of the binding profiles of dozens of TFs in crops, in a timely manner.

NFAT5 Regulated by STUB1, Facilitates Malignant Cell Survival and p38 MAPK Activation by Upregulating AQP5 in Chronic Lymphocytic Leukemia

Chronic lymphocytic leukemia (CLL) is a clonal proliferative disease of mature B lymphocytes. To further improve the prognosis of patients, it is necessary to further elucidate the pathogenesis of CLL and find more effective therapeutic targets. Nuclear Factor of Activated T cells 5 (NFAT5) is the major activated transcription factor (TF) upon osmotic pressure increase in mammalian cells, and it also regulates many target genes to affect various cellular functions. The effects of NFAT5 on tumor growth and metastasis have also been widely revealed. However, the effects of NFAT5 on the progression of CLL are still unclear. In this study, we found abnormally high expression of NFAT5 in human CLL patients. Additionally, NFAT5 depletion suppressed proliferation and stimulated apoptosis of CLL cells. Our data further confirmed NFAT5 regulated AQP5 expression and the phosphorylation of p38 MAPK. We also found that AQP5 overexpression reversed the inhibitory effect of NFAT5 depletion on cell proliferation in CLL cells. Furthermore, we revealed STUB1 directly bound to NFAT5 and promoted its degradation. Taken together, our results indicate the involvement of NFAT5 in CLL progression and suggest that NFAT5 could serve as a promising therapeutic target for CLL treatment.

Chromatin topology reorganization and transcription repression by PML-RARα in acute promyeloid leukemia

BACKGROUND

Acute promyeloid leukemia (APL) is characterized by the oncogenic fusion protein PML-RARα, a major etiological agent in APL. However, the molecular mechanisms underlying the role of PML-RARα in leukemogenesis remain largely unknown.

RESULTS

Using an inducible system, we comprehensively analyze the 3D genome organization in myeloid cells and its reorganization after PML-RARα induction and perform additional analyses in patient-derived APL cells with native PML-RARα. We discover that PML-RARα mediates extensive chromatin interactions genome-wide. Globally, it redefines the chromatin topology of the myeloid genome toward a more condensed configuration in APL cells; locally, it intrudes RNAPII-associated interaction domains, interrupts myeloid-specific transcription factors binding at enhancers and super-enhancers, and leads to transcriptional repression of genes critical for myeloid differentiation and maturation.

CONCLUSIONS

Our results not only provide novel topological insights for the roles of PML-RARα in transforming myeloid cells into leukemia cells, but further uncover a topological framework of a molecular mechanism for oncogenic fusion proteins in cancers.

Cell-Free Protein Synthesis of Plant Transcription Factors

Biotic and abiotic stimuli induce profound transcriptional reprograming in plants through sophisticated regulation of transcription factors (TFs). Recombinant proteins of TFs play an important role in unveiling their molecular functions. Cell-free protein synthesis (CFPS) system from wheat germ has been developed as one of the most efficient protein synthesis platforms. However, preparation of linear DNA templates for in vitro transcription is time-consuming and laborious. Here, we describe a versatile method for in vitro transcription and translation of the wheat germ CFPS system. Our two-step PCR method enables researchers to generate a variety of transcription templates from a single plasmid including fusion proteins of an N- or C-terminal tag and truncated proteins. Thus, this method supports a rapid and high-throughput expression of proteins for a large-scale proteomics analysis.

Multiple transcription factors contribute to inter-chromosomal interaction in yeast

BACKGROUND

Chromatin interactions medicated by genomic elements located throughout the genome play important roles in gene regulation and can be identified with the technologies such as high-throughput chromosome conformation capture (Hi-C), followed by next-generation sequencing. These techniques were wildly used to reveal the relative spatial disposition of chromatins in human, mouse and yeast. Unlike metazoan where CTCF plays major roles in mediating chromatin interactions, in yeast, the transcription factors (TFs) involved in this biological process are poorly known.

RESULTS

Here, we presented two computational approaches to estimate the TFs enriched in the chromatin physical inter-chromosomal interactions in yeast. Through the Chi-square method, we found TFs whose binding data are differentially distributed in different interaction groups, including Cin5, Stp1 and Sut1, whose binding data are negatively correlated with the chromosome spatial distance. A multivariate linear regression model was employed to estimate the potential contribution of different transcription factors against the physical distance of chromosomes. Rlr1, Set12 and Dig1 were found to be top positively participated in these chromosomal interactions. Ste12 was highlighted to be involved in gene reposition. Overall, we found 10 TFs enriched from both computational approaches, potentially to be involved in inter-chromosomal interactions.

CONCLUSIONS

No transcription factor (TF) in our study was found to have a dominant impact on the inter-chromosomal interaction as CTCF did in human or other metazoan, suggesting species without CTCF might have different regulatory systems in mediating inter-chromosomal interactions. In summary, we presented a systematic examination of TFs involved in chromatin interaction in yeast and the results provide candidate TFs for future studies.

OsSND2, a NAC family transcription factor, is involved in secondary cell wall biosynthesis through regulating MYBs expression in rice

BACKGROUND

As one of the most important staple food crops, rice produces huge agronomic biomass residues that contain lots of secondary cell walls (SCWs) comprising cellulose, hemicelluloses and lignin. The transcriptional regulation mechanism underlying SCWs biosynthesis remains elusive.

RESULTS

In this study, we isolated a NAC family transcription factor (TF), OsSND2 through yeast one-hybrid screening using the secondary wall NAC-binding element (SNBE) on the promoter region of OsMYB61 which is known transcription factor for regulation of SCWs biosynthesis as bait. We used an electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation analysis (ChIP) to further confirm that OsSND2 can directly bind to the promoter of OsMYB61 both in vitro and in vivo. OsSND2, a close homolog of AtSND2, is localized in the nucleus and has transcriptional activation activity. Expression pattern analysis indicated that OsSND2 was mainly expressed in internodes and panicles. Overexpression of OsSND2 resulted in rolled leaf, increased cellulose content and up-regulated expression of SCWs related genes. The knockout of OsSND2 using CRISPR/Cas9 system decreased cellulose content and down-regulated the expression of SCWs related genes. Furthermore, OsSND2 can also directly bind to the promoters of other MYB family TFs by transactivation analysis in yeast cells and rice protoplasts. Altogether, our findings suggest that OsSND2 may function as a master regulator to mediate SCWs biosynthesis.

CONCLUSION

OsSND2 was identified as a positive regulator of cellulose biosynthesis in rice. An increase in the expression level of this gene can improve the SCWs cellulose content. Therefore, the study of the function of OsSND2 can provide a strategy for manipulating plant biomass production.

MicroRNA alterations in iPSC-derived dopaminergic neurons from Parkinson disease patients

MicroRNA (miRNA) misregulation in peripheral blood has been linked to Parkinson disease (PD) but its role in the disease progression remains elusive. We performed an explorative genome-wide study of miRNA expression levels in dopaminergic neurons (DAn) from PD patients generated by somatic cell reprogramming and induced pluripotent stem cells differentiation. We quantified expression levels of 377 miRNAs in DAn from 3 sporadic PD patients (sPD), 3 leucine-rich repeat kinase 2-associated PD patients (L2PD) (total 6 PD), and 4 healthy controls. We identified differential expression of 10 miRNA of which 5 were upregulated in PD (miR-9-5p, miR-135a-5p, miR-135b-5p, miR-449a, and miR-449b-5p) and 5 downregulated (miR-141-3p, miR-199a-5p, miR-299-5p, miR-518e-3p, and miR-519a-3p). Changes were similar in sPD and L2PD. Integrative analysis revealed significant correlations between miRNA/mRNA expression. Moreover, upregulation of miR-9-5p and miR-135b-5p was associated with downregulation of transcription factors related to the DNA hypermethylation of enhancer elements in PD DAn (FOXA1 and NR3C1). In summary, miRNA changes are associated with monogenic L2PD and sPD and co-occur with epigenetic changes in DAn from PD patients.

Chromatin Immunoprecipitation Sequencing (ChIP-Seq) for Transcription Factors and Chromatin Factors in Arabidopsis thaliana Roots: From Material Collection to Data Analysis

Chromatin immunoprecipitation combined with next-generation sequencing (ChIP-seq) is a powerful technique to investigate in vivo transcription factor (TF) binding to DNA, as well as chromatin marks. Here we provide a detailed protocol for all the key steps to perform ChIP-seq in Arabidopsis thaliana roots, also working on other A. thaliana tissues and in most non-ligneous plants. We detail all steps from material collection, fixation, chromatin preparation, immunoprecipitation, library preparation, and finally computational analysis based on a combination of publicly available tools.

ABP9, a maize bZIP transcription factor, enhances tolerance to salt and drought in transgenic cotton

ABP9 , encoding a bZIP transcription factor from maize, enhances tolerance to multiple stresses and may participate in the ABA signaling pathway in transgenic cotton by altering physiological and biochemical processes and stress-related gene expression. Abiotic stresses, such as soil salinity and drought, negatively affect growth, development, and yield in cotton. Gene ABP9, which encodes a bZIP transcription factor, binds to the abscisic acid (ABA)-responsive-element (ABRE2) motif of the maize catalase1 gene. Its expression significantly improves tolerance in Arabidopsis to multiple abiotic stresses, but little is known about its role in cotton. In the present study, the ABP9 gene was introduced into upland cotton (Gossypium hirsutum L.) cultivar R15 by Agrobacterium tumefaciens-mediated transformation, and 12 independent transgenic cotton lines were obtained. Cotton plants over-expressing ABP9 have enhanced tolerance to salt and osmotic stress. Under stress, they developed better root systems in a greenhouse and higher germination, reduced stomatal aperture, and stomatal density in a growth chamber. Under drought conditions, survival rate and relative water content (RWC) of transgenic cotton were higher than those of R15 plants. Under salt and osmotic stresses, chlorophyll, proline, and soluble sugar contents significantly increased in transgenic cotton leaves and the malondialdehyde (MDA) content was lower than in R15. Overexpression of ABP9 also enhanced oxidative stress tolerance, reduced cellular levels of reactive oxygen species (ROS) through increased activities of antioxidative enzymes, and alleviated oxidative damage to cell. Interestingly, ABP9 over-expressing cotton was more sensitive to exogenous ABA than R15 at seed germination, root growth, stomatal aperture, and stomatal density. Moreover, ABP9 overexpression upregulated significantly the transcription levels of stress-related genes such as GhDBP2, GhNCED2, GhZFP1, GhERF1, GhHB1, and GhSAP1 under salt treatment. Conjointly, these results showed that overexpression of ABP9 conferred enhanced tolerance to multiple abiotic stresses in cotton. The stress-tolerant transgenic lines provide valuable resources for cotton breeding.

RSAT::Plants: Motif Discovery in ChIP-Seq Peaks of Plant Genomes

In this protocol, we explain how to run ab initio motif discovery in order to gather putative transcription factor binding motifs (TFBMs) from sets of genomic regions returned by ChIP-seq experiments. The protocol starts from a set of peak coordinates (genomic regions) which can be either downloaded from ChIP-seq databases, or produced by a peak-calling software tool. We provide a concise description of the successive steps to discover motifs, cluster the motifs returned by different motif discovery algorithms, and compare them with reference motif databases. The protocol is documented with detailed notes explaining the rationale underlying the choice of options. The interpretation of the results is illustrated with an example from the model plant Arabidopsis thaliana.

An integrated RNA-Seq and network study reveals a complex regulation process of rice embryo during seed germination

Seed germination is a crucial stage for plant development and agricultural production. To investigate its complex regulation process, the RNA-Seq study of rice embryo was conducted at three time points of 0, 12 and 48 h post imbibition (HPI). Dynamic transcriptional alterations were observed, especially in the early stage (0-12 HPI). Seed related genes, especially those encoding desiccation inducible proteins and storage reserves in embryo, decreased drastically after imbibition. The expression profiles of phytohormone related genes indicated distinct roles of abscisic acid (ABA), gibberellin (GA) and brassinosteroid (BR) in germination. Moreover, network analysis revealed the importance of protein phosphorylation in phytohormone interactions. Network and gene ontology (GO) analyses suggested that transcription factors (TFs) played a regulatory role in functional transitions during germination, and the enriched TF families at 0 HPI implied a regulation of epigenetic modification in dry seeds. In addition, 35 germination-specific TF genes in embryo were identified and seven genes were verified by qRT-PCR. Besides, enriched TF binding sites (TFBSs) supported physiological changes in germination. Overall, this study expands our comprehensive knowledge of multiple regulation factors underlying rice seed germination.

In silico-prediction of downstream MYB interacting partners of MAPK3 in Arabidopsis

Protein-Protein interactions (PPIs) are vital to most biological processes thus the identification of PPIs is of primary importance. Here, we endeavor to identify the downstream interacting partners of (AtMAPK3P) in Arabidopsis thaliana using the docking approach. Out of 131 members of MYB transcription factors 41 members are showing interactions with AtMAPK3P while the rest are showing non-interaction. Using minimal sequence motif search as well as through docking approach several novel MYB interacting proteins were also reported in the present study which need to be confirmed by in vitro kinase assay. Together, the results obtained essentially enhance our knowledge of the MAPK interacting protein network and provide a valuable research resource for developing a nearly important link between pathogen-activated MAPK signaling pathways and downstream transcriptional programming.

Chromatin immunoprecipitation and multiplex sequencing (ChIP-Seq) to identify global transcription factor binding sites in the nematode Caenorhabditis elegans

The global identification of transcription factor (TF) binding sites is a critical step in the elucidation of the functional elements of the genome. Several methods have been developed that map TF binding in human cells, yeast, and other model organisms. These methods make use of chromatin immunoprecipitation, or ChIP, and take advantage of the fact that formaldehyde fixation of living cells can be used to cross-link DNA sequences to the TFs that bind them in vivo. In ChIP, the cross-linked TF-DNA complexes are sheared by sonication, size fractionated, and incubated with antibody specific to the TF of interest to generate a library of TF-bound DNA sequences. ChIP-chip was the first technology developed to globally identify TF-bound DNA sequences and involves subsequent hybridization of the ChIP DNA to oligonucleotide microarrays. However, ChIP-chip proved to be costly, labor-intensive, and limited by the fixed number of probes available on the microarray chip. ChIP-Seq combines ChIP with massively parallel high-throughput sequencing (see Explanatory Chapter: Next Generation Sequencing) and has demonstrated vast improvement over ChIP-chip with respect to time and cost, signal-to-noise ratio, and resolution. In particular, multiplex sequencing can be used to achieve a higher throughput in ChIP-Seq analyses involving organisms with genomes of lower complexity than that of human (Lefrançois et al., 2009) and thereby reduce the cost and amount of time needed for each result. The multiplex ChIP-Seq method described in this section has been developed for Caenorhabditis elegans, but is easily adaptable for other organisms.

in silico-prediction of downstream WRKY interacting partners of MAPK3 in Brassica

Protein-Protein interactions (PPIs) are vital to most biological processes thus the identification of PPIs is of primary importance. Here, we endeavor to identify the downstream interacting partners of (BjMPK3P) in Brassica juncea using the docking approach. Out of 63 and 37 members of BrWRKY and BnWRKY transcription factors, 50 and 29 members are showing interactions with BjMPK3P respectively while the rest are showing non-interaction. Twenty two WRKY members are common to both the species. Using minimal sequence motif search as well as through docking approach several novel WRKY interacting proteins were also reported in the present study which need to be confirmed by in vitro kinase assay. Together, the results obtained essentially enhance our knowledge of the MAPK interacting protein network and provide a valuable research resource for developing a nearly important link between pathogen-activated MAPK signaling pathways and downstream transcriptional programming.