Identification of valid reference housekeeping genes for gene expression analysis in tumor neovascularization studies

CDK5R1, GSE1, HSPG2 and WDFY3 as indirect epigenetic-sensitive genes in atrial fibrillation

BACKGROUND

Although mounting evidence supports that aberrant DNA methylation occurs in the hearts of patients with atrial fibrillation (AF), noninvasive epigenetic characterization of AF has not yet been defined.

METHODS

We investigated DNA methylome changes in peripheral blood CD4+ T cells isolated from 10 patients with AF relative to 11 healthy subjects (HS) who were enrolled in the DIANA clinical trial (NCT04371809) via reduced-representation bisulfite sequencing (RRBS).

RESULTS

An atrial-specific PPI network revealed 18 hub differentially methylated genes (DMGs), wherein ROC curve analysis revealed reasonable diagnostic performance of DNA methylation levels found within CDK5R1 (AUC = 0.76; p = 0.049), HSPG2 (AUC = 0.77; p = 0.038), WDFY3 (AUC = 0.78; p = 0.029), USP49 (AUC = 0.76; p = 0.049), GSE1 (AUC = 0.76; p = 0.049), AIFM1 (AUC = 0.76; p = 0.041), CDK5RAP2 (AUC = 0.81; p = 0.017), COL4A1 (AUC = 0.86; p < 0.001), SEPT8 (AUC = 0.90; p < 0.001), PFDN1 (AUC = 0.90; p < 0.01) and ACOT7 (AUC = 0.78; p = 0.032). Transcriptional profiling of the hub DMGs provided a significant overexpression of PSDM6 (p = 0.004), TFRC (p = 0.01), CDK5R1 (p < 0.001), HSPG2 (p = 0.01), WDFY3 (p < 0.001), USP49 (p = 0.004) and GSE1 (p = 0.021) in AF patients vs HS.

CONCLUSIONS

CDK5R1, GSE1, HSPG2 and WDFY3 resulted the best discriminatory genes both at methylation and gene expression level. Our results provide several candidate diagnostic biomarkers with the potential to advance precision medicine in AF.

RIZ2 at the crossroad of the EGF/EGFR signaling in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is the third most deadly and fourth most diagnosed cancer worldwide. Despite the progress in early diagnosis and advanced therapeutic options, CRC shows a poor prognosis with a 5 year survival rate of ~ 45%. PRDM2/RIZ, a member of PR/SET domain family (PRDM), expresses two main molecular variants, the PR-plus isoform (RIZ1) and the PR-minus (RIZ2). The imbalance in their expression levels in favor of RIZ2 is observed in many cancer types. The full length RIZ1 has been extensively investigated in several cancers where it acts as a tumor suppressor, whereas few studies have explored the RIZ2 oncogenic properties. PRDM2 is often target of frameshift mutations and aberrant DNA methylation in CRC. However, little is known about its role in CRC.

METHODS

We combined in-silico investigation of The Cancer Genome Atlas (TCGA) CRC datasets, cellular and molecular assays, transcriptome sequencing and functional annotation analysis to assess the role of RIZ2 in human CRC.

RESULTS

Our in-silico analysis on TCGA datasets confirmed that PRDM2 gene is frequently mutated and transcriptionally deregulated in CRC and revealed that a RIZ2 increase is highly correlated with a significant RIZ1 downregulation. Then, we assayed several CRC cell lines by qRT-PCR analysis for the main PRDM2 transcripts and selected DLD1 cell line, which showed the lowest RIZ2 levels. Therefore, we overexpressed RIZ2 in these cells to mimic TCGA datasets analysis results and consequently to assess the PRDM2/RIZ2 role in CRC. Data from RNA-seq disclosed that RIZ2 overexpression induced profound changes in CRC cell transcriptome via EGF pathway deregulation, suggesting that RIZ2 is involved in the EGF autocrine regulation of DLD1 cell behavior. Noteworthy, the forced RIZ2 expression increased cell viability, growth, colony formation, migration and organoid formation. These effects could be mediated by the release of high EGF levels by RIZ2 overexpressing DLD1 cells.

CONCLUSIONS

Our findings add novel insights on the putative RIZ2 tumor-promoting functions in CRC, although additional efforts are warranted to define the underlying molecular mechanism.

Exploring the putative role of PRDM1 and PRDM2 transcripts as mediators of T lymphocyte activation

BACKGROUND

T cell activation and programming from their naïve/resting state, characterized by widespread modifications in chromatin accessibility triggering extensive changes in transcriptional programs, is orchestrated by several cytokines and transcription regulators. PRDM1 and PRDM2 encode for proteins with PR/SET and zinc finger domains that control several biological processes, including cell differentiation, through epigenetic regulation of gene expression. Different transcripts leading to main protein isoforms with (PR +) or without (PR-) the PR/SET domain have been described. Although many studies have established the critical PRDM1 role in hematopoietic cell differentiation, maintenance and/or function, the single transcript contribution has not been investigated before. Otherwise, very few evidence is currently available on PRDM2. Here, we aimed to analyze the role of PRDM1 and PRDM2 different transcripts as mediators of T lymphocyte activation.

METHODS

We analyzed the transcription signature of the main variants from PRDM1 (BLIMP1a and BLIMP1b) and PRDM2 (RIZ1 and RIZ2) genes, in human T lymphocytes and Jurkat cells overexpressing PRDM2 cDNAs following activation through different signals.

RESULTS

T lymphocyte activation induced an early increase of RIZ2 and RIZ1 followed by BLIMP1b increase and finally by BLIMP1a increase. The "first" and the "second" signals shifted the balance towards the PR- forms for both genes. Interestingly, the PI3K signaling pathway modulated the RIZ1/RIZ2 ratio in favor of RIZ1 while the balance versus RIZ2 was promoted by MAPK pathway. Cytokines mediating different Jak/Stat signaling pathways (third signal) early modulated the expression of PRDM1 and PRDM2 and the relationship of their different transcripts confirming the early increase of the PR- transcripts. Different responses of T cell subpopulations were also observed. Jurkat cells showed that the acute transient RIZ2 increase promoted the balancing of PRDM1 forms towards BLIMP1b. The stable forced expression of RIZ1 or RIZ2 induced a significant variation in the expression of key transcription factors involved in T lymphocyte differentiation. The BLIMP1a/b balance shifted in favor of BLIMP1a in RIZ1-overexpressing cells and of BLIMP1b in RIZ2-overexpressing cells.

CONCLUSIONS

This study provides the first characterization of PRDM2 in T-lymphocyte activation/differentiation and novel insights on PRDM1 and PRDM2 transcription regulation during initial activation phases.

Association Between Circulating CD4+ T Cell Methylation Signatures of Network-Oriented SOCS3 Gene and Hemodynamics in Patients Suffering Pulmonary Arterial Hypertension

