Core Transcription Factors Research Articles

Distinct oncogenic signatures in malignant PEComa and leiomyosarcoma identified by integrative RNA-seq and H3K27ac ChIP-seq analysis.

11552 Background: Malignant perivascular epithelioid cell tumor (PEComa) and leiomyosarcoma (LMS) are two sarcomas with overlapping morphologic and immunophenotypic features which can make their diagnostic distinction challenging. We aimed to characterize the transcriptional and epigenetic landscape of PEComa and LMS to identify distinguishing features. Methods: We performed whole transcriptome RNA-sequencing on 19 PEComas and compared their gene expression profile to 259 sarcomas from The Cancer Genome Atlas (TCGA) including 104 LMS. ChIP-sequencing for H3K27ac, a histone modification associated with activation of nearby genes/open chromatin, was conducted on 9 malignant PEComas and 12 LMS and were compared with publicly available data from 4 other sarcoma subtypes (chordoma; osteosarcoma; undifferentiated pleomorphic sarcoma; rhabdomyosarcoma; n = 29 tumors). Results: Genome-wide epigenetic and transcriptional analyses revealed overlapping patterns between PEComa and LMS, which were distinct from other sarcomas. However, we also identified a set of highly expressed and epigenetically distinct transcripts which may represent diagnostic?biomarkers: e.g., DAPL1, MLANA, SULT1C2, GPR143, and CHI3L1 for PEComa; and MYOCD, WDFC2, DES, MYH11, and CNN1 for LMS; each of which showed >17x fold higher expression for each tumor entity by DESeq2 (FDR<0.0001). Gene Set Enrichment Analyses (GSEA) demonstrated enrichment in the KEGG Lysosome pathway for PEComa (FDR=0.11), whereas myogenesis and smooth muscle contraction pathways were enriched in LMS (FDR=0.09). Integrative transcriptomic and epigenetic analyses revealed a unique set of master core transcription factors for each tumor type including among others MYOCD for LMS; MITF for PEComa, which require further functional investigation. Twelve selected genes including new as well as known and standard diagnostic markers (e.g., DAPL1, MLANA, GPR143, PNL2, CHI3L1, DES, MYH11, ER, CD68, PU.1, pS6 and CNN1) were validated by immunohistochemistry (IHC) in multiple sections from PEComa and LMS (n = 26). The combination of three melanocytic markers (HMB45, MLANA, PNL2) and pS6 can distinguish LMS from PEComas (**** p<0.0001). IHC for CD68 and PU.1 macrophage markers did not show any difference regarding the degree of immune infiltration in PEComa vs. LMS. Conclusions: Our studies revealed novel epigenetic signatures translating into lysosomal and melanocytic proteins for PEComa and myogenic proteins for LMS, which may serve as useful diagnostic biomarkers in the distinction of these two sarcoma subtypes.

Stub1 maintains proteostasis of master transcription factors in embryonic stem cells.

The pluripotency and differentiation states of embryonic stem cells (ESCs) are regulated by a set of core transcription factors, primarily Sox2, Oct4, and Nanog. Although their transcriptional regulation has been studied extensively, the contribution of posttranslational modifications in Sox2, Oct4, and Nanog are poorly understood. Here, using a CRISPR-Cas9 knockout library screen in murine ESCs, we identify the E3 ubiquitin ligase Stub1 as a negative regulator of pluripotency. Manipulation of Stub1 expression in murine ESCs shows that ectopic Stub1 expression significantly reduces the protein half-life of Sox2, Oct4, and Nanog. Mechanistic investigations reveal Stub1 catalyzes the polyubiquitination and 26S proteasomal degradation of Sox2 and Nanog through K48-linked ubiquitin chains and Oct4 via K63 linkage. Stub1 deficiency positively enhances somatic cell reprogramming and delays differentiation, whereas its enforced expression triggers ESC differentiation. The discovery of Stub1 as an integral pluripotency regulator strengthens our understanding of ESC regulation beyond conventional transcriptional control mechanisms.

