Deep mining of the Sequence Read Archive reveals major genetic innovations in coronaviruses and other nidoviruses of aquatic vertebrates.

Virus discovery by genomics and metagenomics empowered studies of viromes, facilitated characterization of pathogen epidemiology, and redefined our understanding of the natural genetic diversity of viruses with profound functional and structural implications. Here we employed a data-driven virus discovery approach that directly queries unprocessed sequencing data in a highly parallelized way and involves a targeted viral genome assembly strategy in a wide range of sequence similarity. By screening more than 269,000 datasets of numerous authors from the Sequence Read Archive and using two metrics that quantitatively assess assembly quality, we discovered 40 nidoviruses from six virus families whose members infect vertebrate hosts. They form 13 and 32 putative viral subfamilies and genera, respectively, and include 11 coronaviruses with bisegmented genomes from fishes and amphibians, a giant 36.1 kilobase coronavirus genome with a duplicated spike glycoprotein (S) gene, 11 tobaniviruses and 17 additional corona-, arteri-, cremega-, nanhypo- and nangoshaviruses. Genome segmentation emerged in a single evolutionary event in the monophyletic lineage encompassing the subfamily Pitovirinae. We recovered the bisegmented genome sequences of two coronaviruses from RNA samples of 69 infected fishes and validated the presence of poly(A) tails at both segments using 3'RACE PCR and subsequent Sanger sequencing. We report a genetic linkage between accessory and structural proteins whose phylogenetic relationships and evolutionary distances are incongruent with the phylogeny of replicase proteins. We rationalize these observations in a model of inter-family S recombination involving at least five ancestral corona- and tobaniviruses of aquatic hosts. In support of this model, we describe an individual fish co-infected with members from the families Coronaviridae and Tobaniviridae. Our results expand the scale of the known extraordinary evolutionary plasticity in nidoviral genome architecture and call for revisiting fundamentals of genome expression, virus particle biology, host range and ecology of vertebrate nidoviruses.

Lipidomic risk scores are independent of polygenic risk scores and can predict incidence of diabetes and cardiovascular disease in a large population cohort.

Type 2 diabetes (T2D) and cardiovascular disease (CVD) represent significant disease burdens for most societies and susceptibility to these diseases is strongly influenced by diet and lifestyle. Physiological changes associated with T2D or CVD, such has high blood pressure and cholesterol and glucose levels in the blood, are often apparent prior to disease incidence. Here we integrated genetics, lipidomics, and standard clinical diagnostics to assess future T2D and CVD risk for 4,067 participants from a large prospective population-based cohort, the Malmö Diet and Cancer-Cardiovascular Cohort. By training Ridge regression-based machine learning models on the measurements obtained at baseline when the individuals were healthy, we computed several risk scores for T2D and CVD incidence during up to 23 years of follow-up. We used these scores to stratify the participants into risk groups and found that a lipidomics risk score based on the quantification of 184 plasma lipid concentrations resulted in a 168% and 84% increase of the incidence rate in the highest risk group and a 77% and 53% decrease of the incidence rate in lowest risk group for T2D and CVD, respectively, compared to the average case rates of 13.8% and 22.0%. Notably, lipidomic risk correlated only marginally with polygenic risk, indicating that the lipidome and genetic variants may constitute largely independent risk factors for T2D and CVD. Risk stratification was further improved by adding standard clinical variables to the model, resulting in a case rate of 51.0% and 53.3% in the highest risk group for T2D and CVD, respectively. The participants in the highest risk group showed significantly altered lipidome compositions affecting 167 and 157 lipid species for T2D and CVD, respectively. Our results demonstrated that a subset of individuals at high risk for developing T2D or CVD can be identified years before disease incidence. The lipidomic risk, which is derived from only one single mass spectrometric measurement that is cheap and fast, is informative and could extend traditional risk assessment based on clinical assays.

Conservation of the HBV RNA element epsilon in nackednaviruses reveals ancient origin of protein-primed reverse transcription.

