Validation of Reference Genes for Normalization of Relative qRT-PCR Studies in Papillary Thyroid Carcinoma

Selection and Validation of Appropriate Reference Genes for qRT-PCR Analysis of Iris germanica L. Under Various Abiotic Stresses

Choosing the appropriate reference genes for quantitative real-time PCR (qRT-PCR) is very important for accurately evaluating expression of target genes. Iris germanica L. is a widely used horticultural plant with high ornamental value, which also shows a strong ability to tolerate abiotic stresses. No comprehensive research has been carried out on optimal reference genes in Iris germanica L. under abiotic stress. In this study, nine reference genes were selected as candidates based on the transcriptome sequencing data of Iris germanica L. The assessment of expression stability under various abiotic stress was conducted using four distinct methods (GeNorm, NormFinder, BestKeeper, and RefFinder). It was found that the optimal reference genes were ACT and F3H for drought and different temperature stresses. EF1α and ACT exhibited superior performance under salt stress. The expression of the IgP5CS gene was evaluated to provide additional validation for the accuracy of the selected optimal reference genes, indicating that inappropriate may lead to significant deviations in the results. This research identified reliable reference genes in I. germanica L. across various abiotic stress conditions, thereby facilitating the investigation into the molecular mechanisms responsible for stress tolerance in I. germanica L.

Expression of TSPAN1 and its link to thyroid nodules: one step forward on the path to thyroid tumorigenesis biomarkers

BACKGROUND

Thyroid cancer is ranked as the most common malignancy within the endocrine system and the seventh most prevalent cancer in women globally. Thyroid malignancies require evaluating biomarkers capable of distinguishing between them for accurate diagnosis. We examined both mRNA and protein levels of TSPAN1 in plasma and tissue samples from individuals with thyroid nodules to aid this endeavour.

METHODS

In this case-control study, TSPAN1 was assessed at both protein and mRNA levels in 90 subjects, including papillary thyroid cancer (PTC; N = 60), benign (N = 30), and healthy subjects (N = 26) using enzyme-linked immunosorbent assay (ELISA) and SYBR-Green Real-Time PCR, respectively.

RESULTS

TSPAN1 plasma levels were decreased in PTC and benign compared to healthy subjects (P = 0.002). TSPAN1 mRNA levels were also decremented in the tumoral compared to the paired normal tissues (P = 0.012); this drop was also observed in PTC patients compared to benign patients (P = 0.001). Further, TSPAN1 had an appropriate diagnostic value for detecting PTC patients from healthy plasma samples with a sensitivity of 76.7% and specificity of 65.4% at the cutoff value < 2.7 (ng/ml).

CONCLUSION

TSPAN1 levels are significantly reduced in patients with benign and PTC, demonstrating its potential value as a diagnostic biomarker. Additionally, the significant reduction in TSPAN1 mRNA expression within PTC tumor tissues may suggest its involvement in tumor progression and development. Further studies, including larger-scale validation studies and mechanistic investigations, are imperative to clarify the molecular mechanisms behind TSPAN1 and, ultimately, its clinical utility for treating thyroid disorders.

Screening and validating the optimal panel of housekeeping genes for 4T1 breast carcinoma and metastasis studies in mice

The 4T1 model is extensively employed in murine studies to elucidate the mechanisms underlying the carcinogenesis of triple-negative breast cancer. Molecular biology serves as a cornerstone in these investigations. However, accurate gene expression analyses necessitate data normalization employing housekeeping genes (HKGs) to avert spurious results. Here, we initially delve into the characteristics of the tumor evolution induced by 4T1 in mice, underscoring the imperative for additional tools for tumor monitoring and assessment methods for tracking the animals, thereby facilitating prospective studies employing this methodology. Subsequently, leveraging various software platforms, we assessed ten distinct HKGs (GAPDH, 18 S, ACTB, HPRT1, B2M, GUSB, PGK1, CCSER2, SYMPK, ANKRD17) not hitherto evaluated in the 4T1 breast cancer model, across tumors and diverse tissues afflicted by metastasis. Our principal findings underscore GAPDH as the optimal HKG for gene expression analyses in tumors, while HPRT1 emerged as the most stable in the liver and CCSER2 in the lung. These genes demonstrated consistent expression and minimal variation among experimental groups. Furthermore, employing these HKGs for normalization, we assessed TNF-α and VEGF expression in tissues and discerned significant disparities among groups. We posit that this constitutes the inaugural delineation of an ideal HKG for experiments utilizing the 4T1 model, particularly in vivo settings.

Evaluation of reference genes suitable for studying mRNAs and microRNAs expression in thyroid neoplasms

