Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Latest Articles
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • Archive
  • Front Matter
  • News
    • For the Press
    • Highlights from Latest Articles
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Purpose and Scope
    • Editorial and Journal Policies
    • Submission Procedures
    • For Reviewers
    • Author FAQ
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Latest Articles
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • Archive
  • Front Matter
  • News
    • For the Press
    • Highlights from Latest Articles
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Purpose and Scope
    • Editorial and Journal Policies
    • Submission Procedures
    • For Reviewers
    • Author FAQ

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology

Changes in transcript abundance relating to colony collapse disorder in honey bees (Apis mellifera)

Reed M. Johnson, Jay D. Evans, Gene E. Robinson, and May R. Berenbaum
PNAS September 1, 2009 106 (35) 14790-14795; https://doi.org/10.1073/pnas.0906970106
Reed M. Johnson
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jay D. Evans
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gene E. Robinson
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
May R. Berenbaum
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Contributed by May R. Berenbaum, July 14, 2009 (received for review February 18, 2009)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Abstract

Colony collapse disorder (CCD) is a mysterious disappearance of honey bees that has beset beekeepers in the United States since late 2006. Pathogens and other environmental stresses, including pesticides, have been linked to CCD, but a causal relationship has not yet been demonstrated. Because the gut acts as a primary interface between the honey bee and its environment as a site of entry for pathogens and toxins, we used whole-genome microarrays to compare gene expression between guts of bees from CCD colonies originating on both the east and west coasts of the United States and guts of bees from healthy colonies sampled before the emergence of CCD. Considerable variation in gene expression was associated with the geographical origin of bees, but a consensus list of 65 transcripts was identified as potential markers for CCD status. Overall, elevated expression of pesticide response genes was not observed. Genes involved in immune response showed no clear trend in expression pattern despite the increased prevalence of viruses and other pathogens in CCD colonies. Microarray analysis revealed unusual ribosomal RNA fragments that were conspicuously more abundant in the guts of CCD bees. The presence of these fragments may be a possible consequence of picorna-like viral infection, including deformed wing virus and Israeli acute paralysis virus, and may be related to arrested translation. Ribosomal fragment abundance and presence of multiple viruses may prove to be useful diagnostic markers for colonies afflicted with CCD.

  • microarray
  • picorna-like virus
  • ribosomal RNA

As the premier managed pollinator in the United States, the western honey bee, Apis mellifera, contributes more than $14 billion to agriculture annually (1). Beginning in fall 2006, the American apiculture industry experienced catastrophic losses of unknown origin. The phenomenon, called colony collapse disorder (CCD), was identified by a set of distinctive characteristics, including the absence of dead bees in or near the colony and the presence of abundant brood, honey, and pollen despite vastly reduced numbers of adult workers (2). Losses were estimated at 23% over the winter of 2006–2007 (3) and at 36% over the winter of 2007–2008 (4).

Using metagenomics, Cox-Foster et al. (3) compared the microbial flora of honey bees in hives diagnosed with CCD and in ostensibly healthy hives, evaluating pathogens with respect to their association with diagnosed CCD. The most predictive pathogen was Israeli acute paralysis virus (IAPV), a picorna-like virus (Dicistroviridae) hitherto unreported in the United States. However, a later study found IAPV in U.S. bees before the appearance of CCD (5), discounting an exclusive causal relationship. The metagenomic analysis identified other pathogens associated with CCD bees, including the picorna-like viruses Kashmir bee virus (KBV) and deformed wing virus (DWV), and 2 species in the microsporidian genus Nosema (3), one of which, N. ceranae, was subsequently linked to collapses of colonies in Spain (6), although this may be unrelated to CCD in the United States. A high prevalence of multiple pathogens in CCD bees suggests that a compromised immune response may be integral to CCD.

Pesticides also have been suspected to play a role in CCD. Sublethal exposures leading to behavioral disruptions consistent with the failure of foragers to return to their hives have been associated with the neonicotinoid pesticides (7). The long-term use of combinations of in-hive pesticides for control of honey bee parasites also may have contributed to otherwise unexplained bee mortality (8).

To differentiate among possible explanations for CCD, we used whole-genome microarray analysis, comparing gut gene expression in adult worker bees from healthy and CCD colonies. We assayed the gut because it is the principal site of pesticide detoxification and an integral component in the immune defense against pathogens in A. mellifera. Samples from colonies varying in CCD severity were collected on the U.S. east and west coasts in the winter of 2006–2007 and compared with healthy (“historical”) controls collected in 2004 and 2005.