Hepadnaviruses, with the human hepatitis B virus as prototype, are small, enveloped hepatotropic DNA viruses which replicate by reverse transcription of an RNA intermediate. Replication is initiated by a unique protein-priming mechanism whereby a hydroxy amino acid side chain of the terminal protein (TP) domain of the viral polymerase (P) is extended into a short DNA oligonucleotide, which subsequently serves as primer for first-strand synthesis. A key component in the priming of reverse transcription is the viral RNA element epsilon, which contains the replication origin and serves as a template for DNA primer synthesis. Here, we show that recently discovered non-enveloped fish viruses, termed nackednaviruses [C. Lauber et al., Cell Host Microbe 22, 387-399 (2017)], employ a fundamentally similar replication mechanism despite their huge phylogenetic distance and major differences in genome organization and viral lifestyle. In vitro cross-priming studies revealed that few strategic nucleotide substitutions in epsilon enable site-specific protein priming by heterologous P proteins, demonstrating that epsilon is functionally conserved since the two virus families diverged more than 400 Mya. In addition, other cis elements crucial for the hepadnavirus-typical replication of pregenomic RNA into relaxed circular double-stranded DNA were identified at conserved positions in the nackednavirus genomes. Hence, the replication mode of both hepadnaviruses and nackednaviruses was already established in their Paleozoic common ancestor, making it a truly ancient and evolutionary robust principle of genome replication that is more widespread than previously thought.

The Human Liver-Expressed Lectin CD302 Restricts Hepatitis C Virus Infection.

C-type lectin domain-containing proteins (CTLDcps) shape host responses to pathogens and infectious disease outcomes. Previously, we identified the murine CTLDcp Cd302 as restriction factor, limiting hepatitis C virus (HCV) infection of murine hepatocytes. In this study, we investigated in detail the human orthologue's ability to restrict HCV infection in human liver cells. CD302 overexpression in Huh-7.5 cells potently inhibited infection of diverse HCV chimeras representing seven genotypes. Transcriptional profiling revealed abundant CD302 mRNA expression in human hepatocytes, the natural cellular target of HCV. Knockdown of endogenously expressed CD302 modestly enhanced HCV infection of Huh-7.5 cells and primary human hepatocytes. Functional analysis of naturally occurring CD302 transcript variants and engineered CD302 mutants showed that the C-type lectin-like domain (CTLD) is essential for HCV restriction, whereas the cytoplasmic domain (CPD) is dispensable. Coding single nucleotide polymorphisms occurring in human populations and mapping to different domains of CD302 did not influence the capacity of CD302 to restrict HCV. Assessment of the anti-HCV phenotype at different life cycle stages indicated that CD302 preferentially targets the viral entry step. In contrast to the murine orthologue, overexpression of human CD302 did not modulate downstream expression of nuclear receptor-controlled genes. Ectopic CD302 expression restricted infection of liver tropic hepatitis E virus (HEV), while it did not affect infection rates of two respiratory viruses, including respiratory syncytial virus (RSV) and the alpha coronavirus HVCoV-229E. Together, these findings suggest that CD302 contributes to liver cell-intrinsic defense against HCV and might mediate broader antiviral defenses against additional hepatotropic viruses. IMPORTANCE The liver represents an immunoprivileged organ characterized by enhanced resistance to immune responses. However, the importance of liver cell-endogenous, noncytolytic innate immune responses in pathogen control is not well defined. Although the role of myeloid cell-expressed CTLDcps in host responses to viruses has been characterized in detail, we have little information about their potential functions in the liver and their relevance for immune responses in this organ. Human hepatocytes endogenously express the CTLDcp CD302. Here, we provide evidence that CD302 limits HCV infection of human liver cells, likely by inhibiting a viral cell entry step. We confirm that the dominant liver-expressed transcript variant, as well as naturally occurring coding variants of CD302, maintain the capacity to restrict HCV. We further show that the CTLD of the protein is critical for the anti-HCV activity and that overexpressed CD302 limits HEV infection. Thus, CD302 likely contributes to human liver-intrinsic antiviral defenses.

Initial HCV infection of adult hepatocytes triggers a temporally structured transcriptional program containing diverse pro- and anti-viral elements.