Pathogenic DNA methylation changes may be involved in pulmonary arterial hypertension (PAH) onset and its progression, but there is no data on potential associations with patient-derived hemodynamic parameters. The reduced representation bisulfite sequencing (RRBS) platform identified N = 631 differentially methylated CpG sites which annotated to N = 408 genes (DMGs) in circulating CD4⁺ T cells isolated from PAH patients vs. healthy controls (CTRLs). A promoter-restricted network analysis established the PAH subnetwork that included 5 hub DMGs (SOCS3, GNAS, ITGAL, NCOR2, NFIC) and 5 non-hub DMGs (NR4A2, GRM2, PGK1, STMN1, LIMS2). The functional analysis revealed that the SOCS3 gene was the most recurrent among the top ten significant pathways enriching the PAH subnetwork, including the growth hormone receptor and the interleukin-6 signaling. Correlation analysis showed that the promoter methylation levels of each network-oriented DMG were associated individually with hemodynamic parameters. In particular, SOCS3 hypomethylation was negatively associated with right atrial pressure (RAP) and positively associated with cardiac index (CI) (|r|≥ 0.6). A significant upregulation of the SOCS3, ITGAL, NFIC, NCOR2, and PGK1 mRNA levels (qRT-PCR) in peripheral blood mononuclear cells from PAH patients vs. CTRLs was found (P ≤ 0.05). By immunoblotting, a significant upregulation of the SOCS3 protein was confirmed in PAH patients vs. CTRLs (P < 0.01). This is the first network-oriented study which integrates circulating CD4⁺ T cell DNA methylation signatures, hemodynamic parameters, and validation experiments in PAH patients at first diagnosis or early follow-up. Our data suggests that SOCS3 gene might be involved in PAH pathogenesis and serve as potential prognostic biomarker.

Gene Expression Studies in Formalin-Fixed Paraffin-Embedded Samples of Cutaneous Cancer: The Need for Reference Genes

Formalin-fixed paraffin-embedded (FFPE) tumour samples may provide crucial data regarding biomarkers for neoplasm progression. Analysis of gene expression is frequently used for this purpose. Therefore, mRNA expression needs to be normalized through comparison to reference genes. In this study, we establish which of the usually reported reference genes is the most reliable one in cutaneous malignant melanoma (MM) and cutaneous squamous cell carcinoma (CSCC). ACTB, TFRC, HPRT1 and TBP expression was quantified in 123 FFPE samples (74 MM and 49 CSCC biopsies) using qPCR. Expression stability was analysed by NormFinder and Bestkeeper softwares, and the direct comparison method between means and SD. The in-silico analysis with BestKeeper indicated that HPRT1 was more stable than ACTB and TFRC in MM (1.85 vs. 2.15) and CSCC tissues (2.09 vs. 2.33). The best option to NormFinder was ACTB gene (0.56) in MM and TFRC (0.26) in CSCC. The direct comparison method showed lower SD means of ACTB expression in MM (1.17) and TFRC expression in CSCC samples (1.00). When analysing the combination of two reference genes for improving stability, NormFinder indicated HPRT1 and ACTB to be the best for MM samples, and HPRT1 and TFRC genes for CSCC. In conclusion, HPRT1 and ACTB genes in combination are the most appropriate choice for normalization in gene expression studies in MM FFPE tissue, while the combination of HPRT1 and TFRC genes are the best option in analysing CSCC FFPE samples. These may be used consistently in forthcoming studies on gene expression in both tumours.

ABCA1, TCF7, NFATC1, PRKCZ, and PDGFA DNA methylation as potential epigenetic-sensitive targets in acute coronary syndrome via network analysis

Acute coronary syndrome (ACS) is the most severe clinical manifestation of coronary heart disease.We performed an epigenome-wide analysis of circulating CD4⁺ and CD8⁺ T cells isolated from ACS patients and healthy subjects (HS), enrolled in the DIANA clinical trial, by reduced-representation bisulphite sequencing (RRBS). In CD4⁺ T cells, we identified 61 differentially methylated regions (DMRs) associated with 57 annotated genes (53% hyper- and 47% hypo-methylated) by comparing ACS patients vs HS. In CD8⁺ T cells, we identified 613 DMRs associated with 569 annotated genes (28% hyper- and 72% hypo-methylated) in ACS patients as compared to HS. In CD4⁺ vs CD8⁺ T cells of ACS patients we identified 175 statistically significant DMRs associated with 157 annotated genes (41% hyper- and 59% hypo-methylated). From pathway analyses, we selected six differentially methylated hub genes (NFATC1, TCF7, PDGFA, PRKCB, PRKCZ, ABCA1) and assessed their expression levels by q-RT-PCR. We found an up-regulation of selected genes in ACS patients vs HS (P < 0.001). ABCA1, TCF7, PDGFA, and PRKCZ gene expression was positively associated with CK-MB serum concentrations (r = 0.75, P = 0.03; r = 0.760, P = 0.029; r = 0.72, P = 0.044; r = 0.74, P = 0.035, respectively).This pilot study is the first single-base resolution map of DNA methylome by RRBS in CD4⁺ and CD8⁺ T cells and provides specific methylation signatures to clarify the role of aberrant methylation in ACS pathogenesis, thus supporting future research for novel epigenetic-sensitive biomarkers in the prevention and early diagnosis of this pathology.

Searching for a Putative Mechanism of RIZ2 Tumor-Promoting Function in Cancer Models

Positive Regulatory Domain (PRDM) gene family members commonly express two main molecular variants, the PR-plus isoform usually acting as tumor suppressor and the PR-minus one functioning as oncogene. Accordingly, PRDM2/RIZ encodes for RIZ1 (PR-plus) and RIZ2 (PR-minus). In human cancers, genetic or epigenetic modifications induce RIZ1 silencing with an expression level imbalance in favor of RIZ2 that could be relevant for tumorigenesis. Additionally, in estrogen target cells and tissues, estradiol increases RIZ2 expression level with concurrent increase of cell proliferation and survival. Several attempts to study RIZ2 function in HeLa or MCF-7 cells by its over-expression were unsuccessful. Thus, we over-expressed RIZ2 in HEK-293 cells, which are both RIZ1 and RIZ2 positive but unresponsive to estrogens. The forced RIZ2 expression increased cell viability and growth, prompted the G2-to-M phase transition and organoids formation. Accordingly, microarray analysis revealed that RIZ2 regulates several genes involved in mitosis. Consistently, RIZ silencing in both estrogen-responsive MCF-7 and -unresponsive MDA-MB-231 cells induced a reduction of cell proliferation and an increase of apoptosis rate. Our findings add novel insights on the putative RIZ2 tumor-promoting functions, although additional attempts are warranted to depict the underlying molecular mechanism.

Determination of reference genes for ovine pulmonary adenocarcinoma infected lung tissues using RNA-seq transcriptome profiling

