Down on the farm: Farm coronaviruses are a large reservoir for new spillover events into humans

Until COVID‐19, coronaviruses were largely overlooked by virologists. Yet, their abundance in mammalian species could cause new spillover events into humans. Subject Categories: Ecology, Microbiology, Virology & Host Pathogen Interaction

By now, coronaviruses have become notorious. Yet, before the SARS coronaviruses started wreaking havoc with human society, coronaviruses were largely overlooked, even amongst virologists. “When I told another virologist at a meeting what I was working on”, recalled Peter Rottier at Utrecht University, the Netherlands, who began his coronavirus research in 1979, “that person said they''d never heard of coronaviruses”. With plenty of other viruses such as influenza and HIV demanding attention, funders saw no reason to support research on the endemic human coronaviruses in the 1980s and 1990s since they caused mostly common cold.

With plenty of other viruses such as influenza and HIV demanding attention, funders saw no reason to support research on the endemic human coronaviruses in the 1980s and 90s…

Clinicians and virologists began taking much notice around 2003 and 2004 after the SARS coronavirus inflicted about 800 deaths before the outbreak receded, but the potential of this virus family to spillover into new species and mutate to more virulent strains had long been recognized by veterinarians, because coronaviruses inflict significant diseases in poultry, cattle and especially pigs. Virologists now are taking a closer look at some of these diseases to better understand coronaviruses and why they may be relevant for potential future pandemics.     

Deep-Sequencing Analysis of the Association between the Quasispecies Nature of the Hepatitis C Virus Core Region and Disease Progression

Variation of core amino acid (aa) 70 of hepatitis C virus (HCV) has been shown recently to be closely correlated with liver disease progression, suggesting that the core region might be present as a quasispecies during persistent infection and that this quasispecies nature might have an influence on the progression of disease. In our investigation, the subjects were 79 patients infected with HCV genotype 1b (25 with chronic hepatitis [CH], 29 with liver cirrhosis [LC], and 25 with hepatocellular carcinoma [HCC]). Deep sequencing of the HCV core region was carried out on their sera by using a Roche 454 GS Junior pyrosequencer. Based on a plasmid containing a cloned HCV sequence (pCV-J4L6S), the background error rate associated with pyrosequencing, including the PCR procedure, was calculated as 0.092 ± 0.005/base. Deep sequencing of the core region in the clinical samples showed a mixture of “mutant-type” Q/H and “wild-type” R at the core aa 70 position in most cases (71/79 [89.9%]), and the ratio of mutant residues to R in the mixture increased as liver disease advanced to LC and HCC. Meanwhile, phylogenetic analysis of the almost-complete core region revealed that the HCV isolates differed genetically depending on the mutation status at core aa 70. We conclude that the core aa 70 mixture ratio, determined by deep sequencing, reflected the status of liver disease, demonstrating a significant association between core aa 70 and disease progression in CH patients infected with HCV genotype 1b.     

A human coronavirus evolves antigenically to escape antibody immunity

There is intense interest in antibody immunity to coronaviruses. However, it is unknown if coronaviruses evolve to escape such immunity, and if so, how rapidly. Here we address this question by characterizing the historical evolution of human coronavirus 229E. We identify human sera from the 1980s and 1990s that have neutralizing titers against contemporaneous 229E that are comparable to the anti-SARS-CoV-2 titers induced by SARS-CoV-2 infection or vaccination. We test these sera against 229E strains isolated after sera collection, and find that neutralizing titers are lower against these “future” viruses. In some cases, sera that neutralize contemporaneous 229E viral strains with titers >1:100 do not detectably neutralize strains isolated 8–17 years later. The decreased neutralization of “future” viruses is due to antigenic evolution of the viral spike, especially in the receptor-binding domain. If these results extrapolate to other coronaviruses, then it may be advisable to periodically update SARS-CoV-2 vaccines.     

Real-Time PCR Investigation of Potential Vectors, Reservoirs, and Shedding Patterns of Feline Hemotropic Mycoplasmas