Results

Gene Expression Differences.

The microarray contained oligonucleotide probes representing 9,867 different genes (with duplicate spots), based on gene predictions and annotation from the honey bee genome sequencing project; in addition, 2,729 probes specific for ESTs matching no existing annotation were represented on the array (9–11), including 21 probes derived from expressed sequence tags (ESTs) specific for rRNA [supporting information (SI) Table S1].

Of the 6,777 probes with expression above background level, 1,305 probes showed significant differences in expression in at least 1 of the 5 contrasts [false discovery rate (FDR), P < .01; and fold-change >2) (Fig. 1]. The west coast versus historical comparison alone generated a list of 948 differentially expressed probes; however, 668 of these were not differentially expressed in any of the other comparisons.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Experimental design for microarray comparisons of CCD and healthy bees. East coast CCD includes bees collected in Florida and Pennsylvania, and west coast CCD includes bees from California. Guts from bees remaining in colonies that were classified as mild or severe were compared with each other and together, through a reference gut RNA sample, with guts of healthy historical bees collected in 2004 and 2005. A total of 22 microarrays were hybridized.

There was little concordance in the lists of differentially expressed genes between the east and west coast comparison of guts taken from “severe” and “mild” colonies. Just 54 probes shared differential expression between the east and west coast severe versus mild comparisons, with 191 probes showing unique differential expression. That so many genes are differentially transcribed in the severe versus mild comparisons suggests these classifications may not be equivalent across geographical samplings.

Reflecting the strong differences in gene expression, hierarchical cluster analysis, using expression values for all probes, grouped the samples broadly by health status and geography (Fig. 2). All healthy bees clustered together, as did all CCD bees, which were further clustered according to geography, with west coast California bees forming one group and east coast Pennsylvania and Florida bees clustering together in another group. Differences in relative colony health at the time of collection, either mild or severe CCD, were not reflected in the cluster analysis.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Heatmap showing expression values for 65 probes demonstrating differential expression in 3 contrasts: east coast CCD versus historical, west coast CCD versus historical, and east coast severe versus east coast mild [P < .01 (FDR); fold-change >2]. Yellow indicates increased transcript abundance and red indicates decreased transcript abundance relative to the mean for all colonies. Probes lacking a description correspond to EST sequences for which no matching gene has been found. More information about the probes is provided in Table S5. CA, California; PA, Pennsylvania; FL, Florida. Numbers correspond to colony number. S, severe CCD; M, mild CCD. Colonywise clustering was performed using expression values for all 6,777 probes. Euclidean distances were calculated and clustered using the “complete” method.

Although geographic variation among the east and west coast CCD samples is evident, if CCD is a distinct phenomenon, then it should be possible to identify a list of CCD-related genes that show differential expression in both east and west coast CCD bees relative to healthy bees. But because these comparisons rely on the historical samples, which may carry their own geographic biases, we further winnowed the gene list by using only genes differentially expressed in the east coast severe versus mild comparison (Fig. 3). This comparison was informative; expression of 65 transcripts was consistently up or down when comparing relatively sick bees (CCD or severe) and relatively healthy bees (historically healthy or mild) (Fig. 2); 29 probes on this list correspond to annotated genes, 35 probes correspond to ESTs, and 1 probe corresponds to a pathogen, DWV.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Venn diagram showing the number of transcripts up-regulated and down-regulated in 3 comparisons: west coast (California) CCD versus historical, east coast (Florida and Pennsylvania) CCD versus historical, and east coast severe CCD versus mild CCD.

Gene Ontology Functional Analyses of Gene Expression Differences.

Genes differentially expressed between groups were categorized according to Gene Ontology (GO) (12). FlyBase orthologs of bee genes that were differentially expressed were used for GO enrichment analyses (Table 1). In the cellular component ontology, genes associated with lipid particle (GO:0005811), cytosolic small ribosomal subunit (GO:0022627), and mitochondrial respiratory chain (GO:0005746) generally demonstrated reduced expression in west coast bee guts relative to historical samples.

View this table:
  • View inline
  • View popup
Table 1.

GO analyses for differentially expressed genes in 3 contrasts