Ovine pulmonary adenocarcinoma (OPA) is a globally occurring tumor of lung epithelium which seriously affects the development of sheep farming. In our research, lung tissues of 3 naturally infected OPA individuals and 3 healthy individuals (2-4 years old) were collected. RNA was extracted for transcriptome analysis and reference gene selection. According to transcriptome analysis, 7 candidate reference genes (eukaryotic translation initiation factor 1, EIF1; glyceraldehyde-3-phosphate dehydrogenase, GAPDH; beta-actin, ACTB; GABA Type A receptor-associated protein, GABARAP; activating transcription factor 4, ATF4; ribosomal protein S15, RPS15; and Y-Box binding protein 1, YBX1) showed fragments per kilobase of transcript per million fragments mapped (FPKM) values > 200.0 and standard errors of the means (SEM) < 20.0. Expression of the above candidate reference genes was evaluated by Real-time quantitative polymerase chain reaction (RT-qPCR) combined with the analysis using GeNorm, NormFinder, and BestKeeper software. Comprehensive analysis of the results showed that ACTB was the most stable one, followed by EIF1 and GABARAP. Then, expression stability of the above three genes were validated, suggesting as suitable reference genes in sheep lung tissue, in additional 30 OPA-affected lung tissues and 10 healthy ovine lung tissues. Finally, our findings will be helpful for the subsequent study on the tumorigenic mechanism of OPA.

Reliable reference genes for the quantification of mRNA in human T-cells and PBMCs stimulated with live influenza virus

Background

Quantitative PCR (qPCR) is a powerful tool that is particularly well-suited to measure mRNA levels in clinical samples, especially those with relatively low cell counts. However, a caveat of this approach is that reliable, stably expressed reference (housekeeping) genes are vital in order to ensure reproducibility and appropriate biological inference. In this study, we evaluated the expression stability of six reference genes in peripheral blood mononuclear cells (PBMCs) and isolated CD3⁺ T-cells from young and old adults (n = 10), following ex vivo stimulation with mock (unstimulated) or live influenza virus. Our genes included: β-actin (ACTB), glyercaldehyde-3-phostphate dehydrogenase (GAPDH), ribosomal protein L13a (RPL13a), ribosomal protein S18 (RPS18), succinate dehydrogenase complex flavoprotein subunit A (SDHA), and ubiquitin-conjugating enzyme E2D2 (UBE2D2).

Results

Reference gene expression varied significantly depending on cell type and stimulation conditions, but not age. Using the comparative ΔCt method, and the previously published software BestKeeper, NormFinder, and geNorm, we show that in PBMCs and T-cells, UBE2D2 and RPS18 were the most stable reference genes, followed by ACTB; however, the expression of UBE2D2 and RPS18 was found to increase with viral stimulation in isolated T-cells, while ACTB expression did not change significantly. No age-related differences in stability were observed for any gene

Conclusions

This study suggests the use of a combination of UBE2D2, RPS18, and ACTB for the study of influenza responses in PBMCs and T-cells, although ACTB alone may be the most optimal choice if choosing to compare target gene expression before and after viral stimulation. Both GAPDH and RPL13a were found to be poor reference genes and should be avoided for studies of this nature.

Postmenopausal Osteoporosis reference genes for qPCR expression assays

Osteoporosis (OP) is a multifactorial disease influenced by genetic factors in more than half of the cases. In spite of the efforts to clarify the relationship among genetic factors and susceptibility to develop OP, many genetic associations need to be further functionally validated. Besides, some limitations as the choice of stably expressed reference genes (RG) should be overcome to ensure the quality and reproducibility of gene expression assays. To our knowledge, a validation study for RG in OP is still missing. We compared the expression levels, using polymerase chain reaction quantitative real time (qPCR) of 10 RG (G6PD, B2M, GUSB, HSP90, EF1A, RPLP0, GAPDH, ACTB, 18 S and HPRT1) to assess their suitability in OP analysis by using GeNorm, Normfinder, BestKeeper and RefFinder programs. A minimal number of two RG was recommended by GeNorm to obtain a reliable normalization. RPLP0 and B2M were identified as the most stable genes in OP studies while ACTB, 18 S and HPRT1 were inadequate for normalization in our data set. Moreover, we showed the dramatic effects of suboptimal RG choice on the quantification of a target gene, highlighting the importance in the identification of the most appropriate reference gene to specific diseases. We suggest the use of RPLP0 and B2M as the most stable reference genes while we do not recommend the use of the least stable reference genes HPRT1, 18 S and ACTB in OP expression assays using PBMC as biological source. Additionally, we emphasize the importance of individualized and careful choice in software and reference genes selection.

Non-nutritional sweeteners effects on endothelial vascular function

AIM

Hyperglycemia status induces endothelial dysfunction, although the underlying pathogenic mechanisms are not fully understood. There are several studies connecting sugar/sweetened beverages to the cardiovascular disease. Currently, many sweeteners have been extensively introduced into lifestyle to normalize blood glucose levels without altering the sweet taste. However, there is growing concern for their impact on metabolic health.

METHODS

Human endothelial cells were treated with Glucose, Fructose, Aspartame, Rebaudioside A, Stevioside, or Steviol. Morphological characteristics, in vitro angiogenesis and array gene expression were analyzed.

RESULTS

High-glucose and fructose concentrations significantly decreased cell features such as angiogenic capability. Interestingly, non-caloric sweeteners did not significantly modified all cell characteristics and they did not compromised cell angiogenic ability. Array gene expression analysis revealed that the chemokine fractalkine (CX3CL1) and the enzyme transferase (HPRT1) were always significantly upregulated and downregulated respectively, after glucose and fructose treatments (P > .05), whereas they were non-differentially expressed with all the other sweeteners. Interestingly, both genes are considered as cardiovascular disease risk biomarkers. Specifically, upregulation of CX3CL1/CX3CR1 occurs in the human placenta and serum levels of the ligand are associated with markers of insulin resistance in GDM.

CONCLUSIONS

Differently from glucose and fructose, steviol glycosides do not damage endothelial cells. Prospective preclinical studies and clinical trials are warranted to confirm the long-term safety of such compounds.

Hybrid 18F-FDG-PET/MRI Measurement of Standardized Uptake Value Coupled with Yin Yang 1 Signature in Metastatic Breast Cancer. A Preliminary Study

