Detection of skunk adenovirus-1 in domestic ferrets (Mustela putorius).

We describe gross pathology, histopathology, ultrastructure, and molecular characterization of skunk adenovirus 1 (SkAdV-1) in 3 ferrets (Mustela putorius) submitted to the Infectious Diseases Laboratory, Zoo and Exotic Animal Pathology Service at the University of Georgia. Case 1 came from a cohort with a 6-week history of multiple ferrets with respiratory disease and mortality. Cases 2 and 3 were submitted 19 days later; they died 3 days after arrival at the vendor's facility. Histopathological evaluation found evidence of suppurative bronchopneumonia, with intralesional gram-positive cocci in case 1. The others had evidence of subacute multifocal histiocytic and lymphoplasmacytic interstitial pneumonia. Viral isolation, polymerase chain reaction (PCR), and DNA sequencing identified SkAdV-1 in lung tissue. In situ hybridization confirmed adenovirus within the lung lesion, and adenovirus particles were visible under an electron microscope. This is the seventh species in which SkAdV-1 has been identified. Cross-species transmission events have important implications for wildlife management and multispecies facilities.

An Exon-Capture System for the Entire Class Ophiuroidea.

Exon-capture studies have typically been restricted to relatively shallow phylogenetic scales due primarily to hybridization constraints. Here, we present an exon-capture system for an entire class of marine invertebrates, the Ophiuroidea, built upon a phylogenetically diverse transcriptome foundation. The system captures approximately 90% of the 1,552 exon target, across all major lineages of the quarter-billion-year-old extant crown group. Key features of our system are 1) basing the target on an alignment of orthologous genes determined from 52 transcriptomes spanning the phylogenetic diversity and trimmed to remove anything difficult to capture, map, or align; 2) use of multiple artificial representatives based on ancestral state reconstructions rather than exemplars to improve capture and mapping of the target; 3) mapping reads to a multi-reference alignment; and 4) using patterns of site polymorphism to distinguish among paralogy, polyploidy, allelic differences, and sample contamination. The resulting data give a well-resolved tree (currently standing at 417 samples, 275,352 sites, 91% data-complete) that will transform our understanding of ophiuroid evolution and biogeography.

An Acute Mortality Event Associated with Novel Macrorhabdus ornithogaster Infection and Underlying Factors in a Newly-Established Captive Group of American White Ibis (Eudocimus albus) Nestlings.

Twenty-four American white ibis (Eudocimus albus) nestlings were collected in Florida (USA) on 17 April 2017 to establish a captive flock. On 7 May 2017, three birds died suddenly, following severe lethargy, hemorrhaging from the mouth and nares, anorexia, and production of bright-green colored feces. An additional ibis with delayed growth and pathological fractures was euthanized 18 May 2017. Severe ventriculitis associated with Macrorhabdus ornithogaster was noted in all four birds, bacterial sepsis was confirmed in one bird by culture and histologic examination, and bacterial endotoxemia was suspected in two birds based on gross and histologic examination, but no bacteria were isolated from these birds. Birds also had vitamin E liver levels consistent with coagulopathy previously described in pelicans. We sampled feces from 91 adult, free-living, healthy ibis in Florida in July 2017 and found 71% were shedding organisms with morphologic characteristics consistent with Macrorhabdus sp. Molecular characterization of the ibis-origin M. ornithogaster showed it was phylogenetically related to numerous M. ornithogaster sequences. It is unknown if M. ornithogaster infection resulted in clinical disease as a result of dietary or stress-related dysbiosis, or other factors. Macrorhabdus-associated disease has not previously been confirmed in wading birds. We discuss potential associations of gastric M. ornithogaster infection with morbidity and mortality in these cases and highlight the need for additional studies on this pathogen in free-living birds.

Chlamydia buteonis in birds of prey presented to California wildlife rehabilitation facilities.

