Using genomic epidemiology of SARS-CoV-2 to support contact tracing and public health surveillance in rural Humboldt County, California.

BACKGROUND: During the COVID-19 pandemic within the United States, much of the responsibility for diagnostic testing and epidemiologic response has relied on the action of county-level departments of public health. Here we describe the integration of genomic surveillance into epidemiologic response within Humboldt County, a rural county in northwest California. METHODS: Through a collaborative effort, 853 whole SARS-CoV-2 genomes were generated, representing ~58% of the 1,449 SARS-CoV-2-positive cases detected in Humboldt County as of March 12, 2021. Phylogenetic analysis of these data was used to develop a comprehensive understanding of SARS-CoV-2 introductions to the county and to support contact tracing and epidemiologic investigations of all large outbreaks in the county. RESULTS: In the case of an outbreak on a commercial farm, viral genomic data were used to validate reported epidemiologic links and link additional cases within the community who did not report a farm exposure to the outbreak. During a separate outbreak within a skilled nursing facility, genomic surveillance data were used to rule out the putative index case, detect the emergence of an independent Spike:N501Y substitution, and verify that the outbreak had been brought under control. CONCLUSIONS: These use cases demonstrate how developing genomic surveillance capacity within local public health departments can support timely and responsive deployment of genomic epidemiology for surveillance and outbreak response based on local needs and priorities.

Single mosquito metatranscriptomics identifies vectors, emerging pathogens and reservoirs in one assay.

Mosquitoes are major infectious disease-carrying vectors. Assessment of current and future risks associated with the mosquito population requires knowledge of the full repertoire of pathogens they carry, including novel viruses, as well as their blood meal sources. Unbiased metatranscriptomic sequencing of individual mosquitoes offers a straightforward, rapid, and quantitative means to acquire this information. Here, we profile 148 diverse wild-caught mosquitoes collected in California and detect sequences from eukaryotes, prokaryotes, 24 known and 46 novel viral species. Importantly, sequencing individuals greatly enhanced the value of the biological information obtained. It allowed us to (a) speciate host mosquito, (b) compute the prevalence of each microbe and recognize a high frequency of viral co-infections, (c) associate animal pathogens with specific blood meal sources, and (d) apply simple co-occurrence methods to recover previously undetected components of highly prevalent segmented viruses. In the context of emerging diseases, where knowledge about vectors, pathogens, and reservoirs is lacking, the approaches described here can provide actionable information for public health surveillance and intervention decisions.

Clinical features, diagnostics, and outcomes of patients presenting with acute respiratory illness: a comparison of patients with and without COVID-19.

BACKGROUND: Emerging data on the clinical presentation, diagnostics, and outcomes of patients with COVID-19 have largely been presented as case series. Few studies have compared these clinical features and outcomes of COVID-19 to other acute respiratory illnesses. METHODS: We examined all patients presenting to an emergency department in San Francisco, California between February 3 and March 31, 2020 with an acute respiratory illness who were tested for SARS-CoV-2. We determined COVID-19 status by PCR and metagenomic next generation sequencing (mNGS). We compared demographics, comorbidities, symptoms, vital signs, and laboratory results including viral diagnostics using PCR and mNGS. Among those hospitalized, we determined differences in treatment (antibiotics, antivirals, respiratory support) and outcomes (ICU admission, ICU interventions, acute respiratory distress syndrome, cardiac injury). FINDINGS: In a cohort of 316 patients, 33 (10%) tested positive for SARS-CoV-2; 31 patients, all without COVID-19, tested positive for another respiratory virus (16%). Among patients with additional viral testing, no co-infections with SARS-CoV-2 were identified by PCR or mNGS. Patients with COVID-19 reported longer symptoms duration (median 7 vs. 3 days), and were more likely to report fever (82% vs. 44%), fatigue (85% vs. 50%), and myalgias (61% vs 27%); p<0.001 for all comparisons. Lymphopenia (55% vs 34%, p=0.018) and bilateral opacities on initial chest radiograph (55% vs. 24%, p=0.001) were more common in patients with COVID-19. Patients with COVID-19 were more often hospitalized (79% vs. 56%, p=0.014). Of 186 hospitalized patients, patients with COVID-19 had longer hospitalizations (median 10.7d vs. 4.7d, p<0.001) and were more likely to develop ARDS (23% vs. 3%, p<0.001). Most comorbidities, home medications, signs and symptoms, vital signs, laboratory results, treatment, and outcomes did not differ by COVID-19 status. INTERPRETATION: While we found differences in clinical features of COVID-19 compared to other acute respiratory illnesses, there was significant overlap in presentation and comorbidities. Patients with COVID-19 were more likely to be admitted to the hospital, have longer hospitalizations and develop ARDS, and were unlikely to have co-existent viral infections. These findings enhance understanding of the clinical characteristics of COVID-19 in comparison to other acute respiratory illnesses. .