Purpose: Detection of breast cancer (BC) metastasis at the early stage is important for the assessment of BC progression status. Image analysis represents a valuable tool for the management of oncological patients. Our preliminary study combined imaging parameters from hybrid ¹⁸F-FDG-PET/MRI and the expression level of the transcriptional factor Yin Yang 1 (YY1) for the detection of early metastases. Methods: The study enrolled suspected n = 217 BC patients that underwent ¹⁸F-FDG-PET/MRI scans. The analysis retrospectively included n = 55 subjects. n = 40 were BC patients and n = 15 imaging-negative female individuals were healthy subjects (HS). Standard radiomics parameters were extracted from PET/MRI image. RNA was obtained from peripheral blood mononuclear cells and YY1 expression level was evaluated by real time reverse transcription polymerase chain reactions (qRT-PCR). An enzyme-linked immuosorbent assay (ELISA) was used to determine the amount of YY1 serum protein. Statistical comparison between subgroups was evaluated by Mann-Whitney U and Spearman’s tests. Results: Radiomics showed a significant positive correlation between Greg-level co-occurrence matrix (GLCM) and standardized uptake value maximum (SUVmax) (r = 0.8 and r = 0.8 respectively) in BC patients. YY1 level was significant overexpressed in estrogen receptor (ER)-positive/progesteron receptor-positive/human epidermal growth factor receptor2-negative (ER+/PR+/HER2-) subtype of BC patients with synchronous metastasis (SM) at primary diagnosis compared to metachronous metastasis (MM) and HS (p < 0.001) and correlating significantly with ¹⁸F-FDG-uptake parameter (SUVmax) (r = 0.48). Conclusions: The combination of functional ¹⁸F-FDG-PET/MRI parameters and molecular determination of YY1 could represent a novel integrated approach to predict synchronous metastatic disease with more accuracy than ¹⁸F-FDG-PET/MRI alone.

Correlation of Circulating miR-765, miR-93-5p, and miR-433-3p to Obstructive Coronary Heart Disease Evaluated by Cardiac Computed Tomography

Epigenetic-sensitive mechanisms may be correlated both to pathogenesis and prognosis of coronary heart disease (CHD). We prospectively investigated some plasma circulating microRNA levels in patients undergoing cardiac computed tomography for suspected CHD (n = 95). We show that let-7c-5p, miR-765, miR-483-5p, miR-31-5p, and miR-206 were upregulated in CHD patients (n = 66) versus healthy subjects HS (n = 29); moreover, let-7c-5p, miR-765, miR- 483-5p showed higher expression in obstructive CHD (n = 36) compared to no obstructive CHD patients (n = 66). Remarkably, miR-765, miR-93-5p, and miR-433-3p showed an upregulation in patients with critical coronary stenosis. Multivariate regression analysis demonstrated that miR-765, miR-31-5p, and miR-206 were independently associated with CHD while circulating levels of miR-765 (p = 0.035), miR-433-3p (p = 0.043), and miR-93-5p (p = 0.041) were significantly higher in critical stenosis patients. Receiver operating characteristic curve analysis revealed a good performance for miR-765, miR-93-5p, and miR-433-3p on predicting CHD severity. In conclusion, our study represents a combined epigenetic/imaging approach useful to support the diagnosis and prediction of CHD.

Variations in contents of hyaluronan in the peritumoral micro-environment of human chondrosarcoma

A concept consolidated in recent years is that tumor growth depends to a great extent on the micro-environment surrounding the tumor, which has a fundamental role in tumor progression and in determining the effectiveness of therapies. Our analysis focuses on chondrosarcoma, the second primary malignant bone tumor, resistant to both chemotherapeutic and radiation therapy. We quantified hyaluronan, one of the main components of the extracellular matrix, with the aim of comparing its amount in the connective tissue surrounding the tumor with intra-tumoral tissue and healthy fascia of the same anatomic district, viewed as a health control. We demonstrate that hyaluronan increased significantly in the peritumoral stroma compared with the healthy fascia, which showed an average amount according to the physical characteristics of body districts by a mean value of 26.9 μg/g. In the peritumoral stroma, the mean hyaluronan content reached 132.6 μg/g (mean value of 63.2 μg/g). The p-value was less than 0.01, showing a highly significant statistical difference. Surprisingly, no significant differences were detected as a function of age, gender, or tumor grade. The levels of hyaluronan were comparable in peritumoral and tumor tissues, although there were differences depending on the state of necrosis. In addition, data on the expression of hyaluronic acid synthetase showed a decrease of about 50% in peritumoral and tumor tissues, indicating alterations in hyaluronan turnover and synthesis. This work demonstrates a variation in hyaluronan contents around the chondrosarcoma, likely correlated with the aggressiveness and resistance to chemotherapy of this tumors. Statement of Clinical Significance: Deeper knowledge about the composition of the peritumoral stroma, rich in extracellular matrix, will enhance better study and understanding of the metastatic potential of tumors and their prognostic indices. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:503-509, 2019.

Evidence of association of circulating epigenetic-sensitive biomarkers with suspected coronary heart disease evaluated by Cardiac Computed Tomography

Circulating biomarkers available in clinical practice do not allow to stratify patients with coronary heart disease (CHD) prior the onset of a clinically relevant event. We evaluated the methylation status of specific genomic segments and gene expression in peripheral blood of patients undergoing Cardiac Computed Tomography (CCT) for CHD (n = 95).

We choose to investigate cholesterol metabolism. Methylation and gene expression of low density lipoprotein receptor (LDLR), sterol regulatory element-binding factor 2 (SREBF2) and ATP-binding cassette transporter 1 (ABCA1) were evaluated by qRT-PCR. Calcium score (CACS), stenosis degree, total plaque volume (TPV), calcified plaque volume (CPV), non-calcified plaque volume (NCPV) and plaque burden (PB) were assessed in all CHD patients (n = 65). The percentage of methylation at the specific analyzed segment of LDLR promoter was higher in CHD patients vs healthy subjects (HS) (n = 30) (p = 0.001). LDLR, SREBF2 and ABCA1 mRNAs were up-regulated in CHD patients vs HS (p = 0.02; p = 0.019; p = 0.008). SREBF2 was overexpressed in patients with coronary stenosis ≥50% vs subjects with stenosis <50% (p = 0.036). After adjustment for risk factors and clinical features, ABCA1 (p = 0.005) and SREBF2 (p = 0.010) gene expression were identified as independent predictors of CHD and severity. ROC curve analysis revealed a good performance of ABCA1 on predicting CHD (AUC = 0.768; p<0.001) and of SREBF2 for the prediction of disease severity (AUC = 0.815; p<0.001). Moreover, adjusted multivariate analysis demonstrated SREBF2 as independent predictor of CPV, NCPV and TPV (p = 0.022; p = 0.002 and p = 0.006) and ABCA1 as independent predictor of NCPV and TPV (p = 0.002 and p = 0.013).

CHD presence and characteristics are related to selected circulating transcriptional and epigenetic-sensitive biomarkers linked to cholesterol pathway. More extensive analysis of CHD phenotypes and circulating biomarkers might improve and personalize cardiovascular risk stratification in the clinical settings.

PR/SET Domain Family and Cancer: Novel Insights from the Cancer Genome Atlas