Three hemotropic mycoplasmas have been identified in pet cats: Mycoplasma haemofelis, “Candidatus Mycoplasma haemominutum,” and “Candidatus Mycoplasma turicensis.” The way in which these agents are transmitted is largely unknown. Thus, this study aimed to investigate fleas, ticks, and rodents as well as saliva and feces from infected cats for the presence of hemotropic mycoplasmas, to gain insight into potential transmission routes for these agents. DNA was extracted from arthropods and from rodent blood or tissue samples from Switzerland and from salivary and fecal swabs from two experimentally infected and six naturally infected cats. All samples were analyzed with real-time PCR, and some positive samples were confirmed by sequencing. Feline hemotropic mycoplasmas were detected in cat fleas and in a few Ixodes sp. and Rhipicephalus sp. ticks collected from animals but not in ticks collected from vegetation or from rodent samples, although the latter were frequently Mycoplasma coccoides PCR positive. When shedding patterns of feline hemotropic mycoplasmas were investigated, “Ca. Mycoplasma turicensis” DNA was detected in saliva and feces at the early but not at the late phase of infection. M. haemofelis and “Ca. Mycoplasma haemominutum” DNA was not amplified from saliva and feces of naturally infected cats, despite high hemotropic mycoplasma blood loads. Our results suggest that besides an ostensibly indirect transmission by fleas, direct transmission through saliva and feces at the early phase of infection could play a role in the epizootiology of feline hemotropic mycoplasmas. Neither the investigated tick nor the rodent population seems to represent a major reservoir for feline hemotropic mycoplasmas in Switzerland.     

Total infectome characterization of respiratory infections in pre-COVID-19 Wuhan,China

At the end of 2019 Wuhan witnessed an outbreak of “atypical pneumonia” that later developed into a global pandemic. Metagenomic sequencing rapidly revealed the causative agent of this outbreak to be a novel coronavirus denoted SARS-CoV-2. To provide a snapshot of the pathogens in pneumonia-associated respiratory samples from Wuhan prior to the emergence of SARS-CoV-2, we collected bronchoalveolar lavage fluid samples from 408 patients presenting with pneumonia and acute respiratory infections at the Central Hospital of Wuhan between 2016 and 2017. Unbiased total RNA sequencing was performed to reveal their “total infectome”, including viruses, bacteria and fungi. We identified 35 pathogen species, comprising 13 RNA viruses, 3 DNA viruses, 16 bacteria and 3 fungi, often at high abundance and including multiple co-infections (13.5%). SARS-CoV-2 was not present. These data depict a stable core infectome comprising common respiratory pathogens such as rhinoviruses and influenza viruses, an atypical respiratory virus (EV-D68), and a single case of a sporadic zoonotic pathogen–Chlamydia psittaci. Samples from patients experiencing respiratory disease on average had higher pathogen abundance than healthy controls. Phylogenetic analyses of individual pathogens revealed multiple origins and global transmission histories, highlighting the connectedness of the Wuhan population. This study provides a comprehensive overview of the pathogens associated with acute respiratory infections and pneumonia, which were more diverse and complex than obtained using targeted PCR or qPCR approaches. These data also suggest that SARS-CoV-2 or closely related viruses were absent from Wuhan in 2016–2017.     

Amplification Biases and Consistent Recovery of Loci in a Double-Digest RAD-seq Protocol

A growing variety of “genotype-by-sequencing” (GBS) methods use restriction enzymes and high throughput DNA sequencing to generate data for a subset of genomic loci, allowing the simultaneous discovery and genotyping of thousands of polymorphisms in a set of multiplexed samples. We evaluated a “double-digest” restriction-site associated DNA sequencing (ddRAD-seq) protocol by 1) comparing results for a zebra finch (Taeniopygia guttata) sample with in silico predictions from the zebra finch reference genome; 2) assessing data quality for a population sample of indigobirds (Vidua spp.); and 3) testing for consistent recovery of loci across multiple samples and sequencing runs. Comparison with in silico predictions revealed that 1) over 90% of predicted, single-copy loci in our targeted size range (178–328 bp) were recovered; 2) short restriction fragments (38–178 bp) were carried through the size selection step and sequenced at appreciable depth, generating unexpected but nonetheless useful data; 3) amplification bias favored shorter, GC-rich fragments, contributing to among locus variation in sequencing depth that was strongly correlated across samples; 4) our use of restriction enzymes with a GC-rich recognition sequence resulted in an up to four-fold overrepresentation of GC-rich portions of the genome; and 5) star activity (i.e., non-specific cutting) resulted in thousands of “extra” loci sequenced at low depth. Results for three species of indigobirds show that a common set of thousands of loci can be consistently recovered across both individual samples and sequencing runs. In a run with 46 samples, we genotyped 5,996 loci in all individuals and 9,833 loci in 42 or more individuals, resulting in <1% missing data for the larger data set. We compare our approach to similar methods and discuss the range of factors (fragment library preparation, natural genetic variation, bioinformatics) influencing the recovery of a consistent set of loci among samples.     