Analysis of gene expression is a pivotal method at the core of biomarkers studies and cancer research. Currently, RT-qPCR is the most commonly used technique to investigate the expression of certain genes. The accurate and reliable result relies on an effective normalization step using suitable reference genes. The present study was designed to evaluate the eligibility of a set of candidate mRNAs and snoRNA as reference genes in the most common human thyroid neoplasms. We tested the expression levels of eleven mRNA and small RNA housekeeping genes in thyroid samples. The stability of the candidate genes was examined in different thyroid lesions and under different experimental conditions. Results were compared to the reported data in the TCGA database. Our results suggested HPRT1 and ACTB as the best mRNA reference genes, SNORD96A, and SNORD95 as the best miRNA reference genes in thyroid tissues. These genes showed the most stable expression pattern among different thyroid lesions as well as different experimental conditions. The findings in this study highlight the effect of reference genes selection on data interpretation and emphasize the importance of testing for suitable reference genes to be used in specific types of cells and experimental conditions to ensure the validity and accuracy of results.

Exploring the potential of myo-inositol in thyroid disease management: focus on thyroid cancer diagnosis and therapy

Thyroid cancer (TC) is a malignancy that is increasing in prevalence on a global scale, necessitating the development of innovative approaches for both diagnosis and treatment. Myo-inositol (MI) plays a crucial role in a wide range of physiological and pathological functions within human cells. To date, studies have investigated the function of MI in thyroid physiology as well as its potential therapeutic benefits for hypothyroidism and autoimmune thyroiditis. However, research in the field of TC is very restricted. Metabolomics studies have highlighted the promising diagnostic capabilities of MI, recognizing it as a metabolic biomarker for identifying thyroid tumors. Furthermore, MI can influence therapeutic characteristics by modulating key cellular pathways involved in TC. This review evaluates the potential application of MI as a naturally occurring compound in the management of thyroid diseases, including hypothyroidism, autoimmune thyroiditis, and especially TC. The limited number of studies conducted in the field of TC emphasizes the critical need for future research to comprehend the multifaceted role of MI in TC. A significant amount of research and clinical trials is necessary to understand the role of MI in the pathology of TC, its diagnostic and therapeutic potential, and to pave the way for personalized medicine strategies in managing this intricate disease.

Assessment of circulating proteins in thyroid cancer: Proteome-wide Mendelian randomization and colocalization analysis

The causality between circulating proteins and thyroid cancer (TC) remains unclear. We employed five large-scale circulating proteomic genome-wide association studies (GWASs) with up to 100,000 participants and a TC meta-GWAS (nCase = 3,418, nControl = 292,703) to conduct proteome-wide Mendelian randomization (MR) and Bayesian colocalization analysis. Protein and gene expressions were validated in thyroid tissue. Through MR analysis, we identified 26 circulating proteins with a putative causal relationship with TCs, among which NANS protein passed multiple corrections (P BH = 3.28e-5, 0.05/1,525). These proteins were involved in amino acids and organic acid synthesis pathways. Colocalization analysis further identified six proteins associated with TCs (VCAM1, LGMN, NPTX1, PLEKHA7, TNFAIP3, and BMP1). Tissue validation confirmed BMP1, LGMN, and PLEKHA7's differential expression between normal and TC tissues. We found limited evidence for linking circulating proteins and the risk of TCs. Our study highlighted the contribution of proteins, particularly those involved in amino acid metabolism, to TCs.

Identification of suitable reference genes for RT-qPCR studies in human parathyroid tissue glandular cells

Identifying a proper reference gene allows us to understand fundamental changes in many biological processes. Normalization during gene expression analyses is essential for every tissue/cell type, including parathyroid tissue glandular cells. Quantitative method of gene expression analyses via qRT-PCR method provides the accurate examination of every target gene. There are limited reports to present commonly used reference genes in human parathyroid tissues rather than for glandular cell types. This study aims to determine and compare the most stable to least stable genes for parathyroid tissue cells. 43 human parathyroid tissue obtained from primary and secondary hyperparathyroidism patients and glandular cells isolated enzymatically by the removal of extracellular matrix components. After extraction of the total RNA, cDNA synthesis was performed, then qRT-PCR evaluated 14 candidate reference genes. Stability was determined by RefFinder software (Delta ct, BestKeeper, Genorm, and NormFinder algorithms), and the outcome was evaluated for five groups. Even if assessed with different groups, the most stable genes were RPLP0 and GAPDH, while the CLTC and RNA 18S were the least stable. We have confirmed the comprehensive ranking of the most stable three genes alone with the NormFinder algorithm to understand intergroup variation and found out that RPLP0>GAPDH>PGK1. Lastly, comparisons of relative target gene (GCM2) expression revealed similar expression patterns for the most stable reference genes. The most stable reference gene is recommended for the stages where stability is evaluated using the results of four different approaches using RefFinder. We aspire for this study to assist future research to conduct thorough assessments of appropriate reference genes before engaging in gene expression analyses for parathyroid tissue.

The role of rs2236242 at SERPINA12 gene and vaspin level on papillary thyroid carcinoma

