A comprehensive overview of genomic imprinting in breast and its deregulation in cancer

A compendium of methods: Searching allele specific expression via RNA sequencing

Diploid mammalian genome has paired alleles for each gene; typically allowing for equal expression of the two alleles within the cell/tissue. However, genetic regulatory elements and epigenetic modifications can disrupt this equality, leading to preferential expression of one allele. Examining high-confidence allele-specific expression (ASE) is vital for understanding genetic variations and their impact on major diseases like cancers and diabetes. ASE analysis not only aids in disease prognosis and diagnosis but also helps to identify regulatory mechanisms operating within the genome. While advances in sequencing technologies have greatly improved our understanding of ASE, challenges remain in estimating it accurately. In this article, we reviewed methods for detecting ASE using both bulk RNASeq and single-cell RNASeq data to provide deeper insights beyond the mere prediction of ASE genes. Fundamentally, ASE detection methods are data-driven and can be classified according to type of data used. Some methods utilize both, DNA genotyping information and RNASeq while others rely solely on RNASeq data. This article offers a comparative analysis of these methods and compilation of repositories providing valuable insights.

Characterization of Loss-of-Imprinting in Breast Cancer at the Cellular Level by Integrating Single-Cell Full-Length Transcriptome with Bulk RNA-Seq Data

Genomic imprinting, the parent-of-origin-specific gene expression, plays a pivotal role in growth regulation and is often dysregulated in cancer. However, screening for imprinting is complicated by its cell-type specificity, which bulk RNA-seq cannot capture. On the other hand, large-scale single-cell RNA-seq (scRNA-seq) often lacks transcript-level detail and is cost-prohibitive. Here, we address this gap by integrating bulk RNA-seq with full-length transcript scRNA-seq to investigate imprinting dynamics in breast cancer. By analyzing scRNA-seq data from 486 cancer cells across subtypes, we identified multiple SNPs in imprinted genes, including HM13, MEST (PEG1), SNHG14 and PEG10, showing consistent biallelic expression. Bulk RNA-seq, however, revealed that this biallelic expression arises from transcript-specific imprinting, rather than loss-of-imprinting (LOI). The imprinted SNPs identified in bulk RNA-seq predominantly demonstrate proper monoallelic expression in scRNA-seq. As a clear exception, an HER2+ breast cancer sample exhibited distinct LOI of MEST. Previous bulk RNA-seq-based observations about MEST LOI in breast cancer could not exclude a non-cancer cell impact, but our results validate that MEST LOI is cancer-specific. This study demonstrates the complementary utility of bulk and scRNA-seq in imprinting studies, confirming MEST LOI as a genuine event in breast cancer.

High expression of HM13 correlates with poor prognosis in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) has a high mortality rate, and the identification of early prognostic markers is crucial for improving patient outcomes. This study aimed to investigate the correlation between the expression of Histocompatibility Minor 13 (HM13) and the prognosis of HCC patients. HM13 protein expression was assessed in HCC tissues and cells through immunohistochemistry (IHC), quantitative reverse transcription PCR (qRT-PCR), and western blot. The relationship between HM13 expression and clinicopathological data of HCC was evaluated. Bioinformatics analyses, including Gene Expression Omnibus (GEO) database, Gene Expression Profiling Interactive Analysis (GEPIA), and Kaplan-Meier plotter (K-M plotter), were employed to analyze HM13 expression and its association with patient survival. HM13 was significantly overexpressed in HCC tissues and cells compared to normal controls. IHC revealed that HM13 protein was primarily localized in the cytoplasm and highly expressed in HCC tissues. Interestingly, patients with high HM13 expression had significantly poorer overall survival (OS), progression-free survival (PFS), recurrence-free survival (RFS), and disease-specific survival (DSS) than those with low expression. HM13 expression was associated with Edmondson grade, metastasis, microvascular invasion, and alpha-fetoprotein (AFP) levels. Multivariate analysis identified HM13 as an independent prognostic factor for poor OS in HCC. HM13 was markedly overexpressed in HCC and correlated with poor prognosis, suggesting its potential as a promising biomarker for early prognostic detection in HCC patients.

Identification of candidate causal variants and target genes at 41 breast cancer risk loci through differential allelic expression analysis

Understanding breast cancer genetic risk relies on identifying causal variants and candidate target genes in risk loci identified by genome-wide association studies (GWAS), which remains challenging. Since most loci fall in active gene regulatory regions, we developed a novel approach facilitated by pinpointing the variants with greater regulatory potential in the disease's tissue of origin. Through genome-wide differential allelic expression (DAE) analysis, using microarray data from 64 normal breast tissue samples, we mapped the variants associated with DAE (daeQTLs). Then, we intersected these with GWAS data to reveal candidate risk regulatory variants and analysed their cis-acting regulatory potential. Finally, we validated our approach by extensive functional analysis of the 5q14.1 breast cancer risk locus. We observed widespread gene expression regulation by cis-acting variants in breast tissue, with 65% of coding and noncoding expressed genes displaying DAE (daeGenes). We identified over 54 K daeQTLs for 6761 (26%) daeGenes, including 385 daeGenes harbouring variants previously associated with BC risk. We found 1431 daeQTLs mapped to 93 different loci in strong linkage disequilibrium with risk-associated variants (risk-daeQTLs), suggesting a link between risk-causing variants and cis-regulation. There were 122 risk-daeQTL with stronger cis-acting potential in active regulatory regions with protein binding evidence. These variants mapped to 41 risk loci, of which 29 had no previous report of target genes and were candidates for regulating the expression levels of 65 genes. As validation, we identified and functionally characterised five candidate causal variants at the 5q14.1 risk locus targeting the ATG10 and ATP6AP1L genes, likely acting via modulation of alternative transcription and transcription factor binding. Our study demonstrates the power of DAE analysis and daeQTL mapping to identify causal regulatory variants and target genes at breast cancer risk loci, including those with complex regulatory landscapes. It additionally provides a genome-wide resource of variants associated with DAE for future functional studies.

Imprinting as Basis for Complex Evolutionary Novelties in Eutherians

The epigenetic phenomenon of genomic imprinting is puzzling. While epigenetic modifications in general are widely known in most species, genomic imprinting in the animal kingdom is restricted to autosomes of therian mammals, mainly eutherians, and to a lesser extent in marsupials. Imprinting causes monoallelic gene expression. It represents functional haploidy of certain alleles while bearing the evolutionary cost of diploidization, which is the need of a complex cellular architecture and the danger of producing aneuploid cells by mitotic and meiotic errors. The parent-of-origin gene expression has stressed many theories. Most prominent theories, such as the kinship (parental conflict) hypothesis for maternally versus paternally derived alleles, explain only partial aspects of imprinting. The implementation of single-cell transcriptome analyses and epigenetic research allowed detailed study of monoallelic expression in a spatial and temporal manner and demonstrated a broader but much more complex and differentiated picture of imprinting. In this review, we summarize all these aspects but argue that imprinting is a functional haploidy that not only allows a better gene dosage control of critical genes but also increased cellular diversity and plasticity. Furthermore, we propose that only the occurrence of allele-specific gene regulation mechanisms allows the appearance of evolutionary novelties such as the placenta and the evolutionary expansion of the eutherian brain.

Loss-of-Imprinting of HM13 Leads to Poor Prognosis in Clear Cell Renal Cell Carcinoma

Genomic imprinting refers to the epigenetic silencing of one of both alleles in a parent-of-origin-specific manner, particularly in genes regulating growth and development. Impaired genomic imprinting leading to the activation of the silenced allele, also called canonical loss-of-imprinting (LOI), is considered an early factor in oncogenesis. As LOI studies in clear cell renal cell carcinoma (ccRCC) are limited to IGF2, we performed a genome-wide analysis in 128 kidney normal solid tissue and 240 stage 1 ccRCC samples (TCGA RNA-seq data) to screen for canonical LOI in early oncogenesis. In ccRCC, we observed LOI (adj. p = 2.74 × 10-3) of HM13 (Histocompatibility Minor 13), a signal peptide peptidase involved in epitope generation. HM13 LOI samples featured HM13 overexpression, both compared to normal solid tissues (p = 3.00 × 10-7) and non-LOI (p = 1.27 × 10-2) samples. Upon adjustment for age and sex, HM13 expression was significantly associated with poor survival (p = 7.10 × 10-5). Moreover, HM13 overexpression consistently exacerbated with increasing tumor stage (p = 2.90 × 10-8). For IGF2, LOI was observed in normal solid tissues, but the prevalence did not increase in cancer. In conclusion, HM13 LOI is an early event in ccRCC, causing overexpression leading to poor prognosis.

