Amy L Kistler*1, Ady Gancz2, Susan Clubb3, Peter Skewes-Cox1,6, Kael Fischer1, Katherine Sorber1, Charles Y Chiu1,4, Avishai Lublin5,
Sara Mechani5, Yigal Farnoushi5, Alexander Greninger1, Christopher C Wen6, Scott B Karlene7, Don Ganem1 and Joseph L DeRisi1
Address: 1 Departments of Biochemistry, Microbiology and Medicine, Howard Hughes Medical Institute and University of California, San Francisco, 94143, USA, 2 The Exotic Clinic, Herzlyia, 46875, Israel, 3 Rainforest Clinic for Birds and Exotics, Loxahatchee, FL, 33470, USA, 4 Division of Infectious Diseases, University of California, San Francisco, 94143, USA, 5 Division of Avian & Fish Diseases, Kimron Veterinary Institute, Bet Dagan, 50250, Israel, 6 Biological and Medical Informatics Program, University of California, San Francisco, 94143, USA and 7 Lahser Interspecies Research Foundation, Bloomfield Hills, MI, 48302, USA
Email: Amy L Kistler* - This email address is being protected from spambots. You need JavaScript enabled to view it.; Ady Gancz - This email address is being protected from spambots. You need JavaScript enabled to view it.; Susan Clubb - This email address is being protected from spambots. You need JavaScript enabled to view it.; Peter Skewes- Cox - This email address is being protected from spambots. You need JavaScript enabled to view it.; Kael Fischer - This email address is being protected from spambots. You need JavaScript enabled to view it.; Katherine Sorber - This email address is being protected from spambots. You need JavaScript enabled to view it.; Charles Y Chiu - This email address is being protected from spambots. You need JavaScript enabled to view it.; Avishai Lublin - This email address is being protected from spambots. You need JavaScript enabled to view it.; Sara Mechani - This email address is being protected from spambots. You need JavaScript enabled to view it.; Yigal Farnoushi - This email address is being protected from spambots. You need JavaScript enabled to view it.; Alexander Greninger - This email address is being protected from spambots. You need JavaScript enabled to view it.; Christopher C Wen - This email address is being protected from spambots. You need JavaScript enabled to view it.; Scott B Karlene- This email address is being protected from spambots. You need JavaScript enabled to view it.; Don Ganem - This email address is being protected from spambots. You need JavaScript enabled to view it.; Joseph L DeRisi - This email address is being protected from spambots. You need JavaScript enabled to view it.
* Corresponding author
Published: 31 July 2008
Virology Journal 2008, 5:88 doi:10.1186/1743-422X-5-88
This article is available from: http://www.virologyj.com/content/5/1/88
© 2008 Kistler et al; licensee BioMed Central Ltd.
Received: 30 June 2008
Accepted: 31 July 2008
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract
Background: Proventricular dilatation disease (PDD) is a fatal disorder threatening domesticated and wild psittacine birds worldwide. It is characterized by lymphoplasmacytic infiltration of the ganglia of the central and peripheral nervous system, leading to central nervous system disorders as well as disordered enteric motility and associated wasting. For almost 40 years, a viral etiology for PDD has been suspected, but to date no candidate etiologic agent has been reproducibly linked to the disease.
Results: Analysis of 2 PDD case-control series collected independently on different continents using a pan-viral microarray revealed a bornavirus hybridization signature in 62.5% of the PDD cases (5/8) and none of the controls (0/8). Ultra high throughput sequencing was utilized to recover the complete viral genome sequence from one of the virus-positive PDD cases. This revealed a bornavirus-like genome organization for this agent with a high degree of sequence divergence from all prior bornavirus isolates. We propose the name avian bornavirus (ABV) for this agent. Further specific ABV PCR analysis of an additional set of independently collected PDD cases and controls yielded a significant difference in ABV detection rate among PDD cases (71%, n = 7) compared to controls (0%, n = 14) (P = 0.01; Fisher's Exact Test). Partial sequence analysis of a total of 16 ABV isolates we have now recovered from these and an additional set of cases reveals at least 5 distinct ABV genetic subgroups.
Conclusion: These studies clearly demonstrate the existence of an avian reservoir of remarkably diverse bornaviruses and provide a compelling candidate in the search for an etiologic agent of PDD.
Background
Proventricular dilatation disease (PDD) is considered by many to be the greatest threat to aviculture of psittacine birds (parrots). This disease has been documented in mul- tiple continents in over 50 different species of psittacines as well as captive and free-ranging species in at least 5 other orders of birds [1-5]. Most, if not all major psittacine collections throughout the world have experienced cases of PDD. It has been particularly devastating in countries like Canada and northern areas of the United States where parrots are housed primarily indoors. However, it is also problematic in warmer regions where birds are typically bred in outdoor aviaries. Moreover, captive breeding efforts for at least one psittacine which is thought to be extinct in the wild, the Spix's macaw (Cyanopsitta spixii), have been severely impacted by PDD.
PDD is an inflammatory disease of birds, first described in the 1970s as Macaw Wasting Disease during an outbreak among macaws (reviewed in [3]). PDD primarily affects the autonomic nerves of the upper and middle digestive tract, including the esophagus, crop, proventriculus, ven- triculus, and duodenum. Microscopically, the disease is recognized by the presence of lymphoplasmacytic infil- trates within myenteric ganglia and nerves. Similar infil- trates may also be present in the brain, spinal cord, peripheral nerves, conductive tissue of the heart, smooth and cardiac muscle, and adrenal glands. Non-suppurative leiomyositis and/or myocarditis may accompany the neu- ral lesions [6-9]. Clinically, PDD cases present with GI tract dysfunction (dysphagia, regurgitation, and passage of undigested food in feces), neurologic symptoms (e.g. ataxia, abnormal gait, proprioceptive defects), or both [3]. Although the clinical course of the disease can vary, it is generally fatal in untreated animals [3].
The cause of PDD is unknown, but several studies have raised the possibility that PDD may be caused by a viral pathogen. Evidence for an infectious etiology stems from the initial outbreaks of Macaw Wasting Disease, and other subsequent outbreaks of PDD [2,10]. Reports of pleomor- phic virus-like particles of variable size (30–250 nm) observed in tissues of PDD affected birds [8] led to the proposal that paramyxovirus (PMV) may cause the dis- ease; however, serological data has shown that PDD affected birds lack detectable antibodies against PMV of serotypes 1–4, 6, and 7, as well as against avian herpes viruses, polyomavirus, and avian encephalitis virus [3]. Similarly, a proposed role for equine encephalitis virus in PDD has been ruled out [11]. Enveloped virus-like parti- cles of approximately 80 nm in diameter derived from the feces of affected birds have been shown to produce cyto- pathic effect in monolayers of macaw embryonic cells [12], but to date no reports confirming these results or identifying this possible agent have been published. Like-wise, adeno-like viruses, enteroviruses, coronaviruses and reoviruses have also been sporadically documented in tis- sues or excretions of affected birds [3,13,14] yet in each case, follow-up evidence for reproducible isolation specif- ically from PDD cases or identification of these candidate agents has not been reported. Thus, the etiology of PDD has remained an open question.
To address this question, we have turned to a comprehen- sive, high throughput strategy to test for the presence of known or novel viruses in PDD affected birds. We employed the Virus chip, a DNA microarray containing representation of all viral taxonomy to interrogate 2 PDD case/control series independently collected on two differ- ent continents for the presence of viral pathogens. We report here the detection of a novel bornavirus signature in 62.5% of the PDD cases and none of the controls. These bornavirus-positive samples were confirmed by virus-spe- cific PCR testing, and the complete genome sequence has been recovered by ultra-high throughput sequencing combined with conventional PCR-based cloning.
Bornaviruses are a family of negative strand RNA viruses whose prototype member is Borna Disease Virus (BDV), an agent of encephalitis whose natural reservoir is prima- rily horses and sheep [15]. Although experimental trans- mission of BDV to many species (including chicks [16]) has been described, there is little information on natural avian infection, and existing BDV isolates are remarkable for their relative sequence homogeneity. The agent reported here, which we designate avian bornavirus (ABV) is highly diverged from all previously identified members of the Bornaviridae family and represents the first full- length bornavirus genome cloned directly from avian tis- sue. Subsequent PCR screening for similar ABVs con- firmed a detection rate of approximately 70% among PDD cases and none among the controls. Sequence anal- ysis of a single complete genome and all of the additional partial sequences that we have recovered directly from the PDD case specimens suggests that the viruses detected in cases of PDD form a new, genetically diverse clade of the Bornaviridae.
Results
Microarray-based detection of a Bornaviridae signature in PDD cases
To identify a possible viral cause of PDD, we applied the Virus chip, a DNA microarray containing 70 mer oligonu- cleotide probes representing all known viral sequences conserved at multiple nodes of the viral taxonomic tree [17,18] to identify viral signatures unique to histologi- cally confirmed cases of PDD. At the outset of this study, specimens from two independently collected PDD case/ control series were available for this investigation (Figure 1, Materials and Methods). The first series (n = 8), from
Figure 1
Clinical presentation of proventricular dilatation disease (PDD) cases and controls. A. Necropsy view of control (left panel) African gray parrot (Psittacus erithacus) that died of other causes. The normal-sized proventriculus is not visible in this view as it lies under the left liver lobe (L). Necropsy view of a great green macaw (Ara ambiguus) with PDD (right panel). The proventriculus (PV) is markedly distended and extends laterally well beyond the left lobe of L. The heart (H) is marked for orientation. B. Contrast fluoroscopy view of control (left panel) African gray parrot (Psittacus erithacus) 1.5 hours after admin- istration of barium sulfate. The kidney (K) is marked for orientation. The outline of both the PV and V is clearly visible, with normal size and shape. Within the intestinal loops (IL), wider and thinner sections represent active peristalsis. Right panel, rep- resentative PDD case, Eclectus parrot (Eclectus roratus) 18 hours after administration of barium. The PV is markedly distended and contains most of the contrast material, with less in the V and within the IL. A large filling defect (*) representing impacted food material. The kidney (K) is shown for orientation. These findings are typical for PDD; however PDD was not confirmed by histology in this case. C. Proventriculus histopathology. Hematoxylin and eosin staining of proventriculus histological sec- tions from a blue and yellow macaw (Ara ararauna) with PDD. Proventricular gland (G) is shown for orientation. Left panel, normal appearing myenteric ganglion detected within the proventriculus of this case (arrow); right panel, marked lymphoplas- macytic infiltration present within a myenteric ganglion (arrows). Right panel inset, higher magnification. D. CNS histopathol- ogy. Hematoxylin and eosin staining of a cerebral section from a control (left panel) African gray parrot (Psittacus erithacus) that died of other causes. Right panel, African gray parrot (Psittacus erithacus) with PDD. Perivascular cuffing is evident around blood vessels (arrows). Inset, higher magnification.
Figure 2
Avian bornavirus (ABV) genome sequence recovery and comparative analysis to Borna disease virus (BDV) genomes. A. Bornaviridae genome schematic. Grey bar at base, non-segmented negative sense viral RNA (vRNA) of Bornaviri- dae genome; coordinates of major sequence landmarks highlighted below. Green bars and dashed lines, transcription initiation sites (TISs); red bars, transcription termination sites. Distinct ORF-encoding transcription products and the gene products they encode are diagrammed above: TIS1 transcripts encoding nucleocapsid (N) gene, pink; TIS2 transcripts encoding phosphopro- tein (P) and X genes, green; TIS3 transcripts encoding the matrix (M), glycoprotein (G) and polymerase (large or 'L') gene, blue. Exons, thick solid black lines; introns, thin solid black lines; dashed black lines, 3'ends of transcripts generated transcription ter- mination read-through; shaded boxes, location of ORFs in transcripts; reading frames for ORFs from multiple genes generated from TIS3 indicated at right. Array probes track, Bornaviridae oligonucleotide 70 mer probes from the Virochip array. PCR primers track, primers generated for PCR follow up and screening of specimens in this study for detection of Bornaviridae spe- cies with expected product diagrammed below. vRNA RT-PCR track, overlapping vRNA clones and RACE products recovered directly from RNA extracted from crop tissue of a histologically confirmed case of PDD. Solexa reads track shows distribution of 33 mer reads with at least 15 bp sequence identity to recovered ABV genome sequence. Sequence identity with BDV genomes track shows scanning average pairwise nucleotide sequence identity (window size of 100 nucleotides, advanced in sin- gle nucleotide steps) shared between ABV and all BDV genome sequences in NCBI. A dashed line on the graph indicates 50% identity threshold for reference. B. Phylogenetic analysis of ABV genome and the 4 representative BDV genome isolates.
Neighbor-joining phylogenetic trees based on nucleotide sequences of the ABV genome sequence [GenBank:EU781967] and the following representative BDV genome sequences: H1766 [GenBank:AJ311523], V/Ref [GenBank:NC_001607], He/80 [GenBank:L27077], and No/98 [GenBank:AJ311524)] Scale bar, genetic distance.
samples originating in the United States, consisted of crop biopsy specimens from 3 histologically confirmed PDD cases and 5 controls that were provided for nucleic acid extraction and follow-up Virus chip analysis. The samples from the second series (n = 8) originated in Israel, where total RNA and DNA from proventriculus, ventriculus and brain specimens were extracted from 5 PDD cases and 3 controls. For each series, total RNA was reverse-tran- scribed with random primers, PCR-amplified, and fluo- rescently labeled and hybridized to the Virus chip microarray as previously described [18].
In these combined PDD case/control series, a Bornaviridae signature was detectable in 62.5% of the cases and none of the controls (Table 1). In the US cohort, which con- tained only GI tract specimens, we detected a bornavirus in 2 of 3 cases. Surprisingly, in samples from the Israeli PDD case/control series for which we had both GI tract and brain specimen RNA for each animal, we detected the Bornaviridae signature in 3 of the cases, but only in sam- ples derived from brain tissue. These signatures were unambiguously confirmed by follow-up PCR and sequence recovery, using primers based on the sequences of the most strongly annealing Bornaviridae oligonucle- otides on the microarray (Figure 2, Array probes and PCR probes tracks). These analyses revealed the presence of a set of surprisingly divergent avian bornaviruses (ABVs) in the PDD cases; the recovered sequences shared less than 70% sequence identity to any of the previously identified mammalian bornavirus isolates in the NCBI database.
Recovery of complete genome sequence of a divergent avian bornavirus (ABV) from a PDD case via ultra high- throughput sequencing and conventional RT-PCR
To determine if the sequence fragments we detected among specimens derived from PDD cases corresponded to the presence of a full-length bornavirus, we performed unbiased deep sequencing on a PCR-confirmed bornavi- rus positive PDD case that contained the highest concen- tration of RNA. To recover both mRNA and vRNA present in the sample, RNA from this specimen was linearly amplified with both oligo(dT) and random hexamer primers, and then PCR-amplified using a modified ran- dom amplification strategy compatible with the Solexa sequencing platform (Materials and Methods). An initial set of 1.4 million 33 mer reads was obtained from this template material. Filtering on read quality, insert pres- ence, and sequence complexity reduced this data set to 600,000 unique reads. Additional ELAND and iterative BLAST analyses ([19] Materials and Methods) of these reads against all avian sequences in NCBI (including ESTs, n = 918,511) identified reads in the dataset with at least 22 nucleotides of sequence identity likely derived from host transcripts randomly amplified during sequencing sample preparation. The 322,790 reads that passed this host filter were next screened for the presence of bornavi- rus sequence through similar ELAND and iterative BLAST analyses (Materials and Methods) using a database gener- ated from all Borna Disease virus (BDV) sequences present in NCBI (n = 207) and the sequences we had recovered from PCR follow-up of the PDD samples that tested positive for bornavirus by Virus chip microarray (n= 5). These analyses provided us with 1400 reads with at least a match of 15 or more nucleotides (blastn) or 7 or more predicted amino acids (tblastx) to known BDV sequences.
Mapping these 1400 reads onto their corresponding posi- tions on a consensus sequence for the 14 publicly availa- ble BDV genome sequences revealed spikes of high read coverage distributed discontinuously across the entire span of the BDV genome consensus. Reads containing blastn scores ? 90% identity to known BDV sequences were used as source sequences for primer design for PCR and sequence recovery of additional bornavirus sequence from both mRNA and vRNA templates present in the PDD specimen. Sequences recovered in this manner facilitated subsequent primer design for recovery of complete genome sequence via RT-PCR of 3 large overlapping frag- ments of the genome and 5'- and 3'-RACE (Figure 2A, vRNA RT-PCR track) directly from negative stranded vRNA present in the total RNA extracted from this clinical specimen.
As our initial PCR results suggested, the bornavirus genome sequence we recovered is quite diverged from all known BDV genomes, including the BDV isolate No/98, a divergent isolate sharing only 81% sequence identity with all other BDV genomes [20]. Overall, this newly recovered bornavirus genome sequence shares only 64% sequence identity at the nucleotide level to each of the complete BDV genomes. Scanning pairwise sequence identity anal- ysis indicates this genetic divergence exists across the entire genome (Figure 2A, Sequence identity shared with BDV genomes track). Given this divergence, we re-exam- ined the depth and distribution of the 322,790 reads from this specimen that passed the host filter to determine if we had missed reads derived from the recovered ABV in our initial screen against all BDV sequences. Not surprisingly, this retrospective BLAST analysis revealed an additional 2600 reads from across the recovered bornavirus genome that were missed in the initial BLAST analyses due to the lack of sequence conservation between the ABV sequence and the available BDV sequences (Figure 2A, Solexa reads track). In total, approximately 1% of all the high through- put shotgun reads could be mapped to the recovered bor- navirus genome.
Despite this sequence divergence, this avian bornavirus genomic sequence possesses all of the hallmarks of a Bornaviridae family member (Figure 2A): six distinct ORFs encoding homologs of the N, X, P, M, G, and L genes are detectable. Likewise, non-coding regulatory sequence ele- ments (the inverted terminal repeat sequences ([21], see Additional file 1: Alignment of bornavirus genomes 5' and 3' termini), the transcription initiation and termina- tion sites ([22], see Additional file 2: Alignment of tran- scription initiation and termination sites in bornavirus genomes), and each of the signals for pre-mRNA splicing ([23], see Additional file 3: Alignment of splice donor and acceptor sequences in bornavirus genomes) are all con- served in sequence and location, with the exception of the splice acceptor site 3 at position 4560 that has been previ- ously found in a subset, but not all BDV genomes [24,25]. Taken together, these data provide evidence that our anal- ysis has uncovered a novel divergent avian bornavirus (ABV) present in cases of PDD.
Phylogenetic and pairwise sequence analyses support this conclusion. Genomic and sub-genomic phylogenetic analyses of nucleotide sequences place the recovered ABV sequence on a branch distant from representative mem- bers of the 4 distinct genetic isolates of BDV for which complete genome sequences are available (Figure 2B, see Additional file 4: Phylogenetic relationships between sub- genomic loci of ABV and representative BDV genomes). Strikingly, the ABV genome sequence segregates to a posi- tion virtually equidistant from both the set of 3 closely related BDV isolates (V/Ref, H1766, and He/80) and the divergent No/98 BDV isolate (Figure 2B). Moreover, in contrast to the previously identified divergent No/98 iso- late, which retains a high level of conservation with other BDV isolates at the amino acid level, the ABV isolate also shows significant sequence divergence in the predicted amino acid sequence of every ORF in the genome (Table 2).
PCR screening of additional PDD cases and controls suggests an association between the presence of ABV and PDD
Recovery of the complete ABV genome sequence con- firmed that the microarray hybridization signature we detected accurately reflected the presence of bornaviruses
Table 1: ABV detection in PDD
a3 crop biopsies from US source and 5 brain and proventriculus/ ventriculus biopsies from Israel source were examined, with ABV detected in 2 of crop specimens and 3 brain specimens. b5 crop biopsies from US source and 3 brain and proventriculus/ventriculus biopsies from Israel source were examined.
in our PDD specimens. With these results in hand, we designed a set of PCR primers to perform ABV-specific PCR screening of an independent set of PDD case and control specimens to investigate the association between the presence of ABV and clinical signs and symptoms of PDD. An additional set of 21 samples derived from upper GI tract specimens (crop, proventriculus or ventriculus) from PDD cases and controls were screened for ABV sequences in a blinded fashion (Materials and Methods). For this analysis, we targeted three regions of the genome:
1) the L gene region of the genome that we used for PCR confirmation of the microarray results, (Figure 2, PCR probes track), 2) a subregion within the N gene and 3) a subregion within the M gene (Materials and Methods). PCR for glyceraldehyde 3 phosphate dehydrogenase (GAPDH) mRNA was performed in parallel with the ABV PCR on all specimens to control for integrity of RNA pro- vided from each specimen. Of the 21 specimens analyzed,
5 were positive for ABV by PCR and confirmed by sequence recovery. Unmasking the clinical status of these samples revealed that 7 of the samples were derived from confirmed PDD cases and 14 samples were derived from PDD controls. Among the PDD cases, we found 71% (5/
7) to be positive by ABV PCR (Table 3). In contrast, all PDD controls were negative by ABV PCR, and positive only for GAPDH mRNA. This PCR analysis provides an independent test of the statistical significance of the asso- ciation between the presence of ABV and histologically confirmed PDD (P = 0.01, Fisher's Exact Test). Although we do not observe ABV in 100% of PDD cases in this series (see Discussion), our results nonetheless indicate a signif- icant association of ABV with PDD.
Additional ABV isolates identified through PCR screening
Because we applied stringent inclusion criteria for the above-described association analysis study, a number of ABV (+) and ABV (-) samples were excluded. From these materials, six additional ABV isolates were detected – 5 derived from cases considered clinically suspicious and a sixth isolate derived from a confirmed PDD case for which only GI content and liver specimens were available. Addi- tional PCR screening of a set of 12 PDD control crop biopsy specimens provided to us unblinded again yielded solely ABV PCR (-) and GAPDH (+) results. These samples were excluded from the association analysis because we knew their clinical status prior to screening. We note that inclusion of these samples in statistical analyses would not diminish the association of ABV with known or suspected PDD.
Sequence analysis of ABV isolates indicates at least 5 divergent isolates in this branch of the Bornaviridae family
Recovery of partial sequence from additional isolates of ABV (from the above PDD case/control specimens as well as an additional samples derived from known or suspected PDD cases (Materials and Methods)) from 3 dis- tinct regions of the ABV genome provided the opportunity to further investigate the genetic diversity within this new branch of the Bornaviridae. Here, our description of results is restricted to comparison with representative members of the 4 major isolates of BDV, but virtually identical results were obtained when all available BDV sequences were analyzed.
As we observed for the complete ABV genome sequence, phylogenetic analysis of the recovered subgenomic ABV sequences revealed that each of the ABV isolates we recov- ered resides on a branch distant from the BDV isolates (Figure 3). PCR with the L gene consensus primers detected 14 isolates corresponding to 4 genetic subgroups of ABV. Each of these isolates were also detected with at least one of primer sets corresponding to the more highly expressed N gene and more conserved M gene regions of the genome; however, PCR with these two additional primer sets identified 2 additional ABV isolates that segre- gate to a genetically distinct 5th subgroup among the ABVs (ABV5, Figure 3B and 3C). Although these 5 distinct branches correlate largely according to the geographic ori- gin of the isolates, the genetic diversity we detect cannot be ascribed solely to differences in geographic origin of the isolates, since one of the branches (ABV4) is com- prised of isolates derived from both the U.S. and Israel. Likewise, we did not detect an obvious correlation between host species and genetic subgroup of ABV among the recovered isolates.
Pairwise sequence analyses of the nucleotide and pre- dicted amino acid sequence from the L region of the genome provide additional evidence for surprising genetic diversity among the ABV branches compared to that seen among the BDV branches (Table 4). Although derived from coding sequences of one of the more diver- gent genes of the bornavirus genome (Table 2, L gene), the region of the L gene we have used for PCR screening is rel- atively conserved among the BDV isolates, ranging from 81–98% at the nucleotide level, and 96–99% at the amino acid level (Table 4). In contrast, the sequence identity shared across this region of the genome among the ABV branches of the tree ranges from 77–83% at the nucle- otide level and 86–94% at the amino acid level. Taken together with the phylogenetic analysis described above, these data provide evidence that these ABV isolates form a new, genetically diverse branch of the Bornaviridae phyl- ogeny that is significantly diverged from the founder BDV isolates.
Discussion
It has been almost 40 years since the first description of PDD. Although a viral etiology has long been suspected, a convincing lead for a responsible viral pathogen has been lacking. By combining veterinary clinical investigation with genomics and molecular biology, we have identified a genetically diverse set of novel avian bornaviruses (ABVs) that are likely to play a significant role in this dis- ease. Through microarray analysis and follow-up PCR, we detected ABV sequences in 62.5% of the PDD cases in a set of specimens from two carefully collected PDD case/con- trol series originating from two different continents. We confirmed that these assays faithfully reflect the presence of full-length bornavirus in ABV PCR positive specimens through cloning of the complete ABV vRNA sequence directly from RNA extracted from one of these ABV PCR positive PDD case specimens. We next found evidence for a significant association between the presence of ABV and clinically confirmed PDD in follow-up blinded PCR screening of a set of additional PDD cases and controls, with ABV was detected in 71% of PDD cases and none of the controls (P = 0.01, Fisher's Exact Test).
Almost all prior sightings of bornaviruses in nature have been among mammals, and the mammalian isolates have been remarkably homogeneous at the sequence level (Table 2 and [15]). The latter is a surprising feature for RNA viruses, whose RNA-dependent RNA polymerases typically have high error rates. By contrast, the ABV iso- lates reported here are quite diverged from their mammalian counterparts, and show substantial heterogeneity among themselves. We note with interest that a single ear- lier report suggesting a potential avian reservoir for bornaviruses has been presented [26]. In that study, RT-PCR based on mammalian BDV sequences was used to recover partial sequences from stool collected near duck ponds where wild waterfowl congregate. However, the resulting sequences shared ca. 98% amino acid sequence homology to the mammalian BDVs, raising the possibility that these putative avian sequences might have resulted from possi- ble environmental or laboratory contamination [15]. Our ABV isolates, which are unequivocally of avian origin, are clearly very different from these sequences; it remains to be seen if other wild birds can indeed harbor BDV-like agents. The expanded sequence diversity of the bornavi- ruses discovered here should facilitate design of PCR primers that will enable expanded detection of diverse bornaviral types in future epidemiological studies.
Table 2: Predicted amino acid sequence similarity between ABV, the divergent BDV-No/98 and other BDV genomes
The known neurotropism of bornaviruses makes them attractive and biologically plausible candidate etiologic agents in PDD, since (i) PDD cases have well-described neurological symptoms such as ataxia, proproceptive defects and motor abnormalities; and (ii) the central GI tract pathology in the disorder results from inflammation and destruction of the myenteric ganglia that control per- istaltic activity. However, despite our success in ABV detection in PDD, we did not observe ABV in every PDD case analyzed. There are several possible explanations for this result. First, we do not know the tissue distribution (tropism) of ABV infection, or how viral copy number may vary at different sites as a function of the stage of the disease. By weighting our sample collection towards clin- ically overt PDD, we may have biased specimen accrual towards advanced disease. At this stage, where destruction of myenteric ganglial elements is often extensive, loss of infected cells may have contributed to detection difficul- ties (We note with interest in this context that in one of our case collections from Israel, virus detection occurred preferentially in CNS rather than in GI specimens). There are many precedents for such temporal variation in clini- cal virology – for example, in chronic hepatitis B viral loads typically decline by several orders of magnitude over the long natural history of the infection [27]. It is also pos- sible that our detection rate may merely reflect subopti- mal selection of PCR primers employed for screening; after all, our consensus primer selection was based on sequences we had recovered (L gene consensus primers) or sequence homology between the first fully sequenced
Table 3: Analysis of significance of ABV detection rate in PDD
ABV genome we recovered and a set of highly related mammalian BDV genome sequences (N and M gene con- sensus primers). We now recognize that there is substan- tial sequence variation within the ABVs (see Fig. 3); as more sequence diversity is recognized, better choices for more highly conserved primers will become apparent and could impact upon these prevalence estimates.
Finally, there could actually be multiple etiologic agents in PDD, with ABV infection accounting for only ~70% of the cases. Certainly both human and veterinary medicine are replete with examples of multiple agents that can trig- ger the same clinical syndrome – for example, at least 5 genetically unrelated viruses (hepatitis viruses A-E) are associated with acute hepatitis, and at least 3 of these can be implicated in chronic liver injury; similarly, several agents (RSV, rhinoviruses and occasionally influenza viruses) are implicated in bronchiolitis. To investigate this possibility, further high-throughput sequencing analysis of PDD cases that were negative for bornaviruses by PCR screening is currently underway.
Although ABV is clearly a leading candidate etiologic agent in PDD, formally establishing a causal role for ABV in PDD will require further experimentation. Such exper- iments could include (i) attempts to satisfy Koch's postu- lates via cultivation of ABV, followed by experimental transmission of infection and disease in inoculees, (ii) examination of seroprevalence rates in flocks with high and low PDD incidences, (iii) documentation of serocon- version accompanying development of PDD-like illnesses and (iv) examination of PDD cases by immunohisto- chemistry or in situ hybridization for evidence of colocal- ization of ABV infection at sites of histopathology. The recovery and characterization of a complete ABV genome and multiple isolates from this diverse new branch of the Bornaviridae family now opens the door to such investigations.
Conclusion
By combining clinical veterinary medical investigation with comprehensive pan-viral microarray and high throughput sequence analyses, we have identified a highly diverged set of avian bornaviruses directly from tissues of PDD cases, but not controls. These results are significant for a number of reasons. First, they provide a compelling lead in the long-standing search for a viral etiology of PDD, and pave the way for further investigations to assess the link between ABV and PDD. Second, these results also unambiguously demonstrate the existence of an avian res- ervoir of bornaviruses, expanding our understanding of the bornavirus host range. Finally, these results also pro- vide the first evidence that the Bornaviridae family is not confined to a set of genetically homogeneous species as was previously thought, but actually encompasses a set of heretofore unanticipated genetically diverse viral species.
Methods
PDD case and control definitions, specimen collection, and RNA extraction for pan-viral microarray screening
Two independent sets of PDD case and control specimens collected from two distinct geographic locations were independently prepared for pan-viral microarray screening and subsequent PCR screening. Sampling collection and inclusion criteria for each set are described below. Detailed information on each sample, along with results from histology, microarray, and PCR assays are provided in Additional file 5: Summary of clinical and molecular data for specimens provided in this study.
United States PDD case/control series
Specimen collection
All specimens provided for initial screening were crop tissue biopsies obtained from live psittacine birds to be used as normal controls or multiple tissue samples collected from clinically diseased birds at the time of euthanasia. Specimens were collected from client-owned birds from approximately August 2006 to May 2008 (All samples collected by S. Clubb). All of these samples originated from the southeast region of Florida. Crop biopsy tissue was collected from live birds under isoflurane anaesthesia. Following routine surgical preparation and sterile technique, the skin was incised over the center of the crop. The crop tissue was exposed and a section of tissue removed taking care to include large visible blood vessels. Fresh crop biopsy tissue was trimmed into tissue slices < 5 mm thick and submersed in RNAlater (Qiagen, Inc., USA, Valencia, CA) solution immediately upon harvest and frozen within 2 minutes of collection at -20°C to -80°C according to manufacturer's protocol, and held in this manner until shipped. A duplicate sample was fixed in 10% buffered formalin for routine histological examination with H & E stain. Time of frozen storage varied (2 weeks to 12 months) as samples were accumulated prior to shipping frozen. Clinically affected birds submitted as positives were euthanized under isoflurane anaesthesia and mixed tissues (proventriculus, ventriculus, heart, liver, spleen, kidneys, brain) were placed into RNAlater within 1 minute of death and frozen within 2 minutes of death. Duplicate samples were collected for histopathologic diagnosis of PDD.
Inclusion criteria
PDD-positive cases were required to meet the following criteria 1) Clinical history of characteristic wasting/malabsorption syndrome with dilation of the proventriculus and/or ventriculus and presence of undigested food in the stool and in most cases, a clinical history of ataxia or other CNS signs consistent with clinical PDD, and 2) histopathology confirming the presence of moderate to extensive lymphoplasmacytic ganglioneuritis affecting crop tissue and at least one of the following additional areas: proventriculus, ventriculus, brain, adrenal gland, or myocar- dium. PDD-negative controls were required to be from birds with no evidence of lymphoplasmacytic neurogan- gliitis on histopathology derived either from 1) normal birds with no clinical history of PDD or no known expo- sure to PDD or 2) birds which died of other causes. Crop biopsies from samples from living birds classified as suspicious cases were also submitted. Suspicious cases were defined histologically as having lymphocytes and plasma cells surrounding neurons but not infiltrating into the neurons. An additional specimen derived from a live bird raised with two necropsy-confirmed PDD birds in Virginia was also collected for analysis. Here, only cloacal swab and blood specimens were available and the lack of histopathological confirmation and crop tissue excluded this specimen from the ABV-PDD association analysis. How- ever, we did perform ABV PCR on these clinically suspicious specimens and the resulting viral sequences isolated were included in the subsequent comparative sequence analyses.
Table 4: Average pairwise sequence identity shared between ABV and BDV isolates*
RNA extractions
For RNA extractions, specimens were thawed in RNALater, sliced into 0.5 mm × 0.5 mm pieces, transferred to 2 ml of RNABee solution (Tel-Test, Inc., Friendswood, TX), homogenized with freeze thawing and scapel mincing, then extracted in the presence of chloroform according to manufacterer's instructions. Resulting RNA was next incu- bated with DNase (DNA-free, Applied Biosystems/ Ambion, Austin, TX) to remove any potential contaminat- ing DNA present in the specimen.
Israeli case/control series Specimen collection
Tissue samples were obtained from psittacine birds sub- mitted to the Division of Avian and Fish Diseases, Kimron Veterinary Institute (KVI) Bet Dagan, Israel, for diagnostic necropsy between July 2004 and March 2008. A few addi- tional specimens were obtained through private veterinar- ians. Some tissues were kept for nearly 4 years frozen either at -20°C or -80°C prior to testing, while others were fresh tissues from recent cases. The types of banked frozen tissue varied from case to case, while for some of the older cases only gastrointestinal content had been banked. Clin- ical histories for these birds were available from the sub- mission forms or through communication with the submitting veterinarians. The results of ancillary tests per- formed at the KVI were available through the KVI compu- terized records.
Inclusion criteria
Only cases for which appropriate histological sections were available for inspection were considered for this study. These had to include brain and at least two of the following tissues: crop, proventriculus, ventriculus. The tissue-types examined for each bird for which specimens were provided are listed in Data File S1. PDD-positive cases were required to have evidence of lymphoplasma- cytic infiltration of myenteric nerves and/or ganglia within one or more of the upper GI tract tissues men- tioned above. These were all derived from birds that had been suspected to have PDD based on their clinical case histories and/or necropsy findings. PDD-negative con- trols had no detectable lesions and no evidence of non- suppurative encephalitis. For most birds in the PDD-neg- ative group, a cause of death (other than PDD) has been determined. Two birds that came from a known PDD out- break, but showed only cerebral lymphoplasmacytic perivascular cuffing, were classified as 'suspicious'. These were excluded from the statistical analysis, as were all other birds for which a PDD status could not be clearly determined and classified as 'inconclusive' (e.g. due to poor tissue preservation, poor section quality, or scarcity of myenteric nerves within the tissues examined).
RNA extraction
When possible, a sample of brain as well as a combined proventricular/ventricular sample was prepared for RNA extraction for each bird. If not available, other tissues and/ or gastrointestinal content were used (see Additional file 5: Summary of clinical and molecular data for specimens provided in this study). Frozen samples were allowed to thaw for 1–2 hours at room temperature prior to han- dling. Then, under a laminar flow biohazard hood and using aseptic technique, approximately 1 cm3 of each tis- sue was macerated by two passages through a 2.5 ml ster- ile syringe and transferred into sterile test tubes containing 4 ml nuclease-free PBS. The content of the tubes was mixed by vortex for 30 sec, and the tubes were placed overnight at 4°C. RNA extraction was performed on the following day, using either the QIAamp® viral RNA kit (Qiagen, Valencia, CA; batch1&2, specimens 1–8) or the TRI Reagent® kit (Molecular Research Center, Cincin- nati, OH; all other specimens), following the manufactur- ers' instructions. The end product was either provided lyophilized (batches 1 and 2, samples 1–9) as a dry pellet, or re-suspended in 40 ul nuclease-free water.
Virus chip hybridization experiments
Microarray analysis of specimens was carried out as previ- ously described [18]. Briefly, 50–200 ng of DNAse-treated total RNA from each sample was amplified and labelled using a random-primed amplification protocol and hybridized to the Virochip. Microarrays (NCBI GEO plat- form GPL3429) were scanned with an Axon 4000B scan- ner (Axon Instruments). Virochip results were analyzed using E-Predict [28] and vTaxi (K. Fischer et al., in prepa- ration).
PCR primers for detection of avian bornaviruses
Microarray-based Bornaviridae PCR primers
Initial PCR primers were generated based on two of the 70 mer microarray probes with hybridization signal in the Bornaviridae positive arrays that localize to positions 3676–3745 and 4201–4270 of the Bornaviridae reference sequence [GenBank:NC_001607]. Subsequences within each of these probes (BDV_LconsensusF: 5'- CCTCGCGAGGAGGAGACGCCTC-3' and BDV_LconsensusR: 5' CTGCTCTTGGCTGTGTCTGCTGC-
3'; positions 3710–3729 and 4252–4230, respectively of the NCBI Bornaviridae reference sequence) that are 100% conserved across the 12 other fully sequenced bornavirus genome isolates in NCBI (huP2br [GenBank:AB258389], Bo/04w [GenBank:AB246670], No/98 [Gen- Bank:AJ311524], H1766 [GenBank:AJ311523], He/80/FR [GenBank:AJ311522], V/FR [GenBank:AJ311522], virus rescue plasmid pBRT7-HrBDVc [GenBank:AY05791], CRNP5 [GenBank:AY114163], CRP3B [Gen- Bank:AY114162], CRP3A [GenBank:AY114161], He/80
[GenBank:L27077], and V [GenBank:U040608]) were utilized for initial follow-up PCR and sequence confirmation of microarray screening results. Briefly, 1 ul of the ran- domly amplified nucleic acid prepared for microarray hybridization from all specimens was utilized as template for 35 cycles of PCR, under the following conditions: 94°C, 30 seconds; 50°C, 30 seconds; 72°C, 30 seconds. Resulting PCR products were gel purified, subcloned into the TOPO TA cloning vector pCR2.1 (Invitrogen, USA, Carlsbad CA) and sequenced with M13F and M13R prim- ers.
Generation of ABV consensus PCR primers
Sequences recovered from BDV_LconsensusF and BDV_LconsensusR PCR products were aligned, and an additional set of ABV consensus primers biased towards the ABV sequences were identified: ABV_LconsensusF, 5'- CGCCTCGGAAGGTGGTCGG-3' (maps to positions aligning with residues 3724–3742 of BDV reference genome) and ABV_LconsensusR, 5'-GGCAYCAYCK- ACTCTTRAYYGTRTCAGC-3' (maps to positions aligning with residues 4233–4257 of BDV reference genome). Using identical PCR cycling conditions as described above for the microarray-based Bornaviridae PCR assay, these ABV consensus primers were found to be > 100X more sensitive for ABV detection compared to BDV_LconsensusF and BDV_LconsensusR primers, and were thus utilized to re-screen the initial set of PDD case and control samples provided for microarray analysis (no additional positives identified) and all subsequently pro- vided samples. Two additional PCR primers in the N (ABV_NconsensusF: 5'-CCHCATGAGGCTATWGATT-
GGATTAACG-3' and ABV_NconsensusR: 5'-GCMCGG- TAGCCNGCCATTGTDGG-3') and M (ABV_MconsensusF: 5'-GGRCAAGGTAATYGTYCCT-
GGATGGCC-3' and ABV_PconsensusR: 5'-CCAACAC- CAATGTTCCGAAGMCG-3') that mapped to conserved sequences shared between the complete ABV genome sequence and the 14 other fully sequenced BDV genomes in the NCBI database were also employed for PCR screen- ing of PDD cases and controls.
Ultra high-throughput sequencing
Sample preparation and sequencing
500 ng of total RNA derived from one of the PDD case specimens was linearly amplified via modification of the MesssageAmp aRNA kit (Applied Biosystems/Ambion, Austin, TX). To ensure the amplification of both mRNA and vRNA present in the specimen, T7-tailed random nonamer was mixed in an equimolar ratio with the man- ufacturer-provided T7-oligo(dT) primer during the 1st strand synthesis step. The resulting aRNA was next used as input for a modified version of Genomic DNA sample preparation protocol for ultra high-throughput Solexa sequencing (Illumina, Hayward, CA). 400 ng of the input aRNA was reverse-transcribed with reverse transcriptase (Clontech Laboratories, Inc., Mountain View, CA) using a random nonamer tailed with 19 bp of the Solexa Long (5'- CACGACGCTCTTCCGATCTNNNNNNNNN-3') primer sequence (Illumina, Hayward CA). Following termination of reaction, first strand cDNA products were purified from the reaction with Qiagen MinElute spin column (Qiagen USA, Valencia CA). To ensure stringent separation from primers, the MinElute eluate was then filtered through a Microcon YM30 centrifugal filter (Millipore Corp., Biller- ica, MA). The resulting eluate served as template for 2nd strand synthesis in a standard Sequenase 2.0 (USB, Cleve- land, OH) reaction primed with a random nonamer tailed with 22 bp (5'-GGCATACGA GCTCTTC-CGATCTNNNNNNNNN-3') of the Solexa Short primer sequence (Illumina, Hayward CA). Double-stranded DNA products were separated from primers and very short products through a second Qiagen MinElute spin column run followed by a Microcon YM50 centrifugal filter. This eluate was used as template for 10 cycles of PCR amplifi- cation with the full length Solexa L and S primers using KlenTaq LA DNA polymerase mix (Sigma-Aldrich, St. Louis, MO). PCR product was purified from the reaction with a MinElute spin column. Following cluster genera- tion, Solexa sequencing primer was annealed to the flow cell, and 36 cycles of single base pair extensions were per- formed with image capture using a 1G Genome Analyzer (Illumina, Hayward, CA). The Solexa Pipeline software suite version 0.2.2.6 (Illumina, Hayward, CA) was utilized for base calling from these images. Using software default quality filters, cycles 4–36 were deemed high quality, resulting in a total of 1.4 million 33 mer reads for down- stream sequence analyses.
Identification of Bornaviridae reads
Reads sharing 100% identity to each other or the Solexa amplification primers were filtered, reducing our initial set of 1.4 million reads to a working set of 600,000 unique reads. In order to quickly assess the homology of this set of reads to different sequence databases, we employed an iterative strategy using ELAND (Efficient Local Alignment of Nucleotide Data) and BLAST analyses. To filter reads from our analysis potentially derived from psittacine host tissue, the working set of reads were aligned to a database of all Aves sequences from NCBI (n = 918,511) using ELAND, which tolerates no more than 2 base mismatches, and discards both low quality reads and reads with low sequence complexity. Reads that did not align to the Aves database by ELAND analysis were next re-aligned to the Aves database for high stringency blastn analysis (e = 10-7, word size = 11), followed by progressively lower stringen- cies (down to e = 10-2, word size = 8), corresponding to reads containing only 22 nucleotide identities to sequences in the Aves database. To identify reads with some homology to Bornaviridae sequences in the resulting set of 322,790 host-filtered reads, we re-implemented the ELAND/iterative blastn analysis strategy (down to ? 15 nucleotides identity) using a database of all NCBI BDV sequences (n = 207) augmented by our previously recov- ered ABV sequences (n = 5). An additional iterative tblastx analysis was incorporated to capture distantly related reads that shared similarity to the known BDV sequences only at the level of predicted amino acid sequence (down to ? 6 amino acid identity).
Complete ABV vRNA genome sequence recovery by RT- PCR
Initial genome sequence recovery
Sequences from 33 mer reads from the deep sequencing with a minimum of 91% sequence identity with known BDV sequences present in the NCBI database were utilized to generate a set of primers for additional cloning and sequence recovery by RT-PCR of both mRNA and vRNA present in the clinical specimen. In this manner, we gen- erated a hybrid assembly derived from multiple overlap- ping clones and 5' RACE products encompassing the ABV genome sequence.
vRNA genome sequence recovery
To ensure recovery of accurate sequence across the ABV genome, especially at splice junctions and transcription initiation and termination sites, we utilized the sequence from ABV hybrid assembly to design primers for recovery of 3 overlapping products by RT-PCR directed against the vRNA present in the specimen. Aliquots of 500 ng of DNAse-treated total RNA extracted from the clinical spec- imen were annealed with 3 primers complementary to the predicted vRNA sequences: ABV1r, 5'-ATGACCAGGAC- GAGGAGATG-3' (maps to residues 8831-8812 of vRNA), ABV2r, 5'-CCTGTGAATGTCTCGTTTCTG-3' (maps to resi- dues 5754-5733 of vRNA), and ABV3r 5-TTCTTTCAG-
CAACCACTGACG-3' (maps to residues 2563-2543 of vRNA). Reverse transcription was carried out at 50°C for 1 hr with SuperScriptIII (Invitrogen, Carlsbad CA) accord- ing to manufacturer's instructions. Following RNase H treatment, PCR was performed on the resulting cDNA with Phusion polymerase (NEB, Ipswich, MA) with the primers used for reverse transcription and the following primers: ABV1f: 5'-GGATCATTCCTTGATGATGTATTAGC-
3', (maps to residues 5567-5589) ABV2f: 5'-CAAATGGA- GAGCCTGATTGG-3' (maps to residues 2378-2397) ABV3f: 5'-AATCGGTAAGTCCAGAGTCAAGG-3' (maps to residues 155-177). All products were amplified for 35 cycles under the following conditions: 98°C, 3 minutes; 98°C, 10 seconds, 50°C, 30 seconds, 72°C 3 minutes. Resulting products were gel purified, and subcloned into the TOPO T/A cloning vector pCR2.1 after incubation with Taq polymerase and dATP for 10 minutes at 72°C. For each product, 4 independent transformants were pre- pared for standard dideoxy sequencing on an ABI3730 sequencer (ElimBio, Hayward CA). Forward and reverse reads spanning each clone were generated using M13F and M13R and additional overlapping primers spaced at 600–800 bp intervals across the each of the clones.
5' and 3' RACE to sequence at vRNA termini
vRNA RT-PCR products containing uncapped vRNA ter- mini were captured using the First Choice RLM RACE kit (Ambion, Austin TX) with the following modifications to the standard protocol: 1) tobacco acid phosphotase treat- ment was omitted, 2) a phosphorylated RNA, RNAligate, 5'-p-GUUAUCACUUUCACCC-3' (gift of J. Shock, DeRisi lab) was substituted for the 3' RNA ligation-mediated RACE primer provided in the kit and ligated to 3' ends as per manufacterer's 5' RACE protocol, and 3) in the 3' RACE reverse transcription reactions, two reverse tran- scription reactions were performed and carried forward in parallel: one with random decamers and one with a DNA oligo complementary to oJSmer utilized in the RNA liga- tion step (ligateRC, 5'-p-GGGTGAAAGTGATAAC-3'). For 5' RACE, a single round of PCR was sufficient to generate a product using the vRNA specific primer ABV5RaceOuter, 5'-CAGTCGGTTCTTGGACTTGAAGTATCTAGG-3' (maps to residues 346-317 of vRNA) and manufacturer provided outer PCR primer. For 3' RACE, nested PCR was required to recover detectable PCR product of expected size using outer PCR primers oJSmerRC and the gene specific primer ABV3RaceOuter, 5'-CCCGTCTACTGTTCTTTCGCCG-3' (maps to residues 8479-8497 of vRNA), followed by inner PCR using Tailed_RNAligateRC, 5'- AAGCAGTGGTAACAACGCAGAGTACGGGTGAAAGT- GATAAC-3' and the gene specific primer, ABV3RaceInner, 5'-GCAATCCAGGAATAAGCAAGCACAAA-3' (maps to residues 8595-8620 of vRNA). Both of the RACE PCR reactions were carried out with Platinum Taq polymerase (Invitrogen, Carlsbad, CA) in 35 cycles of gradient PCR (with varying annealing temperature): 94°C, 30 seconds; 55–58°C, 30 seconds; 72°C, 30 seconds. Resulting PCR products were gel purified and subcloned into TOPO T/A cloning vector pCR 4.0. For the 5' RACE products, 7 inde- pendent transformants from 3 independently generated PCR products were subcloned and sequenced with M13F and M13R primers. For the 3' RACE products, 6 independ- ent transformants from 4 independently generated PCR products were subcloned and sequenced with M13F and M13R primers. Terminal sequences reported here reflect the majority consensus sequence obtained from these reads.
Genome sequence assembly
Genome sequence assemblies from both initial genome sequence recovery and vRNA genome sequence recovery were generated using Consed, version 16.0 software [29]. All bases from the resulting vRNA genome sequence assembly are covered at least 4× with a minimum Phred value of 20.
Blinded PCR screening of additional PDD cases and controls
Beyond the initial set of 16 specimens provided for micro- array analysis, specimens from a total of 38 additional PDD cases, PDD controls, and PDD suspicious birds with varied clinical histories were provided to us blinded by our 2 collaborators (see Additional file 5: Summary of clinical and molecular data for specimens provided in this study).
Sample processing
For specimens provided in tissue form from the US collab- orators, total RNA was extracted as described above with RNABee, DNase treated, then reverse-transcribed and PCR-amplified according to our random amplification protocol for microarray sample preparation (Materials and Methods). Specimens provided from Israel in the form of extracted RNA were similarly DNAse-treated and amplified prior to PCR screening.
PCR screening
1 ul of the randomly amplified material generated from these RNA samples was used as input template for ABV consensus PCRs as described above. In parallel, as an independent control for input specimen RNA integrity, PCR for glyceraldehyde 3-phosphate dehydrogenase (GAPDH) mRNA was performed on all specimens using designed based on Friedman-Einat et al [30] and Gallu gal- lus GAPDH sequence: Gg_GAPDHf: 5'-AGTCATCCCT- GAGCTSAAYGG*GAAGC-3' (bp708-733 in Gallus gallus cDNA (NCBI accession NM_204305), * indicates the junction of GAPDH exon 8 and 9 spanned by this primer), Gg_GAPDHr 5'-ACCATCAAGTC-
CACAACACGG-3' (Spans bp 1037-1017 in Gallus gallus GAPDH cDNA (NCBI accession NM_204305), maps to GAPDH exon 12). After PCR results were tallied, clinical information on all specimens tested was unmasked. A complete accounting of ABV, GAPDH PCR results, speci- men type and clinical status is provided in Additional file 5: Summary of clinical and molecular data for specimens provided in this study.
Sample inclusion for association analysis
To reduce potential confounding due to differences in viral detection resulting from specimen tissue source, only specimens derived from upper GI tract tissue (crop, prov- entriculus/ventriculus) that tested positive by GAPDH mRNA PCR were included in association analysis pre- sented in Table 3. This consisted of a total of 21 speci- mens, 7 of which were derived from histologically confirmed PDD cases and 14 derived from histologically negative control specimens.
Samples excluded from association analysis
The remaining 17 samples were excluded from the analy- sis because they were either 1) GAPDH-positive or GAPDH-negative samples derived from specimen other than upper GI tract tissue (GI content, brain, liver, or intestine) or 2) derived from cases that were histologically or clinically 'suspicious', but unconfirmed PDD cases. Six additional ABV PCR positives were identified among this set of samples excluded from the statistical analyses: 1 derived from GI content from a confirmed PDD case, and 5 derived from a variety of tissues from the PDD suspi- cious cases.
Phylogenetic and comparative sequence analysis
Multiple sequence alignments of complete genome sequences or partial sequences derived from PCR screen- ing studies were generated with ClustalW [31] version
1.83. Resulting alignments were used for scanning pair- wise sequence analysis (window size, 100; step size 1 nucleotide steps). Additional ClustalW alignments and neighbor-joining phylogenetic trees were generated using Mega software, version 4.0.2 [32].
List of abbreviations
ABV: Avian bornavirus; BDV: Borna diseae virus; PDD: Proventricular dilatation disease.
Competing interests
Sequence information obtained here has been disclosed for patenting purposes. ALK, AG, SC, PS-C, KF, KS, CYC, AL, AG, SKB, DG, and JLD were all party to this disclosure in conjunction with UCSF Office of Technology Manage- ment.
Authors' contributions
ALK participated in the conception, design, and coordina- tion of the study, performed specimen extraction of spec- imens from Florida case/control study, array analyses for both sets of PDD case/control series, follow-up PCR screening and sequencing of samples and wrote the man- uscript, AG orchestrated and collected the PDD case/con- trol specimens from Israel and coordinated the clinical and histopathology analyses, and nucleic acid extraction for samples from Israel, and participated in revising the manuscript, SC orchestrated and collected the Florida PDD case/control specimens and oversaw the clinical and histopathologic analyses of these samples from Florida, and participated in revising the manuscript, PS-C carried out filtering and iterative BLAST analysis of ultra high throughput sequence data for ABV genome sequence recovery, participated in primer design and complete genome sequence recovery, and drafting the manuscript, KF participated in array analysis, developed pipeline for ultra high throughput sequence analysis, and participated in design of filtering and iterative BLAST analysis, KS per-
formed modified library preparation for ultra high throughput sequencing and participated in revising the manuscript, CYC performed ultra high throughput sequencing and participated in revising the manuscript, AL, SM, and YF participated in clinical evaluation, speci- men collection and extraction of samples from Israel, AG participated in extraction of specimens from Florida and follow-up microarray analysis and high throughput sequencing, CCW developed additional primers for PCR follow-up studies, SBK assisted in the selection of the PDD case/control specimens from Florida and partici- pated in review of clinical and histological status of cases and controls included in the study, DG and JLD oversaw the overall conception and design of the project and supervised all phases of its execution and the drafting and revision of the manuscript.
Additional material
Additional file 1
Alignment of bornavirus genomes 5' and 3' termini. Bornavirus genome organization overview diagrammed as in Figure 2. Sequences in alignments shown are complementary to vRNA sequence, genome isolate names shown at left. 3' end sequence recovered for ABV genome and other BDV genomes is shown in left panel, 5' end sequence recovered for ABV genome and other BDV genomes is shown in right panel. Accession num- bers for genomes aligned: hu2Pbr [GenBank:AB258389], Bo/04w [Gen- Bank:AB246670], H1766 [GenBank:AJ311523], Ref [GenBank:NC_001607], V [GenBank:U04608], V/FR [Gen- Bank:AJ311521], CRNP5 [GenBank:AY114163], CRP3B [Gen- Bank:AY114162], CRP3A [GenBank:AY114161], He/80/FR [GenBank:AJ311522], He/80 [GenBank:L27077], pBRT7-HrBDVc [GenBank:AY705791], No/98 [GenBank:AJ311524], ABV [Gen- Bank:EU781967].
Click here for file [http://www.biomedcentral.com/content/supplementary/1743-422X-5-88-S1.pdf422X-5-88-S1.pdf]
Additional file 2
Alignment of transcription initiation and termination sites in borna- virus genomes. Panel A, alignment of the 3 bornavirus transcription ini- tiation sites (TIS) and 6 nucleotides of flanking sequences. Panel B, alignment of the 4 bornavirus transcription termination sites. Source genomes for alignments are shown at left. Black trianges highlight ABV sequences.
Click here for file [http://www.biomedcentral.com/content/supplementary/1743-422X-5-88-S2.pdf 422X-5-88-S2.pdf]
Additional file 3
Alignment of splice donor and splice acceptor sequences in bornavirus genomes. Panel A, alignment of splice donor 1 and splice acceptor 1 sequences; Panel B, alignment of splice donor 2 and splice acceptor 2 sequences; Panel C, alignment of splice acceptor 3 sequences. Source genomes for alignments are shown at left.
Click here for file [http://www.biomedcentral.com/content/supplementary/1743-422X-5-88-S3.pdf 422X-5-88-S3.pdf]
Additional file 4
Phylogenetic relationships between sub-genomic loci of ABV and rep- resentative BDV genomes. Neighbor-joining trees generated for the indi- cated nucleotide sequences of ABV and a representative set of BDV genomes are shown for each ORF in the bornavirus genome. Accession numbers of representative BDV genomes are: Ref/V [Gen- Bank:NC_001607], H1766 [GenBank:AJ311523], He/80 [Gen- Bank:AY705791], No/98 [GenBank:AJ311524].
Click here for file [http://www.biomedcentral.com/content/supplementary/1743-422X-5-88-S4.pdf 422X-5-88-S4.pdf]
Additional file 5
Summary of clinical and molecular data for specimens provided in this study. Microsoft Excel file containing two spreadsheet (US specimens and Israel specimens) summarizing clinical and epidemiologic information available for each specimen, as well as the associated results from the described microarray/PCR/sequence experiments.
Click here for file [http://www.biomedcentral.com/content/supplementary/1743-422X-5-88-S5.xls422X-5-88-S5.xls]
Acknowledgements
Jenny Shock (DeRisi lab, UCSF) for providing RNA oligos for 3' RACE exper- iments; Prof. Shmuel Perl, head of the Division of Pathology (KVI), for allowing us access to the KVI histopathology specimen collection; Dr. Asaf Berkovich.
(KVI) for assistance with specimen preparation and retrieval, Dr. Uri Bend- heim, Dr. Revital Harari, and Dr. Anthony Poutous for submitting case mate- rial from their practices; and the Lahser Interspecies Research Foundation for providing funding for US specimen collection and veterinary care. The remain- der of this work was supported by HHMI grants to JLD and DG.
References
1. Daoust PY, Julian RJ, Yason CV, Artsob H: Proventricular impac- tion associated with nonsuppurative encephalomyelitis and ganglioneuritis in two Canada geese. J Wildl Dis 1991, 27:513-517.
2. Doneley RJ, Miller RI, Fanning TE: Proventricular dilatation dis- ease: an emerging exotic disease of parrots in Australia. Aust Vet J 2007, 85:119-123.
3. Gregory C, Latimer KS, Niagro F, Ritchie BW, Campagnoli RP, Nor- ton TM, Greenacre CB: A review of proventricular dilatation syndrome. J Assoc Avian Vet 1994, 8:69-75.
4. Perpinan D, Fernandez-Bellon H, Lopez C, Ramis A: Lymphoplas- macytic myenteric, subepicardial, and pulmonary ganglione- uritis in four nonpsittacine birds. J Avian Med Surg 2007, 21:210-214.
5. Sullivan ND, Mackie JT, Miller RI, Giles A: First case of psittacine proventricular dilatation syndrome (macaw wasting disease) in Australia. Aust Vet J 1997, 75:674.
6. Berhane YSD, Newman S, Taylor M, Nagy E, Binnington B, Hunter B: Peripheral neuritis in psittacine birds with proventricular dil- atation disease. Avian Pathol 2001, 30:563-570.
7. Lutz ME, Wilson RB: Psittacine proventricular dilatation syn- drome in an umbrella cockatoo. J Am Vet Med Assoc 1991, 198:1962-1964.
8. Mannl A, Gerlach H, Leipold R: Neuropathic gastric dilatation in psittaciformes. Avian Dis 1987, 31:214-221.
9. Vice CA: Myocarditis as a component of psittacine proven- tricular dilatation syndrome in a Patagonian conure. Avian Dis 1992, 36:1117-1119.
10. Lublin A, Mechani S, Farnoushi I, Perl S, Bendheim U: An outbreak of proventricular dilatation disease in psittacine breeding farm in Israel. Israel Journal of Veterinary Medicine 2006, 61:16-19 [http://www.isrvma.org/journal.htm].
11. Gregory CR, Latimer KS, Niagro FD, Roberts AW, Campagnoli RP, Pesti DA, Ritchie BW, Lukert PD: Investigations of Eastern Equine Encephalomyelitis Virus as the Causative Agent of Psittacine Proventricular Dilatation Syndrome. J Avian Medi- cine and Surgery 1997, 11:187-193.
12. Gough RE, Drury SE, Harcourt-Brown NH, Higgins RJ: Virus-like particles associated with macaw wasting disease. Vet Rec 1996, 139:24.
13. Gough RE, Drury SE, Culver F, Britton P, Cavanagh D: Isolation of a coronavirus from a green-cheeked Amazon parrot (Ama- zon viridigenalis Cassin). Avian Pathol 2006, 35:122-126.
14. Ritchie B: Avian Viruses: Function and Control Lake Worth: Wingers Publishing; 1995.
15. Durrwald R, Kolodziejek J, Muluneh A, Herzog S, Nowotny N: Epi- demiological pattern of classical Borna disease and regional genetic clustering of Borna disease viruses point towards the existence of to-date unknown endemic reservoir host popu- lations. Microbes Infect 2006, 8:917-929.
16. Rott R, Becht H: Natural and experimental Borna disease in animals. Curr Top Microbiol Immunol 1995, 190:17-30.
17. Chiu CY, Alizadeh AA, Rouskin S, Merker JD, Yeh E, Yagi S, Schnurr D, Patterson BK, Ganem D, DeRisi JL: Diagnosis of a critical res- piratory illness caused by human metapneumovirus by use of a pan-virus microarray. J Clin Microbiol 2007, 45:2340-2343.
18.Chiu CY, Rouskin S, Koshy A, Urisman A, Fischer K, Yagi S, Schnurr D, Eckburg PB, Tompkins LS, Blackburn BG, et al.: Microarray detection of human parainfluenzavirus 4 infection associated with respiratory failure in an immunocompetent adult. Clin Infect Dis 2006, 43:e71-76.
19. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lip- man DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25:3389-3402.
20.Nowotny N, Kolodziejek J, Jehle CO, Suchy A, Staeheli P, Schwemmle M: Isolation and characterization of a new subtype of Borna disease virus. J Virol 2000, 74:5655-5658.
21. Schneider U, Schwemmle M, Staeheli P: Genome trimming: a unique strategy for replication control employed by Borna disease virus. Proc Natl Acad Sci USA 2005, 102:3441-3446.
22. Schneemann A, Schneider PA, Kim S, Lipkin WI: Identification of signal sequences that control transcription of borna disease virus, a nonsegmented, negative-strand RNA virus. J Virol 1994, 68:6514-6522.
23. Schneider PA, Schneemann A, Lipkin WI: RNA splicing in Borna disease virus, a nonsegmented, negative-strand RNA virus. J Virol 1994, 68:5007-5012.
24. Cubitt B, Ly C, de la Torre JC: Identification and characteriza- tion of a new intron in Borna disease virus. J Gen Virol 2001, 82:641-646.
25. Tomonaga K, Kobayashi T, Lee BJ, Watanabe M, Kamitani W, Ikuta K: Identification of alternative splicing and negative splicing activity of a nonsegmented negative-strand RNA virus, Borna disease virus. Proc Natl Acad Sci USA 2000, 97:12788-12793.
26. Berg MJM, Montell H, Berg A-L: Wild birds as a possible natural reservoir of Borna Disease Virus. Epidemiology and Infection 2001, 127:173-178.
27. Ganem D, Prince AM: Hepatitis B virus infection – natural his- tory and clinical consequences. N Engl J Med 2004, 350:1118-1129.
28. Urisman A, Fischer KF, Chiu CY, Kistler AL, Beck S, Wang D, DeRisi JL: E-Predict: a computational strategy for species identifica- tion based on observed DNA microarray hybridization patterns. Genome Biol 2005, 6:R78.
29. Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res 1998, 8:195-202.
30. Friedman-Einat M, Boswell T, Horev G, Girishvarma G, Dunn IC, Tal- bot RT, Sharp PJ: The chicken leptin gene: has it been cloned? Gen Comp Endocrinol 1999, 115:354-363.
31. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994, 22:4673-4680.
32. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolu- tionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 24:1596-1599.