Several studies showed that adipokines are associated with types of cancer which are documented to be effective in cancer biology. This study aimed to determine the relationship between vaspin rs2236242 polymorphism and the vaspin level with papillary thyroid carcinoma (PTC), and multinodular goiter (MNG). In this cross-sectional study, we recruited 170 candidates. Ninety patients with newly diagnosed (PTC 60 patients and MNG 30 patients), and 80 participants as a control group referred to Shariati Hospital, Tehran, Iran, were enrolled in the study. The vaspin hormone measurements were conducted utilizing the Elisa Kit. Using Tetra amplification resistant-mutation system polymerase chain reaction (T-ARMS-PCR), the genotype of single nucleotide polymorphism (SNP) rs2236242 was determined. The statistical analysis was performed using SPSS software version 20. Our findings showed significant age and genotype frequency differences in three groups (p-value < 0.05). There was no significant difference in vaspin levels between PTC, and control groups. The level of vaspin in MNG compared to the control group had significantly different, but there were no differences after adjustment for age. Results showed the genotypes of vaspin rs2236242 polymorphism are not associated with the level of vaspin. The genotypes and allele frequencies of vaspin rs2236242 in the PTC and MNG groups were significant compared to the control group. We have found vaspin rs2236242 gene polymorphism as a potential marker of papillary thyroid cancer. The A allele of the vaspin SNP rs2236242 plays a protective role against PTC and MNG. SNP at rs2236242 was not significantly associated with vaspin levels.

Evaluation of reference genes for gene expression analysis by real-time quantitative PCR (qPCR) in different tissues from mice infected by Ascaris suum

Human ascariasis is the most prevalent helminth infection, affecting 445 million people worldwide. To better understand the impact of the immune system on the pathophysiology of individuals infected with Ascaris suum, mice have been used as experimental models. The RT-qPCR technique is a critical auxiliary tool of investigation used to quantify mRNA levels. However, proper normalization using reference genes is essential to ensure reliable outcomes to avoid analytical errors and false results. Despite the importance of reference genes for experimental A. suum infection studies, no specific reference genes have been identified yet. Therefore, we conducted a study to assess five potential reference genes (GAPDH, 18s, ACTB, B2M, and HPRT1) in different tissues (liver, lungs, small and large intestines) affected by A. suum larval migration in C57BL/6j mice. Tissue collection was carried out to analyze parasite burden and confirm the presence of larvae during the peak of migration in each tissue. Upon confirmation, we analyzed different genes in the tissues and found no common gene with stable expression. Our results highlight the importance of analyzing different genes and using different software programs to ensure reliable relative expression results. Based on our findings, B2M was ranked as the ideal reference gene for the liver, while 18S was the most stable gene in the lung and small intestine. ACTB, or a combination of ACTB with GAPDH, was deemed suitable as reference genes for the large intestine due to their stable expression and less variation between the control and infected groups. To further demonstrate the impact of using different reference genes, we normalized the expression of a chemokine gene (CXCL9) in all tissues. Significant differences in CXCL9 expression levels were observed between different groups in all tissues except for the large intestine. This underscores the importance of selecting appropriate reference genes to avoid overestimating target gene expression levels and encountering normalization-related issues that can lead to false results. In conclusion, our study highlights the significance of using reliable reference genes for accurate RT-qPCR analysis, especially in the context of A. suum infection studies in different tissues. Proper normalization is crucial to ensure the validity of gene expression data and avoid potential pitfalls in interpreting results.

Selection and Validation of Reference Genes in Sudan Grass (Sorghum sudanense (Piper) Stapf) under Various Abiotic Stresses by qRT-PCR

Quantitative reverse transcription PCR (qRT-PCR) can screen applicable reference genes of species, and reference genes can be used to reduce experimental errors. Sudan grass (Sorghum sudanense (Piper) Stapf) is a high-yield, abiotic-tolerant annual high-quality forage with a wide range of uses. However, no studies have reported reference genes suitable for Sudan grass. Therefore, we found eight candidate reference genes, including UBQ10, HIS3, UBQ9, Isoform0012931, PP2A, ACP2, eIF4α, and Actin, under salt stress (NaCl), drought stress (DR), acid aluminum stress (AlCl3), and methyl jasmonate treatment (MeJA). By using geNorm, NormFinder, BestKeeper, and RefFinder, we ranked eight reference genes on the basis of their expression stabilities. The results indicated that the best reference gene was PP2A under all treatments. eIF4α can be used in CK, MeJA, NaCl, and DR. HIS3 can serve as the best reference gene in AlCl3. Two target genes (Isoform0007606 and Isoform0002387) belong to drought-stress-response genes, and they are highly expressed in Sudan grass according to transcriptome data. They were used to verify eight candidate reference genes under drought stress. The expression trends of the two most stable reference genes were similar, but the trend in expression for Actin showed a significant difference. The reference genes we screened provided valuable guidance for future research on Sudan grass.

Analyzing the postulated inhibitory effect of Manumycin A on farnesyltransferase