The PR/SET domain gene family (PRDM) encodes 19 different transcription factors that share a subtype of the SET domain [Su(var)3-9, enhancer-of-zeste and trithorax] known as the PRDF1-RIZ (PR) homology domain. This domain, with its potential methyltransferase activity, is followed by a variable number of zinc-finger motifs, which likely mediate protein⁻protein, protein⁻RNA, or protein⁻DNA interactions. Intriguingly, almost all PRDM family members express different isoforms, which likely play opposite roles in oncogenesis. Remarkably, several studies have described alterations in most of the family members in malignancies. Here, to obtain a pan-cancer overview of the genomic and transcriptomic alterations of PRDM genes, we reanalyzed the Exome- and RNA-Seq public datasets available at The Cancer Genome Atlas portal. Overall, PRDM2, PRDM3/MECOM, PRDM9, PRDM16 and ZFPM2/FOG2 were the most mutated genes with pan-cancer frequencies of protein-affecting mutations higher than 1%. Moreover, we observed heterogeneity in the mutation frequencies of these genes across tumors, with cancer types also reaching a value of about 20% of mutated samples for a specific PRDM gene. Of note, ZFPM1/FOG1 mutations occurred in 50% of adrenocortical carcinoma patients and were localized in a hotspot region. These findings, together with OncodriveCLUST results, suggest it could be putatively considered a cancer driver gene in this malignancy. Finally, transcriptome analysis from RNA-Seq data of paired samples revealed that transcription of PRDMs was significantly altered in several tumors. Specifically, PRDM12 and PRDM13 were largely overexpressed in many cancers whereas PRDM16 and ZFPM2/FOG2 were often downregulated. Some of these findings were also confirmed by real-time-PCR on primary tumors.

Splicing regulators in endothelial cell differentiation

AIMS

Alternative splicing represents a key mechanism of gene regulation. Despite its role in regulating cell pluripotency and differentiation being well known, the underlining mechanisms are still poorly studied. Here, we investigated the possible involvement of splicing regulators during the different steps of endothelial cell differentiation through expression studies on human circulating progenitors.

METHODS

Total RNAs were extracted from all cells and reverse-transcribed. Semiquantitative and real-time RT-PCR was performed using selective oligonucleotides. Differences between group means were considered significant at P value less than 0.05 and more significant at P value less than 0.01. Protein extracts were incubated with an antibody directed against MED23. Immunoprecipitation of supernatants and pellets was probed with both anti-Muscleblind-like splicing regulator (MBNL)1 and anti-MBNL2 antibodies.

RESULTS

Several clinical trials demonstrated the safety and efficacy of progenitor cells in regenerative therapy of the cardiovascular system. Particularly, we analyzed the expression of genes belonging to muscleblind family members and MED complex subunits, which are known to be involved during differentiation in other models. This study shows that MED23, MBNL1 and MBNL2 were all expressed at high levels only in differentiated cells. Moreover, immunoprecipitation assays indicated that MED23 is able to bind MBNLs in endothelial cells.

CONCLUSION

Our data suggest that MED23, MBNL1 and MBNL2 could regulate alternative splicing events activated during differentiation through a common mechanism. Hence, these observations corroborate previous evidence that splicing regulators may have an essential role in the basic apparatus required for cell pluripotency and reprogramming, allowing identification of novel biomarkers to use for early diagnosis in cardiovascular diseases.

Identification of suitable reference genes for real-time quantitative PCR analysis of hydrogen peroxide-treated human umbilical vein endothelial cells

Background

Oxidative stress can induce cell injury in vascular endothelial cells, which is the initial event in the development of atherosclerosis. Although quantitative real-time polymerase chain reaction (qRT-PCR) has been widely used in gene expression studies in oxidative stress injuries, using carefully validated reference genes has not received sufficient attention in related studies. The objective of this study, therefore, was to select a set of stably expressed reference genes for use in qRT-PCR normalization in oxidative stress injuries in human umbilical vein endothelial cells (HUVECs) induced by hydrogen peroxide (H₂O₂).

Results

Using geNorm analysis, we found that five stably expressed reference genes were sufficient for normalization in qRT-PCR analysis in HUVECs treated with H₂O₂. Genes with the most stable expression according to geNorm were U6, TFRC, RPLP0, GAPDH, and ACTB, and according to NormFinder were ALAS1, TFRC, U6, GAPDH, and ACTB.

Conclusion

Taken together, our study demonstrated that the expression stability of reference genes may differ according to the statistical program used. U6, TFRC, RPLP0, GAPDH, and ACTB was the optimal set of reference genes for studies on gene expression performed by qRT-PCR assays in HUVECs under oxidative stress study.

Characterizing biomarkers in osteosarcoma metastasis based on an ego-network

OBJECTIVES

To characterize biomarkers that underlie osteosarcoma (OS) metastasis based on an ego-network.

RESULTS

From the microarray data, we obtained 13,326 genes. By combining PPI data and microarray data, 10,520 shared genes were found and constructed into ego-networks. 17 significant ego-networks were identified with p < 0.05. In the pathway enrichment analysis, seven ego-networks were identified with the most significant pathway.

CONCLUSIONS

These significant ego-modules were potential biomarkers that reveal the potential mechanisms in OS metastasis, which may contribute to understanding cancer prognoses and providing new perspectives in the treatment of cancer.

In Vivo and In Vitro Analysis in Coronary Artery Disease Related to Type 2 Diabetes

AIM

The leading cause of morbidity and mortality in patients with type 2 diabetes mellitus (DM) is coronary artery disease (CAD), a condition often asymptomatic but severe in these patients. Although glucose metabolism impairment and oxidative stress are known actors in the endothelial dysfunction/remodeling that occurs in diabetic patients, the relationship between cardiovascular disorders and DM is not fully understood. We have performed both an in vivo imaging and in vitro molecular analysis to investigate diabetic-specific CAD alterations.

METHODS

Computed tomography coronary angiography (CTCA) was performed in a group of 20 diabetic patients with CAD (DM⁺CAD⁺), 20 non-diabetic with CAD (DM^-CAD⁺), 10 diabetic non-CAD patients (DM⁺CAD^-), and 20 non-diabetic healthy subjects (HS). Imaging quantitative parameters such as calcium score (Cascore), calcified plaque volume (CPV), non-calcified plaque volume (NCPV), total plaque volume (TPV), remodeling index (RI), and plaque burden were extracted for each CAD subject. Moreover, the expression levels of superoxide dismutase 2 (SOD2) and liver X receptor alpha (LXRα) genes were analyzed in the peripheral blood mononuclear cells, whereas hyaluronan (HA) concentrations were evaluated in the plasma of each subject.

RESULTS