Biallelic ATOH1 Gene Variant in Siblings With Pontocerebellar Hypoplasia, Developmental Delay, and Hearing Loss.

Background and ObjectivesTo report on the novel association of biallelic variant in atonal basic helix-loop-helix transcription factor 1 (ATOH1) gene and pontocerebellar hypoplasia (PCH), severe global developmental delay, intellectual disability, and hearing loss in a family with 2 affected siblings.MethodsA detailed clinical assessment and exome sequencing of peripheral blood sample were performed. Segregation analysis with Sanger sequencing and structural modeling of the variant was performed to support the pathogenicity of the variant.ResultsA homozygous missense variant (NM_005172.1:c.481C>G) in the ATOH1 gene was identified in the proband and his affected sister. The segregation analysis subsequently confirmed its segregation with an apparently recessive PCH in this family. ATOH1 encodes for the atonal basic helix-loop-helix (bHLH) transcription factor 1, a core transcription factor in the developing cerebellum, brainstem, and dorsal spinal cord, and in the ear. The identified variant results in the p.(Arg161Gly) amino acid substitution in the evolutionarily conserved DNA-binding bHLH domain of the ATOH1 protein. Biallelic missense variants in this domain were previously reported to result in disordered cerebellar development and hearing loss in animal models. In silico homology modeling revealed that p.Arg161Gly in ATOH1 protein probably disrupts a salt bridge with DNA backbone phosphate and increases the flexibility of the bHLH helix–both of which together affect the binding capability of the bHLH domain to the DNA.DiscussionBased on the sequencing results and evidence from structural modeling of the identified variant, as well as with previous reports of ATOH1 gene disruption, we conclude that ATOH1 may represent a novel candidate gene associated with the phenotype of PCH, global developmental delay, and hearing loss in humans.

GATA2-Mediated Transcriptional Activation of Notch3 Promotes Pancreatic Cancer Liver Metastasis.

The liver is the predominant metastatic site for pancreatic cancer. However, the factors that determine the liver metastasis and the specific molecular mechanisms are still unclear. In this study, we used human pancreatic cancer cell line Hs766T to establish Hs766T-L3, a subline of Hs766T with stable liver metastatic ability. We performed RNA sequencing of Hs766T-L3 and its parental cell line Hs766T, and revealed huge differences in gene expression patterns and pathway activation between these two cell lines. We correlated the difference in pathway activation with the expression of the four core transcriptional factors including STAT1, NR2F2, GATA2, and SMAD4. Using the TCGA database, we examined the relative expression of these transcription factors (TFs) in pan-cancer and their relationship with the prognosis of the pancreatic cancer. Among these TFs, we considered GATA2 is closely involved in tumor metastasis and may serve as a potential metastatic driver. Further in vitro and in vivo experiments confirmed that GATA2-mediated transcriptional activation of Notch3 promotes the liver metastasis of Hs766T-L3, and knockdown of either GATA2 or Notch3 reduces the metastatic ability of Hs766T-L3. Therefore, we claim that GATA2 may serve as a metastatic driver of pancreatic cancer and a potential therapeutic target to treat liver metastasis of pancreatic cancer.

CBL-Interacting Protein Kinase OsCIPK18 Regulates the Response of Ammonium Toxicity in Rice Roots.