Imprinted gene detection effectively improves the diagnostic accuracy for papillary thyroid carcinoma

BACKGROUND

Papillary thyroid carcinoma (PTC) is the most frequent histological type of thyroid carcinoma. Although an increasing number of diagnostic methods have recently been developed, the diagnosis of a few nodules is still unsatisfactory. Therefore, the present study aimed to develop and validate a comprehensive prediction model to optimize the diagnosis of PTC.

METHODS

A total of 152 thyroid nodules that were evaluated by postoperative pathological examination were included in the development and validation cohorts recruited from two centres between August 2019 and February 2022. Patient data, including general information, cytopathology, imprinted gene detection, and ultrasound features, were obtained to establish a prediction model for PTC. Multivariate logistic regression analysis with a bidirectional elimination approach was performed to identify the predictors and develop the model.

RESULTS

A comprehensive prediction model with predictors, such as component, microcalcification, imprinted gene detection, and cytopathology, was developed. The area under the curve (AUC), sensitivity, specificity, and accuracy of the developed model were 0.98, 97.0%, 89.5%, and 94.4%, respectively. The prediction model also showed satisfactory performance in both internal and external validations. Moreover, the novel method (imprinted gene detection) was demonstrated to play a role in improving the diagnosis of PTC.

CONCLUSION

The present study developed and validated a comprehensive prediction model for PTC, and a visualized nomogram based on the prediction model was provided for clinical application. The prediction model with imprinted gene detection effectively improves the diagnosis of PTCs that are undetermined by the current means.

Tumor Microenvironment and Epigenetic Implications in Breast Cancer Progression

Breast cancer (BC) poses significant challenges, driven by its diverse nature and intricate dynamics. Epigenetic modifications, such as DNA methylation, histone modifications, and noncoding RNAs, have emerged as key regulators of gene expression and BC metastasis plasticity or therapeutic resistance. Targeting epigenetic regulators and pathways associated with therapeutic resistance holds promise for overcoming treatment obstacles and enhancing treatment efficacy.

Identification and validation of a novel necroptosis-related molecular signature to evaluate prognosis and immune features in breast cancer

Necroptosis has been shown to play an important role in the development of tumors. However, the characteristics of the necroptosis-related subtypes and the associated immune cell infiltration in the tumor microenvironment (TME) of breast cancer (BRCA) remain unclear. In this study, we identified three clusters related to necroptosis using the expression patterns of necroptosis-relevant genes (NRGs), and found that these three clusters had different clinicopathological features, prognosis and immune cell infiltration in the TME. Cluster 2 was characterized by less infiltration of immune cells in the TME and was associated with a worse prognosis. Then, a necroptosis risk score (NRS) composed of 14 NRGs was constructed using the least absolute shrinkage and selection operator regression (LASSO) Cox regression method. Based on NRS, all BRCA patients in the TCGA datasets were classified into a low-risk group and a high-risk group. Patients in the low-risk group were characterized by longer overall survival (OS), lower mutation burden, and higher infiltration level of immune cells in the TME. Moreover, the NRS was significantly associated with chemotherapeutic drug sensitivity. Finally, the knockdown of VDAC1 reduced the proliferation and migration of BRCA cells, and promoted cell death induced by necroptosis inducer. This study identified a novel necroptosis-related subtype of BRCA, and a comprehensive analysis of NRGs in BRCA revealed its potential roles in prognosis, clinicopathological features, TME, chemotherapy, tumor proliferation, and tumor necroptosis. These results may improve our understanding of NRGs in BRCA and provide a reference for developing individualized therapeutic strategies.

Novel epigenetic molecular therapies for imprinting disorders

Genomic imprinting disorders are caused by the disruption of genomic imprinting processes leading to a deficit or increase of an active allele. Their unique molecular mechanisms underlying imprinted genes offer an opportunity to investigate epigenetic-based therapy for reactivation of an inactive allele or reduction of an active allele. Current treatments are based on managing symptoms, not targeting the molecular mechanisms underlying imprinting disorders. Here, we highlight molecular approaches of therapeutic candidates in preclinical and clinical studies for individual imprinting disorders. These include the significant progress of discovery and testing of small molecules, antisense oligonucleotides, and CRISPR mediated genome editing approaches as new therapeutic strategies. We discuss the significant challenges of translating these promising therapies from the preclinical stage to the clinic, especially for genome editing based approaches.

Imprinted Genes: Genomic Conservation, Transcriptomic Dynamics and Phenomic Significance in Health and Diseases

Since its discovery in 1991, genomic imprinting has been the subject of numerous studies into its mechanisms of establishment and regulation, evolution and function, and presence in multiple genomes. Disturbance of imprinting has been implicated in a range of diseases, ranging from debilitating syndromes to cancers to fetal deficiencies. Despite this, studies done on the prevalence and relevance of imprinting on genes have been limited in scope, tissue types available, and focus, by both availability and resources. This has left a gap in comparative studies. To address this, we assembled a collection of imprinted genes available in current literature covering five species. Here we sought to identify trends and motifs in the imprinted gene set (IGS) in three distinct arenas: evolutionary conservation, across-tissue expression, and health phenomics. Overall, we found that imprinted genes displayed less conservation and higher proportions of non-coding RNA while maintaining synteny. Maternally expressed genes (MEGs) and paternally expressed genes (PEGs) occupied distinct roles in tissue expression and biological pathway use, while imprinted genes collectively showed a broader tissue range, notable preference for tissue specific expression and limited gene pathways than comparable sex differentiation genes. Both human and murine imprinted genes showed the same clear phenotypic trends, that were distinct from those displayed by sex differentiation genes which were less involved in mental and nervous system disease. While both sets had representation across the genome, the IGS showed clearer clustering as expected, with PEGs significantly more represented than MEGs.

Co-Occurrence of Beckwith-Wiedemann Syndrome and Early-Onset Colorectal Cancer

CRC is an adult-onset carcinoma representing the third most common cancer and the second leading cause of cancer-related deaths in the world. EO-CRC (<45 years of age) accounts for 5% of the CRC cases and is associated with cancer-predisposing genetic factors in half of them. Here, we describe the case of a woman affected by BWSp who developed EO-CRC at age 27. To look for a possible molecular link between BWSp and EO-CRC, we analysed her whole-genome genetic and epigenetic profiles in blood, and peri-neoplastic and neoplastic colon tissues. The results revealed a general instability of the tumor genome, including copy number and methylation changes affecting genes of the WNT signaling pathway, CRC biomarkers and imprinted loci. At the germline level, two missense mutations predicted to be likely pathogenic were found in compound heterozygosity affecting the Cystic Fibrosis (CF) gene CFTR that has been recently classified as a tumor suppressor gene, whose dysregulation represents a severe risk factor for developing CRC. We also detected constitutional loss of methylation of the KCNQ1OT1:TSS-DMR that leads to bi-allelic expression of the lncRNA KCNQ1OT1 and BWSp. Our results support the hypothesis that the inherited CFTR mutations, together with constitutional loss of methylation of the KCNQ1OT1:TSS-DMR, initiate the tumorigenesis process. Further somatic genetic and epigenetic changes enhancing the activation of the WNT/beta-catenin pathway likely contributed to increase the growth advantage of cancer cells. Although this study does not provide any conclusive cause-effect relationship between BWSp and CRC, it is tempting to speculate that the imprinting defect of BWSp might accelerate tumorigenesis in adult cancer in the presence of predisposing genetic variants.

Differential haplotype expression in class I MHC genes during SARS-CoV-2 infection of human lung cell lines

Introduction

