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: 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.

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.

Immunosignaturing shows promise as a general approach to diagnosis. It has been shown to detect immunological signs of infection early during the course of disease and to distinguish Alzheimer’s disease from healthy controls. Here we test whether immunosignatures correspond to clinical classifications of disease using samples from people with brain tumors. Blood samples from patients undergoing craniotomies for therapeutically naïve brain tumors with diagnoses of astrocytoma (23 samples), Glioblastoma multiforme (22 samples), mixed oligodendroglioma/astrocytoma (16 samples), oligodendroglioma (18 samples), and 34 otherwise healthy controls were tested by immunosignature. Because samples were taken prior to adjuvant therapy, they are unlikely to be perturbed by non-cancer related affects. The immunosignaturing platform distinguished not only brain cancer from controls, but also pathologically important features about the tumor including type, grade, and the presence or absence of O⁶-methyl-guanine-DNA methyltransferase methylation promoter (MGMT), an important biomarker that predicts response to temozolomide in Glioblastoma multiformae patients.

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.

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: 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.

Epigenetic changes enable genomes to respond to changes in the environment, such as altered nutrition, activity, or social setting. Epigenetic modifications, thereby, provide a source of phenotypic plasticity in many species. The honey bee (Apis mellifera) uses nutritionally sensitive epigenetic control mechanisms in the development of the royal caste (queens) and the workers. The workers are functionally sterile females that can take on a range of distinct physiological and/or behavioral phenotypes in response to environmental changes. Honey bees have a wide repertoire of epigenetic mechanisms which, as in mammals, include cytosine methylation, hydroxymethylated cytosines, together with the enzymatic machinery responsible for these cytosine modifications. Current data suggests that honey bees provide an excellent system for studying the “social repertoire” of the epigenome. In this review, we elucidate what is known so far about the honey bee epigenome and its mechanisms. Our discussion includes what may distinguish honey bees from other model animals, how the epigenome can influence worker behavioral task separation, and how future studies can answer central questions about the role of the epigenome in social behavior.

Dynamical processes occurring on the edges in complex networks are relevant to a variety of real-world situations. Despite recent advances, a framework for edge controllability is still required for complex networks of arbitrary structure and interaction strength. Generalizing a previously introduced class of processes for edge dynamics, the switchboard dynamics, and exploit- ing the exact controllability theory, we develop a universal framework in which the controllability of any node is exclusively determined by its local weighted structure. This framework enables us to identify a unique set of critical nodes for control, to derive analytic formulas and articulate efficient algorithms to determine the exact upper and lower controllability bounds, and to evaluate strongly structural controllability of any given network. Applying our framework to a large number of model and real-world networks, we find that the interaction strength plays a more significant role in edge controllability than the network structure does, due to a vast range between the bounds determined mainly by the interaction strength. Moreover, transcriptional regulatory networks and electronic circuits are much more strongly structurally controllable (SSC) than other types of real-world networks, directed networks are more SSC than undirected networks, and sparse networks are typically more SSC than dense networks.

There is an increasing awareness that health care must move from post-symptomatic treatment to presymptomatic intervention. An ideal system would allow regular inexpensive monitoring of health status using circulating antibodies to report on health fluctuations. Recently, we demonstrated that peptide microarrays can do this through antibody signatures (immunosignatures). Unfortunately, printed microarrays are not scalable. Here we demonstrate a platform based on fabricating microarrays (~10 M peptides per slide, 330,000 peptides per assay) on silicon wafers using equipment common to semiconductor manufacturing. The potential of these microarrays for comprehensive health monitoring is verified through the simultaneous detection and classification of six different infectious diseases and six different cancers. Besides diagnostics, these high-density peptide chips have numerous other applications both in health care and elsewhere.