Parallel Epigenomic and Transcriptomic Responses to Viral Infection in Honey Bees ()
Honey bees are a critical pollinator of a wide variety of agricultural crops, but beekeepers experience heavy annual losses of honey bee colonies. Several factors are associated with colony losses, in particular infection with Israeli Acute Paralysis Virus (IAPV). Despite the importance of viruses to honey bee health, our understanding of the molecular mechanisms mediating host-pathogen interactions is limited. Here, we characterized the epigenomic and transcriptomic responses of honey bees to short term (<24 hour) IAPV infection. We found significant expression differences in 753 genes between infected and control bees, including genes involved in immune and epigenetic pathways. However, IAPV-regulated genes did not overlap significantly with genes that respond to bacterial or microsporidian infection, suggesting that honey bees utilize distinct pathways for different immune challenges. Additionally, DNA methylation status of 156 genes changed significantly, including several genes that have been linked to antiviral immune responses in humans. Interestingly, there was no significant overlap between the differentially methylated and differentially expressed genes, suggesting that honey bees may possess a parallel genomic mechanisms to respond to viral infection.
Published in the journal:
Parallel Epigenomic and Transcriptomic Responses to Viral Infection in Honey Bees (). PLoS Pathog 11(3): e32767. doi:10.1371/journal.ppat.1004713
Category:
Research Article
doi:
https://doi.org/10.1371/journal.ppat.1004713
Summary
Honey bees are a critical pollinator of a wide variety of agricultural crops, but beekeepers experience heavy annual losses of honey bee colonies. Several factors are associated with colony losses, in particular infection with Israeli Acute Paralysis Virus (IAPV). Despite the importance of viruses to honey bee health, our understanding of the molecular mechanisms mediating host-pathogen interactions is limited. Here, we characterized the epigenomic and transcriptomic responses of honey bees to short term (<24 hour) IAPV infection. We found significant expression differences in 753 genes between infected and control bees, including genes involved in immune and epigenetic pathways. However, IAPV-regulated genes did not overlap significantly with genes that respond to bacterial or microsporidian infection, suggesting that honey bees utilize distinct pathways for different immune challenges. Additionally, DNA methylation status of 156 genes changed significantly, including several genes that have been linked to antiviral immune responses in humans. Interestingly, there was no significant overlap between the differentially methylated and differentially expressed genes, suggesting that honey bees may possess a parallel genomic mechanisms to respond to viral infection.
Introduction
Honey bee populations are in decline throughout the world [1]. A recent survey found that US beekeepers lose 30% of their colonies annually [2]. Because honey bees are critical pollinators of our agricultural crops and over 70% of major global food crops, including fruits, vegetables and nuts, benefit from honey bees and other pollinators [3], these losses have substantial implications for the sustainability of our agricultural industry. Several factors are thought to contribute to these declines, including pathogens, parasites, habitat loss, poor nutrition due to monocropping systems, and pesticide use [1]. Furthermore, several of these factors appear to interact synergistically; for example, the titers and impacts of viruses increase in bees infested with Varroa mites [4], fed poor diets [5] or exposed to pesticides [6]. Honey bees have several viruses, with more than 20 identified thus far [7,8]. These viruses can have a wide range of impacts, from no obvious symptoms to developmental defects [9], altered learning and memory [10], and loss of muscle coordination and premature death [11]. In several studies, increasing numbers of viruses or viral titers have been linked to colony losses [12–14]. There is currently a wealth of information on the pleiotropic effects of viral infection in honey bees, but information about the molecular and physiological responses of honey bees to viral infections, especially on the genome-wide scale, is surprisingly limited.
One of the viruses that has been implicated in the decline of honey bees is Israeli Acute Paralysis Virus (IAPV). IAPV is a positive sense RNA virus in the family Dicistroviridae [15]. Infection of honey bee pupae results in cessation of development and altered expression of several candidate genes involved in ribosomal biogenesis, though the impacts of infection can be quite variable [16]. Infection of caged adult bees results in decreased locomotion, muscle spasms, and increased mortality, with 80% of infected bees dying prematurely[11]. Though IAPV infection was originally associated with the phenomenon known as "Colony Collapse Disorder" [17,18], subsequent studies failed to provide a strong association between IAPV infection and colony losses associated with this syndrome [15,19,20]. However, a recent study demonstrated that IAPV is more prevalent in weak colonies which fail to survive the winter [21]. Thus, though the role of IAPV in colony losses is not well understood, it seems likely that presence of this virus reduces colony strength and survival.
The molecular mechanisms involved in the antiviral responses in insects have been best studied in Drosophila melanogaster. Antiviral response appear to be mediated by two primary pathways, one involving JAK-STAT (Janus kinase—signal transducers and activators of transcription) and the other RNA interference (RNAi) [22]. The JAK-STAT pathway uses cytokine molecules to activate the signaling cascade upon binding to receptors on the cell surface. The activated receptor then causes a dimerization of Stat transcription factors, which enter the nucleus and binds to downstream viral response genes with Stat binding sites in their promoter regions. The RNAi pathway functions by cleaving double stranded RNA (dsRNA) into small fragments, which are used to target other mRNA transcripts with the same sequence, and thereby preventing the translation from mRNA to protein and effectively decreasing the activity of the particular gene. Mutant D. melanogaster that do not possess a functional RNAi system have increased mortality and higher viral titers, suggesting that this pathway is an essential part of the antiviral immunity in flies [23]. Previous studies have suggested that the RNAi pathway may also be involved in the honey bee antiviral response, since honey bees show decreased viral titers after treatment with virus specific dsRNA and non-specific dsRNA[11,24–26]. However, though these studies demonstrate the principle that viral titers can be reduced through RNAi, they did not demonstrate that transcriptional regulation of the RNAi pathway is an antiviral response in honey bees. Other recent studies have examined transcriptional responses to Deformed Wing Virus in pre-adult stages[27], three days post infection with a model virus (Sindbis) in young worker bees[24], and in older bees chronically infected with IAPV[21], but the global transcriptional and epigenetic responses to acute viral infections in adult bees has not been investigated.
Epigenetic mechanisms, such as RNAi and DNA methylation, have been shown to play a role in mediating antiviral defenses in vertebrates and plants through post translational modification (protein folding, RNAi, and alternative splicing) [28–30]. Epigenetic mechanisms enable modification of function or expression of DNA, RNA, or proteins without changing the original DNA sequence. In humans and plants, viral infection (of both integrating and non-integrating viruses) has the ability to alter the host DNA methylation patterns[31–36]. Though the majority of insect immune studies have involved D. melanogaster, it does not possess a full complement of DNA methylation enzymes, and thus while DNA methylation is present, it is at very low levels[37,38]. However, many other insect species, including honey bees, do possess a fully functional DNA methylation system [39], and thus DNA methylation may represent an important antiviral response in insects as well.
In insects, most methylation occurs in the gene body (from the transcriptional start site to the end of transcription, encompassing both exons and introns only), and how this impacts gene function is still under investigation[40,41]. DNA methylation involves the addition of a methyl group (CH3) to the 5th carbon in the cytosine pyrimidine ring and primarily occurs within CpG (cytosine phosphate guanine) dinucleotides, but it also occurs to a much lesser extent in CHG and CHH sequence contexts (H = A, C, or T) [42]. Methylated DNA can recruit histone modifying proteins that then alter the structure of the chromatin [42]. This change in DNA packaging can have different effects depending on the genomic context in which it occurs. If DNA methylation occurs within or near a promoter region, the outcome is typically transcriptional silencing [42]. The function of DNA methylation within the gene body is less clear; it has been proposed to function in alternative splicing [43], phenotypic plasticity [44], and suppression of transcriptional noise [45,46]. More research is needed to definitively characterize the function of intragenic DNA methylation.
Here, we examine the transcriptional and epigenetic (DNA methylation) responses of young worker honey bees to acute (20–24 h) infection with IAPV. In control and virus infected workers, we (1) characterized genome-wide gene expression patterns in the fat body tissues, the major sites of metabolism, detoxification and immune responses in insects (2) characterized associated genome-wide DNA methylation patterns, and (3) determined if there is a correlation between gene expression and DNA methylation patterns. Our studies provide insight into the transcriptional and epigenetic mechanisms underlying both antiviral responses in honey bees and the pleiotropic effects of these viruses, while also allowing us to examine in detail the relationship between gene expression and DNA methylation.
Methods
Identification of IAPV-free experimental colonies
Colonies of honey bees were maintained at apiaries at Penn State using standard commercial apiculture practices. Colonies headed by single-drone inseminated (SDI) queens (obtained from Glenn Apiaries, Fallbrook, CA) were screened to identify colonies that were free of detectable IAPV infections and had, overall, the lowest numbers and titers of viruses. RNA was isolated from the pooled abdomens of five individuals from each SDI colony using an RNeasy (Qiagen, Valencia, CA) protocol. Using a standard protocol[47], 150 ng of RNA was converted to cDNA and quantitative real time PCR (qRT-PCR) was performed using primers specific for the following viruses: Black Queen Cell Virus (BQCV), Deformed Wing Virus (DWV), Israeli Acute Paralysis Virus (IAPV), and Kashmir Bee Virus (KBV) (see S1B Table for a listing of the primers used and the associated references). Workers from the colony found to be IAPV-free and with the lowest total viral loads were selected and used in subsequent analyses. Note that the lack of IAPV viral infection and reduced levels of other viruses in control bees was confirmed by qRT-PCR and by scanning for viral transcripts in the transcriptomic data (S1C&D Table).
Treatments
Honeycomb frames of emerging brood were collected placed in a 34°C, 50%RH incubator overnight. Emerged bees (<24 hours old) were collected. Individual workers were fed extracts of virus-infected bees (which contained IAPV, DWV, BQCV, KBV and SBV; bees were obtained from D. Cox-Foster, Penn State) in a 50%-sucrose solution while control bees were fed a 50%-sucrose control (4 μL/bee). These honey bee viruses currently cannot be produced as pure strains under laboratory conditions, so we modified the protocol from Li et al. [48] to extract the viruses from live bees known to be infected with multiple viruses. The viral extracts were obtained by homogenizing two whole bees in 2ml of molecular grade water. Prior to feeding, the inoculum was dialyzed through a 0.20 μm cell filter to purify the virus extract, as in Hunter et al. and Liu et al [25,49]. Sucrose was added to produce a 50% sucrose solution.
Bees were maintained in Plexiglas cages (10 x 10 x 7 cm) in groups of 45 under red light in a 34.5°C and ∼50% RH incubator. Bees were fed 50% sucrose/water and a 50% MegaBee diet/honey mixture ad libitum. Each cage was also provided with 0.1 queen equivalents of queen mandibular pheromone (QMP) to approximate normal colony conditions [50]. Bees were monitored hourly for symptoms of viral infection (trembling, twitching) starting 15 hours post infection. Samples for analysis were taken from control and virus-infected bees that began exhibiting symptoms 20–24 hours after treatment (when they were <48 hours old). Only bees exhibiting the neurological symptoms of a successful viral infection were sampled from the virus-infected treatments.
Confirmation of viral infection
A subset of control and infected individuals was used to confirm viral infection (or lack thereof, in the case of the controls) using qRT-PCR, as above. RNA was isolated from eviscerated abdomens with attached fat bodies of individual bees using TRIzol Reagent (Invitrogen, Carlsbad, CA). Bees were individually screened for infections with Chronic Bee Paralysis Virus (CBPV), Acute Bee Paralysis Virus (ABPV), Sacbrood Virus (SBV) and IAPV (see S1B Table for primers and associated references). Furthermore, we confirmed the presence of the negative strand of IAPV in the treated bees, demonstrating that the virus was actively replicating (S1 Fig.).
Sample preparation for RNA-Seq and BS-Seq
Eviscerated abdomens with attached fat bodies from individual bees were dissected in ice-cold RNAlater (Qiagen, Valencia, CA) and excess RNAlater was removed from the samples. RNA was isolated using TRIzol Reagent (Invitrogen, Carlsbad, CA). Samples were then treated with TURBO DNase (Ambion, Austin, TX) to remove DNA contamination. The samples were then pooled into three groups containing three individuals for both treatment and control groups. The quality and concentrations of the RNA pools was verified using a 2100 Bioanalyzer (Agilent, Santa Clara, CA).
DNA was isolated from the fat bodies from a second set of bees using a Gentra Puregene kit (Qiagen, Valencia, CA). DNA quality was assessed by electrophoresis, using 1% Agarose gel, and the concentration was determined using a Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA). Due to the large amount of DNA required for BS-Seq, the DNA was pooled to produce a single control (n = 9 bees) and a single treatment sample (n = 9 bees).
RNA-Seq and BS-Seq
The RNA and DNA samples were sent to Beijing Genomics Institute (BGI)(Shenzhen, China) for library preparation and sequencing. The library preparation (each library was individually tagged for each pool) and Illumina paired end sequencing of the RNA samples and bisulfite conversion and Illumina paired end sequencing of DNA samples was performed by BGI. RNA and DNA Samples were sequenced on two separate lanes of the Illumina Hiseq2000 platform, resulting in paired end reads that were 2x100bp in length for the RNA samples and 2x90bp for the DNA samples. The total number of reads for the RNA-seq samples were between 51,231,170 and 53,823,222, and for the BS-seq, the total number of reads for each sample were between 101,423,954 and 102,213,608. Data generated from these studies has been deposited in NCBI’s Gene Expression Omnibus[51] and are accessible through the GEO Series accession number GSE65659.
Analysis of RNA-Seq data
Transcriptome sequencing reads were preprocessed using Illumina HCS 1.1 software. We removed adaptor sequences, reads composed of more than 5% unknown nucleotides, and reads with more than 20% of the base qualities under 10. The transcriptome sequencing reads were aligned to the most recent honey bee genome assembly (Amel_4.5) [52] using Tophat[53]. Aligned read counts were imported into R statistical software (http://www.r-project.org). Genes with low read counts (less than 5 reads per gene) were removed. The data was normalized using a trimmed mean of M-values (TMM) method [54]. The DESeq package in R statistical program was used to identify significantly differentially expressed genes [55].
For the gene ontology (GO) analysis, D. melanogaster orthologs of the differentially expressed genes were uploaded to DAVID Bioinformatics Resources 6.7[56]. The overrepresented GO terms of differentially expressed genes during viral infection were determined using a background list of all genes expressed among the samples. To refine the list of GO terms, the significant terms (FDR < 0.05) that were generated using DAVID Bioinformatics Resources were imported into REVIGO [57], which clusters the GO terms using semantic similarity measures.
qRT-PCR validation of expression differences
To validate the results of the transcriptomic analysis, we repeated the experiment using a new set of bees. Emerging worker bees were collected and fed either viral extracts in a 50% sucrose solution or a 50% sucrose control solution as before. As an additional control, we also included a bee lysate control, using healthy, non-symptomatic bees, that was extensively washed to remove as much of the virus as possible. The groups were reared under the same conditions as before, and monitored hourly for symptoms of viral infection starting 15 hours post infection. All bees were collected 24 hours after treatment. RNA was extracted as above and qRT-PCR was then used to validate the RNA-seq results (see S1B Table for a listing of the primers and associated references).
Analysis of BS-Seq data
The bisulfite-converted DNA reads in both control and treated group are aligned to honey bee genome Apimel4.5 using BSMAP [58] Ambiguously mapped reads were filtered with parameter-u and potential PCR duplications were removed by parameter-r. To avoid erroneous annotation, we restricted our analyses to genes that fit the following criteria: 1) mapped by RNA-seq data, 2) have orthologs in D. melanogaster, and 3) with corresponding assembly in the previous assembly 2.0. Following this filtering step, we examined a total of 12,446 genes.
Methylated CpG sites are determined by the binomial test [59]. The error rate (non-conversion rate and sequencing error) is estimated from the sequencing result of the unmethylated λ phage mixed in our sample [60]. Error rates estimated from these were 0.0041 and 0.0037 from control and treatment groups, respectively. P-values were subsequently corrected for multiple testing [61], and sites with adjusted p-value < 0.1 were considered methylated.
Fractional methylation values were calculated for each CpG site as mCG/CG, where mCG is the number of reads with a methylated cytosine at a CpG site (according to non-conversion) and CG is the total number of reads mapped to the site [41]. Fractional methylation for specific genomic regions (for example, gene body, exon or intron) was calculated as the mean of all CpG fractional methylation values within that region. The degree of CpG depletion, or CpGO/E, is calculated as in Elango et al. [62].
Differential methylation
The differential methylation level of each CpG site was assessed by a Fisher’s exact test with P-values adjusted for multiple testing [61]. Differential methylation of specific genomic regions (such as gene bodies) was assessed by a generalized linear model with binomial family, similar to [59,63]. In the generalized model, the response vector is the number of methylated and nonmethylated CpGs in specific genomic regions, modeled by two categorical factors, treatment (control or infected) and CpGs locations [59,63].
Alternative splicing
Transcriptome sequencing reads were mapped to honeybee assembly 4.5 and processed by TopHat[53]. Upon detecting a splice junction, two types of alternative splicing events were summarized by comparing them to the known annotation (S2 Fig.): I) Intron retention if a pair of junctions occur inside an annotated exon (intron is spliced in our sample, and retained in the reference sequence) and II) Exon skipping if the region between a pair of junction covers one or more annotated exons (exon is skipped in our sample).
Results and Discussion
Confirmation of viral infection levels
We used qRT-PCR to examine titers of viruses in our control and treated samples. For DWV, CBPV, and BQCV, in all samples (including both treatment and control) viral RNA was either undetectable or showed very late amplification (average Ct > 35), while IAPV titers were detected in the treatment samples only (see S1C Table). We also confirmed the presence of the negative strand of IAPV, demonstrating viral replication in the treatment samples (see S1 Fig.). As another measure to confirm successful viral infection, the reads from the RNA-seq data set were aligned to a panel of viral genomes (see S1D Table for accession numbers). The number of reads mapping to IAPV in the treatment samples were substantially higher than those in the controls (S1D Table). Reads for additional viruses (BQCV, DWV, and KBV) were also present in the samples, but at much lower abundance (over 2000x lower than the IAPV reads in the treated group) and were equally abundant in both control and treated samples. While we screened the transcriptome for the common viruses that are pathogenic to honey bees in North America (S1D Table), it is possible that additional viruses were present that we did not account for.
Transcriptional response to viral infection
Using RNA-seq to monitor genome-wide expression patterns in the fat bodies of control and treated bees, we identified 753 genes that were significantly differentially expressed (FDR p-value < 0.05; S1A Table), with 607 genes up-regulated and 146 down-regulated in response to viral infection (Fig. 1A). IAPV-infected and control bees clearly had distinct expression profiles and were clearly separated according to a hierarchical clustering analysis (Fig. 1B). Gene ontology analysis revealed that the genes that were down-regulated were related to response to toxic substances, cofactor metabolic process, glucose catabolic process, and alcohol catabolic process (S1E Table), while genes that were up-regulated were significantly enriched for categories including regulation of MAPK cascade, metabolism, morphogenesis, phosphorylation, autophagic cell death, nucleic acid metabolism, and negative regulation of transcription from RNA polymerase II promoter (S1F Table). Interestingly, transcriptional pausing of RNA polymerase II has been suggested as an innate immune response in D. melanogaster, priming the host genome by increasing the accessibility of promoter regions associated with virally induced genes, allowing for a more rapid response[64]. Our data also indicate that several core insect immune pathways that have been previously associated with antiviral defense in other insects [22,65,66], including genes from two immune signaling pathways (Toll and Jak-STAT) and the RNAi pathway, are involved in the acute antiviral immune response in honey bees. We found a significant upregulation of multiple representatives of the RNAi pathway (Fig. 2A), including Argonaute-2 and Dicer-like using qRT-PCR. Honey bees encode multiple Dicer proteins (Dicer-1 and Dicer-like), but Dicer-like is primarily involved in RNAi[67]. Using qRT-PCR, we confirmed upregulation of genes involved in the RNAi and Toll pathways (Fig. 2B). All but one of the candidate genes (tarbp2-like) we tested matched the expression pattern that was seen in the RNA-seq data. This is the first time the RNAi pathway has been implicated in a transcriptional antiviral response in honey bees, though previous studies have demonstrated that feeding with double-stranded RNA constructs, both virus specific and non-specific, (which induces an RNAi response) can reduce viral titers[11,24–26,49].
Next, we performed a series of comparative studies with previous-genome wide analyses of transcriptional responses to parasite and pathogen infections in honey bees and D. melanogaster. A comparison of our suite of regulated genes with the 177 genes identified as part of the honey bee immune response pathways in Evans et al. [68] revealed an overlap of only 19 genes (S1G Table), demonstrating that there are many other genes that may play a role in antiviral responses (Fig. 3A). Comparing our gene expression results to an antiviral study in D. melanogaster by Xu et al. [64], we found an overlap of 29 genes (S1H Table). Interestingly, Xu et al. also found that the transcriptional pausing pathway was involved in the D. melanogaster acute antiviral response, including upregulation of the gene encoding negative elongation factor (NELF), which plays an important role in transcriptional pausing and was also found to be upregulated in our study. Indeed, of the 29 genes overlapping between our study and Xu et al, there were two nonsignificantly overrepresented GO terms, regulation of transcription and autophagic cell death, further suggesting a link between transcriptional pausing and antiviral response.
Previous studies identified genes differentially regulated in the fat bodies of worker bees injected with saline, beads or Escherichia coli bacteria [69] or of worker bees chronically infected with microsporidian (Nosema spp.) intestinal parasites [70]. We found little overlap among the genes identified in these studies and those responding to viral infection in our current study (38 for N. spp. and 27 for E. coli), and only 5 were present in all three treatments (Fig. 3B and S1I Table). The lack of a significant overlap (p > 0.05, Fisher Exact Test) with either of these two studies suggests that honey bees use different molecular mechanisms to defend themselves against different pathogens. When we compared our gene expression results to three other studies (Fig. 3C, S1I Table) examining transcriptional responses to viral infection in honey bees (DWV, IAPV, and Sindbis Virus (a model virus)), we found a greater degree of overlap, however, there is still variability among the antiviral responses. This variability could be due to differences in antiviral response across life stages or to the different viruses. Notably, of all the studies we compared, the greatest overlap with our results was found with a study examining the effects of chronic IAPV infection on older adult workers[21].
Characteristics of methylated genes in honey bees
In agreement with previous studies[62], we found that DNA methylation is mostly targeted to gene bodies (exons and introns, Fig. 4). Fractional methylation levels of gene bodies exhibit a distinctive bimodal pattern where one group of genes shows little DNA methylation (fractional methylation level peak at 0.03) while the other group of genes shows substantial methylation peak at 0.3 (Fig. 4). This pattern is similar to the results of a previous analysis of experimental data [71] as well as computational predictions [62,72,73]. The bimodal pattern of DNA methylation of genic regions is consistently observed when we restrict our analyses to CpGs occurring only in exons or introns, or in the first 500 bps of a gene (Fig. 4). Based upon this observation, we tentatively refer to genes with fractional methylation levels of exons < 0.03 as ‘non-methylated genes (nonMGs)’, and those greater than 0.03 as ‘methylated genes (MGs)’. Defined this way, most genes remained in the same group (MGs or nonMGs) in the control and treatment groups (see below).
Within a methylated gene, exons are much more heavily methylated than introns, similar to the pattern in Nasonia vitripennis [74]. DNA methylation levels peak in the second exon (Fig. 5A) and sharply decrease outside of the coding regions of methylated genes (Fig. 5B). Non-methylated genes are significantly longer than methylated genes (P < 10–16; Mann-Whitney test), as previously inferred based upon nucleotide composition [75]. Moreover, within methylated genes, gene length is negatively correlated with fractional methylated levels of each exon (Spearman’s correlation coefficient-0.34, P < 10–16). To exclude the possibility that longer genes appear to be less methylated since DNA methylation is mostly targeted to the 5’ end of genes, we examined the relationship between gene length and DNA methylation using the fractional methylation levels of only the first 500 bps of genes. Gene length is still significantly negatively correlated with methylation level in this analysis (Spearman’s correlation coefficient-0.10, P < 10–15), indicating that the relationship between DNA methylation and gene (and exon) length is likely to have true biological basis.
DNA methylation and gene expression
We next examined the relationship between DNA methylation and gene expression. Overall, we found that methylated genes are expressed at significantly higher levels than non-methylated genes (P<10–16, Mann-Whitney test, for both control and treatment, Fig. 6A&B), consistent with the previous studies demonstrating that gene body methylation promotes expression [62,71,76,77]. However, as shown above, DNA methylation is tightly correlated with gene length (longer genes are less methylated). It is also known that gene expression levels co-vary with gene length, in that long genes tend to exhibit reduced expression [78,79]. Fig. 6 demonstrates the interactions between gene length, DNA methylation and gene expression levels, by illustrating the relationship between gene length and gene expression in different methylation bins. First, genes with extremely low levels of DNA methylation are expressed at significantly lower levels than others. For example, genes with methylation < 0.02 have significantly lower expression compared to the rest of the genes for both control and treatment (P < 10–16; Mann-Whitney test). Second, gene expression level is strongly associated with gene length. Longer genes generally were expressed at lower levels than shorter genes, except the case of the lowest methylated quantile. Thus, both gene length and DNA methylation affect gene expression, which is also supported by an ANOVA (Table 1). Both control and infected groups exhibited similar trends.
Differential DNA methylation as a response to viral infection
We examined the changes of DNA methylation as a response to viral infection by using a generalized linear model where CpG positions are explicitly modeled (Materials and Methods). In a moderately stringent condition of FDR < 0.01, only 40 genes are identified as differentially methylated (S2A Table). Even when we relax the statistical criterion to FDR < 0.1 (a cutoff selected based upon the empirical distributions of P-values to better reflect the data), the number of differentially methylated genes (DMGs) only increases to 156 (S2A Table). Among these 156 putatively differentially methylated genes, 33 of them are more highly methylated in controls while 123 of them are more highly methylated in treatment. Several genes involved in cellular signaling were found to be differentially methylated including five genes that are associated with the T-cell activation pathway in vertebrates, a major immune pathway (S2B&C Table). These genes are inositol 1,4,5,-tris-phosphate receptor (IP3R), phosphotidylinositol 3 kinase (PI3K), protein kinase C-delta (PKC-delta), Vav-like (Vav), and son of sevenless (Sos). All of these genes interact with RNA viruses in human systems. For example, IP3R is required for the release of human immunodeficiency virus (HIV) particles[80]. Inhibition of PI3K negatively affects cell entry of influenza virus and hepatitis C virus[81,82] thereby affecting viral replication. Inhibition of PKC-delta has a negative effect on HIV post cell entry replication[83]. Dominant negative Vav constructs decrease intercellular trafficking of a critical HIV protein[84]. Finally, inactivation of Sos results in an inhibition of protein kinase (PKR) activity, ultimately leading to a resistance to reovirus infection[85]. Unfortunately, information on how these genes respond to viral infection in other species is severely understudied. Though our results and other recent studies have indicated that DNA methylation can be involved in rapid processes—indeed, even learning[86]—it is possible that there are significantly more differentially methylated genes at later time points.
Differential DNA methylation does not lead to differential gene expression
We then directly examined whether significantly differentially methylated genes (DMGs) are also differentially expressed genes (DEGs), as expected if differential DNA methylation drives the transcriptomic responses to viral infection. In contrast to this prediction, we found that there is little overlap between DMGs and DEGs—only 4 genes were found to be both differentially expressed and differentially methylated, which is not a significant association according to the Fisher’s exact test (P = 0.068; 5 genes would be expected to be shared by chance). Among the 4 genes that were both differentially expressed and differentially methylated, three were up-regulated in response to increased DNA methylation, while the other was down-regulated in response to increased DNA methylation. The lack of correspondence between DMGs and DEGs may have been a result of using different sets of bees for each analysis, which was necessary due to the high quantities of RNA and DNA required for these analyses. However, the lack of overlap between DMGs and DEGs was further supported by single CpG level analysis where neither enrichment nor reduction of differentially methylated CpG sites was observed in the 753 differential expressed genes (P = 0.38; Chi-square test).
In fact, DEGs and DMGs had distinctive characteristics at many levels (Fig. 7). DMGs tend to be highly methylated, and consequently, highly expressed genes (P < 10–10; Mann-Whitney test). On the other hand, DEGs are significantly less methylated than non-DEGs (mean fractional methylation levels of DEGs is 0.063 compared to 0.197 of non-DEGs, P < 10–16; Mann-Whitney test). Among the 753 DEGs, 610 of them were classified as non-methylated in both controls and virus-infected bees. DEGs also were highly expressed (P < 10–16; t test), however, they were not significantly longer or shorter than other genes (median coding region length of DEGs was 1261 as compared to 1203 of non-DEGs, P = 0.1, t-test). But coding regions of DEGs are significantly shorter than DMGs (P < 10–16; Mann-Whitney test). Thus our study suggests that different types of genes respond to biological stimuli in different ways: some genes tend to be differentially expressed, while others tend to be differentially methylated. These two types of genes exhibit different, almost completely opposite, genomic characteristics. However, transcriptional regulation seems to occur at a much greater level in response to viral infection than methylation, with 753 DEGs (FDR<0.05) and only 156 DMGs (FDR <0.1).
DNA methylation and alternative splicing
Recent studies have demonstrated significant roles of DNA methylation on alternative splicing [43,87–89]. While alternative splicing indicates usage of alternative transcripts, it can also interfere with normal gene function. For example, retaining a non-coding portion of a gene (by intron retention) or missing a crucial exon (exon skipping) in the final transcript may lead to misfolded or non-functional proteins. Thus, DNA methylation changes may reduce gene function via splicing, rather than transcription.
Here we examined the relationship between DNA methylation and alternative splicing (Materials and Methods). We identified 1345 genes with exon skipping (with 90 genes unique to controls, and 43 unique to treatment; S2D&E Table) and 1578 with intron retention (with 81 genes unique to controls, and 156 unique to treatment; S2F&G Table). GO analysis found that genes with exon skipping are enriched in glycerophospholipid metabolic process. Interestingly, glycerophospholipid metabolism plays a key role in flock house virus (FHV, which, like IAPV, is a positive sense RNA virus) replication in D. melanogaster, where it is involved in FHV RNA replication complex assembly and function[90].
We found that exons involved in intron retention are significantly more methylated than other exons in the same genes, in both biological conditions (P < 10–3, paired Mann-Whitney test). Exons that are skipped in an exon-skipping event show also higher methylation level compared to other exons of that gene (P <10–9, paired Mann-Whitney test, Fig. 8). However, no significant overlap was observed between those differentially spliced genes and DMGs. Moreover, exons involved in alternative splicing between the control and treatment exhibit little difference in DNA methylation between control and treatment (Table 2).
Conclusion
Here we characterized, for the first time, global gene expression and DNA methylation patterns associated with acute viral infection in an insect species, the honey bee. We demonstrated that several transcriptionally regulated genes are associated with previously identified viral response pathways in insects, including the JAK-STAT pathway, the RNAi pathway, and transcriptional pausing. Note that these last two pathways have not previously been observed to be transcriptionally upregulated in response to viral infection in honey bees, and thus represent intriguing new candidates to consider in the development of viral-resistant honey bee strains or management techniques to mitigate the impact of viruses. We found very little overlap of our gene expression patterns with those identified in previous studies of immunostimulation of adult bees with bacteria[69] and microsprodia[70], and though overlap was greater with studies of viral infection[24,27] (particularly chronic infection of older honey bees with IAPV[21]), these results suggest that immune responses can vary substantially across pathogens, parasites and life stages. We have also demonstrated that epigenetic processes, namely DNA methylation, are altered in response to viral infection in honey bees. Changes in DNA methylation patterns have not previously been described to be an antiviral response in insects, but this is likely due to the fact that most insect molecular immunity studies involve D. melanogaster, which does not have a functional de novo DNA methylation system and has very low global DNA methylation levels. Several of the differentially methylated genes identified in our study are associated with antiviral responses in humans, but have not been associated with immune responses in insects, and thus these represent key genes for future studies of insect immunity. There was no significant overlap in the sets of differentially expressed and differentially methylated genes, and these sets of genes have distinct characteristics (genes that were differentially methylated were longer, more methylated, and less expressed than genes whose methylation patterns did not change). The lack of a significant overlap between differentially expressed and differentially methylated genes suggested honey bees, and perhaps other insects, possess parallel transcriptomic and epigenomic response pathways to viral infection.
Supporting Information
Zdroje
1. Potts SG, Biesmeijer JC, Kremen C, Neumann P, Schweiger O, Kunin WE. Global pollinator declines: trends, impacts and drivers. Trends Ecol Evol. Elsevier Ltd; 2010;25: 345–53. doi: 10.1016/j.tree.2010.01.007 20188434
2. Steinhauer N a, Rennich K, Wilson ME, Caron DM, Lengerich EJ, Pettis JS, et al. A national survey of managed honey bee 2012–2013 annual colony losses in the USA: results from the Bee Informed Partnership. J Apic Res. 2014;53: 1–18. doi: 10.3896/IBRA.1.53.1.01
3. Klein A-M, Vaissière BE, Cane JH, Steffan-Dewenter I, Cunningham S a, Kremen C, et al. Importance of pollinators in changing landscapes for world crops. Proc Biol Sci. 2007;274: 303–13. doi: 10.1098/rspb.2006.3721 17164193
4. Yang X, Cox-Foster DL. Impact of an ectoparasite on the immunity and pathology of an invertebrate: evidence for host immunosuppression and viral amplification. Proc Natl Acad Sci U S A. 2005;102: 7470–5. doi: 10.1073/pnas.0501860102 15897457
5. DeGrandi-Hoffman G, Chen Y, Huang E, Huang MH. The effect of diet on protein concentration, hypopharyngeal gland development and virus load in worker honey bees (Apis mellifera L.). J Insect Physiol. Elsevier Ltd; 2010;56: 1184–91. doi: 10.1016/j.jinsphys.2010.03.017 20346950
6. Di Prisco G, Cavaliere V, Annoscia D, Varricchio P, Caprio E, Nazzi F, et al. Neonicotinoid clothianidin adversely affects insect immunity and promotes replication of a viral pathogen in honey bees. Proc Natl Acad Sci U S A. 2013;110: 18466–71. doi: 10.1073/pnas.1314923110 24145453
7. Chen YP, Siede R. Honey bee viruses. Adv Virus Res. 2007;70: 33–80. doi: 10.1016/S0065-3527(07)70002-7 17765703
8. Runckel C, Flenniken ML, Engel JC, Ruby JG, Ganem D, Andino R, et al. Temporal analysis of the honey bee microbiome reveals four novel viruses and seasonal prevalence of known viruses, Nosema, and Crithidia. PLoS One. 2011;6: e20656. doi: 10.1371/journal.pone.0020656 21687739
9. Bailey L, Ball B V. Honey Bee Pathology. 2nd ed. London: Academic Press INC.; 1991.
10. Iqbal J, Mueller U. Virus infection causes specific learning deficits in honeybee foragers. Proc Biol Sci. 2007;274: 1517–21. doi: 10.1098/rspb.2007.0022 17439851
11. Maori E, Paldi N, Shafir S, Kalev H, Tsur E, Glick E, et al. IAPV, a bee-affecting virus associated with Colony Collapse Disorder can be silenced by dsRNA ingestion. Insect Mol Biol. 2009;18: 55–60. doi: 10.1111/j.1365-2583.2009.00847.x 19196347
12. Dainat B, Evans JD, Chen YP, Gauthier L, Neumann P. Predictive markers of honey bee colony collapse. PLoS One. 2012;7: e32151. doi: 10.1371/journal.pone.0032151 22384162
13. Genersch E, von Der Ohe W, Kaatz H, Schroeder A, Otten C, Buchler R, et al. The German bee monitoring project: a long term study to understand periodically high winter losses of honey bee colonies. Apidologie. 2010;41: 332–352.
14. Van Dooremalen C, Gerritsen L, Cornelissen B, van der Steen JJM, van Langevelde F, Blacquière T. Winter survival of individual honey bees and honey bee colonies depends on level of Varroa destructor infestation. PLoS One. 2012;7: e36285. doi: 10.1371/journal.pone.0036285 22558421
15. De Miranda JR, Cordoni G, Budge G. The Acute bee paralysis virus-Kashmir bee virus-Israeli acute paralysis virus complex. J Invertebr Pathol. Elsevier Inc.; 2010;103 Suppl: S30–47. doi: 10.1016/j.jip.2009.06.014 19909972
16. Boncristiani HF, Evans JD, Chen Y, Pettis J, Murphy C, Lopez DL, et al. In vitro infection of pupae with Israeli acute paralysis virus suggests disturbance of transcriptional homeostasis in honey bees (Apis mellifera). PLoS One. 2013;8: e73429. doi: 10.1371/journal.pone.0073429 24039938
17. Cox-Foster DL, Conlan S, Holmes EC, Palacios G, Evans JD, Moran N a, et al. A metagenomic survey of microbes in honey bee colony collapse disorder. Science. 2007;318: 283–7. doi: 10.1126/science.1146498 17823314
18. Hou C, Rivkin H, Slabezki Y, Chejanovsky N. Dynamics of the presence of israeli acute paralysis virus in honey bee colonies with colony collapse disorder. Viruses. 2014;6: 2012–27. doi: 10.3390/v6052012 24800677
19. Vanengelsdorp D, Evans JD, Saegerman C, Mullin C, Haubruge E, Nguyen BK, et al. Colony collapse disorder: a descriptive study. PLoS One. 2009;4: e6481. doi: 10.1371/journal.pone.0006481 19649264
20. Cornman RS, Tarpy DR, Chen Y, Jeffreys L, Lopez D, Pettis JS, et al. Pathogen webs in collapsing honey bee colonies. PLoS One. 2012;7: e43562. doi: 10.1371/journal.pone.0043562 22927991
21. Chen YP, Pettis JS, Corona M, Chen WP, Li CJ, Spivak M, et al. Israeli Acute Paralysis Virus: Epidemiology, Pathogenesis and Implications for Honey Bee Health. PLoS Pathog. 2014;10: e1004261. doi: 10.1371/journal.ppat.1004261 25079600
22. Merkling SH, van Rij RP. Beyond RNAi: antiviral defense strategies in Drosophila and mosquito. J Insect Physiol. 2013;59: 159–70. doi: 10.1016/j.jinsphys.2012.07.004 22824741
23. Wang X, Aliyari R, Li W-X, Li H-W, Kim K, Carthew R, et al. RNA interference directs innate immunity against viruses in adult Drosophila. Science. 2006;312: 452–4. doi: 10.1126/science.1125694 16556799
24. Flenniken ML, Andino R. Non-specific dsRNA-mediated antiviral response in the honey bee. PLoS One. 2013;8: e77263. doi: 10.1371/journal.pone.0077263 24130869
25. Liu X, Zhang Y, Yan X, Han R. Prevention of Chinese sacbrood virus infection in Apis cerana using RNA interference. Curr Microbiol. 2010;61: 422–8. doi: 10.1007/s00284-010-9633-2 20379718
26. Desai SD, Eu Y-J, Whyard S, Currie RW. Reduction in deformed wing virus infection in larval and adult honey bees (Apis mellifera L.) by double-stranded RNA ingestion. Insect Mol Biol. 2012;21: 446–55. doi: 10.1111/j.1365-2583.2012.01150.x 22690671
27. Ryabov E V, Wood GR, Fannon JM, Moore JD, Bull JC, Chandler D, et al. A virulent strain of deformed wing virus (DWV) of honeybees (Apis mellifera) prevails after Varroa destructor-mediated, or in vitro, transmission. PLoS Pathog. 2014;10: e1004230. doi: 10.1371/journal.ppat.1004230 24968198
28. Finnegan EJ, Genger RK, Peacock WJ, Dennis ES. Dna methylation in plants. Annu Rev Plant Physiol Plant Mol Biol. 1998; 223–247.
29. Goodwin J, Chapman K, Swaney S, Parks TD, Wernsman E, Dougherty W. Genetic and Biochemical Dissection of Transgenic RNA-Mediated Virus Resistance. Plant Cell. 1996;8: 95–105. 8597662
30. Robertson KD. DNA methylation and human disease. Nat Rev Genet. 2005;6: 597–610. doi: 10.1038/nrg1655 16136652
31. Fang J, Hao Q, Liu L, Li Y, Wu J, Huo X, et al. Epigenetic changes mediated by microRNA miR29 activate cyclooxygenase 2 and lambda-1 interferon production during viral infection. J Virol. 2012;86: 1010–20. doi: 10.1128/JVI.06169-11 22072783
32. Kathiria P, Sidler C, Golubov A, Kalischuk M, Kawchuk LM, Kovalchuk I. Tobacco mosaic virus infection results in an increase in recombination frequency and resistance to viral, bacterial, and fungal pathogens in the progeny of infected tobacco plants. Plant Physiol. 2010;153: 1859–70. doi: 10.1104/pp.110.157263 20498336
33. Li W, Sun W, Liu L, Yang F, Li Y, Chen Y, et al. IL-32: a host proinflammatory factor against influenza viral replication is upregulated by aberrant epigenetic modifications during influenza A virus infection. J Immunol. 2010;185: 5056–65. doi: 10.4049/jimmunol.0902667 20889550
34. Mukherjee S, Vipat VC, Chakrabarti AK. Infection with influenza A viruses causes changes in promoter DNA methylation of inflammatory genes. Influenza Other Respi Viruses. 2013;1: 1–8. doi: 10.1111/irv.12127
35. Takacs M, Segesdi J, Banati F, Koroknai A, Wolf H, Niller HH, et al. The importance of epigenetic alterations in the development of epstein-barr virus-related lymphomas. Mediterr J Hematol Infect Dis. 2009;1: e2009012. doi: 10.4084/MJHID.2009.012 21416002
36. Tang B, Zhao R, Sun Y, Zhu Y, Zhong J, Zhao G, et al. Interleukin-6 expression was regulated by epigenetic mechanisms in response to influenza virus infection or dsRNA treatment. Mol Immunol. Elsevier Ltd; 2011;48: 1001–8. doi: 10.1016/j.molimm.2011.01.003 21353307
37. Urieli-Shoval S, Gruenbaum Y, Sedat J, Razin A. The absence of detectable methylated bases in Drosophila melanogaster DNA. 1982;146: 148–152. 6814955
38. Capuano F, Kok R, Blom HJ, Ralser M. Cytosine DNA Methylation Is Found in Drosophila melanogaster but Absent in Saccharomyces cerevisiae, Schizosaccharomyces pombe, and Other Yeast Species. Anal Chem. 2014; 3697–3702. doi: 10.1021/ac500447w 24640988
39. Wang Y, Jorda M, Jones PL, Maleszka R, Ling X, Hugh M, et al. Functional System CpG Methylation in a Social Insect. 2006;314: 645–647. 17068262
40. Suzuki MM, Bird A. DNA methylation landscapes: provocative insights from epigenomics. Nat Rev Genet. 2008;9: 465–76. doi: 10.1038/nrg2341 18463664
41. Zemach A, McDaniel IE, Silva P, Zilberman D. Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science. 2010;328: 916–9. doi: 10.1126/science.1186366 20395474
42. Klose RJ, Bird AP. Genomic DNA methylation: the mark and its mediators. Trends Biochem Sci. 2006;31: 89–97. doi: 10.1016/j.tibs.2005.12.008 16403636
43. Foret S, Kucharski R, Pellegrini M, Feng S, Jacobsen SE, Robinson GE, et al. DNA methylation dynamics, metabolic fluxes, gene splicing, and alternative phenotypes in honey bees. Proc Natl Acad Sci U S A. 2012;109: 4968–73. doi: 10.1073/pnas.1202392109 22416128
44. Kucharski R, Maleszka J, Foret S, Maleszka R. Nutritional control of reproductive status in honeybees via DNA methylation. Science. 2008;319: 1827–30. doi: 10.1126/science.1153069 18339900
45. Bird AP. Gene number, noise reduction and biological complexity. Trends Genet. 1995;11: 94–100. doi: 10.1016/S0168-9525(00)89009-5 7732579
46. Huh I, Zeng J, Park T, Yi S V. DNA methylation and transcriptional noise. Epigenetics Chromatin. 2013;6: 9. doi: 10.1186/1756-8935-6-9 23618007
47. Fussnecker BL, McKenzie AM, Grozinger CM. cGMP modulates responses to queen mandibular pheromone in worker honey bees. J Comp Physiol A Sensory, Neural, Behav Physiol. 2013;197: 939–948. doi: 10.1007/s00359-011-0654-5 cGMP
48. Li Z, Chen Y, Zhang S, Chen S, Li W, Yan L, et al. Viral infection affects sucrose responsiveness and homing ability of forager honey bees, Apis mellifera L. PLoS One. 2013;8: e77354. doi: 10.1371/journal.pone.0077354 24130876
49. Hunter W, Ellis J, Vanengelsdorp D, Hayes J, Westervelt D, Glick E, et al. Large-scale field application of RNAi technology reducing Israeli acute paralysis virus disease in honey bees (Apis mellifera, Hymenoptera: Apidae). PLoS Pathog. 2010;6: e1001160. doi: 10.1371/journal.ppat.1001160 21203478
50. Richard F-J, Tarpy DR, Grozinger CM. Effects of insemination quantity on honey bee queen physiology. PLoS One. 2007;2: e980. doi: 10.1371/journal.pone.0000980 17912357
51. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30: 207–10. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=99122&tool=pmcentrez&rendertype=abstract 11752295
52. Elsik CG, Worley KC, Bennett AK, Beye M, Camara F, Childers CP, et al. Finding the missing honey bee genes: lessons learned from a genome upgrade. BMC Genomics. BMC Genomics; 2014;15: 86. doi: 10.1186/1471-2164-15-86 24479613
53. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25: 1105–11. doi: 10.1093/bioinformatics/btp120 19289445
54. Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11: R25. doi: 10.1186/gb-2010-11-3-r25 20196867
55. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. BioMed Central Ltd; 2010;11: R106. doi: 10.1186/gb-2010-11-10-r106 20979621
56. Huang DW, Sherman BT, Lempicki R a. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4: 44–57. doi: 10.1038/nprot.2008.211 19131956
57. Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6: e21800. doi: 10.1371/journal.pone.0021800 21789182
58. Xi Y, Li W. BSMAP: whole genome bisulfite sequence MAPping program. BMC Bioinformatics. 2009;10: 232. doi: 10.1186/1471-2105-10-232 19635165
59. Lyko F, Foret S, Kucharski R, Wolf S, Falckenhayn C, Maleszka R. The honey bee epigenomes: differential methylation of brain DNA in queens and workers. PLoS Biol. 2010;8: e1000506. doi: 10.1371/journal.pbio.1000506 21072239
60. Lister R, Pelizzola M, Dowen RH, Hawkins RD, Hon G, Tonti-Filippini J, et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. Nature Publishing Group; 2009;462: 315–22. doi: 10.1038/nature08514 19829295
61. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B …. 1995;
62. Elango N, Hunt BG, Goodisman M a D, Yi S V. DNA methylation is widespread and associated with differential gene expression in castes of the honeybee, Apis mellifera. Proc Natl Acad Sci U S A. 2009;106: 11206–11. doi: 10.1073/pnas.0900301106 19556545
63. Lockett G a, Kucharski R, Maleszka R. DNA methylation changes elicited by social stimuli in the brains of worker honey bees. Genes Brain Behav. 2012;11: 235–42. doi: 10.1111/j.1601-183X.2011.00751.x 22098706
64. Xu J, Grant G, Sabin LR, Gordesky-Gold B, Yasunaga A, Tudor M, et al. Transcriptional pausing controls a rapid antiviral innate immune response in Drosophila. Cell Host Microbe. Elsevier Inc.; 2012;12: 531–43. doi: 10.1016/j.chom.2012.08.011 23084920
65. Kingsolver MB, Huang Z, Hardy RW. Insect antiviral innate immunity: pathways, effectors, and connections. J Mol Biol. Elsevier Ltd; 2013;425: 4921–36. doi: 10.1016/j.jmb.2013.10.006 24120681
66. Mukherjee S, Hanley K a. RNA interference modulates replication of dengue virus in Drosophila melanogaster cells. BMC Microbiol. 2010;10: 127. doi: 10.1186/1471-2180-10-127 20420715
67. Niu J, Meeus I, Cappelle K, Piot N, Smagghe G. The immune response of the small interfering RNA pathway in the defense against bee viruses. Curr Opin Insect Sci. Elsevier Inc; 2014; 22–27. doi: 10.1016/j.cois.2014.09.014 25401084
68. Evans JD, Aronstein K, Chen YP, Hetru C, Imler J-L, Jiang H, et al. Immune pathways and defence mechanisms in honey bees Apis mellifera. Insect Mol Biol. 2006;15: 645–56. doi: 10.1111/j.1365-2583.2006.00682.x 17069638
69. Richard F-J, Holt HL, Grozinger CM. Effects of immunostimulation on social behavior, chemical communication and genome-wide gene expression in honey bee workers (Apis mellifera). BMC Genomics. BMC Genomics; 2012;13: 558. doi: 10.1186/1471-2164-13-558 23072398
70. Holt HL, Aronstein K a, Grozinger CM. Chronic parasitization by Nosema microsporidia causes global expression changes in core nutritional, metabolic and behavioral pathways in honey bee workers (Apis mellifera). BMC Genomics. 2013;14: 799. doi: 10.1186/1471-2164-14-799 24245482
71. Sarda S, Zeng J, Hunt BG, Yi S V. The evolution of invertebrate gene body methylation. Mol Biol Evol. 2012;29: 1907–16. doi: 10.1093/molbev/mss062 22328716
72. Wang Y, Leung FCC. In silico prediction of two classes of honeybee genes with CpG deficiency or CpG enrichment and sorting according to gene ontology classes. J Mol Evol. 2009;68: 700–5. doi: 10.1007/s00239-009-9244-3 19466376
73. Foret S, Kucharski R, Pittelkow Y, Lockett G a, Maleszka R. Epigenetic regulation of the honey bee transcriptome: unravelling the nature of methylated genes. BMC Genomics. 2009;10: 472. doi: 10.1186/1471-2164-10-472 19828049
74. Wang X, Wheeler D, Avery A, Rago A, Choi J-H, Colbourne JK, et al. Function and evolution of DNA methylation in Nasonia vitripennis. PLoS Genet. 2013;9: e1003872. doi: 10.1371/journal.pgen.1003872 24130511
75. Zeng J, Yi S V. DNA methylation and genome evolution in honeybee: gene length, expression, functional enrichment covary with the evolutionary signature of DNA methylation. Genome Biol Evol. 2010;2: 770–80. doi: 10.1093/gbe/evq060 20924039
76. Hellman A, Chess A. Gene body-specific methylation on the active X chromosome. Science. 2007;315: 1141–3. doi: 10.1126/science.1136352 17322062
77. Jones P a. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. Nature Publishing Group; 2012;13: 484–92. doi: 10.1038/nrg3230 22641018
78. Duret L, Mouchiroud D. Expression pattern and, surprisingly, gene length shape codon usage in Caenorhabditis, Drosophila, and Arabidopsis. Proc Natl Acad Sci U S A. 1999;96: 4482–7. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=16358&tool=pmcentrez&rendertype=abstract 10200288
79. Park J, Xu K, Park T, Yi S V. What are the determinants of gene expression levels and breadths in the human genome? Hum Mol Genet. 2012;21: 46–56. doi: 10.1093/hmg/ddr436 21945885
80. Ehrlich LS, Medina GN, Khan MB, Powell MD, Mikoshiba K, Carter C a. Activation of the inositol (1,4,5)-triphosphate calcium gate receptor is required for HIV-1 Gag release. J Virol. 2010;84: 6438–51. doi: 10.1128/JVI.01588-09 20427533
81. Hrincius ER, Dierkes R, Anhlan D, Wixler V, Ludwig S, Ehrhardt C. Phosphatidylinositol-3-kinase (PI3K) is activated by influenza virus vRNA via the pathogen pattern receptor Rig-I to promote efficient type I interferon production. Cell Microbiol. 2011;13: 1907–19. doi: 10.1111/j.1462-5822.2011.01680.x 21899695
82. Liu Z, Tian Y, Machida K, Lai MMC, Luo G, Foung SKH, et al. Transient activation of the PI3K-AKT pathway by hepatitis C virus to enhance viral entry. J Biol Chem. 2012;287: 41922–30. doi: 10.1074/jbc.M112.414789 23095753
83. Contreras X, Mzoughi O, Gaston F, Peterlin MB, Bahraoui E. Protein kinase C-delta regulates HIV-1 replication at an early post-entry step in macrophages. Retrovirology. BioMed Central Ltd; 2012;9: 37. doi: 10.1186/1742-4690-9-37 22554282
84. Xu W, Santini P a, Sullivan JS, He B, Shan M, Ball SC, et al. HIV-1 evades virus-specific IgG2 and IgA responses by targeting systemic and intestinal B cells via long-range intercellular conduits. Nat Immunol. Nature Publishing Group; 2009;10: 1008–17. doi: 10.1038/ni.1753 19648924
85. Strong JE, Coffey MC, Tang D, Sabinin P, Lee PW. The molecular basis of viral oncolysis: usurpation of the Ras signaling pathway by reovirus. EMBO J. 1998;17: 3351–62. doi: 10.1093/emboj/17.12.3351 9628872
86. Lockett G a, Helliwell P, Maleszka R. Involvement of DNA methylation in memory processing in the honey bee. Neuroreport. 2010;21: 812–6. doi: 10.1097/WNR.0b013e32833ce5be 20571459
87. Park J, Peng Z, Zeng J, Elango N, Park T, Wheeler D, et al. Comparative analyses of DNA methylation and sequence evolution using Nasonia genomes. Mol Biol Evol. 2011;28: 3345–54. doi: 10.1093/molbev/msr168 21693438
88. Shukla S, Kavak E, Gregory M, Imashimizu M, Shutinoski B, Kashlev M, et al. CTCF-promoted RNA polymerase II pausing links DNA methylation to splicing. Nature. Nature Publishing Group; 2011;479: 74–9. doi: 10.1038/nature10442 21964334
89. Herb BR, Wolschin F, Hansen KD, Aryee MJ, Langmead B, Irizarry R, et al. Reversible switching between epigenetic states in honeybee behavioral subcastes. Nat Neurosci. Nature Publishing Group; 2012;15: 1371–3. doi: 10.1038/nn.3218 22983211
90. Castorena KM, Stapleford KA, Miller DJ. targeted functional genetic analyses in cultured Drosophila cells highlight the role of glycerophospholipid metabolism in Flock House virus RNA replication. 2010;
Štítky
Hygiena a epidemiológia Infekčné lekárstvo LaboratóriumČlánok vyšiel v časopise
PLOS Pathogens
2015 Číslo 3
- Parazitičtí červi v terapii Crohnovy choroby a dalších zánětlivých autoimunitních onemocnění
- Očkování proti virové hemoragické horečce Ebola experimentální vakcínou rVSVDG-ZEBOV-GP
- Koronavirus hýbe světem: Víte jak se chránit a jak postupovat v případě podezření?
Najčítanejšie v tomto čísle
- Bacterial Immune Evasion through Manipulation of Host Inhibitory Immune Signaling
- Antimicrobial-Induced DNA Damage and Genomic Instability in Microbial Pathogens
- Is Antigenic Sin Always “Original?” Re-examining the Evidence Regarding Circulation of a Human H1 Influenza Virus Immediately Prior to the 1918 Spanish Flu
- An 18 kDa Scaffold Protein Is Critical for Biofilm Formation