ASU Scholarship Showcase

Honeybee workers are essentially sterile female helpers that make up the majority of individuals in a colony. Workers display a marked change in physiology when they transition from in-nest tasks to foraging. Recent technological advances have made it possible to unravel the metabolic modifications associated with this transition. Previous studies have revealed extensive remodeling of brain, thorax, and hypopharyngeal gland biochemistry. However, data on changes in the abdomen is scarce. To narrow this gap we investigated the proteomic composition of abdominal tissue in the days typically preceding the onset of foraging in honeybee workers.

In order to get a broader representation of possible protein dynamics, we used workers of two genotypes with differences in the age at which they initiate foraging. This approach was combined with RNA interference-mediated downregulation of an insulin/insulin-like signaling component that is central to foraging behavior, the insulin receptor substrate (irs), and with measurements of glucose and lipid levels.
Our data provide new insight into the molecular underpinnings of phenotypic plasticity in the honeybee, invoke parallels with vertebrate metabolism, and support an integrated and irs-dependent association of carbohydrate and lipid metabolism with the transition from in-nest tasks to foraging.

Background: Gene bodies are the most evolutionarily conserved targets of DNA methylation in eukaryotes. However, the regulatory functions of gene body DNA methylation remain largely unknown. DNA methylation in insects appears to be primarily confined to exons. Two recent studies in Apis mellifera (honeybee) and Nasonia vitripennis (jewel wasp) analyzed transcription and DNA methylation data for one gene in each species to demonstrate that exon-specific DNA methylation may be associated with alternative splicing events. In this study we investigated the relationship between DNA methylation, alternative splicing, and cross-species gene conservation on a genome-wide scale using genome-wide transcription and DNA methylation data.

Results: We generated RNA deep sequencing data (RNA-seq) to measure genome-wide mRNA expression at the exon- and gene-level. We produced a de novo transcriptome from this RNA-seq data and computationally predicted splice variants for the honeybee genome. We found that exons that are included in transcription are higher methylated than exons that are skipped during transcription. We detected enrichment for alternative splicing among methylated genes compared to unmethylated genes using fisher’s exact test. We performed a statistical analysis to reveal that the presence of DNA methylation or alternative splicing are both factors associated with a longer gene length and a greater number of exons in genes. In concordance with this observation, a conservation analysis using BLAST revealed that each of these factors is also associated with higher cross-species gene conservation.

Conclusions: This study constitutes the first genome-wide analysis exhibiting a positive relationship between exon-level DNA methylation and mRNA expression in the honeybee. Our finding that methylated genes are enriched for alternative splicing suggests that, in invertebrates, exon-level DNA methylation may play a role in the construction of splice variants by positively influencing exon inclusion during transcription. The results from our cross-species homology analysis suggest that DNA methylation and alternative splicing are genetic mechanisms whose utilization could contribute to a longer gene length and a slower rate of gene evolution.

Insect immune systems can recognize specific pathogens and prime offspring immunity. High specificity of immune priming can be achieved when insect females transfer immune elicitors into developing oocytes. The molecular mechanism behind this transfer has been a mystery. Here, we establish that the egg-yolk protein vitellogenin is the carrier of immune elicitors. Using the honey bee, Apis mellifera, model system, we demonstrate with microscopy and western blotting that vitellogenin binds to bacteria, both Paenibacillus larvae – the gram-positive bacterium causing American foulbrood disease – and to Escherichia coli that represents gram-negative bacteria. Next, we verify that vitellogenin binds to pathogen-associated molecular patterns; lipopolysaccharide, peptidoglycan and zymosan, using surface plasmon resonance. We document that vitellogenin is required for transport of cell-wall pieces of E. coli into eggs by imaging tissue sections. These experiments identify vitellogenin, which is distributed widely in oviparous species, as the carrier of immune-priming signals. This work reveals a molecular explanation for trans-generational immunity in insects and a previously undescribed role for vitellogenin.

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.

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: Phosphatase and TENsin (PTEN) homolog is a negative regulator that takes part in IIS (insulin/insulin-like signaling) and Egfr (epidermal growth factor receptor) activation in Drosophila melanogaster. IIS and Egfr signaling events are also involved in the developmental process of queen and worker differentiation in honey bees (Apis mellifera). Here, we characterized the bee PTEN gene homologue for the first time and begin to explore its potential function during bee development and adult life.

Results: Honey bee PTEN is alternatively spliced, resulting in three splice variants. Next, we show that the expression of PTEN can be down-regulated by RNA interference (RNAi) in the larval stage, when female caste fate is determined. Relative to controls, we observed that RNAi efficacy is dependent on the amount of PTEN dsRNA that is delivered to larvae. For larvae fed queen or worker diets containing a high amount of PTEN dsRNA, PTEN knockdown was significant at a whole-body level but lethal. A lower dosage did not result in a significant gene down-regulation. Finally, we compared same-aged adult workers with different behavior: nursing vs. foraging. We show that between nurses and foragers, PTEN isoforms were differentially expressed within brain, ovary and fat body tissues. All isoforms were expressed at higher levels in the brain and ovaries of the foragers. In fat body, isoform B was expressed at higher level in the nurse bees.