Evolutionary insights into the ecology of coronaviruses

Although many novel members of the Coronaviridae have recently been recognized in different species, the ecology of coronaviruses has not been established. Our study indicates that bats harbor a much wider diversity of coronaviruses than any other animal species. Dating of different coronavirus lineages suggests that bat coronaviruses are older than those recognized in other animals and that the human severe acute respiratory syndrome (SARS) coronavirus was directly derived from viruses from wild animals in wet markets of southern China. Furthermore, the most closely related bat and SARS coronaviruses diverged in 1986, an estimated divergence time of 17 years prior to the outbreak, suggesting that there may have been transmission via an unknown intermediate host. Analysis of lineage-specific selection pressure also indicated that only SARS coronaviruses in civets and humans were under significant positive selection, also demonstrating a recent interspecies transmission. Analysis of population dynamics revealed that coronavirus populations in bats have constant population growth, while viruses from all other hosts show epidemic-like increases in population. These results indicate that diverse coronaviruses are endemic in different bat species, with repeated introductions to other animals and occasional establishment in other species. Our findings suggest that bats are likely the natural hosts for all presently known coronavirus lineages and that all coronaviruses recognized in other species were derived from viruses residing in bats. Further surveillance of bat and other animal populations is needed to fully describe the ecology and evolution of this virus family.     

COVID‐19 and the gain of function debates: Improving biosafety measures requires a more precise definition of which experiments would raise safety concerns

The COVID‐19 pandemic has rekindled debates about gain‐of‐function experiments. This is an opportunity to clearly define safety risks and appropriate countermeasures. Subject Categories: Economics, Law & Politics, Microbiology, Virology & Host Pathogen Interaction, Science Policy & Publishing

The so‐called “gain of function” research has been recently debated in the context of viral research on coronaviruses and whether it is too risky to undertake such experiments. However, the meaning of “gain of function” or “GOF” in a science policy context has changed over time. The term was originally coined to describe two controversial research projects on H5N1 avian influenza virus and was later applied to specific experiments on coronavirus. Subsequent policies and discussions have attempted to define GOF in different ways, but no single definition has been widely accepted by the community. The fuzzy and imprecise nature of the term has led to misunderstandings and has hampered discussions on how to properly assess the benefit of such experiments and biosafety measures.

The fuzzy and imprecise nature of the term GOF has led to misunderstandings and has hampered discussions on how to properly assess the benefit of such experiments and biosafety measures

     

Evaluating the effectiveness of ensemble voting in improving the accuracy of consensus signals produced by various DTWA algorithms from step-current signals generated during nanopore sequencing

Nanopore sequencing device analysis systems simultaneously generate multiple picoamperage current signals representing the passage of DNA or RNA nucleotides ratcheted through a biomolecule nanopore array by motor proteins. Squiggles are a noisy and time-distorted representation of an underlying nucleotide sequence, “gold standard model”, due to experimental and algorithmic artefacts. Other research fields use dynamic time warped-space averaging (DTWA) algorithms to produce a consensus signal from multiple time-warped sources while preserving key features distorted by standard, linear-averaging approaches. We compared the ability of DTW Barycentre averaging (DBA), minimize mean (MM) and stochastic sub-gradient descent (SSG) DTWA algorithms to generate a consensus signal from squiggle-space ensembles of RNA molecules Enolase, Sequin R1-71-1 and Sequin R2-55-3 without knowledge of their associated gold standard model. We propose techniques to identify the leader and distorted squiggle features prior to DTWA consensus generation. New visualization and warping-path metrics are introduced to compare consensus signals and the best estimate of the “true” consensus, the study’s gold standard model. The DBA consensus was the best match to the gold standard for both Sequin studies but was outperformed in the Enolase study. Given an underlying common characteristic across a squiggle ensemble, we objectively evaluate a novel “voting scheme” that improves the local similarity between the consensus signal and a given fraction of the squiggle ensemble. While the gold standard is not used during voting, the increase in the match of the final voted-on consensus to the underlying Enolase and Sequin gold standard sequences provides an indirect success measure for the proposed voting procedure in two ways: First is the decreased least squares warped distance between the final consensus and the gold model, and second, the voting generates a final consensus length closer to the known underlying RNA biomolecule length. The results suggest considerable potential in marrying squiggle analysis and voted-on DTWA consensus signals to provide low-noise, low-distortion signals. This will lead to improved accuracy in detecting nucleotides and their deviation model due to chemical modifications (a.k.a. epigenetic information). The proposed combination of ensemble voting and DTWA has application in other research fields involving time-distorted, high entropy signals.     