Cell entry of SARS-CoV-2 causes genome-wide disruption of the transcriptional profiles of genes and biological pathways involved in the pathogenesis of COVID-19. Expression allelic imbalance is characterized by a deviation from the Mendelian expected 1:1 expression ratio and is an important source of allele-specific heterogeneity. Expression allelic imbalance can be measured by allele-specific expression analysis (ASE) across heterozygous informative expressed single nucleotide variants (eSNVs). ASE reflects many regulatory biological phenomena that can be assessed by combining genome and transcriptome information. ASE contributes to the interindividual variability associated with the disease. We aim to estimate the transcriptome-wide impact of SARS-CoV-2 infection by analyzing eSNVs.

Methods

We compared ASE profiles in the human lung cell lines Calu-3, A459, and H522 before and after infection with SARS-CoV-2 using RNA-Seq experiments.

Results

We identified 34 differential ASE (DASE) sites in 13 genes (HLA-A, HLA-B, HLA-C, BRD2, EHD2, GFM2, GSPT1, HAVCR1, MAT2A, NQO2, SUPT6H, TNFRSF11A, UMPS), all of which are enriched in protein binding functions and play a role in COVID-19. Most DASE sites were assigned to the MHC class I locus and were predominantly upregulated upon infection. DASE sites in the MHC class I locus also occur in iPSC-derived airway epithelium basal cells infected with SARS-CoV-2. Using an RNA-Seq haplotype reconstruction approach, we found DASE sites and adjacent eSNVs in phase (i.e., predicted on the same DNA strand), demonstrating differential haplotype expression upon infection. We found a bias towards the expression of the HLA alleles with a higher binding affinity to SARS-CoV-2 epitopes.

Discussion

Independent of gene expression compensation, SARS-CoV-2 infection of human lung cell lines induces transcriptional allelic switching at the MHC loci. This suggests a response mechanism to SARS-CoV-2 infection that swaps HLA alleles with poor epitope binding affinity, an expectation supported by publicly available proteome data.

Phasing DNA Methylation

Haplotyping enables the study of allele-specific events. Heterozygous variants, primarily single nucleotide variants (SNVs), enable the assignment of the paternal and maternal origin of the chromosomes and are widely employed to phase sequencing reads to their haplotype of origin. Certain long-read technologies enable the detection of both the DNA sequence and DNA modifications. These long reads and their inherent methylation information are suitable for genome-wide haplotyping and allele-specific DNA methylation analysis. Here, we describe the workflow to phase reads and DNA methylation using nanopore sequencing.

Maximizing Small Biopsy Patient Samples: Unified RNA-Seq Platform Assessment of over 120,000 Patient Biopsies

Despite its wide-ranging benefits, whole-transcriptome or RNA exome profiling is challenging to implement in a clinical diagnostic setting. The Unified Assay is a comprehensive workflow wherein exome-enriched RNA-sequencing (RNA-Seq) assays are performed on clinical samples and analyzed by a series of advanced machine learning-based classifiers. Gene expression signatures and rare and/or novel genomic events, including fusions, mitochondrial variants, and loss of heterozygosity were assessed using RNA-Seq data generated from 120,313 clinical samples across three clinical indications (thyroid cancer, lung cancer, and interstitial lung disease). Since its implementation, the data derived from the Unified Assay have allowed significantly more patients to avoid unnecessary diagnostic surgery and have played an important role in guiding follow-up decisions regarding treatment. Collectively, data from the Unified Assay show the utility of RNA-Seq and RNA expression signatures in the clinical laboratory, and their importance to the future of precision medicine.

Increased copy number of imprinted genes in the chromosomal region 20q11-q13.32 is associated with resistance to antitumor agents in cancer cell lines

BACKGROUND

Parent of origin-specific allelic expression of imprinted genes is epigenetically controlled. In cancer, imprinted genes undergo both genomic and epigenomic alterations, including frequent copy number changes. We investigated whether copy number loss or gain of imprinted genes in cancer cell lines is associated with response to chemotherapy treatment.

RESULTS

We analyzed 198 human imprinted genes including protein-coding genes and noncoding RNA genes using data from tumor cell lines from the Cancer Cell Line Encyclopedia and Genomics of Drug Sensitivity in Cancer datasets. We examined whether copy number of the imprinted genes in 35 different genome locations was associated with response to cancer drug treatment. We also analyzed associations of pretreatment expression and DNA methylation of imprinted genes with drug response. Higher copy number of BLCAP, GNAS, NNAT, GNAS-AS1, HM13, MIR296, MIR298, and PSIMCT-1 in the chromosomal region 20q11-q13.32 was associated with resistance to multiple antitumor agents. Increased expression of BLCAP and HM13 was also associated with drug resistance, whereas higher methylation of gene regions of BLCAP, NNAT, SGK2, and GNAS was associated with drug sensitivity. While expression and methylation of imprinted genes in several other chromosomal regions was also associated with drug response and many imprinted genes in different chromosomal locations showed a considerable copy number variation, only imprinted genes at 20q11-q13.32 had a consistent association of their copy number with drug response. Copy number values among the imprinted genes in the 20q11-q13.32 region were strongly correlated. They were also correlated with the copy number of cancer-related non-imprinted genes MYBL2, AURKA, and ZNF217 in that chromosomal region. Expression of genes at 20q11-q13.32 was associated with ex vivo drug response in primary tumor samples from the Beat AML 1.0 acute myeloid leukemia patient cohort. Association of the increased copy number of the 20q11-q13.32 region with drug resistance may be complex and could involve multiple genes.

CONCLUSIONS

Copy number of imprinted and non-imprinted genes in the chromosomal region 20q11-q13.32 was associated with cancer drug resistance. The genes in this chromosomal region may have a modulating effect on tumor response to chemotherapy.

Overexpressed Histocompatibility Minor 13 was Associated with Liver Hepatocellular Carcinoma Progression and Prognosis

Among primary liver carcinoma cases, the proportion of liver hepatocellular carcinoma (LIHC) cases is 75%-85%. Current treatments for LIHC include chemotherapy, surgical excision, and liver transplantation, which are effective for early LIHC treatment. Nevertheless, the early symptoms of liver carcinoma are atypical, so a large proportion of LIHC patients are diagnosed at an advanced stage. Histocompatibility minor 13 (HM13), located in the endoplasmic reticulum, is responsible for catalysing the hydrolysis of some signal peptides after cleavage from the precursor protein. Here, we studied the role of HM13 in LIHC development through bioinformatics analysis. Database analysis showed that HM13 was of great significance for LIHC tumorigenesis. Compared to normal liver tissues, HM13 expression was increased to a greater extent in LIHC tissues. After analysis of Kaplan‒Meier plotter and Gene Expression Profiling Interactive Analysis (GEPIA) datasets, we discovered that highly expressed HM13 exhibited an association with shorter overall survival (OS), disease-free survival (DFS), and disease-specific survival (DSS). We conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses to analyse HM13-related genes, and the data indicated that these genes obviously participated in rRNA processing, ribosome biogenesis, spliceosome, Huntington's disease, and ATP-dependent helicase activity. The Cell Counting Kit-8 (CCK-8) assay and Transwell assay showed that reducing HM13 expression hindered LIHC cell proliferation, migration, and invasion. In conclusion, these findings indicate that HM13 is a biomarker and is related to the poor prognosis of LIHC. Our results are conducive to discovering new targets for LIHC treatment.

Pan-cancer analysis suggests histocompatibility minor 13 is an unfavorable prognostic biomarker promoting cell proliferation, migration, and invasion in hepatocellular carcinoma