In the molecular function ontology, genes related to transcription factor activity (GO:0003700) were more highly expressed in both east coast and west coast CCD bees and in severe CCD bees. West coast bees also showed reduced expression of both heme binding (GO:0020037) and carboxylesterase activity (GO:0004091) categorizations, which include cytochrome P450 monooxygenases (P450s) and carboxylesterases (COEs), the principal xenobiotic metabolizing enzymes. One P450, CYP4G11, was expressed at a lower level in all 3 contrasts tested.

In the biological process ontology, genes related to development (GO:0016348 and GO:0007509) and the peptidoglycan catabolic process (GO:0009253) were overrepresented but did not demonstrate consistent directional differences. The peptidoglycan recognition proteins (PGRPs) contribute to immunity by detecting bacteria (13); PGRP-S2 was up-regulated and PGRP-S3 was down-regulated in CCD bees. Transcripts encoding sallimus, a large complex protein that mediates muscle elasticity (14), were more abundant in west coast bees. This protein also may contribute to immunity, because it contains Ig repeats and its transcripts are overexpressed in Anopheles gambiae mosquitoes after bacterial infection (15).

The gene mblk-1 was more highly expressed in CCD bees; the D. melanogaster ortholog of this gene plays a role in hormone-triggered cell death (16) and larval midgut histolysis (17). Thus, overexpression of mblk-1 may indicate apoptosis in the guts of CCD bees. However, of the 39 apoptosis-related (GO:0006915) genes on the microarray, only mblk-1 and GB14659, which is similar to apoptotic peptidase activating factor 1, were more highly expressed in the guts of CCD bees.

Pathogens in CCD Bees.

Twenty-two probes specific for 8 bee pathogens were present on the microarray [KBV, DWV, black queen cell virus (BQCV), acute bee paralysis virus, sacbrood virus, Ascophaera apis, Nosema apis, and Paenibacillus larvae] (Fig. S2), providing an opportunity to survey the samples for pathogens. In some cases, the number of colonies sampled was small with respect to standard pathogen surveys, so the results are best interpreted in a qualitative rather than quantitative manner. Three pathogens were detectable at different levels among the sampled bees. Chalkbrood (A. apis) RNA was more abundant in both east coast CCD bees and in severe CCD bees. BQCV RNA was more abundant in CCD bees, as was DWV, but DWV also was more abundant in east coast severe CCD bees. The DWV probe on the microarray also may have detected the closely related Kakugo virus or Varroa destructor virus 1 (18), with which it shares 91% identity. A probe specific for IAPV (19) was not included on the array and likely is not detected by the probe for KBV, with which it shares just 60% identity. The presence of Nosema, IAPV, and other viruses was surveyed with quantitative PCR (qPCR), indicating a generally higher viral load in CCD colonies (Table 3).

Validation of Microarray Results with qPCR.

Expression changes of 8 transcripts were verified by qPCR. Five of 8 transcripts had similar expression profiles when measured using the same RNA samples analyzed by both methods (Table S4). Biological validation of microarray results was performed using qPCR on gut samples from healthy and CCD bees collected in California in 2007, a year after the samples used for microarray analysis were obtained. Differential expression of 3 of 8 transcripts was confirmed: CCD-associated up-regulation of the ESTs QW33 and jdeC15, corresponding to 28S and 5.8S rRNA subunits. DWV also was more abundant in CCD gut samples (Table 2).

View this table:
  • View inline
  • View popup
Table 2.

Comparison of fold-change expression differences as measured by microarray and qPCR

Further qPCR validation of markers was performed with abdomens of bees from 147 colonies in healthy and CCD apiaries in 6 states. Only an rRNA-coding EST, QW33, was differentially expressed between CCD and healthy abdomens (Table 2). These CCD bees also were more likely to harbor both a larger number of viruses (as detected by these methods) and nonviral pathogens, as in the original samples (Table 3).

View this table:
  • View inline
  • View popup
Table 3.

Proportion of colonies with detectable levels of viral and microsporidial pathogens, as measured by qPCR

Discussion

No simple explanations for the cause of CCD emerge from the microarray analysis. Expression of detoxification and immune gene transcripts, which would indicate toxins or disease as the cause of CCD, was largely unchanged. However, considerable geographic variation existed among CCD bees, with west coast bees seemingly more severely affected. Transcripts of genes related to basic cellular processes involving ribosomal and mitochondrial function were generally less abundant in west coast bees.

