Multi-Tiered Assessment of Gene Expression Provides Evidence for Mechanisms That Underlie Risk for Type 2 Diabetes

Introduction

Integrated transcriptome and microRNA differential gene expression (DEG) analyses may help to explain type 2 diabetes (T2D) pathogenesis in at-risk populations. The purpose of this study was to characterize DEG in banked biospecimens from underactive adult participants who responded to a randomized clinical trial measuring the effects of lifestyle interventions on T2D risk factors. DEGs were further examined within the context of annotated biological pathways.

Methods

Participants (n = 52) in a previously completed clinical trial that assessed a 12-week behavioural intervention for T2D risk reduction were included. Participants who showed >6mg/dL decrease in fasting blood glucose were identified as responders. Gene expression was measured by RNASeq, and overrepresentation analysis within KEGG pathways and weighted gene correlation network analysis (WGCNA) were performed.

Results

No genes remained significantly differentially expressed after correction for multiple comparisons. One module derived by WGCNA related to body mass index was identified, which contained genes located in KEGG pathways related to known mechanisms underlying risk for T2D as well as pathways related to neurodegeneration and protein misfolding. A network analysis showed indirect connections between genes in this module and islet amyloid polypeptide (IAPP), which has previously been hypothesized as a mechanism for T2D.

Discussion

We validated prior studies that showed pathways related to metabolism, inflammation/immunity, and endocrine/hormone function are related to risk for T2D. We identified evidence for new potential mechanisms that include protein misfolding. Additional studies are needed to determine whether these are potential therapeutic targets to decrease risk for T2D.

MicroRNA biomarkers target genes and pathways associated with type 2 diabetes

AIMS/HYPOTHESIS

Our prior analysis of the Diabetes Prevention Program study identified a subset of five miRNAs that predict incident type 2 diabetes. The purpose of this study was to identify mRNAs and biological pathways targeted by these five miRNAs to elucidate potential mechanisms of risk and responses to the tested interventions.

METHODS

Using experimentally validated data from miRTarBase version 8.0 and R (2021), we identified mRNAs with strong evidence to be regulated by individual or combinations of the five predictor miRNAs. Overrepresentation of the mRNA targets was assessed in pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation database.

RESULTS

The five miRNAs targeted 167 pathways and 122 mRNAs. Nine of the pathways have known associations with type 2 diabetes: Insulin signaling, Insulin resistance, Diabetic cardiomyopathy, Type 2 diabetes, AGE-RAGE signaling in diabetic complications, HIF-1 signaling, TGF-beta signaling, PI3K/Akt signaling, and Adipocytokine signaling pathways. Vascular endothelial growth factor A (VEGFA) has prior genetic associations with risk for type 2 diabetes and was the most commonly targeted mRNA for this set of miRNAs.

CONCLUSIONS/INTERPRETATION

These findings show that miRNA predictors of incident type 2 diabetes target mRNAs and pathways known to underlie risk for type 2 diabetes. Future studies should evaluate miRNAs as potential therapeutic targets for preventing and treating type 2 diabetes.

Review of databases for experimentally validated human microRNA-mRNA interactions

MicroRNAs (miRs) may contribute to disease etiology by influencing gene expression. Numerous databases are available for miR target prediction and validation, but their functionality is varied, and outputs are not standardized. The purpose of this review is to identify and describe databases for cataloging validated miR targets. Using Tools4miRs and PubMed, we identified databases with experimentally validated targets, human data, and a focus on miR-messenger RNA (mRNA) interactions. Data were extracted about the number of times each database was cited, the number of miRs, the target genes, the interactions per database, experimental methodology and key features of each database. The search yielded 10 databases, which in order of most cited to least were: miRTarBase, starBase/The Encyclopedia of RNA Interactomes, DIANA-TarBase, miRWalk, miRecords, miRGator, miRSystem, miRGate, miRSel and targetHub. Findings from this review suggest that the information presented within miR target validation databases can be enhanced by adding features such as flexibility in performing queries in multiple ways, downloadable data, ongoing updates and integrating tools for further miR-mRNA target interaction analysis. This review is designed to aid researchers, especially those new to miR bioinformatics tools, in database selection and to offer considerations for future development and upkeep of validation tools. Database URL http://mirtarbase.cuhk.edu.cn/.