Manumycin A is postulated to be a specific inhibitor against the farnesyltransferase (FTase) since this effect has been shown in 1993 for yeast FTase. Since then, plenty of studies investigated Manumycin A in human cells as well as in model organisms like Caenorhabditis elegans. Some studies pointed to additional targets and pathways involved in Manumycin A effects like apoptosis. Therefore, these studies created doubt whether the main mechanism of action of Manumycin A is FTase inhibition. For some of these alternative targets half maximal inhibitory concentrations (IC₅₀) of Manumycin A are available, but not for human and C. elegans FTase. So, we aimed to 1) characterize missing C. elegans FTase kinetics, 2) elucidate the IC₅₀ and K_i values of Manumycin A on purified human and C. elegans FTase 3) investigate Manumycin A dependent expression of FTase and apoptosis genes in C. elegans. C. elegans FTase has its temperature optimum at 40°C with K_M of 1.3 µM (farnesylpyrophosphate) and 1.7 µM (protein derivate). Whilst other targets are inhibitable by Manumycin A at the nanomolar level, we found that Manumycin A inhibits cell-free FTase in micromolar concentrations (K_{i human} 4.15 μM; K_{i C. elegans} 3.16 μM). Furthermore, our gene expression results correlate with other studies indicating that thioredoxin reductase 1 is the main target of Manumycin A. According to our results, the ability of Manumycin A to inhibit the FTase at the micromolar level is rather neglectable for its cellular effects, so we postulate that the classification as a specific FTase inhibitor is no longer valid.

Identification and Validation of Common Reference Genes for Normalization of Esophageal Squamous Cell Carcinoma Gene Expression Profiles

Esophageal squamous cell carcinoma (ESCC) is one of the subtypes of esophageal cancer with Chinese characteristics, and its five-year survival rate is less than 20%. Early diagnosis is beneficial to improving the survival rate of ESCC significantly. Quantitative Real-Time Polymerase Chain Reaction is a high-throughput technique that can quantify tumor-related genes for early diagnosis. Its accuracy largely depends on the stability of the reference gene. There is no systematic scientific basis to demonstrate which reference gene expression is stable in ESCC and no consensus on the selection of internal reference. Therefore, this research used four software programs (The comparative delta-Ct method, GeNorm, NormFinder, and BestKeeper) to evaluate the expression stability of eight candidate reference genes commonly used in other tumor tissues and generated a comprehensive analysis by RefFinder. Randomly selected transcriptome sequencing analysis confirmed the SPP1 gene is closely related to ESCC. It was found that the expression trend of SPP1 obtained by RPS18 and PPIA as internal reference genes were the same as that of sequencing. The results show that RPS18 and PPIA are stable reference genes, and PPIA + RPS18 are a suitable reference gene combination. This is a reference gene report that combines transcriptome sequencing analysis and only focuses on ESCC, which makes the quantification more precise, systematic, and standardized, and promotes gene regulation research and the early diagnosis of ESCC in the future.

Amelioration for an ignored pitfall in reference gene selection by considering the mean expression and standard deviation of target genes

Routine tissue-specific reference genes are often used in expression studies, but target genes are not taken into account. Using the relative RT-qPCR approach, we evaluated the expression of three target genes. At the same time, meta-analyses were conducted in various ethnic groups, genders, and thyroid cancer subtypes. When eight common reference genes were examined, it was discovered that some of them not only lacked consistent expression but also had considerable expression variance. It is worth noting that while choosing a reference gene, the mean gene expression and its standard deviation should be carefully addressed. An equation was developed based on this, and it was used to perform statistical analysis on over 25,000 genes. According to the subtype of thyroid cancer and, of course, the target genes in this investigation, appropriate reference genes were proposed. The intuitive choice of GAPDH as a common reference gene caused a major shift in the quantitative expression data of target genes, inverting the relative expression values. As a result, choosing the appropriate reference gene(s) for quantification of transcription data, and especially for relative studies of the expression of target gene(s), is critical and should be carefully considered during the study design.

Thyroid hormone economy in mice overexpressing iodothyronine deiodinases

In peripheral tissues, triiodothyronine (T3) production and consequent thyroid hormone actions are mainly regulated by iodothyronine deiodinases (DIOs) classified into 3 types: D1, D2, and D3. We aimed to investigate the effects of peripheral DIOs on thyroid hormone economy independent of the hypothalamus-pituitary-thyroid axis. We cloned coding sequences of human DIOs with FLAG-tag and HiBiT-tag sequences into a pcDNA3.1 vector. To obtain full-length proteins, we modified these vectors by cloning the selenocysteine insertion sequence of each DIO (SECIS vectors). Western blot analyses and HiBiT lytic assay using HEK293T cells revealed that SECIS vectors expressed full-length proteins with substantial activity. Subsequently, in vivo transfections of pLIVE-based SECIS vectors into male C57BL/6J mice were performed by hydrodynamic gene delivery to generate mice overexpressing DIOs predominantly in the liver (D1, D2, and D3 mice). After 7 days from transfections, mice were analyzed to clarify phenotypes. To summarize, serum thyroid hormone levels did not change in D1 mice but D2 mice had higher serum free T3 levels. D3 mice developed hypothyroidism with higher serum reverse T3 (rT3) levels. Transfections with levothyroxine administration suggested that thyroid hormone action was upregulated in D2 mice. Our DIO-overexpressing mice provided insights on the physiological properties of upregulated DIOs: D2 augments local thyroid hormone action and recruits T3 into the circulation: D3 decreases circulating T3 and T4 levels with elevated rT3, leading to consumptive hypothyroidism. As D3 mice are expected to be a novel hypothyroidism model, they can contribute to progress in the field of thyroid hormone economy and action.

