The invention is in the field of virology and relates to the deformed wing virus (DWV) and the use of inhibitors of DWV to prevent and/or treat DWV infection in bees.
Honeybees are of great importance to the global economy, far surpassing their contribution in terms of honey production. In all, 52 of the world's 115 leading agricultural crops rely on honeybee pollination to some extent. These crops represent approximately 35% of the human diet.
Honeybee numbers have decreased in recent years. One factor contributing to this decrease is disease. In particular, honeybees are susceptible to a host of picorna-like viruses, including the closely related Acute Bee Paralysis Virus (ABPV), Kashmir Bee Virus (KBV), and Israeli Acute Paralysis Virus (IAPV). Viral infection can have a devastating effect on the bee population, resulting in high mortality rates and a decrease in the number of bees and colonies.
The mite, Varroa destructor, which feeds on honeybee haemolymph, originated in Asia and arrived in the UK in the 1980s. The Varroa mites are known to transmit bee viruses. Research has identified several viral pathogens of honeybees that are transmitted by Varroa mites, including the deformed wing virus (DWV). The DWV is so called because the virus causes the appearance of honeybees with characteristic wing deformities together with other developmental defects such as abdominal stunting within bee colonies. Infestation of honeybee colonies with Varroa results in a dramatic increase in DWV levels.
The present inventors have identified a novel strain of deformed wing virus (DWV) that is surprisingly predominant in bees infested with Varroa mites. This novel strain of the invention comprises a recombinant genome containing the Varroa Desctructor Virus (VDV-1)-derived structural genes and the DWV-derived non-structural genes. The inventors have identified two predominant genomic sequences for the DWV strain. These two sequences differ only in the 5′ non-coding regions (NCRs). Thus, the regions of the two genomic sequences which code for all the viral proteins are identical. The two genomic sequences are set out in SEQ ID NOs: 1 and 2.
The novel strain of the invention is found at high concentrations within individual bees in Varroa-infested colonies, with other, highly divergent, strains of DWV present at much lower concentrations. Typically the newly identified strain is present at 1,000-10,000 times the concentration of the other DWV strains. This is in contrast to non-Varroa infested colonies, where, although a high diversity of DWV strains is observed, individual bees exhibit a much lower viral load. The particular strain of DWV identified by the inventors can be used in diagnostics to identify those colonies at risk of Varroa-transmitted deformed wing disease. Inhibitors of the particular strain may be used in the treatment and/or prevention of DWV-induced disease.
Also, transgenic bees can be generated which are resistant to the particular strain of DWV and hence resistant to Varroa mite-induced deformed wing disease.
Accordingly, the present invention provides a polynucleotide comprising (a) the sequence shown in SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62, (b) a variant sequence having at least 98% homology to SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 based on nucleotide identity over its entire length or (c) a sequence that is complementary to the sequence in (a) or (b).
The invention also provides an oligonucleotide which specifically hybridises to a part of the polynucleotide of the invention.
The invention also provides an oligonucleotide which comprises 50 or fewer consecutive nucleotides from a polynucleotide of the invention.
The invention also provides a polynucleotide comprising at least two oligonucleotides of the invention.
The invention also provides an isolated strain of deformed wing virus (DWV) which comprises the varroa destructor virus 1 (VDV-1) capsid proteins (CP) and the DWV non-structural proteins (NS). This isolated strain of deformed wing virus (DWV) is referred to herein as the isolated DWV strain of the invention, or the isolated strain of the invention. This isolated strain may comprise a polynucleotide comprising (a) the sequence shown in SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62, (b) a variant sequence having at least 98% homology to SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 based on nucleotide identity over its entire length or (c) a sequence that is complementary to the sequence in (a) or (b).
The invention also provides an antibody which specifically binds to an isolated strain of the invention.
The invention also provides a method of determining the presence or absence of an isolated strain of the invention in a sample, comprising (a) contacting the sample with an oligonucleotide of the invention or an antibody of the invention and (b) detecting specific hybridisation of the oligonucleotide or specific binding of the antibody and thereby determining the presence or absence of the isolated strain.
The invention also provides a method of determining the concentration of an isolated strain of the invention in a sample, comprising (a) contacting the sample with an oligonucleotide of the invention or an antibody of the invention and (b) detecting the amount of specific hybridisation of the oligonucleotide or the amount of specific binding of the antibody and thereby determining the concentration of the isolated strain.
The invention further provides a vector comprising a polynucleotide of the invention or an oligonucleotide of the invention, wherein said polynucleotide or oligonucleotide is operably linked to a promoter.
The invention further provides a composition comprising a polynucleotide of the invention, an oligonucleotide of the invention, an antibody of the invention and/or a vector of the invention and a delivery vehicle.
The invention also provides a method of treating or preventing in a bee or bee colony an infection with an isolated strain of the invention, comprising contacting the bee or bee colony with an inhibitor of the isolated strain.
The invention also a method of treating or preventing deformed wing disease in a Varroa mite-infested bee or bee colony, comprising contacting the bee or bee colony with an inhibitor of an isolated strain of the invention.
The invention also provides a method of diagnosing in a bee or bee colony infection with an isolated strain of the invention, comprising determining the presence or absence of the isolated strain, wherein the presence of the isolated strain is indicative of the presence of infection with the isolated strain and wherein the absence of the isolated strain is indicative of the absence of infection with the isolated strain.
The invention also provides a method of diagnosing deformed wing disease in a Varroa mite-infested bee or bee colony, comprising determining the presence or absence of an isolated strain of the invention, wherein the presence of the isolated strain is indicative of the presence of deformed wing and wherein the absence of the isolated strain is indicative of the absence of deformed wing disease.
The invention also provides a transgenic bee that is resistant to infection by an isolated strain of the invention. The invention also provides a transgenic bee that is resistant to Varroa-mite induced deformed wing disease, wherein at least one cell of the bee expresses an oligonucleotide of the invention.
The invention further provides a method of generating a transgenic queen bee that is resistant to infection by an isolated strain of the invention, comprising (a) incorporating an oligonucleotide of the invention or a polynucleotide of the invention into the genome of one or more bee germ cells; and (b) generating the queen bee from said one or more germ cells.
The invention also provides a method of generating a transgenic bee of the invention, comprising using a transgenic queen bee of the invention or a transgenic queen bee produced using a method of the invention to generate the bee.
The invention also provides a method of preventing in a bee or bee colony an infection with an isolated strain of the invention, comprising using a transgenic queen bee of the invention or a transgenic queen bee produced using a method of the invention to generate the bee or bee colony.
The invention further provides a method of preventing deformed wing disease in a Varroa mite-infested bee or bee colony, comprising using a transgenic queen bee of the invention or a transgenic queen bee produced using a method of the invention to generate the bee or bee colony.
The invention also provides a method of producing a bee or bee colony that is resistant to infection from an isolated strain of the invention, comprising using a transgenic queen bee of the invention or a transgenic queen bee produced using a method of the invention to generate the bee or bee colony.
The invention also provides one or more bee germ cells comprising an oligonucleotide of the invention, a polynucleotide of the invention or a vector of the invention.
The invention also provides a bee colony wherein at least 50% of the bees within the colony are resistant to infection by an isolated strain of the invention.
The invention further provides a bee colony wherein at least 50% of the bees within the colony are resistant to Varroa mite-induced deformed wing disease.
SEQ ID NOs: 1 and 2 are the genomic sequences of the DWV strain of the invention. As discussed herein, the regions encoding the structural proteins and non-structural proteins are identical between SEQ ID NOs: 1 and 2. Only the sequence of the 5′ non-coding region (NCR) differs between SEQ ID NOs: 1 and 2.
SEQ ID NO: 3 is the coding sequence for the Varroa-destructor virus-1 (VDV-1) capsid proteins (CP).
SEQ ID NO: 4 is the coding sequence for the DWV non-structural proteins (NS).
SEQ ID NOs: 5 and 6 are the sequences comprised in the preferred polynucleotides of the invention. Preferred oligonucleotides of the invention may be derived from these sequences.
SEQ ID NO: 7 is the sequence of a preferred transducing vector, pTYF-mCMW/SYN-EGFP.
SEQ ID NO: 8 is the sequence of the newly-identified honeybee heat shock protein 70 (hsp70) promoter.
SEQ ID NO: 9 is the sequence of a preferred honeybee transcription terminator.
SEQ ID NO: 10 is the sequence of a transducing vector derived from pTYF-mCMW/SYN-EGFP which comprises a luciferase reporter gene operably linked to the honeybee hsp70 promoter of SEQ ID NO: 7.
SEQ ID NOs: 11 to 51 are sequences of exemplary oligonucleotides of the invention (Table 1).
SEQ ID NO: 52 is the sequence of Clone-202-complete-KJ437447.
SEQ ID NO: 53 is the sequence of Clone-complete-204
SEQ ID NO: 54 is the sequence of Clone-partial-21
SEQ ID NO: 55 is the sequence of Clone-23-partial
SEQ ID NO: 56 is the sequence of NGS-Library-E7
SEQ ID NO: 57 is the sequence of NGS-Library-E8
SEQ ID NO: 58 is the sequence of NGS-Library-E10
SEQ ID NO: 59 is the sequence of NGS-Library-E11
SEQ ID NO: 60 is the sequence of NGS-Library-INJ4
SEQ ID NO: 61 is the sequence of NGS-Library-INJ5
SEQ ID NO: 62 is the sequence of NGS-Library-INJ6
SEQ ID NOs: 52 to 55 are PCR amplified cDNAs that have been sequenced by standard Sanger sequencing.
SEQ ID NOs: 56 to 62 are consensus sequences determined from next generation sequencing, either of naturally infected bees (SEQ ID NOs: 56 to 59) or honeybee pupae injected with virus in the laboratory (SEQ ID NOs: 60 to 62).
SEQ ID NOs: 63 to 72 are sequences of exemplary oligonucleotides of the invention (Table 5).
SEQ ID NO: 73 is the sequence of a preferred transducing vector, pTYF-mCMW/SYN-EGFP.
It is to be understood that different applications of the disclosed products and methods may be tailored to the specific needs in the art. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments of the invention only, and is not intended to be limiting.
In addition as used in this specification and the appended claims, the singular forms “a”, “an”, and “the” include plural referents unless the content clearly dictates otherwise. Thus, for example, reference to “a polynucleotide” includes two or more such polynucleotides, reference to “an oligonucleotide” includes two or more such oligonucleotides, reference to “a virus” includes two or more such viruses, and the like.
All publications, patents and patent applications cited herein, whether supra or infra, are hereby incorporated by reference in their entirety.
Deformed wing virus (DWV), Varroa destructor virus-1 (VDV-1) (Genbank Accession No. NC—006494, Version No. NC—006494.1 GI:56121875, Ongus et al 2004 J. Gen. Virol. 85: 3747-3755), their recombinants, and Kakugo virus (KV) (Genbank Accession No. NC—005876, Version No. NC—005876.1 GI:47177088, Fujiyuki et al 2004 J. Virol. 78(3): 1093-1100) are closely related picorna-like RNA viruses infecting the European honeybee (Apis mellifera). Some of these viruses also infect the ectoparasitic mite Varroa destructor. The overall nucleotide homology between the viruses related to DWV is no less than 84%. Therefore, for the purposes of the present invention, these viruses are considered as strains of the same virus, DWV.
Picorna-like viruses, including DWV, possess genomes that are essentially ‘modular’. These modular genomes consist of 4 modules. From the 5′ end (i.e. the direction of translation and the direction that they are always represented in print) the four modules are:
The individual modules can be recombined to create recombinant viral genomes.
Recombinant viruses must have all four modules in the correct order, and the modules must be compatible with each other. Module compatibility is achieved through selection, as only functional viral genomes will allow viral growth and the generation of progeny viruses.
DWV is a single-stranded, positive sense RNA virus. Prior to the present invention, DWV has been relatively poorly characterised. The DWV genome is approximately 10150 nucleotides in length and has a modular architecture as discussed above. In particular, the four DWV modules have been previously defined as follows (Genbank Accession No. NC—004830; Version No. NC—004830.2 GI:71480055, Lanzi et al. 2006 J. Virol. 80(10): 4998-5009):
In the DWV genome there is an additional protein, the L protein (a non-structural protein) which is coded before the first of the structural proteins. Although other strains of DWV are known in the art, the present strain is newly identified by the inventors.
DWV is present in the majority of honeybee colonies, but in the absence of its ectoparasite, the mite V. destructor which feeds on the bee haemolymph, DWV generally causes asymptomatic infection and accumulates to low concentrations in infected bees. Conversely, in Varroa-infested colonies bees show impaired development and increased mortality associated with very high DWV concentrations.
The present inventors have investigated the effect of Varroa infestation on the DWV strain composition and DWV concentration in DWV infected colonies. Pupae infected orally but remaining mite-free showed significant changes in DWV strain composition, but only a slight increase in overall DWV concentrations. Varroa-infestation resulted in further changes in DWV strain compositions with concentrations of DWV showing bimodal distribution, being either similar to the orally infested or 1000-times higher. Surprisingly, the present inventors found that in bees with the very highest concentrations of DWV-like viruses, which comprise about 80% of the Varroa mite-infested bee pupae, there was a single recombinant strain of DWV/VDV. The predominant (up to 99.9% of the virus genomes present) strain in the Varroa-infested bees had recombinant genomes containing the VDV-1-derived structural genes and the DWV-derived non-structural genes.
Accordingly, the present invention provides an isolated strain of DWV comprising a recombinant genome containing the VDV-1-derived structural genes and the DWV-derived non-structural genes. In one embodiment, in the isolated strain of the invention the precise recombination junction between the VDV-1 derived sequence and the DWV-derived sequence matches exactly to the junction of the VDV-1-derived structural genes and the DWV-derived non-structural genes. In another embodiment, the precise recombination junction between the VDV-1 derived sequence and the DWV-derived sequence does not matches exactly to the junction of the VDV-1-derived structural genes and the DWV-derived non-structural genes. For example, the recombination junction may lie within 100 nucleotides, within 50 nucleotides, within 20 nucleotides, within 10 nucleotides, within 9 nucleotides, within 8 nucleotides, within 7 nucleotides, within 6 nucleotides, within 5 nucleotides, within 4 nucleotides, within 3 nucleotides, within 2 nucleotides or within 1 nucleotide of the junction of the structural genes and the non-structural genes.
The present inventors have sequenced the genome of the isolated DWV strain and have identified two predominant sequences. The sequences of the isolated strain are set out in SEQ ID NOs: 1 and 2. SEQ ID NO: 1 comprises (1) the DWV 5′ NCR module; (2) a structural protein module encoding the varroa destructor virus 1 (VDV-1) capsid proteins (CP); (3) a non-structural protein module encoding the DWV non-structural proteins (NS); and (4) a 3′ NCR. SEQ ID NO: 2 comprises (1) the VDV-1 5′ NCR module; (2) a structural protein module encoding the VDV-1 CP; (3) a non-structural protein module encoding the DWV NS; and (4) a 3′ NCR.
The VDV-1 CP are encoded by SEQ ID NO: 3, which is present in both SEQ ID NO: 1 and 2. The DWV NS are encoded by SEQ ID NO: 4, which is present in both SEQ ID NO:1 and 2. SEQ ID NO: 3 is found at nucleotides 1118 to 4594 of SEQ ID NO: 1 and at nucleotides 1145 to 4621 of SEQ ID NO: 2. SEQ ID NO: 4 is found at nucleotides 4595 to 9799 of SEQ ID NO: 1 and at nucleotides 4622 to 9826 of SEQ ID NO: 2. Corresponding regions in SEQ ID NO: 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 can be determined using sequence alignments. Methods of aligning sequences are discussed in more detail below.
The sequences of SEQ ID NOs: 1 and 2 differ in their 5′ NCR. The regions of the two genomic sequences which code for all the viral proteins (approximately 90% of the genome) are identical. Thus, the isolated strain of the invention can be thought of as possessing a single defined protein coding region with one of two alternative 5′ NCRs. It is the particular combination of protein coding regions of the isolate strain, i.e. modules 2 and 3 identified above, which give the isolated strain of the invention its unique characteristics. The 5′ NCR is not determinative. Therefore, the strain of the invention may comprise the NCR of either DWV or VDV-1, such that translation of the strain of the invention is mediated by the 5′ NCR of either DWV or VDV-1.
In one embodiment the isolated strain of the invention comprises a polynucleotide comprising a non-structural protein module encoding the DWV NS and/or a structural protein module encoding the VDV-1 CP. In a preferred embodiment the isolated strain comprises a polynucleotide comprising both a non-structural protein module encoding the DWV NS and a structural protein module encoding the VDV-1 CP.
The VDV-1 CP is the capsid protein region of the DWV strain of the invention. The capsid proteins form the outer proteins of the virus particle and are therefore considered to be structural proteins (Lanzi et al. 2006 J. Virol. 80(10): 4998-5009). The VDV-1 CP comprises 4 capsid proteins, VP1, VP2, VP3 and VP4 which are expressed as part of a polyprotein, preceded by protein L.
The DWV NS is the non-structural protein region of the DWV strain of the invention. This region comprises: (1) a helicase that is predicted to have RNA structure unwinding activity; (2) a viral protein genome linked (VPg) protein which covalently attaches to the viral genome during replication; (3) a second protease similar to the 3C protease of other picornaviruses; and (4) an RNA-dependent RNA polymerase (RdRp).
In one embodiment, the isolated strain of the invention comprises a polynucleotide sequence of SEQ ID NO: 3, which encodes the VDV-1 CP and/or a polynucleotide sequence of SEQ ID NO: 4, which encodes the DWV NS.
In a preferred embodiment, the isolated strain of the invention comprises a polynucleotide sequence of the invention as discussed below. The polynucleotide sequence of the isolated strain is preferably SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62. The polynucleotide sequence of the isolated strain may be any of the variant sequences discussed below.
The strain of the invention is “isolated” in the sense that it is isolated (or separated) from its natural state and the molecules with which it is found in nature. The virus may be mixed with other molecules, carriers or diluents, such as those disclosed below, which will not interfere with its intended use. Typically, the isolated virus is not present in a bee or a Varroa mite.
An isolated strain of the virus may be produced by expressing a polynucleotide of the invention are discussed below.
The invention provides a polynucleotide comprising (a) the sequence of the genome of the isolated strain of the invention, namely SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62, (b) a variant sequence having at least 98% homology to SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 based on nucleotide identity over its entire length or (c) a sequence which is complementary to SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 or a variant sequence as defined in (b).
A polynucleotide, such as a nucleic acid, is a polymer comprising two or more nucleotides. The nucleotides can be naturally occurring or artificial.
A nucleotide typically contains a nucleobase, a sugar and at least one linking group, such as a phosphate, 2′O-methyl, 2′ methoxy-ethyl, phosphoramidate, methylphosphonate or phosphorothioate group. The nucleobase is typically heterocyclic. Nucleobases include, but are not limited to, purines and pyrimidines and more specifically adenine (A), guanine (G), thymine (T), uracil (U) and cytosine (C). The sugar is typically a pentose sugar. Nucleotide sugars include, but are not limited to, ribose and deoxyribose. The nucleotide is typically a ribonucleotide or deoxyribonucleotide. The nucleotide typically contains a monophosphate, diphosphate or triphosphate. Phosphates may be attached on the 5′ or 3′ side of a nucleotide.
Nucleotides include, but are not limited to, adenosine monophosphate (AMP), adenosine diphosphate (ADP), adenosine triphosphate (ATP), guanosine monophosphate (GMP), guanosine diphosphate (GDP), guanosine triphosphate (GTP), thymidine monophosphate (TMP), thymidine diphosphate (TDP), thymidine triphosphate (TTP), uridine monophosphate (UMP), uridine diphosphate (UDP), uridine triphosphate (UTP), cytidine monophosphate (CMP), cytidine diphosphate (CDP), cytidine triphosphate (CTP), 5-methylcytidine monophosphate, 5-methylcytidine diphosphate, 5-methylcytidine triphosphate, 5-hydroxymethylcytidine monophosphate, 5-hydroxymethylcytidine diphosphate, 5-hydroxymethylcytidine triphosphate, cyclic adenosine monophosphate (cAMP), cyclic guanosine monophosphate (cGMP), deoxyadenosine monophosphate (dAMP), deoxyadenosine diphosphate (dADP), deoxyadenosine triphosphate (dATP), deoxyguanosine monophosphate (dGMP), deoxyguanosine diphosphate (dGDP), deoxyguanosine triphosphate (dGTP), deoxythymidine monophosphate (dTMP), deoxythymidine diphosphate (dTDP), deoxythymidine triphosphate (dTTP), deoxyuridine monophosphate (dUMP), deoxyuridine diphosphate (dUDP), deoxyuridine triphosphate (dUTP), deoxycytidine monophosphate (dCMP), deoxycytidine diphosphate (dCDP) and deoxycytidine triphosphate (dCTP), 5-methyl-2′-deoxycytidine monophosphate, 5-methyl-2′-deoxycytidine diphosphate, 5-methyl-2′-deoxycytidine triphosphate, 5-hydroxymethyl-2′-deoxycytidine monophosphate, 5-hydroxymethyl-2′-deoxycytidine diphosphate and 5-hydroxymethyl-2′-deoxycytidine triphosphate. The nucleotides are preferably selected from AMP, TMP, GMP, UMP, dAMP, dTMP, dGMP or dCMP.
The nucleotides may contain additional modifications. In particular, suitable modified nucleotides include, but are not limited to, 2′amino pyrimidines (such as 2′-amino cytidine and 2′-amino uridine), 2′-hyrdroxyl purines (such as, 2′-fluoro pyrimidines (such as 2′-fluorocytidine and 2′fluoro uridine), hydroxyl pyrimidines (such as 5′-α-P-borano uridine), 2′-O-methyl nucleotides (such as 2′-O-methyl adenosine, 2′-O-methyl guanosine, 2′-O-methyl cytidine and 2′-O-methyl uridine), 4′-thio pyrimidines (such as 4′-thio uridine and 4′-thio cytidine) and nucleotides have modifications of the nucleobase (such as 5-pentynyl-2′-deoxy uridine, 5-(3-aminopropyl)-uridine and 1,6-diaminohexyl-N-5-carbamoylmethyl uridine).
One or more nucleotides in the polynucleotide can be oxidized or methylated. One or more nucleotides in the polynucleotide may be damaged. For instance, the polynucleotide may comprise a pyrimidine dimer. Such dimers are typically associated with damage by ultraviolet light.
One or more nucleotides in the polynucleotide may be modified, for instance with a label or a tag. The label may be any suitable label which allows the polynucleotide to be detected. Suitable labels include, but are not limited to fluorescent molecules, radioisotopes, e.g. 125I, 35S, enzymes, antibodies, antigens, other polynucleotides and ligands such as biotin.
The nucleotides in the polynucleotide may be attached to each other in any manner. The nucleotides may be linked by phosphate, 2′O-methyl, 2′ methoxy-ethyl, phosphoramidate, methylphosphonate or phosphorothioate linkages. The nucleotides are typically attached by their sugar and phosphate groups as in nucleic acids. The nucleotides may be connected via their nucleobases as in pyrimidine dimers.
The polynucleotide may be double stranded. The polynucleotide is preferably single stranded. The polynucleotide may be one strand from a double stranded polynucleotide.
The polynucleotide can be nucleic acids, such as deoxyribonucleic acid (DNA) or a ribonucleic acid (RNA). The polynucleotide may be any synthetic nucleic acid known in the art, such as peptide nucleic acid (PNA), glycerol nucleic acid (GNA), threose nucleic acid (TNA), locked nucleic acid (LNA), morpholino nucleic acid or other synthetic polymers with nucleotide side chains. The polynucleotide may comprise any of the nucleotides discussed above, including the modified nucleotides. The polynucleotide of the invention is preferably RNA.
SEQ ID NOs: 1 and 2 encode an isolated strain of DWV as discussed above. The VDV-1 CP are encoded by SEQ ID NOs: 3 and 4 within SEQ ID NO: 1 and 2 as discussed above. SEQ ID NOs: 52 to 55 are PCR amplified cDNAs that have been sequenced by standard Sanger sequencing. SEQ ID NOs: 56 to 62 are consensus sequences determined from next generation sequencing, either of naturally infected bees (SEQ ID NOS: 56 to 59) or honeybee pupae injected with virus in the laboratory (SEQ ID NOS: 60 to 62).
The polynucleotide of the invention may comprise a variant sequence based on SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62. A variant sequence is a polynucleotide that has a nucleotide sequence which varies from that of SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 and which typically retains its ability to encode an isolated strain of the invention. The variant sequence preferably retains the ability to encode an isolated strain which comprises the varroa destructor virus 1 (VDV-1) capsid proteins (CP) and the DWV non-structural proteins (NS). The ability of a variant sequence to encode such a virus can be assayed using any method known in the art. For instance, the polynucleotide may be expressed as discussed below.
The variant sequence may comprise any of the nucleotides discussed above, including the modified nucleotides. The variant sequence is typically the same length as SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62, but may be longer or shorter.
Over the entire length of the polynucleotide sequence of SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62, a variant sequence is at least 98% homologous to that sequence based on nucleotide identity. The variant sequence may be at least 98% homologous to the region of SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 that comprises the structural protein module and/or the non-structural protein module of the DWV strain of the invention. Typically, the variant sequence may be at least 98% homologous to the region of SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 that comprises the structural protein module and the non-structural protein module of the DWV strain of the invention.
More preferably, the variant sequence may be at least 98.1%, at least 98.2%, at least 98.3%, at least 98.4%, at least 98.5%, at least 98.6%, at least 98.7%, at least 98.8%, at least 98.9%, at least 98.95%, at least 98.96% at least 98.97% at least 98.98% at least 98.99%, at least 99%, at least 99.1%, at least 99.2%, at least 99.3%, at least 99.4%, at least 99.5%, at least 99.6%, at least 99.7%, at least 99.8%, at least 99.9%, at least 99.95%, at least 99.96%, at least 99.97%, at least 99.98% or at least 99.99% homologous based on nucleotide identity to SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 over the region of SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 that comprises the structural protein module and the non-structural protein module of the DWV strain of the invention. The variant sequence may be at least 98.1%, at least 98.2%, at least 98.3%, at least 98.4%, at least 98.5%, at least 98.6%, at least 98.7%, at least 98.8%, at least 98.9%, at least 98.95%, at least 98.96% at least 98.97% at least 98.98% at least 98.99%, at least 99%, at least 99.1%, at least 99.2%, at least 99.3%, at least 99.4%, at least 99.5%, at least 99.6%, at least 99.7%, at least 99.8%, at least 99.9%, at least 99.95%, at least 99.96%, at least 99.97%, at least 99.98% or at least 99.99% homologous based on nucleotide identity to SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 over its entire sequence (i.e. over its entire length). There may be at least 99%, for example at least 99.5%, at least 99.8%, at least 99.9% or 100%, nucleotide identity over a stretch of 1000 or more, for example 2000, 3000, 4000, 5000, 6000, 7000, 8000 or 9000 or more, contiguous nucleotides of SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 (“hard homology”).
Methods of measuring polynucleotide homology or identity are known in the art. For example the UWGCG Package provides the BESTFIT program which can be used to calculate homology (e.g. used on its default settings) (Devereux et al (1984) Nucleic Acids Research 12, p 387-395).
The PILEUP and BLAST algorithms can also be used to calculate homology or line up sequences (typically on their default settings), for example as described in Altschul S. F. (1993) J Mol Evol 36:290-300; Altschul, S, F et al (1990) J Mol Biol 215:403-10.
Software for performing BLAST analysis is publicly available through the National Centre for Biotechnology Information (http://www.ncbi.nlm.nih.gov/). This algorithm involves first identifying high scoring sequence pair (HSPs) by identifying short words of length W in the query sequence that either match or satisfy some positive-valued threshold score T when aligned with a word of the same length in a database sequence. T is referred to as the neighbourhood word score threshold (Altschul et al, supra). These initial neighbourhood word hits act as seeds for initiating searches to find HSPs containing them. The word hits are extended in both directions along each sequence for as far as the cumulative alignment score can be increased. Extensions for the word hits in each direction are halted when: the cumulative alignment score goes to zero or below, due to the accumulation of one or more negative-scoring residue alignments; or the end of either sequence is reached. The BLAST algorithm parameters W, T and X determine the sensitivity and speed of the alignment. The BLAST program uses as defaults a word length (W) of 11, the BLOSUM62 scoring matrix (see Henikoff and Henikoff (1992) Proc. Natl. Acad. Sci. USA 89:10915-10919) alignments (B) of 50, expectation (E) of 10, M=5, N=4, and a comparison of both strands.
The BLAST algorithm performs a statistical analysis of the similarity between two sequences; see e.g., Karlin and Altschul (1993) Proc. Natl. Acad. Sci. USA 90:5873-5787. One measure of similarity provided by the BLAST algorithm is the smallest sum probability (P(N)), which provides an indication of the probability by which a match between two nucleotide or amino acid sequences would occur by chance. For example, a sequence is considered similar to another sequence if the smallest sum probability in comparison of the first sequence to the second sequence is less than about 1, preferably less than about 0.1, more preferably less than about 0.01, and most preferably less than about 0.001.
One of skill in the art is capable of generating sequences which are complementary to SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 or a variant sequence using canonical base pairing. The complementary sequence may comprise any of the nucleotides discussed above, including the modified nucleotides.
The invention also provides an oligonucleotide which specifically hybridises to a part of a polynucleotide of the invention, i.e. an oligonucleotide which specifically hybridises to part of (a) SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62, (b) a variant sequence as defined above or (c) a sequence which is complementary to (a) or b). Oligonucleotides are short nucleotide polymers which typically have 50 or fewer nucleotides, such 40 or fewer, 30 or fewer, 22 or fewer, 21 or fewer, 20 or fewer, 10 or fewer or 5 or fewer nucleotides. The oligonucleotide of the invention is preferably 21 or 22 nucleotides in length. The oligonucleotide may comprise any of the nucleotides discussed above, including the modified nucleotides. The oligonucleotide may be double stranded. The oligonucleotide is preferably single stranded. The oligonucleotide is preferably RNA.
An oligonucleotide of the invention specifically hybridises to a part of a polynucleotide of the invention, hereafter called the target sequence. The length of the target sequence typically corresponds to the length of the oligonucleotide. For instance, a 21 or 22 nucleotide oligonucleotide typically specifically hybridises to a 21 or 22 nucleotide target sequence. The target sequence may therefore be any of the lengths discussed above with reference to the length of the oligonucleotide. The target sequence is typically consecutive nucleotides within the polynucleotide of the invention.
The target sequence is preferably present within the region of SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 that codes for structural proteins or non-structural proteins. In particular, the target sequence is preferably present in the region of SEQ ID NO: 1 or 2 that encodes the VDV-1 CP or the DWV NS, i.e. within SEQ ID NO: 3 or SEQ ID NO: 4. The target sequence is preferably in a corresponding region in SEQ ID NO: 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62. These can be determined using standard sequence alignment.
An oligonucleotide “specifically hybridises” to a target sequence when it hybridises with preferential or high affinity to the target sequence but does not substantially hybridise, does not hybridise or hybridises with only low affinity to other sequences, for example other sequences in SEQ ID NO: 1 or sequences from other strains of DWV (i.e. strains outside the scope of the present invention). For instance, an oligonucleotide “specifically hybridises” to the VDV-1 CP coding sequence when it hybridises with preferential or high affinity to that sequence (SEQ ID NO: 3) but does not substantially hybridise, does not hybridise or hybridises with only low affinity to other sequences in SEQ ID NO: 1, such as the DWV NS coding sequence in SEQ ID NO: 1 (SEQ ID NO: 4).
An oligonucleotide “specifically hybridises” if it hybridises to the target sequence with a melting temperature (Tm) that is at least 2° C., such as at least 3° C., at least 4° C., at least 5° C., at least 6° C., at least 7° C., at least 8° C., at least 9° C. or at least 10° C., greater than its Tm for other sequences in the HCV genome. More preferably, the oligonucleotide hybridises to the target sequence with a Tm that is at least 2° C., such as at least 3° C., at least 4° C., at least 5° C., at least 6° C., at least 7° C., at least 8° C., at least 9° C., at least 10° C., at least 20° C., at least 30° C. or at least 40° C., greater than its Tm for other nucleic acids. Preferably, the portion hybridises to the target sequence with a Tm that is at least 2° C., such as at least 3° C., at least 4° C., at least 5° C., at least 6° C., at least 7° C., at least 8° C., at least 9° C., at least 10° C., at least 20° C., at least 30° C. or at least 40° C., greater than its Tm for a sequence which differs from the target sequence by one or more nucleotides, such as by 1, 2, 3, 4 or 5 or more nucleotides. The portion typically hybridises to the target sequence with a Tm of at least 90° C., such as at least 92° C. or at least 95° C. Tm can be measured experimentally using known techniques, including the use of DNA microarrays, or can be calculated using publicly available Tm calculators, such as those available over the internet.
More preferably, the oligonucleotide does not hybridise to sequences from other DWV strains or sequences from picorna-like viruses to which bees are susceptible, even under high stringency conditions. Most preferably, the oligonucleotide does not hybridise to any other nucleic acid even under high stringency conditions or sequences from picorna-like viruses to which bees are susceptible.
Conditions that permit the hybridisation are well-known in the art (for example, Sambrook et al., 2001, Molecular Cloning: a laboratory manual, 3rd edition, Cold Spring Harbour Laboratory Press; and Current Protocols in Molecular Biology, Chapter 2, Ausubel et al., Eds., Greene Publishing and Wiley-Interscience, New York (1995)). Hybridisation can be carried out under low stringency conditions, for example in the presence of a buffered solution of 30 to 35% formamide, 1 M NaCl and 1% SDS (sodium dodecyl sulfate) at 37° C. followed by a 20 wash in from 1× (0.1650 M Na+) to 2× (0.33 M Na+) SSC (standard sodium citrate) at 50° C. Hybridisation can be carried out under moderate stringency conditions, for example in the presence of a buffer solution of 40 to 45% formamide, 1 M NaCl, and 1% SDS at 37° C., followed by a wash in from 0.5× (0.0825 M Na+) to 1× (0.1650 M Na+) SSC at 55° C. Hybridisation can be carried out under high stringency conditions, for example in the presence of a buffered solution of 50% formamide, 1 M NaCl, 1% SDS at 37° C., followed by a wash in 0.1× (0.0165 M Na+) SSC at 60° C.
The oligonucleotide of the invention may comprise a sequence which is substantially complementary to the target sequence. Typically, the oligonucleotides are 100% complementary. However, lower levels of complementarity may also be acceptable, such as 95%, 90%, 85% and even 80%. Complementarity below 100% is acceptable as long as the oligonucleotides specifically hybridise to the target sequence. An oligonucleotide may therefore have 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 or more mismatches across a region of 5, 10, 15, 20, 21, 22, 30, 40 or 50 nucleotides. Preferably 100% complementarity is present at positions in the target sequence that are unique to the isolated strain of the invention.
The oligonucleotide of the invention preferably comprises a sequence derived from SEQ ID NO: 5 or 6.
The invention also provides an oligonucleotide which comprises 50 or fewer consecutive nucleotides from a polynucleotide of the invention, i.e. an oligonucleotide which comprises 50 or fewer consecutive nucleotides from (a) SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62, (b) a variant sequence as defined above or (c) a sequence which is complementary to (a) or (b). The oligonucleotide may be any of the lengths discussed above. It is preferably 21 or 22 nucleotides in length. The oligonucleotide may comprise any of the nucleotides discussed above, including the modified nucleotides. The oligonucleotide preferably comprises or consists of any one of SEQ ID NOs: 11 to 51 to 63 to 72. The oligonucleotide may be double stranded. The oligonucleotide is preferably single stranded. The oligonucleotide is preferably RNA.
The oligonucleotide may derived from any region of a polynucleotide of the invention. In particular, the oligonucleotide of the invention may be any of the target sequences discussed above. The oligonucleotide is preferably derived from within the region of SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 that codes for structural proteins or non-structural proteins. The oligonucleotide of the invention is more preferably derived from the region of SEQ ID NO: 1 that encodes that encodes the VDV-1 CP or the DWV NS, i.e. from within SEQ ID NO: 3 or SEQ ID NO: 4. The oligonucleotide of the invention is more preferably derived from the region of SEQ ID NO: 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 that encodes that encodes the VDV-1 CP or the DWV NS.
The invention also provides a polynucleotide which comprise two or more oligonucleotides of the invention. The polynucleotide may be used to generate multiple oligonucleotides of the invention within a bee cell. Typically the polynucleotide is processed by the bee cell to produce two or more oligonucleotides of the invention. The polynucleotide may be processed by the bee cell to produce two or more oligonucleotides of the invention without maintenance of the polynucleotide within the bee cell. These resultant oligonucleotides may then generate a small interfering RNA (siRNA or RNAi) effect as discussed below.
Alternatively, the polynucleotide of the invention may be maintained within a bee cell. For example, the polynucleotide of the invention may be reversed-transcribed within the bee cell and the resulting DNA molecule either imported into the nucleus and integrated into the bee genome or maintained as extra-chromosomal DNA. The DNA can then be transcribed by the bee cell to produce two or more oligonucleotides of the invention, which may be useful in generating a RNAi effect as discussed below.
In a preferred embodiment, the polynucleotide of the invention comprises at least 2, at least 3, at least 4, at least 5, at least 10 or more oligonucleotides of the invention. The polynucleotides of the invention may be used to reduce the likelihood of resistance occurring within the isolated DWV strain as a result of the selection of single point mutations. The use of a polynucleotide comprising more than one oligonucleotide of the invention may give rise to an increased RNAi immune response compared with the use of a single oligonucleotide. This is discussed in more detail below.
In a preferred embodiment, the polynucleotide of the invention may comprise the sequence of SEQ ID NO: 5 and/or SEQ ID NO: 6.
As will be clear from the discussion below, the polynucleotides and oligonucleotides of the invention may be useful in generating an RNAi immune response against the isolated strain of the invention. The polynucleotide or oligonucleotide of the invention is preferably capable of generating an RNAi immune response to the isolated strain of the invention in a bee cell. This RNAi immune response typically leads to gene silencing within the isolated strain. Gene silencing in bees using oligonucleotides has been demonstrated in the art in relation to other bee viruses (Maori et al., Insect Molecular Biology, 2009; 18(1): 55-60; and Hunter et al., PLoS Pathogens, 2010; 6(12): e1001160). Typically an RNAi immune response involves the generation of dsRNA molecules or hairpin RNA molecules. The mechanism of action of such RNA molecules is well characterised in the art. The polynucleotide or oligonucleotide of the invention may be derived from any region to be silenced, including but not limited to, the coding sequence of the gene to be silenced, the 5′ untranslated region, the 3′ untranslated region, the promoter of the gene to be silenced, or any combination thereof.
More preferably the generated RNAi immune response is specific to the isolated strain of the invention. An RNAi immune response may be obtained when the polynucleotide or oligonucleotide is applied exogenously to a bee or bee cell. Exogenous application of the polynucleotide or oligonucleotide may or may not lead to maintenance of the polynucleotide or oligonucleotide in the bee or bee cell. Specifically, exogenous application of the polynucleotide or oligonucleotide may or may not lead to integration of the polynucleotide or oligonucleotide into the bee genome. A polynucleotide of the invention may be delivered to a bee or bee cell. Preferably the polynucleotide will be processed to generate two or more dsRNA molecules comprising oligonucleotides of the invention. These dsRNA molecules may then generate an RNAi immune response.
Alternatively, upon delivery the polynucleotide or oligonucleotide of the invention may be maintained in a bee or bee cell. The polynucleotide or oligonucleotide of the invention may be maintained as a stable genetic element outside the bee genome, i.e. the polynucleotide or oligonucleotide may be maintained within the bee or bee cell without integration into the bee genome. The polynucleotide or oligonucleotide of the invention may be maintained by being integrated into the bee genome. Preferably, the polynucleotide or oligonucleotide of the invention may be capable of acting as described in Goic et al. (Nature Immunology 2013: 14(4): 396-403). Specifically, in a preferred embodiment, the polynucleotide or oligonucleotide is capable of being reverse-transcribed in a bee cell, resulting in DNA forms embedded in retrotransposon sequences. The virus-retrotransposon DNA chimeras may then produce transcripts that are processed by the RNAi machinery within the bee cell, which in turn inhibits replication of the isolated strain. This is discussed in more detail below.
Thus, in generating an RNAi immune response, the polynucleotide or oligonucleotide of the invention is preferably capable of reducing the level of the isolated strain of the invention in a bee cell, bee or bee colony.
The invention also provides a method of generating a polynucleotide or oligonucleotide of the invention derived from SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62. Typically the method involves the design and/or production of primers specific to the target sequence in SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 and using the specific primers to amplify the target sequence. Specific primers can be designed and generated using standard techniques known in the art. Similarly, methods for amplifying target sequences using such primers are also known in the art. A polynucleotide or oligonucleotide of the invention derived from SEQ ID NO: 1, 2, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61 or 62 may also be generated synthetically using known techniques.
The invention provides an antibody which specifically binds to the isolated strain of the invention.
Screening methods are well known to those of skill in the art which may be used to generate and identify antibodies that are capable of specifically binding to the isolated strain of the invention (e.g. Making monoclonals: A practical beginners' guide to the production and characterization of monoclonal antibodies against bacteria and viruses, Newell, Public Health Laboratory Service (1988), ISBN-13: 978-0901144232). Antibodies of the invention can be tested for specific binding to an isolated strain of the invention by, for example, standard ELISA or Western blotting. An ELISA assay can also be used to screen for hybridomas that show positive reactivity with the strain. The binding specificity of an antibody may also be determined by monitoring binding of the antibody to the virus, for example by flow cytometry.
Antibodies of the invention will specifically bind to antigens and epitopes within the isolated strain of the invention. The antigens and epitopes may be identified and used to prepare additional antibodies of the invention.
An antibody “specifically binds” or “specifically recognises” an isolated strain of the invention when it binds with preferential or high affinity to the strain for which it is specific but does not substantially bind, or binds with low affinity, to other viruses or proteins. The specificity of an antibody of the invention for the isolated strain may be further studied by determining whether or not the antibody binds to other related strains or whether it discriminates between them.
An antibody of the invention binds with preferential or high affinity if it binds with a Kd of 1×10−7 M or less, more preferably 5×10−8 M or less, more preferably 1×10−8 M or less or more preferably 5×10−9 M or less. An antibody binds with low affinity if it binds with a Kd of 1×10−6 M or more, more preferably 1×10−5 M or more, more preferably 1×10−4 M or more, more preferably 1×10−3 M or more, even more preferably 1×10−2 M or more. A variety of protocols for competitive binding or immunoradiometric assays to determine the specific binding capability of antibodies are well known in the art (see for example Maddox et al, J. Exp. Med. 158, 1211-1226, 1993).
The term “antibody” as referred to herein includes whole antibodies and any antigen binding fragment (i.e., “antigen-binding portion”) or single chains thereof. An antibody refers to a glycoprotein comprising at least two heavy (H) chains and two light (L) chains inter-connected by disulfide bonds, or an antigen binding portion thereof. Each heavy chain is comprised of a heavy chain variable region (abbreviated herein as VH) and a heavy chain constant region. Each light chain is comprised of a light chain variable region (abbreviated herein as VL) and a light chain constant region. The variable regions of the heavy and light chains contain a binding domain that interacts with an antigen. The VH and VL regions can be further subdivided into regions of hypervariability, termed complementarity determining regions (CDR), interspersed with regions that are more conserved, termed framework regions (FR). The constant regions of the antibodies may mediate the binding of the immunoglobulin to host tissues or factors, including various cells of the immune system (e.g., effector cells) and the first component (Clq) of the classical complement system.
An antibody of the invention may be a monoclonal antibody or a polyclonal antibody, and will preferably be a monoclonal antibody. An antibody of the invention may be a chimeric antibody, a CDR-grafted antibody, a nanobody, a human or humanised antibody or an antigen binding portion of any thereof. For the production of both monoclonal and polyclonal antibodies, the experimental animal is typically a nonhuman mammal such as a goat, rabbit, rat or mouse but may also be raised in other species such as camelids. Polyclonal antibodies may be produced by routine methods such as immunisation of a suitable animal, with the antigen of interest. Blood may be subsequently removed from the animal and the IgG fraction purified. Monoclonal antibodies (mAbs) of the invention can be produced by a variety of techniques, including conventional monoclonal antibody methodology e.g., the standard somatic cell hybridization technique of Kohler and Milstein. The preferred animal system for preparing hybridomas is the murine system. Hybridoma production in the mouse is a very well-established procedure and can be achieved using techniques well known in the art.
The term “antigen-binding portion” of an antibody refers to one or more fragments of an antibody that retain the ability to specifically bind to an antigen. It has been shown that the antigen-binding function of an antibody can be performed by fragments of a full-length antibody. Examples of binding fragments encompassed within the term “antigen-binding portion” of an antibody include a Fab fragment, a F(ab′)2 fragment, a Fab′ fragment, a Fd fragment, a Fv fragment, a dAb fragment and an isolated complementarity determining region (CDR). Single chain antibodies such as scFv antibodies are also intended to be encompassed within the term “antigen-binding portion” of an antibody. These antibody fragments may be obtained using conventional techniques known to those of skill in the art, and the fragments may be screened for utility in the same manner as intact antibodies.
An antibody of the invention may be prepared, expressed, created or isolated by recombinant means, such as (a) antibodies isolated from an animal (e.g., a mouse) that is transgenic or transchromosomal for the immunoglobulin genes of interest or a hybridoma prepared therefrom, (b) antibodies isolated from a host cell transformed to express the antibody of interest, e.g., from a transfectoma, (c) antibodies isolated from a recombinant, combinatorial antibody library, and (d) antibodies prepared, expressed, created or isolated by any other means that involve splicing of immunoglobulin gene sequences to other DNA sequences.
Once a suitable antibody has been identified and selected, the amino acid sequence of the antibody may be identified by methods known in the art. The genes encoding the antibody can be cloned using degenerate primers. The antibody may be recombinantly produced by routine methods.
The invention also provides a method of determining the presence or absence of an isolated strain of the invention. The method of the invention involves detecting the presence or absence of an isolated strain of the invention. In other words, the method involves determining whether or not the isolated strain of the invention is present in the sample. The method may give a positive result, i.e. where the isolated strain of the invention is present in the sample. The method may alternatively give a negative result, i.e. where the isolated strain of the invention is not present in the sample. The method may be helpful to identify diagnose bee colonies as discussed below.
The method comprises contacting the sample with an oligonucleotide of the invention or an antibody of the invention and (b) detecting specific hybridisation of the oligonucleotide or specific binding of the antibody. The presence of specific hybridization or specific binding indicates the presence of the isolated strain in the sample. The absence of specific hybridization or specific binding indicates the absence of the isolated strain in the sample.
The invention also provides a method of determining the concentration of an isolated strain of the invention in a sample, comprising (a) contacting the sample with an oligonucleotide or antibody of the invention and (b) detecting the level of specific hybridisation of the oligonucleotide or specific binding of the antibody and thereby determining the concentration of the isolated strain.
Determining the presence, absence or concentration of the isolated strain allows the detection of the isolated strain of the invention in any bee or bee tissue as described herein. Also, determining the presence, absence or concentration of the isolated strain allows the detection of sequences derived from the strain of the invention, including a polynucleotide and/or oligonucleotide of the invention, that are integrated within the genome of a bee or bee.
The method of determining the concentration of the isolated strain may be used to quantify the isolate strain of the invention by other measures, for example the level, amount or number of viral particles of the isolated strain of the invention.
The sample may be any suitable sample. The invention is typically carried out on a sample that is known to contain or suspected of containing the isolated strain of the invention. The invention may be carried out on a sample that contains one or more viral strains whose identity is unknown. Alternatively, the invention may be carried out on a sample to confirm the identity of the isolated strain of the invention whose presence in the sample is known or expected.
The sample may be a biological sample. The invention may be carried out in vitro on a sample obtained from or extracted from a bee or a bee colony. The sample may be take from any bee tissue. For example, the sample may be taken from tissues from the abdomen, thorax or head of a bee. In particular, a sample may be from brain tissue or gut tissue of a bee. The sample is preferably a fluid sample. Preferably the sample is haemolymph.
The sample may be obtained from a bee at any stage of its lifecycle, including bee eggs, particularly freshly laid eggs, bee larvae, bee pupae and adult bees of any age. The bee may be of any type, including a worker bee, drone bee or queen bee,
The sample to be tested may also be obtained from or extracted from waste material produced during the development of a bee, for example the pupal case that is shed during development, or faeces.
The sample is typically processed prior to being assayed, for example by centrifugation or by passage through a membrane that filters out unwanted molecules or cells. The sample may be measured immediately upon being taken. The sample may also be typically stored prior to assay, preferably below −70° C.
The level of the isolated strain present in a sample encompasses all appropriate measures for quantifying a virus, and includes, for example, the amount, the concentration, the number of virus particles given in any appropriate units.
The presence or level of the isolated strain may be measured direct or indirect means. Any method that allows for the detecting of the isolated strain and the quantification, or relative quantification of the isolated strain may be used. Appropriate technique are known in the art, including, but not limited to, plaque assays, 50% tissue culture infective dose (TCID50) assays, fluorescent focus assays (FFA), protein assays such as the single radial immunodiffusion assay, flow cytometry, quantitative polymerase chain reaction (qPCR) and antibody-based assays.
The level of the isolated strain in a sample of interest may be compared with the level of the isolated strain in another sample. Alternatively, the level of the isolated strain in a sample of interest may be compared with the level of the other DWV strains in the same or another sample.
Any method may be used to detect and quantify specific hybridisation of the oligonucleotide of the invention. The oligonucleotide is typically tagged (or “labelled”) with a molecular marker of either radioactive or fluorescent molecules. Suitable tags (or labels) include, but are not limited to, any of those discussed above, 32P (a radioactive isotope of phosphorus incorporated into the phosphodiester bond in the probe DNA) an digoxigenin, which is a non-radioactive, antibody-based marker.
Any method may be used to detect and quantify specific binding of the antibody of the invention. Methods of quantitatively measuring the binding of an antibody to an antigen are well known in the art. For example, when an isolated strain of the invention is present in the sample, it may bind or substantially bind with the antibody to form antibody-virus complexes, which may then be detected or quantitatively measured. Detection of such complexes is typically carried out using a secondary antibody which recognises general features in the antibody of the invention. The secondary antibody is typically labelled with a detectable label. This facilitates identification of the autoantibody-antigen complex. Any detectable label may be used. Suitable labels include those discussed above with reference to the polynucleotides of the invention.
For instance, the secondary antibody may be conjugated to an enzyme such as, for example, horseradish peroxidise (HRP), so that detection of an antibody-virus complexes is achieved by addition of an enzyme substrate and subsequent colorimetric, chemiluminescent or fluorescent detection of the enzymatic reaction products, or it may be conjugated to a fluorescent or luminescent signal. Alternatively, the secondary antibody may be labelled with a reporter molecule such as a heavy metal or a radioactive tag. Preferably, the intensity of the signal from the secondary antibody is indicative of the relative amount of the antigen-autoantibody complex in the sample when compared to a positive or negative control, and using different dilutions of the samples.
The binding of antibodies to viruses may be detected by any immunological assay technique, of which many are well known in the art. Examples of suitable techniques include enzyme-linked immunosorbent assays (ELISAs), radioimmunoassays, competition assays, inhibition assays, sandwich assays, fluorescent microscopy, microarrays (such as a protein microarray), fluorescence-activated cell sorting (FACS) or the like.
The polynucleotides and oligonucleotides of the present invention may be provided in the form of a vector which includes a promoter operably linked to the inserted sequence, thus allowing for expression of the polynucleotide or oligonucleotide in a bee cell. Any suitable vector may be used which enables the expression of the polynucleotide or oligonucleotide of the invention.
Accordingly, the present invention thus provides a vector that comprises a polynucleotide or oligonucleotide of the invention. The vector is routinely constructed in the art of molecular biology and may, for example, involve the use of plasmid DNA and appropriate initiators, promoters, enhancers and other elements, such as for example polyadenylation signals which may be necessary, and which are positioned in the correct orientation, in order to allow for expression of the polynucleotide or oligonucleotide of the invention. Other suitable vectors would be apparent to persons skilled in the art. By way of further example in this regard refer to Sambrook et al., 2001, Molecular Cloning: a laboratory manual, 3rd edition, Cold Spring Harbour Laboratory Press; and Current Protocols in Molecular Biology, Ausubel et al., Eds., Greene Publishing and Wiley-Interscience, New York (1995)). Therefore, the vector of the invention may be for example, a plasmid, virus or phage vector provided with an origin of replication, a promoter for the expression of the polynucleotide or oligonucleotide and optionally a regulator of the promoter.
In a preferred embodiment vector is a plasmid vector. A plasmid is an autonomously replicating, extrachromosomal circular or linear polynucleotide. The plasmid may include additional elements, such as an origin of replication, or selector genes. Such elements are known in the art and can be included using standard techniques. Numerous suitable expression plasmids are known in the art. For example, one suitable plasmid is phGFP-S65T (Robinson et al. 2000 Insect Molecular Biology 9(6):625-634), pHCMW-G (Kurita et al. 2004 PANS 101(5): 1263-1267). A plasmid vector of the invention may be modified to include retrotrasposon sequences from bees or related insects, for example from Drosophila. Such retrotransposons are known in the art, with 60-450 such sequences being known in honeybees. Any functional retrotrasnposon sequence may be used, including, but not limited to, the mariner element known in the art.
In one embodiment, the plasmid vector of the invention is a bacterial plasmid, i.e. a DNA plasmid comprising T7 promoters operably linked to the polynucleotide or oligonucleotide of the invention. The T7 promoters allow the expression of dsRNA molecules from the polynucleotide or oligonucleotide of the invention, which are useful in generating an RNAi immune response as described herein.
In another preferred embodiment, the vector of the invention is a recombinant viral vector. Suitable recombinant viral vectors include but are not limited to adenovirus vectors, adeno-associated viral (AAV) vectors, herpes-virus vectors, a retroviral vector, lentiviral vectors, baculoviral vectors, pox viral vectors or parvovirus vectors. In the case of viral vectors, administration of the polynucleotide or oligonucleotide is mediated by viral infection of a bee cell.
For example, the polynucleotide or oligonucleotide can be inserted into a viral vector and packaged as retroviral particles using techniques known in the art. The isolated strain of the invention may be produced synthetically be inserted the polynucleotide of the invention into a suitable viral vector or viral envelope. Production of the isolated strain can be carried out be introducing the virus particle into an appropriate host cell in which the virus can replicate. This can be carried out using standard techniques known in the art. The virus particles can then be isolated and/or concentrated using known techniques.
The recombinant virus comprising the polynucleotide or oligonucleotide of the invention can then be isolated and delivered to a bee cell. Preferably the virus is delivered to a bee cell, most preferably a bee germ cell. The virus may be delivered either in vivo or ex vivo, preferably in vitro.
Viral vectors may be based on any suitable virus. For example, retroviral vectors may be based upon the Moloney murine leukaemia virus (Mo-MLV) or human-immunodeficiency virus (HIV). A number of adenovirus vectors are known. Adenovirus subgroup C serotypes 2 and 5 are commonly used as vectors. The wild type adenovirus genome is approximately 35 kb of which up to 30 kb can be replaced with foreign DNA. Suitable adenoviral vectors include Ad5 vectors and simian adenovirus vectors. Viral vectors may also be derived from the pox family of viruses, including vaccinia viruses and avian poxvirus such as fowlpox vaccines. Additional types of virus such as adeno-associated virus (AAV) and herpes simplex virus (HSV) may also be used to develop suitable vector systems.
In a preferred embodiment, the polynucleotide, oligonucleotide or vector of the invention is encapsulated in a viral envelope. In a most preferred embodiment, the viral envelope is a vesicular stomatitis virus envelope (pseudotype virus).
Methods for the production of viral vectors are well known in the art. One non-limiting example of such an approach is a viral vector comprising three separate plasmids: (1) an envelope vector comprising a plasmid encoding a viral coat (capsid) or envelope operably linked to a promoter; (2) a plasmid encoding the non-structural viral genes required for viral replication operably linked to a promoter (the so-called “packaging” vector); and (3) a transducing vector comprising a plasmid in which the sequence of interest operably linked to a promoter.
One such viral vector is described in Coleman et al. (Physiol. Genomics (2002) 12: 221-228). Colmen et al. discloses a viral vector comprising three separate plasmids: (1) the pHEF-VSVG plasmid which encodes the Vesicular Stomatitis Virus (VSV) glycoprotein under the control of an EF1 promoter (commercially available from Addgene, http://www.addgene.org/22501/); (2) the pNHP plasmid which encodes lentiviral gag and pol proteins (commercially available from Addgene, http://www.addgene.org/22500/); and (3) a pTYF-mCMV/SYN-EGFP which encodes a green fluorescent protein under the control of a cytomegalovirus promoter (commercially available from Addgene, http://www.addgene.org/199751).
In the present invention, the transducing vector comprises a polynucleotide or oligonucleotide of the invention. The transducing vector may be derived from any suitable vector, such as the pTYF-mCMV/SYN-EGFP of Colmen et al. to comprise a polynucleotide or oligonucleotide of the invention under the control of a suitable promoter, such as the particular promoters described below.
As an alternative to viral vectors, liposomal preparations can alternatively be used to deliver the polynucleotide or oligonucleotide of the invention. Useful liposomal preparations include cationic (positively charged), anionic (negatively charged) and neutral preparations, with cationic liposomes particularly preferred. Cationic liposomes may mediate intracellular delivery of the polynucleotide or oligonucleotide.
As another alternative to viral vector systems, the polynucleotide or oligonucleotide may be encapsulated, adsorbed to, or associated with, particulate carriers. Suitable particulate carriers include those derived from polymethyl methacrylate polymers, as well as PLG microparticles derived from poly(lactides) and poly(lactide-co-glycolides). Other particulate systems and polymers can also be used, for example, polymers such as polylysine, polyarginine, polyornithine, spermine, spermidine, as well as conjugates of these molecules.
The vector of the invention comprises the polynucleotide or oligonucleotide operably linked to a promoter. Operably linked refers to an arrangement of elements wherein the components so described are configured so as to perform their usual function. Thus, a promoter operably linked to the polynucleotide or oligonucleotide sequence is capable of effecting the expression of that sequence when the proper enzymes are present. Promoters for use in the present invention are discussed further below.
In addition to the polynucleotide or oligonucleotide sequence operably linked to a promoter, additional sequences may be present. These additional sequences are discussed further below.
The vector may be used in vitro, for example to transfect or transform a bee cell to produce a bee cell comprising the vector of the invention. Preferably a vector comprising the oligonucleotide of the invention enables incorporation of the oligonucleotide into the genome of the bee cell. Preferably the bee cell is a bee germ cell. The vectors may also be adapted to be used in vivo, for example to allow in vivo expression of the polynucleotide or oligonucleotide. This is discussed in more detail below.
The polynucleotide, oligonucleotide or vector of the present invention may be administered directly in a “naked” form, meaning that they are not carried in any delivery vehicle. Alternatively, as discussed below, the polynucleotide, oligonucleotide or vector may be administered in a composition comprising a delivery vehicle, such as a sugar syrup.
A “promoter” is a nucleotide sequence which initiates and regulates transcription of a polypeptide-encoding polynucleotide. Promoters can include inducible promoters (where expression of a polynucleotide or oligonucleotide sequence operably linked to the promoter is induced by an analyte, cofactor, regulatory protein, etc.), repressible promoters (where expression of a polynucleotide sequence or oligonucleotide operably linked to the promoter is repressed by an analyte, cofactor, regulatory protein, etc.), and constitutive promoters. It is intended that the term “promoter” includes full-length promoter regions and functional (e.g., controls transcription or translation) segments of these regions.
Promoters for use in the invention are selected to be compatible with the bee cell for which expression is designed. In a preferred embodiment, the promoter is the honeybee heat shock protein 70 (hsp70) promoter. In a most preferred embodiment the honeybee hsp70 promoter comprises the sequence of SEQ ID NO: 8. The honeybee hsp70 promoter is located at nucleotide positions 7306390 to 7307483 of the honeybee genome (Genbank Accession No. NC—007070, Version No. NC—007070.3 GI:323388987). The present inventors have demonstrated that the honeybee hsp70 promoter is function using experiments in bumblebees (Bombus mori), by attaching the promoter in front of a luciferase reporter gene and demonstrating expression of the luciferase reporter gene upon injection of bumblebees with the construct.
The promoter need not be contiguous with the polynucleotide or oligonucleotide sequence, so long as it functions to direct the expression thereof. Thus, for example, intervening untranslated yet transcribed sequences can be present between the promoter sequence and the polynucleotide or oligonucleotide sequence and the promoter sequence can still be considered “operably linked” to the polynucleotide or oligonucleotide sequence.
The vector of the invention may comprise one or more additional sequence in addition to the polynucleotide/oligonucleotide and promoter of the invention. These additional sequences may regulate the expression of the polynucleotide or oligonucleotide.
A polynucleotide, oligonucleotide or vector of the invention may comprise an untranslated leader sequence. In general the untranslated leader sequence has a length of from about 10 to about 200 nucleotides, for example from about 15 to 150 nucleotides, preferably 15 to about 130 nucleotides. Leader sequences comprising, for example, 15, 50, 75 or 100 nucleotides may be used. Generally a functional untranslated leader sequence is one which is able to provide a translational start site for expression of a coding sequence in operable linkage with the leader sequence.
Typically, transcription termination and polyadenylation sequences will also be present, located 3′ to the translation stop codon. Preferably, a sequence for optimization of initiation of translation, located 5′ to the coding sequence, is also present. Examples of transcription terminator/polyadenylation signals include those derived from SV40, as described in Sambrook et al., supra. In a preferred embodiment, a vector of the invention comprises a honeybee transcription terminator sequence of SEQ ID NO: 9. This honeybee transcription terminator is located at nucleotide positions 7311340 to 7311727 of the honeybee genome (Genbank Accession No. NC—007070, Version No. NC—007070.3 GI:323388987).
The polynucleotide, oligonucleotide vector may include transcriptional modulator elements, referred to as “enhancers”. Enhancers are broadly defined as a cis-acting agent, which when operably linked to a promoter and polynucleotide/oligonucleotide sequence, will increase transcription of that polynucleotide or oligonucleotide sequence. Enhancers can function from positions that are much further away from the polynucleotide or oligonucleotide sequence of interest than other expression control elements (e.g. promoters), and may operate when positioned in either orientation relative to the sequence of interest. Enhancers have been identified from a number of viral sources, including polyoma virus, BK virus, cytomegalovirus (CMV), adenovirus, simian virus 40 (SV40), Moloney sarcoma virus, bovine papilloma virus and Rous sarcoma virus. Examples of suitable enhancers include the SV40 early gene enhancer, the enhancer/promoter derived from the long terminal repeat (LTR) of the Rous Sarcoma Virus, and elements derived from human or murine CMV, for example, elements included in the CMV intron A sequence.
The vector may contain one or more selectable marker genes, for example an ampicillin resistance gene in the case of a bacterial plasmid or a resistance gene for a fungal vector.
The vector may include a retroviral polynucleotide or oligonucleotide which facilitate the insertion of the oligonucleotide of the invention into the bee genome. Typically plasmid vectors of the invention include such retroviral sequences.
As discussed above, the inventors have identified a single strain of DWV, the isolated strain of the invention, that is predominant in bees and bee colonies infested with Varroa mites and which is present at high levels in the infested bees and colonies.
Therefore, inhibitors of the isolated strain of the invention may be used to treat and/or prevent infection by the strain, and so reduce the impact of infection on the bee or bee colony. Inhibitors of the isolated strain of the invention may also be used to treat and/or prevent deformed wing in a bee or bee colony infested with Varroa mites.
Accordingly, the present invention provides methods and compositions for inhibiting the isolated strain of the invention. The inhibitors for use in the invention are able to reduce the level, amount or concentration of the isolated strain within a bee. The compositions and methods of the invention can be used in individual cells, cells or tissue in culture, or in vivo in bees, or in organs or other portions of bees. The compositions and methods of the invention can be used to treat part of or entire bee colonies.
Inhibition of the isolated strain of the invention may be measured by any suitable means. For example, the Examples herein disclose methods for determining infection or viral load. Infection or viral load may be used as an indicator of the effectiveness of the inhibitors of the invention: the more effective an inhibitor, the greater the decrease in infection or viral load.
According to the present invention, inhibitors of the isolated strain of the invention can be used for reducing the presence and/or level of viral particles in a bee cell, bee or bee colony. The inhibitors may achieve this reduction by reducing the expression of or cause silencing or knockdown of any gene of the isolated strain of the invention. The inhibitors of the invention may cause reduction of expression of at least one viral protein. Typically, the reduction in level of viral particles in the cell, bee or bee colony may be at least 50% of the amount of viral particles in the absence of the inhibitor. Preferably, the reduction is at least 55%, at least 60%, at least 65%, at least 70%, at least 75% or at least 80% or more preferably, at least 85%, at least 90% or at least 95%. A method for determining the relative amount of viral particles may be any suitable method known in the art.
Any suitable inhibitor may be used according to the invention, for example peptides and peptidomimetics, antibodies, small molecule inhibitors, double-stranded and antisense RNA, aptamers and ribozymes. Preferred inhibitors are polynucleotide or oligonucleotide inhibitors, more preferably oligonucleotide inhibitors, and most preferably RNA oligonucleotide inhibitors.
Inhibitors for use in the invention include the antibodies, polynucleotides, oligonucleotides, vectors, and compositions of the invention. The most preferred inhibitors are the polynucleotides and oligonucleotides of the invention.
The inhibitors of the invention may be delivered to a bee or bee cell by any appropriate means. In particular, the polynucleotides and oligonucleotides of the invention may be delivered as naked sequences, preferably in the form of dsRNA or hairpin RNA structures, as DNA molecules operably linked to a suitable promoter, in a vector, particularly a plasmid or viral vector as discussed herein.
In some embodiments, the polynucleotide, oligonucleotide, antibody or vector of the invention is contained in a composition comprising a delivery vehicle. The delivery vehicle is preferably a sugar syrup. This allows the composition to be fed to a target bee and allows uptake of the molecule or inhibitor by the target bee. In a preferred embodiment the sugar syrup comprises a 1:1 w/v mix of granulated sugar and water, a 2:1 w/v mix of granulated sugar and water, invertase treated sugar or sugar solutions. Fondant, particularly bakers fondant may also be used a delivery vehicle. Alternatively, or in addition to these sugars-based delivery vehicles, real or artificial pollen patties may be used.
The formulated compositions will include an appropriate amount of the polynucleotide, oligonucleotide, antibody or vector of the invention which is sufficient achieve the desired result. For example, the composition may comprise an amount of the oligonucleotide of the invention that is sufficient to generate an RNAi immune response against the isolated strain of the invention when the composition is contacted with a bee. An appropriate effective amount can be readily determined by one of skill in the art. Such an amount will fall in a relatively broad range that can be determined through routine trials. The compositions may can be administered directly to a bee or, alternatively, delivered ex vivo, to bee cells, using methods known to those skilled in the art.
The invention concerns treating and diagnosing bees. The bee may be any type of bee, including, without limitation, honeybees (of the genus Apis) and bumblebees (of the genus Bombus). Preferably the bee is Apis mellifera, Apis cerana, Bombus terrestris or Bombus impatiens. More preferably the bee is Apis mellifera. The term “bee” as used herein covers all stages of the bee life cycle. Thus, egg, larval, pupal and adult bees are all covered by the term “bee”.
Bee cells, including bee germ cells, may also be used according to the present invention. A bee germ cell is a cell that is responsible for transferring the bee genome to the next generation of bee. Examples include bee gametes, such as bee spermatozoa or bee eggs (oocytes).
The present invention may also be applied to a bee colony. A bee colony typically comprises at least one queen bee and two or more, preferably many, worker and/or drone bees. Preferably the bee colony also comprises at least one bee egg and/or at least one bee larva. The bee colony may also comprise at least one bee pupa.
The compositions are administered to a bee in an amount that is compatible with the dosage formulation and that will be prophylactically and/or therapeutically effective. An appropriate effective amount will fall in a relatively broad range but can be readily determined by one of skill in the art by routine trials.
As used herein, the term “prophylactically or therapeutically effective dose” typically means a dose in an amount sufficient to reduce the level of an isolated strain of the invention as discussed above. The term “prophylactically or therapeutically effective dose” may mean a dose in an amount sufficient to alleviate, reduce, cure or at least partially arrest symptoms and/or complications from infection with an isolated strain of the invention. The term “prophylactically or therapeutically effective dose” may mean a dose in an amount sufficient to alleviate, reduce, cure or at least partially arrest symptoms and/or complications from Varroa mite-induced deformed wing.
Prophylaxis or therapy can be accomplished by a single direct administration at a single time point or by multiple administrations, optionally at multiple time points. Administration can also be delivered to a single or to multiple sites within a bee. Further, administration may be to an individual bee or to more than one bee. Administration may be to the whole or part of one or more bee colony. Those skilled in the art can adjust the dosage and concentration to suit the particular route of delivery. In one embodiment, a single dose is administered to a bee or a bee colony on a single occasion. In one embodiment, multiple doses are administered to a bee or bee colony on multiple occasions. Any combination of such administration regimes may be used.
Different administrations may be performed on the same occasion, on the same day, one, two, three, four, five or six days apart, one, two, three, four or more weeks apart. Preferably, administrations are 1 to 5 weeks apart, more preferably 2 to 4 weeks apart, such as 2 weeks, 3 weeks or 4 weeks apart. The schedule and timing of such multiple administrations can be optimised for a particular composition or compositions by one of skill in the art by routine trials.
For the administration of inhibitors of the invention, particularly exogenous polynucleotides or oligonucleotides of the invention (for the generation of an RNAi immune response) via a feeding delivery route, the preferred times for delivery are early spring, typically March or April in the UK, or late summer/early autumn, typically August or September in the UK. Bee keepers typically feed bees at both these times, such that either solid or liquid food may be used as a delivery vehicle as described herein. Typically bee keepers feed bees for three to six weeks within these time periods. The feeding periods are weather-dependent and may be altered as necessary according to the specific conditions.
Further, administration may take place at any other time a bee keeper routinely feeds the bee. For example, bee keepers typically feed bees after hiving a swarm. Therefore, an inhibitor of the invention may be administered on feeding after hiving a swarm. This is a particularly preferred time for administration, because it would allow the treatment and/or prevention of the isolated strain of the invention or Varroa mite-induced deformed wing disease according to the invention in a newly formed bee colony.
Dosages for administration will depend upon a number of factors including the nature of the composition, the route of administration and the schedule and timing of the administration regime. The dose will also vary according to the severity of the infection with the isolated strain and/or the deformed wing in a Varroa mite-infested bee or bee colony. The dose will also vary according to number of bees or bee colonies to be treated. Optimum dosages may vary depending on the relative potency of the oligonucleotide or antibody of the invention, and can generally be estimated based on EC50s found to be effective in vitro and in in vivo animal models.
As discussed above, the isolated strain of the invention has been identified at high concentrations and with low diversity in bee colonies infested with Varroa mites. DWV infection and Varroa infestation are coincident with the appearance of pathology in the bees, including the appearance of wing deformities. In contrast, bee colonies not infested with Varroa mites have low-concentration, covert DWV infections with a high diversity of DWV strains.
Accordingly, the present invention provides a method of treating or preventing in a bee or bee colony infection with the isolated strain of the invention, comprising contacting the bee or bee colony with an inhibitor of the isolated strain. The invention also provides an inhibitor of the isolated strain of the invention for use in a method of treating and/or preventing in a bee or bee colony infection with the isolated strain, comprising contacting the bee or bee colony with an inhibitor of the isolated strain. The invention also provides use of an inhibitor of the isolated strain of the invention in the manufacture of a medicament for treating and/or preventing in a bee or bee colony infection with the isolated strain
The invention also provides a method of treating or preventing deformed wing disease in a Varroa mite-infested bee or bee colony infection, comprising contacting the bee or bee colony with an inhibitor of the isolated strain. The invention also provides an inhibitor of the isolated strain of the invention for use in a method of treating and/or preventing deformed wing disease in a Varroa mite-infested bee or bee colony, comprising contacting the bee or bee colony with an inhibitor of the isolated strain of the invention. The invention also provides use of an inhibitor of the isolated strain of the invention in the manufacture of a medicament for treating and/or preventing deformed wing disease in a Varroa mite-infested bee or bee colony.
In a preferred embodiment the inhibitor is specific for the isolated strain of the invention. An inhibitor is specific for the isolated strain of the invention if is inhibits the isolated strain of the invention, but inhibits other viruses, such as other strains of DWV, to a lower degree or not at all. In a more preferred embodiment, the inhibitor is a polynucleotide of the invention, an oligonucleotide of the invention, an antibody of the invention, a vector of the invention or a composition of the invention.
The therapeutic or prophylactic method of the invention may comprise contacting at least one bee with an oligonucleotide of the invention to stimulate the bee's immune system to generate resistance to the isolated strain of the invention.
The methods of the invention preferably involve generating bees that are resistant to the isolated strain of the invention by contacting at least one bee with a polynucleotide or an oligonucleotide of the invention and thereby stimulating the at least one bee's immune system to generate resistance to the isolates strain. In a preferred embodiment, the polynucleotide or oligonucleotide is fed to the at least one bee. Typically the polynucleotide or oligonucleotide is operably linked to a promoter, contained in a vector and encapsulated in a viral envelop (preferably from vesicular stomatitis virus) to form a pseudotyped retrovirus particle. The pseudotyped retrovirus particles are then typically added to a sugar syrup, facilitating their ingestion by the at least one bee. Adult bees and larval bees will ingest the syrup.
Alternatively, the pseudotyped retrovirus particles may be mixed with drone semen and used during instrumental insemination of virgin bee queens to produce transgenic fertilised bee queens that are capable of transmitting the RNAi-expressing DNA vertically to progeny bees.
Alternatively, a non-integrating or integrating plasmid vector of the invention capable of being maintained and vertically transmitted may be introduced to drone semen for use in instrumental insemination.
The at least one bee is contacted with a polynucleotide or an oligonucleotide of the invention, either in its naked form or comprised in a vector or composition of the invention. On contact, the bee will ingest the polynucleotide or oligonucleotide. After ingestion, the polynucleotide or oligonucleotide is taken up by cells of the at least one bee. Once inside a bee cell, the polynucleotide or oligonucleotide is processed by the bee's intracellular machinery to generate an RNAi immune response. In particular, the polynucleotide or oligonucleotide may be expressed in the bee cell in a hairpin configuration. These hairpin sequences do not integrate into the bee genome, but may be processed by the endogenous enzyme (Dicer), and associate with Argonaute protein (Ago). Dicer and Ago are involved in the early stages of RNAi processing and expression. Dicer cleaves the dsRNA into siRNAs. Ago directs the resulting siRNAs protein complex, the RNA-induced silencing complex (RISC), from the bee cell to complementary RNA molecules from the isolated strain within the bee cell. Contact of the RISC to the viral RNAs targets the viral RNAs for degradation using cellular ribonucleases or inhibits translation of the RNA. Thus, the ability of the isolated strain to replicate in the bee cell is reduced. Consequently, the bee itself is resistant to the isolated strain.
Entry of polynucleotides or oligonucleotides delivered in a viral vector of the invention to a bee cell is mediated by the viral coat or capsid, which is preferably the VSV capsid as described herein. Once the bee cell is infected with the viral vector, the viral vector migrates to the bee cell nucleus and integrates. The vector payload, i.e. the polynucleotide or oligonucleotide of the invention operably linked to a suitable promoter is then expressed in a form that generates dsRNA molecules. These dsRNA molecules are processed by Dicer, leading to the Ago and RISC-mediated cascade described above.
Expression of one or more oligonucleotides of the invention as hairpin sequences allows a highly focused means of suppressing DWV replication. Alternatively, multiple oligonucleotides of the invention or a more extensive region of the polynucleotide of the invention may be expressed, enabling the generation of longer dsRNA for Dicer-mediated processing. Such an embodiment may be used to reduce the risk of resistant forms of DWV emerging, because a greater number of RNAs from the isolated strain can be targeted by the bee protein complexes. Such multiple oligonucleotides or more extensive polynucleotides may be expressed with two opposing promoters, preferably honeybee promoters.
The present invention also provides a transgenic bee that is resistant to the isolated strain of the invention. The transgenic bee may be an adult bee or a larval bee. Typically the bee is engineered to express an oligonucleotide of the invention, making the bee resistant to the isolated strain of the invention. Typically the transgenic bee comprises at least one cell that has been modified in a way that confers resistance to the isolated strain of the invention. The transgenic bee may be modified such that at least one cell expresses at least one oligonucleotide of the invention. The transgenic bee may be modified such that an oligonucleotide or polynucleotide of the invention has been incorporated into the genome of at least one cell, such that at least one oligonucleotide of the invention is synthesised within the at least one cell. The invention also provides a transgenic bee that is resistant to Varroa mite-induced deformed wing disease, wherein at least one cell of the bee has been modified in a way that confers resistance to deformed wing disease. The transgenic bee may be modified such that at least one cell expresses at least one oligonucleotide of the invention. The transgenic bee may be modified such that an oligonucleotide or polynucleotide of the invention has been incorporated into the genome of at least one cell, such that at least one oligonucleotide of the invention is synthesised within the at least one cell.
Typically a sufficient number of bee cells have been modified as described above such that the bee is resistant to the isolated strain or to deformed wing disease. The polynucleotide or at least one oligonucleotide may be stably incorporated into the bee genome. Alternatively, the polynucleotide or at least one oligonucleotide, or the vector, particularly viral or plasmid vector, containing the polynucleotide or at least one oligonucleotide may be retained within the bee cell(s) without incorporation into the genome.
In a preferred embodiment, the transgenic bee of the invention is a queen bee. The queen bee may be an adult queen bee or a larval queen bee. A transgenic queen bee allows the production of a bee colony which is resistant to the isolated strain of the invention.
The invention also provides one or more bee germ cells which expresses at least one oligonucleotide of the invention. Preferably the oligonucleotide is incorporated into the genome of the bee germ cell(s). The one or more bee germ cell may comprise a polynucleotide of the invention incorporated into the genome of at the bee germ cell(s), such that at least one oligonucleotide of the invention is synthesised within the at least one cell. Thus, the invention provides one or more bee germ cells comprising an oligonucleotide, poly nucleotide or vector of the invention. The one or more bee germ cells are preferably one or more bee spermatozoa.
The transgenic bee of the invention may be designed to express regions of the genome of the isolated strain of the invention, such a way that they induce small interfering RNA (siRNA or RNAi) molecules. In a preferred embodiment, these oligonucleotides are engineered into the genome of the bee, such that the resulting resistance to the isolate strain is stable and may be passed on to subsequent generations.
Bees can be considered a superorganism with a single fertilized queen, several thousand female (unfertilized) workers and a small number of male (fertile) drones within a single colony. Queens mate within 2 weeks of emergence and never leave the colony (other than when swarming). Workers cannot mate and only under exceptional circumstances lay unfertilized eggs (which develop into drones). Fertilized eggs develop into queens or workers depending upon how they are maintained in the first 3 days after hatching. Beekeepers have developed ways to raise queens from any newly hatched larvae.
To introduce polynucleotides or oligonucleotides of the invention into the germline of bees it is necessary to manipulate a queen bee to lay fertilized eggs containing the polynucleotides or oligonucleotides. Sperm mediated gene transfer (SMGT) in honeybees is known in the art (Robinson et al. 2000 Insect Molecular Biology 9(6): 625-634).
Accordingly, the invention provides a method of generating a transgenic queen bee that is resistant to infection by an isolated strain of the invention, comprising (a) incorporating a polynucleotide or an oligonucleotide of the invention into the genome of one or more bee germ cells, such as one or more spermatozoa; and (b) generating the queen bee from said one or more germ cells. This method may also be used to generate a transgenic queen bee that is resistant to Varroa mite-induced deformed wing disease. Bee spermatozoa may be collected from one or more drone bee.
The polynucleotide or oligonucleotide of the invention may be contacted with the bee germ cell using methods known in the art (Robinson et al 2000), such that the polynucleotide or oligonucleotide is incorporated into the genome of the bee germ cell. The polynucleotide or oligonucleotide may be contacted with the bee germ cell as a naked polynucleotide or oligonucleotide, or in a vector or composition of the invention. In particular, the polynucleotide or oligonucleotide may be incorporated into the bee germ cell using retroviral mediated infection. The retroviral vector used in said method may be as described herein.
The bee germ cell, particularly the bee spermatozoa, comprising the polynucleotide or oligonucleotide of the invention may then be transferred to a virgin queen bee. In the case of bee spermatozoa, the spermatozoa will migrate to the spermatheca where they are stored for months to years and can be used to fertilise the eggs of the virgin queen bee to generate a transgenic bee. Preferably the transgenic bee is a queen bee.
Once the polynucleotide or oligonucleotide of the invention has been contacted with the bee germ cell it will first be taken up by that cell. After uptake of the polynucleotide or oligonucleotide by the bee germ cell, the polynucleotide or oligonucleotide will be incorporated into the genome of the bee germ cell. This can be mediated by retroviral sequences provided in a vector of the invention as discussed above. Incorporation of the polynucleotide or oligonucleotide into the genome allows the oligonucleotide to be transmitted vertically to a progeny bee, thus enabling the production of bee that are resistant to the isolated strain of the invention or to Varroa-mite induced deformed wing.
As discussed above, the polynucleotides and oligonucleotides of the invention may be useful in generating an small interfering RNA (siRNA or RNAi) immune response against the isolated strain of the invention. As a result of this RNAi immune response, the transgenic bee will be resistant to the isolated strain of the invention or to Varroa-mite induced deformed wing disease. In a preferred embodiment, the polynucleotide or oligonucleotide is reverse-transcribed in a bee cell, resulting in DNA forms embedded in retrotransposon sequences to generate the RNAi immune response.
Alternatively, transgenic bees of the invention may be generated by other suitable techniques, such as feeding developing drones with the retroviral vectors of the invention, or microinjecting developing drone larvae or pupae with the retroviral vectors of the invention.
Once a transgenic queen bee of the invention has been produced, this transgenic queen can be used produce transgenic progeny bees that also comprise the oligonucleotide of the invention and so are resistant to the isolated strain of the invention or to Varroa-mite induced deformed wing disease. The transgenic queen bee of the invention can be used to initiate new a new bee colony, the worker and drone bees of which will all transgenic and hence resistant to the isolated strain of the invention or to Varroa-mite induced deformed wing disease. Accordingly, the invention also provides a method of generating a transgenic bee of the invention, comprising producing the transgenic bee from a transgenic queen bee of the invention or from a transgenic queen bee producing using a method of the invention.
The invention further provides a method of preventing in a bee or bee colony an infection with an isolated strain of the invention, comprising using a transgenic queen bee of the invention or from a transgenic queen bee produced using a method of the invention to generate the bee or the bee colony.
The invention further provides a method of preventing deformed wing disease in a Varroa mite-infested bee or bee colony, comprising using a transgenic queen bee of the invention or from a transgenic queen bee produced using a method of the invention to generate the bee or bee colony.
The invention further provides a method of producing a bee or bee colony that is resistant to infection from an isolated strain the invention, comprising using a transgenic queen bee of the invention or from a transgenic queen bee produced using a method of the invention to generate the bee or bee colony.
As discussed herein, the inventors have identified a sub-population of bees present in Varroa mite-infested colonies that have sequences from the strain of the present invention integrated into their genomes. Such bees are naturally resistant to the isolated strain. However, whilst a sub-population of bees in Varroa mite-infested colonies might be naturally resistant to the isolated strain, the presence of such a naturally resistant sub-population is not sufficient for colony survival. In practice, a bee colony will only survive if a sufficient number of individual bees are resistant to the isolated strain.
Bees comprising such incorporated sequences may be identified according to the present invention, for example using the oligonucleotides or antibodies of the invention in the methods of determining the presence or absence of the isolated strain of the invention. These naturally resistant bees may be used to generate bees or bee colonies that are also resistant to the isolated strain of the invention, or that are resistant to Varroa mite-induced deformed wing disease.
The invention provides for the first time a bee colony comprising sufficient bees that are individually resistant to infection by the isolated strain of the invention, such that the colony as a whole is resistant infection by the isolated strain or to Varroa mite-induced deformed wing disease. Typically at least 50%, at least 60%, at least 70%, at least 80, at least 90% or more, up to 100% of the individual bees within a colony must be resistant to the isolated strain in order for the colony to be resistant to the isolated strain or to Varroa mite-induced deformed wing disease.
Accordingly, the invention provides a bee colony wherein at least 50% of the bees within the colony are resistant to the isolated strain of the invention. The invention also provides a bee colony wherein at least 50% of the bees within the colony are resistant to Varroa mite-induced deformed wing disease. Typically the at least 50% of the bees comprise an oligonucleotide of the invention, a polynucleotide of the invention comprising two or more oligonucleotides of the invention, or a vector of the invention. In a preferred embodiment, the individual bees that are resistant to the isolated strain of the invention have an oligonucleotide of the invention or a polynucleotide of the invention comprising two or more oligonucleotides of the invention incorporated into their genome.
The present invention also provides a method of producing a bee or bee colony that is resistant to the isolated strain of the invention comprising: identifying a bee or bee larva that have sequences from the strain of the invention incorporated into the bee or bee larva genome; generating a queen bee from said bee or bee larva; and using said queen bee to generate the bee or bee colony. The present invention also provides a method of producing a bee or bee colony that is resistant to Varroa mite-induced deformed wing disease comprising: identifying a bee or bee larva that have sequences from the strain of the invention incorporated into the bee or bee larva genome; generating a queen bee from said bee or bee larva; and using said queen bee to generate the bee or bee colony.
The invention provides a method of diagnosing in a bee or bee colony infection with an isolated strain of the invention, comprising determining the presence or absence of the isolated strain, wherein the presence of the isolated strain is indicative of the presence of infection with the isolated strain and wherein the absence of the isolated strain is indicative of the absence of infection with the isolated strain. Typically the determination step takes place in vitro in sample from the bee or colony. Suitable samples are discussed above.
The invention also provides a method of diagnosing deformed wing disease in a Varroa mite-infested bee or bee colony, comprising determining the presence or absence of an isolated strain of the invention, wherein the presence of the isolated strain is indicative of the presence of deformed wing disease and wherein the absence of the isolated strain is indicative of the absence of deformed wing disease. Typically the determination step takes place in vitro in sample from the bee or colony. Suitable samples are discussed above.
In a preferred embodiment, an oligonucleotide of the invention or an antibody of the invention is used to identify the presence or absence of the isolated strain as described above.
Methods of determining the presence or absence of the isolated strain of the invention and methods of determining the level of the isolated strain as disclosed herein may be used in the diagnosis in a bee or bee colony of infection with the isolated strain of the invention or the diagnosis of deformed wing disease in a Varroa mite-infested bee or bee colony.
The invention is illustrated by the following Examples:
There are indications that in the absence of the Varroa mite DWV is transmitted vertically through eggs and male sperm. This route of transmission favours covert infection with low pathogenicity and DWV levels, because it depends on the survival of the host. Moreover, it cannot be ruled out that DWV may be associated with favourable bee traits. For example KV (a Japanese strain of DWV) is associated with aggressive behaviour traits, which could be beneficial for the colony survival. With the spread of Varroa destructor mites, which reached the UK in the 1990s, overt DWV infections became common. There is a possibility that pathogenicity of DWV increased in Varroa-infested colonies due to changes in DWV populations facilitated by the mite. The Varroa mite, which is a host to at least some DWV variants, provides a route for effective horizontal transmission of DWV, therefore the virus was no longer dependent on the survival of infected individuals. Moreover, aggressive strains of DWV accumulating to higher levels faster could be more easily acquired from infected bees and transmitted by the mite. Changes in DWV population composition on a large temporal and spatial scale following Varroa infestation were reported for Hawaiian honeybees, although this study did not provide enough DWV sequence data to pinpoint the identity of the winning virus strains. There is no agreement on whether the feeding by Varroa mites suppresses bee innate immunity, which would make possible the replication of DWV to high levels. Earlier reports showed that bees exposed to Varroa feeding are immuno-compromised, but later transcriptome analysis studies found little or no evidence of negative effects of Varroa on the expression of innate immunity genes. It should be noted, however, that analysis of the impact of Varroa mite feeding on antivirus defences in bees is difficult, because these defences remain poorly understood in all insects, including the much better studied Drosophila and mosquito. Although RNA interference (RNAi) is considered to be the major defence against RNA viruses in insects, there is evidence suggesting that Toll, Imd and Jak-Stat signalling pathways are also involved in antivirus resistance. It is possible that Varroa destructor interferes with antivirus pathways in bees and that contributes to activation of DWV.
The present inventors sought to explore both hypotheses on the role of Varroa in dramatic increase of DWV pathogenicity and the replication levels in bees, namely: (i) Varroa infestation results in the introduction of pathogenic strains of DWV, which were selected for their ability to be effectively transmitted by Varroa mites, or (ii) Varroa mites disrupt bee antivirus defence mechanisms, including immune defence pathways, RNAi machinery, or tissue barriers. To test these distinct, but not mutually exclusive hypotheses, the inventors devised a system which modelled encounter of naïve Varroa-free bees with Varroa mites and the mite-associated DWV strains. We infected bees from a Varroa-free honeybee colony with the DWV strains circulating in a Varroa-infested colony, either solely orally or orally and via injection into bee haemolymph by the mites. The inventors analysed the virus population dynamics, RNAi responses, and the global gene expression in the bees subjected to the above treatments and in the control bees. Analysis of the resulting changes in DWV populations and in the bee transcriptome provided data indicating that both introduction of pathogenic strains of DWV and suppression of antivirus defences play their parts in activation of DWV replication in Varroa-infested bees. The results suggest that pathogenic strains of DWV are effectively suppressed when delivered orally, although oral infection results in changes in DWV diversity. Varroa mite feeding on developing pupae, which allowed direct injection of pathogenic DWV strains into haemolymph as well as suppression of innate immunity, were the likely causes of DWV activation in the mite associated pupae.
Frame Transfer System for Honeybee Infestation with DWV and Varroa Mites
Infestation of bees included transfer of a brood frame with newly hatched larvae from Varroa-free colony to Varroa-infested colony and subsequent sampling the capped pupae. A Varroa-free colony of black British honeybee, Apis mellifera mellifera with a naturally-mated one year old queen was imported from Colonsay, Scotland, an island with no historic reports of Varroa incidence and no imports of bees from Varroa-infested areas. The presence of DWV strains associated with Varroa mite infestation could thus be excluded. As a source of Varroa mites and the mite-associated DWV strains, a Warwickshire honeybee colony, heavily infested with Varroa and having high DWV levels in bees and mites was selected. The Varroa-free and Varroa-infested colonies were contained in separate mesh flight cages and were maintained on an artificial diet of sugar syrup and pollen. The pollen was imported from Varroa-free Australia to exclude contamination with Varroa-associated viruses through foraged food. In order to minimise possible effects on bee gene expression due through differences in nutrition, both the control Varroa-free and the Varroa-infested colonies were maintained in flight cages in the same apiary (at the University of Warwick, UK) and were fed on the artificial diet for two months before the start of the frame transfer experiments.
The experimental infestation, summarised in
In this study we carried out comprehensive analysis of virus diversity and gene expression in the same individual bees which were subjected to oral infection with the strains of DWV associated with Varroa destructor infestation, or both oral DWV in, as well as analysis siRNA and miRNA.
The whole individual honeybee pupae were ground to fine powder in liquid nitrogen, and a half of ground bee was used for RNA extraction was carried out using 1 mL of Trizol Reagent (Invitrogen) according to the manufacturer's instructions. Total RNA extraction from mites was carried out using RNeasy spin columns (Qiagen RNeasy Plant Mini kit). Real-time reverse transcription PCR was carried out essentially as in Moore et al 2011, in brief, RNA extracts were treated with DNAse, then purified DNA-free total RNA preparations were used to as a template to produce cDNA using random primer and Superscript III reverse transcriptase (Invitrogen). The cDNA samples produced were used for real-time PCR quantification of the DWV or host transcripts using SYBR green mix (Agilent Technologies).
Individual bee samples were used to take into account differences in pathogen load and diversity between individuals, as well as to have an opportunity to consider effect of genetic differences on transcriptional responses to the treatments. Bees were analysed bees at the same developmental stage, narrowing to one day, to reduce developmental differences in gene expression which would allow to detect changes induced by parasite and/or pathogens.
Amplification of the cDNA fragments corresponding to the central region of DWV genomic RNA was carried out by nested PCR using GoTaq PCR mix (Promega) and appropriate primers in the cDNA to total RNA samples produced from bee or Varroa mite samples in the pooled bee sample from the Varroa-infested colony used in the transfer experiment.
The outside PCR primers were designed to amplify known DWV strains. For each first round reaction four second round amplification reactions were carried out using VDV-1 or DWV specific primers which allowed VDV-1-type and DWV-type CP and NS to be distinguished, thus amplifying all potential combinations even those which were present at very low levels. The PCR fragments were cloned into plasmid vector, pGemT-Easy (Promega) and sequenced using Sanger dideoxy method.
Worker honeybees from a Varroa-free colony (sourced in a region with no historic contacts with Varroa), were infected orally with the Varroa-associated DWV strains at the larval stage (from day four until capping on day nine). A proportion of the pupae capped with Varroa mites were also subjected to the mites feeding on their haemolymph during the pupal stage until sampling on day fifteen, five days after capping. Feeding of the mites (adult females) on pupae was confirmed by the presence of at least one protonymph. Sampling capped bee pupae allowed not only identification of the individual bees on which Varroa mites were fed, but also provided bees at the developmental stage when DWV and Varroa mite introduction results in highly pronounced symptoms, such as deformed wings.
The total levels of DWV viruses in each collected pupae was assessed by qRT-PCR using a pair of primers for the conserved polymerase-coding region, designed to detect all known DWV strains, including DWV, VDV-1 and KV. The real-time PCR Ct values showed bimodal distribution, with distinct low DWV level and high DWV level groups (P<0.0001) (
A sample of eight randomly selected bee pupae was taken from each of the four treatment groups for further in-depth analysis of DWV diversity (which included real time qPCR quantification of DWV strains, sequencing of the cDNA clones and next generation sequencing of genomic RNA), transcriptional responses (which included microarray transcriptional profiling and high-throughput micro RNA sequencing), and antiviral RNAi responses (which included high throughput siRNA sequencing).
The levels of DWV, VDV-1 and their recombinants were quantified in individual pupae and the mites associated with the pupae. Sets of primers designed to distinguish between the DWV-type and VDV-1-type sequences coding for the capsid proteins (CP) and the non-structural proteins (NS) were used (Table 1 below).
The level of the DWV-type and VDV-1-type CP and NS was determined in each of the 32 bees (eight bees in each of the four treatment groups) and in the 15 Varroa mite samples associated with the VH and VL pupae. The 32 bees were also used for the whole genome expression microarray analysis.
The strain-specific quantification showed that the bees in each of the four treatment groups had unique and significantly distinct combinations of the levels of each of four tested target RNA sequences. The most pronounced was an increase of the number of DWV-like genomes with the VDV-1 CP and the DWV NS sequences in the group VH pupae compared to those in other treatment groups. In comparison with the control group C, the group VH had about a 6,000-fold increase in the VDV-1 CP and about a 26,000-fold increase in the DWV NS. When compared with the group NV, the group VH showed lower relative increases (312-fold for VDV-1 CP and 2500-fold for DWV NS, both with P<0.0001) because an increase of VDV-1 CP and DWV NS also took place in group NV compared to the group C (
Varroa β-actin, R
Varroa β-actin F
Primer descriptions are given as follows: target (position in DWV or VDV-1 nucleotide sequence), polarity (F, forward; R, reverse). GenBank accession numbers used to express primer positions are AJ489744 (DWV), NC—006494 (VDV-1), and AB242568 (Varroa destructor β-actin mRNA).
The strain-specific real-time PCR also showed that the oral acquisition of the DWV-like viruses present in a Varroa-infested hive results in changes in virus diversity. The group NV bees had statistically significant 27-fold increase of the VDV-1 CP (P=0.0217) and 10-fold increase of DWV NS (P=0.0314) compared to the Varroa-free control group C.
A significant difference in DWV diversity was observed between the groups NV and VL. Although both these groups contained similar overall levels of the DWV-like viruses (
To determine if increase of the VDV-1 CP/DWV NS recombinants in the Varroa-infested pupae was linked to the increased levels of these recombinants in associated Varroa mites, the strain-specific real-time PCR was carried out in the mite samples. Surprisingly, the viral levels in the mites did not show strong correlation with viral levels in the corresponding bee pupae. The highest statistically significant (P<0.05) correlations were found for the total DWV levels determined using universal NS primers across all samples of the VL and VH groups (r=0.6025) and for the VDV-1 NS levels across the samples of both VH and VL groups (r=0.5342) and for the samples of the VH group (r=0.7850). Notably, no significant correlation was found between the levels of VDV-1 CP and DWV NS in bees and corresponding mites.
The real-time PCR quantification of DWV and VDV-1 targets suggested that in the bees of group VH the predominant DWV variants were recombinant genomes with VDV-1 CP and the DWV NS, but this approach did not yield information about the type of these recombinants.
Previously it was found that the recombination points could be located in different positions in the 4500-5500 nt region of the DWV genome. Therefore, sequence analysis of this region was required to provide data on phylogenetic relationships between the DWV strains in the treatment groups, including information on the types of the VDV-1/DWV recombinants. The cDNA fragments corresponding to the central regions of the DWV genome were amplified by nested RT-PCR in the bee or Varroa mite sample pools for each treatment group, as well as in the pooled bee sample from the Varroa-infested colony used in the transfer experiment. The outside primers were designed to amplify known DWV strains including VDV-1 and DWV. For each first round reaction four second round amplification reactions were carried out using the specific primers, which allowed VDV-1-type and DWV-type CP and NS to be distinguished, allowing all potential combinations to be amplified.
Recombinant fragments amplified with the VDV-1 CP and the DWV NS primers were detected in the pupae of all treatment groups as well as in the mite samples, being most abundant in the group VH and the Varroa-infected hive sample. The full-length VDV-1 and DWV sequences were detected in treatment groups C and NV and the Varroa-infested colony sample. Notably, no recombinant sequences with DWV CP and VDV-1 NS parts were detected in any sample. The PCR fragments were cloned into a plasmid vector and 8 to 18 individual clones per treatment group were sequenced. In total 93 individual clones were sequenced, which were aligned with 12 published full-length and nearly full-length sequences of DWV strains including DWV, KV, VDV-1 and DWV-VDV-1 recombinants and used for phylogenetic analysis. The 1330 nt alignment, corresponding to the region 4926-6255 nt of the type DWV strain (GenBank Accession No. AJ489744) had 27% divergent positions (23% of divergent positions among the sequences produced in this study). The high proportion of divergent positions allowed us to generate a robust phylogenetic tree of the central genomic regions of DWV strains.
The sequences derived from the bees of the treatment groups C, NV, VL and the Varroa-infested colony appeared in all major clusters of the phylogenetic tree, including DWV, VDV-1 and VDV-1/DWV recombinant clusters,
Surprisingly, we found that Varroa mites harboured wider diversity of DWV variants than the pupae on which they were feeding. Although all sequences from the group VH pupae were in a single clade “VDV-1/DWV Recombinant-C” (
Genomic RNA of DWV and DWV-related viruses was amplified using nested PCR to generate overlapping DNA fragments spanning the entire viral genome. Primers used were directed against conserved regions of the DWV genome and were designed to amplify all known DWV-like viruses, including DWV, VDV-1, KV. Uncloned cDNA populations were sequenced commercially (GATC Biotech.) using standard library preparation methodologies. Paired end sequence data (˜5 million reads per cDNA fragment) was assembled using tview and Bowtie by alignment with reference sequences from Genbank (DWV reference NC—004830) and previous recombinant cDNAs sequenced by the inventors (Moore et al. 2011 J. Gen. Virol. 92: 156-161, Genbank HM063437 and HM063438).
Sequence reads were used to generate a consensus sequence for the majority DWV-like virus present in the preparation. Two predominant forms exist, differing only in their 5′ non-coding regions. The predominant forms were both recombinants in which the structural proteins were more closely related to VDV and the non-structural proteins to DWV. The two recombinant forms had either DVD-like or VDV-like 5′ non coding regions. The sequences for the two predominant forms are given in SEQ ID NOs: 1 and 2.
Sequence variation was determined by analysing all aligned reads showing mismatches to the consensus and templates sequences. ˜81% of all reads in the library were DWV-like. The overall variation of the predominant virus sequence was determined as 0.99051%.
The small RNA (15 to 40 nt) fractions were isolated from the bee total RNA samples (which were used for real-time PCR), pooled according to the treatment groups and used for small RNA library preparations and small RNA sequencing. Libraries which contained 15 to 25 million reads were aligned to the reference DWV and VDV-1 sequences using Bowtie (Langmead B et al. Genome Biol 10:R25.).
DWV and VDV-1 specific siRNAs of both polarities were present in all treatment groups. The most abundant siRNA species were 21 and 22 nt in length, consistent with previous reports on siRNA in insect and suggesting that Dicer operated bees in all treatment groups. It is generally believed that in insects siRNAs are produced by processing of dsRNA replication intermediates and no amplification of siRNA takes place due to the lack of host-encoded RNA-dependent polymerase gene(s). The presence of the DWV-specific 21-22 nt siRNA derived from both positive and negative strands was an indication that replication of these viruses takes place in developing bees of all groups, including those with the low levels of DWV-like viruses (“Control” C, and “No Varroa” NV).
Most of the reads were 21 or 22 nt long, with the three to four times more reads being produced from positive RNA strand than from negative RNA strand (
Pools of siRNA were selected which were unambiguously matching either DWV or VDV-1 reference sequences and profiles of the DWV-type and VDV-1-type siRNA for the libraries derived from the pooled samples of the groups C, NV, VH, and VL were created and analysed (
Honeybees (Apis mellifera) can be engineered to express foreign sequences, stably integrated into their genome, by retroviral vector-mediated infection of spermatozoa and instrumental insemination. Transgenic honeybees, resistant to acute infection with pathogenic viruses such as deformed wing virus (DWV), may be created that express regions of the virus genome engineered in such a way that they induce small interfering RNA (siRNA or RNAi) molecules. For this resistance to be inheritable the virus sequences will need to be engineered into the genome of the honeybee.
To introduce novel genetic material into the germline of honeybees it will be necessary to manipulate a queen to lay fertilized eggs containing the novel genetic material. Sperm mediated gene transfer (SMGT) in honeybees has been achieved in the art by mixing plasmid DNA with spermatozoa used for instrumental insemination. The sequences introduced were stably maintained for several months but did not integrate into the genome. The introduced sequences were also passed on to progeny raised from the manipulated queen. It appears as though a limitation to this approach is the absence of genome integration.
In the present application, to circumvent this limitation, retroviral mediated infection of honeybee spermatozoa is used. This approach has been used to generate transgenic animals from other species, such as zebrafish. Retroviruses are RNA viruses that integrate into DNA during their replication cycle. Virologists have developed methods to create novel retroviruses carrying different ‘payloads’ e.g. foreign genes. To ensure expression of these foreign genes a suitable promoter must be present.
A functional honeybee promoter from a heat shock protein (hsp70) has been identified and cloned (SEQ ID NO: 8). This promoter has been demonstrated to function upon microinjection of bumblebee (Bombus) pupae. A recombinant retrovirus genome (designated plasmid 1) has been generated based upon a standard lentivirus (such as human immunodeficiency virus (HIV), several of which are available commercially) which contains the hsp70 promoter adjacent to a luciferase reporter gene. The sequence of this recombinant retrovirus genome is given in SEQ ID NO: 10. Whole body transfection of bee larvae with the recombinant retrovirus has been shown to result in luciferase expression within the bee larvae (
To generate recombinant retrovirus particles suitable cells e.g. mammalian HeLa cells, are co-tranfected with plasmid 1 together with a plasmid (designated plasmid 2) encoding a suitable envelope glycoprotein (for example, the envelope glycoprotein from vesicular stomatitis virus (VSV) which (a) is usually used to create such recombinant retrovirus particles and (b) is known to attach to and allow infection of insect cells) and a plasmid (plasmid 3) encoding the gag and pol genes of a suitable retrovirus. With the exception of the newly-identified honeybee promoter, which is disclosed herein for the first time, this technology is all known in the art and available commercially.
The resulting recombinant retrovirus is mixed with donor spermatozoa harvested from suitable drones and used for retroviral vector-mediated infection of spermatozoa and instrumental insemination of virgin honeybee queens. The resulting fertilized queens are used to establish small colonies and the progeny screened for the presence of introduced foreign DNA (e.g. luciferase).
White eye pupae (day 12-13 of development) maintained in vitro (as described in Mockel et al. 2011 J Gen Virol 92: 370-377) were directly injected with virus particles purified from groups C, NV and VH pupae according to Moore et al (2011). The proportion of the DWV- and VDV-1-type CP coding regions in the inocula and injected pupae (following incubation to the blue eye stage for 3 days) were determined by qRT-PCR using strain-specific primers to the CP and universal primers to the NS region. All preparations contained higher and broadly similar levels of VDV-1-like CP coding regions. In contrast, the amount of DWV-like CP coding regions was much higher in the virus preparation from the group C pupae (where it accounted for ˜12% of the population) than from either the NV or VH group pupae (
Next generation sequencing (Illumina paired end reads) was conducted to more completely characterise the group C inocula and the viruses present in pupae injected with the group C virus. The composition of the inoculum, as determined by qRT-PCR and subsequent MosaicSolver analysis of the NGS reads, were in close agreement and consisted of 12.5% DWV, 42% VDV-1 with the remainder being a VDV-1 CP-encoding recombinants with a DWV-like NS region (
Pupae inoculated with group C virus exhibited a marked reduction in the DWV content (1%) and a concomitant statistically significant increase in recombinant forms of the virus (70% of the total) that were characterised by the presence of VDV-1 CP coding region and DWV-like NS regions (
Host-pathogen interactions can be broadly divided into asymptomatic or symptomatic infections [1]. In the former, the absence of symptomatic disease is typically due to restricted pathogen replication, which reduces the opportunities for horizontal transmission within its host population. Conversely, prolonged survival of the infected host increases the likelihood of vertical transmission of the pathogen[2]. In contrast, symptomatic infections are typically characterized by high levels of pathogen replication, with consequent enhanced virulence, thereby maximizing horizontal transmission [1-4]. The ‘lifestyle choice’ of asymptomatic or symptomatic infection is determined by multiple factors including the duration of host-pathogen co-evolution, host physiology and anti-pathogen responses, routes of transmission and environmental factors. Evolutionary changes in pathogen virulence may be triggered by changes in pathogen-host assemblages [5]. In the case of multi-host pathogens with interspecies transmission, a pathogen's virulence may increase following introduction of a second host, when the constraint on pathogen virulence in a given host is removed [6].
The European honeybee (Apis mellifera) is the predominant managed pollinating insect and delivers economically important pollination services for agriculture which are estimated to add ˜$40 bn globally to crop value/annum [7]. Factors that influence colony health and viability are therefore important for colony survival and pollination performance. In addition to the bacterial foulbroods, the most important diseases of A. mellifera are caused by a range of viruses many of which are vectored by the ectoparasitic mite Varroa destructor when feeding on honeybee haemolymph. Varroa is believed to have expanded its host range from Apis cerana to A. mellifera during the first half of the 20th century and subsequently spread to all beekeeping regions of the world with the exception of Australia [8-11].
Deformed wing virus (DWV), a picorna-like single-stranded, positive-sense, RNA virus [12,13], is present in the majority of honeybee colonies [10]. DWV is closely related to Varroa destructor virus type 1 (VDV-1) [14]. Their recombinants [15,16] and Kakugo virus (KV) [17], which together exhibit at least 84% nucleotide identity, can be considered as strains of the same virus (henceforth we use the term DWV to refer to this related group of viruses). In the absence of Varroa, DWV generally causes asymptomatic infection and is present at low levels in honeybees. In contrast, in Varroa-infested colonies, mite-exposed pupae can exhibit very high DWV levels which may result in impaired development of the general adult honeybee and increased mortality [10,13]. The mechanisms underlying the transition of DWV from a relatively benign virus to a major honeybee pathogen in the presence of Varroa remain unclear. Two possibilities, not mutually exclusive, have been proposed: suppression of honeybee antivirus defences by Varroa mites which allows the virus to proliferate [18,19], and a Varroa-driven selection of particular DWV genotypes, potentially due to replication in the mite [15,20].
Previous studies using functional or gene expression analysis have produced contradictory conclusions on the impact of Varroa on the immune responses of honeybees. Initial reports indicated that Varroa-exposed honeybees were immuno-compromised [18,19], although later transcriptome analysis found little or no effect on genes implicated in insect immunity [21,22]. Additional studies have shown down-regulation of a honeybee NF-κB transcription factor [23]. Recent reports have implicated the Drosophila Toll, Imd and Jak-Stat signalling pathways in controlling RNA virus infection [24] and RNA interference (RNAi), which has long been considered the major antiviral mechanism in insects [25], has recently been associated with controlling the persistence of RNA virus infections in Drosophila [26]. It was therefore possible that high levels of DWV in Varroa-exposed honeybees could be the result of a suppression of these antivirus responses and so warranted further analysis.
We have previously demonstrated that Varroa infestation is associated with the accumulation in mite-exposed pupae of a particular subset of DWV-like viruses [15]. These recombinant forms (RF) are predominantly comprised of genomes with structural and non-structural coding regions that most closely align with VDV-1 and DWV respectively. The organisation of these recombinants suggests that, as with other picorna-like viruses, DWV likely has a modular genome, with a 5′ untranslated region (5′-UTR) driving translation of the structural or capsid (CP) and non-structural (NS) ‘modules’ [15]. We hypothesised that such recombinants were transmitted more efficiently between Varroa and honeybees, resulting in their amplification to the markedly elevated levels observed in Varroa-parasitized pupae (about 1000 times higher than in unexposed pupae). In recent complementary studies, changes in the composition of the DWV population over a large temporal and spatial scale following Varroa infestation were reported for honeybees colonies following accidental introduction of Varroa into the Hawaiian islands [20]. The introduction of Varroa was associated with a marked restriction in DWV diversity measured in the pooled honeybee samples collected from the Varroa-infested colonies, although the precise identity of the dominant virus was not determined [20].
In the present study we devised a novel experimental system to specifically test two hypotheses on the role of Varroa in the development of high-level DWV infection in the honeybee, namely that the mite (i) amplifies and transmits virulent genotypes of DWV, and (ii) suppresses antiviral responses, including immune signalling pathways and RNA interference. The experimental procedure included exposure of Varroa-naïve honeybees to mites and their associated DWV payload together with the per os in-hive horizontal transmission. The use of Varroa-naïve honeybees from a Varroa-free region allowed us to monitor changes in DWV diversity and loads, as well as potential antivirus responses in the honeybee responses, following exposure to the viral genotypes associated with Varroa infestation. Importantly, we analysed immune responses and viral load/diversity in individual mite-exposed and -unexposed pupae, rather than in pooled samples. This allowed us to stratify individual responses into four distinct experimental groups, characterised by Varroa exposure and viral load, that clearly correlated with characteristic changes in the transcriptome and virus population diversity. In addition, we recapitulated the exposure of Varroa-free honeybees to DWV by direct injection and analysed virus diversity in bees of a colony with long-established Varroa infestation.
Our results indicate that a virulent recombinant form of DWV, while transmissible orally, only replicates to high levels when directly inoculated into honeybee haemolymph—by Varroa or experimental injection. This results in massive reduction of DWV diversity in bees with high virus levels, both in the Varroa-exposed pupae and newly emerged bees with symptomatic deformed wing disease. Significantly, the same virulent recombinant form of DWV reached the highest levels in mite-exposed pupae and in adult bees exhibiting characteristic deformed wing symptoms. Although exposure to Varroa resulted in changes in expression of a number of immune-related genes, the roles of which should be further explored, we demonstrate that it is the route of virus acquisition that is responsible for the amplification of this virulent form of DWV in a Varroa-infested colony.
Worker honeybee larvae from a Varroa-free colony (sourced from a region with no historic contacts with or presence of Varroa) were moved in a frame transfer experiment to a Varroa-infested colony. The larvae were subsequently exposed through feeding to DWV strains circulating in the infested colony from day 4 until the cells were capped at day 9 (all times relative to egg laying;
We assessed the total levels of DWV viruses in 80 individual pupae by qRT-PCR using a primer pair for a conserved polymerase-coding region, designed to detect all known DWV strains, including DWV, VDV-1 and KV (Tables 1 and 5). The real-time PCR Ct values showed a clear bimodal distribution, with distinct low- and high-levels of DWV (p<10−16;
We sampled eight honeybee pupae selected at random from each of the four groups (C, NV, VH and VL;
We used a two-colour dye-balanced loop design microarray [30,31] to determined the genome-wide transcriptional profile using RNA extracted from the 32 samples defined above, (8 pupae from each experimental group). The oligonucleotide expression array contained probes to all protein-coding transcripts of A. mellifera [32], as well as probes to all known viral and fungal pathogens of honeybees, including distinct DWV and VDV-1 probes. After array normalization, differentially expressed (DE) genes were determined for each contrast between experimental groups (
There were high levels of commonality and additivity for DE genes in the contrasts considered (
To obtain insight into the functional consequences of DE gene expression we carried out Gene Ontology (GO) analysis, focusing on the GO Biological Process (BP) [33]. A number of overrepresented GO BP terms related to cell division were associated with DE genes in the C to NV contrast, while those related to regulation of various cellular processes were associated with the DE genes in the NV to VL contrast (
Although no significant changes in expression of genes associated with the RNAi response (e.g. Argonaute, Dicer) [25] were observed in the microarray analysis, there could be post-transcriptional effects on RNAi generation. We therefore analysed the DWV-related RNAi population and compared it with the levels and identity of virus in pupae from the four experimental groups (
All RNA libraries analysed contained similar proportions of host-encoded miRNA reads, 12 to 18% of the total (Table 6), indicating both successful isolation of small RNA libraries and broad equivalence of the pooled sample sets. DWV- and VDV-1-specific siRNAs of both polarities were present in all treatment groups. DWV- and VDV-1-specific siRNAs could originate from either DWV or VDV-1, or from the previously reported [15] recombinants between these parental viruses (
We and others have previously reported changes in virus diversity at the population level [20] and the predominance of particular virus recombinants (i.e. a reduction in diversity) in honeybee pupae exhibiting high viral loads [15]. To quantify both viral load and diversity in individual honeybee pupae and their associated Varroa mites we used generic NS qRT-PCR primers or primer pairs specific for DWV or VDV-1 CP or NS coding regions (Table 1). We quantified the total virus levels (
As already indicated (
DWV- or VDV-1 specific qRT-PCR analysis demonstrated only a weak correlation with virus levels in the corresponding honeybee pupae (
The dominance of recombinant viruses bearing VDV-1 CP and DWV NS coding regions in the VH group was strongly suggested by qRT-PCR (
The resulting dendrogram contained six distinct clusters, one each for non-recombinant DWV- or VDV-1-like sequences, together with four different VDV-1/DWV recombinant forms (designated RF1-RF4;
One interpretation of the near-clonality of viral sequences in the VH group was that these were the only ones carried, and hence transmitted, by the mite. However, with the exception of the non-recombinant DWV cluster, which was not detected in the mite, the 32 viral sequences obtained from Varroa were widely distributed within the dendrogram (open symbols in
Inoculation of Honeybees by Injection into Haemolymph Results in Preferential Amplification of Specific VDV-1/DWV Recombinants
We directly injected white eye pupae (day 12-13 of development) maintained in vitro (as described in [39]) with virus particles purified from groups C, NV and VH pupae as described previously [15]. As before, we determined the proportion of the DWV- and VDV-1-type CP coding regions in the inocula and injected pupae (following incubation to the purple-eye stage for 3 days) by qRT-PCR using strain-specific primers to the CP and universal primers to the NS region. Virus preparations from groups NV and VH contained higher and broadly similar levels of VDV-1-like CP coding regions. The amount of DWV-like CP coding regions was much higher in the virus preparation from the group C pupae (where it accounted for ˜12% of the population) than from either the NV or VH group pupae (
We additionally conducted next generation sequencing (Illumina paired-end reads) to comprehensively characterise the group C inocula and the viruses present in pupae injected with the group C virus. The composition of the inoculum, as determined by qRT-PCR and subsequent MosaicSolver [40] analysis of the NGS reads, were in close agreement and consisted of 12.5% DWV (in excellent agreement with the qRT-PCR-determined level, see above), 42% VDV-1 with the remainder being VDV-1 CP-encoding recombinants with a DWV-like NS region (
Independent Verification of DWV Diversity Reduction by Deep Sequencing of the Honeybees from a Varroa-Infested Colony.
The sequence analysis of DWV in Varroa-exposed pupae (
The NGS reads were aligned to reference DWV and VDV-1 sequences (GenBank Accession numbers GU109335 and AY251269 respectively), and the pileup profiles were analysed. The proportions of DWV and VDV-1 reads in the libraries (each containing about 10 million reads) showed a bimodal distribution and were either very high (from 7.41% to 83.87%,
The introduction and global spread of the parasitic mite Varroa destructor in honeybees (Apis mellifera) has had significant impact on the health and survival of infested colonies [8,13]. Colony losses associated with Varroa are predominantly attributed to the RNA virus payload vectored by the parasite and transmitted when the mite feeds on honeybee haemolymph [10,45]. Although Varroa were reported to vector at least 5 RNA viruses, the picoma-like Iflavirus deformed wing virus (DWV) is of particular interest and importance; deformed wing disease is associated with mite infestation [9] and high levels of DWV exacerbate overwintering colony losses [10]. Furthermore, of the viruses analysed, only DWV levels increased upon introduction of Varroa to Hawaii [20]. Here we show that virus levels exhibited a bimodal distribution in developing honeybee pupae—low in the absence but generally high following mite exposure—with high viral levels associated with emergence of a distinct recombinant form of the virus. This complements observations of the effects of the introduction of Varroa at the regional scale (e.g. Hawaii), which is associated with a dramatic reduction in DWV variation and the emergence of dominant strains over a scale of months or years [20]. Our study also strongly suggests that DWV is widespread and present in all UK honeybees, in Varroa-free and Varroa-infested colonies. By using different molecular approaches including RNA-seq, cloned cDNA sequencing and qRT-PCR, we demonstrated that, contrary to the recently reported figure of 36% DWV presence in the UK honeybees [46], DWV was present in all n=250 tested honeybees, including those from a Varroa-free region (n=47) (
To better understand the interactions of the honeybee host, Varroa and DWV that account for the observed emergence of presumed pathogenic, near-clonal strains of the virus we analysed events at the level of individual honeybee pupae. We reasoned that the large scale global, temporal and population changes observed reflect the cumulative outcome of interactions that occur at the level of individual pupae within the colony. We further reasoned that analysis at this level would allow two hypotheses accounting for the high levels of overt deformed wing disease beekeepers associate with heavy Varroa infestation to be tested. These hypotheses, themselves not mutually exclusive, are that; a) Varroa suppresses the antiviral honeybee defences so allowing unrestricted DWV replication and b) the presence of Varroa results in the selection and transmission of particular pathogenic variants of DWV, resulting in serial amplification within a mite infested colony. We tested these hypotheses by transferring honeybee larvae from a Varroa-free to an experimental Varroa-infested environment, stratified pupae of a standardised age in terms of mite-exposure and viral load, and investigated transcriptome and RNAi responses of the host and the virus population in individual mites and associated pupae.
We demonstrated that Varroa-naïve honeybee larvae (group C) transferred into a Varroa-infested colony (effectively mimicking the exposure to oral and mite-transmitted DWV during larval and pupal development) can, after incubation for 6 days, be stratified into three distinct pupal groups by Varroa exposure (presence or absence) and DWV level (high or low). High virus levels were not observed in the absence of Varroa. Group NV comprised pupae from capped cells free of Varroa. As Varroa enters the cell immediately prior to capping [27], we assume this group only acquired viruses per os during larval feeding. The VH and VL groups contained Varroa within the capped cell with evidence of Varroa feeding on the pupae—including nymphal forms present and signs of abdominal piercing. These groups of pupae harbour strikingly different DWV populations: the VL group having low viral levels and high diversity that are not significantly different from the C and NV groups, whereas the VH group carries 1,000-10,000 times the viral loads of a single phylogenetic type (
Previous analyses of transcriptome or immune response changes in Varroa-exposed honeybees have produced contradictory results, perhaps due to the analysis of pooled individuals and/or pupae of different ages. These confounding influences may have obscured the marked changes in gene expression that we observed in response to either mite or viral pathogens, as emphasised by the transcriptome differences observed in the Varroa-associated pupae in groups VL and VH, which respectively exhibited 493 (˜5% of transcriptome) and 951 (˜9%) significantly differentially expressed (DE) genes when compared with the control group C (
Gene Ontology (GO) analysis gave additional insights into the transcriptional responses in honeybees following experimental treatments (
We analysed changes in expression of known and presumed immune-related genes (Table 3,
Varroa exposure (VL or VH groups) resulted in down-regulation of several putative components of the honeybee Toll signalling pathway [51], including two Toll receptor orthologs (GB10640, GB17781), CLIP-domain protease spirit (GB14044) and the Toll receptor ligand spatzle (GB15688). In addition, spatzle was down-regulated when the VH group was compared against the other experimental groups, suggesting down-regulation of this gene may be a response to the elevated virus levels in the group VH, rather than the presence of Varroa per se. Toll signaling pathways are implicated in antiviral resistance to the RNA virus Drosophila virus X [54], possibly controlling proliferation of haemocytes which, because of their involvement in phagocytosis, play a central role in insect immunity [54,55]. We also observed down-regulation of a Tetraspanin 68C (Tsp68C) ortholog (GB16002, GB13670), a cell surface membrane scaffolding protein previously implicated in receptor modulation during hemopoiesis [56], an ortholog of pannier (GB 19895), a GATA transcription factor required for differentiation of plasmatocytes (which resemble the mammalian macrophage lineage [57]), and a serrate ortholog (GB15155), a membrane ligand for the Notch receptor implicated in differentiation of haemocyte-related crystal cell precursors which function in pathogen defence via melanisation [58]. These transcriptome changes may help explain functional studies in which salivary secretions from Varroa mites damage moth caterpillar haemocytes [59] and suggest that Varroa-mediated depletion of haemocytes, a key component of the immune response of insects [60-62], may contribute to enhanced susceptibility to DWV and other viruses. Interestingly, we also observed suppression of the Friends-of-GATA transcription factor U-shaped (ush) ortholog (GB 16457), in the C to VL and NV to VL contrasts. Drosophila ush is reported to antagonise crystal cell development [63,64], implying that the low level of virus accumulation in the VL group may be due to elevated numbers of crystal cells resulting from ush down-regulation.
Although by definition descriptive, transcriptome analysis of pupae stratified according to Varroa and virus-exposure, also provides direct insights into possible pathogenic mechanisms. In the contrast C to VH we observed differential expression of orthologs of five Drosophila homeobox genes (summarised in Table 8) encoding transcription factors which are involved in insect development [65]. Most of these DE genes are reported to be expressed at early pupal stages and involved in abdomen (Abdominal B), appendage (apterous) or brain development (extradenticle). This may explain previously reported developmental abnormalities in the honeybee that are associated with high DWV levels at the pupal stage [13] and warrant further investigation to potentially determine the molecular mechanism underlying DWV pathogenesis.
Notwithstanding the absence of significant changes in gene expression of key components, such as Dicer and Argonaute, of the RNAi response—the major antiviral mechanism in insects [25]—we explored the relationship between DWV-like virus levels and the corresponding siRNA populations. In particular, we sought to investigate if high levels of DWV in VH honeybees was associated with the limited accumulation of virus-derived siRNA, implying the virus may express an siRNA suppressor as, for example, demonstrated in Alphavirus infection of mosquitos [66]. Although DWV- and VDV-1 specific siRNAs were recently detected in adult honeybees [67,68], these studies could not show if RNAi is involved in suppression of the virus, because viral genomic RNA levels were not quantified. Analysis of siRNAs in the honeybees of the frame transfer experiment showed that the predominant DWV- and VDV-1-specific siRNAs were 22 nt in length with genome sense strand-specific siRNAs present at a 3-4 fold excess over antisense. This was consistent both with the presence of replicating DWV-like viruses in all experimental groups and with the known activity of Dicer in other insects [25] and strongly suggests normal functioning of Dicer in honeybees [69]. As insects do not amplify siRNA populations [25], it was unsurprising that virus-specific siRNA levels were broadly proportional to the level of viral genomic RNA determined by qRT-PCR; the C and NV groups exhibited ˜103 times less viral genomic RNA than the VH group and normalised siRNA levels differed by ˜770 times (Table 6). The exception to this was the siRNA response in the VL group which was ˜5 times higher proportionally than the level of VL virus genomic RNA (Table 6). The relationship between the levels and compositions of the viral genomic RNA and virus-derived siRNA may be altered by differences in targeting of the individual components of DWV-like virus population by the honeybee RNAi machinery, as observed during West Nile virus infection of mosquitos [70]. Although the presence of virus-specific siRNAs does not necessarily correlate with effective silencing—viruses may encode late-acting suppressors such as the Argonaute-inhibiting 1A protein of cricket paralysis virus [71]—the robust siRNA response in the VL group may contribute to suppression of DWV replication and the differences between this response and that observed in the VH group may be a fruitful area for further analysis.
Genetic Diversity of the DWV Population is Determined by Route of Transmission Rather than Preferential Amplification of Virus in Varroa
The introduction of the parasitic Varroa mite elevates the level of DWV-like viruses [20], amplifies particular recombinant forms (RF) bearing the capsid determinants of VDV-1 and non-structural genome region from DWV [15,16] and dramatically reduces the diversity of DWV-like viruses in a population [20]. Using complementary approaches including strain-specific qRT-PCR and sequencing together with next generation sequencing of the virus genome and host siRNA response to infection, we analysed individual pupae exposed to DWV during larval feeding and following mite exposure, and recapitulated horizontal transmission of virus by Varroa using direct injection.
The C, NV and VL exposure groups all carried low viral loads and exhibited high virus diversity (
We reasoned that there were two possibilities that might account for the marked amplification of DWVV in the VH group pupae. Firstly, the mite may have delivered a high dose of one specific recombinant form, perhaps reflecting its preferential replication in the ectoparasite. Secondly, we considered that DWVV might have a growth advantage when inoculated into haemolymph by Varroa (potentially in addition to the preferential amplification in the mite). To distinguish between these possibilities we sequenced qRT-PCR amplified viral RNA from mites co-isolated from capped cells containing group VL and VH pupae. We also investigated the consequences of Varroa-independent mechanical virus transmission by direct injection of mixed virus preparations to Varroa-naïve pupae and subsequent monitoring of virus levels and diversity.
Individual Varroa mites contained a diversity of DWV-like sequences that were well distributed throughout the phylogenetic tree of virus sequences from pupae (square symbols,
Although further studies will be required to determine whether sampling stochasticity accounts for the apparent absence of non-recombinant DWV in Varroa, together these results suggest that at least at the level of the entire mite—there is no selection, either by absolute presence or replication capability, for the DWVV RF that accumulates to high levels after mite exposure in VH group pupae.
Since the diversity of virus present in Varroa indicates that the near-clonal virus population in the VH group is not due to the mite delivering either a restricted virus type or to elevated levels of DWVV in the mite we went on to inoculate pupae with a mixed virus population prepared from group C pupae and characterised the resulting virus population after three days. Recipient purple-eye stage pupae contained high virus loads which, compared with the inocula, had markedly reduced levels of DWV-like virus and elevated levels of a VDV-1/DWV recombinant (
We propose that the strikingly elevated levels and associated restricted diversity of DWVV (RF4-type;
a Honeybee gene ID according to the Apis mellifera Official Gene Set 1 [30]. Honeybee immune-related genes included in the analysis were either those described in [33] or the honeybee homologues of the Drosophila melanogaster genes associated with Gene Ontology Biological Process term “Immune System Process” GO:0002376.
b Drosophila melanogaster homologue showing highest similarity in BLAST.
c Drosophila melanogaster gene name according to FlyBase.
d Pathway or category of gene if known (Toll-Toll signalling pathway; SP-serine protease; Scav-Scavenger receptor A; Notch-Notch signalling pathway; JNK-JNK signalling pathway; J-ST-JAK-STAT signalling pathway; IG-IG Superfamily Genes; GATA-GATA transcription factor), Immune-Immune system process gene.
The single read libraries were aligned using Bowtie to the Apis mellifera miRNA and to the reference full-length DWV and VDV-1 sequences, GeneBank Accession numbers GU109335 and AY251269 respectively.
Without proper management Varroa has a devastating effect on honeybee colony viability and consequent honey production and pollination services. We show here that the markedly elevated levels of DWV-like viruses in Varroa-exposed honeybee pupae are likely attributable to the direct inoculation of a specific virus, DWVV, by Varroa to haemolymph. Repeated cycles of Varroa-replication within an infested colony would preferentially amplify DWVV, potentially resulting in it becoming the predominant virus present, transferred both by Varroa and per os. Further studies will be required to determine whether such a virus, if sufficient were ingested, would also cause symptomatic infection. Oral susceptibility to a virulent form of DWV may also explain reported cases of deformed wing disease symptoms seen in Varroa-free colonies in Hawaii [20] and Scotland (Andrew Abrahams, pers. comm.), but may also reflect genetic variation and the presence of particularly susceptible pupae in the colony.
Our study demonstrates that a proportion of Varroa-exposed pupae (the VL group) do not exhibit elevated levels of the near-clonal DWVV recombinant (
This study was based around an experiment in which a brood frame containing newly hatched larvae from a Varroa-free colony was introduced into a Varroa-infested colony. The larvae were left to develop within the Varroa-infested colony, and pupae were collected 11 days later from capped brood cells at the purple eye stage and analysed using a range of molecular methods. The Varroa-free honeybee (Apis mellifera) colony with a naturally-mated one-year-old queen was imported from Colonsay, Scotland, an island with no historic reports of Varroa incidence and no imports of honeybees from Varroa-infested areas. This allowed us to exclude the presence of DWV strains associated with Varroa mite infestation. As a source of Varroa mites and the mite-associated DWV strains, we selected a Warwickshire honeybee colony, heavily infested with Varroa and having high DWV levels in honeybees and mites. The Varroa-free and Varroa-infested colonies were contained in separate mesh flight cages (dimensions: 6 meters long, 2.5 meters wide, 2 meters high) and maintained on an artificial diet of sugar syrup and pollen. The pollen was imported from Varroa-free Australia to exclude possible contamination with Varroa-associated viruses through foraged food and was pre-screened by PCR before use for DWV-like viruses. In order to minimise possible effects on honeybee gene expression due to the differences in nutrition, both the control Varroa-free and the Varroa-infested colonies were maintained in flight cages in the same apiary (at the University of Warwick, UK) and were fed on the artificial diet for two months before the start of the frame transfer experiments. Neither colony was treated with miticides.
The experimental infestation, summarised in
The pupae and the Varroa mites associated with each infested pupa were individually snap frozen in liquid nitrogen immediately after being removed from brood cells and stored at −80° C. prior to total RNA extraction. For total RNA extraction, whole individual honeybee pupae were ground to fine powder in liquid nitrogen, and half of the powder used for RNA extraction, carried out using 1 mL of Trizol Reagent (Invitrogen) according to the manufacturer's instructions. Total RNA extraction from Varroa mites was carried out using RNeasy spin columns (Qiagen RNeasy Plant Mini kit).
Virus purification from honeybee material and extraction of the viral genomic RNA from virus particles were carried out as described previously [15].
For genome-wide analysis of the honeybee transcriptome total RNA preparations from eight individual honeybees from each of the four experimental groups (32 honeybees in total) were purified further using RNeasy Plant Mini kit spin columns (Qiagen). RNA concentration, purity and integrity were assessed using a 2100 Bioanalyzer and an RNA 6000 LabChip (Agilent Technologies). The probe preparation, hybridization and scanning were carried out according to the Agilent instructions, essentially as in [75]. Total RNA extracts from an individual honeybee were used to produce Cy3- and Cy5-labelled aRNA samples using a Low Input RNA fluorescent linear amplification kit (Agilent Technologies), according to the manufacturer's instructions. The Cy3- and Cy5-labelled samples were used in a two-colour dye-balanced loop design [30,31] for a genome-wide analysis of the honeybee transcriptome with the custom expression oligonucleotide microarray. Four slides, each with eight two-channel arrays were employed, allowing two replicates per sample, one green and one red. Different treatment groups were allocated to the green and red channels in each array; the loop design ensured that each sample was indirectly compared with all other samples. The array, in 60K×8 format, included 60 nt oligonucleotides specific to 10,157 transcripts of the Apis mellifera Official Gene Set 1, OGS1 [32]; the array also contained probes to all the honeybee RNA viruses known to date. Each probe was replicated five times to enable robust statistical analysis. Sequences of the probes of the honeybee whole genome expression microarray, a 60-mer oligonucleotidide array based on the Apis mellifera transcriptome (OGS1) and Apis mellifera fungal and viral pathogens (Agilent ID: 027104, SurePrint G3 Custom GE 8×60K), are available in the ArrayExpress database (www.ebi.ac.uk/arrayexpress) under accession number A-MEXP-2251. Following hybridisation, the microarrays were scanned using Agilent Technologies G2565CA Scanner and the fluorescence intensity data were processed using feature extraction software (Agilent Technologies). Cy3 and Cy5 fluorescence intensities for each spot were measured as values of green and red pixels respectively. The details of the array experiment design, sample description, and microarray data are available in the ArrayExpress database (www.ebi.ac.uk/arrayexpress) [76] under accession number E-MTAB-1285. One array failed (assigned to VL green and NV red) leaving 62 channels for final analysis.
For additional confirmation we conducted qRT-PCR analysis with primers specific to Paenibacillus larvae ssp and Melissococcus plutonius (Table 5), the causal agents of American foulbrood and European foulbrood, respectively, which showed that the samples were free of detectable levels of these bacterial pathogens.
The unprocessed intensity scanning values were both within-array and between-array normalized using the linear model based Limma R package [77]. Differentially expressed (DE) genes in all six possible contrasts were found using Limma (via function “lmscFit” incorporating intraspot correlation) and also the R GaGa package for gamma-gamma Bayesian hierarchical modeling [78-80]. A gene was considered as differentially expressed (DE) in a given contrast (using a t-statistic moderated across genes) when the average expression exceeded 6.0, the fold change exceeded 14%, the Limma analysis p-value adjusted for multiple genes was less than 0.05 and the posterior probability determined by GaGa was above 0.6. Microarray results were validated by qRT-PCR using a set of primers for certain honeybee genes and DWV (Tables 1 and 5).
For Gene Ontology (GO) analysis a three-stage process was used. Genes in the latest A. mellifera genome annotation, Amel—4.0 (http://hymenopteragenome.org/beebase), corresponding to genes in A. mellifera OGS1 were found using protein blast. GO terms associated with Amel—4.0 genes were then obtained using Blast2GO [81] with the SwissProt database option. Finally, over- and under-represented GO terms in the sets of DE honeybee OGS1 genes in each contrast were obtained with BiNGO, using a hypergeometric test, a Benjamini and Hochberg FDR correction and a significance level of 0.05 [82].
For Principal Component Analysis (PCA), the significant DE genes in all six contrasts were pooled and ranked by their adjusted p-value. The 60 with the lowest adjusted p-value were selected, all of which appeared in the contrast C to VH; the other contrasts' contributions were 35 (C to VL), 21 (NV to VH), 19 (C to NV), 4 (NV to VL) and 11 (VL to VH) genes. Principal components of the expression profiles across the 62 microarray channels were found and (the first two) plotted using the princomp and biplot functions in R [83].
For high throughput sequencing of small RNA, we pooled equal amounts of the Trizol-extracted total RNA from individual honeybees and isolated the 15 to 40 nt RNA fraction, which was separated using denaturing polyacrylamide gel. The RNA pools were ligated to the oligonucleotide adapters, reverse-transcribed and amplified using the TruSeq Small RNA Sample Prep Kit (Illumina small RNA kit). The libraries were sequenced using the Illumina HighSeq 2000 platform, producing 15-25 million reads per libraries (GATC-Biotech, Germany). The small RNA NGS sequencing data are available in the ArrayExpress database (www.ebi.ac.uk/arrayexpress) [76] under accession number E-MTAB-1671. The reads were cropped to remove adapter sequences and aligned to reference viral and miRNA sequences using Bowtie [37]. Samtools mpileup was used to produce the siRNA and miRNA coverage profiles.
Real-time reverse transcription PCR was carried out essentially as in [15]. In brief, RNA extracts were treated with DNAse, then purified DNA-free total RNA preparations were used as a template to produce cDNA using random primer and Superscript III reverse transcriptase (Invitrogen). The cDNA samples produced were used for real-time PCR quantification of the DWV or host transcripts using SYBR green mix (Agilent Technologies). Oligonucleotide primers are summarized in Tables 1 and 5.
For strand-specific quantification of viral RNA of DWV and VDV-1 types reverse transcription was carried out at 50° C. using Superscript III reverse transcriptase (Invitrogen) and the tagged primers designed to anneal to the negative strands RNA of DWV or VDV-1, primers 389 and 391 respectively (SEQ ID NOs: 50 and 51 in Table 1). The qPCR step was carried out using corresponding DWV or VDV-1 specific primers in negative polarity (Table 1, SEQ ID NOs: 15 and 16) and primer 388 identical to the sequence of the tag (Table 1, SEQ ID NO: 49).
Amplification of the cDNA fragments corresponding to the central region of DWV genomic RNA was carried out by nested PCR using GoTaq PCR mix (Promega) and primers 155 and 156 (SEQ ID NOs: 23 and 24 in Table 1) using the cDNA extracted from the honeybees and the mites, pooled according to their treatment groups. The outside PCR primers were designed to amplify all known DWV-like sequences. For each first round reaction we carried four second round amplification reactions using VDV-1- or DWV-specific primers, 151-154 (SEQ ID NOs: 23 to 26 in Table 1), which allowed distinction of VDV-1-type and DWV-type CP and NS regions, thereby enabling amplification of all potential combinations, even those present at very low levels. The PCR fragments were cloned into pGemT-Easy (Promega) and sequenced using the Sanger dideoxy method. GenBank accession numbers for the reported sequences are JX661628-JX661712 and KC249926-KC249933. The full-length cDNA of DWV, GenBank accession number KJ437447, was amplified by RT-PCR using primers specific to the published termini of DWV and VDV-1 RNA and cloned into the pCR-TOPO-XL vector (Invitrogen) as described in [15]. The sequences were aligned using CLUSTAL X [84], and phylogenetic analysis of the sequences was carried out using the PHYLIP package [85].
For the next generation sequencing of RNA, a series of overlapping cDNA fragments were produced using viral RNA or total RNA preparations using the set of primers designed to the sequences of the genomic RNA conserved among DWV, VDV-1 and KV (Table 1). The fragments were pooled and libraries of paired-end reads (101 nt.), about 5 million per sample, were generated using an Illumina HiSeq 2000 (GATC-Biotech). The virus genomic RNA NGS sequencing data are available in the ArrayExpress database (www.ebi.ac.uk/arrayexpress)[76] under accession number E-MTAB-1675.
The next generation sequencing of the poly(A) RNA fraction (RNA-seq) of the total RNA preparations isolated from the honeybees was carried out using Illumina HiSeq 2000 (GATC-Biotech) protocol, with about 10 million 101 nucleotide-long reads generated for each sample. The RNA-seq sequencing data are available in the EBI Sequence Read Archive [86] under accession number PRJEB5249. This RNA-seq dataset was used to calculate Shannon's diversity index values of DWV populations using the following procedure. First we selected the reads aligning to the reference DWV and VDV-1 sequences (GenBank Accession numbers GU109335 and AY251269 respectively) from the original RNA-seq libraries using Bowtie. To take into account the effect of difference in coverage of low virus levels and high virus level RNA-seq libraries we used two approaches, (i) correction for NGS error for complete libraries ([44]) and (ii) multiple sampling. For the latter we produced five samples of 3285 reads (the lowest number of the viral reads among the libraries), which were aligned using Bowtie to the reference DWV and VDV-1 sequences, and the NGS nucleotide pileups were then generated for each nucleotide position of the reference sequences using samtools. Shannon's diversity index of the aligned nucleotides was calculated for each position in the reference sequence. Then, the average Shannon's index values were calculated for the selected regions in the reference genomes for each sample. The averages values of and standard deviation of five samples were used in the statistical analysis.
Number | Date | Country | Kind |
---|---|---|---|
1306812.7 | Apr 2013 | GB | national |
1314063.7 | Aug 2013 | GB | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/GB2014/051176 | 4/15/2014 | WO | 00 |