Abstract 601: Serological responses to tumor-associated adducts in pancreatic cancer

The incidence of pancreatic cancer has increased over the past several decades and is now the fourth leading cause of cancer death in the USA. Only 8.5% of the patients live beyond 5 years after their initial diagnosis. This dismal survival rate is partly explained by the fact that most pancreatic cancer cases are detected late when treatment options are limited. Novel biomarkers are crucially needed. It has long been recognized that cancer cells abnormally accumulate chemical groups covalently linked to proteins, DNA and other molecules. Humoral responses to these “adducts” have not been thoroughly explored. To investigate these responses, we developed an ELISA platform for the detection of serum IgG reactive to 93 adducts, including post-translational modifications known to contribute to tumorigenesis, progression, and metastasis. Our panel also includes oxidation-related modification, advanced glycation end products, and certain co-enzymes that qualify as adducts based on their binding properties. Using this assay, we measured anti-adduct antibodies in the serum of adult healthy donors (N=24; age range 40-80) as well as patients with pancreatic cancer (N=31; age range 50-90). Importantly, all patient specimens were collected before any treatment. Reactivity to all adducts was analyzed using a random forest predictive model with 10,000 bootstrapped decision trees to distinguish between the two groups. Results revealed a distinctive anti-adduct IgG reactivity profile for pancreatic cancer patients when compared to controls. Applying feature selection using the Boruta algorithm identified a reduced set of 19 target adducts recognized by IgG that most efficiently discriminated between healthy donors and pancreatic cancer cases. All 19 adducts were confirmed as IgG targets in cancer patients using multiple unpaired t-tests (p≤0,05). Aside from their predictive value, we reasoned that the detection of adduct-specific antibodies in cancer patients reflected their accumulation in the tumor cells. To test this hypothesis, we assessed the level of two representative target adducts by immunofluorescence staining in tumor tissue compared to control adjacent non-tumoral tissue or healthy tonsil. Results unequivocally show higher levels of these two adducts in cancer cells compared to non-cancer cells. Taken together, these preliminary findings support our hypothesis of the development of anti-adduct antibodies in pancreatic cancer patients. Furthermore, the specific antibody signature detected in these patients is consistent with the abnormal presence of corresponding adducts in the transformed cells. Our findings pave the way for the development of diagnostic tests based on the detection of serum IgG specific to tumor-associated adducts.

Citation Format: Talita Aguiar, Shunya Mashiko, Poulomi Roy, Max Dietzel, Kesava Asam, Bradley Aouizerat, Jeanine Genkinger, Helen Remotti, Emmanuel Zorn. Serological responses to tumor-associated adducts in pancreatic cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 601.

Neurotrophin Pathway Receptors NGFR and TrkA Control Perineural Invasion, Metastasis, and Pain in Oral Cancer

Oral squamous cell carcinoma (OSCC) patients suffer from poor survival due to metastasis or locoregional recurrence, processes that are both facilitated by perineural invasion (PNI). OSCC has higher rates of PNI than other cancer subtypes, with PNI present in 80% of tumors. Despite the impact of PNI on oral cancer prognosis and pain, little is known about the genes that drive PNI, which in turn drive pain, invasion, and metastasis. In this study, clinical data, preclinical, and in vitro models are leveraged to elucidate the role of neurotrophins in OSCC metastasis, PNI, and pain. The expression data in OSCC patients with metastasis, PNI, or pain demonstrate dysregulation of neurotrophin genes. TrkA and nerve growth factor receptor (NGFR) are focused, two receptors that are activated by NGF, a neurotrophin expressed at high levels in OSCC. It is demonstrated that targeted knockdown of these two receptors inhibits proliferation and invasion in an in vitro and preclinical model of OSCC, and metastasis, PNI, and pain. It is further determined that TrkA knockdown alone inhibits thermal hyperalgesia, whereas NGFR knockdown alone inhibits mechanical allodynia. Collectively the results highlight the ability of OSCC to co-opt different components of the neurotrophin pathway in metastasis, PNI, and pain.

Schwann Cells Induce Phenotypic Changes in Oral Cancer Cells