Bioengineering the ameloblastoma tumour to study its effect on bone nodule formation

Ameloblastoma is a benign, epithelial cancer of the jawbone, which causes bone resorption and disfigurement to patients affected. The interaction of ameloblastoma with its tumour stroma drives invasion and progression. We used stiff collagen matrices to engineer active bone forming stroma, to probe the interaction of ameloblastoma with its native tumour bone microenvironment. This bone-stroma was assessed by nano-CT, transmission electron microscopy (TEM), Raman spectroscopy and gene analysis. Furthermore, we investigated gene correlation between bone forming 3D bone stroma and ameloblastoma introduced 3D bone stroma. Ameloblastoma cells increased expression of MMP-2 and -9 and RANK temporally in 3D compared to 2D. Our 3D biomimetic model formed bone nodules of an average surface area of 0.1 mm² and average height of 92.37 [Formula: see text] 7.96 μm over 21 days. We demonstrate a woven bone phenotype with distinct mineral and matrix components and increased expression of bone formation genes in our engineered bone. Introducing ameloblastoma to the bone stroma, completely inhibited bone formation, in a spatially specific manner. Multivariate gene analysis showed that ameloblastoma cells downregulate bone formation genes such as RUNX2. Through the development of a comprehensive bone stroma, we show that an ameloblastoma tumour mass prevents osteoblasts from forming new bone nodules and severely restricted the growth of existing bone nodules. We have identified potential pathways for this inhibition. More critically, we present novel findings on the interaction of stromal osteoblasts with ameloblastoma.

HMBS is the most suitable reference gene for RT-qPCR in human HCC tissues and blood samples

Reverse transcription-quantitative (RT-q) PCR is the most feasible and useful technique for identifying and evaluating cancer biomarkers; however, the method requires suitable reference genes for gene expression analysis. The aim of the present study was to identify the most suitable reference gene for the normalization of relative gene expression in human hepatocellular carcinoma (HCC) tissue and blood samples. First, 14 candidate reference genes were selected through a systematic literature search. The expression levels of these genes (ACTB, B2M, GAPDH, GUSB, HMBS, HPRT1, PGK1, PPIA, RPLP0, RPL13A, SDHA, TBP, TFRC and YWHAZ) were evaluated using human multistage HCC transcriptome data (dataset GSE114564), which included normal liver (n=15), chronic hepatitis (n=20), liver cirrhosis (n=10), and early (n=18) and advanced HCC (n=45). From the 14 selected genes, five genes, whose expression levels were stable in all liver disease statuses (ACTB, GAPDH, HMBS, PPIA and RPLP0), were further assessed using RT-qPCR in 40 tissues (20 paired healthy tissues and 20 tissues from patients with HCC) and 40 blood samples (20 healthy controls and 20 samples from patients with HCC). BestKeeper statistical algorithms were used to identify the most stable reference genes, of which HMBS was found to be the most stable in both HCC tissues and blood samples. Therefore, the results of the present study suggest HMBS as a promising reference gene for the normalization of relative RT-qPCR techniques in HCC-related studies.

The gene expression of GPER1 is low in fresh samples of papillary thyroid carcinoma (PTC), and in silico analysis

Papillary thyroid cancer (PTC), whose incidence has been increasing in the last years, occurs more frequently in women. Experimental studies suggested that estrogen could be an important risk factor for the higher female incidence. In fact, it has been demonstrated that 17β-estradiol (E2) could increase proliferation and dedifferentiation in thyroid follicular cells. Genomic estrogen responses are typically mediated through classical estrogen receptors, the α and β isoforms, which have been described in normal and abnormal human thyroid tissue. Nevertheless, effects mediated through G protein estrogen receptor 1 (GPR30/GPER/GPER1), described in some thyroid cancer cell lines, could be partially responsible for the regulation of growth in normal cells. In this study, GPER1 gene and protein expression are described in non-malignant and in papillary thyroid cancer (PTC), as well as its association with clinical features of patients with PTC. The GPER1 expression was lower in PTC as compared to paired non-malignant thyroid tissues in fresh samples of PTC and in silico analysis of GEO and TCGA databases. In PTC cases of TCGA database, low GPER1 mRNA expression was independently associated with metastatic lymph nodes, female gender, and BRAF mutation. Besides, GPER1 mRNA levels were positively correlated with mRNA levels of thyroid differentiation genes. These results support the hypothesis that GPER1 have a role in PTC tumorigenesis and might be a potential target for its therapy. Further studies are needed to determine the functionality of these receptors in normal and diseased thyroid.