P450s and COEs play multiple metabolic roles and are particularly important in the detoxification of natural and synthetic toxins in bees and other insects (20). However, the few genes in these superfamilies differentially transcribed with CCD likely serve functions other than detoxification. Expression of CYP4G11 was elevated in CCD bees, but the function of this P450 remains uncharacterized. Its ortholog in D. melanogaster, CYP4G15, is expressed only in the central nervous system, where it may be associated with ecdysteroid metabolism (21). The single COE overexpressed in CCD bees, GB10854, differs from other insect COEs in terms of its catalytic site, such that it may be incapable of detoxificative carboxylester hydrolysis (22).

Although 88 known immune-related genes were included on the microarray, these also were mostly unchanged in expression in CCD bees. Given the association between pathogens and CCD (3), the finding that so few immune genes were differentially expressed in CCD bees was unexpected. Insect guts can produce antimicrobial peptides, including the apidaecins in honey bees (23), yet CCD-related changes in transcription of these genes were not detected. This lack of response does not rule out viral infection, however, because antimicrobial proteins are ineffective against viruses. In D. melanogaster, fungal or bacterial infection induces antimicrobial peptides, but flies infected with Drosophila C virus (DCV), another picorna-like virus, did not produce more antimicrobial peptides (24) and demonstrated elevated transcription of only a small number of immune genes relative to fungal and bacterial infection (25).

Many of the transcripts differentially expressed in CCD bees were detected by probes based on ESTs matching no existing annotation (10, 11). Of the 2,729 geneless ESTs on the array, 28 were differentially expressed in CCD bees. Of all of the ESTs on which probes were based, 13 were discovered to contain rRNA sequence, and 5 of these were more highly expressed in CCD bee guts. These rRNA-coding ESTs were not identified as rRNA previously, because the A. mellifera rRNA sequence is not available in online repositories, and identification of these probes was accomplished only by manually transcribing the published rRNA sequence from printed text (26). Had these ESTs been previously identified as rRNA, the associated probes likely would have been excluded from the microarray.

The presence of rRNA on these microarrays is surprising, because both the EST projects on which the probes were based and the microarray hybridization rely on a polyadenylated 3′ tail for the initial priming step in the reverse-transcription reaction. Ribosomal RNA is transcribed by RNA polymerase I, and the resulting transcript is not polyadenylated; as such, oligo(dT)-based reverse-transcription reactions should not reverse-transcribe rRNA. Yet transcripts of rRNA have long appeared in EST projects, and their reverse transcription has been attributed to poly(A)-rich internal sequences hybridizing to the oligo(dT) or to genomic DNA contamination (27). This differential expression of polyadenylated rRNA transcripts in insects currently defies technical explanation; differential expression of these transcripts associated with age and caste in a termite (Reticulotermes flavipes) EST project (28) suggests an underlying biological explanation.

One possible explanation for the presence of poly(A) rRNA sequences in bees is that they are degradation intermediates. Polyadenylation is a well-known marker of RNA degradation in bacteria and organelles (29), and poly(A) addition is known to contribute to eukaryotic rRNA degradation (30). Polyadenylated 25S rRNA comprises ≈0.02% of all of the 25S rRNA transcripts in fission yeast Schizosaccharomyces pombe (31). Another explanation for these sequences is the “ribosomal filter hypothesis,” in which conventional mRNA transcripts include tens to hundreds of nucleotides similar to rRNA sequence in either the sense or antisense direction (32). Transcripts of genes containing rRNA sequence can bind to ribosomal subunits and alter translation efficiency (33). Short sections of rRNA function as an internal ribosome entry site (IRES) to interact with the ribosome and alter translation. Indeed, functions have been assigned to some mRNAs containing large sections of rRNA sequence (34), although none of the predicted genes annotated from the honey bee genome contains sufficient rRNA sequence to hybridize to any of the rRNA probes on the array (Table S1). A related hypothesis posits that long sequences of rRNA-like transcript may aid in protein folding, taking on part of the enzymatic activity of rRNA (35).

The picorna-like virus DWV was more abundant in CCD bees analyzed in the microarray experiment and a broader sampling of bee colonies revealed that CCD bees carry a larger number of different picorna-like viruses (Table 3). Picorna-like viruses, including IAPV and KBV, that have been associated with CCD (3) initiate transcription of viral proteins by the ribosome through an IRES rather than through the 5′-methylation cap that initiates translation of most mRNA (36). The IRES sequences from both DWV and Varroa destructor virus 1 effectively enhance translation of a reporter gene (36). Picorna-like viruses in mammals also halt translation of host mRNA through cleavage of the translation initiation factor eIF4G and the poly(A)-binding protein (37), leaving the ribosomes incapable of binding host mRNA and giving the viral RNA little competition for its IRES-mediated binding to host ribosomes (38).