Head and neck cancer (HNC) is the seventh most common cancer worldwide, the majority being oral squamous cell carcinoma. Despite advances in cancer diagnosis and treatment, the survival rate of patients with HNC remains stagnant. The cancer-nerve interaction has been recognized as an important driver of cancer progression. Schwann cells, a type of peripheral glia, have been implicated in promoting cancer cell growth, migration, dispersion, and invasion into the nerve in many cancers. Here, it is demonstrated that the presence of Schwann cells makes oral cancer cells more aggressive by promoting their proliferation, extracellular matrix breakdown, and altering cell metabolism. Furthermore, oral cancer cells became larger, more circular, with more projections and nuclei following co-culturing with Schwann cells. RNA-sequencing analysis in oral cancer cells following exposure to Schwann cells shows corresponding changes in genes involved in the hallmarks of cancer and cell metabolism; the enriched KEGG pathways are spliceosome, RNA transport, cell cycle, axon guidance, signaling pathways regulating pluripotency of stem cells, cAMP signaling, WNT signaling, proteoglycans in cancer and PI3K-Akt signaling. Taken together, these results suggest a significant role for Schwann cells in facilitating oral cancer progression, highlighting their potential as a target to treat oral cancer progression.

Co‑expressed microRNAs, target genes and pathways related to metabolism, inflammation and endocrine function in individuals at risk for type 2 diabetes

MicroRNAs (miRNAs) may be considered important regulators of risk for type 2 diabetes (T2D). The aim of the present study was to identify novel sets of miRNAs associated with T2D risk, as well as their gene and pathway targets. Circulating miRNAs (n=59) were measured in plasma from participants in a previously completed clinical trial (n=82). An agnostic statistical approach was applied to identify novel sets of miRNAs with optimal co-expression patterns. In silico analyses were used to identify the messenger RNA and biological pathway targets of the miRNAs within each factor. A total of three factors of miRNAs were identified, containing 18, seven and two miRNAs each. Eight biological pathways were revealed to contain genes targeted by the miRNAs in all three factors, 38 pathways contained genes targeted by the miRNAs in two factors, and 55, 18 and two pathways were targeted by the miRNAs in a single factor, respectively (all q<0.05). The pathways containing genes targeted by miRNAs in the largest factor shared a common theme of biological processes related to metabolism and inflammation. By contrast, the pathways containing genes targeted by miRNAs in the second largest factor were related to endocrine function and hormone activity. The present study focused on the pathways uniquely targeted by each factor of miRNAs in order to identify unique mechanisms that may be associated with a subset of individuals. Further exploration of the genes and pathways related to these biological themes may provide insights about the subtypes of T2D and lead to the identification of novel therapeutic targets.

Brush swab as a noninvasive surrogate for tissue biopsies in epigenomic profiling of oral cancer

BACKGROUND

Oral squamous cell carcinoma (OSCC) has poor survival rates. There is a pressing need to develop more precise risk assessment methods to tailor clinical treatment. Epigenome-wide association studies in OSCC have not produced a viable biomarker. These studies have relied on methylation array platforms, which are limited in their ability to profile the methylome. In this study, we use MethylCap-Seq (MC-Seq), a comprehensive methylation quantification technique, and brush swab samples, to develop a noninvasive, readily translatable approach to profile the methylome in OSCC patients.

METHODS

Three OSCC patients underwent collection of cancer and contralateral normal tissue and brush swab biopsies, totaling 4 samples for each patient. Epigenome-wide DNA methylation quantification was performed using the SureSelectXT Methyl-Seq platform. DNA quality and methylation site resolution were compared between brush swab and tissue samples. Correlation and methylation value difference were determined for brush swabs vs. tissues for each respective patient and site (i.e., cancer or normal). Correlations were calculated between cancer and normal tissues and brush swab samples for each patient to determine the robustness of DNA methylation marks using brush swabs in clinical biomarker studies.

RESULTS

There were no significant differences in DNA yield between tissue and brush swab samples. Mapping efficiency exceeded 90% across all samples, with no differences between tissue and brush swabs. The average number of CpG sites with at least 10x depth of coverage was 2,716,674 for brush swabs and 2,903,261 for tissues. Matched tissue and brush swabs had excellent correlation (r = 0.913 for cancer samples and r = 0.951 for normal samples). The methylation profile of the top 1000 CpGs was significantly different between cancer and normal samples (mean p-value = 0.00021) but not different between tissues and brush swabs (mean p-value = 0.11).