Transcriptional profiling provides global snapshots of virus-mediated cellular reprogramming, which can simultaneously encompass pro- and antiviral components. To determine early transcriptional signatures associated with HCV infection of authentic target cells, we performed ex vivo infections of adult primary human hepatocytes (PHHs) from seven donors. Longitudinal sampling identified minimal gene dysregulation at six hours post infection (hpi). In contrast, at 72 hpi, massive increases in the breadth and magnitude of HCV-induced gene dysregulation were apparent, affecting gene classes associated with diverse biological processes. Comparison with HCV-induced transcriptional dysregulation in Huh-7.5 cells identified limited overlap between the two systems. Of note, in PHHs, HCV infection initiated broad upregulation of canonical interferon (IFN)-mediated defense programs, limiting viral RNA replication and abrogating virion release. We further find that constitutive expression of IRF1 in PHHs maintains a steady-state antiviral program in the absence of infection, which can additionally reduce HCV RNA translation and replication. We also detected infection-induced downregulation of ∼90 genes encoding components of the EIF2 translation initiation complex and ribosomal subunits in PHHs, consistent with a signature of translational shutoff. As HCV polyprotein translation occurs independently of the EIF2 complex, this process is likely pro-viral: only translation initiation of host transcripts is arrested. The combination of antiviral intrinsic and inducible immunity, balanced against pro-viral programs, including translational arrest, maintains HCV replication at a low-level in PHHs. This may ultimately keep HCV under the radar of extra-hepatocyte immune surveillance while initial infection is established, promoting tolerance, preventing clearance and facilitating progression to chronicity.IMPORTANCEAcute HCV infections are often asymptomatic and therefore frequently undiagnosed. We endeavored to recreate this understudied phase of HCV infection using explanted PHHs and monitored host responses to initial infection. We detected temporally distinct virus-induced perturbations in the transcriptional landscape, which were initially narrow but massively amplified in breadth and magnitude over time. At 72 hpi, we detected dysregulation of diverse gene programs, concurrently promoting both virus clearance and virus persistence. On the one hand, baseline expression of IRF1 combined with infection-induced upregulation of IFN-mediated effector genes suppresses virus propagation. On the other, we detect transcriptional signatures of host translational inhibition, which likely reduces processing of IFN-regulated gene transcripts and facilitates virus survival. Together, our data provide important insights into constitutive and virus-induced transcriptional programs in PHHs, and identifies simultaneous antagonistic dysregulation of pro-and anti-viral programs which may facilitate host tolerance and promote viral persistence.

Evidence for Internal Initiation of RNA Synthesis by the Hepatitis C Virus RNA-Dependent RNA Polymerase NS5B In Cellulo.

Initiation of RNA synthesis by the hepatitis C virus (HCV) RNA-dependent RNA polymerase (RdRp) NS5B has been extensively studied in vitro and in cellulo Intracellular replication is thought to rely exclusively on terminal de novo initiation, as it conserves all genetic information of the genome. In vitro, however, additional modes of initiation have been observed. In this study, we aimed to clarify whether the intracellular environment allows for internal initiation of RNA replication by the HCV replicase. We used a dual luciferase replicon harboring a terminal and an internal copy of the viral genomic 5' untranslated region, which was anticipated to support noncanonical initiation. Indeed, a shorter RNA species was detected by Northern blotting with low frequency, depending on the length and sequence composition upstream of the internal initiation site. By introducing mutations at either site, we furthermore established that internal and terminal initiation shared identical sequence requirements. Importantly, lethal point mutations at the terminal site resulted exclusively in truncated replicons. In contrast, the same mutations at the internal site abrogated internal initiation, suggesting a competitive selection of initiation sites, rather than recombination or template-switching events. In conclusion, our data indicate that the HCV replicase is capable of internal initiation in its natural environment, although functional replication likely requires only terminal initiation. Since many other positive-strand RNA viruses generate subgenomic messenger RNAs during their replication cycle, we surmise that their capability for internal initiation is a common and conserved feature of viral RdRps.IMPORTANCE Many aspects of viral RNA replication of hepatitis C virus (HCV) are still poorly understood. The process of RNA synthesis is driven by the RNA-dependent RNA polymerase (RdRp) NS5B. Most mechanistic studies on NS5B so far were performed with in vitro systems using isolated recombinant polymerase. In this study, we present a replicon model, which allows the intracellular assessment of noncanonical modes of initiation by the full HCV replicase. Our results add to the understanding of the biochemical processes underlying initiation of RNA synthesis by NS5B by the discovery of internal initiation in cellulo Moreover, they validate observations made in vitro, showing that the viral polymerase acts very similarly in isolation and in complex with other viral and host proteins. Finally, these observations provide clues about the evolution of RdRps of positive-strand RNA viruses, which might contain the intrinsic ability to initiate internally.