Ammonium () is one of the major nitrogen sources for plants. However, excessive ammonium can cause serious harm to the growth and development of plants, i.e., ammonium toxicity. The primary regulatory mechanisms behind ammonium toxicity are still poorly characterized. In this study, we showed that OsCIPK18, a CBL-interacting protein kinase, plays an important role in response to ammonium toxicity by comparative analysis of the physiological and whole transcriptome of the T-DNA insertion mutant (cipk18) and the wild-type (WT). Root biomass and length of cipk18 are less inhibited by excess compared with WT, indicating increased resistance to ammonium toxicity. Transcriptome analysis reveals that OsCIPK18 affects the uptake by regulating the expression of OsAMT1;2 and other transporters, but does not affect ammonium assimilation. Differentially expressed genes induced by excess in WT and cipk18 were associated with functions, such as ion transport, metabolism, cell wall formation, and phytohormones signaling, suggesting a fundamental role for OsCIPK18 in ammonium toxicity. We further identified a transcriptional regulatory network downstream of OsCIPK18 under stress that is centered on several core transcription factors. Moreover, OsCIPK18 might function as a transmitter in the auxin and abscisic acid (ABA) signaling pathways affected by excess ammonium. These data allowed us to define an OsCIPK18-regulated/dependent transcriptomic network for the response of ammonium toxicity and provide new insights into the mechanisms underlying ammonium toxicity.

Cellular heterogeneity and transcriptomic profiles during intrahepatic cholangiocarcinoma initiation and progression.

Intrahepatic cholangiocarcinoma (ICC) is not fully investigated, and how stromal cells contribute to ICC formation is poorly understood. We aimed to uncover ICC origin, cellular heterogeneity, and critical modulators during ICC initiation/progression, and to decipher how fibroblast and endothelial cells in the stromal compartment favor ICC progression. We performed single-cell RNA sequencing (scRNA-seq) using AKT/Notch intracellular domain-induced mouse ICC tissues at early, middle, and late stages. We analyzed the transcriptomic landscape, cellular classification and evolution, and intercellular communication during ICC initiation/progression. We confirmed the findings using quantitative real-time PCR, western blotting, immunohistochemistry or immunofluorescence, and gene knockout/knockdown analysis. We identified stress-responding and proliferating subpopulations in late-stage mouse ICC tissues and validated them using human scRNA-seq data sets. By integrating weighted correlation network analysis and protein-protein interaction through least absolute shrinkage and selection operator regression, we identified zinc finger, MIZ-type containing 1 (Zmiz1) and Y box protein 1 (Ybx1) as core transcription factors required by stress-responding and proliferating ICC cells, respectively. Knockout of either one led to the blockade of ICC initiation/progression. Using two other ICC mouse models (YAP/AKT, KRAS/p19) and human ICC scRNA-seq data sets, we confirmed the orchestrating roles of Zmiz1 and Ybx1 in ICC occurrence and development. In addition, hes family bHLH transcription factor 1, cofilin 1, and inhibitor of DNA binding 1 were identified as driver genes for ICC. Moreover, periportal liver sinusoidal endothelial cells could differentiate into tip endothelial cells to promote ICC development, and this was Dll4-Notch4-Efnb2 signaling-dependent. Stress-responding and ICC proliferating subtypes were identified, and Zmiz1 and Ybx1 were revealed as core transcription factors in these subtypes. Fibroblast-endothelial cell interaction promotes ICC development.

Hsa-miR-422a Originated from Short Interspersed Nuclear Element Increases ARID5B Expression by Collaborating with NF-E2.

MicroRNAs (miRNAs) are a class of small non-coding RNAs that regulate the expression of target messenger RNA (mRNA) complementary to the 3’ untranslated region (UTR) at the post-transcriptional level. Hsa-miR-422a, which is commonly known as miRNA derived from transposable element (MDTE), was derived from short interspersed nuclear element (SINE). Through expression analysis, hsa-miR-422a was found to be highly expressed in both the small intestine and liver of crab-eating monkey. AT-Rich Interaction Domain 5 B (ARID5B) was selected as the target gene of hsa-miR-422a, which has two binding sites in both the exon and 3’UTR of ARID5B. To identify the interaction between hsa-miR-422a and ARID5B, a dual luciferase assay was conducted in HepG2 cell line. The luciferase activity of cells treated with the hsa-miR-422a mimic was upregulated and inversely downregulated when both the hsa-miR-422a mimic and inhibitor were administered. Nuclear factor erythroid-2 (NF-E2) was selected as the core transcription factor (TF) via feed forward loop analysis. The luciferase expression was downregulated when both the hsa-miR-422a mimic and siRNA of NF-E2 were treated, compared to the treatment of the hsa-miR-422a mimic alone. The present study suggests that hsa-miR-422a derived from SINE could bind to the exon region as well as the 3’UTR of ARID5B. Additionally, hsa-miR-422a was found to share binding sites in ARID5B with several TFs, including NF-E2. The hsa-miR-422a might thus interact with TF to regulate the expression of ARID5B, as demonstrated experimentally. Altogether, hsa-miR-422a acts as a super enhancer miRNA of ARID5B by collaborating with TF and NF-E2.

