ASU Scholarship Showcase

Background: Healthy individuals on the lower end of the insulin sensitivity spectrum also have a reduced gene expression response to exercise for specific genes. The goal of this study was to determine the relationship between insulin sensitivity and exercise-induced gene expression in an unbiased, global manner.

Methods and Findings: Euglycemic clamps were used to measure insulin sensitivity and muscle biopsies were done at rest and 30 minutes after a single acute exercise bout in 14 healthy participants. Changes in mRNA expression were assessed using microarrays, and miRNA analysis was performed in a subset of 6 of the participants using sequencing techniques. Following exercise, 215 mRNAs were changed at the probe level (Bonferroni-corrected P<0.00000115). Pathway and Gene Ontology analysis showed enrichment in MAP kinase signaling, transcriptional regulation and DNA binding. Changes in several transcription factor mRNAs were correlated with insulin sensitivity, including MYC, r=0.71; SNF1LK, r=0.69; and ATF3, r= 0.61 (5 corrected for false discovery rate). Enrichment in the 5’-UTRs of exercise-responsive genes suggested regulation by common transcription factors, especially EGR1. miRNA species of interest that changed after exercise included miR-378, which is located in an intron of the PPARGC1B gene.

Conclusions: These results indicate that transcription factor gene expression responses to exercise depend highly on insulin sensitivity in healthy people. The overall pattern suggests a coordinated cycle by which exercise and insulin sensitivity regulate gene expression in muscle.

Although insulin resistance in skeletal muscle is well-characterized, the role of circulating whole blood in the metabolic syndrome phenotype is not well understood. We set out to test the hypothesis that genes involved in inflammation, insulin signaling and mitochondrial function would be altered in expression in the whole blood of individuals with metabolic syndrome. We further wanted to examine whether similar relationships that we have found previously in skeletal muscle exist in peripheral whole blood cells. All subjects (n=184) were Latino descent from the Arizona Insulin Resistance registry. Subjects were classified based on the metabolic syndrome phenotype according to the National Cholesterol Education Program’s Adult Treatment Panel III. Of the 184 Latino subjects in the study, 74 were classified with the metabolic syndrome and 110 were without. Whole blood gene expression profiling was performed using the Agilent 4x44K Whole Human Genome Microarray. Whole blood microarray analysis identified 1,432 probes that were altered in expression ≥1.2 fold and P<0.05 after Benjamini-Hochberg in the metabolic syndrome subjects. KEGG pathway analysis revealed significant enrichment for pathways including ribosome, oxidative phosphorylation and MAPK signaling (all Benjamini-Hochberg P<0.05). Whole blood mRNA expression changes observed in the microarray data were confirmed by quantitative RT-PCR. Transcription factor binding motif enrichment analysis revealed E2F1, ELK1, NF-kappaB, STAT1 and STAT3 significantly enriched after Bonferroni correction (all P<0.05). The results of the present study demonstrate that whole blood is a useful tissue for studying the metabolic syndrome and its underlying insulin resistance although the relationship between blood and skeletal muscle differs.

Background: Although the effect of the fat mass and obesity-associated (FTO) gene on adiposity is well established, there is a lack of evidence whether physical activity (PA) modifies the effect of FTO variants on obesity in Latino populations. Therefore, the purpose of this study was to examine PA influences and interactive effects between FTO variants and PA on measures of adiposity in Latinos.

Results: After controlling for age and sex, participants who did not engage in regular PA exhibited higher BMI, fat mass, HC, and WC with statistical significance (P < 0.001). Although significant associations between the three FTO genotypes and adiposity measures were found, none of the FTO genotype by PA interaction assessments revealed nominally significant associations. However, several of such interactive influences exhibited considerable trend towards association.

Conclusions: These data suggest that adiposity measures are associated with PA and FTO variants in Latinos, but the impact of their interactive influences on these obesity measures appear to be minimal. Future studies with large sample sizes may help to determine whether individuals with specific FTO variants exhibit differential responses to PA interventions.

Background: Obesity is a metabolic disease caused by environmental and genetic factors. However, the epigenetic mechanisms of obesity are incompletely understood. The aim of our study was to investigate the role of skeletal muscle DNA methylation in combination with transcriptomic changes in obesity.