Histocompatibility Minor 13 (HM13) encoding the signal peptide peptidase plays an important role in maintaining protein homeostasis but its role in tumors remains unclear. In this study, 33 tumor RNA-seq datasets were extracted from The Cancer Genome Atlas (TCGA) database, and the pan-cancer expression profile of HM13 was evaluated in combination with The Genotype-Tissue Expression (GTEx) datasets. The prognostic significance of abnormal HM13 pan-cancer expression was evaluated by univariate Cox regression and Kaplan-Meier analyses. Co-expression analysis was performed to examine the correlation between abnormal pan-cancer expression of HM13 and immune cell infiltration, immune checkpoint, molecules related to RNA modification, tumor mutational burden (TMB), microsatellite instability (MSI), and other related molecules. CellMiner database was used to evaluate the relationship between the expression of HM13 and drug sensitivity. The results showed overexpression of HM13 in almost all tumors except kidney chromophobe (KICH). Abnormally high expression of HM13 in adrenocortical carcinoma (ACC), kidney renal papillary cell carcinoma (KIRP), uveal melanoma (UVM), liver hepatocellular carcinoma (LIHC), brain lower grade glioma (LGG), head and neck squamous cell carcinoma (HNSC), and kidney renal clear cell carcinoma (KIRC) was associated with poor prognosis. Expression of HM13 correlated strongly with pan-cancer immune checkpoint gene expression and immune cell infiltration. Drug sensitivity analysis indicated that the expression of HM13 was an excellent predictor of drug sensitivity. We verified that both mRNA and protein levels of HM13 were abnormally upregulated in HCC tissues, and were independent risk factors for poor prognosis. Furthermore, interference with HM13 expression in Huh-7 and HCCLM3 cells significantly inhibited proliferation, migration, and invasion. Therefore, our findings demonstrate that HM13 is a potential pan-cancer prognostic marker, thus providing a new dimension for understanding tumor development.

An Interplay between Epigenetics and Translation in Oocyte Maturation and Embryo Development: Assisted Reproduction Perspective

Germ cell quality is a key prerequisite for successful fertilization and early embryo development. The quality is determined by the fine regulation of transcriptomic and proteomic profiles, which are prone to alteration by assisted reproduction technology (ART)-introduced in vitro methods. Gaining evidence shows the ART can influence preset epigenetic modifications within cultured oocytes or early embryos and affect their developmental competency. The aim of this review is to describe ART-determined epigenetic changes related to the oogenesis, early embryogenesis, and further in utero development. We confront the latest epigenetic, related epitranscriptomic, and translational regulation findings with the processes of meiotic maturation, fertilization, and early embryogenesis that impact the developmental competency and embryo quality. Post-ART embryo transfer, in utero implantation, and development (placentation, fetal development) are influenced by environmental and lifestyle factors. The review is emphasizing their epigenetic and ART contribution to fetal development. An epigenetic parallel among mouse, porcine, and bovine animal models and human ART is drawn to illustrate possible future mechanisms of infertility management as well as increase the awareness of the underlying mechanisms governing oocyte and embryo developmental complexity under ART conditions.

Genome-wide detection of imprinted differentially methylated regions using nanopore sequencing

Imprinting is a critical part of normal embryonic development in mammals, controlled by defined parent-of-origin (PofO) differentially methylated regions (DMRs) known as imprinting control regions. Direct nanopore sequencing of DNA provides a means to detect allelic methylation and to overcome the drawbacks of methylation array and short-read technologies. Here, we used publicly available nanopore sequencing data for 12 standard B-lymphocyte cell lines to acquire the genome-wide mapping of imprinted intervals in humans. Using the sequencing data, we were able to phase 95% of the human methylome and detect 94% of the previously well-characterized, imprinted DMRs. In addition, we found 42 novel imprinted DMRs (16 germline and 26 somatic), which were confirmed using whole-genome bisulfite sequencing (WGBS) data. Analysis of WGBS data in mouse (Mus musculus), rhesus monkey (Macaca mulatta), and chimpanzee (Pan troglodytes) suggested that 17 of these imprinted DMRs are conserved. Some of the novel imprinted intervals are within or close to imprinted genes without a known DMR. We also detected subtle parental methylation bias, spanning several kilobases at seven known imprinted clusters. At these blocks, hypermethylation occurs at the gene body of expressed allele(s) with mutually exclusive H3K36me3 and H3K27me3 allelic histone marks. These results expand upon our current knowledge of imprinting and the potential of nanopore sequencing to identify imprinting regions using only parent-offspring trios, as opposed to the large multi-generational pedigrees that have previously been required.

GALNT8 suppresses breast cancer cell metastasis potential by regulating EGFR O-GalNAcylation

Breast cancer represents the most lethal malignancy that threatens the health of females. Metastasis is the fatal hallmark of breast cancer, and current effective therapeutic targets of metastasis are still lacking. Aberrant O-GalNAcylation, which is attributed to alteration of polypeptide N-acetylgalactosaminyl transferases (GALNTs), has been implicated in cancer metastasis. However, GALNTs that drive metastasis in breast cancer and their underlying mechanisms are largely unclear. In the present study, a negative correlation between GALNT8 and the prognosis of breast cancer patients was observed in multiple groups of Gene Expression Omnibus (GEO) datasets. We then constructed a stable GALNT8 knockdown MCF7 cell line and performed transcriptome analysis using RNA sequencing, which revealed that the expression of multiple migration-related genes was changed. GALNT8 was identified as a regulator of epithelial-mesenchymal transition (EMT) markers, including E-cadherin, N-cadherin, ZO-1 and vimentin. Moreover, loss- and gain-of-function GALNT8 assays demonstrated that this glycosyltransferase inhibited the metastatic potential of breast cancer cells. Interestingly, the O-GalNAcylation of EGFR, which is the key factor related to the metastasis cascade, was impacted by GALNT8. Furthermore, our results suggested that the GALNT8-mediated O-GalNAcylation led to the suppression of the EGFR signaling pathway and metastatic potential in breast cancer cells. These results suggested that GALNT8 acts as a tumor suppressor, represses tumor metastasis and inhibits the EMT process through the EGFR signaling pathway. This finding may provide insight into the mechanism by which aberrant O-glycosylation modulates breast cancer metastasis.

Genomic Imprinting in the New Omics Era: A Model for Systems-Level Approaches

Genomic imprinting represents a noteworthy inheritance mechanism leading to allele-specific regulations dependent of the parental origin. Imprinted loci are especially involved in essential mammalian functions related to growth, development and behavior. In this mini-review, we first offer a summary of current representations associated with genomic imprinting through key results of the three last decades. We then outline new perspectives allowed by the spread of new omics technologies tackling various interacting levels of imprinting regulations, including genomics, transcriptomics and epigenomics. We finally discuss the expected contribution of new omics data to unresolved big questions in the field.

BrewerIX enables allelic expression analysis of imprinted and X-linked genes from bulk and single-cell transcriptomes

Genomic imprinting and X chromosome inactivation (XCI) are two prototypical epigenetic mechanisms whereby a set of genes is expressed mono-allelically in order to fine-tune their expression levels. Defects in genomic imprinting have been observed in several neurodevelopmental disorders, in a wide range of tumours and in induced pluripotent stem cells (iPSCs). Single Nucleotide Variants (SNVs) are readily detectable by RNA-sequencing allowing the determination of whether imprinted or X-linked genes are aberrantly expressed from both alleles, although standardised analysis methods are still missing. We have developed a tool, named BrewerIX, that provides comprehensive information about the allelic expression of a large, manually-curated set of imprinted and X-linked genes. BrewerIX does not require programming skills, runs on a standard personal computer, and can analyze both bulk and single-cell transcriptomes of human and mouse cells directly from raw sequencing data. BrewerIX confirmed previous observations regarding the bi-allelic expression of some imprinted genes in naive pluripotent cells and extended them to preimplantation embryos. BrewerIX also identified misregulated imprinted genes in breast cancer cells and in human organoids and identified genes escaping XCI in human somatic cells. We believe BrewerIX will be useful for the study of genomic imprinting and XCI during development and reprogramming, and for detecting aberrations in cancer, iPSCs and organoids. Due to its ease of use to non-computational biologists, its implementation could become standard practice during sample assessment, thus raising the robustness and reproducibility of future studies.

Epigenetic imprinting alterations as effective diagnostic biomarkers for early-stage lung cancer and small pulmonary nodules

Background

Early lung cancer detection remains a clinical challenge for standard diagnostic biopsies due to insufficient tumor morphological evidence. As epigenetic alterations precede morphological changes, expression alterations of certain imprinted genes could serve as actionable diagnostic biomarkers for malignant lung lesions.

Results