The dynamic alteration of transcriptional regulation by crucial TFs during tumorigenesis of gastric cancer

BackgroundThe mechanisms of Gastric cancer (GC) initiation and progression are complicated, at least partly owing to the dynamic changes of gene regulation during carcinogenesis. Thus, investigations on the changes in regulatory networks can improve the understanding of cancer development and provide novel insights into the molecular mechanisms of cancer.MethodsDifferential co-expression analysis (DCEA), differential gene regulation network (GRN) modeling and differential regulation analysis (DRA) were integrated to detect differential transcriptional regulation events between gastric normal mucosa and cancer samples based on GSE54129 dataset. Cytological experiments and IHC staining assays were used to validate the dynamic changes of CREB1 regulated targets in different stages.ResultsA total of 1955 differentially regulated genes (DRGs) were identified and prioritized in a quantitative way. Among the top 1% DRGs, 14 out of 19 genes have been reported to be GC relevant. The four transcription factors (TFs) among the top 1% DRGs, including CREB1, BPTF, GATA6 and CEBPA, were regarded as crucial TFs relevant to GC progression. The differentially regulated links (DRLs) around the four crucial TFs were then prioritized to generate testable hypotheses on the differential regulation mechanisms of gastric carcinogenesis. To validate the dynamic alterations of gene regulation patterns of crucial TFs during GC progression, we took CREB1 as an example to screen its differentially regulated targets by using cytological and IHC staining assays. Eventually, TCEAL2 and MBNL1 were proved to be differentially regulated by CREB1 during tumorigenesis of gastric cancer.ConclusionsBy combining differential networking information and molecular cell experiments verification, testable hypotheses on the regulation mechanisms of GC around the core TFs and their top ranked DRLs were generated. Since TCEAL2 and MBNL1 have been reported to be potential therapeutic targets in SCLC and breast cancer respectively, their translation values in GC are worthy of further investigation.

Characterization of nanog in Nile tilapia (Oreochromis niloticus) and its spatiotemporal expression patterns during embryonic and gonadal development.

Nanog is one of the well-characterized core transcription factors in pluripotency maintenance network. So far, studies on fishes indicate that the Nanog expression occurs from embryonic 1-cell stage to blastula stage, and is restricted to the gonadal germline cells in adult tissues, which is strikingly different from that in mammals. However, whether this expression profile is conservative in fishes remains to be investigated. Here Nile tilapia (Oreochromis niloticus) nanog (named as Ong) was identified and its spatiotemporal expression patterns during embryonic and gonadal development were investigated. The Ong cDNA contains an open reading frame of 678bp, encoding 226 amino acids. The anti-Ong antibody was prepared through prokaryotic protein expression and its specificity was validated. The Ong expression in embryonic 1-cell stage did not appear until the early stage of blastocyst and continued to the late stage of blastocyst. In adult tissues, its expression was limited to gonads. Its expression patterns during gonadal development were further investigated by in situ hybridization and immunohistochemical staining. In testis, Ong was not expressed at 30 dah (days after hatching), but highly expressed in spermatogonia and spermatocytes at 150 dah; in ovaries at 30 and 150 dah, it was not expressed in germline cells but in all somatic cells. This expression profile is strikingly different from reports in fishes to date. Our study firstly indicates that the Nanog expression profile is not conservative in fishes. This study is valuable for further functional and evolutionary study of this gene.

DeepCAGE: Incorporating Transcription Factors in Genome-wide Prediction of Chromatin Accessibility

