ASU Scholarship Showcase

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: Malignant brain tumors affect people of all ages and are the second leading cause of cancer deaths in children. While current treatments are effective and improve survival, there remains a substantial need for more efficacious therapeutic modalities. The ketogenic diet (KD) - a high-fat, low-carbohydrate treatment for medically refractory epilepsy - has been suggested as an alternative strategy to inhibit tumor growth by altering intrinsic metabolism, especially by inducing glycopenia.

Methods: Here, we examined the effects of an experimental KD on a mouse model of glioma, and compared patterns of gene expression in tumors vs. normal brain from animals fed either a KD or a standard diet.

Results: Animals received intracranial injections of bioluminescent GL261-luc cells and tumor growth was followed in vivo. KD treatment significantly reduced the rate of tumor growth and prolonged survival. Further, the KD reduced reactive oxygen species (ROS) production in tumor cells. Gene expression profiling demonstrated that the KD induces an overall reversion to expression patterns seen in non-tumor specimens. Notably, genes involved in modulating ROS levels and oxidative stress were altered, including those encoding cyclooxygenase 2, glutathione peroxidases 3 and 7, and periredoxin 4.

Conclusions: Our data demonstrate that the KD improves survivability in our mouse model of glioma, and suggests that the mechanisms accounting for this protective effect likely involve complex alterations in cellular metabolism beyond simply a reduction in glucose.

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.

Background: Microarray image analysis processes scanned digital images of hybridized arrays to produce the input spot-level data for downstream analysis, so it can have a potentially large impact on those and subsequent analysis. Signal saturation is an optical effect that occurs when some pixel values for highly expressed genes or peptides exceed the upper detection threshold of the scanner software (2¹⁶ - 1 = 65, 535 for 16-bit images). In practice, spots with a sizable number of saturated pixels are often flagged and discarded. Alternatively, the saturated values are used without adjustments for estimating spot intensities. The resulting expression data tend to be biased downwards and can distort high-level analysis that relies on these data. Hence, it is crucial to effectively correct for signal saturation.

Results: We developed a flexible mixture model-based segmentation and spot intensity estimation procedure that accounts for saturated pixels by incorporating a censored component in the mixture model. As demonstrated with biological data and simulation, our method extends the dynamic range of expression data beyond the saturation threshold and is effective in correcting saturation-induced bias when the lost information is not tremendous. We further illustrate the impact of image processing on downstream classification, showing that the proposed method can increase diagnostic accuracy using data from a lymphoma cancer diagnosis study.

Conclusions: The presented method adjusts for signal saturation at the segmentation stage that identifies a pixel as part of the foreground, background or other. The cluster membership of a pixel can be altered versus treating saturated values as truly observed. Thus, the resulting spot intensity estimates may be more accurate than those obtained from existing methods that correct for saturation based on already segmented data. As a model-based segmentation method, our procedure is able to identify inner holes, fuzzy edges and blank spots that are common in microarray images. The approach is independent of microarray platform and applicable to both single- and dual-channel microarrays.

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.

One of the gravest dangers facing cancer patients is an extended symptom-free lull between tumor initiation and the first diagnosis. Detection of tumors is critical for effective intervention. Using the body’s immune system to detect and amplify tumor-specific signals may enable detection of cancer using an inexpensive immunoassay. Immunosignatures are one such assay: they provide a map of antibody interactions with random-sequence peptides. They enable detection of disease-specific patterns using classic train/test methods. However, to date, very little effort has gone into extracting information from the sequence of peptides that interact with disease-specific antibodies. Because it is difficult to represent all possible antigen peptides in a microarray format, we chose to synthesize only 330,000 peptides on a single immunosignature microarray. The 330,000 random-sequence peptides on the microarray represent 83% of all tetramers and 27% of all pentamers, creating an unbiased but substantial gap in the coverage of total sequence space. We therefore chose to examine many relatively short motifs from these random-sequence peptides. Time-variant analysis of recurrent subsequences provided a means to dissect amino acid sequences from the peptides while simultaneously retaining the antibody–peptide binding intensities. We first used a simple experiment in which monoclonal antibodies with known linear epitopes were exposed to these random-sequence peptides, and their binding intensities were used to create our algorithm. We then demonstrated the performance of the proposed algorithm by examining immunosignatures from patients with Glioblastoma multiformae (GBM), an aggressive form of brain cancer. Eight different frameshift targets were identified from the random-sequence peptides using this technique. If immune-reactive antigens can be identified using a relatively simple immune assay, it might enable a diagnostic test with sufficient sensitivity to detect tumors in a clinically useful way.