Imaging parameters, such as Cascore, CPV, RI, and plaque burden, were significantly higher in DM⁺CAD⁺ group, compared to DM^-CAD⁺ (P = 0.019; P = 0.014; P < 0.001, P < 0.001, respectively). SOD2 mRNA was downregulated, while LXRα gene expression was upregulated in DM⁺CAD^-, DM⁺CAD⁺, and DM^-CAD⁺ groups compared to HS (P = 0.001, P = 0.03, and P = 0.001 for SOD2 and P = 0.006, P = 0.008, and P < 0.001 for LXRα, respectively). Plasmatic levels of HA were higher in DM^-CAD⁺, DM⁺CAD^-, and DM⁺CAD⁺ groups, compared to HS (P = 0.001 for the three groups). When compared to DM^-CAD⁺, HA concentration was higher in DM⁺CAD^- (P = 0.008) and DM⁺CAD⁺ (P < 0.001) with a significant difference between the two diabetic groups (P = 0.003). Moreover, HA showed a significant association with diabetes (P = 0.01) in the study population, and the correlation between HA levels and glycemia was statistically significant (ρ = 0.73, P < 0.001).

CONCLUSION

In our population, imaging parameters highlight a greater severity of CAD in diabetic patients. Among molecular parameters, HA is modulated by diabetic CAD-related alterations while SOD2 and LXRα are found to be more associated with CAD but do not discriminate between diabetic and non-diabetic subgroups.

Hormone receptor expression in human fascial tissue

Many epidemiologic, clinical, and experimental findings point to sex differences in myofascial pain in view of the fact that adult women tend to have more myofascial problems with respect to men. It is possible that one of the stimuli to sensitization of fascial nociceptors could come from hormonal factors such as estrogen and relaxin, that are involved in extracellular matrix and collagen remodeling and thus contribute to functions of myofascial tissue. Immunohistochemical and molecular investigations (real-time PCR analysis) of relaxin receptor 1 (RXFP1) and estrogen receptor-alpha (ERα) localization were carried out on samples of human fascia collected from 8 volunteers patients during orthopedic surgery (all females, between 42 and 70 yrs, divided into pre- and post-menopausal groups), and in fibroblasts isolated from deep fascia, to examine both protein and RNA expression levels. We can assume that the two sex hormone receptors analyzed are expressed in all the human fascial districts examined and in fascial fibroblasts culture cells, to a lesser degree in the post-menopausal with respect to the pre-menopausal women. Hormone receptor expression was concentrated in the fibroblasts, and RXFP1 was also evident in blood vessels and nerves. Our results are the first demonstrating that the fibroblasts located within different districts of the muscular fasciae express sex hormone receptors and can help to explain the link between hormonal factors and myofascial pain. It is known, in fact, that estrogen and relaxin play a key role in extracellular matrix remodeling by inhibiting fibrosis and inflammatory activities, both important factors affecting fascial stiffness and sensitization of fascial nociceptors.

Mining for genes related to choroidal neovascularization based on the shortest path algorithm and protein interaction information

BACKGROUND

Choroidal neovascularization (CNV) is a serious eye disease that may cause visual loss, especially for older people. Many factors have been proven to induce this disease including age, gender, obesity, and so on. However, until now, we have had limited knowledge on CNV's pathogenic mechanism. Discovering the genes that underlie this disease and performing extensive studies on them can help us to understand how CNV occurs and design effective treatments.

METHODS

In this study, we designed a computational method to identify novel CNV-related genes in a large protein network constructed using the protein-protein interaction information in STRING. The candidate genes were first extracted from the shortest paths connecting any two known CNV-related genes and then filtered by a permutation test and using knowledge of their linkages to known CNV-related genes.

RESULTS

A list of putative CNV-related candidate genes was accessed by our method. These genes are deemed to have strong relationships with CNV.

CONCLUSIONS

Extensive analyses of several of the putative genes such as ANK1, ITGA4, CD44 and others indicate that they are related to specific biological processes involved in CNV, implying they may be novel CNV-related genes.

GENERAL SIGNIFICANCE

The newfound putative CNV-related genes may provide new insights into CNV and help design more effective treatments. This article is part of a Special Issue entitled "System Genetics" Guest Editor: Dr. Yudong Cai and Dr. Tao Huang.

Selection of reliable reference genes in eutopic and ectopic endometrium for quantitative expression studies

PURPOSE

Physiological changes during menstrual cycle cause the endometrium and endometriosis to develop specific kind of tissues, especially in regard to the gene expression profiles, which may include also housekeeping genes, commonly used as reference genes (RGs) in quantitative studies. Reverse transcription, followed by quantitative polymerase chain reaction (RT-qPCR) is the most precise and commonly used method in gene expression studies. In order to reduce effects of technical approaches and biological variability of gene's expression level, the studies often employ RGs in experimental data normalization. However, the expression of RGs is not always stable and depends on several variables. Thus, the selection of appropriate RG is one of the most significant steps to obtain reliable results in RT-qPCR-based methods.

MATERIAL AND METHODS

With the usage of RT-qPCR, we researched the expression of seven genes (ACTB, B2M, G6PD, GAPD, GUSB, HPRT and PPIA) as reliable reference genes in eutopic and ectopic endometrial tissue specimens obtained during standard surgery of women of reproductive age. Stability of expression level was analyzed by the most universal MS Excel plug-ins including: geNorm, NormFinder and BestKeeper. The descriptive statistics were evaluated using Statistica software.

RESULTS

The distribution of threshold (Ct) values was not equal. We identified genes with higher expression level (referring to Ct values) such as ACTB and B2M, medium e.g., GAPD and low expression level, e.g., G6PD and HPRT. We demonstrated that the stability of the analyzed reference genes was not homogenous, and different algorithms pointed to PPIA, GAPD and B2M as the most stable ones in eutopic and ectopic endometrium. On the contrary to these, GUSB and G6PD were the most unstable ones.

CONCLUSIONS

In RT-qPCR-based analyses of gene expression level in eutopic and ectopic endometrium, we strongly recommend that a minimum of two reference genes are to be used and we determined that the most suitable seem to be PPIA and GAPD.

Doubling Throughput of a Real-Time PCR

The invention of polymerase chain reaction (PCR) in 1983 revolutionized many areas of science, due to its ability to multiply a number of copies of DNA sequences (known as amplicons). Here we report on a method to double the throughput of quantitative PCR which could be especially useful for PCR-based mass screening. We concurrently amplified two target genes using only single fluorescent dye. A FAM probe labelled olionucleotide was attached to a quencher for one amplicon while the second one was without a probe. The PCR was performed in the presence of the intercalating dye SYBR Green I. We collected the fluorescence amplitude at two points per PCR cycle, at the denaturation and extension steps. The signal at denaturation is related only to the amplicon with the FAM probe while the amplitude at the extension contained information from both amplicons. We thus detected two genes within the same well using a single fluorescent channel. Any commercial real-time PCR systems can use this method doubling the number of detected genes. The method can be used for absolute quantification of DNA using a known concentration of housekeeping gene at one fluorescent channel.

Community-Based HIV-1 Early Diagnosis and Risk Behavior Analysis of Men Having Sex with Men in Hong Kong