Although computational approaches have been complementing high-throughput biological experiments for the identification of functional regions in the human genome, it remains a great challenge to systematically decipher interactions between transcription factors (TFs) and regulatory elements to achieve interpretable annotations of chromatin accessibility across diverse cellular contexts. To solve this problem, we propose DeepCAGE, a deep learning framework that integrates sequence information and binding statuses of TFs, for the accurate prediction of chromatin accessible regions at a genome-wide scale in a variety of cell types. DeepCAGE takes advantage of a densely connected deep convolutional neural network architecture to automatically learn sequence signatures of known chromatin accessible regions and then incorporates such features with expression levels and binding activities of human core TFs to predict novel chromatin accessible regions. In a series of systematic comparisons with existing methods, DeepCAGE exhibits superior performance in not only the classification but also the regression of chromatin accessibility signals. In a detailed analysis of TF activities, DeepCAGE successfully extracts novel binding motifs and measures the contribution of a TF to the regulation with respect to a specific locus in a certain cell type. When applied to whole-genome sequencing data analysis, our method successfully prioritizes putative deleterious variants underlying a human complex trait and thus provides insights into the understanding of disease-associated genetic variants. DeepCAGE can be downloaded from https://github.com/kimmo1019/DeepCAGE.

Comprehensive 3D epigenomic maps define limbal stem/progenitor cell function and identity

The insights into how genome topology couples with epigenetic states to govern the function and identity of the corneal epithelium are poorly understood. Here, we generate a high-resolution Hi-C interaction map of human limbal stem/progenitor cells (LSCs) and show that chromatin multi-hierarchical organisation is coupled to gene expression. By integrating Hi-C, epigenome and transcriptome data, we characterize the comprehensive 3D epigenomic landscapes of LSCs. We find that super-silencers mediate gene repression associated with corneal development, differentiation and disease via chromatin looping and/or proximity. Super-enhancer (SE) interaction analysis identified a set of SE interactive hubs that contribute to LSC-specific gene activation. These active and inactive element-anchored loop networks occur within the cohesin-occupied CTCF-CTCF loops. We further reveal a coordinated regulatory network of core transcription factors based on SE-promoter interactions. Our results provide detailed insights into the genome organization principle for epigenetic regulation of gene expression in stratified epithelia.

E2F4 may be a core transcription factor in the lncRNA-TF regulatory network in cervical cancer.

BackgroundCervical cancer is the most common gynecological cancer worldwide and is associated with high morbidity and mortality. Despite improvements in therapeutic strategies, the network regulation mechanism remains unclear and the treatment effect is not satisfactory. Therefore, there is a need to continue studying the mechanism of cervical cancer to explore effective gene targets and precise targeted therapy drugs.MethodsFirst, three paired tissues (cancer tissues and noncancerous tissues) from patients with cervical squamous cell carcinoma were collected, grouped, and analyzed by microarray. Second, differentially expressed mRNAs (DEMs) and differentially expressed lncRNAs (DELs) (|fold change| ≥ 2 and p < 0.05) between the two groups were screened. For DEMs, functional annotation and pathway analysis were performed using DAVID. Functional prediction of DELs was then performed and their cis‐regulatory and trans‐regulatory networks were explored.ResultsFunction prediction of DELs (both up‐regulated and down‐regulated) shows that the highest frequency Cellular Component (CC) item is cytosol, the highest frequency Molecular function (MF) item is mitotic cell cycle and the highest frequency Biological Process (BP) item is protein binding. Through cis‐regulation analysis of DELs, the cis‐regulatory relationship of 96 DELs was predicted. The lncRNA‐trans‐regulation network analysis suggested that E2F4 may be the core transcription factor in the lncRNA‐TF regulatory network in cervical cancer.ConclusionsThe lncRNA‐TF regulatory network plays an important role in the occurrence and progression of cervical cancer, and E2F4 may be a critical transcription factor in the regulatory network.