While there is no direct evidence that elevated poly(A) rRNA is a result of picorna-like viral infection, one consequence of the viral disabling of ribosomal function may well be increased ribosomal degradation. Picornavirus infection in mammals both reduces protein production and causes strings of translating polyribosomes to break down, and these idle ribosomal subunits may be more susceptible to degradation (39). Ribosomal degradation also could occur in CCD bees independent of viral infection; for example, starving bacteria demonstrate elevated ribosomal degradation (39). Regardless of the cause, the abundance of poly(A) rRNA may be a useful diagnostic marker for determining CCD status, although the presence of poly(A) rRNA does not in itself indicate CCD, because healthy bees, including those used to generate EST libraries, contain these fragments.

Although gene transcript analysis did not clearly identify a specific cause for CCD, our study documents several patterns suggestive of a causal mechanism. The reduced protein synthetic capabilities that would accompany ribosomal hijacking by multiple picorna-like viruses would leave bees unable to respond to additional stresses from pesticides, nutrition, or pathogens. Although any interpretation of the presence of these rRNA fragments is speculative, the reported interaction between bee picorna-like viruses and rRNA is suggestive of a possible root cause of CCD.

To establish a causal relationship, the quantitative association between multiple picorna-like virus infections and polyadenylated rRNA fragment abundance merits further exploration. In addition, the consequences of viral infection and CCD on the function of ribosomes should be explored through assays of translational efficiency. Because of the potential for translational interference, studies on immune suppression should focus on bioassays or protein abundance rather than on immune gene transcripts.

Our results also indicate an unappreciated variation in gene expression patterns and pathogen loads with geography. This variation provides insight into the different stresses facing bees and clearly demonstrates that diagnostic surveys must sample extensively across numerous bee populations. Colony surveillance via assay of rRNA-like transcript abundance may provide an earlier indication of CCD status than has hitherto been available and allow beekeepers to take actions to reduce losses.

Materials and Methods

Biological Material and RNA Extraction.

Remaining adult bees in colonies diagnosed with CCD at apiaries in Florida, California, and Pennsylvania were collected during winter 2006–2007. CCD colony health was scored at the time of collection as either severe or mild, depending on the apparent strength of the colony. Historical bees, collected before the appearance of CCD and hence ostensibly healthy, were collected in 2004 and 2005 from colonies set up on new equipment with no miticide treatments in apiaries of Pennsylvania State University near State College, Pennsylvania (Table S3). Bees for the reference (ostensibly healthy) sample were collected from 5 hives near Urbana, Illinois in July 2007. To date, there have been no confirmed cases of CCD in central Illinois. All bees were immediately frozen on dry ice and kept at −80 °C until dissection.

Dissection, RNA extraction, microarray, statistical and qPCR validation analyses are described in SI Text. Pooled samples of RNA from 6 guts collected from each CCD colony were hybridized to 2 arrays, incorporating a dye swap (Fig. S1). Two separate pools of RNA from 6 guts each were created for both historical samplings, and each was hybridized against the reference on a single microarray. Experiments were designed to meet Minimum Information About a Microarray Experiment standards, and all microarray data obtained in these studies were deposited at ArrayExpress [www.ebi.ac.uk/arrayexpress (accession no. E-MEXP-2292)].

Acknowledgments

We thank D. vanEngelsdorp, J. Pettis, and D. Cox-Foster for sample collection and colony health scores, and D. Lopez and T. Newman for laboratory assistance. We also thank I. Lipkin, J. Willis, and J. Gillespie for an insightful and constructive review. This work was supported by U.S. Department of Agriculture Grants USDAG 2007–37610-18499 (to S. Ratcliffe, G.E.R., and M.B.), USDAG 2008–3532-18831 (to M.R.B.), and USDA-NRI 2004–35604-14277 (to G.E.R.).

Footnotes

  • 2To whom correspondence should be addressed. E-mail: maybe{at}illinois.edu
  • Author contributions: R.M.J., J.D.E., G.E.R., and M.R.B. designed research; R.M.J. and J.D.E. performed research; J.D.E. and G.E.R. contributed new reagents/analytic tools; R.M.J., J.D.E., G.E.R., and M.R.B. analyzed data; and R.M.J., J.D.E., G.E.R., and M.R.B. wrote the paper.

  • The authors declare no conflict of interest.

  • Data deposition: The data reported in this paper have been deposited in the ArrayExpress database, www.ebi.ac.uk/arrayexpress (accession no. E-MEXP-2292).

  • This article contains supporting information online at www.pnas.org/cgi/content/full/0906970106/DCSupplemental.