CONCLUSIONS

Our results demonstrate that MC-Seq is an efficient platform for epigenome profiling in cancer biomarker studies, with broader methylome coverage than array-based platforms. Brush swab biopsy provides adequate DNA yield for MC-Seq, and taken together, our findings set the stage for development of a non-invasive methylome quantification technique for oral cancer with high translational potential.

R-Spondin 3 Regulates Mammalian Dental and Craniofacial Development

Development of the teeth requires complex signaling interactions between the mesenchyme and the epithelium mediated by multiple pathways. For example, canonical WNT signaling is essential to many aspects of odontogenesis, and inhibiting this pathway blocks tooth development at an early stage. R-spondins (RSPOs) are secreted proteins, and they mostly augment WNT signaling. Although RSPOs have been shown to play important roles in the development of many organs, their role in tooth development is unclear. A previous study reported that mutating Rspo2 in mice led to supernumerary lower molars, while teeth forming at the normal positions showed no significant anomalies. Because multiple Rspo genes are expressed in the orofacial region, it is possible that the relatively mild phenotype of Rspo2 mutants is due to functional compensation by other RSPO proteins. We found that inactivating Rspo3 in the craniofacial mesenchyme caused the loss of lower incisors, which did not progress beyond the bud stage. A simultaneous deletion of Rspo2 and Rspo3 caused severe disruption of craniofacial development from early stages, which was accompanied with impaired development of all teeth. Together, these results indicate that Rspo3 is an important regulator of mammalian dental and craniofacial development.

The REASON score: an epigenetic and clinicopathologic score to predict risk of poor survival in patients with early stage oral squamous cell carcinoma

BACKGROUND

Oral squamous cell carcinoma (OSCC) is a capricious cancer with poor survival rates, even for early-stage patients. There is a pressing need to develop more precise risk assessment methods to appropriately tailor clinical treatment. Genome-wide association studies have not produced a viable biomarker. However, these studies are limited by using heterogeneous cohorts, not focusing on methylation although OSCC is a heavily epigenetically-regulated cancer, and not combining molecular data with clinicopathologic data for risk prediction. In this study we focused on early-stage (I/II) OSCC and created a risk score called the REASON score, which combines clinicopathologic characteristics with a 12-gene methylation signature, to predict the risk of 5-year mortality.

METHODS

We combined data from an internal cohort (n = 515) and The Cancer Genome Atlas (TCGA) cohort (n = 58). We collected clinicopathologic data from both cohorts to derive the non-molecular portion of the REASON score. We then analyzed the TCGA cohort DNA methylation data to derive the molecular portion of the risk score.

RESULTS

5-year disease specific survival was 63% for the internal cohort and 86% for the TCGA cohort. The clinicopathologic features with the highest predictive ability among the two the cohorts were age, race, sex, tobacco use, alcohol use, histologic grade, stage, perineural invasion (PNI), lymphovascular invasion (LVI), and margin status. This panel of 10 non-molecular features predicted 5-year mortality risk with a concordance (c)-index = 0.67. Our molecular panel consisted of a 12-gene methylation signature (i.e., HORMAD2, MYLK, GPR133, SOX8, TRPA1, ABCA2, HGFAC, MCPH1, WDR86, CACNA1H, RNF216, CCNJL), which had the most significant differential methylation between patients who survived vs. died by 5 years. All 12 genes have already been linked to survival in other cancers. Of the genes, only SOX8 was previously associated with OSCC; our study was the first to link the remaining 11 genes to OSCC survival. The combined molecular and non-molecular panel formed the REASON score, which predicted risk of death with a c-index = 0.915.

CONCLUSIONS

The REASON score is a promising biomarker to predict risk of mortality in early-stage OSCC patients. Validation of the REASON score in a larger independent cohort is warranted.

Molecular patterning of the embryonic cranial mesenchyme revealed by genome-wide transcriptional profiling

