Identification of the DNA binding element of the human ZNF333 protein

Identifying genetic loci that are associated with changes in gene expression in PTSD in a South African cohort

The molecular mechanisms underlying posttraumatic stress disorder (PTSD) are yet to be fully elucidated, especially in underrepresented population groups. Expression quantitative trait loci (eQTLs) are DNA sequence variants that influence gene expression, in a local (cis-) or distal (trans-) manner, and subsequently impact cellular, tissue, and system physiology. This study aims to identify genetic loci associated with gene expression changes in a South African PTSD cohort. Genome-wide genotype and RNA-sequencing data were obtained from 32 trauma-exposed controls and 35 PTSD cases of mixed-ancestry, as part of the SHARED ROOTS project. The first approach utilised 108 937 single-nucleotide polymorphisms (SNPs) (MAF > 10%) and 11 312 genes with Matrix eQTL to map potential eQTLs, while controlling for covariates as appropriate. The second analysis was focused on 5638 SNPs related to a previously calculated PTSD polygenic risk score for this cohort. SNP-gene pairs were considered eQTLs if they surpassed Bonferroni correction and had a false discovery rate <0.05. We did not identify eQTLs that significantly influenced gene expression in a PTSD-dependent manner. However, several known cis-eQTLs, independent of PTSD diagnosis, were observed. rs8521 (C > T) was associated with TAGLN and SIDT2 expression, and rs11085906 (C > T) was associated with ZNF333 expression. This exploratory study provides insight into the molecular mechanisms associated with PTSD in a non-European, admixed sample population. This study was limited by the cross-sectional design and insufficient statistical power. Overall, this study should encourage further multi-omics approaches towards investigating PTSD in diverse populations.

Genetic variation influencing DNA methylation provides insights into molecular mechanisms regulating genomic function

We determined the relationships between DNA sequence variation and DNA methylation using blood samples from 3,799 Europeans and 3,195 South Asians. We identify 11,165,559 SNP-CpG associations (methylation quantitative trait loci (meQTL), P < 10-14), including 467,915 meQTL that operate in trans. The meQTL are enriched for functionally relevant characteristics, including shared chromatin state, High-throuhgput chromosome conformation interaction, and association with gene expression, metabolic variation and clinical traits. We use molecular interaction and colocalization analyses to identify multiple nuclear regulatory pathways linking meQTL loci to phenotypic variation, including UBASH3B (body mass index), NFKBIE (rheumatoid arthritis), MGA (blood pressure) and COMMD7 (white cell counts). For rs6511961 , chromatin immunoprecipitation followed by sequencing (ChIP-seq) validates zinc finger protein (ZNF)333 as the likely trans acting effector protein. Finally, we used interaction analyses to identify population- and lineage-specific meQTL, including rs174548 in FADS1, with the strongest effect in CD8+ T cells, thus linking fatty acid metabolism with immune dysregulation and asthma. Our study advances understanding of the potential pathways linking genetic variation to human phenotype.

Essential role of ALKBH5-mediated RNA demethylation modification in bile acid-induced gastric intestinal metaplasia

Bile acid reflux and subsequent caudal-related homeobox 2 (CDX2) activation contribute to gastric intestinal metaplasia (IM), a precursor of gastric cancer; however, the mechanism underlying this phenomenon is unclear. Here, we demonstrate that alkylation repair homolog protein 5 (ALKBH5), a major RNA N⁶-adenosine demethylase, is required for bile acid-induced gastric IM. Mechanistically, we revealed the N⁶-methyladenosine (m⁶A) modification profile in gastric IM for the first time and identified ZNF333 as a novel m⁶A target of ALKBH5. ALKBH5 was shown to demethylate ZNF333 mRNA, leading to enhanced ZNF333 expression by abolishing m⁶A-YTHDF2-dependent mRNA degradation. In addition, ALKBH5 activated CDX2 and downstream intestinal markers by targeting the ZNF333/CYLD axis and activating NF-κB signaling. Reciprocally, p65, the key transcription factor of the canonical NF-κB pathway, enhanced the transcription activity of ALKBH5 in the nucleus, thus forming a positive feedforward circuit. Furthermore, ALKBH5 levels were positively correlated with ZNF333 and CDX2 levels in IM tissues, indicating significant clinical relevance. Collectively, our findings suggest that an m⁶A modification-associated positive feedforward loop between ALKBH5 and NF-κB signaling is involved in generating the IM phenotype of gastric epithelial cells. Targeting the ALKBH5/ZNF333/CYLD/CDX2 axis may be a useful therapeutic strategy for gastric IM in patients with bile regurgitation.

Apak competes with p53 for direct binding to intron 1 of p53AIP1 to regulate apoptosis

The KRAB-type zinc-finger protein Apak was recently identified as a negative regulator of p53-mediated apoptosis. However, the mechanism of this selective regulation is not fully understood. Here, we show that Apak recognizes the TCTTN2−30TTGT consensus sequence through its zinc-fingers. This sequence is specifically found in intron 1 of the proapoptotic p53 target gene p53AIP1 and largely overlaps with the p53-binding sequence. Apak competes with p53 for binding to this site to inhibit p53AIP1 expression. Upon DNA damage, Apak dissociates from the DNA, which abolishes its inhibitory effect on p53-mediated apoptosis.