Crystal Structure of the C-Terminal Cytoplasmic Domain of Non-Structural Protein 4 from Mouse Hepatitis Virus A59

Background

The replication of coronaviruses takes place on cytoplasmic double membrane vesicles (DMVs) originating in the endoplasmic reticulum (ER). Three trans-membrane non-structural proteins, nsp3, nsp4 and nsp6, are understood to be membrane anchors of the coronavirus replication complex. Nsp4 is localized to the ER membrane when expressed alone but is recruited into the replication complex in infected cells. It is revealed to contain four trans-membrane regions and its N- and C-termini are exposed to the cytosol.

Methodology/Principal Findings

We have determined the crystal structures of the C-terminal hydrophilic domain of nsp4 (nsp4C) from MHV strain A59 and a C425S site-directed mutant. The highly conserved 89 amino acid region from T408 to Q496 is shown to possess a new fold. The wild-type (WT) structure features two monomers linked by a Cys425-Cys425 disulfide bond in one asymmetric unit. The monomers are arranged with their N- and C-termini in opposite orientations to form an “open” conformation. Mutation of Cys425 to Ser did not affect the monomer structure, although the mutant dimer adopts strikingly different conformations by crystal packing, with the cross-linked C-termini and parallel N-termini of two monomers forming a “closed” conformation. The WT nsp4C exists as a dimer in solution and can dissociate easily into monomers in a reducing environment.

Conclusions/Significance

As nsp4C is exposed in the reducing cytosol, the monomer of nsp4C should be physiological. This structure may serve as a basis for further functional studies of nsp4.     

Genotype Imputation Reference Panel Selection Using Maximal Phylogenetic Diversity

The recent dramatic cost reduction of next-generation sequencing technology enables investigators to assess most variants in the human genome to identify risk variants for complex diseases. However, sequencing large samples remains very expensive. For a study sample with existing genotype data, such as array data from genome-wide association studies, a cost-effective approach is to sequence a subset of the study sample and then to impute the rest of the study sample, using the sequenced subset as a reference panel. The use of such an internal reference panel identifies population-specific variants and avoids the problem of a substantial mismatch in ancestry background between the study population and the reference population. To efficiently select an internal panel, we introduce an idea of phylogenetic diversity from mathematical phylogenetics and comparative genomics. We propose the “most diverse reference panel”, defined as the subset with the maximal “phylogenetic diversity”, thereby incorporating individuals that span a diverse range of genotypes within the sample. Using data both from simulations and from the 1000 Genomes Project, we show that the most diverse reference panel can substantially improve the imputation accuracy compared to randomly selected reference panels, especially for the imputation of rare variants. The improvement in imputation accuracy holds across different marker densities, reference panel sizes, and lengths for the imputed segments. We thus propose a novel strategy for planning sequencing studies on samples with existing genotype data.     

Middle East Respiratory Syndrome Coronavirus Intra-Host Populations Are Characterized by Numerous High Frequency Variants