BET proteins are essential for the specification and maintenance of the epiblast lineage in mouse preimplantation embryos

BackgroundDuring mammalian preimplantation development, as the fertilized egg develops and differentiates, three cell lineages become specified: trophectoderm (TE), epiblast, and primitive endoderm (PrE). Through two steps of cell fate decisions, 16-cell blastomeres develop into TE and an inner cell mass (ICM), and thereafter, the latter differentiates into pluripotent epiblast and PrE. Although bromodomain and extra-terminal domain (BET) proteins, such as BRD4, are necessary for the transcriptional activation of genes involved in the maintenance of mouse embryonic stem cells by occupying their enhancers, their roles in the development of mouse preimplantation are unknown.ResultsTo evaluate the effect of BET protein deficiency on cell lineage formation, we cultured preimplantation embryos in the presence of JQ1, which blocks the binding of BET bromodomains to acetylated-histones. We found BET inhibition blocked the transcriptional activation of genes, such as Nanog, Otx2, and Sox2, important for the formation of the epiblast lineage in blastocysts. Expression studies with lineage-specific markers in morulae and blastocysts revealed BET proteins were essential for the specification and maintenance of the epiblast lineage but were dispensable for the formation of primarily extraembryonic TE and PrE lineages. Additional Ingenuity Pathway Analysis and expression studies with a transcriptionally active form of signal transducer and activator of the transcription 3 (STAT3) suggested BET-dependent activation was partly associated with the STAT3-dependent pathway to maintain the epiblast lineage. To identify BET proteins involved in the formation of the epiblast lineage, we analyzed mutant embryos deficient in Brd4, Brd2, and double mutants. Abolishment of NANOG-positive epiblast cells was only evident in Brd4/Brd2 double-deficient morulae. Thus, the phenotype of JQ1-treated embryos is reproduced not by a Brd4- or Brd2-single deficiency, but only Brd4/Brd2-double deficiency, demonstrating the redundant roles of BRD2 and BRD4 in the specification of the epiblast lineage.ConclusionsBET proteins are essential to the specification and maintenance of the epiblast lineage by activating lineage-specific core transcription factors during mouse preimplantation development. Among BET proteins, BRD4 plays a central role and BRD2 a complementary role in the specification and maintenance of epiblast lineages. Additionally, BET-dependent maintenance of the epiblast lineage may be partly associated with the STAT3-dependent pathway.

UPF1 contributes to the maintenance of endometrial cancer stem cell phenotype by stabilizing LINC00963

Endometrial cancer stem cells (ECSCs) play a vital role in endometrial cancer (EC) metastasis, relapse, and chemoresistance. However, the molecular mechanisms that sustain ECSCs remain elusive. Here, we showed that the expression of UPF1 was upregulated in EC tissues and ECSCs and correlated with poor clinicopathological characteristics. UPF1 silencing suppressed ECSC hallmarks, such as sphere formation ability, carboplatin resistance, migration and invasion, and cell cycle progression. UPF1 regulated the behavior and fate of ECSCs by stabilizing LINC00963. LINC00963 further shares the same miRNA response element with the core transcription factor SOX2 and relieved the suppression of SOX2 by miR-508-5p in self-renewing ECSCs. Notably, inhibition of UPF1 and LINC00963 in combination severely impaired the in vivo tumorigenic potential of ECSCs. We demonstrate that the UPF1/LINC00963/miR-508-5p/SOX2 axis has potential value in modulating ECSC maintenance, chemoresistance, and tumorigenesis in EC, which highlights a novel promising target for EC treatment.

Etiological roles of core promoter variation in triple-negative breast cancer