Using the previously established quantitative chromogenic imprinted gene in situ hybridization (QCIGISH) method, elevated aberrant allelic expression of imprinted genes GNAS, GRB10, SNRPN and HM13 was observed in lung cancers over benign lesions and normal controls, which were pathologically confirmed among histologically stained normal, paracancerous and malignant tissue sections. Based on the differential imprinting signatures, a diagnostic grading model was built on 246 formalin-fixed and paraffin-embedded (FFPE) surgically resected lung tissue specimens, tested against 30 lung cytology and small biopsy specimens, and blindly validated in an independent cohort of 155 patients. The QCIGISH diagnostic model demonstrated 99.1% sensitivity (95% CI 97.5–100.0%) and 92.1% specificity (95% CI 83.5–100.0%) in the blinded validation set. Of particular importance, QCIGISH achieved 97.1% sensitivity (95% CI 91.6–100.0%) for carcinoma in situ to stage IB cancers with 100% sensitivity and 91.7% specificity (95% CI 76.0–100.0%) noted for pulmonary nodules with diameters ≤ 2 cm.

Conclusions

Our findings demonstrated the diagnostic value of epigenetic imprinting alterations as highly accurate translational biomarkers for a more definitive diagnosis of suspicious lung lesions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01203-5.

The Role of Long Non-Coding RNAs (lncRNAs) in Female Oriented Cancers

Recent advances in molecular biology have discovered the mysterious role of long non-coding RNAs (lncRNAs) as potential biomarkers for cancer diagnosis and targets for advanced cancer therapy. Studies have shown that lncRNAs take part in the incidence and development of cancers in humans. However, previously they were considered as mere RNA noise or transcription byproducts lacking any biological function. In this article, we present a summary of the progress on ascertaining the biological functions of five lncRNAs (HOTAIR, NEAT1, H19, MALAT1, and MEG3) in female-oriented cancers, including breast and gynecological cancers, with the perspective of carcinogenesis, cancer proliferation, and metastasis. We provide the current state of knowledge from the past five years of the literature to discuss the clinical importance of such lncRNAs as therapeutic targets or early diagnostic biomarkers. We reviewed the consequences, either oncogenic or tumor-suppressing features, of their aberrant expression in female-oriented cancers. We tried to explain the established mechanism by which they regulate cancer proliferation and metastasis by competing with miRNAs and other mechanisms involved via regulating genes and signaling pathways. In addition, we revealed the association between stated lncRNAs and chemo-resistance or radio-resistance and their potential clinical applications and future perspectives.

Prognostic prediction and gene regulation network of EIF2S2 in hepatocellular carcinoma based on data mining

BACKGROUND

Hepatocellular carcinoma (HCC) is a malignant tumor with a high fatality rate, predicting poor prognosis and therapeutic effect. Screening potential prognostic genes in HCC could be a creative way to advance clinical treatment. Eukaryotic translation initiation factor 2 subunit beta (EIF2S2) has reportedly been linked to several tumors, including liver cancer, but the prognostic predictions remain unknown. Therefore, we aimed to clarify the prognostic role and interaction network of EIF2S2 in HCC using bioinformatics data.

METHODS

We screened EIF2S2 using the Oncomine, Ualcan, and TCGA databases. R software was used to analyze the mRNA level and clinicopathological characteristics of hepatocellular carcinoma. Evaluation of the correlations between EIF2S2 and patients' survival was made using the Kaplan-Meier curves and Cox proportional hazards regression model. Then, the influence of EIF2S2 gene mutations on the prognosis of patients was explored by cBioPortal. The protein-protein interaction network of 50 similar genes related to EIF2S2 was implemented by GEPIA2 and Metascape. The LinkedOmics database allowed us to carry out Gene Set Enrichment Analysis. Finally, we constructed the EIF2S2 kinase, miRNA, and transcription factor target networks using GeneMANIA.

RESULTS

EIF2S2 mRNA was overexpressed in HCC and was closely associated with clinicopathological features, including gender, age, race, tumor grade, and stage. There was no correlation between EIF2S2 genetic mutations and prognostic survival. Combining Cox proportional hazards regression model analyses, high-expressed EIF2S2 predicted poor prognosis in HCC patients. Additionally, we screened the top three EIF2S2-related genes (PFDN4, HM13, and SNRPD1), the 50 similar genes, and then constructed a 50-similar-gene protein-protein interaction network identified by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways using Metascape. EIF2S2 target networks in HCC were identified in kinase, miRNA, and transcription factor networks, including the mitogen-activated protein kinase 1 (MAPK1), miRNAs (Mir-144), and transcription factors (GGAANCGGAANY_UNKNOWN) using GeneMANIA.

CONCLUSIONS

EIF2S2 plays a crucial role in the gene-regulating network of HCC and may be a potential prognostic marker or therapeutic target for HCC patients.

The Role of Long Non-coding RNAs in Human Imprinting Disorders: Prospective Therapeutic Targets

Genomic imprinting is a term used for an intergenerational epigenetic inheritance and involves a subset of genes expressed in a parent-of-origin-dependent way. Imprinted genes are expressed preferentially from either the paternally or maternally inherited allele. Long non-coding RNAs play essential roles in regulating this allele-specific expression. In several well-studied imprinting clusters, long non-coding RNAs have been found to be essential in regulating temporal- and spatial-specific establishment and maintenance of imprinting patterns. Furthermore, recent insights into the epigenetic pathological mechanisms underlying human genomic imprinting disorders suggest that allele-specific expressed imprinted long non-coding RNAs serve as an upstream regulator of the expression of other protein-coding or non-coding imprinted genes in the same cluster. Aberrantly expressed long non-coding RNAs result in bi-allelic expression or silencing of neighboring imprinted genes. Here, we review the emerging roles of long non-coding RNAs in regulating the expression of imprinted genes, especially in human imprinting disorders, and discuss three strategies targeting the central long non-coding RNA UBE3A-ATS for the purpose of developing therapies for the imprinting disorders Prader-Willi syndrome and Angelman syndrome. In summary, a better understanding of long non-coding RNA-related mechanisms is key to the development of potential therapeutic targets for human imprinting disorders.

ASMdb: a comprehensive database for allele-specific DNA methylation in diverse organisms

DNA methylation is known to be the most stable epigenetic modification and has been extensively studied in relation to cell differentiation, development, X chromosome inactivation and disease. Allele-specific DNA methylation (ASM) is a well-established mechanism for genomic imprinting and regulates imprinted gene expression. Previous studies have confirmed that certain special regions with ASM are susceptible and closely related to human carcinogenesis and plant development. In addition, recent studies have proven ASM to be an effective tumour marker. However, research on the functions of ASM in diseases and development is still extremely scarce. Here, we collected 4400 BS-Seq datasets and 1598 corresponding RNA-Seq datasets from 47 species, including human and mouse, to establish a comprehensive ASM database. We obtained the data on DNA methylation level, ASM and allele-specific expressed genes (ASEGs) and further analysed the ASM/ASEG distribution patterns of these species. In-depth ASM distribution analysis and differential methylation analysis conducted in nine cancer types showed results consistent with the reported changes in ASM in key tumour genes and revealed several potential ASM tumour-related genes. Finally, integrating these results, we constructed the first well-resourced and comprehensive ASM database for 47 species (ASMdb, www.dna-asmdb.com).

Molecular epigenetic dynamics in breast carcinogenesis

Breast cancer has become one of the most common dreadful diseases that target women across the globe. The most obvious reasons we associate with it are either genetic mutations or dysregulation of pathways. However, there is yet another domain that has a significant role in influencing the genetic mutations and pathways. Epigenetic mechanisms influence these pathways either independently or in association with genetic mutations, thereby expediting the process of breast carcinogenesis. Breast cancer is governed by various transduction pathways such as PI3K/AKT/mTOR, NOTCH, β Catenin, NF-kB, Hedgehog, etc. There are many proteins as well that serve to be tumor suppressors but somehow lose their ability to function. This may be because of either genetic mutation or a process that represses their function. Apart from these, there are a lot of individual factors like puberty, breastfeeding, abortion, parity, circadian rhythm, alcohol consumption, pollutants, and obesity that drive these mutations and hence alter the pathways. Epigenetic mechanisms like DNA methylation, histone modifications, and lncRNAs directly or indirectly bring alterations in the proteins that are involved in the pathways. They do this by either promoting the transcription of genes or by repressing it at the ground genetic level that advances breast carcinogenesis. Epigenetics precedes genetic mutation in driving carcinogenesis and so, it needs to be explored further to diversify the possibilities of target specific treatments. In this review, the general role of DNA methylation, histone modification, and lncRNAs in breast cancer and their role in influencing the oncogenic signaling pathways along with the various factors governing them have been discussed for a better understanding of the role of epigenetics in breast carcinogenesis.