References

  1. ↵
    1. National Academy of Sciences
    (2007) Status of Pollinators in North America (National Academies Press, Washington, DC).
  2. ↵
    1. vanEngelsdorp D,
    2. Underwood R,
    3. Caron D,
    4. Hayes J
    (2007) An estimate of managed colony losses in the winter of 2006–2007: A report commissioned by the apiary inspectors of America. Am Bee J 147:599–603.
    OpenUrl
  3. ↵
    1. Cox-Foster DL,
    2. et al.
    (2007) A metagenomic survey of microbes in honey bee colony collapse disorder. Science 318:283–287.
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. vanEngelsdorp D,
    2. Hayes J, Jr,
    3. Underwood RM,
    4. Pettis J
    (2008) A survey of honey bee colony losses in the U.S., fall 2007 to spring 2008. PLoS One 3:e4071.
    OpenUrlCrossRefPubMed
  5. ↵
    1. The Honeybee Genome Sequencing Consortium
    (2006) Insights into social insects from the genome of the honeybee, Apis mellifera. Nature 443:931–949.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Chen YP,
    2. Evans JD
    (2007) Historical presence of Israeli acute paralysis virus in the United States. Am Bee J 147:1027–1028.
    OpenUrl
  7. ↵
    1. Higes M,
    2. et al.
    (2009) Honeybee colony collapse due to Nosema ceranae in professional apiaries. Environ Microbiol Rep 1:110–113.
    OpenUrlCrossRef
  8. ↵
    1. Decourtye A,
    2. Devillers J,
    3. Cluzeau S,
    4. Charreton M,
    5. Pham-Delègue M
    (2004) Effects of imidacloprid and deltamethrin on associative learning in honeybees under semi-field and laboratory conditions. Ecotoxicol Environ Saf 57:410–419.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Johnson RM,
    2. Pollock HS,
    3. Berenbaum MR
    (2009) Synergistic interactions between in-hive miticides in Apis mellifera. J Econ Entomol 102:474–479.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Evans JD,
    2. Wheeler DE
    (2000) Expression profiles during honeybee caste determination. Genome Biol 2:research0001.1–0001.6.
  11. ↵
    1. Whitfield CW,
    2. et al.
    (2002) Annotated expressed sequence tags and cDNA microarrays for studies of brain and behavior in the honey bee. Genome Res 12:555–566.
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Ashburner M,
    2. et al.
    (2000) Gene Ontology: Tool for the unification of biology. Nat Genet 25:25–29.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Werner T,
    2. Borge-Renberg K,
    3. Mellroth P,
    4. Steiner H,
    5. Hultmark D
    (2003) Functional diversity of the Drosophila PGRP-LC gene cluster in the response to lipopolysaccharide and peptidoglycan. J Biol Chem 278:26319–26322.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Burkart C,
    2. et al.
    (2007) Modular proteins from the Drosophila sallimus (sls) gene and their expression in muscles with different extensibility. J Mol Biol 367:953–969.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Garver LS,
    2. Xi Z,
    3. Dimopoulos G
    (2008) Immunoglobulin superfamily members play an important role in the mosquito immune system. Dev Comp Immunol 32:519–531.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Lee T,
    2. Marticke S,
    3. Sung C,
    4. Robinow S,
    5. Luo L
    (2000) Cell-autonomous requirement of the USP/EcR-B ecdysone receptor for mushroom body neuronal remodeling in Drosophila. Neuron 28:807–818.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Buszczak M,
    2. Segraves WA
    (2000) Insect metamorphosis: Out with the old, in with the new. Curr Biol 10:830–833.
    OpenUrlCrossRef
  18. ↵
    1. Lanzi G,
    2. et al.
    (2006) Molecular and biological characterization of deformed wing virus of honeybees (Apis mellifera L.) J Virol 80:4998–5009.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Maori E,
    2. et al.
    (2007) Isolation and characterization of Israeli acute paralysis virus, a dicistrovirus affecting honeybees in Israel: Evidence for diversity due to intra- and inter-species recombination. J Gen Virol 88:3428–3438.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Iatrou K,
    2. Gilbert L,
    3. Gill S
    1. Feyereisen R
    (2005) in Comprehensive Molecular Insect Science, Insect cytochrome P450, eds Iatrou K, Gilbert L, Gill S (Elsevier Pergamon), Vol 4, pp 1–77.
    OpenUrlCrossRef
  21. ↵
    1. Maïbèche-Coisne M,
    2. Monti-Dedieu L,
    3. Aragon S,
    4. Dauphin-Villemant C
    (2000) A new cytochrome P450 from Drosophila melanogaster, CYP4G15, expressed in the nervous system. Biochem Biophys Res Comm 273:1132–1137.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Claudianos C,
    2. et al.
    (2006) A deficit of detoxification enzymes: Pesticide sensitivity and environmental response in the honeybee. Insect Mol Biol 15:615–636.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Li W,
    2. Ma G,
    3. Zhou X
    (2006) Apidaecin-type peptides: Biodiversity, structure–function relationships and mode of action. Peptides 27:2350–2359.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Sabatier L,
    2. et al.
    (2003) Pherokine-2 and −3. Eur J Biochem 270:3398–3407.
    OpenUrlPubMed
  25. ↵
    1. Roxström-Lindquist K,
    2. Terenius O,
    3. Faye I
    (2004) Parasite-specific immune response in adult Drosophila melanogaster: A genomic study. EMBO Rep 5:207–212.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Gillespie JJ,
    2. Johnston JS,
    3. Cannone JJ,
    4. Gutell RR
    (2006) Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): Structure, organization, and retrotransposable elements. Insect Mol Biol 15:657–686.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Gonzalez IL,
    2. Sylvester JE
    (1997) Incognito rRNA and rDNA in databases and libraries. Genome Res 7:65–70.
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Scharf ME,
    2. Wu-Scharf D,
    3. Zhou X,
    4. Pittendrigh BR,
    5. Bennett GW
    (2005) Gene expression profiles among immature and adult reproductive castes of the termite Reticulitermes flavipes. Insect Mol Biol 14:31–44.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Slomovic S,
    2. Portnoy V,
    3. Liveanu V,
    4. Schuster G
    (2006) RNA polyadenylation in prokaryotes and organelles: Different tails tell different tales. Crit Rev Plant Sci 25:65–77.
    OpenUrlCrossRef
  30. ↵
    1. Slomovic S,
    2. Laufer D,
    3. Geiger D,
    4. Schuster G
    (2006) Polyadenylation of ribosomal RNA in human cells. Nucleic Acids Res 34:2966–2975.
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Win TZ,
    2. et al.
    (2006) Requirement of fission yeast Cid14 in polyadenylation of rRNAs. Mol Cell Biol 26:1710–1721.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Mauro VP,
    2. Edelman GM
    (2002) The ribosome filter hypothesis. Proc Natl Acad Sci USA 99:12031–12036.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Tranque P,
    2. Hu MC,
    3. Edelman GM,
    4. Mauro VP
    (1998) rRNA complementarity within mRNAs: A possible basis for mRNA–ribosome interactions and translational control. Proc Natl Acad Sci USA 95:12238–12243.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Chang Y,
    2. Stockinger MP,
    3. Tashiro H,
    4. Lin CG
    (2008) A novel noncoding RNA rescues mutant SOD1-mediated cell death. FASEB J 22:691–702.
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Kong Q,
    2. Stockinger MP,
    3. Chang Y,
    4. Tashiro H,
    5. Lin CG
    (2008) The presence of rRNA sequences in polyadenylated RNA and its potential functions. Biotech J 3:1041–1046.
    OpenUrlCrossRef
  36. ↵
    1. Ongus JR,
    2. Roode EC,
    3. Pleij CWA,
    4. Vlak JM,
    5. van Oers MM
    (2006) The 5′ non-translated region of Varroa destructor virus 1 (genus Iflavirus): Structure prediction and IRES activity in Lymantria dispar cells. J Gen Virol 87:3397–3407.
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Joachims M,
    2. Van Breugel PC,
    3. Lloyd RE
    (1999) Cleavage of poly(A)-binding protein by enterovirus proteases concurrent with inhibition of translation in vitro. J Virol 73:718–727.
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Hellen CT,
    2. Sarnow P
    (2001) Internal ribosome entry sites in eukaryotic mRNA molecules. Genes Dev 15:1593–1612.
    OpenUrlFREE Full Text
  39. ↵
    1. Zundel MA,
    2. Basturea GN,
    3. Deutscher MP
    (2009) Initiation of ribosome degradation during starvation in Escherichia coli. RNA 15:977–983.
    OpenUrlAbstract/FREE Full Text
