4
views
0
recommends
+1 Recommend
0 collections
    0
    shares
      • Record: found
      • Abstract: found
      • Article: found
      Is Open Access

      First evidence of a new simian adenovirus clustering with Human mastadenovirus F viruses

      research-article

      Read this article at

      Bookmark
          There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

          Abstract

          Background

          Adenoviruses play an important role as human pathogens, though most infections are believed to be asymptomatic. The over 100 human adenovirus types are classified into seven species (A-G), some of which include simian adenoviruses. Recent findings have highlighted that simian adenoviruses have a zoonotic potential and that some human adenoviruses are likely the result of relatively recent spillover events.

          Methods

          In order to evaluate the risks associated with primates hunted and sold as bushmeat, multiple samples from 24 freshly killed monkeys were collected in the Republic of the Congo and tested for adenovirus DNA by PCRs targeting the conserved DNA polymerase and hexon genes.

          Results

          The DNA of a novel simian adenovirus was detected in a moustached monkey ( Cercopithecus cephus) by the DNA polymerase PCR, but not by the hexon PCR. The 275 nucleotide amplicon was most closely related to members of the Human mastadenovirus F species (93% HAdV-40 and 89% HAdV-41 amino acid identity), rather than to other known simian adenoviruses.

          Conclusions

          The phylogenetic clustering with Human mastadenovirus F sequences suggests a common ancestor, more recent than the last common ancestor of humans and moustached monkeys. The findings increase concerns about the zoonotic potential of simian adenoviruses and highlight the need for more research and surveillance on the issue.

          Related collections

          Most cited references18

          • Record: found
          • Abstract: found
          • Article: not found

          Phylogeny and geography predict pathogen community similarity in wild primates and humans.

          In natural systems, host species are often co-infected by multiple pathogen species, and recent work has suggested that many pathogens can infect a wide range of host species. An important question therefore is what determines the host range of a pathogen and the community of pathogens found within a given host species. Using primates as a model, we show that infectious diseases are more often shared between species that are closely related and inhabit the same geographical region. We find that host relatedness is the best overall predictor of whether two host species share the same pathogens. A higher frequency of pathogen host shifts between close relatives or inheritance of pathogens from a common ancestor may explain this result. For viruses, geographical overlap among neighbouring primate hosts is more important in determining host range. We suggest this is because rapid evolution within viral lineages allows host jumps across larger evolutionary distances. We also show that the phylogenetic pattern of pathogen sharing with humans is the same as that between wild primates. For humans, this means we share a higher proportion of pathogens with the great apes, including chimpanzees and gorillas, because these species are our closest relatives.
            Bookmark
            • Record: found
            • Abstract: found
            • Article: not found

            New adenovirus species found in a patient presenting with gastroenteritis.

            An unidentified agent was cultured in primary monkey cells at the Los Angeles County Public Health Department from each of five stool specimens submitted from an outbreak of gastroenteritis. Electron microscopy and an adenovirus-specific monoclonal antibody confirmed this agent to be an adenovirus. Since viral titers were too low, complete serotyping was not possible. Using the DNase-sequence-independent viral nucleic acid amplification method, we identified several nucleotide sequences with a high homology to human adenovirus 41 (HAdV-41) and simian adenovirus 1 (SAdV-1). However, using anti-SAdV-1 sera, it was determined that this virus was serologically different than SAdV-1. Genomic sequencing and phylogenetic analysis confirmed that this new adenovirus was so divergent from the known human adenoviruses that it was not only a new type but also represented a new species (human adenovirus G). In a retrospective clinical study, this new virus was detected by PCR in one additional patient from a separate gastroenteritis outbreak. This study suggests that HAdV-52 may be one of many agents causing gastroenteritis of unknown etiology.
              Bookmark
              • Record: found
              • Abstract: found
              • Article: found
              Is Open Access

              Cross-Species Transmission of a Novel Adenovirus Associated with a Fulminant Pneumonia Outbreak in a New World Monkey Colony

              Introduction Adenoviruses, first isolated in the 1950s from explanted adenoid tissue, are double-stranded nonenveloped DNA viruses that naturally infect many vertebrates, including humans and nonhuman primates. The human adenoviruses in the Mastadenovirus genus, comprised of all mammalian adenoviruses, are classified into 7 species A-G, and at least 51 different serotypes (and 5 proposed types, HAdV-52 to HAdV-56) have been described to date [1], [2]. Adenoviruses are the cause of an estimated 5–10% of febrile illnesses in children worldwide [3]. Some serotypes, such as human adenovirus type 14 (HAdV-14), have been associated with severe and potentially fatal outbreaks of pneumonia in residential facilities and military bases [4]. Adenoviruses have also been associated with other clinical syndromes including conjunctivitis, hepatitis, and diarrhea [5]. In nonhuman primates, most epidemiologic studies of adenoviruses have focused on their identification in fecal samples from asymptomatic animals [6], [7], [8]. Overt respiratory disease associated with simian adenoviruses has also been observed [9]. Although adenoviruses are significant pathogens, genetically modified strains are being actively explored as potential vectors for vaccines and gene therapy [10]. Infection by adenoviruses has generally been thought to be species-specific. Human adenoviruses do not usually replicate in monkey cells in the absence of helper viruses [11], and do not productively infect rodents (and vice versa) [12]. Studies of sera from animal handlers and zoo workers exposed to chimpanzees in captivity fail to detect antibodies to chimpanzee adenoviruses [13], [14]. However, recent serological surveys have found antibodies to New World and Old World monkey adenoviruses in donor human sera from regions where the monkeys are endemic [14], [15]. In addition, phylogenetic analyses of adenoviruses from greater apes reveal that they fall precisely into “human” adenoviral species B, C, and E [7]. The high degree of sequence relatedness within members of each species suggests that at least some adenoviral strains may be capable of infecting both nonhuman primates and humans. Beginning in May of 2009, a deadly outbreak of fulminant pneumonia and hepatitis occurred in a closed colony of New World titi monkeys of the Callicebus genus at the California National Primate Research Center (CNPRC). Routine microbiological testing for an infectious etiology was negative. We previously developed the Virochip (University of California, San Francisco) as a broad-spectrum surveillance assay for identifying viral causes of unknown acute and chronic illnesses [16], [17], [18], [19], [20], [21], [22]. The Virochip, a pan-viral microarray containing ∼19,000 probes derived from all viral species in GenBank (n∼2500) [21], [23], has been previously successful in detection of novel outbreak viruses such as the SARS coronavirus [22], [24] and the 2009 pandemic H1N1 influenza virus [23]. Here we apply the Virochip to identify a novel and highly divergent adenovirus as the cause of the titi monkey outbreak. In addition, we present clinical and serological evidence that this virus may have infected a researcher at the CNPRC and a family member, thus demonstrating for the first time the potential for cross-species infection by adenoviruses. Results An outbreak of fulminant pneumonia in a titi monkey colony In early 2009, the CNPRC housed 65 titi monkeys in one quadrant of an animal building. The index case, a healthy adult titi monkey, presented on May 14, 2009 with cough, lethargy, and decreased appetite (Fig. 1A, T1). Despite aggressive treatment with intravenous fluids and antibiotics, the animal developed severe respiratory distress and was humanely euthanized 5 days later. A second case presented 4 weeks later near the entrance to the building (Fig. 1A, T54). In the interim period, 3 healthy titi monkeys had been relocated from a separate building (Fig. 1A, T2, T3, and T19), with 2 of the 3 monkeys placed into the cage formerly occupied by the index case, reflecting a total at-risk population of 68. Over the ensuing 2 months, 21 additional monkeys, including one of the relocated monkeys, presented with clinical signs similar to those shown by the index case (attack rate  = 23/68, or 34%) (Figs. 1A and 1B). Clinical signs in affected animals included cough, lethargy, poor appetite, tachypnea, and abdominal breathing. These symptoms progressed to overt respiratory distress and death or humane euthanasia within an average of 8 days. Chest radiographs typically revealed diffuse interstitial pulmonary changes and bronchoalveolar consolidation indicative of pneumonia, with right middle lobe predominance (Fig. 1C). Animals displaying clinical signs were quarantined and aggressively treated by veterinarians with supplemental oxygen, anti-inflammatory medications, bronchodilators (nebulized albuterol), broad-spectrum antibiotics, and antivirals (oseltamivir and/or ribavirin). In total, 19 animals died or were euthanized due to the illness during the outbreak (case fatality rate  = 19/23, or 83%). Only 4 monkeys survived, even though the majority of sick animals (17/23, or 74%) consisted of apparently healthy adults and juveniles. Interestingly, none of the 133 rhesus macaques (Macaca mulatta) housed in the same building became sick during the outbreak, and neither did any of the Old World monkeys from surrounding outdoor colonies of rhesus and cynomolgus macaques (Macaca fascicularis). 10.1371/journal.ppat.1002155.g001 Figure 1 Clinical and epidemiologic features of the titi monkey outbreak. (A) Map of the titi monkey cages situated in one quadrant of a building, showing the locations of asymptomatic, at-risk monkeys (brown or green), affected surviving monkeys (black), and monkeys who died from their illness (skeleton). 3 monkeys were moved into the building (arrows pointing down and to the left) and 4 monkeys out of the building (arrows pointing up and to the right) during the 3rd week of the outbreak. The upper left photograph shows an image of an adult male titi monkey and his infant. The upper right inset shows the location of the titi monkey cages relative to other rhesus monkey cages in the building. Asymptomatic monkeys with positive serum antibody titers to TMAdV 4 months after the outbreak are shown in green. (B) Epidemic curve of the outbreak, with the number of cases in blue and cumulative attack rate in red. (C) Anteroposterior chest radiograph of an affected titi monkey, showing bilateral basilar infiltrates and a prominent right middle lobe consolidation. (D) 1 – gross photograph of lungs at necropsy; the lungs failed to fully collapse upon opening the chest, and a single ∼1.5 cm focus of dark red discoloration (hemorrhage) can be seen in the left caudal lobe. 2 – photomicrograph of H&E stained lung tissue showing a severe diffuse necrotizing bronchopneumonia characterized by the presence of hemorrhage and intranuclear inclusions (arrows). 3 – photomicrograph of H&E stained liver tissue showing a multifocal necrotizing hepatitis with numerous intranuclear inclusions (arrows). 4 – transmission electron micrograph of an affected lung alveolus (scale bar  = 1 µm) filled with adenovirus-like particles (inset, scale bar  = 0.1 µm). Gross necropsy findings were similar in all titi monkeys and were characterized primarily by diffuse, consolidated pneumonias, with occasional evidence of fibrinous pleuritis, pericardial/pleural edema, and hemorrhage (Fig. 1D-1). Some livers, spleens, and lymph nodes were found to be abnormally enlarged. Hepatic necrosis and hemorrhage, along with ascites, were occasionally appreciated. On histologic examination, the normal cellular architecture of the lung and trachea was destroyed, and prominent intranuclear inclusion bodies were observed in the liver, lung, and trachea (Figs. 1D-2 and 1D-3). A routine microbiological workup for infectious causes of the outbreak, including bacterial, mycoplasma, and fungal cultures, was negative. Respiratory viral testing failed to detect evidence of respiratory syncytial virus, adenovirus, influenza virus A and B, human metapneumovirus, and parainfluenza virus types 1, 2, and 3. Virochip identification, PCR screening and electron microscopic (EM) confirmation of TMAdV Given the clinical presentation of a severe acute viral respiratory illness and the appearance of intranuclear inclusion bodies on histological examination, we strongly suspected that a virus that had eluded detection by conventional assays was the cause of the titi monkey outbreak. Nasal, lung, and liver swab samples collected during necropsy were analyzed using the Virochip [21], [23]. Microarrays were analyzed using ranked Z-scores to assess the highest-intensity viral probes [18]. From a lung swab sample from an affected monkey, 4 of the top 80 probes on the Virochip corresponded to adenoviruses. Other viruses or viral families with ≥4 probes among the top 80, including chimpanzee herpesvirus (Herpesviridae), bovine viral diarrhea virus (Flaviviridae), and endogenous retroviruses (Retroviridae), were regarded as less likely to cause fulminant pneumonia and hepatitis, so were not pursued any further. The 4 adenovirus probes mapped to 2 different gene regions corresponding to the DNA polymerase and penton base (Fig. 2A). Interestingly, the 4 viral probes were derived from 2 different Adenoviridae genera (SAdV-23, simian adenovirus 23, PAdV-A, porcine adenovirus A, and HAdV-5, human adenovirus 5, in the Mastadenovirus genus; FAdV-D, fowl adenovirus D, in the Aviadenovirus genus), suggesting the presence of a divergent adenovirus that was not a member of any previously known species. 10.1371/journal.ppat.1002155.g002 Figure 2 Discovery and whole-genome characterization of the novel adenovirus TMAdV. (A) The locations of the 4 Virochip probes derived from adenovirus sequences and used to detect TMAdV are mapped onto the ∼37 kB genome. The 4 Virochip probe sequences are also aligned with the corresponding sequence in the TMAdV genome, with mismatches highlighted in pink. (B) Coverage map of deep sequencing reads corresponding to TMAdV using BLASTN (blue) and TBLASTX (transparent blue) alignments to SAdV-18. The actual coverage achieved by deep sequencing as determined by alignments to the fully sequenced genome of TMAdV is much higher (red). (C) Genome organization of TMAdV. Predicted protein coding regions are shown as boxes. Boxes above the central black line represent open reading frames (ORFs) that are encoded on the forward strand, while boxes underneath the black line represent reverse-strand encoded ORFs. Early region ORFs are shaded in gray. The x-axis refers to the nucleotide position along the ∼37k genome of TMAdV. Abbreviations: FAdV, fowl adenovirus; SAdV, simian adenovirus; PAdV, porcine adenovirus; HAdV, human adenovirus, TMAdV, titi monkey adenovirus. To confirm the Virochip finding of an adenovirus, we used consensus primers to amplify a 301 bp fragment from the hexon gene by PCR [25]. The fragment shared ∼86% nucleotide identity with its closest phylogenetic relatives in GenBank, SAdV-18, an Old World vervet monkey adenovirus, and the human species D adenoviruses. The newly identified adenovirus was designated TMAdV, or titi monkey adenovirus. Specific PCR for TMAdV was then used to screen body fluids and tissues from affected monkeys (Table 1). PCR results were positive from post-necropsy liver and lung tissues as well as from sera, conjunctival swabs, oral swabs, and nasal swabs collected at time of quarantine in 8 different affected monkeys, but were negative from a throat swab from an asymptomatic animal whose other 5 cage mates had become sick. In addition, nasal swabs were negative in 3 asymptomatic, minimal-risk titi monkeys housed in a separate building. To confirm the presence of virus in diseased tissues, we examined lung tissue from affected monkeys by transmission electron microscopy, revealing abundant icosahedral particles characteristic of adenovirus filling the alveoli (Fig. 1D-4). 10.1371/journal.ppat.1002155.t001 Table 1 PCR screening for TMAdV. Sample Sample Type PCR Result Date Presenting with Clinical Signs Date of Necropsy Affected, at-risk titi monkeys (died) T1 serum§ − 5/14/2009 5/19/2009 T26 serum¶ + 7/23/2009 7/30/2009 T28 conjunctival swab¶ + 7/16/2009 7/25/2009 nasal swab¶ + liver swab¶ + lung swab¶ + T29 serum¶ + 7/26/2009 7/31/2009 T30 serum¶ − 7/25/2009 7/30/3009 T33 lung swab¶ + 6/23/2009 6/29/2009 nasal swab¶ + T36 lung swab¶ + 7/7/2009 7/14/2009 lung swab¶ + T60 serum¶ − 7/15/2009 7/22/2009 T63 serum¶ − 6/20/2009 8/1/2009 T67 nasal swab# − 7/7/2009 8/13/2009 nasal swab¶ + Affected, at-risk titi monkeys (survived) T31 serum* − 7/10/2009 N/A T32 serum* − 7/12/2009 N/A T34 serum* − 6/23/2009 N/A T62 serum* − 7/8/2009 N/A Asymptomatic at-risk and minimal-risk titi monkeys T27 throat swab (n = 1)¶ − N/A N/A at-risk titi stool from cages (n = 14)* − N/A N/A at-risk titi serum (n = 29)* − N/A N/A minimal-risk titi oral swab (n = 3)¶ − N/A N/A minimal-risk titi stool from cages (n = 5)* − N/A N/A minimal-risk titi serum (n = 8)* − N/A N/A minimal-risk titi stool from cages (n = 8)* − N/A N/A minimal-risk titi serum (n = 8)* − N/A N/A Other rhesus rectal swabs (n = 26)* − N/A N/A human serum (n = 15)∞ − N/A N/A rodent droppings (n = 2)* − N/A N/A For titi monkey cage designations (TXX), please refer to Fig. 1. §: initial case. # collected prior to outbreak. ¶: collected during outbreak. *collected 2 months after outbreak. ∞: collected 4 months after outbreak. Next, to assess persistent subclinical infection from TMAdV, we analyzed serum samples from at-risk asymptomatic or affected surviving monkeys 2 months after the outbreak (n = 41). All post-outbreak serum samples were negative for TMAdV by PCR (Table 1). To assess potential TMAdV shedding, stool samples collected from all cages housing titi monkeys 2 months post-outbreak were analyzed by PCR (n = 27), and were found to be negative. In addition, we checked for TMAdV in rectal swab samples from rhesus macaques housed in the same building as the titi monkeys (n = 26) and in pooled droppings from wild rodents (n = 2) living near the titi monkey cages. All macaque and rodent fecal samples were negative for TMAdV by PCR. We also sought to determine whether PCR assays commonly used to detect human adenoviruses in clinical or public health settings could detect TMAdV. Adenovirus PCR was performed on a TMAdV-positive clinical sample, a TMAdV culture, and a human adenovirus B culture (as a positive control) using an additional 5 pairs of primers, according to previously published protocols [26], [27], [28] Three of the 5 primer pairs, designed to detect human respiratory adenoviruses B, C, and E, failed to amplify TMAdV [27]. The remaining 2 pairs of primers, both derived from highly conserved sequences in the hexon gene [26], [28], were able to detect TMAdV in culture as well as directly from clinical material. Whole-genome sequencing, features, and phylogenetic analysis of TMAdV To facilitate whole-genome sequencing of TMAdV, deep sequencing of a lung swab from one affected titi monkey and lung tissue from another affected monkey was performed. Out of ∼11.9 million high-quality reads, 2,782 reads and 3,767 reads aligned to the SAdV-18 genome by BLASTN (Fig. 2B, blue) and TBLASTX (Fig. 2B, transparent blue), respectively, with reads mapping to sites across the genome. De novo assembly of the complete TMAdV genome from reads that aligned to SAdV-18 was not possible due to insufficient sequence coverage ( 1∶512, while 1 monkey exhibited a positive but much weaker response. Affected titi monkeys who died during the outbreak exhibited a wide range of neutralizing Ab titers, from 1∶512 (those without Ab likely died before mounting a response). 10.1371/journal.ppat.1002155.g006 Figure 6 Seroprevalence of TMAdV in humans and monkeys. Sera from titi monkeys (circles), rhesus macaques (squares), and humans (triangles) were tested for antibodies to TMAdV by virus neutralization. Arrows designate pre-outbreak and post-outbreak serum samples from the same individual monkey. Pre-outbreak serum samples were previously banked in 2007. Sera from CNPRC personnel and close contacts (orange triangles) were collected 4 months post-outbreak, except for the two family members of the clinically ill researcher, whose sera were collected 1 year post-outbreak. *, clinically ill researcher; **, family member of the researcher, who was also sick. Abbreviations: CNPRC, California National Primate Research Center; NEG, negative. To investigate the possibility of subclinical infection by TMAdV, we also examined serum samples from asymptomatic titi monkeys (n = 40) and nearby rhesus macaques (n = 36), collected 2 months after the outbreak. Fourteen of 40 asymptomatic titi monkeys tested (35%) had antibody to TMAdV, indicating that the incidence of subclinical infection was significant (Fig. 1A; Fig 6). In fact, one of the 14 asymptomatic titi monkeys with positive Ab titers was located in the minimal-risk building. In contrast, only 1 of 36 rhesus macaque samples was positive, with an Ab titer of 1∶16. The 1 antibody-positive rhesus serum sample was negative by specific PCR for TMAdV (data not shown), as was stool from the cage in which the rhesus monkey was housed (Table 1). Approximately 4 months after the outbreak, serum samples were collected from CNPRC personnel in direct contact with the titi monkeys. Serum samples were also collected from the two family members of the clinically ill CNPRC researcher 1 year after the outbreak. Only two samples were found positive for neutralizing Abs to TMAdV: (1) Ab titers for the clinically ill researcher were 1∶32, and (2) Ab titers for one of the family members of the clinically ill researcher were 1∶8. Among 81 random blood donors from the Western United States, 2 individuals (2/81, 2.5%) had positive Ab titers of 1∶16 and 1∶8. Pooled rabbit sera containing antibodies to human adenovirus serotypes 1 through 35, representing species A–E, were unable to neutralize TMAdV (data not shown). Thus, the results of our serological survey appear unlikely to be due to nonspecific cross-reactivity from prior exposure to known human adenoviruses. Discussion In this study, we employed a pan-viral microarray assay, the Virochip, to identify a novel adenovirus associated with a fulminant pneumonia outbreak in a colony of New World titi monkeys. Despite the absence of an animal model, which precludes a strict fulfillment of Koch's postulates, there are several lines of evidence implicating this novel adenovirus, TMAdV, as the cause of the outbreak. First, conventional testing for other pathogens, including other viruses by Virochip, was negative, and affected monkeys did not respond to empiric therapy with antibiotics or antivirals (ribavirin and oseltamivir in anecdotal use are not effective against adenoviral infections) [30]. Second, the clinical presentation of pneumonia and hepatitis is consistent with the known spectrum of disease associated with adenoviral infections. Third, TMAdV sequence was recovered by PCR in various body fluids and tissues from affected monkeys, including blood, respiratory secretions, and lung/liver tissue (Table 1). Fourth, the finding of intranuclear inclusions in diseased tissues, as well as direct visualization of adenoviral-like particles (TMAdV) in lung alveoli by electron microscopy (Figs. 1D-2 to 1D-4), support a primary role for TMAdV in the pathogenesis of tissue injury in affected monkeys. Finally, there was a significant neutralizing Ab response in surviving animals, with 2 monkeys having titers undetectable prior to the outbreak but rising to >1∶512 at convalescence (Fig. 6). Although TMAdV retains the core genomic features common to all adenoviruses (Fig. 2C), phylogenetic analysis clearly places TMAdV within a separate branch, with no closely related neighbors (Figs. 3 and S1). A phylogenetic distance of >10% combined with the lack of cross-neutralization defines TMAdV as a new species [31]. Since emerging adenovirus strains such as HAdV-14 and HAdV-D22/H8 (otherwise known as HAdV-D53) are known to arise from recombination events among related ancestral strains [32], [33], we performed bootscanning analysis to look for such events in TMAdV. The bootscanning analysis, however, failed to show evidence of recombination, likely because closely related and/or ancestral strains to TMAdV have not yet been identified. Entry of adenoviruses into cells involves an initial attachment of the fiber knob to the cell receptor, followed by internalization via a secondary interaction of the penton base with αv integrins [34], [35]. The presence of an RGD motif in the TMAdV penton base implies that the virus uses αv integrins for internalization [35]. However, the high sequence divergence in the fiber protein (Table 2), as well as the absence of fiber motifs conserved among adenoviruses that bind CAR [36], [37] (coxsackievirus-adenovirus receptor) or CD46 [38], [39], [40] (data not shown), suggest that neither of these two human adenoviral receptors may be the attachment receptor for TMAdV. Further studies will be necessary to identify the preferred cellular attachment and internalization receptors for TMAdV. Despite its isolation from affected titi monkeys, we were unable to propagate TMAdV in both established (BSC-1) and primary (PMK) monkey kidney cells (Fig. 4). The virus, however, grew efficiently in a human A549 lung adenocarcinoma cell line. One explanation for this finding is that TMAdV may be unable to productively infect cells derived from Old World monkeys (e.g. rhesus and African green monkeys). An alternative possibility is that successful propagation of TMAdV may depend on infection of a specific host cell type, such as A549 lung, and not BSC-1 or PMK kidney cells. Nevertheless, after 10 passages in human A549 cells, the fully adapted strain of TMAdV exhibits an extended host range with the ability to productively infect both monkey and human cells. This observation implies that TMAdV possesses an inherent capacity to cross the species barrier and infect both humans and nonhuman primates. Efforts to identify host range and cell tropism of TMAdV, as well as the specific sequence changes responsible for adaptation to growth in cell culture, are currently underway. The virulence of TMAdV in healthy and apparently immunocompetent titi monkeys (83% case fatality rate) is highly unusual for infections by adenovirus. In humans, deaths due to adenovirus infections or outbreaks are generally low (up to 18% for pneumonia associated with HAdV-14 [4]). Furthermore, severe infections from human adenoviruses are more commonly associated with older age, immunosuppression, and chronic underlying conditions such as kidney failure [4], [41]. Young, healthy individuals are in general much less likely to succumb to adenoviral-related illness. The severity of TMAdV-related illness in affected titi monkeys suggests that this species of monkey may not be the natural host for the virus. The failure to detect fecal shedding of TMAdV in convalescent or asymptomatic animals also suggests that the virus does not normally infect titi monkeys (Table 1). Although the exact origin of TMAdV remains unclear, we can speculate on several possibilities. One possibility is that a cross-species “jump” from captive macaques to a susceptible colony of titi monkeys precipitated the outbreak. As there have been no new introductions of monkeys into the closed colony for the past 2 years, this conjecture relies on asymptomatic infection and transmission of TMAdV in the captive rhesus/cynomolgus macaque population at the CNPRC. CNPRC personnel who visited macaque areas would occasionally enter titi rooms with no change in personal protective equipment, thus providing a potential route of transmission for the virus. In addition, specific antibodies were detected in 1 of 36 (2.8%) asymptomatic rhesus macaques housed in the same building (Fig. 6), indicating that TMAdV has the capacity to infect this species of Old World monkey. Notably, the closest identified phylogenetic relative to TMAdV among the complete genomic sequences available in GenBank is a rhesus monkey adenovirus, SAdV-3 (Fig. 4; Table 2). Furthermore, serological evidence for cross-species adenoviral transmission events between different nonhuman primate species has been previously reported in the literature [42]. Although we failed to detect TMAdV in rodent droppings found near titi monkey cages (Table 2), it is still possible that the virus arose from an unknown animal reservoir. In this regard, the high sequence divergence of TMAdV relative to the known human/simian adenoviruses (Fig. 3), and comparable sequence similarity in the polymerase gene to a porcine adenovirus (Figs. 3 and S1) are striking. The four-week interval between the index case and the second case appears overly long given a typical incubation period for adenovirus infections of no more than 1 week [43]. This may be explained by our finding of a high rate of subclinical infection by TMAdV in asymptomatic titi monkeys (35%), but may also be due to separate introductions of TMAdV into the colony from an as-yet unidentified reservoir. Our study data also support the potential for cross-species transmission of TMAdV between monkeys and humans. The researcher's fever, cough, and pleuritic symptoms (“burning sensation in the lungs”) are consistent with the development of a prolonged viral respiratory illness. Interestingly, pleurisy has been specifically reported in association with certain human adenovirus infections [44]. The clinical presentation, time of illness concurrent with the onset of the outbreak, and presence of neutralizing Abs in convalescent serum all strongly point to primary infection of the researcher by TMAdV. The detection of weakly neutralizing Abs (1∶8) in a serum sample from a sick family member of the researcher also suggests that TMAdV may be capable of human-to-human transmission. The decreased levels of neutralizing Abs to TMAdV in the researcher (1∶32) and a family member (1∶8) relative to those in infected titi monkeys (up to >1∶512) are consistent with a recent study showing much higher levels of neutralizing antibodies in chimpanzees than in humans with adenovirus infections, possibly due to more robust adenovirus-specific T-cell responses in humans than in monkeys [45]. Several lines of evidence support the contention that the direction of TMAdV transmission was zoonotic (monkeys to humans) rather than anthroponotic (humans to monkeys). First, the closest known relative to TMAdV in GenBank is SAdV-3, an Old World monkey adenovirus (Fig. 3; Table 2). Second, our results show that PCR assays for human adenoviruses in common use are capable of detecting TMAdV. Although sequencing of PCR amplicons for human adenoviruses is not performed routinely in diagnostic virology, TMAdV would presumably have been detected previously in large-scale studies of hexon sequencing of Ad field isolates if it were circulating in the community [46], [47]. Finally, the available sequence data in GenBank is heavily biased towards human adenoviruses, and much less is known about the potential diversity of the simian adenoviruses. We also cannot formally exclude the possibility that the outbreak arose from anthroponotic transmission. In our study, 2 of 81, or 2.5% of random adult blood donors exhibited borderline titers of neutralizing antibody to TMAdV, indicating either a low prevalence of TMAdV in the human population or cross-reactivity to a related virus (although no evidence of cross-reactivity was found with HAdV serotypes 1 through 35). Future large-scale studies of TMAdV seroepidemiology will be needed to better understand transmission of TMAdV between monkeys and humans. Nevertheless, our discovery of TMAdV, a novel adenovirus with the capacity to cross species barriers, highlights the need to monitor adenoviruses closely for outbreak or even pandemic potential. Materials and Methods Ethics statement This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. The use and care of all animals followed policies and guidelines established by the University of California, Davis Institutional Animal Care and Use Committee (IACUC) and CNPRC (Animal Welfare Assurance #A3433-01). The protocol for the maintenance and breeding of the titi monkey colony was approved by the University of California, Davis IACUC (Protocol #15730). No specific animal research protocol was drafted for this study, as only excess clinical samples were analyzed for diagnostic purposes. Animals in extreme respiratory distress were humanely euthanized by veterinarians. Extensive veterinary care was provided to all animals affected by the outbreak in order to minimize pain and distress. Serum samples from staff at the CNPRC, close contacts, and random adult blood donors were collected under protocols approved by institutional review boards of the University of California, Davis (Protocol #200917650-1) and University of California, San Francisco (Protocol #H49187-35245-01). Specifically, written informed consent was obtained from staff at the CNPRC and close contacts for analysis of their samples. Any potentially identifying information has been provided with the explicit permission of the individuals involved. Sera from random blood donors were obtained from the Blood Systems Research Institute (San Francisco, CA); sera were derived from affiliated donor banks in California (Blood Centers of the Pacific, San Francisco, CA), Nevada (United Blood Service, Reno, NV), and Wyoming (United Blood Services, Cheyenne, Wyoming) and de-identified prior to analysis. The California National Primate Research Center (CNPRC) The California National Primate Research Center (CNPRC), which houses over 5,000 nonhuman primates, is a part of the National Primate Research Centers Program and is accredited by the Association for the Assessment and Accreditation of Laboratory Animal Care, International (AAALAC). At the beginning of 2009, the CNPRC maintained a colony of 74 titi monkeys (Callicebus cupreus) and a colony of over 4,500 rhesus macaques (Macaca mulatta). No new animals have been introduced into either colony for over 2 years. All titi monkeys are maintained in small social groups, while rhesus macaques are maintained in small or large social groups. All animal facilities are maintained in compliance with United States Department of Agriculture specifications. Eighty-eight percent of the titi monkey population (n = 65) were housed in 1 quadrant of an indoor animal building, and all titi monkeys demonstrating clinical signs originated from this building (i.e. “at-risk” room) (Fig. 1A). Rhesus macaques (n = 133) were housed in the other 3 quadrants of this same building, and surrounding the building were outdoor housing units with rhesus macaques and cynomolgus macaques (Macaca fasicularis). Three additional titi monkeys were moved into the at-risk room less than 2 weeks after presentation of the index case, reflecting a total at-risk population of 68 animals. The remaining 6 titi monkeys were housed in an indoor animal building greater than 500 yards from the at-risk population (i.e. “minimal-risk” room). Outbreak investigation and microbiological testing The outbreak lasted approximately 3 months from May to August of 2009. Affected titi monkeys died from 3–24 days after appearance of clinical signs, with an average time to death or humane euthanasia of 8 days. Clinical and epidemiological data, including daily census reports, were tracked and recorded by veterinary and management staff. All personnel entering the titi monkey rooms (both at-risk rooms and minimal-risk rooms) needed to pass within approximately 20 feet of macaque enclosures prior to entry. CNPRC personal protective equipment (PPE) policy requires a change of PPE between entrance/exit of animal rooms housing different species. Staff compliance of this policy may have been compromised. Measures have since been taken by CNPRC management to ensure compliance with existing policies. Bacterial, mycoplasma, and fungal cultures were performed at the CNPRC. Clinical samples were also sent to an outside laboratory (Focus Diagnostics, Cypress, CA) for respiratory viral testing by centrifugation-enhanced shell vial culture followed by direct fluorescent antibody staining for 8 viruses (respiratory syncytial virus, adenovirus, influenza virus A and B, parainfluenza virus types 1, 2, and 3, and human metapneumovirus). Gross, histopathological, and ultrastructural analyses Gross and histopathological analyses of post-mortem tissues were performed by a board-certified veterinary pathologist specializing in nonhuman primate/laboratory animal medicine, a branch of Primate Services at the CNPRC. At necropsy, tissue samples from the trachea, lung, and liver were collected and fixed in 10% formalin. Tissues were routinely processed and embedded in paraffin. 3-µm sections were stained with hematoxylin and eosin (HE) and examined by light microscopy. For transmission electron microscopy, tissue fragments (2×2 mm) were excised from paraffin blocks of lung, deparaffinized, and processed as previously described [48]. Nucleic acid extraction and cDNA library preparation Total nucleic acid was extracted from body fluid or swab samples using commercially available kits (Qiagen, Valencia, CA). 200 µL of sample were passed through a 0.22 µm filter (Millipore, Temecula, CA) to remove bacteria and cellular debris and then treated with Turbo DNase (Ambion, Culver City, CA) to degrade host genomic DNA prior to extraction. For tissue samples, lung or liver tissue was homogenized in a 15 mL Eppendorf tube using a disposable microtube pestle (Eppendorf, San Diego, CA) and scalpel, and RNA extraction was then performed using TRIzol LS (Invitrogen, Carlsbad, CA), followed by isopropanol precipitation and two washes in 70% ethanol. Extracted nucleic acid was amplified using a random PCR method to generate cDNA libraries for Virochip and deep sequencing analyses as previously described [18], [21]. Virochip analysis The current 8×60 k Virochip microarrays used in this study contain 19,058 70mer probes representing all viral species in GenBank, and combine probes from all previous Virochip designs [17], [18], [21], [23]. Four probes derived from 2 different Adenoviridae genera (SAdV-23, PAdV-A, HAdV-5, and FAdV-D) yielded an adenovirus signature on the Virochip that was found to be TMAdV. With the exception of SAdV-23, these highly conserved probes are part of the core Virochip design and were derived from all available adenoviral sequences in GenBank as of 2002 [21]. One explanation why more high-intensity probes to simian adenoviruses were not seen by Virochip analysis is that the genomes of many simian Ads, including SAdV-3 and SAdV-18 (the two closest phylogenetic relatives to TMAdV in GenBank), were not sequenced until after 2004 [7], [49], and thus their genomes are not as well-represented on the Virochip microarray. Virochip analysis was performed as previously described [21], [23]. Briefly, samples were labeled with Cy3 or Cy5 fluorescent dye, normalized to 10 pmol of incorporated dye, and hybridized overnight using the Agilent Gene Expression Hybridization kit (Agilent Technologies, Santa Clara, California). Slides were scanned at 3 µm resolution using an Agilent DNA Microarray Scanner. Virochip microarrays were analyzed with Z-score analysis [18], hierarchical cluster analysis [50], and E-Predict, an automated computational algorithm for viral species prediction from microarrays [51]. Only Z-score analysis, a method for assessing the statistical significance of individual Virochip probes, yielded a credible viral signature on the microarray. PCR screening We initially used consensus primers derived from a highly conserved region of the hexon gene to confirm the Virochip finding of an adenovirus by PCR [25]. After recovering the full genome sequence, we then designed a set of specific PCR primers for detection of TMAdV, TMAdV-F (5′-GTGACGTCATAGTTGTGGTC) and TMAdV-R (5′-CTTCGAAGGCAACTACGC). The TMAdV-specific quantitative real-time PCR was performed on a Stratagene MX3005P real-time PCR system using a 25 µL master mix consisting of 18 µL of water, 2.5 µL of 10X Taq buffer, 1 µL of MgCl2 (50 mM), 0.5 µL of deoxynucleoside triophosphates (dNTPs; 12.5 mM), 0.5 µL of each primer (10 µM), 0.5 µL SYBR green, 0.5 µL of Taq polymerase (Invitrogen, Carlsbad, CA), and 1 µL of extracted nucleic acid. Conditions for the PCR reaction were 40 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 30 s. Amplicons were purified on a 2% agarose gel, cloned into plasmid vectors using TOPO TA (Invitrogen, Carlsbad, CA), and sent to an outside company (Elim Biopharmaceuticals, Hayward, CA) for Sanger sequencing in both directions using vector primers M13F and M13R. To assess linearity and limits of sensitivity for the TMAdV PCR assay, 12 serial log dilutions were made of a standard plasmid constructed by cloning the 157-bp TMAdV amplicon into a TOPO plasmid vector. Purified plasmid clones at each serial dilution were quantified using a Nanodrop spectrophotometer and then spiked into negative serum, stool, or oral swab sample matrix, each matrix consisting of a pool of 10 sera, 10 stool samples, or 3 oral swabs, respectively. For each sample type, a standard curve for the TMAdV PCR assay was calculated from 3 PCR replicates at each dilution of nucleic acid extracted from the spiked matrix (data not shown). To determine limits of sensitivity for the assay, probit analysis of results from 6 PCR replicates of 7 serial log dilutions (from a starting concentration of ∼1.2×105 copies/mL) was performed using SPSS 16.0 (SPSS Inc., Chicago, IL). By probit analysis, the 95% limit of detection for TMAdV was 781, 377, or 35 viral genome equivalents/mL for serum, stool, or oral swab samples, respectively (data not shown). Whole-genome sequencing To facilitate whole-genome sequencing of TMAdV, we prepared amplified cDNA/DNA libraries for deep sequencing from lung tissue and a lung swab sample from 2 different monkeys using previously published protocols [23], [52]. Briefly, randomly amplified libraries were cleaved with a Type IIs restriction endonuclease (GsuI) and truncated adapters were ligated on the resulting strand ends. Full-length adapters containing strict 6-nt barcodes were added via an additional 15 cycles of PCR. Amplified libraries were size-selected on a 2% agarose gel at approximately 350 bp average length and then sent to an outside company (Elim Biopharmaceuticals, Hayward, CA) for deep sequencing on an Illumina Genome Analyzer IIx (Illumina, San Diego, CA). Paired-end reads were sequenced for 73 cycles in each direction. Paired-end reads were subsequently filtered to eliminate low-complexity sequences with a Lempel-Ziv-Welch (LZW) compression ratio below 0.4 [53], split into individual reads, classified by barcode, and stripped of any remaining primer sequences using BLASTN alignments (word size = 11, E-value = 1×10−5). After low-complexity filtering and barcode/primer trimming, 11,950,557 sequence reads remained, with each read consisting of 67 nucleotides, for a total of ∼800 megabases of sequence. Reads were then aligned using BLASTN (word size = 11, E-value = 1×10−5) and TBLASTX (word size = 11, E-value = 1×10−5) to the genome sequence of SAdV-18 (Fig. 2B). Overlapping reads aligning to SAdV-18 were used to assemble portions of the TMAdV genome with Geneious software (version 3.6.5) [54], employing the SAdV-18 genome as a reference sequence and requiring a 20-bp minimum overlap and 95% overlap identity. Aligning reads were also used to design PCR primers to close remaining gaps in the TMAdV genome. Amplicons derived from specific TMAdV PCR primers were gel-purified, cloned, and sequenced as described above. The 5′ end corresponding to the inverted terminal repeat (ITR) of TMAdV was obtained by PCR using a forward degenerate consensus primer and a reverse TMAdV-specific primer. The 3′ end was recovered using a forward primer near the 3′ end of the genome and a reverse primer derived from 5′-ITR sequence. Structural features and phylogenetic analysis To identify predicted coding regions in the TMAdV genome, we used the fully annotated genome sequence of SAdV-21 in GenBank as a reference. First, we aligned the two genomes and identified all ORFs that were present with Geneious [54]. We then selected the candidate ORF that best matched the corresponding ORF in the annotated reference genome. For adenovirus genes that are spliced (e.g. E1A), the identification of a GT-AG intron start-stop signal was used to pinpoint the correct ORF. To confirm the accuracy of the coding sequence, the sequence of each identified ORF was aligned to a database containing all adenoviral proteins in GenBank by BLASTX. To generate whole-genome and individual gene nucleotide phylogeny trees, all 95 fully sequenced unique adenovirus genomes were first downloaded from GenBank. Multiple sequence alignments were then performed on a 48-core Linux system using ClustalW-MPI [55]. Trees were constructed after bootstrapping to 1000 replicates by the neighbor-joining method (based on Jukes-Cantor distances) in Geneious [54], [56]. Pairwise alignments were calculated using Shuffle-LAGAN (window size, 400 bp; step size 40 bp; translated anchoring), a glocal alignment algorithm that is able to calculate optimal alignments by using both local alignments and global maps of sequence rearrangements (e.g. duplications of the fiber gene in adenovirus genomes with 2 fibers) [57]. Sliding window analysis of the Shuffle-LAGAN pairwise alignments was performed using the online mVISTA platform [58]. More accurate alignments were obtained with Shuffle-LAGAN than with either ClustalW-MPI or Geneious (data not shown). Bootscanning analysis was performed according to the Kimura 2-parameter method using 1000 replicates with Simplot (version 3.5.1) [59]. Pairwise amino acid amino acid alignments between predicted TMAdV proteins and corresponding proteins in other adenoviruses (Table 2) were performed using Geneious [54]. Virus cultivation A549 (human lung adenocarcinoma) and BSC-1 (African green monkey kidney epithelial) cell lines as well as PMK (primary rhesus monkey kidney) cells are routinely maintained at the Viral and Rickettsial Disease Laboratory (VRDL) branch of the California Department of Public Health. Media consisting of Hank's medium (for A549 cells) or Dulbecco's modified Eagle's medium (DMEM) (for BSC-1 cells) were supplemented with 1×nonessential amino acids (Invitrogen, Carlsbad, CA), 10% fetal bovine serum, 100 U of penicillin/mL and 100 µg of streptomycin/mL. PMK cells were maintained in tubes containing growth media and antibodies to SV-40 and SV-5 polyomaviruses (Viromed, Pasadena, CA). Clinical samples were clarified by centrifugation for 10 min×4000 g and passaged through a 0.2-µm filter. Cell culture passages were subjected to 3 freeze-thaw cycles and clarified as above. After achieving 80–90% confluency, cell culture media were changed to maintenance media with 2% FBS and were inoculated with 200 µL of clinical sample or 100 µL of passaged viral supernatant. Viral replication was monitored over 14 days by visual inspection under light microscopy for cytopathic effect (CPE). To confirm the generation of infectious virus, viral supernatants were quantitated by an end-point dilution assay. Virus neutralization assay (human and monkey sera) A virus stock of TMAdV (passage 7) was produced on human A549 cells, aliquoted, and quantitated by end-point dilution. To perform the virus neutralization assay, 55 µL of viral supernatant at a concentration of 100 TCID50 and 55 µL of serum (starting at a 1∶8 dilution) were mixed and incubated for 1 hour at 37°C. As a control for each serum sample, 55 µL of culture media and 55 µL of serum were mixed and treated in an identical fashion. While mixtures were incubating, A549 cells grown in T-25 plates were trypsinized and 4,000 cells in 100 µL of media were added to each well of a 96-well plate. After incubation, 100 µL of mixture were inoculated into appropriate wells containing 4,000 cells per well and the entire plate was placed in a 37°C 5% CO2 incubator. Cells in the plate wells were observed for evidence of CPE every other day for 1 week. For wells that showed inhibition of viral CPE, the corresponding serum samples were diluted in six 2-fold steps and then retested. The reciprocal of the highest dilution that completely inhibited viral CPE was taken as the neutralizing antibody titer. Virus neutralization assay (rabbit typing sera) To assess cross-neutralization of TMAdV by known human adenoviruses, 7 pools of in-house reference sera at the VRDL (rabbit hyperimmune sera) and collectively containing antibodies to human adenovirus serotypes 1 through 35 were available for testing. For each pool, 55 µL of rabbit sera and 55 µL of viral supernatant at a concentration of 100 TCID50 were mixed, incubated for 1 hour at 37°C, and inoculated onto A549 cells in wells of a 96-well plate as described above. Cells in the plate wells were observed for evidence of CPE every other day for 1 week. Microarray and nucleotide sequence accession numbers All Virochip microarrays used in this study have been submitted to the NCBI GEO database (study accession number GSE26898; microarray accession numbers GSM662370-GSM662391; microarray design accession number GPL11662). The annotated, whole-genome sequence of TMAdV has been submitted to GenBank (accession number HQ913600). Deep sequencing reads have been submitted to the NCBI Sequence Read Archive (accession number SRA031285). Supporting Information Figure S1 Phylogenetic analysis of the hexon, polymerase, penton base, and fiber genes of TMAdV. A multiple sequence alignment of selected genes from all 95 unique, fully-sequenced adenovirus genomes in GenBank and TMAdV is performed and the results displayed as a radial phylogenetic tree. The branch corresponding to TMAdV is highlighted in boldface red. Abbreviations: HAdV, human adenovirus, SAdV, simian adenovirus; PAdV, porcine adenovirus; FAdV, fowl adenovirus. (TIF) Click here for additional data file. Figure S2 Bootscanning recombination analysis of TMAdV. Bootscanning analysis was initially performed with all 95 unique, fully-sequenced adenovirus genomes in GenBank (data not shown). After removal of similar viral genomes, bootscan plots of the whole genome and individual genes from a subset representing human/simian adenoviruses in species A–G and all non-primate vertebrate adenoviruses were generated. The window size is 400 bp with a step size of 40 bp for the whole genome, and 200 bp with a step size of 20 bp for the individual genes. The x-axis refers to the nucleotide position. For definition of abbreviations, please refer to Fig. 3. (TIF) Click here for additional data file.
                Bookmark

                Author and article information

                Contributors
                clange_virology@gmx.de
                Journal
                Virol J
                Virol. J
                Virology Journal
                BioMed Central (London )
                1743-422X
                27 November 2019
                27 November 2019
                2019
                : 16
                : 147
                Affiliations
                [1 ]Metabiota Inc, Nanaimo, BC Canada
                [2 ]National Laboratory of Public Health, Brazzaville, Republic of the Congo
                [3 ]Wildlife Conversation Society, Bronx, NY USA
                [4 ]Unites States Fish and Wildlife Service, Crossroads VA, Bailey’s, USA
                [5 ]Ministry of Agriculture and Livestock, Brazzaville, Republic of the Congo
                [6 ]ISNI 0000 0004 1936 9684, GRID grid.27860.3b, One Health Institute, School of Veterinary Medicine, , University of California, ; Davis, CA USA
                [7 ]ISNI 0000 0004 0409 4702, GRID grid.420826.a, EcoHealth Alliance, ; New York, NY USA
                [8 ]Ministry of National Defense, Brazzaville, Republic of the Congo
                [9 ]British Columbia Ministry of Environment and Climate Change Strategy, Victoria, BC Canada
                Article
                1248
                10.1186/s12985-019-1248-z
                6880561
                31775793
                0905122f-f5c7-4bec-a6f8-f5c2d25b517c
                © The Author(s). 2019

                Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License ( http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

                History
                : 16 September 2019
                : 5 November 2019
                Funding
                Funded by: USAID
                Award ID: GHN-A-OO-09-00010-00
                Categories
                Research
                Custom metadata
                © The Author(s) 2019

                Microbiology & Virology
                adenovirus,primate,zoonosis,africa,evolution,bushmeat
                Microbiology & Virology
                adenovirus, primate, zoonosis, africa, evolution, bushmeat

                Comments

                Comment on this article