DNA methylation and histone variants in aging and cancer

Aging-related diseases such as cancer can be traced to the accumulation of molecular disorder including increased DNA mutations and epigenetic drift. We provide a comprehensive review of recent results in mice and humans on modifications of DNA methylation and histone variants during aging and in cancer. Accumulated errors in DNA methylation maintenance lead to global decreases in DNA methylation with relaxed repression of repeated DNA and focal hypermethylation blocking the expression of tumor suppressor genes. Epigenetic clocks based on quantifying levels of DNA methylation at specific genomic sites is proving to be a valuable metric for estimating the biological age of individuals. Histone variants have specialized functions in transcriptional regulation and genome stability. Their concentration tends to increase in aged post-mitotic chromatin, but their effects in cancer are mainly determined by their specialized functions. Our increased understanding of epigenetic regulation and their modifications during aging has motivated interventions to delay or reverse epigenetic modifications using the epigenetic clocks as a rapid readout for efficacity. Similarly, the knowledge of epigenetic modifications in cancer is suggesting new approaches to target these modifications for cancer therapy.

The RNA Atlas expands the catalog of human non-coding RNAs

Existing compendia of non-coding RNA (ncRNA) are incomplete, in part because they are derived almost exclusively from small and polyadenylated RNAs. Here we present a more comprehensive atlas of the human transcriptome, which includes small and polyA RNA as well as total RNA from 300 human tissues and cell lines. We report thousands of previously uncharacterized RNAs, increasing the number of documented ncRNAs by approximately 8%. To infer functional regulation by known and newly characterized ncRNAs, we exploited pre-mRNA abundance estimates from total RNA sequencing, revealing 316 microRNAs and 3,310 long non-coding RNAs with multiple lines of evidence for roles in regulating protein-coding genes and pathways. Our study both refines and expands the current catalog of human ncRNAs and their regulatory interactions. All data, analyses and results are available for download and interrogation in the R2 web portal, serving as a basis for future exploration of RNA biology and function.

Megabase-scale methylation phasing using nanopore long reads and NanoMethPhase

The ability of nanopore sequencing to simultaneously detect modified nucleotides while producing long reads makes it ideal for detecting and phasing allele-specific methylation. However, there is currently no complete software for detecting SNPs, phasing haplotypes, and mapping methylation to these from nanopore sequence data. Here, we present NanoMethPhase, a software tool to phase 5-methylcytosine from nanopore sequencing. We also present SNVoter, which can post-process nanopore SNV calls to improve accuracy in low coverage regions. Together, these tools can accurately detect allele-specific methylation genome-wide using nanopore sequence data with low coverage of about ten-fold redundancy.

The Prognostic Value and Mechanisms of TMEM16A in Human Cancer

As a calcium ion-dependent chloride channel transmembrane protein 16A (TMEM16A) locates on the cell membrane. Numerous research results have shown that TMEM16A is abnormally expressed in many cancers. Mechanically, TMEM16A participates in cancer proliferation and migration by affecting the MAPK and CAMK signaling pathways. Additionally, it is well documented that TMEM16A exerts a regulative impact on the hyperplasia of cancer cells by interacting with EGFR in head and neck squamous cell carcinoma (HNSCC), an epithelial growth factor receptor in head and neck squamous cell carcinoma respectively. Meanwhile, as an EGFR activator, TMEM16A is considered as an oncogene or a tumor-promoting factor. More and more experimental data showed that down-regulation of TMEM16A or gene targeted therapy may be an effective treatment for cancer. This review summarized its role in various cancers and research advances related to its clinical application included treatment and diagnosis.

High heterogeneity undermines generalization of differential expression results in RNA-Seq analysis

BACKGROUND

RNA sequencing (RNA-Seq) has been widely applied in oncology for monitoring transcriptome changes. However, the emerging problem that high variation of gene expression levels caused by tumor heterogeneity may affect the reproducibility of differential expression (DE) results has rarely been studied. Here, we investigated the reproducibility of DE results for any given number of biological replicates between 3 and 24 and explored why a great many differentially expressed genes (DEGs) were not reproducible.

RESULTS

Our findings demonstrate that poor reproducibility of DE results exists not only for small sample sizes, but also for relatively large sample sizes. Quite a few of the DEGs detected are specific to the samples in use, rather than genuinely differentially expressed under different conditions. Poor reproducibility of DE results is mainly caused by high variation of gene expression levels for the same gene in different samples. Even though biological variation may account for much of the high variation of gene expression levels, the effect of outlier count data also needs to be treated seriously, as outlier data severely interfere with DE analysis.

CONCLUSIONS

High heterogeneity exists not only in tumor tissue samples of each cancer type studied, but also in normal samples. High heterogeneity leads to poor reproducibility of DEGs, undermining generalization of differential expression results. Therefore, it is necessary to use large sample sizes (at least 10 if possible) in RNA-Seq experimental designs to reduce the impact of biological variability and DE results should be interpreted cautiously unless soundly validated.

Causes, effects, and clinical implications of perturbed patterns within the cancer epigenome

Somatic mutations accumulating over a patient's lifetime are well-defined causative factors that fuel carcinogenesis. It is now clear, however, that epigenomic signature is also largely perturbed in many malignancies. These alterations support the transcriptional program crucial for the acquisition and maintenance of cancer hallmarks. Epigenetic instability may arise due to the genetic mutations or transcriptional deregulation of the proteins implicated in epigenetic signaling. Moreover, external stimulation and physiological aging may also participate in this phenomenon. The epigenomic signature is frequently associated with a cell of origin, as well as with tumor stage and differentiation, which all reflect its high heterogeneity across and within various tumors. Here, we will overview the current understanding of the causes and effects of the altered and heterogeneous epigenomic landscape in cancer. We will focus mainly on DNA methylation and post-translational histone modifications as the key regulatory epigenetic signaling marks. In addition, we will describe how this knowledge is translated into the clinic. We will particularly concentrate on the applicability of epigenetic alterations as biomarkers for improved diagnosis, prognosis, and prediction. Finally, we will also review current developments regarding epi-drug usage in clinical and experimental settings.

Both Epimutations and Chromosome Aberrations Affect Multiple Imprinted Loci in Aggressive Wilms Tumors

The embryonal renal cancer Wilms tumor (WT) accounts for 7% of all children's malignancies. Its most frequent molecular defect is represented by DNA methylation abnormalities at the imprinted 11p15.5 region. Multiple imprinted methylation alterations dictated by chromosome copy-number variations have been recently demonstrated in adult cancers, raising the question of whether multiple imprinted loci were also affected in WT. To address this issue, we analyzed DNA methylation and chromosome profiles of 7 imprinted loci in 48 WT samples. The results demonstrated that methylation abnormalities of multiple imprinted loci occurred in 35% of the cases, but that they were associated with either chromosome aberrations or normal chromosome profiles. Multiple imprinted methylation changes were correlated with tumor stage and presence of metastasis, indicating that these epimutations were more frequent in highly aggressive tumors. When chromosome profiles were affected, these alterations were extended to flanking cancer driver genes. Overall, this study demonstrates the presence of multiple imprinted methylation defects in aggressive WTs and suggests that the mechanism by which they arise in embryonal and adult cancers is different.

IGF2 loss of imprinting enhances colorectal cancer stem cells pluripotency by promoting tumor autophagy

Cancer stem cells (CSCs) are believed to be the driving force behind the tumor growth. We performed this study to further explore the role of IGF2 epigenetic on CRC stem cells pluripotency which showed that IGF2 LOI CRC cells usually had a higher CD133 expression and sphere forming efficiency than MOI cells. IGF2 LOI CSCs were also found to have a higher level of autophagy than MOI CSCs. Moreover, IGF2/IR-A signal was determined to play a more important role in CSCs formation than IGF2/IGF1R. At last, by using miRNA-195 mimics, we fortunately found the increased IR-A expression might be due to the degradation of miRNA-195 in CRC. In conclusion, our results might reveal that IGF2 LOI could promote CRC stem cells pluripotency by promoting CSCs autophagy. For the degradation of miRNA-195, IGF2 showed a higher ability in interacting with overexpressed IR-A rather than IGF1R which would further activate CSCs autophagy. All these findings might provide a novel mechanistic insight into CRC diagnosis and therapy.

