Search

Abstract Background Estriol (E3) is a steroid hormone formed only during pregnancy in primates including humans. Although E3 is synthesized at large amounts through a complex pathway involving the fetus and placenta, it is not required for the maintenance of pregnancy and has classically been considered virtually inactive due to associated very weak canonical estrogen signaling. However, estrogen exposure during pregnancy may have an effect on organs both within and outside the reproductive system, and compounds with binding affinity for estrogen receptors weaker than E3 have been found to impact reproductive organs and the brain. Here, we explore potential effects of E3 on fetal development using mouse as a model system. Results We administered E3 to pregnant mice, exposing the fetus to E3. Adult females exposed to E3 in utero (E3-mice) had increased fertility and superior pregnancy outcomes. Female and male E3-mice showed decreased anxiety and increased exploratory behavior. The expression levels and DNA methylation patterns of multiple genes in the uteri and brains of E3-mice were distinct from controls. E3 promoted complexing of estrogen receptors with several DNA/histone modifiers and their binding to target genes. E3 functions by driving epigenetic change, mediated through epigenetic modifier interactions with estrogen receptors rather than through canonical nuclear transcriptional activation. Conclusions We identify an unexpected functional role for E3 in fetal reproductive system and brain. We further identify a novel mechanism of estrogen action, through recruitment of epigenetic modifiers to estrogen receptors and their target genes, which is not correlated with the traditional view of estrogen potency.

Abstract Chromatin accessibility is essential for transcriptional activation of genomic regions. It is well established that transcription factors (TFs) and histone modifications (HMs) play critical roles in chromatin accessibility regulation. However, there is a lack of studies that quantify these relationships. Here we constructed a two-layer model to predict chromatin accessibility by integrating DNA sequence, TF binding, and HM signals. By applying the model to two human cell lines (GM12878 and HepG2), we found that DNA sequences had limited power for accessibility prediction, while both TF binding and HM signals predicted chromatin accessibility with high accuracy. According to the HM model, HM features determined chromatin accessibility in a cell line shared manner, with the prediction power attributing to five core HM types. Results from the TF model indicated that chromatin accessibility was determined by a subset of informative TFs including both cell line-specific and generic TFs. The combined model of both TF and HM signals did not further improve the prediction accuracy, indicating that they provide redundant information in terms of chromatin accessibility prediction. The TFs and HM models can also distinguish the chromatin accessibility of proximal versus distal transcription start sites with high accuracy.

Abstract Background Neuroblastoma (NB) is the most common extracranial solid tumor found in children. The frequent gain/loss of many chromosome bands in tumor cells and absence of mutations found at diagnosis suggests that NB is a copy number-driven cancer. Despite the previous work, a systematic analysis that investigates the relationship between such frequent gain/loss of chromosome bands and patient prognosis has yet to be implemented. Methods First, we analyzed two NB CNV datasets to select chromosomal bands with a high frequency of gain or loss. Second, we applied a computational approach to infer sample-specific CNVs for each chromosomal band selected in step 1 based on gene expression data. Third, we applied univariate Cox proportional hazards models to examine the association between the resulting inferred copy number values (iCNVs) and patient survival. Finally, we applied multivariate Cox proportional hazards models to select chromosomal bands that remained significantly associated with prognosis after adjusting for critical clinical variables, including age, stage, gender, and MYCN amplification status. Results Here, we used a computational method to infer the copy number variations (CNVs) of sample-specific chromosome bands from NB patient gene expression profiles. The resulting inferred CNVs (iCNVs) were highly correlated with the experimentally determined CNVs, demonstrating CNVs can be accurately inferred from gene expression profiles. Using this iCNV metric, we identified 58 frequent gain/loss chromosome bands that were significantly associated with patient survival. Furthermore, we found that 7 chromosome bands were still significantly associated with patient survival even when clinical factors, such as MYCN status, were considered. Particularly, we found that the chromosome band chr11p14 has high potential as a novel candidate cytogenetic biomarker for clinical use. Conclusion Our analysis resulted in a comprehensive list of prognostic chromosome bands supported by strong statistical evidence. In particular, the chr11p14 gain event provided additional prognostic value in addition to well-established clinical factors, including MYCN status, and thereby represents a novel candidate cytogenetic biomarker with high clinical potential. Additionally, this computational framework could be readily extended to other cancer types, such as leukemia.

Abstract Background Intra-tumor heterogeneity stems from genetic, epigenetic, functional, and environmental differences among tumor cells. A major source of genetic heterogeneity comes from DNA sequence differences and/or whole chromosome and focal copy number variations (CNVs). Whole chromosome CNVs are caused by chromosomal instability (CIN) that is defined by a persistently high rate of chromosome mis-segregation. Accordingly, CIN causes constantly changing karyotypes that result in extensive cell-to-cell genetic heterogeneity. How the genetic heterogeneity caused by CIN influences gene expression in individual cells remains unknown. Methods We performed single-cell RNA sequencing on a chromosomally unstable glioblastoma cancer stem cell (CSC) line and a control normal, diploid neural stem cell (NSC) line to investigate the impact of CNV due to CIN on gene expression. From the gene expression data, we computationally inferred large-scale CNVs in single cells. Also, we performed copy number adjusted differential gene expression analysis between NSCs and glioblastoma CSCs to identify copy number dependent and independent differentially expressed genes. Results Here, we demonstrate that gene expression across large genomic regions scales proportionally to whole chromosome copy number in chromosomally unstable CSCs. Also, we show that the differential expression of most genes between normal NSCs and glioblastoma CSCs is largely accounted for by copy number alterations. However, we identify 269 genes whose differential expression in glioblastoma CSCs relative to normal NSCs is independent of copy number. Moreover, a gene signature derived from the subset of genes that are differential expressed independent of copy number in glioblastoma CSCs correlates with tumor grade and is prognostic for patient survival. Conclusions These results demonstrate that CIN is directly responsible for gene expression changes and contributes to both genetic and transcriptional heterogeneity among glioblastoma CSCs. These results also demonstrate that the expression of some genes is buffered against changes in copy number, thus preserving some consistency in gene expression levels from cell-to-cell despite the continuous change in karyotype driven by CIN. Importantly, a gene signature derived from the subset of genes whose expression is buffered against copy number alterations correlates with tumor grade and is prognostic for patient survival that could facilitate patient diagnosis and treatment.