The increasing prevalence of HIV-1 among men having sex with men (MSM) calls for an investigation of HIV-1 prevalence and incidence in MSM by early diagnosis to assist with early preventive interventions in Hong Kong. The participants were recruited randomly from MSM communities within a one-year period. Rapid HIV Test (RHT) and real-time dried blood spot (DBS)-based quantitative polymerase chain reaction (DBS-qPCR) were used for the early diagnosis of 474 participants. Risk behavior analysis was performed by studying information obtained from the participants during the study period. The HIV-1 prevalence and incident rates in the studied MSM population were 4.01% (19/474) and 1.47% (7/474), respectively. Three infected participants were found at the acute phase of infection by DBS-qPCR. Only 46.4% (220/474) MSM were using condoms regularly for anal sex. HIV infection significantly correlated with unprotected receptive anal sex and syphilis infection. An increased number of infections was found among foreign MSM in Hong Kong. This study is the first to use DBS-qPCR to identify acutely infected individuals in a community setting and to provide both the prevalence and incident rates of HIV-1 infection among MSM in Hong Kong. The risk analysis provided evidence that behavior intervention strengthening is necessary to fight against the increasing HIV-1 epidemic among MSM in Hong Kong and surrounding regions in Asia.

Evaluation and application of the strand-specific protocol for next-generation sequencing

Next-generation sequencing (NGS) has become a powerful sequencing tool, applied in a wide range of biological studies. However, the traditional sample preparation protocol for NGS is non-strand-specific (NSS), leading to biased estimates of expression for transcripts overlapped at the antisense strand. Strand-specific (SS) protocols have recently been developed. In this study, we prepared the same RNA sample by using the SS and NSS protocols, followed by sequencing with Illumina HiSeq platform. Using real-time quantitative PCR as a standard, we first proved that the SS protocol more precisely estimates gene expressions compared with the NSS protocol, particularly for those overlapped at the antisense strand. In addition, we also showed that the sequence reads from the SS protocol are comparable with those from conventional NSS protocols in many aspects. Finally, we also mapped a fraction of sequence reads back to the antisense strand of the known genes, originally without annotated genes located. Using sequence assembly and PCR validation, we succeeded in identifying and characterizing the novel antisense genes. Our results show that the SS protocol performs more accurately than the traditional NSS protocol and can be applied in future studies.

Steroid responsive regulation of IFNγ2 alternative splicing and its possible role in germ cell proliferation in medaka

Interferon gamma (IFNγ) is an active player in estrogen dependent immuno-regulation of fish. The present work was aimed to characterize the alternatively spliced isoforms of IFNγ2 in the gonadal sex development in medaka. Phylogenetic analysis demonstrated that IFNγ2a and 2b were clustered with fish specific interferon gamma. Our in vitro promoter and mini-genome analysis data confirmed that alternative splicing of IFNγ2 is regulated by estrogens and androgens. Tissue distribution, quantitative PCR and ISH data demonstrated ubiquitous expression of IFNγ2a, while IFNγ2b was only expressed predominantly in female germ cells than males. This was further confirmed by germ cell specific GFP signals in the IFNγ2b-GFP over-expressed embryos and specific induction of IFNγ2b expression in the BrdU positive cells. All together our data suggest that steroid responsive alternatively spliced IFNγ2b isoforms might have some indirect roles in germ cell proliferation and thus can be an important candidate for immuno-reproductive interaction studies.

Surgical procedures and postsurgical tissue processing significantly affect expression of genes and EGFR-pathway proteins in colorectal cancer tissue

An understanding of tissue data variability in relation to processing techniques during and postsurgery would be desirable when testing surgical specimens for clinical diagnostics, drug development, or identification of predictive biomarkers. Specimens of normal and colorectal cancer (CRC) tissues removed during colon and liver resection surgery were obtained at the beginning of surgery and postsurgically, tissue was fixed at 10, 20, and 45 minutes. Specimens were analyzed from 50 patients with primary CRC and 43 with intrahepatic metastasis of CRC using a whole genome gene expression array. Additionally, we focused on the epidermal growth factor receptor pathway and quantified proteins and their phosphorylation status in relation to tissue processing timepoints. Gene and protein expression data obtained from colorectal and liver specimens were influenced by tissue handling during surgery and by postsurgical processing time. To obtain reliable expression data, tissue processing for research and diagnostic purposes needs to be highly standardized.

RNA-Seq for the identification of novel Mediator transcripts in endothelial progenitor cells

Mediator (MED) complex is a multiprotein playing a key role in the eukaryotic transcription. Alteration of MED function may have enormous pathophysiological consequences and several MED genes have been implicated in human diseases. Here, we have combined computational and experimental approaches to identify and characterize, new transcripts generated by alternative splicing (AS) for all MED genes, through the analysis of our recently published RNA-Sequencing datasets of endothelial progenitor cells (EPCs). This combined strategy allowed us to identify novel transcripts for MED4, MED9, MED11, MED14, MED27 and CDK8 most of them generated by AS. All the newly identified transcripts, except MED11, are predicted to encode novel protein isoforms. The identification of novel MED variants could lead to the finding of other MED complexes with different functions depending on their subunit composition. Finally, the expression profile of all MED genes, together with an extensive gene expression analysis, may be useful to better classify the diverse subsets of cell populations that contribute to neovascularization.

Identification of stable reference genes for gene expression analysis of three-dimensional cultivated human bone marrow-derived mesenchymal stromal cells for bone tissue engineering

The principles of tissue engineering (TE) are widely used for bone regeneration concepts. Three-dimensional (3D) cultivation of autologous human mesenchymal stromal cells (MSCs) on porous scaffolds is the basic prerequisite to generate newly formed bone tissue. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is a specific and sensitive analytical tool for the measurement of mRNA-levels in cells or tissues. For an accurate quantification of gene expression levels, stably expressed reference genes (RGs) are essential to obtain reliable results. Since the 3D environment can affect a cell's morphology, proliferation, and gene expression profile compared with two-dimensional (2D) cultivation, there is a need to identify robust RGs for the quantification of gene expression. So far, this issue has not been adequately investigated. The aim of this study was to identify the most stably expressed RGs for gene expression analysis of 3D-cultivated human bone marrow-derived MSCs (BM-MSCs). For this, we analyzed the gene expression levels of n=31 RGs in 3D-cultivated human BM-MSCs from six different donors compared with conventional 2D cultivation using qRT-PCR. MSCs isolated from bone marrow aspirates were cultivated on human cancellous bone cube scaffolds for 14 days. Osteogenic differentiation was assessed by cell-specific alkaline phosphatase (ALP) activity and expression of osteogenic marker genes. Expression levels of potential reference and target genes were quantified using commercially available TaqMan(®) assays. mRNA expression stability of RGs was determined by calculating the coefficient of variation (CV) and using the algorithms of geNorm and NormFinder. Using both algorithms, we identified TATA box binding protein (TBP), transferrin receptor (p90, CD71) (TFRC), and hypoxanthine phosphoribosyltransferase 1 (HPRT1) as the most stably expressed RGs in 3D-cultivated BM-MSCs. Notably, genes that are routinely used as RGs, for example, beta actin (ACTB) and ribosomal protein L37a (RPL37A), were among the least stable genes. We recommend the combined use of TBP, TFRC, and HPRT1 for the accurate and robust normalization of qRT-PCR data of 3D-cultivated human BM-MSCs.