View Abstract
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Changes in transcript abundance relating to colony collapse disorder in honey bees (Apis mellifera)
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
Citation Tools
Changes in transcript abundance relating to colony collapse disorder in honey bees (Apis mellifera)
Reed M. Johnson, Jay D. Evans, Gene E. Robinson, May R. Berenbaum
Proceedings of the National Academy of Sciences Sep 2009, 106 (35) 14790-14795; DOI: 10.1073/pnas.0906970106

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Changes in transcript abundance relating to colony collapse disorder in honey bees (Apis mellifera)
Reed M. Johnson, Jay D. Evans, Gene E. Robinson, May R. Berenbaum
Proceedings of the National Academy of Sciences Sep 2009, 106 (35) 14790-14795; DOI: 10.1073/pnas.0906970106
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 116 (8)
Current Issue

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

News Feature: Cities serve as testbeds for evolutionary change
Urban living can pressure flora and fauna to adapt in intriguing ways. Biologists are starting to take advantage of this convenient laboratory of evolution.
Image credit: Kristin Winchell (Washington University in St. Louis, St. Louis).
Several aspects of the proposal, which aims to expand open access, require serious discussion and, in some cases, a rethink.
Opinion: “Plan S” falls short for society publishers—and for the researchers they serve
Several aspects of the proposal, which aims to expand open access, require serious discussion and, in some cases, a rethink.
Image credit: Dave Cutler (artist).
Featured Profile
PNAS Profile of NAS member and biochemist Hao Wu
 Nonmonogamous strawberry poison frog (Oophaga pumilio).  Image courtesy of Yusan Yang (University of Pittsburgh, Pittsburgh).