A Toolkit for Rapid Modular Construction of Biological Circuits in Mammalian Cells.

The ability to rapidly assemble and prototype cellular circuits is vital for biological research and its applications in biotechnology and medicine. Current methods for the assembly of mammalian DNA circuits are laborious, slow, and expensive. Here we present the Mammalian ToolKit (MTK), a Golden Gate-based cloning toolkit for fast, reproducible, and versatile assembly of large DNA vectors and their implementation in mammalian models. The MTK consists of a curated library of characterized, modular parts that can be assembled into transcriptional units and further weaved into complex circuits. We showcase the capabilities of the MTK by using it to generate single-integration landing pads, create and deliver libraries of protein variants and sgRNAs, and iterate through dCas9-based prototype circuits. As a biological proof of concept, we demonstrate how the MTK can speed the generation of noninfectious viral circuits to enable rapid testing of pharmacological inhibitors of emerging viruses that pose a major threat to human health.

JC polyoma viruria associates with protection from chronic kidney disease independently from apolipoprotein L1 genotype in African Americans.

Background: Viral infections can trigger chronic kidney disease (CKD) and the urine virome may inform risk. The Natural History of APOL1-Associated Nephropathy Study (NHAANS) reported that urine JC polyomavirus (JCPyV) associated with a lower risk of APOL1-associated nephropathy in African Americans. Herein, association was assessed between urine JCPyV with CKD in African Americans independent from the APOL1 genotype. Methods: Quantitative polymerase chain reaction was performed for urinary detection of JCPyV and BK polyoma virus (BKPyV) in 200 newly recruited nondiabetic African Americans. A combined analysis was performed in these individuals plus 300 NHAANS participants. Results: In the 200 new participants, urine JCPyV was present in 8.8% of CKD cases and 45.8% of nonnephropathy controls (P = 3.0 × 10-8). In those with APOL1 renal-risk genotypes, JCPyV was detected in 5.1% of cases and 40.0% of controls (P = 0.0002). In those lacking APOL1 renal-risk genotypes, JCPyV was detected in 12.2% of cases and 48.8% of controls (P = 8.5 × 10-5). BKPyV was detected in 1.3% of cases and 0.8% of controls (P = 0.77). In a combined analysis with 300 NHAANS participants (n = 500), individuals with urine JCPyV had a 63% lower risk of CKD compared with those without urine JCPyV (odds ratio 0.37; P = 4.6 × 10-6). RNA fluorescence in situ hybridization confirmed the presence of JCPyV genomic DNA and JCPyV messenger RNA (mRNA) in nondiseased kidney. Conclusions: Inverse relationships exist between JCPyV viruria and non-diabetic CKD. Future studies should determine whether renal inflammation associated with CKD is less permissive for JCPyV reactivation/replication or whether JCPyV is a marker of reduced host immune responsiveness that diminishes immune pathologic contributions to CKD.

Identification of a previously undescribed divergent virus from the Flaviviridae family in an outbreak of equine serum hepatitis.

Theiler's disease is an acute hepatitis in horses that is associated with the administration of equine blood products; its etiologic agent has remained unknown for nearly a century. Here, we used massively parallel sequencing to explore samples from a recent Theiler's disease outbreak. Metatranscriptomic analysis of the short sequence reads identified a 10.5-kb sequence from a previously undescribed virus of the Flaviviridae family, which we designate "Theiler's disease-associated virus" (TDAV). Phylogenetic analysis clusters TDAV with GB viruses of the recently proposed Pegivirus genus, although it shares only 35.3% amino acid identity with its closest relative, GB virus D. An epidemiological survey of additional horses from three separate locations supports an association between TDAV infection and acute serum hepatitis. Experimental inoculation of horses with TDAV-positive plasma provides evidence that several weeks of viremia preceded liver injury and that liver disease may not be directly related to the level of viremia. Like hepatitis C virus, the best characterized Flaviviridae species known to cause hepatitis, we find TDAV is capable of efficient parenteral transmission, engendering acute and chronic infections associated with a diversity of clinical presentations ranging from subclinical infection to clinical hepatitis.

Identification, characterization, and in vitro culture of highly divergent arenaviruses from boa constrictors and annulated tree boas: candidate etiological agents for snake inclusion body disease.