Abnormal gene expression plays key role in cancer development. A core promoter is located around the transcriptional start site. Through interaction between core promoter sequences and transcriptional factors, core promoter controls transcriptional initiation. We hypothesized that in cancer, core promoter sequences could be mutated to interfere the interaction with transcriptional factors, resulting in altered transcriptional initiation and abnormal gene expression and cancer development. We used triple-negative breast cancer (TNBC) as a model to test our hypothesis. We collected genome-wide core promoter variants from 279 TNBC genomes. After extensive filtering of normal genomic polymorphism, we identified 19,427 recurrent somatic variants in 1,238 core promoters of 1,274 genes and 1,694 recurrent germline variants in 272 core promoters of 294 genes. Many of the affected genes were oncogenes and tumor suppressors. Analysis of RNA-seq data from the same patient cohort identified increased or decreased gene expression in 439 somatic and 85 germline variants-affected genes, and the results were validated by luciferase reporter assay. By comparing with the core promoter variation data from 610 unclassified breast cancer, we observed that core promoter variants in TNBC were highly TNBC-specific. We further identified the drugs targeting the genes with core promoter variation. Our study demonstrates that core promoter is highly mutable in cancer, and can play etiological roles in TNBC and other types of cancer through influencing transcriptional initiation.

Deep neural network prediction of genome-wide transcriptome signatures \u2013 beyond the Black-box

Prediction algorithms for protein or gene structures, including transcription factor binding from sequence information, have been transformative in understanding gene regulation. Here we ask whether human transcriptomic profiles can be predicted solely from the expression of transcription factors (TFs). We find that the expression of 1600 TFs can explain >95% of the variance in 25,000 genes. Using the light-up technique to inspect the trained NN, we find an over-representation of known TF-gene regulations. Furthermore, the learned prediction network has a hierarchical organization. A smaller set of around 125 core TFs could explain close to 80% of the variance. Interestingly, reducing the number of TFs below 500 induces a rapid decline in prediction performance. Next, we evaluated the prediction model using transcriptional data from 22 human diseases. The TFs were sufficient to predict the dysregulation of the target genes (rho = 0.61, P < 10−216). By inspecting the model, key causative TFs could be extracted for subsequent validation using disease-associated genetic variants. We demonstrate a methodology for constructing an interpretable neural network predictor, where analyses of the predictors identified key TFs that were inducing transcriptional changes during disease.

Genetic effect of ischemia-reperfusion injury upon primary graft dysfunction and chronic lung allograft dysfunction in lung transplantation: evidence based on transcriptome data

The unclear mechanism that ischemia–reperfusion injury (IRI) contributes to the development of primary graft dysfunction (PGD) and chronic lung allograft dysfunction (CLAD) remains a major issue in lung transplantation. Differentially expressed PGD-related genes and CLAD-related genes during IRI (IRI-PGD common genes and IRI-CLAD common genes) were identified using GEO datasets (GSE127003, GSE8021, GSE9102) and GeneCards datasets. Enrichment analysis and four network analyses, namely, protein-protein interaction, microRNA (miRNA)-gene, transcription factor (TF)-gene, and drug-gene networks, were then performed. Moreover, GSE161520 was analyzed to identify the differentially expressed core miRNAs during IRI in rats. Finally, Pearson correlation analysis and ROC analysis were performed. Eight IRI-PGD common genes (IL6, TNF, IL1A, IL1B, CSF3, CXCL8, SERPINE1, and PADI4) and 10 IRI-CLAD common genes (IL1A, ICAM1, CCL20, CCL2, IL1B, TNF, PADI4, CXCL8, GZMB, and IL6) were identified. Enrichment analysis showed that both IRI-PGD and IRI-CLAD common genes were significantly enriched in “AGE-RAGE signaling pathway in diabetic complication” and “IL-17 signaling pathway”. Among the core miRNAs, miR-1-3p and miR-335 were differentially expressed in IRI rats. Among core TFs, CEBPB expression had a significant negative correlation with P/F ratio (r = −0.33, P = 0.021). In the reperfused lung allografts, the strongest positive correlation was exhibited between PADI4 expression and neutrophil proportion (r = 0.76, P < 0.001), and the strongest negative correlation was between PADI4 expression and M2 macrophage proportion (r = −0.74, P < 0.001). In lung allografts of PGD recipients, IL6 expression correlated with activated dendritic cells proportion (r = 0.86, P < 0.01), and IL1B expression correlated with the neutrophils proportion(r = 0.84, P < 0.01). In whole blood of CLAD recipients, GZMB expression correlated with activated CD4+ memory T cells proportion (r = 0.76, P < 0.001).Our study provides the novel insights into the molecular mechanisms by which IRI contributes to PGD and CLAD and potential targets for therapeutic intervention.