New evidence on tumor suppressor activity of PTEN and KLLN in papillary thyroid carcinoma

This study aimed to address the hypothesis that the expression of PTEN and KLLN tumor suppressor genes could diminish in papillary thyroid cancer (PTC) compared to paired normal tissue (PNT) and multinodular goiter (MNG). PTEN and KLLN expressions were assessed at both mRNA and protein levels in 82 tissue samples, including 30 PTC, 30 PNT, and 26 MNG using SYBR-Green Real-Time PCR and enzyme-linked immunosorbent assay (ELISA), respectively. Bioinformatics studies were performed to evaluate the genomic location and the genes promoter region. The mRNA expression of PTEN and KLLN in PTC was significantly lower than PNT (PTEN, P = 0.0033; KLLN, P = 0.0005). A significant decrease in the mRNA level of KLLN was also observed in PTC than MNG (P = 0.0304). Decreased level of PTEN mRNA (odds ratio=0.391; P = 0.013) or KLLN mRNA (odds ratio=0.023; P = 0.025) was associated with an increased risk of PTC tumorigenesis. Areas under the ROC curve for PTEN and KLLN were 0.69 and 0.78, respectively. PTEN and KLLN protein expressions in PTC compared to PNT or MNG were not significantly different. The bioinformatics studies revealed the sequence near the promoter region is lowly conserved across species. Four GC boxes were found upstream of the PTEN transcription start site (TSS), and one TATA box and one GC box were found upstream of KLLN TSS. The results suggest PTEN and KLLN are the two tumor suppressor genes that decreasing or loss of both of them occurs in sporadic PTC tumorigenesis. It appears they could have a promising application in both diagnostic and therapeutic areas.

G Protein-coupled Receptors in Radioiodine-refractory Thyroid Cancer in the Era of Precision Medicine

CONTEXT

Radioiodine-refractory thyroid cancers have poor outcomes and limited therapeutic options (tyrosine kinase inhibitors) due to transient efficacy and toxicity of treatments. Therefore, combinatorial treatments with new therapeutic approaches are needed. Many studies link G protein-coupled receptors (GPCRs) to cancer cell biology.

OBJECTIVE

To perform a specific atlas of GPCR expression in progressive and refractory thyroid cancer to identify potential targets among GPCRs aiming at drug repositioning.

METHODS

We analyzed samples from tumor and normal thyroid tissues from 17 patients with refractory thyroid cancer (12 papillary thyroid cancers [PTCs] and 5 follicular thyroid cancers [FTCs]). We assessed GPCR mRNA expression using NanoString technology with a custom panel of 371 GPCRs. The data were compared with public repositories and pharmacological databases to identify eligible drugs. The analysis of prognostic value of genes was also performed with TCGA datasets.

RESULTS

With our transcriptomic analysis, 4 receptors were found to be downregulated in FTC (VIPR1, ADGRL2/LPHN2, ADGRA3, and ADGRV1). In PTC, 24 receptors were deregulated, 7 of which were also identified by bioinformatics analyses of publicly available datasets on primary thyroid cancers (VIPR1, ADORA1, GPRC5B, P2RY8, GABBR2, CYSLTR2, and LPAR5). Among all the differentially expressed genes, 22 GPCRs are the target of approved drugs and some GPCRs are also associated with prognostic factors.

DISCUSSION

For the first time, we performed GPCR mRNA expression profiling in progressive and refractory thyroid cancers. These findings provide an opportunity to identify potential therapeutic targets for drug repositioning and precision medicine in radioiodine-refractory thyroid cancer.

Validation of reference genes for use in untreated bovine fibroblasts

Proper normalization of RT-qPCR data is pivotal to the interpretation of results and accuracy of scientific conclusions. Though different approaches may be taken, normalization against multiple reference genes is now standard practice. Genes traditionally used and deemed constitutively expressed have demonstrated variability in expression under different experimental conditions, necessitating the proper validation of reference genes prior to utilization. Considering the wide use of fibroblasts in research and scientific applications, it is imperative that suitable reference genes for fibroblasts of different animal origins and conditions be elucidated. Previous studies on bovine fibroblasts have tested limited genes and/or samples. Herein, we present an extensive study investigating the expression stability of 16 candidate reference genes across 7 untreated bovine fibroblast cell lines subjected to controlled conditions. Data were analysed using various statistical tools and algorithms, including geNorm, NormFinder, BestKeeper, and RefFinder. A combined use of GUSB and RPL13A was determined to be the best approach for data normalization in untreated bovine fibroblasts.

GAPDH and PUM1: Optimal Housekeeping Genes for Quantitative Polymerase Chain Reaction-Based Analysis of Cancer Stem Cells and Epithelial-Mesenchymal Transition Gene Expression in Rectal Tumors

Background