UNLABELLED: Inclusion body disease (IBD) is an infectious fatal disease of snakes typified by behavioral abnormalities, wasting, and secondary infections. At a histopathological level, the disease is identified by the presence of large eosinophilic cytoplasmic inclusions in multiple tissues. To date, no virus or other pathogen has been definitively characterized or associated with the disease. Using a metagenomic approach to search for candidate etiologic agents in snakes with confirmed IBD, we identified and de novo assembled the complete genomic sequences of two viruses related to arenaviruses, and a third arenavirus-like sequence was discovered by screening an additional set of samples. A continuous boa constrictor cell line was established and used to propagate and isolate one of the viruses in culture. Viral nucleoprotein was localized and concentrated within large cytoplasmic inclusions in infected cells in culture and tissues from diseased snakes. In total, viral RNA was detected in 6/8 confirmed IBD cases and 0/18 controls. These viruses have a typical arenavirus genome organization but are highly divergent, belonging to a lineage separate from that of the Old and New World arenaviruses. Furthermore, these viruses encode envelope glycoproteins that are more similar to those of filoviruses than to those of other arenaviruses. These findings implicate these viruses as candidate etiologic agents of IBD. The presence of arenaviruses outside mammals reveals that these viruses infect an unexpectedly broad range of species and represent a new reservoir of potential human pathogens. IMPORTANCE: Inclusion body disease (IBD) is a common infectious disease of captive snakes. IBD is fatal and can cause the loss of entire animal collections. The cause of the disease has remained elusive, and no treatment exists. In addition to being important to pet owners, veterinarians, breeders, zoological parks, and aquariums, the study of animal disease is significant since animals are the source of virtually every emerging infectious human disease. We searched for candidate causative agents in snakes diagnosed with IBD and found a group of novel viruses distantly related mainly to arenaviruses but also to filoviruses, both of which can cause fatal hemorrhagic fevers when transmitted from animals to humans. In addition to providing evidence that strongly suggests that these viruses cause snake IBD, this discovery reveals a new and unanticipated domain of virus biology and evolution.

Avian bornavirus is present in many tissues of psittacine birds with histopathologic evidence of proventricular dilatation disease.

Proventricular dilatation disease (PDD) is a neurologic disease of psittacine birds suspected to be caused by a recently identified Avian bornavirus (ABV). In the current report, data supporting the causal association of ABV with PDD are presented. Immunohistochemistry (IHC) with rabbit polyclonal antiserum raised against ABV nucleocapsid protein was used to identify cell and organ distribution of viral antigen. The ABV antigen was most consistently detected in brain, spinal cord, adrenal gland, pancreas, and kidney. Histopathologic evaluation was correlated with ABV-specific polymerase chain reaction (PCR) and immunohistochemical tests in multiple tissues from 16 psittacine birds with and without PDD. Using histopathologic diagnosis as the gold standard, the sensitivity and specificity of IHC for ABV antigens were found to be 100% and 100%, respectively. Many more tissues were positive for ABV RNA by reverse transcription PCR than were positive for pathologic changes or viral antigens by IHC, indicating the presence of subclinical or asymptomatic infection at many sites.

Analysis of naturally occurring avian bornavirus infection and transmission during an outbreak of proventricular dilatation disease among captive psittacine birds.

A proventricular dilatation disease (PDD) outbreak provided the opportunity to investigate the transmissibility of avian Bornavirus (ABV) and its linkage to PDD under natural conditions. Upon exposure to a bird with a fatal case of PDD, 10 birds became symptomatic and died. ABV2 RNA was recovered from available tissues. Further screening revealed that 12/46 exposed birds were ABV2(+). Three chicks boarded at this aviary developed PDD. They harbored the same ABV2 isolate and transmitted it to five of eight chicks in their home aviary. These findings demonstrate that ABV infection precedes the development of PDD. ABV-specific Western blotting and reverse transcription-PCR indicate that ABV2 is not strictly neurotropic.

Experimental induction of proventricular dilatation disease in cockatiels (Nymphicus hollandicus) inoculated with brain homogenates containing avian bornavirus 4.