Putative signature of monogamy
A study suggests a putative gene-expression hallmark common to monogamous male vertebrates of some species, namely cichlid fishes, dendrobatid frogs, passeroid songbirds, common voles, and deer mice, and identifies 24 candidate genes potentially associated with monogamy.
Image courtesy of Yusan Yang (University of Pittsburgh, Pittsburgh).
Active lifestyles. Image courtesy of Pixabay/MabelAmber.
Meaningful life tied to healthy aging
Physical and social well-being in old age are linked to self-assessments of life worth, and a spectrum of behavioral, economic, health, and social variables may influence whether aging individuals believe they are leading meaningful lives.
Image courtesy of Pixabay/MabelAmber.

More Articles of This Classification

Biological Sciences

  • Unnatural verticilide enantiomer inhibits type 2 ryanodine receptor-mediated calcium leak and is antiarrhythmic
  • Unleashing floret fertility in wheat through the mutation of a homeobox gene
  • Arg302 governs the pKa of Glu325 in LacY
Show more

Agricultural Sciences

  • Crop variety management for climate adaptation supported by citizen science
  • A salivary effector enables whitefly to feed on host plants by eliciting salicylic acid-signaling pathway
  • Small-scale urban agriculture results in high yields but requires judicious management of inputs to achieve sustainability
Show more

Related Content

  • In This Issue
  • Scopus
  • PubMed
  • Google Scholar

Cited by...

  • Honey Bee Deformed Wing Virus Structures Reveal that Conformational Changes Accompany Genome Release
  • Early gut colonizers shape parasite susceptibility and microbiota composition in honey bee workers
  • Evidence of Apis ceranaSacbrood virus Infection in Apis mellifera
  • Systemic Spread and Propagation of a Plant-Pathogenic Virus in European Honeybees, Apis mellifera
  • Microsporidian genome analysis reveals evolutionary strategies for obligate intracellular growth
  • Global Honey Bee Viral Landscape Altered by a Parasitic Mite
  • CYP9Q-mediated detoxification of acaricides in the honey bee (Apis mellifera)
  • Clarity on Honey Bee Collapse?
  • Scopus (147)
  • Google Scholar

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Latest Articles
  • Archive

PNAS Portals

  • Classics
  • Front Matter
  • Teaching Resources
  • Anthropology
  • Chemistry
  • Physics
  • Sustainability Science

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Press
  • Site Map

Feedback    Privacy/Legal

Copyright © 2019 National Academy of Sciences. Online ISSN 1091-6490