In the head of an embryo, a layer of mesenchyme surrounds the brain underneath the surface ectoderm. This cranial mesenchyme gives rise to the meninges, the calvaria (top part of the skull), and the dermis of the scalp. Abnormal development of these structures, especially the meninges and the calvaria, is linked to significant congenital defects in humans. It has been known that different areas of the cranial mesenchyme have different fates. For example, the calvarial bone develops from the cranial mesenchyme on the baso-lateral side of the head just above the eye (supraorbital mesenchyme, SOM), but not from the mesenchyme apical to SOM (early migrating mesenchyme, EMM). However, the molecular basis of this difference is not fully understood. To answer this question, we compared the transcriptomes of EMM and SOM using high-throughput sequencing (RNA-seq). This experiment identified a large number of genes that were differentially expressed in EMM and SOM, and gene ontology analyses found very different terms enriched in each region. We verified the expression of about 40 genes in the head by RNA in situ hybridization, and the expression patterns were annotated to make a map of molecular markers for 6 subdivisions of the cranial mesenchyme. Our data also provided insights into potential novel regulators of cranial mesenchyme development, including several axon guidance pathways, lectin complement pathway, cyclic-adenosine monophosphate (cAMP) signaling pathway, and ZIC family transcription factors. Together, information in this paper will serve as a unique resource to guide future research on cranial mesenchyme development.

Anti-osteogenic function of a LIM-homeodomain transcription factor LMX1B is essential to early patterning of the calvaria

The calvaria (upper part of the skull) is made of plates of bone and fibrous joints (sutures and fontanelles), and the proper balance and organization of these components are crucial to normal development of the calvaria. In a mouse embryo, the calvaria develops from a layer of head mesenchyme that surrounds the brain from shortly after mid-gestation. The mesenchyme just above the eye (supra-orbital mesenchyme, SOM) generates ossification centers for the bones, which then grow toward the apex gradually. In contrast, the mesenchyme apical to SOM (early migrating mesenchyme, EMM), including the area at the vertex, does not generate an ossification center. As a result, the dorsal midline of the head is occupied by sutures and fontanelles at birth. To date, the molecular basis for this regional difference in developmental programs is unknown. The current study provides vital insights into the genetic regulation of calvarial patterning. First, we showed that osteogenic signals were active in both EMM and SOM during normal development, which suggested the presence of an anti-osteogenic factor in EMM to counter the effect of these signals. Subsequently, we identified Lmx1b as an anti-osteogenic gene that was expressed in EMM but not in SOM. Furthermore, head mesenchyme-specific deletion of Lmx1b resulted in heterotopic ossification from EMM at the vertex, and craniosynostosis affecting multiple sutures. Conversely, forced expression of Lmx1b in SOM was sufficient to inhibit osteogenic specification. Therefore, we conclude that Lmx1b plays a key role as an anti-osteogenic factor in patterning the head mesenchyme into areas with different osteogenic competence. In turn, this patterning event is crucial to generating the proper organization of the bones and soft tissue joints of the calvaria.

Eicosanoyl-5-hydroxytryptamide (EHT) prevents Alzheimer's disease-related cognitive and electrophysiological impairments in mice exposed to elevated concentrations of oligomeric beta-amyloid

Soluble forms of oligomeric beta-amyloid (Aβ) are thought to play a central role in Alzheimer's disease (AD). Transgenic manipulation of methylation of the serine/threonine protein phosphatase, PP2A, was recently shown to alter the sensitivity of mice to AD-related impairments resulting from acute exposure to elevated levels of Aβ. In addition, eicosanoyl-5-hydroxytryptamide (EHT), a naturally occurring component from coffee beans that modulates PP2A methylation, was shown to confer therapeutic benefits in rodent models of AD and Parkinson's disease. Here, we tested the hypothesis that EHT protects animals from the pathological effects of exposure to elevated levels of soluble oligomeric Aβ. We treated mice with EHT-containing food at two different doses and assessed the sensitivity of these animals to Aβ-induced behavioral and electrophysiological impairments. We found that EHT administration protected animals from Aβ-induced cognitive impairments in both a radial-arm water maze and contextual fear conditioning task. We also found that both chronic and acute EHT administration prevented Aβ-induced impairments in long-term potentiation. These data add to the accumulating evidence suggesting that interventions with pharmacological agents, such as EHT, that target PP2A activity may be therapeutically beneficial for AD and other neurological conditions.