Comparative analysis of histologically classified oligodendrogliomas reveals characteristic molecular differences between subgroups.

BACKGROUND: Molecular data of histologically classified oligodendrogliomas are available offering the possibility to stratify these human brain tumors into clinically relevant molecular subtypes. METHODS: Gene copy number, mutation, and expression data of 193 histologically classified oligodendrogliomas from The Cancer Genome Atlas (TCGA) were analyzed by well-established computational approaches (unsupervised clustering, statistical testing, network inference). RESULTS: We applied hierarchical clustering to tumor gene copy number profiles and revealed three molecular subgroups within histologically classified oligodendrogliomas. We further screened these subgroups for molecular glioma markers (1p/19q co-deletion, IDH mutation, gain of chromosome 7 and loss of chromosome 10) and found that our subgroups largely resemble known molecular glioma subtypes. We excluded glioblastoma-like tumors (7a10d subgroup) and derived a gene expression signature distinguishing histologically classified oligodendrogliomas with concurrent 1p/19q co-deletion and IDH mutation (1p/19q subgroup) from those with predominant IDH mutation alone (IDHme subgroup). Interestingly, many signature genes were part of signaling pathways involved in the regulation of cell proliferation, differentiation, migration, and cell-cell contacts. We further learned a gene regulatory network associated with the gene expression signature revealing novel putative major regulators with functions in cytoskeleton remodeling (e.g. APBB1IP, VAV1, ARPC1B), apoptosis (CCNL2, CREB3L1), and neural development (e.g. MYTIL, SCRT1, MEF2C) potentially contributing to the manifestation of differences between both subgroups. Moreover, we revealed characteristic expression differences of several HOX and SOX transcription factors suggesting the activity of different glioma stemness programs in both subgroups. CONCLUSIONS: We show that gene copy number profiles alone are sufficient to derive molecular subgroups of histologically classified oligodendrogliomas that are well-embedded into general glioma classification schemes. Moreover, our revealed novel putative major regulators and characteristic stemness signatures indicate that different developmental programs might be active in these subgroups, providing a basis for future studies.

ICTV Virus Taxonomy Profile: Polyomaviridae.

The Polyomaviridae is a family of small, non-enveloped viruses with circular dsDNA genomes of approximately 5 kbp. The family includes four genera whose members have restricted host range, infecting mammals and birds. Polyomavirus genomes have also been detected recently in fish. Merkel cell polyomavirus and raccoon polyomavirus are associated with cancer in their host; other members are human and veterinary pathogens. Clinical manifestations are obvious in immunocompromised patients but not in healthy individuals. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the taxonomy of the Polyomaviridae, which is available at www.ictv.global/report/polyomaviridae.

Interspecific adaptation by binary choice at de novo polyomavirus T antigen site through accelerated codon-constrained Val-Ala toggling within an intrinsically disordered region.

It is common knowledge that conserved residues evolve slowly. We challenge generality of this central tenet of molecular biology by describing the fast evolution of a conserved nucleotide position that is located in the overlap of two open reading frames (ORFs) of polyomaviruses. The de novo ORF is expressed through either the ALTO protein or the Middle T antigen (MT/ALTO), while the ancestral ORF encodes the N-terminal domain of helicase-containing Large T (LT) antigen. In the latter domain the conserved Cys codon of the LXCXE pRB-binding motif constrains codon evolution in the overlapping MT/ALTO ORF to a binary choice between Val and Ala codons, termed here as codon-constrained Val-Ala (COCO-VA) toggling. We found the rate of COCO-VA toggling to approach the speciation rate and to be significantly accelerated compared to the baseline rate of chance substitution in a large monophyletic lineage including all viruses encoding MT/ALTO and three others. Importantly, the COCO-VA site is located in a short linear motif (SLiM) of an intrinsically disordered region, a typical characteristic of adaptive responders. These findings provide evidence that the COCO-VA toggling is under positive selection in many polyomaviruses, implying its critical role in interspecific adaptation, which is unprecedented for conserved residues.