Results: Muscle biopsies were obtained basally from lean (n = 12; BMI = 23.4 ± 0.7 kg/m[superscript 2]) and obese (n = 10; BMI = 32.9 ± 0.7 kg/m[superscript 2]) participants in combination with euglycemic-hyperinsulinemic clamps to assess insulin sensitivity. We performed reduced representation bisulfite sequencing (RRBS) next-generation methylation and microarray analyses on DNA and RNA isolated from vastus lateralis muscle biopsies. There were 13,130 differentially methylated cytosines (DMC; uncorrected P < 0.05) that were altered in the promoter and untranslated (5' and 3'UTR) regions in the obese versus lean analysis. Microarray analysis revealed 99 probes that were significantly (corrected P < 0.05) altered. Of these, 12 genes (encompassing 22 methylation sites) demonstrated a negative relationship between gene expression and DNA methylation. Specifically, sorbin and SH3 domain containing 3 (SORBS3) which codes for the adapter protein vinexin was significantly decreased in gene expression (fold change −1.9) and had nine DMCs that were significantly increased in methylation in obesity (methylation differences ranged from 5.0 to 24.4 %). Moreover, differentially methylated region (DMR) analysis identified a region in the 5'UTR (Chr.8:22,423,530–22,423,569) of SORBS3 that was increased in methylation by 11.2 % in the obese group. The negative relationship observed between DNA methylation and gene expression for SORBS3 was validated by a site-specific sequencing approach, pyrosequencing, and qRT-PCR. Additionally, we performed transcription factor binding analysis and identified a number of transcription factors whose binding to the differentially methylated sites or region may contribute to obesity.

Conclusions: These results demonstrate that obesity alters the epigenome through DNA methylation and highlights novel transcriptomic changes in SORBS3 in skeletal muscle.

Introduction: Decreased insulin sensitivity blunts the normal increase in gene expression from skeletal muscle after exercise. In addition, chronic inflammation decreases insulin sensitivity. Chronic kidney disease (CKD) is an inflammatory state. How CKD and, subsequently, kidney transplantation affects skeletal muscle gene expression after exercise are unknown.

Methods: Study cohort: non-diabetic male/female 4/1, age 52±2 years, with end-stage CKD who underwent successful kidney transplantation. The following were measured both pre-transplant and post-transplant and compared to normals: Inflammatory markers, euglycemic insulin clamp studies determine insulin sensitivity, and skeletal muscle biopsies performed before and within 30 minutes after an acute exercise protocol. Microarray analyses were performed on the skeletal muscle using the 4x44K Whole Human Genome Microarrays. Since nuclear factor of activated T cells (NFAT) plays an important role in T cell activation and calcineurin inhibitors are mainstay immunosuppression, calcineurin/NFAT pathway gene expression was compared at rest and after exercise. Log transformation was performed to prevent skewing of data and regression analyses comparing measures pre- and post-transplant performed.

Result: Markers of inflammation significantly improved post-transplantation. Insulin infusion raised glucose disposal slightly lower post-transplant compared to pre-transplant, but not significantly, thus concluding differences in insulin sensitivity were similar. The overall pattern of gene expression in response to exercise was reduced both pre-and post-transplant compared to healthy volunteers. Although significant changes were observed among NFAT/Calcineurin gene at rest and after exercise in normal cohort, there were no significant differences comparing NFAT/calcineurin pathway gene expression pre- and post-transplant.

Conclusions: Despite an improvement in serum inflammatory markers, no significant differences in glucose disposal were observed post-transplant. The reduced skeletal muscle gene expression, including NFAT/calcineurin gene expression, in response to a single bout of exercise was not improved post-transplant. This study suggests that the improvements in inflammatory mediators post-transplant are unrelated to changes of NFAT/calcineurin gene expression.

We have previously hypothesized a biological pathway of activity-dependent synaptic plasticity proteins that addresses the dual genetic and environmental contributions to schizophrenia. Accordingly, variations in the immediate early gene EGR3, and its target ARC, should influence schizophrenia susceptibility. We used a pooled Next-Generation Sequencing approach to identify variants across these genes in U.S. populations of European (EU) and African (AA) descent. Three EGR3 and one ARC SNP were selected and genotyped for validation, and three SNPs were tested for association in a replication cohort. In the EU group of 386 schizophrenia cases and 150 controls EGR3 SNP rs1877670 and ARC SNP rs35900184 showed significant associations (p = 0.0078 and p = 0.0275, respectively). In the AA group of 185 cases and 50 controls, only the ARC SNP revealed significant association (p = 0.0448). The ARC SNP did not show association in the Han Chinese (CH) population. However, combining the EU, AA, and CH groups revealed a highly significant association of ARC SNP rs35900184 (p = 2.353 x 10^-7; OR [95% CI] = 1.54 [1.310–1.820]). These findings support previously reported associations between EGR3 and schizophrenia. Moreover, this is the first report associating an ARC SNP with schizophrenia and supports recent large-scale GWAS findings implicating the ARC complex in schizophrenia risk. These results support the need for further investigation of the proposed pathway of environmentally responsive, synaptic plasticity-related, schizophrenia genes.

Background: Immunosignaturing is a new peptide microarray based technology for profiling of humoral immune responses. Despite new challenges, immunosignaturing gives us the opportunity to explore new and fundamentally different research questions. In addition to classifying samples based on disease status, the complex patterns and latent factors underlying immunosignatures, which we attempt to model, may have a diverse range of applications.