Middle East respiratory syndrome coronavirus (MERS-CoV) is an emerging human pathogen related to SARS virus. In vitro studies indicate this virus may have a broad host range suggesting an increased pandemic potential. Genetic and epidemiological evidence indicate camels serve as a reservoir for MERS virus but the mechanism of cross species transmission is unclear and many questions remain regarding the susceptibility of humans to infection. Deep sequencing data was obtained from the nasal samples of three camels that had been experimentally infected with a human MERS-CoV isolate. A majority of the genome was covered and average coverage was greater than 12,000x depth. Although only 5 mutations were detected in the consensus sequences, 473 intrahost single nucleotide variants were identified. Many of these variants were present at high frequencies and could potentially influence viral phenotype and the sensitivity of detection assays that target these regions for primer or probe binding.     

Intrinsic resistance of feline peritoneal macrophages to coronavirus infection correlates with in vivo virulence.

Cats infected with virulent feline coronavirus strains develop feline infectious peritonitis, an invariably fatal, immunologically mediated disease; avirulent strains cause either clinically inapparent infection or mild enteritis. Four virulent coronavirus isolates and five avirulent isolates were assessed by immunofluorescence and virus titration for their ability to infect and replicate in feline peritoneal macrophages in vitro. The avirulent coronaviruses infected fewer macrophages, produced lower virus titers, were less able to sustain viral replication, and spread less efficiently to other susceptible macrophages than the virulent coronaviruses. Thus, the intrinsic resistance of feline macrophages may play a pivotal role in the outcome of coronavirus infection in vivo.     

In Vitro HIV-1 Evolution in Response to Triple Reverse Transcriptase Inhibitors & In Silico Phenotypic Analysis

Background

Effectiveness of ART regimens strongly depends upon complex interactions between the selective pressure of drugs and the evolution of mutations that allow or restrict drug resistance.

Methods

Four clinical isolates from NRTI-exposed, NNRTI-naive subjects were passaged in increasing concentrations of NVP in combination with 1 µM 3 TC and 2 µM ADV to assess selective pressures of multi-drug treatment. A novel parameter inference procedure, based on a stochastic viral growth model, was used to estimate phenotypic resistance and fitness from in vitro combination passage experiments.

Results

Newly developed mathematical methods estimated key phenotypic parameters of mutations arising through selective pressure exerted by 3 TC and NVP. Concentrations of 1 µM 3 TC maintained the M184V mutation, which was associated with intrinsic fitness deficits. Increasing NVP concentrations selected major NNRTI resistance mutations. The evolutionary pathway of NVP resistance was highly dependent on the viral genetic background, epistasis as well as stochasticity. Parameter estimation indicated that the previously unrecognized mutation L228Q was associated with NVP resistance in some isolates.

Conclusion

Serial passage of viruses in the presence of multiple drugs may resemble the selection of mutations observed among treated individuals and populations in vivo and indicate evolutionary preferences and restrictions. Phenotypic resistance estimated here “in silico” from in vitro passage experiments agreed well with previous knowledge, suggesting that the unique combination of “wet-” and “dry-lab” experimentation may improve our understanding of HIV-1 resistance evolution in the future.     

Genome-wide mapping of SARS-CoV-2 RNA structures identifies therapeutically-relevant elements

SARS-CoV-2 is a betacoronavirus with a linear single-stranded, positive-sense RNA genome, whose outbreak caused the ongoing COVID-19 pandemic. The ability of coronaviruses to rapidly evolve, adapt, and cross species barriers makes the development of effective and durable therapeutic strategies a challenging and urgent need. As for other RNA viruses, genomic RNA structures are expected to play crucial roles in several steps of the coronavirus replication cycle. Despite this, only a handful of functionally-conserved coronavirus structural RNA elements have been identified to date. Here, we performed RNA structure probing to obtain single-base resolution secondary structure maps of the full SARS-CoV-2 coronavirus genome both in vitro and in living infected cells. Probing data recapitulate the previously described coronavirus RNA elements (5′ UTR and s2m), and reveal new structures. Of these, ∼10.2% show significant covariation among SARS-CoV-2 and other coronaviruses, hinting at their functionally-conserved role. Secondary structure-restrained 3D modeling of these segments further allowed for the identification of putative druggable pockets. In addition, we identify a set of single-stranded segments in vivo, showing high sequence conservation, suitable for the development of antisense oligonucleotide therapeutics. Collectively, our work lays the foundation for the development of innovative RNA-targeted therapeutic strategies to fight SARS-related infections.