Identification of reference genes for qRT-PCR in human lung squamous-cell carcinoma by RNA-Seq

Although the accuracy of quantitative real-time polymerase chain reaction (qRT-PCR) is highly dependent on the reliable reference genes, many commonly used reference genes are not stably expressed and as such are not suitable for quantification and normalization of qRT-PCR data. The aim of this study was to identify novel reliable reference genes in lung squamous-cell carcinoma. We used RNA sequencing (RNA-Seq) to survey the whole genome expression in 5 lung normal samples and 44 lung squamous-cell carcinoma samples. We evaluated the expression profiles of 15 commonly used reference genes and identified five additional candidate reference genes. To validate the RNA-Seq dataset, we used qRT-PCR to verify the expression levels of these 20 genes in a separate set of 100 pairs of normal lung tissue and lung squamous-cell carcinoma samples, and then analyzed these results using geNorm and NormFinder. With respect to 14 of the 15 common reference genes (B2M, GAPDH, GUSB, HMBS, HPRT1, IPO8, PGK1, POLR2A, PPIA, RPLP0, TBP, TFRC, UBC, and YWHAZ), the expression levels were either too low to be easily detected, or exhibited a high degree of variability either between lung normal and squamous-cell carcinoma samples, or even among samples of the same tissue type. In contrast, 1 of the 15 common reference genes (ACTB) and the 5 additional candidate reference genes (EEF1A1, FAU, RPS9, RPS11, and RPS14) were stably and constitutively expressed at high levels in all the samples tested. ACTB, EEF1A1, FAU, RPS9, RPS11, and RPS14 are ideal reference genes for qRT-PCR analysis of lung squamous-cell carcinoma, while 14 commonly used qRT-PCR reference genes are less appropriate in this context.

Gene expression profile of the whole Mediator complex in human osteosarcoma and normal osteoblasts

Mediator complex (MED) is an essential multi-subunit component of the transcription apparatus and plays a key role in the transcription regulation of many genes involved in several diseases, including cancer. Recently, numerous MED subunits have been implicated in cancer development and metastasis, and specific alterations in their coding genes have been found to correlate with some malignancies. It is conceivable that a specific MED alteration pattern can characterize each cancer type. However, to date, no study has reported the complete picture of MED subunits in a specific tumor. Thus, the aim of this study was to investigate for the first time the gene expression profile of the whole MED complex in human osteosarcoma (OS). To this purpose, we have examined all the MED subunit genes in three OS cell lines compared to normal osteoblasts by real-time RT-PCR. Interestingly, our findings indicate that the expression of most of the MED genes is altered in OS. Moreover, a very high overexpression of MED20 and MED31 can be observed in all the analyzed OS cells, thus suggesting for the first time a potential role of these subunits in human malignancies. Overall, this study may open the way to other functional studies exploring the role of the whole complex in cancer development and progression. These findings may lead to the identification of novel biomarkers, which can be used also in combination with imaging techniques for early detection, and/or to develop novel targets for innovative therapeutic approaches.

Selecting housekeeping genes as references for the normalization of quantitative PCR data in breast cancer

OBJECTIVE

The common reference genes of choice in relative gene expression studies based on quantitative real time polymerase chain reaction, ACTB and B2M, were shown to be regulated differently in respect to tissue type. In this study, the stability of the selected housekeeping genes for normalizing the qPCR data were identified in the tumor and its adjacent tissues in invasive breast cancer, and the variability of their levels according to the stages and the histopathologic subtypes was analyzed.

METHODS

Four housekeeping genes: PUM1, RPL13A, B2M, and ACTB were analyzed in 99 surgically excised tissue specimens (50 tumor, 45 tumor adjacent and 4 normal breast tissues). Three of the most common softwares (GeNorm, NormFinder, and BestKeeper) were used for calculation purposes.

RESULTS

When all of the tissue samples were included in analyses, PUM1 was the most stable gene according to calculations made with both NormFinder and BestKeeper; while PUM1/RPL13A combination was the most stable by GeNorm software. The PUM1 gene was also identified as the most stable gene among the four in all sample groups (in both Estrogen Receptor positive and Estrogen Receptor negative subgroups of invasive breast carcinoma and in normal breast tissue) according to calculations made using the NormFinder software.

CONCLUSION

While suggesting PUM1 is one of the most stable single gene and the PUM1/RPL13A pair as one of the best housekeeping genes for the normalization of expression studies in invasive breast tumor studies, it will be more practical to evaluate stability once more and decide upon the reference gene accordingly within the sample group itself.

Identification of optimal reference genes for RT-qPCR in the rat hypothalamus and intestine for the study of obesity

BACKGROUND

Obesity has a complicated metabolic pathology, and defining the underlying mechanisms of obesity requires integrative studies with molecular end points. Real-time quantitative PCR (RT-qPCR) is a powerful tool that has been widely utilized. However, the importance of using carefully validated reference genes in RT-qPCR seems to have been overlooked in obesity-related research. The objective of this study was to select a set of reference genes with stable expressions to be used for RT-qPCR normalization in rats under fasted vs re-fed and chow vs high-fat diet (HFD) conditions.

DESIGN

Male long-Evans rats were treated under four conditions: chow/fasted, chow/re-fed, HFD/fasted and HFD/re-fed. Expression stabilities of 13 candidate reference genes were evaluated in the rat hypothalamus, duodenum, jejunum and ileum using the ReFinder software program. The optimal number of reference genes needed for RT-qPCR analyses was determined using geNorm.

RESULTS

Using geNorm analysis, we found that it was sufficient to use the two most stably expressed genes as references in RT-qPCR analyses for each tissue under specific experimental conditions. B2M and RPLP0 in the hypothalamus, RPS18 and HMBS in the duodenum, RPLP2 and RPLP0 in the jejunum and RPS18 and YWHAZ in the ileum were the most suitable pairs for a normalization study when the four aforementioned experimental conditions were considered.

CONCLUSIONS

Our study demonstrates that gene expression levels of reference genes commonly used in obesity-related studies, such as ACTB or RPS18, are altered by changes in acute or chronic energy status. These findings underline the importance of using reference genes that are stable in expression across experimental conditions when studying the rat hypothalamus and intestine, because these tissues have an integral role in the regulation of energy homeostasis. It is our hope that this study will raise awareness among obesity researchers on the essential need for reference gene validation in gene expression studies.