BACKGROUND: Proventricular dilatation disease (PDD) is a fatal disorder of psittacine birds worldwide. The disease is characterized by lymphoplasmacytic infiltration of the central and peripheral nervous systems, leading to gastrointestinal motility and/or central nervous system dysfunction. Recently, we detected a significant association between avian bornavirus (ABV) infection and clinical signs of PDD in psittacines. However, it remains unclear whether ABV infection actually causes PDD. To address this question, we examined the impact of ABV inoculation on the cockatiel (Nymphicus hollandicus). RESULTS: Five cockatiels were inoculated via multiple routes (intramuscular, intraocular, intranasal, and oral) with a brain homogenate derived from either a PDD(+) avian bornavirus 4 (ABV4) (+) case (n = 3 inoculees) or from a PDD(-) ABV(-) control (n = 2 inoculees). The control birds remained free of clinical or pathological signs of PDD, and tested ABV(-) by RT-PCR and immunohistochemistry (IHC). In contrast, all three cockatiels inoculated with ABV4(+) brain homogenate developed gross and microscopic PDD lesions, and two exhibited overt clinical signs. In numerous tissues, ABV RT-PCR and sequence analysis demonstrated the presence of ABV4 RNA nearly identical to that in the inoculum. ABV was detected in the central nervous system of the three ABV-inoculees by IHC. Pyrosequencing to investigate the viral flora in the ABV4(+) inoculum uncovered 7 unique reads sharing 73-100% nucleotide sequence identity with previously identified ABV sequences and 24 reads sharing 40-89% amino acid sequence identity with viruses in the Retroviridae and Astroviridae families. Of these candidate viral species, only ABV RNA was recovered from tissues of the inoculated birds. CONCLUSION: In this study, the clinical and pathological manifestations of PDD were induced by inoculation of cockatiels with brain homogenates containing avian bornavirus 4. By using high throughput pyrosequencing an in-depth view of the viral content of the inoculum was achieved, revealing that of 3 candidate virus families detected, only the presence of ABV RNA correlated with the development of PDD. This study provides evidence of a causal association between ABV4 infection and PDD in cockatiels.

Recovery of divergent avian bornaviruses from cases of proventricular dilatation disease: identification of a candidate etiologic agent.

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.

Genome-wide diversity and selective pressure in the human rhinovirus.

BACKGROUND: The human rhinoviruses (HRV) are one of the most common and diverse respiratory pathogens of humans. Over 100 distinct HRV serotypes are known, yet only 6 genomes are available. Due to the paucity of HRV genome sequence, little is known about the genetic diversity within HRV or the forces driving this diversity. Previous comparative genome sequence analyses indicate that recombination drives diversification in multiple genera of the picornavirus family, yet it remains unclear if this holds for HRV. RESULTS: To resolve this and gain insight into the forces driving diversification in HRV, we generated a representative set of 34 fully sequenced HRVs. Analysis of these genomes shows consistent phylogenies across the genome, conserved non-coding elements, and only limited recombination. However, spikes of genetic diversity at both the nucleotide and amino acid level are detectable within every locus of the genome. Despite this, the HRV genome as a whole is under purifying selective pressure, with islands of diversifying pressure in the VP1, VP2, and VP3 structural genes and two non-structural genes, the 3C protease and 3D polymerase. Mapping diversifying residues in these factors onto available 3-dimensional structures revealed the diversifying capsid residues partition to the external surface of the viral particle in statistically significant proximity to antigenic sites. Diversifying pressure in the pleconaril binding site is confined to a single residue known to confer drug resistance (VP1 191). In contrast, diversifying pressure in the non-structural genes is less clear, mapping both nearby and beyond characterized functional domains of these factors. CONCLUSION: This work provides a foundation for understanding HRV genetic diversity and insight into the underlying biology driving evolution in HRV. It expands our knowledge of the genome sequence space that HRV reference serotypes occupy and how the pattern of genetic diversity across HRV genomes differs from other picornaviruses. It also reveals evidence of diversifying selective pressure in both structural genes known to interact with the host immune system and in domains of unassigned function in the non-structural 3C and 3D genes, raising the possibility that diversification of undiscovered functions in these essential factors may influence HRV fitness and evolution.

Rapid DNA mapping by fluorescent single molecule detection.

DNA mapping is an important analytical tool in genomic sequencing, medical diagnostics and pathogen identification. Here we report an optical DNA mapping strategy based on direct imaging of individual DNA molecules and localization of multiple sequence motifs on the molecules. Individual genomic DNA molecules were labeled with fluorescent dyes at specific sequence motifs by the action of nicking endonuclease followed by the incorporation of dye terminators with DNA polymerase. The labeled DNA molecules were then stretched into linear form on a modified glass surface and imaged using total internal reflection fluorescence (TIRF) microscopy. By determining the positions of the fluorescent labels with respect to the DNA backbone, the distribution of the sequence motif recognized by the nicking endonuclease can be established with good accuracy, in a manner similar to reading a barcode. With this approach, we constructed a specific sequence motif map of lambda-DNA. We further demonstrated the capability of this approach to rapidly type a human adenovirus and several strains of human rhinovirus.

E-Predict: a computational strategy for species identification based on observed DNA microarray hybridization patterns.

DNA microarrays may be used to identify microbial species present in environmental and clinical samples. However, automated tools for reliable species identification based on observed microarray hybridization patterns are lacking. We present an algorithm, E-Predict, for microarray-based species identification. E-Predict compares observed hybridization patterns with theoretical energy profiles representing different species. We demonstrate the application of the algorithm to viral detection in a set of clinical samples and discuss its relevance to other metagenomic applications.