Methods: We investigate the utility of a number of statistical methods to determine model performance and address challenges inherent in analyzing immunosignatures. Some of these methods include exploratory and confirmatory factor analyses, classical significance testing, structural equation and mixture modeling.

Results: We demonstrate an ability to classify samples based on disease status and show that immunosignaturing is a very promising technology for screening and presymptomatic screening of disease. In addition, we are able to model complex patterns and latent factors underlying immunosignatures. These latent factors may serve as biomarkers for disease and may play a key role in a bioinformatic method for antibody discovery.

Conclusion: Based on this research, we lay out an analytic framework illustrating how immunosignatures may be useful as a general method for screening and presymptomatic screening of disease as well as antibody discovery.

Background: Glioblastoma is the most aggressive primary central nervous tumor and carries a very poor prognosis. Invasion precludes effective treatment and virtually assures tumor recurrence. In the current study, we applied analytical and bioinformatics approaches to identify a set of microRNAs (miRs) from several different human glioblastoma cell lines that exhibit significant differential expression between migratory (edge) and migration-restricted (core) cell populations. The hypothesis of the study is that differential expression of miRs provides an epigenetic mechanism to drive cell migration and invasion.

Results: Our research data comprise gene expression values for a set of 805 human miRs collected from matched pairs of migratory and migration-restricted cell populations from seven different glioblastoma cell lines. We identified 62 down-regulated and 2 up-regulated miRs that exhibit significant differential expression in the migratory (edge) cell population compared to matched migration-restricted (core) cells. We then conducted target prediction and pathway enrichment analysis with these miRs to investigate potential associated gene and pathway targets. Several miRs in the list appear to directly target apoptosis related genes. The analysis identifies a set of genes that are predicted by 3 different algorithms, further emphasizing the potential validity of these miRs to promote glioblastoma.

Conclusions: The results of this study identify a set of miRs with potential for decreased expression in invasive glioblastoma cells. The verification of these miRs and their associated targeted proteins provides new insights for further investigation into therapeutic interventions. The methodological approaches employed here could be applied to the study of other diseases to provide biomedical researchers and clinicians with increased opportunities for therapeutic interventions.

Background: The discovery of genetic associations is an important factor in the understanding of human illness to derive disease pathways. Identifying multiple interacting genetic mutations associated with disease remains challenging in studying the etiology of complex diseases. And although recently new single nucleotide polymorphisms (SNPs) at genes implicated in immune response, cholesterol/lipid metabolism, and cell membrane processes have been confirmed by genome-wide association studies (GWAS) to be associated with late-onset Alzheimer's disease (LOAD), a percentage of AD heritability continues to be unexplained. We try to find other genetic variants that may influence LOAD risk utilizing data mining methods.

Methods: Two different approaches were devised to select SNPs associated with LOAD in a publicly available GWAS data set consisting of three cohorts. In both approaches, single-locus analysis (logistic regression) was conducted to filter the data with a less conservative p-value than the Bonferroni threshold; this resulted in a subset of SNPs used next in multi-locus analysis (random forest (RF)). In the second approach, we took into account prior biological knowledge, and performed sample stratification and linkage disequilibrium (LD) in addition to logistic regression analysis to preselect loci to input into the RF classifier construction step.

Results: The first approach gave 199 SNPs mostly associated with genes in calcium signaling, cell adhesion, endocytosis, immune response, and synaptic function. These SNPs together with APOE and GAB2 SNPs formed a predictive subset for LOAD status with an average error of 9.8% using 10-fold cross validation (CV) in RF modeling. Nineteen variants in LD with ST5, TRPC1, ATG10, ANO3, NDUFA12, and NISCH respectively, genes linked directly or indirectly with neurobiology, were identified with the second approach. These variants were part of a model that included APOE and GAB2 SNPs to predict LOAD risk which produced a 10-fold CV average error of 17.5% in the classification modeling.

Conclusions: With the two proposed approaches, we identified a large subset of SNPs in genes mostly clustered around specific pathways/functions and a smaller set of SNPs, within or in proximity to five genes not previously reported, that may be relevant for the prediction/understanding of AD.

We study the so-called Descent, or [bar over Q], Equation for the null polygonal supersymmetric Wilson loop in the framework of the pentagon operator product expansion. To properly address this problem, one requires to restore the cyclicity of the loop broken by the choice of OPE channels. In the course of the study, we unravel a phenomenon of twist enhancement when passing to a cyclically shifted channel. Currently, we focus on the consistency of the all-order Descent Equation for the particular case relating the NMHV heptagon to MHV hexagon. We find that the equation establishes a relation between contributions of different twists and successfully verify it in perturbation theory making use of available bootstrap predictions for elementary pentagons.