ZNF424, a novel human KRAB/C2H2 zinc finger protein, suppresses NFAT and p21 pathway

Zinc finger-containing transcription factors are the largest single family of transcriptional regulators in mammals, which play an essential role in cell differentiation, cell proliferation, apoptosis, and neoplastic transformation. Here we have cloned a novel KRAB-related zinc finger gene, ZNF424, encoding a protein of 555aa. ZNF424 gene consisted of 4 exons and 3 introns, and mapped to chromosome 19p13.3. ZNF424 gene was ubiquitously expressed in human embryo tissues by Northern blot analysis. ZNF424 is conserved across species in evolution. Using a GFP-labeled ZNF424 protein, we demonstrate that ZNF424 localizes mostly in the nucleus. Transcriptional activity assays shows ZNF424 suppresses transcriptional activity of L8G5-luciferase. Overexpression of ZNF424 in HEK- 293 cells inhibited the transcriptional activity of NFAT and p21, which may be silenced by siRNA. The results suggest that ZNF424 protein may act as a transcriptional repressor that suppresses NFAT and p21 pathway to mediate cellular functions.

DNA light-strand preferential recognition of human mitochondria transcription termination factor mTERF

Transcription termination of the human mitochondrial genome requires specific binding to termination factor mTERF. In this study, mTERF was produced in E. coli and purified by two-step chromatography. mTERF-binding DNA sequences were isolated from a pool of randomized sequences by the repeated selection of bound sequences by gel-mobility shift assay and polymerase chain reaction. Sequencing and comparison of the 23 isolated clones revealed a 16-bp consensus sequence of 5\'-GTGTGGCA GANCCNGG-3\' in the light-strand (underlined residues were absolutely conserved), which nicely matched the genomic 13-bp terminator sequence 5\'-TGGCAGAGCCC GG-3\'. Moreover, mTERF binding assays of heteroduplex and single-stranded DNAs showed mTERF recognized the light strand in preference to the heavy strand. The preferential binding of mTERF with the light-strand may explain its distinct orientation-dependent termination activity.

The metastasis efficiency modifier ribosomal RNA processing 1 homolog B (RRP1B) is a chromatin-associated factor

There is accumulating evidence for a role of germ line variation in breast cancer metastasis. We have recently identified a novel metastasis susceptibility gene, Rrp1b (ribosomal RNA processing 1 homolog B). Overexpression of Rrp1b in a mouse mammary tumor cell line induces a gene expression signature that predicts survival in breast cancer. Here we extend the analysis of RRP1B function by demonstrating that the Rrp1b activation gene expression signature accurately predicted the outcome in three of four publicly available breast carcinoma gene expression data sets. In addition, we provide insights into the mechanism of RRP1B. Tandem affinity purification demonstrated that RRP1B physically interacts with many nucleosome binding factors, including histone H1X, poly(ADP-ribose) polymerase 1, TRIM28 (tripartite motif-containing 28), and CSDA (cold shock domain protein A). Co-immunofluorescence and co-immunoprecipitation confirmed these interactions and also interactions with heterochromatin protein-1alpha and acetyl-histone H4 lysine 5. Finally, we investigated the effects of ectopic expression of an RRP1B allelic variant previously associated with improved survival in breast cancer. Gene expression analyses demonstrate that, compared with ectopic expression of wild type RRP1B in HeLa cells, the variant RRP1B differentially modulates various transcription factors controlled by TRIM28 and CSDA. These data suggest that RRP1B, a tumor progression and metastasis susceptibility candidate gene, is potentially a dynamic modulator of transcription and chromatin structure.

Brain proteins interacting with the tetramerization region of non-erythroid alpha spectrin

The N-terminal region of non-erythroid alpha spectrin (SpαII) is responsible for interacting with its binding partner, beta spectrin, to form functional spectrin tetramers. We used a yeast-two-hybrid system, with an N-terminal segment of alpha spectrin representing the functional tetramerization site, as a bait to screen human brain c-DNA library for proteins that interact with the alpha spectrin segment. In addition to several beta spectrin isoforms, we identified 14 proteins that interact with SpαII. Seven of the 14 were matched to 6 known proteins: Duo protein, Lysyl-tRNA synthetase, TBP associated factor 1, two isoforms (b and c) of a protein kinase A interacting protein and Zinc finger protein 333 (2 different segments). Four of the 6 proteins are located primarily in the nucleus, suggesting that spectrin plays important roles in nuclear functions. The remaining 7 proteins were unknown to the protein data base. Structural predictions show that many of the 14 proteins consist of a large portion of unstructured regions, suggesting that many of these proteins fold into a rather flexible conformation. It is interesting to note that all but 3 of the 14 proteins are predicted to consist of one to four coiled coils (amphiphilic helices). A mutation in SpαII, V22D, which interferes with the coiled coil bundling of SpαII with beta spectrin, also affects SpαII interaction with Duo protein, TBP associated factor 1 and Lysyl-tRNA synthetase, suggesting that they may compete with beta spectrin for interaction with SpαII. Future structural and functional studies of these proteins to provide interaction mechanisms will no doubt lead to a better understanding of brain physiology and pathophysiology.