Epigenetic Regulation of Megaspore Mother Cell Formation.

In flowering plants, the female gametophyte (FG) initiates from the formation of the megaspore mother cell (MMC). Among a pool of the somatic cells in the ovule primordium, only one hypodermal cell undergoes a transition of cell fate to become the MMC. Subsequently, the MMC undergoes a series of meiosis and mitosis to form the mature FG harboring seven cells with eight nuclei. Although SPL/NZZ, the core transcription factor for MMC formation, was identified several decades ago, which and why only one somatic cell is chosen as the MMC have long remained mysterious. A growing body of evidence reveal that MMC formation is associated with epigenetic regulation at multiple layers, including dynamic distribution of histone variants and histone modifications, small RNAs, and DNA methylation. In this review, we summarize the progress of epigenetic regulation in the MMC formation, emphasizing the roles of chromosome condensation, histone variants, histone methylation, small RNAs, and DNA methylation.

ZNF32 promotes the self-renewal of colorectal cancer cells by regulating the LEPR-STAT3 signaling pathway

Due to the self-renewal characteristics and tumorigenic abilities of cancer stem cells (CSCs), CSCs have been demonstrated to play vital roles in carcinogenesis and antitumor therapy. Our previous report found that Krüppel-like family members (KLFs) and zinc finger protein 32 (ZNF32) play oncogenic roles in carcinogenesis. However, the roles and mechanism of ZNF32 in CSCs are still unknown. Our study demonstrated that ZNF32 was highly expressed in colorectal CSCs, which promoted their self-renewal capacity and tumorigenicity. Overexpression of ZNF32 in colorectal cancer (CRC) cells increased their self-renewal capacity. Furthermore, we identified the leptin receptor (LEPR) as the downstream target gene of ZNF32 and verified that the ZNF32-mediated regulation of CRC self-renewal is achieved via the LEPR- signal transducer and activator of transcription 3 (STAT3) pathway. Moreover, ZNF32 regulated the expression of SOX2, a core transcription factor in stem cells. Finally, we demonstrated that ZNF32 and LEPR were positively correlated in CRC tissues. ZNF32 expression was negatively correlated with the prognosis of CRC patients. Therefore, therapeutically targeting the ZNF32-LEPR-STAT3 pathway in the clinic is tempting.

FOXL2 regulates the expression of the Col4a1 collagen gene in chicken granulosa cells.

Forkhead box L2 (FOXL2), one member in the superfamily of forkhead transcription factors, is a core transcription factor specifically expressed in ovarian granulosa cells and is essential for the development of follicles. FOXL2 has been shown to regulate the transcription of genes encoding enzymes that synthesize steroid hormones and estrogen receptors and regulate the expression of collagen genes in granulosa cells. This study explored the effect of FOXL2 on collagen gene expression in granulosa cells by overexpressing Foxl2 in pregranulosa cells, prehierarchical follicles and preovulation follicle granulosa cells. The results showed that FOXL2 regulated the expression of several genes encoding collagens in chicken granulosa cells and that overexpression of Foxl2 significantly reduced the messenger RNA and protein levels of Col4a1 in different granulosa cells. Moreover, luciferase reporter and chromatin immunoprecipitation assays were performed to study how FOXL2 regulates the expression of collagen genes, and the results showed that FOXL2 directly regulated the expression of Col4a1 by binding to the motif of CAGCAGCACCAGCAG between -640 and -625 bp upstream of the coding region. The results indicated that FOXL2 could regulate the components of the extracellular matrix; however, the biological significance of this regulation needs further clarification.