Chlamydial infections, caused by a group of obligate, intracellular, gram-negative bacteria, have health implications for animals and humans. Due to their highly infectious nature and zoonotic potential, staff at wildlife rehabilitation centers should be educated on the clinical manifestations, prevalence, and risk factors associated with Chlamydia spp. infections in raptors. The objectives of this study were to document the prevalence of chlamydial DNA shedding and anti-chlamydial antibodies in raptors admitted to five wildlife rehabilitation centers in California over a one-year period. Chlamydial prevalence was estimated in raptors for each center and potential risk factors associated with infection were evaluated, including location, species, season, and age class. Plasma samples and conjunctiva/choana/cloaca swabs were collected for serology and qPCR from a subset of 263 birds of prey, representing 18 species. Serologic assays identified both anti-C. buteonis IgM and anti-chlamydial IgY antibodies. Chlamydial DNA and anti-chlamydial antibodies were detected in 4.18% (11/263) and 3.14% (6/191) of patients, respectively. Chamydial DNA was identified in raptors from the families Accipitridae and Strigidae while anti-C.buteonis IgM was identified in birds identified in Accipitridae, Falconidae, Strigidae, and Cathartidae. Two of the chlamydial DNA positive birds (one Swainson's hawk (Buteo swainsoni) and one red-tailed hawk (Buteo jamaicensis)) were necropsied, and tissues were collected for culture. Sequencing of the cultured elementary bodies revealed a chlamydial DNA sequence with 99.97% average nucleotide identity to the recently described Chlamydia buteonis. Spatial clusters of seropositive raptors and raptors positive for chlamydial DNA were detected in northern California. Infections were most prevalent during the winter season. Furthermore, while the proportion of raptors testing positive for chlamydial DNA was similar across age classes, seroprevalence was highest in adults. This study questions the current knowledge on C. buteonis host range and highlights the importance of further studies to evaluate the diversity and epidemiology of Chlamydia spp. infecting raptor populations.

Adapterama I: universal stubs and primers for 384 unique dual-indexed or 147,456 combinatorially-indexed Illumina libraries (iTru & iNext).

Massively parallel DNA sequencing offers many benefits, but major inhibitory cost factors include: (1) start-up (i.e., purchasing initial reagents and equipment); (2) buy-in (i.e., getting the smallest possible amount of data from a run); and (3) sample preparation. Reducing sample preparation costs is commonly addressed, but start-up and buy-in costs are rarely addressed. We present dual-indexing systems to address all three of these issues. By breaking the library construction process into universal, re-usable, combinatorial components, we reduce all costs, while increasing the number of samples and the variety of library types that can be combined within runs. We accomplish this by extending the Illumina TruSeq dual-indexing approach to 768 (384 + 384) indexed primers that produce 384 unique dual-indexes or 147,456 (384 × 384) unique combinations. We maintain eight nucleotide indexes, with many that are compatible with Illumina index sequences. We synthesized these indexing primers, purifying them with only standard desalting and placing small aliquots in replicate plates. In qPCR validation tests, 206 of 208 primers tested passed (99% success). We then created hundreds of libraries in various scenarios. Our approach reduces start-up and per-sample costs by requiring only one universal adapter that works with indexed PCR primers to uniquely identify samples. Our approach reduces buy-in costs because: (1) relatively few oligonucleotides are needed to produce a large number of indexed libraries; and (2) the large number of possible primers allows researchers to use unique primer sets for different projects, which facilitates pooling of samples during sequencing. Our libraries make use of standard Illumina sequencing primers and index sequence length and are demultiplexed with standard Illumina software, thereby minimizing customization headaches. In subsequent Adapterama papers, we use these same primers with different adapter stubs to construct amplicon and restriction-site associated DNA libraries, but their use can be expanded to any type of library sequenced on Illumina platforms.

Systems biology of the qa gene cluster in Neurospora crassa.

An ensemble of genetic networks that describe how the model fungal system, Neurospora crassa, utilizes quinic acid (QA) as a sole carbon source has been identified previously. A genetic network for QA metabolism involves the genes, qa-1F and qa-1S, that encode a transcriptional activator and repressor, respectively and structural genes, qa-2, qa-3, qa-4, qa-x, and qa-y. By a series of 4 separate and independent, model-guided, microarray experiments a total of 50 genes are identified as QA-responsive and hypothesized to be under QA-1F control and/or the control of a second QA-responsive transcription factor (NCU03643) both in the fungal binuclear Zn(II)2Cys6 cluster family. QA-1F regulation is not sufficient to explain the quantitative variation in expression profiles of the 50 QA-responsive genes. QA-responsive genes include genes with products in 8 mutually connected metabolic pathways with 7 of them one step removed from the tricarboxylic (TCA) Cycle and with 7 of them one step removed from glycolysis: (1) starch and sucrose metabolism; (2) glycolysis/glucanogenesis; (3) TCA Cycle; (4) butanoate metabolism; (5) pyruvate metabolism; (6) aromatic amino acid and QA metabolism; (7) valine, leucine, and isoleucine degradation; and (8) transport of sugars and amino acids. Gene products both in aromatic amino acid and QA metabolism and transport show an immediate response to shift to QA, while genes with products in the remaining 7 metabolic modules generally show a delayed response to shift to QA. The additional QA-responsive cutinase transcription factor-1ß (NCU03643) is found to have a delayed response to shift to QA. The series of microarray experiments are used to expand the previously identified genetic network describing the qa gene cluster to include all 50 QA-responsive genes including the second transcription factor (NCU03643). These studies illustrate new methodologies from systems biology to guide model-driven discoveries about a core metabolic network involving carbon and amino acid metabolism in N. crassa.

Systems biology of the clock in Neurospora crassa.

A model-driven discovery process, Computing Life, is used to identify an ensemble of genetic networks that describe the biological clock. A clock mechanism involving the genes white-collar-1 and white-collar-2 (wc-1 and wc-2) that encode a transcriptional activator (as well as a blue-light receptor) and an oscillator frequency (frq) that encodes a cyclin that deactivates the activator is used to guide this discovery process through three cycles of microarray experiments. Central to this discovery process is a new methodology for the rational design of a Maximally Informative Next Experiment (MINE), based on the genetic network ensemble. In each experimentation cycle, the MINE approach is used to select the most informative new experiment in order to mine for clock-controlled genes, the outputs of the clock. As much as 25% of the N. crassa transcriptome appears to be under clock-control. Clock outputs include genes with products in DNA metabolism, ribosome biogenesis in RNA metabolism, cell cycle, protein metabolism, transport, carbon metabolism, isoprenoid (including carotenoid) biosynthesis, development, and varied signaling processes. Genes under the transcription factor complex WCC ( = WC-1/WC-2) control were resolved into four classes, circadian only (612 genes), light-responsive only (396), both circadian and light-responsive (328), and neither circadian nor light-responsive (987). In each of three cycles of microarray experiments data support that wc-1 and wc-2 are auto-regulated by WCC. Among 11,000 N. crassa genes a total of 295 genes, including a large fraction of phosphatases/kinases, appear to be under the immediate control of the FRQ oscillator as validated by 4 independent microarray experiments. Ribosomal RNA processing and assembly rather than its transcription appears to be under clock control, suggesting a new mechanism for the post-transcriptional control of clock-controlled genes.

Genome-wide expression analysis of genetic networks in Neurospora crassa.

The products of five structural genes and two regulatory genes of the qa gene cluster of Neurospora crassa control the metabolism of quinic acid (QA) as a carbon source. A detailed genetic network model of this metabolic process has been reported. This investigation is designed to expand the current model of the QA reaction network. The ensemble method of network identification was used to model RNA profiling data on the qa gene cluster. Through microarray and cluster analysis, genome-wide identification of RNA transcripts associated with quinic acid metabolism in N. crassa is described and suggests a connection to other metabolic circuits. More than 100 genes whose products include carbon metabolism, protein degradation and modification, amino acid metabolism and ribosome synthesis appear to be connected to quinic acid metabolism. The core of the qa gene cluster network is validated with respect to RNA profiling data obtained from microarrays.

Class II transactivator is required for maximal expression of HLA-DOB in B cells.

HLA-DO, encoded by the HLA-DOA and HLA-DOB genes, has been shown to function as a modulator of Ag presentation. DNA microarray comparisons between B cells wild-type and mutant for the master regulator of MHC class II transcription, class II transactivator (CIITA), identified HLA-DOA and HLA-DOB as being up-regulated by CIITA. Although HLA-DOA had been shown previously to be regulated by CIITA, HLA-DOB expression was suggested to be independent of CIITA. A series of assays including quantitative RT-PCR, promoter-reporter assays, chromatin immunoprecipitations, and intracellular staining were performed to corroborate the DNA microarray analysis. The combined data demonstrate that HLA-DOB levels are increased by CIITA, and that this difference has an impact on the overall level of HLA-DO expression. Additionally, unlike the classical MHC class II genes, HLA-DOB expression was present in the absence of CIITA, indicating that additional factors mediate HLA-DOB expression in B cells.