Conclusion: Our results suggest that PTEN plays a central role during growth and development in queen- and worker-destined honey bees. In adult workers, moreover, tissue-specific patterns of PTEN isoform expression are correlated with differences in complex division of labor between same-aged individuals. Therefore, we propose that knowledge on the roles of IIS and Egfr activity in developmental and behavioral control may increase through studies of how PTEN functions can impact bee social phenotypes.

In vitro rearing is an important and useful tool for honey bee (Apis mellifera L.) studies. However, it often results in intercastes between queens and workers, which are normally are not seen in hive-reared bees, except when larvae older than three days are grafted for queen rearing. Morphological classification (queen versus worker or intercastes) of bees produced by this method can be subjective and generally depends on size differences. Here, we propose an alternative method for caste classification of female honey bees reared in vitro, based on weight at emergence, ovariole number, spermatheca size and size and shape, and features of the head, mandible and basitarsus. Morphological measurements were made with both traditional morphometric and geometric morphometrics techniques. The classifications were performed by principal component analysis, using naturally developed queens and workers as controls. First, the analysis included all the characters. Subsequently, a new analysis was made without the information about ovariole number and spermatheca size. Geometric morphometrics was less dependent on ovariole number and spermatheca information for caste and intercaste identification. This is useful, since acquiring information concerning these reproductive structures requires time-consuming dissection and they are not accessible when abdomens have been removed for molecular assays or in dried specimens. Additionally, geometric morphometrics divided intercastes into more discrete phenotype subsets. We conclude that morphometric geometrics are superior to traditional morphometrics techniques for identification and classification of honey bee castes and intermediates.

Background: Juvenile hormone (JH) has been demonstrated to control adult lifespan in a number of non-model insects where surgical removal of the corpora allata eliminates the hormone’s source. In contrast, little is known about how juvenile hormone affects adult Drosophila melanogaster. Previous work suggests that insulin signaling may modulate Drosophila aging in part through its impact on juvenile hormone titer, but no data yet address whether reduction of juvenile hormone is sufficient to control Drosophila life span. Here we adapt a genetic approach to knock out the corpora allata in adult Drosophila melanogaster and characterize adult life history phenotypes produced by reduction of juvenile hormone. With this system we test potential explanations for how juvenile hormone modulates aging.

Results: A tissue specific driver inducing an inhibitor of a protein phosphatase was used to ablate the corpora allata while permitting normal development of adult flies. Corpora allata knockout adults had greatly reduced fecundity, inhibited oogenesis, impaired adult fat body development and extended lifespan. Treating these adults with the juvenile hormone analog methoprene restored all traits toward wildtype. Knockout females remained relatively long-lived even when crossed into a genotype that blocked all egg production. Dietary restriction further extended the lifespan of knockout females. In an analysis of expression profiles of knockout females in fertile and sterile backgrounds, about 100 genes changed in response to loss of juvenile hormone independent of reproductive state.

Conclusions: Reduced juvenile hormone alone is sufficient to extend the lifespan of Drosophila melanogaster. Reduced juvenile hormone limits reproduction by inhibiting the production of yolked eggs, and this may arise because juvenile hormone is required for the post-eclosion development of the vitellogenin-producing adult fat body. Our data do not support a mechanism for juvenile hormone control of longevity simply based on reducing the physiological costs of egg production. Nor does the longevity benefit appear to function through mechanisms by which dietary restriction extends longevity. We identify transcripts that change in response to juvenile hormone independent of reproductive state and suggest these represent somatically expressed genes that could modulate how juvenile hormone controls persistence and longevity.

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.

Recent advancements in genomics provide new tools for evolutionary ecological research. The paper wasp genus Polistes is a model for social insect evolution and behavioral ecology. We developed RNA interference (RNAi)-mediated gene silencing to explore proposed connections between expression of hexameric storage proteins and worker vs. gyne (potential future foundress) castes in naturally-founded colonies of P. metricus. We extended four fragments of putative hexamerin-encoding P. metricus transcripts acquired from a previous study and fully sequenced a gene that encodes Hexamerin 2, one of two proposed hexameric storage proteins of P. metricus. MALDI-TOF/TOF, LC-MSMS, deglycosylation, and detection of phosphorylation assays showed that the two putative hexamerins diverge in peptide sequence and biochemistry. We targeted the hexamerin 2 gene in 5th (last)-instar larvae by feeding RNAi-inducing double-stranded hexamerin 2 RNA directly to larvae in naturally-founded colonies in the field. Larval development and adult traits were not significantly altered in hexamerin 2 knockdowns, but there were suggestive trends toward increased developmental time and less developed ovaries, which are gyne characteristics. By demonstrating how data acquisition from 454/Roche pyrosequencing can be combined with biochemical and proteomics assays and how RNAi can be deployed successfully in field experiments on Polistes, our results pave the way for functional genomic research that can contribute significantly to learning the interactions of environment, development, and the roles they play in paper wasp evolution and behavioral ecology.