Gene expression of adipokines and adipokine receptors in the tumor microenvironment: associations of lower expression with more aggressive breast tumor features

PURPOSE

Limited epidemiologic data are available on the expression of adipokines leptin (LEP) and adiponectin (ADIPOQ) and adipokine receptors (LEPR, ADIPOR1, ADIPOR2) in the breast tumor microenvironment (TME). The associations of gene expression of these biomarkers with tumor clinicopathology are not well understood.

METHODS

NanoString multiplexed assays were used to assess the gene expression levels of LEP, LEPR, ADIPOQ, ADIPOR1, and ADIPOR2 within tumor tissues among 162 Black and 55 White women with newly diagnosed breast cancer. Multivariate mixed effects models were used to estimate associations of gene expression with breast tumor clinicopathology (overall and separately among Blacks).

RESULTS

Black race was associated with lower gene expression of LEPR (P = 0.002) and ADIPOR1 (P = 0.01). Lower LEP, LEPR, and ADIPOQ gene expression were associated with higher tumor grade (P = 0.0007, P < 0.0001, and P < 0.0001, respectively) and larger tumor size (P < 0.0001, P = 0.0005, and P < 0.0001, respectively). Lower ADIPOQ expression was associated with ER- status (P = 0.0005), and HER2-enriched (HER2-E; P = 0.0003) and triple-negative (TN; P = 0.002) subtypes. Lower ADIPOR2 expression was associated with Ki67+ status (P = 0.0002), ER- status (P < 0.0001), PR- status (P < 0.0001), and TN subtype (P = 0.0002). Associations of lower adipokine and adipokine receptor gene expression with ER-, HER2-E, and TN subtypes were confirmed using data from The Cancer Genome Atlas (P-values < 0.005).

CONCLUSION

These findings suggest that lower expression of ADIPOQ, ADIPOR2, LEP, and LEPR in the breast TME might be indicators of more aggressive breast cancer phenotypes. Validation of these findings are warranted to elucidate the role of the adipokines and adipokine receptors in long-term breast cancer prognosis.

The Roles of Imprinted SLC22A18 and SLC22A18AS Gene Overexpression Caused by Promoter CpG Island Hypomethylation as Diagnostic and Prognostic Biomarkers for Non-Small Cell Lung Cancer Patients

Genomic imprinting is a process that involves one gene copy turned-off in a parent-of-origin-dependent manner. The regulation of imprinted genes is broadly dependent on promoter methylation marks, which are frequently associated with both oncogenes and tumor suppressors. The purpose of this study was to assess the DNA methylation patterns of the imprinted solute-carrier family 22 member 18 (SLC22A18) and SLC22A18 antisense (SLC22A18AS) genes in non-small cell lung cancer (NSCLC) patients to study their relevance to the disease. We found that both genes were hypomethylated in adenocarcinoma and squamous cell carcinoma patients. Due to this imprinting loss, SLC22A18 and SLC22A18AS were found to be overexpressed in NSCLC tissues, which is significantly more evident in lung adenocarcinoma patients. These results were validated through analyses of public databases of NSCLC patients. The reversed gene profile of both genes was achieved in vitro by treatment with ademetionine. We then showed that high SLC22A18 and SLC22A18AS expression levels were significantly associated with worsening disease progression. In addition, low levels of SLC22A18AS were also correlated with better overall survival for lung adenocarcinoma patients. We found that SLC22A18 and SLC22A18AS knockdown inhibits cell proliferation in vitro. All these results suggest that both genes may be useful as diagnostic and prognostic biomarkers in NSCLC, revealing novel therapeutic opportunities.

LncRNA MEG3 rs3087918 was associated with a decreased breast cancer risk in a Chinese population: a case-control study

BACKGROUND

LncRNA MEG3 expressed abnormally in various cancers including breast cancer, but no studies reported the correlation between MEG3 SNPs and breast cancer susceptibility among Chinese women.

METHODS

This study is aimed to explore the association between three SNPs of MEG3 (rs3087918, rs7158663, rs11160608) and breast cancer. The study is a population-based case-control study including 434 breast cancer patients and 700 healthy controls. Genotyping was performed using Sequenom MassArray technique. Function prediction of rs3087918 were based on RNAfold and lncRNASNP2 databases.

RESULTS

Pooled analysis indicated that rs3087918 was related to a decreased risk of breast cancer [GG vs. TT: OR (95%) = 0.67(0.45-0.99), P = 0.042; GG vs. TT + TG: OR (95%) = 0.69(0.48-0.99), P = 0.046], especially for women aged <=49 [GG vs. TT: OR (95%) = 0.40(0.22-0.73), P = 0.02]. Comparison between case groups showed genotype GG and TG/GG of rs3087918 were associated with her-2 receptor expression [GG vs. TT: OR (95%) = 2.37(1.24-4.63), P = 0.010; TG + GG vs. TT: OR (95%) = 1.50(1.01-2.24), P = 0.045]. We didn't find statistical significance for rs11160608, rs7158663 and breast cancer. Structure prediction based on RNAfold found rs3087918 may influence the secondary structure of MEG3. The results based on lncRNASNP2 indicated that rs3087918 may gain the targets of hsa-miR-1203 to MEG3, while loss the target of hsa-miR-139-3p and hsa-miR-5091 to MEG3.

CONCLUSIONS

MEG3 rs3087918 was associated with a decreased risk of breast cancer. MEG3 haplotype TCG may increase the risk of breast cancer.

The influence of DNA methylation on monoallelic expression

Monoallelic gene expression occurs in diploid cells when only one of the two alleles of a gene is active. There are three main classes of genes that display monoallelic expression in mammalian genomes: (1) imprinted genes that are monoallelically expressed in a parent-of-origin dependent manner; (2) X-linked genes that undergo random X-chromosome inactivation in female cells; (3) random monoallelically expressed single and clustered genes located on autosomes. The heritability of monoallelic expression patterns during cell divisions implies that epigenetic mechanisms are involved in the cellular memory of these expression states. Among these, methylation of CpG sites on DNA is one of the best described modification to explain somatic inheritance. Here, we discuss the relevance of DNA methylation for the establishment and maintenance of monoallelic expression patterns among these three groups of genes, and how this is intrinsically linked to development and cellular states.

Novel visualized quantitative epigenetic imprinted gene biomarkers diagnose the malignancy of ten cancer types

BACKGROUND

Epigenetic alterations are involved in most cancers, but its application in cancer diagnosis is still limited. More practical and intuitive methods to detect the aberrant expressions from clinical samples using highly sensitive biomarkers are needed. In this study, we developed a novel approach in identifying, visualizing, and quantifying the biallelic and multiallelic expressions of an imprinted gene panel associated with cancer status. We evaluated the normal and aberrant expressions measured using the imprinted gene panel to formulate diagnostic models which could accurately distinguish the imprinting differences of normal and benign cases from cancerous tissues for each of the ten cancer types.

RESULTS

The Quantitative Chromogenic Imprinted Gene In Situ Hybridization (QCIGISH) method developed from a 1013-case study which provides a visual and quantitative analysis of non-coding RNA allelic expressions identified the guanine nucleotide-binding protein, alpha-stimulating complex locus (GNAS), growth factor receptor-bound protein (GRB10), and small nuclear ribonucleoprotein polypeptide N (SNRPN) out of five tested imprinted genes as efficient epigenetic biomarkers for the early-stage detection of ten cancer types. A binary algorithm developed for cancer diagnosis showed that elevated biallelic expression (BAE), multiallelic expression (MAE), and total expression (TE) measurements for the imprinted gene panel were associated with cell carcinogenesis, with the formulated diagnostic models achieving consistently high sensitivities (91-98%) and specificities (86-98%) across the different cancer types.

CONCLUSIONS

The QCIGISH method provides an innovative way to visually assess and quantitatively analyze individual cells for cancer potential extending from hyperplasia and dysplasia until carcinoma in situ and invasion, which effectively supplements standard clinical cytologic and histopathologic diagnosis for early cancer detection. In addition, the diagnostic models developed from the BAE, MAE, and TE measurements of the imprinted gene panel GNAS, GRB10, and SNRPN could provide important predictive information which are useful in early-stage cancer detection and personalized cancer management.