The footprint of genome architecture in the largest genome expansion in RNA viruses.

The small size of RNA virus genomes (2-to-32 kb) has been attributed to high mutation rates during replication, which is thought to lack proof-reading. This paradigm is being revisited owing to the discovery of a 3'-to-5' exoribonuclease (ExoN) in nidoviruses, a monophyletic group of positive-stranded RNA viruses with a conserved genome architecture. ExoN, a homolog of canonical DNA proof-reading enzymes, is exclusively encoded by nidoviruses with genomes larger than 20 kb. All other known non-segmented RNA viruses have smaller genomes. Here we use evolutionary analyses to show that the two- to three-fold expansion of the nidovirus genome was accompanied by a large number of replacements in conserved proteins at a scale comparable to that in the Tree of Life. To unravel common evolutionary patterns in such genetically diverse viruses, we established the relation between genomic regions in nidoviruses in a sequence alignment-free manner. We exploited the conservation of the genome architecture to partition each genome into five non-overlapping regions: 5' untranslated region (UTR), open reading frame (ORF) 1a, ORF1b, 3'ORFs (encompassing the 3'-proximal ORFs), and 3' UTR. Each region was analyzed for its contribution to genome size change under different models. The non-linear model statistically outperformed the linear one and captured >92% of data variation. Accordingly, nidovirus genomes were concluded to have reached different points on an expansion trajectory dominated by consecutive increases of ORF1b, ORF1a, and 3'ORFs. Our findings indicate a unidirectional hierarchical relation between these genome regions, which are distinguished by their expression mechanism. In contrast, these regions cooperate bi-directionally on a functional level in the virus life cycle, in which they predominantly control genome replication, genome expression, and virus dissemination, respectively. Collectively, our findings suggest that genome architecture and the associated region-specific division of labor leave a footprint on genome expansion and may limit RNA genome size.

Characterization of hepatitis C virus intra- and intergenotypic chimeras reveals a role of the glycoproteins in virus envelopment.

Hepatitis C virus (HCV) is highly variable and associated with chronic liver disease. Viral isolates are grouped into seven genotypes (GTs). Accumulating evidence indicates that viral determinants in the core to NS2 proteins modulate the efficiency of virus production. However, the role of the glycoproteins E1 and E2 in this process is currently poorly defined. Therefore, we constructed chimeric viral genomes to explore the role of E1 and E2 in HCV assembly. Comparison of the kinetics and efficiency of particle production by intragenotypic chimeras highlighted core and p7 as crucial determinants for efficient virion release. Glycoprotein sequences, however, had only a minimal impact on this process. In contrast, in the context of intergenotypic HCV chimeras, HCV assembly was profoundly influenced by glycoprotein genes. On the one hand, insertion of GT1a-derived (H77) E1-E2 sequences into a chimeric GT2a virus (Jc1) strongly suppressed virus production. On the other hand, replacement of H77 glycoproteins within the GT1a-GT2a chimeric genome H77/C3 by GT2a-derived (Jc1) E1-E2 increased infectious particle production. Thus, within intergenotypic chimeras, glycoprotein features strongly modulate virus production. Replacement of Jc1 glycoprotein genes by H77-derived E1-E2 did not grossly affect subcellular localization of core, E2, and NS2. However, it caused an accumulation of nonenveloped core protein and increased abundance of nonenveloped core protein structures with slow sedimentation. These findings reveal an important role for the HCV glycoproteins E1 and E2 in membrane envelopment, which likely depends on a genotype-specific interplay with additional viral factors.

Major Histocompatibility Complex class IIb polymorphism influences gut microbiota composition and diversity.

Animals harbour diverse communities of symbiotic bacteria, which differ dramatically among host individuals. This heterogeneity poses an immunological challenge: distinguishing between mutualistic and pathogenic members of diverse and host-specific microbial communities. We propose that Major Histocompatibility class II (MHC) genotypes contribute to recognition and regulation of gut microbes, and thus, MHC polymorphism contributes to microbial variation among hosts. Here, we show that MHC IIb polymorphism is associated with among-individual variation in gut microbiota within a single wild vertebrate population of a small fish, the threespine stickleback. We sampled stickleback from Cedar Lake, on Vancouver Island, and used next-generation sequencing to genotype the sticklebacks' gut microbiota (16S sequencing) and their MHC class IIb exon 2 sequences. The presence of certain MHC motifs was associated with altered relative abundance (increase or decrease) of some microbial Families. The effect sizes are modest and entail a minority of microbial taxa, but these results represent the first indication that MHC genotype may affect gut microbiota composition in natural populations (MHC-microbe associations have also been found in a few studies of lab mice). Surprisingly, these MHC effects were frequently sex-dependent. Finally, hosts with more diverse MHC motifs had less diverse gut microbiota. One implication is that MHC might influence the efficacy of therapeutic strategies to treat dysbiosis-associated disease, including the outcome of microbial transplants between healthy and diseased patients. We also speculate that macroparasite-driven selection on MHC has the potential to indirectly alter the host gut microbiota, and vice versa.

Analysis of Replication, Cell Division-Mediated Spread, and HBV Envelope Protein-Dependent Pseudotyping of Three Mammalian Delta-like Agents.

The human hepatitis delta virus (HDV) is a satellite RNA virus that depends on hepatitis B virus (HBV) surface proteins (HBsAg) to assemble into infectious virions targeting the same organ (liver) as HBV. Until recently, the evolutionary origin of HDV remained largely unknown. The application of bioinformatics on whole sequence databases lead to discoveries of HDV-like agents (DLA) and shed light on HDV's evolution, expanding our understanding of HDV biology. DLA were identified in heterogeneous groups of vertebrates and invertebrates, highlighting that the evolution of HDV, represented by eight distinct genotypes, is broader and more complex than previously foreseen. In this study, we focused on the characterization of three mammalian DLA discovered in woodchuck (Marmota monax), white-tailed deer (Odocoileus virginianus), and lesser dog-like bat (Peropteryx macrotis) in terms of replication, cell-type permissiveness, and spreading pathways. We generated replication-competent constructs expressing 1.1-fold over-length antigenomic RNA of each DLA. Replication was initiated by transfecting the cDNAs into human (HuH7, HeLa, HEK293T, A549) and non-human (Vero E6, CHO, PaKi, LMH) cell lines. Upon transfection and replication establishment, none of the DLA expressed a large delta antigen. A cell division-mediated viral amplification assay demonstrated the capability of non-human DLA to replicate and propagate in hepatic and non-hepatic tissues, without the requirement of envelope proteins from a helper virus. Remarkably L-HDAg but not S-HDAg from HDV can artificially mediate envelopment of WoDV and DeDV ribonucleoproteins (RNPs) by HBsAg to form infectious particles, as demonstrated by co-transfection of HuH7 cells with the respective DLA expression constructs and a plasmid encoding HBV envelope proteins. These chimeric viruses are sensitive to HDV entry inhibitors and allow synchronized infections for comparative replication studies. Our results provide a more detailed understanding of the molecular biology, evolution, and virus-host interaction of this unique group of animal viroid-like agents in relation to HDV.

Drug repurposing screen identifies lonafarnib as respiratory syncytial virus fusion protein inhibitor.

Respiratory syncytial virus (RSV) is a common cause of acute lower respiratory tract infection in infants, older adults and the immunocompromised. Effective directly acting antivirals are not yet available for clinical use. To address this, we screen the ReFRAME drug-repurposing library consisting of 12,000 small molecules against RSV. We identify 21 primary candidates including RSV F and N protein inhibitors, five HSP90 and four IMPDH inhibitors. We select lonafarnib, a licensed farnesyltransferase inhibitor, and phase III candidate for hepatitis delta virus (HDV) therapy, for further follow-up. Dose-response analyses and plaque assays confirm the antiviral activity (IC50: 10-118 nM). Passaging of RSV with lonafarnib selects for phenotypic resistance and fixation of mutations in the RSV fusion protein (T335I and T400A). Lentiviral pseudotypes programmed with variant RSV fusion proteins confirm that lonafarnib inhibits RSV cell entry and that these mutations confer lonafarnib resistance. Surface plasmon resonance reveals RSV fusion protein binding of lonafarnib and co-crystallography identifies the lonafarnib binding site within RSV F. Oral administration of lonafarnib dose-dependently reduces RSV virus load in a murine infection model using female mice. Collectively, this work provides an overview of RSV drug repurposing candidates and establishes lonafarnib as a bona fide fusion protein inhibitor.

From Stockholm to Malawi: recent developments in studying human polyomaviruses.

Until a few years ago the polyomavirus family (Polyomaviridae) included a dozen viruses identified in avian and mammalian hosts. Two of these, the JC and BK-polyomaviruses isolated a long time ago, are known to infect humans and cause severe illness in immunocompromised hosts. Since 2007 an unprecedented number of eight novel polyomaviruses were discovered in humans. Among them are the KI- and WU-polyomaviruses identified in respiratory samples, the Merkel cell polyomavirus found in skin carcinomas and the polyomavirus associated with trichodysplasia spinulosa, a skin disease of transplant patients. Another four novel human polyomaviruses were identified, HPyV6, HPyV7, HPyV9 and the Malawi polyomavirus, so far not associated with any disease. In the same period several novel mammalian polyomaviruses were described. This review summarizes the recent developments in studying the novel human polyomaviruses, and touches upon several aspects of polyomavirus virology, pathogenicity, epidemiology and phylogeny.

Partitioning the genetic diversity of a virus family: approach and evaluation through a case study of picornaviruses.

The recent advent of genome sequences as the only source available to classify many newly discovered viruses challenges the development of virus taxonomy by expert virologists who traditionally rely on extensive virus characterization. In this proof-of-principle study, we address this issue by presenting a computational approach (DEmARC) to classify viruses of a family into groups at hierarchical levels using a sole criterion-intervirus genetic divergence. To quantify genetic divergence, we used pairwise evolutionary distances (PEDs) estimated by maximum likelihood inference on a multiple alignment of family-wide conserved proteins. PEDs were calculated for all virus pairs, and the resulting distribution was modeled via a mixture of probability density functions. The model enables the quantitative inference of regions of distance discontinuity in the family-wide PED distribution, which define the levels of hierarchy. For each level, a limit on genetic divergence, below which two viruses join the same group, was objectively selected among a set of candidates by minimizing violations of intragroup PEDs to the limit. In a case study, we applied the procedure to hundreds of genome sequences of picornaviruses and extensively evaluated it by modulating four key parameters. It was found that the genetics-based classification largely tolerates variations in virus sampling and multiple alignment construction but is affected by the choice of protein and the measure of genetic divergence. In an accompanying paper (C. Lauber and A. E. Gorbalenya, J. Virol. 86:3905-3915, 2012), we analyze the substantial insight gained with the genetics-based classification approach by comparing it with the expert-based picornavirus taxonomy.

Toward genetics-based virus taxonomy: comparative analysis of a genetics-based classification and the taxonomy of picornaviruses.

Virus taxonomy has received little attention from the research community despite its broad relevance. In an accompanying paper (C. Lauber and A. E. Gorbalenya, J. Virol. 86:3890-3904, 2012), we have introduced a quantitative approach to hierarchically classify viruses of a family using pairwise evolutionary distances (PEDs) as a measure of genetic divergence. When applied to the six most conserved proteins of the Picornaviridae, it clustered 1,234 genome sequences in groups at three hierarchical levels (to which we refer as the "GENETIC classification"). In this study, we compare the GENETIC classification with the expert-based picornavirus taxonomy and outline differences in the underlying frameworks regarding the relation of virus groups and genetic diversity that represent, respectively, the structure and content of a classification. To facilitate the analysis, we introduce two novel diagrams. The first connects the genetic diversity of taxa to both the PED distribution and the phylogeny of picornaviruses. The second depicts a classification and the accommodated genetic diversity in a standardized manner. Generally, we found striking agreement between the two classifications on species and genus taxa. A few disagreements concern the species Human rhinovirus A and Human rhinovirus C and the genus Aphthovirus, which were split in the GENETIC classification. Furthermore, we propose a new supergenus level and universal, level-specific PED thresholds, not reached yet by many taxa. Since the species threshold is approached mostly by taxa with large sampling sizes and those infecting multiple hosts, it may represent an upper limit on divergence, beyond which homologous recombination in the six most conserved genes between two picornaviruses might not give viable progeny.

Discovery of the first insect nidovirus, a missing evolutionary link in the emergence of the largest RNA virus genomes.

Nidoviruses with large genomes (26.3-31.7 kb; 'large nidoviruses'), including Coronaviridae and Roniviridae, are the most complex positive-sense single-stranded RNA (ssRNA+) viruses. Based on genome size, they are far separated from all other ssRNA+ viruses (below 19.6 kb), including the distantly related Arteriviridae (12.7-15.7 kb; 'small nidoviruses'). Exceptionally for ssRNA+ viruses, large nidoviruses encode a 3'-5'exoribonuclease (ExoN) that was implicated in controlling RNA replication fidelity. Its acquisition may have given rise to the ancestor of large nidoviruses, a hypothesis for which we here provide evolutionary support using comparative genomics involving the newly discovered first insect-borne nidovirus. This Nam Dinh virus (NDiV), named after a Vietnamese province, was isolated from mosquitoes and is yet to be linked to any pathology. The genome of this enveloped 60-80 nm virus is 20,192 nt and has a nidovirus-like polycistronic organization including two large, partially overlapping open reading frames (ORF) 1a and 1b followed by several smaller 3'-proximal ORFs. Peptide sequencing assigned three virion proteins to ORFs 2a, 2b, and 3, which are expressed from two 3'-coterminal subgenomic RNAs. The NDiV ORF1a/ORF1b frameshifting signal and various replicative proteins were tentatively mapped to canonical positions in the nidovirus genome. They include six nidovirus-wide conserved replicase domains, as well as the ExoN and 2'-O-methyltransferase that are specific to large nidoviruses. NDiV ORF1b also encodes a putative N7-methyltransferase, identified in a subset of large nidoviruses, but not the uridylate-specific endonuclease that - in deviation from the current paradigm - is present exclusively in the currently known vertebrate nidoviruses. Rooted phylogenetic inference by Bayesian and Maximum Likelihood methods indicates that NDiV clusters with roniviruses and that its branch diverged from large nidoviruses early after they split from small nidoviruses. Together these characteristics identify NDiV as the prototype of a new nidovirus family and a missing link in the transition from small to large nidoviruses.

Viroid-like RNA-dependent RNA polymerase-encoding ambiviruses are abundant in complex fungi.

Ambiviruses are hybrid infectious elements encoding the hallmark gene of RNA viruses, the RNA-dependent RNA polymerase, and self-cleaving RNA ribozymes found in many viroids. Ambiviruses are thought to be pathogens of fungi, although the majority of reported genomes have been identified in metatranscriptomes. Here, we present a comprehensive screen for ambiviruses in more than 46,500 fungal transcriptomes from the Sequence Read Archive (SRA). Our data-driven virus discovery approach identified more than 2,500 ambiviral sequences across the kingdom Fungi with a striking expansion in members of the phylum Basidiomycota representing the most complex fungal organisms. Our study unveils a large diversity of unknown ambiviruses with as little as 27% protein sequence identity to known members and sheds new light on the evolution of this distinct class of infectious agents with RNA genomes. No evidence for the presence of ambiviruses in human microbiomes was obtained from a comprehensive screen of respective metatranscriptomes available in the SRA.

Evolutionary Insight into the Association between New Jersey Polyomavirus and Humans.

Advances in viral discovery techniques have led to the identification of numerous novel viruses in human samples. However, the low prevalence of certain viruses in humans raises doubts about their association with our species. To ascertain the authenticity of a virus as a genuine human-infecting agent, it can be useful to investigate the diversification of its lineage within hominines, the group encompassing humans and African great apes. Building upon this rationale, we examined the case of the New Jersey polyomavirus (NJPyV; Alphapolyomavirus terdecihominis), which has only been detected in a single patient thus far. In this study, we obtained and analyzed sequences from closely related viruses infecting all African great ape species. We show that NJPyV nests within the diversity of these viruses and that its lineage placement is compatible with an ancient origin in humans, despite its apparent rarity in human populations.