Background: High-throughput technologies such as DNA, RNA, protein, antibody and peptide microarrays are often used to examine differences across drug treatments, diseases, transgenic animals, and others. Typically one trains a classification system by gathering large amounts of probe-level data, selecting informative features, and classifies test samples using a small number of features. As new microarrays are invented, classification systems that worked well for other array types may not be ideal. Expression microarrays, arguably one of the most prevalent array types, have been used for years to help develop classification algorithms. Many biological assumptions are built into classifiers that were designed for these types of data. One of the more problematic is the assumption of independence, both at the probe level and again at the biological level. Probes for RNA transcripts are designed to bind single transcripts. At the biological level, many genes have dependencies across transcriptional pathways where co-regulation of transcriptional units may make many genes appear as being completely dependent. Thus, algorithms that perform well for gene expression data may not be suitable when other technologies with different binding characteristics exist. The immunosignaturing microarray is based on complex mixtures of antibodies binding to arrays of random sequence peptides. It relies on many-to-many binding of antibodies to the random sequence peptides. Each peptide can bind multiple antibodies and each antibody can bind multiple peptides. This technology has been shown to be highly reproducible and appears promising for diagnosing a variety of disease states. However, it is not clear what is the optimal classification algorithm for analyzing this new type of data.

Results: We characterized several classification algorithms to analyze immunosignaturing data. We selected several datasets that range from easy to difficult to classify, from simple monoclonal binding to complex binding patterns in asthma patients. We then classified the biological samples using 17 different classification algorithms. Using a wide variety of assessment criteria, we found ‘Naïve Bayes’ far more useful than other widely used methods due to its simplicity, robustness, speed and accuracy.

Conclusions: ‘Naïve Bayes’ algorithm appears to accommodate the complex patterns hidden within multilayered immunosignaturing microarray data due to its fundamental mathematical properties.

We have previously shown that the diversity of antibodies in an individual can be displayed on chips on which 130,000 peptides chosen from random sequence space have been synthesized. This immunosignature technology is unbiased in displaying antibody diversity relative to natural sequence space, and has been shown to have diagnostic and prognostic potential for a wide variety of diseases and vaccines. Here we show that a global measure such as Shannon’s entropy can be calculated for each immunosignature. The immune entropy was measured across a diverse set of 800 people and in 5 individuals over 3 months. The immune entropy is affected by some population characteristics and varies widely across individuals. We find that people with infections or breast cancer, generally have higher entropy values than non-diseased individuals. We propose that the immune entropy as measured from immunosignatures may be a simple method to monitor health in individuals and populations.

Honey bees as other insects rely on the innate immune system for protection against diseases. The innate immune system includes the circulating hemocytes (immune cells) that clear pathogens from hemolymph (blood) by phagocytosis, nodulation or encapsulation. Honey bee hemocyte numbers have been linked to hemolymph levels of vitellogenin. Vitellogenin is a multifunctional protein with immune-supportive functions identified in a range of species, including the honey bee. Hemocyte numbers can increase via mitosis, and this recruitment process can be important for immune system function and maintenance. Here, we tested if hemocyte mediated phagocytosis differs among the physiologically different honey bee worker castes (nurses, foragers and winter bees), and study possible interactions with vitellogenin and hemocyte recruitment. To this end, we adapted phagocytosis assays, which—together with confocal microscopy and flow cytometry—allow qualitative and quantitative assessment of hemocyte performance. We found that nurses are more efficient in phagocytic uptake than both foragers and winter bees. We detected vitellogenin within the hemocytes, and found that winter bees have the highest numbers of vitellogenin-positive hemocytes. Connections between phagocytosis, hemocyte-vitellogenin and mitosis were worker caste dependent. Our results demonstrate that the phagocytic performance of immune cells differs significantly between honey bee worker castes, and support increased immune competence in nurses as compared to forager bees. Our data, moreover, provides support for roles of vitellogenin in hemocyte activity.