Landscape of genomic imprinting and its functions in the mouse mammary gland

Genomic imprinting is an epigenetic modification of DNA, whereby gene expression is restricted to either maternally or paternally inherited alleles. Imprinted genes (IGs) in the placenta and embryo are essential for growth regulation and nutrient supply. However, despite being an important nutrition delivery organ, studies on mammary gland genomic imprinting remain limited. In this study, we found that both the number of IGs and their expression levels decreased during development of the mouse mammary gland. IG expression was lineage-specific and related to mammary gland development and lactation. Meta-analysis of single-cell RNA sequencing data revealed that mammary gland IGs were co-expressed in a network that regulated cell stemness and differentiation, which was confirmed by our functional studies. Accordingly, our data indicated that IGs were essential for the self-renewal of mammary gland stem cells and IG decline was correlated with mammary gland maturity. Taken together, our findings revealed the importance of IGs in a poorly studied nutrition-related organ, i.e. the mammary gland, thus providing a reference for further studies on genomic imprinting.

Long non-coding RNA HUMT hypomethylation promotes lymphangiogenesis and metastasis via activating FOXK1 transcription in triple-negative breast cancer

BACKGROUND

Triple-negative breast cancer (TNBC) is the most malignant subtype of breast cancer with highly invasive ability and metastatic nature to the lymph nodes. Long non-coding RNAs (lncRNAs) have been widely explored in cancer tumorigenesis and progression. However, their roles in TNBC lymph node metastasis remains rarely studied.

METHODS

The expression of lncRNA highly upregulated in metastatic TNBC (HUMT) in cell lines and tissues was detected by quantitative real-time PCR (qRT-PCR) and in situ hybridization (ISH). RNA immunoprecipitation (RIP) and RNA pulldown were used to verify the interaction between lncRNA and protein. Chromatin immunoprecipitation (CHIP) and dCas9-gRNA-guided chromatin immunoprecipitation (dCas9-CHIP) were conducted to identify the specific binding site of HUMT-YBX1 complex. Western blot was used to detect the downstream of HUMT.

RESULTS

HUMT was significantly upregulated in lymph node invasive cells and predicted poorer clinical prognosis. Functional study indicated that HUMT promoted lymphangiogenesis and lymph node metastasis. Bioinformatic analysis and qRT-PCR showed that the high expression of HUMT was correlated with the hypomethylation status of its promoter region. Further, HUMT recruited Y-box binding protein 1 (YBX1) to form a novel transcription complex and activated the expression of forkhead box k1 (FOXK1), thus enhancing the expression of vascular endothelial growth factor C (VEGFC). The therapeutic value was further validated in patient-derived xenograft (PDX) models, and a combined marker panel exhibited a better prognostic value for TNBC in receiver operating characteristic (ROC) analysis.

CONCLUSIONS

Our study identified a novel TNBC lymph node metastasis-associated lncRNA, which promoted TNBC progression and indicated a novel biomarker and potential therapeutic target for TNBC lymph node metastasis.

The Emerging Role of Ten-Eleven Translocation 1 in Epigenetic Responses to Environmental Exposures

Mounting evidence from epidemiological studies and animal models has linked exposures to environmental factors to changes in epigenetic markers, especially in DNA methylation. These epigenetic changes may lead to dysregulation of molecular processes and functions and mediate the impact of environmental exposures in complex diseases. However, detailed molecular events that result in epigenetic changes following exposures remain unclear. Here, we review the emerging evidence supporting a critical role of ten-eleven translocation 1 (TET1) in mediating these processes. Targeting TET1 and its associated pathways may have therapeutic potential in alleviating negative impacts of environmental exposures, preventing and treating exposure-related diseases.

The evolution of genomic imprinting: Epigenetic control of mammary gland development and postnatal resource control

Genomic imprinting is an epigenetically regulated process leading to gene expression according to its parental origin. Imprinting is essential for prenatal growth and development, regulating nutritional resources to offspring, and contributing to a favored theory about the evolution of imprinting being due to a conflict between maternal and paternal genomes for the control of prenatal resources-the so-called kinship hypothesis. Genomic imprinting has been mainly studied during embryonic and placental development; however, maternal nutrient provisioning is not restricted to the prenatal period. In this context, the mammary gland acts at the maternal-offspring interface providing milk to the newborn. Maternal care including lactation supports the offspring, delivering nutrients and bioactive molecules protecting against infections and contributing to healthy organ development and immune maturation. The normal developmental cycle of the mammary gland-pregnancy, lactation, involution-is vital for this process, raising the question of whether genomic imprinting might also play a role in postnatal nutrient transfer by controlling mammary gland development. Characterizing the function and epigenetic regulation of imprinted genes in the mammary gland cycle may therefore provide novel insights into the evolution of imprinting since the offspring's paternal genome is absent from the mammary gland, in addition to increasing our knowledge of postnatal nutrition and its relation to life-long health. This article is categorized under: Developmental Biology > Developmental Processes in Health and Disease.

The phenotypic variations of multi-locus imprinting disturbances associated with maternal-effect variants of NLRP5 range from overt imprinting disorder to apparently healthy phenotype

BACKGROUND

A subset of individuals affected by imprinting disorders displays multi-locus imprinting disturbances (MLID). MLID has been associated with maternal-effect variants that alter the maintenance of methylation at germline-derived differentially methylated regions (gDMRs) in early embryogenesis. Pedigrees of individuals with MLID also include siblings with healthy phenotype. However, it is unknown if these healthy individuals have MLID themselves or if their methylation patterns differ from those associated with imprinting disorders, and in general, if MLID affects the clinical phenotype.

METHODS

We have investigated gDMR methylation by locus-specific and whole-genome analyses in a family with multiple pregnancy losses, a child with Beckwith-Wiedemann syndrome (BWS) and a further child with no clinical diagnosis of imprinting disorder or other pathologies.

RESULTS

We detected MLID with different methylation profiles in the BWS-affected and healthy siblings. Whole-exome sequencing demonstrated the presence of novel loss-of-function variants of NLRP5 in compound heterozygosity in the mother. The methylation profiles of the two siblings were compared with those of other cases with MLID and control groups by principal component analysis and unsupervised hierarchical clustering, but while their patterns were clearly separated from those of controls, we were unable to cluster those associated with specific clinical phenotypes among the MLID cases.

CONCLUSION

The identification of two novel maternal-effect variants of NLRP5 associated with poly-abortivity and MLID adds further evidence to the role of this gene in the maintenance of genomic imprinting in early embryos. Furthermore, our results demonstrate that within these pedigrees, MLID can also be present in the progeny with healthy phenotype, indicating that some sort of compensation occurs between altered imprinted loci in these individuals. The analysis of larger cohorts of patients with MLID is needed to formulate more accurate epigenotype-phenotype correlations.

Effects of reprogramming on genomic imprinting and the application of pluripotent stem cells

Pluripotent stem cells are considered to be the ideal candidates for cell-based therapies in humans. In this regard, both nuclear transfer embryonic stem (ntES) cells and induced pluripotent stem (iPS) cells are particularly advantageous because patient-specific autologous ntES and iPS cells can avoid immunorejection and other side effects that may be present in the allogenic pluripotent stem cells derived from unrelated sources. However, they have been found to contain deleterious genetic and epigenetic changes that may hinder their therapeutic applications. Indeed, deregulation of genomic imprinting has been frequently observed in reprogrammed ntES and iPS cells. We will survey the recent studies on genomic imprinting in pluripotent stem cells, particularly in iPS cells. In a previous study published about six years ago, genomic imprinting was found to be variably lost in mouse iPS clones. Intriguingly, de novo DNA methylation also occurred at the previously unmethylated imprinting control regions (ICRs) in a high percentage of iPS clones. These unexpected results were confirmed by a recent independent study with a similar approach. Since dysregulation of genomic imprinting can cause many human diseases including cancer and neurological disorders, these recent findings on genomic imprinting in reprogramming may have some implications for therapeutic applications of pluripotent stem cells.