The overwhelming majority of published articles have taken colon and rectal cancer as a single group, i.e., colorectal cancer, when normalizing gene expression data with housekeeping genes (HKG) in quantitative polymerase chain reaction (qPCR) experiments though there are published reports that suggest the differential expression pattern of genes between the colon and rectal cancer groups and hence the current experiment was attempted to find out the optimal set of housekeeping genes from the list of common HKG for rectal tumor gene expression analysis.

Methods

The expression of five potential housekeeping genes GAPDH, RPNI, PUM1, B2M, and PMM1 was analyzed through qPCR and Bestkeeper software (http://www.wzw.tum.de/gene-quantification/bestkeeper.html) in 20 stage II-IV rectal cancer samples to check for uniformity in their expression pattern. Cancer stem cell (CSC) marker ALDH1 and epithelial-mesenchymal transition marker (EMT) markers E cadherin, vimentin, Twist, and SNAI2 expression were evaluated in conjunction with the two optimal reference genes in 10 rectal cancers as part of validation.

Results

The standard deviation of the cycle threshold value of GAPDH was found the lowest at 0.65 followed by RPN1 at 0.88, PUM1 at 0.94, PMM1 at 0.94, and B2M at 1.21 when analyzed with BestKeeper software. Using GAPDH and PUM1 as the reference gene for the validation phase, rectal cancer patients with stage III/IV showed a 4.79-fold change (P=0.006) in ALDH1 expression, and an 11.76-fold change in Twist expression (P=0.003) with respect to stage II rectal tumor when normalized with GAPDH and PUM1.

Conclusion

GAPDH and PUM1 can be used as an optimal set of housekeeping genes for gene expression-related experiments in rectal tumors. ALDH1 and Twist were found significantly overexpressed in stage III/IV rectal tumors in comparison to stage II rectal cancer. Genes associated with cancer stem cells and EMT markers could be optimally analyzed by normalizing them with GAPDH and PUM1 as housekeeping genes.

The pseudogene problem and RT-qPCR data normalization; SYMPK: a suitable reference gene for papillary thyroid carcinoma

In RT-qPCR, accuracy requires multiple levels of standardization, but results could be obfuscated by human errors and technical limitations. Data normalization against suitable reference genes is critical, yet their observed expression can be confounded by pseudogenes. Eight reference genes were selected based on literature review and analysis of papillary thyroid carcinoma (PTC) microarray data. RNA extraction and cDNA synthesis were followed by RT-qPCR amplification in triplicate with exon-junction or intron-spanning primers. Several statistical analyses were applied using Microsoft Excel, NormFinder, and BestKeeper. In normal tissues, the least correlation of variation (CqCV%) and the lowest maximum fold change (MFC) were respectively recorded for PYCR1 and SYMPK. In PTC tissues, SYMPK had the lowest CqCV% (5.16%) and MFC (1.17). According to NormFinder, the best reference combination was SYMPK and ACTB (stability value = 0.209). BestKeeper suggested SYMPK as the best reference in both normal (r = 0.969) and PTC tissues (r = 0.958). SYMPK is suggested as the best reference gene for overcoming the pseudogene problem in RT-qPCR data normalization, with a stability value of 0.319.

Selection of endogenous control genes for normalising gene expression data derived from formalin-fixed paraffin-embedded tumour tissue

Quantitative real time polymerase chain reaction (qPCR) data are normalised using endogenous control genes. We aimed to: (1) demonstrate a pathway to identify endogenous control genes for qPCR analysis of formalin-fixed paraffin-embedded (FFPE) tissue using bladder cancer as an exemplar; and (2) examine the influence of probe length and sample age on PCR amplification and co-expression of candidate genes on apparent expression stability. RNA was extracted from prospective and retrospective samples and subject to qPCR using TaqMan human endogenous control arrays or single tube assays. Gene stability ranking was assessed using coefficient of variation (CoV), GeNorm and NormFinder. Co-expressed genes were identified from The Cancer Genome Atlas (TCGA) using the on-line gene regression analysis tool GRACE. Cycle threshold (Ct) values were lower for prospective (19.49 ± 2.53) vs retrospective (23.8 ± 3.32) tissues (p < 0.001) and shorter vs longer probes. Co-expressed genes ranked as the most stable genes in the TCGA cohort by GeNorm when analysed together but ranked lower when analysed individually omitting co-expressed genes indicating bias. Stability values were < 1.5 for the 20 candidate genes in the prospective cohort. As they consistently ranked in the top ten by CoV, GeNorm and Normfinder, UBC, RPLP0, HMBS, GUSB, and TBP are the most suitable endogenous control genes for bladder cancer qPCR.

Selection and validation of reference genes for measuring gene expression in Toona ciliata under different experimental conditions by quantitative real-time PCR analysis

BACKGROUND

Before studying gene expression of different organisms, it is important to determine the best reference gene. At present, the most accurate method of detecting gene expression is quantitative real-time PCR (RT-qPCR). With this method, reference genes that are stable in different biological systems and under different conditions can be obtained. Toona ciliata Roem (T. ciliata). is a valuable and fast-growing timber specie. In this study, 20 reference genes were identified using RT-qPCR, as a primary prerequisite for future gene expression analysis. Four different methods, geNorm, NormFinder, BestKeeper, and RankAggreg were used to evaluate the expression stability of the 20 candidate reference genes in various tissues under different conditions.

RESULTS

The experimental results showed that TUB-α was the most stably expressed reference gene across all samples and UBC17 was the most stable in leaves and young stems under Hypsipyla robusta (H. robusta) and methyl jasmonate (MeJA) treatments. In addition, PP2C59 and UBC5B were the best-performing genes in leaves under H. robusta treatment, while HIS1 and ACT7 were the best reference genes in young stems. The two best reference genes were 60S-18 and TUB-α after treatment at 4 °C. The expression of HIS6 and MUB1 was the most stable under PEG6000 treatment. The accuracy of the selected reference genes was verified using the transcription factor MYB3 (TcMYB3) gene.

CONCLUSIONS

This is the first report to verify the best reference genes for normalizing gene expression in T. ciliata under different conditions, which will facilitate future elucidation of gene regulations in this species.

Plasma Metabolic Profiling of Human Thyroid Nodules by Gas Chromatography-Mass Spectrometry (GC-MS)-Based Untargeted Metabolomics

One of the challenges in the area of diagnostics of human thyroid cancer is a preoperative diagnosis of thyroid nodules with indeterminate cytology. Herein, we report an untargeted metabolomics analysis to identify circulating thyroid nodule metabolic signatures, to find new novel metabolic biomarkers. Untargeted gas chromatography-quadrupole-mass spectrometry was used to ascertain the specific plasma metabolic changes of thyroid nodule patients, which consisted of papillary thyroid carcinoma (PTC; n = 19), and multinodular goiter (MNG; n = 16), as compared to healthy subjects (n = 20). Diagnostic models were constructed using multivariate analyses such as principal component analysis, orthogonal partial least squares-discriminant analysis, and univariate analysis including One-way ANOVA and volcano plot by MetaboAnalyst and SIMCA software. Because of the multiple-testing issue, false discovery rate p-values were also computed for these functions. A total of 60 structurally annotated metabolites were subjected to statistical analysis. A combination of univariate and multivariate statistical analyses revealed a panel of metabolites responsible for the discrimination between thyroid nodules and healthy subjects, with variable importance in the projection (VIP) value greater than 0.8 and p-value less than 0.05. Significantly altered metabolites between thyroid nodules versus healthy persons are those associated with amino acids metabolism, the tricarboxylic acid cycle, fatty acids, and purine and pyrimidine metabolism, including cysteine, cystine, glutamic acid, α-ketoglutarate, 3-hydroxybutyric acid, adenosine-5-monophosphate, and uracil, respectively. Further, sucrose metabolism differed profoundly between thyroid nodule patients and healthy subjects. Moreover, according to the receiver operating characteristic (ROC) curve analysis, sucrose could discriminate PTC from MNG (area under ROC curve value = 0.92). This study enhanced our understanding of the distinct metabolic pathways associated with thyroid nodules, which enabled us to distinguish between patients and healthy subjects. In addition, our study showed extensive sucrose metabolism in the plasma of thyroid nodule patients, which provides a new metabolic signature of the thyroid nodule’s tumorigenesis. Accordingly, it suggests that sucrose can be considered as a circulating biomarker for differential diagnosis between malignancy and benignity in indeterminate thyroid nodules.

Catalogue of stage-specific transcripts in Ixodes ricinus and their potential functions during the tick life-cycle

Background

The castor bean tick Ixodes ricinus is an important vector of several clinically important diseases, whose prevalence increases with accelerating global climate changes. Characterization of a tick life-cycle is thus of great importance. However, researchers mainly focus on specific organs of fed life stages, while early development of this tick species is largely neglected.

Methods

In an attempt to better understand the life-cycle of this widespread arthropod parasite, we sequenced the transcriptomes of four life stages (egg, larva, nymph and adult female), including unfed and partially blood-fed individuals. To enable a more reliable identification of transcripts and their comparison in all five transcriptome libraries, we validated an improved-fit set of five I. ricinus-specific reference genes for internal standard normalization of our transcriptomes. Then, we mapped biological functions to transcripts identified in different life stages (clusters) to elucidate life stage-specific processes. Finally, we drew conclusions from the functional enrichment of these clusters specifically assigned to each transcriptome, also in the context of recently published transcriptomic studies in ticks.

Results

We found that reproduction-related transcripts are present in both fed nymphs and fed females, underlining the poorly documented importance of ovaries as moulting regulators in ticks. Additionally, we identified transposase transcripts in tick eggs suggesting elevated transposition during embryogenesis, co-activated with factors driving developmental regulation of gene expression. Our findings also highlight the importance of the regulation of energetic metabolism in tick eggs during embryonic development and glutamate metabolism in nymphs.

Conclusions

Our study presents novel insights into stage-specific transcriptomes of I. ricinus and extends the current knowledge of this medically important pathogen, especially in the early phases of its development.