Optimized quantification of intra-host viral diversity in SARS-CoV-2 and influenza virus sequence data.

High error rates of viral RNA-dependent RNA polymerases lead to diverse intra-host viral populations during infection. Errors made during replication that are not strongly deleterious to the virus can lead to the generation of minority variants. However, accurate detection of minority variants in viral sequence data is complicated by errors introduced during sample preparation and data analysis. We used synthetic RNA controls and simulated data to test seven variant-calling tools across a range of allele frequencies and simulated coverages. We show that choice of variant caller and use of replicate sequencing have the most significant impact on single-nucleotide variant (SNV) discovery and demonstrate how both allele frequency and coverage thresholds impact both false discovery and false-negative rates. When replicates are not available, using a combination of multiple callers with more stringent cutoffs is recommended. We use these parameters to find minority variants in sequencing data from SARS-CoV-2 clinical specimens and provide guidance for studies of intra-host viral diversity using either single replicate data or data from technical replicates. Our study provides a framework for rigorous assessment of technical factors that impact SNV identification in viral samples and establishes heuristics that will inform and improve future studies of intra-host variation, viral diversity, and viral evolution. IMPORTANCE When viruses replicate inside a host cell, the virus replication machinery makes mistakes. Over time, these mistakes create mutations that result in a diverse population of viruses inside the host. Mutations that are neither lethal to the virus nor strongly beneficial can lead to minority variants that are minor members of the virus population. However, preparing samples for sequencing can also introduce errors that resemble minority variants, resulting in the inclusion of false-positive data if not filtered correctly. In this study, we aimed to determine the best methods for identification and quantification of these minority variants by testing the performance of seven commonly used variant-calling tools. We used simulated and synthetic data to test their performance against a true set of variants and then used these studies to inform variant identification in data from SARS-CoV-2 clinical specimens. Together, analyses of our data provide extensive guidance for future studies of viral diversity and evolution.

Optimized Quantification of Intrahost Viral Diversity in SARS-CoV-2 and Influenza Virus Sequence Data.

High error rates of viral RNA-dependent RNA polymerases lead to diverse intra-host viral populations during infection. Errors made during replication that are not strongly deleterious to the virus can lead to the generation of minority variants. However, accurate detection of minority variants in viral sequence data is complicated by errors introduced during sample preparation and data analysis. We used synthetic RNA controls and simulated data to test seven variant calling tools across a range of allele frequencies and simulated coverages. We show that choice of variant caller, and use of replicate sequencing have the most significant impact on single nucleotide variant (SNV) discovery and demonstrate how both allele frequency and coverage thresholds impact both false discovery and false negative rates. We use these parameters to find minority variants in sequencing data from SARS-CoV-2 clinical specimens and provide guidance for studies of intrahost viral diversity using either single replicate data or data from technical replicates. Our study provides a framework for rigorous assessment of technical factors that impact SNV identification in viral samples and establishes heuristics that will inform and improve future studies of intrahost variation, viral diversity, and viral evolution. IMPORTANCE: When viruses replicate inside a host, the virus replication machinery makes mistakes. Over time, these mistakes create mutations that result in a diverse population of viruses inside the host. Mutations that are neither lethal to the virus, nor strongly beneficial, can lead to minority variants that are minor members of the virus population. However, preparing samples for sequencing can also introduce errors that resemble minority variants, resulting in inclusion of false positive data if not filtered correctly. In this study, we aimed to determine the best methods for identification and quantification of these minority variants by testing the performance of seven commonly used variant calling tools. We used simulated and synthetic data to test their performance against a true set of variants, and then used these studies to inform variant identification in data from clinical SARS-CoV-2 clinical specimens. Together, analyses of our data provide extensive guidance for future studies of viral diversity and evolution.

Increased risk of influenza among vaccinated adults who are obese.

BACKGROUND: Influenza infects 5-15% of the global population each year, and obesity has been shown to be an independent risk factor for increased influenza-related complications including hospitalization and death. However, the risk of developing influenza or influenza-like illness (ILI) in a vaccinated obese adult population has not been addressed. OBJECTIVE: This study evaluated whether obesity was associated with increased risk of influenza and ILI among vaccinated adults. SUBJECTS AND METHODS: During the 2013-2014 and 2014-2015 influenza seasons, we recruited 1042 subjects to a prospective observational study of trivalent inactivated influenza vaccine (IIV3) in adults. A total of 1022 subjects completed the study. Assessments of relative risk for laboratory confirmed influenza and ILI were determined based on body mass index. Seroconversion and seroprotection rates were determined using prevaccination and 26-35 days post vaccination serum samples. Recruitment criteria for this study were adults 18 years of age and older receiving the seasonal trivalent inactivated influenza vaccine (IIV3) for the years 2013-2014 and 2014-2015. Exclusion criteria were immunosuppressive diseases, use of immunomodulatory or immunosuppressive drugs, acute febrile illness, history of Guillain-Barre syndrome, use of theophylline preparations or use of warfarin. RESULTS: Among obese, 9.8% had either confirmed influenza or influenza-like-illness compared with 5.1% of healthy weight participants. Compared with vaccinated healthy weight, obese participants had double the risk of developing influenza or ILI (relative risk=2.01, 95% CI 1.12, 3.60, P=0.020). Seroconversion or seroprotection rates were not different between healthy weight and obese adults with influenza or ILI. CONCLUSIONS: Despite robust serological responses, vaccinated obese adults are twice as likely to develop influenza and ILI compared with healthy weight adults. This finding challenges the current standard for correlates of protection, suggesting use of antibody titers to determine vaccine effectiveness in an obese population may provide misleading information.

Presence of influenza viruses in backyard poultry and swine in El Yali wetland, Chile.

In South America little is known regarding influenza virus circulating in backyard poultry and swine populations. Backyard productive systems (BPS) that breed swine and poultry are widely distributed throughout Chile with high density in the central zone, and several BPS are located within the "El Yali" (EY) ecosystem, which is one of the most important wetlands in South America. Here, 130 different wild bird species have been described, of them, at least 22 species migrate yearly from North America for nesting. For this reason, EY is considered as a high-risk zone for avian influenza virus. This study aims to identify if backyard poultry and swine bred in the EY ecosystem have been exposed to influenza A virus and if so, to identify influenza virus subtypes. A biosecurity and handling survey was applied and samples were collected from BPS in two seasons (spring 2013 and fall 2014) for influenza seroprevalence, and in one season (fall 2014) for virus presence. Seroprevalence at BPS level was 42% (95% CI:22-49) during spring 2013 and 60% (95% CI 43-72) in fall 2014. rRT-PCR for the influenza A matrix gene indicated a viral prevalence of 27% (95% CI:14-39) at BPS level in fall 2014. Eight farms (73% of rRT-PCR positive farms) were also positive to the Elisa test at the same time. One BPS was simultaneously positive (rRT-PCR) in multiple species (poultry, swine and geese) and a H1N2 virus was identified from swine, exemplifying the risk that these BPS may pose for generation of novel influenza viruses.

Influenza gain-of-function experiments: their role in vaccine virus recommendation and pandemic preparedness.

In recent years, controversy has arisen regarding the risks and benefits of certain types of gain-of-function (GOF) studies involving avian influenza viruses. In this article, we provide specific examples of how different types of data, including information garnered from GOF studies, have helped to shape the influenza vaccine production process-from selection of candidate vaccine viruses (CVVs) to the manufacture and stockpiling of safe, high-yield prepandemic vaccines for the global community. The article is not written to support a specific pro- or anti-GOF stance but rather to inform the scientific community about factors involved in vaccine virus selection and the preparation of prepandemic influenza vaccines and the impact that some GOF information has had on this process.

Inhibition of Vaccinia virus entry by a broad spectrum antiviral peptide.

Concerns about the possible use of Variola virus, the causative agent of smallpox, as a weapon for bioterrorism have led to renewed efforts to identify new antivirals against orthopoxviruses. We identified a peptide, EB, which inhibited infection by Vaccinia virus with an EC(50) of 15 microM. A control peptide, EBX, identical in composition to EB but differing in sequence, was inactive (EC50>200 microM), indicating sequence specificity. The inhibition was reversed upon removal of the peptide, and EB treatment had no effect on the physical integrity of virus particles as determined by electron microscopy. Viral adsorption was unaffected by the presence of EB, and the addition of EB post-entry had no effect on viral titers or on early gene expression. The addition of EB post-adsorption resulted in the inhibition of beta-galactosidase expression from an early viral promoter with an EC(50) of 45 microM. A significant reduction in virus entry was detected in the presence of the peptide when the number of viral cores released into the cytoplasm was quantified. Electron microscopy indicated that 88% of the virions remained on the surface of cells in the presence of EB, compared to 37% in the control (p<0.001). EB also blocked fusion-from-within, suggesting that virus infection is inhibited at the fusion step. Analysis of EB derivatives suggested that peptide length may be important for the activity of EB. The EB peptide is, to our knowledge, the first known small molecule inhibitor of Vaccinia virus entry.

Examining the cellular pathways involved in influenza virus induced apoptosis.

Apoptosis is essential in many physiological processes including wound healing and development of the immune response. Apoptosis also plays an important role in the pathogenesis of many infectious diseases including those caused by viruses. Influenza viruses induce apoptosis in cells that are permissive for viral replication and cells that do not support viral replication. The cellular pathways involved in influenza virus induced apoptosis are currently ill defined. Previous studies suggest that influenza virus infection increased the expression of the Fas antigen in HeLa cells, and that Fas antigen is partially involved in apoptosis. In these studies we examined the cellular pathways involved in avian influenza virus induced apoptosis in two cell lines that support productive viral replication: Madin-Darby canine kidney cells (MDCK) and mink lung epithelial (Mv1Lu) cells.

Recombinant paramyxovirus type 1-avian influenza-H7 virus as a vaccine for protection of chickens against influenza and Newcastle disease.

Current vaccines to prevent avian influenza rely upon labor-intensive parenteral injection. A more advantageous vaccine would be capable of administration by mass immunization methods such as spray or water vaccination. A recombinant vaccine (rNDV-AIV-H7) was constructed by using a lentogenic paramyxovirus type 1 vector (Newcastle disease virus [NDV] B1 strain) with insertion of the hemagglutinin (HA) gene from avian influenza virus (AIV) A/chicken/NY/13142-5/94 (H7N2). The recombinant virus had stable insertion and expression of the H7 AIV HA gene as evident by detection of HA expression via immunofluorescence in infected Vero cells. The rNDV-AIV-H7 replicated in 9-10 day embryonating chicken eggs and exhibited hemagglutinating activity from both NDV and AI proteins that was inhibited by antisera against both NDV and AIV H7. Groups of 2-week-old white Leghorn chickens were vaccinated with transfectant NDV vector (tNDV), rNDV-AIV-H7, or sterile allantoic fluid and were challenged 2 weeks later with viscerotropic velogenic NDV (vvNDV) or highly pathogenic (HP) AIV. The sham-vaccinated birds were not protected from vvNDV or HP AIV challenge. The transfectant NDV vaccine provided 70% protection for NDV challenge but did not protect against AIV challenge. The rNDV-AIV-H7 vaccine provided partial protection (40%) from vvNDV and HP AIV challenge. The serologic response was examined in chickens that received one or two immunizations of the rNDV-AIV-H7 vaccine. Based on hemagglutination inhibition and enzyme-linked immunosorbent assay (ELISA) tests, chickens that received a vaccine boost seroconverted to AIV H7, but the serologic response was weak in birds that received only one vaccination. This demonstrates the potential for NDV for use as a vaccine vector in expressing AIV proteins.

Localization of astrovirus in experimentally infected turkeys as determined by in situ hybridization.

Twenty-one 3-day-old turkey poults from British United Turkeys of America were orally inoculated with a recently characterized astrovirus, TAstV-2, isolated from turkeys with poult enteritis and mortality syndrome. At 1, 2, 3, 4, 5, 7, and 9 days postinfection (dpi), three inoculated birds were euthanatized, and tissues (intestines, spleen, bursa, and thymus) were collected immediately into 10% neutral buffered formalin. Inoculated birds were diarrheic by 3 dpi, and frothy feces persisted throughout the experimental period. Histologically, there was only slight evidence of enteric damage, which was characterized by mild epithelial necrosis, lamina propria infiltrates, minimal villus atrophy, and mild crypt hyperplasia. In situ hybridization, using a negative sense digoxigenin-labeled riboprobe to the capsid gene of TAstV-2, revealed viral RNA in intestinal epithelial cells at the basal margins of the villi, in distal small intestine, and in cecum at 2 dpi, with subsequent extension to epithelium of the large intestine and proximal small intestine (3-5 dpi). Minimal virus remained by 9 dpi.

Influenza virus ns1 protein induces apoptosis in cultured cells.

The importance of influenza viruses as worldwide pathogens in humans, domestic animals, and poultry is well recognized. Discerning how influenza viruses interact with the host at a cellular level is crucial for a better understanding of viral pathogenesis. Influenza viruses induce apoptosis through mechanisms involving the interplay of cellular and viral factors that may depend on the cell type. However, it is unclear which viral genes induce apoptosis. In these studies, we show that the expression of the nonstructural (NS) gene of influenza A virus is sufficient to induce apoptosis in MDCK and HeLa cells. Further studies showed that the multimerization domain of the NS1 protein but not the effector domain is required for apoptosis. However, this mutation is not sufficient to inhibit apoptosis using whole virus.

Inactivation of an astrovirus associated with poult enteritis mortality syndrome.

Outbreaks of poult enteritis mortality syndrome (PEMS) continue to cause financial losses to the turkey industry. Clinically, PEMS is defined by mortality profiles, diarrhea, flock unevenness, and immunosuppression. PEMS is a very difficult disease to control and prevent. Depopulation of PEMS-affected flocks and thorough cleaning of the contaminated housing have failed to prevent infection (disease) in subsequent flock placements. The relationship of PEMS to other enteric disease complexes of young turkeys is unknown, partly because the causative agent of PEMS remains unknown. Recently, we isolated a unique astrovirus strain from the thymus and intestines of PEMS-infected poults. This strain is molecularly and serologically distinct from the astrovirus that circulated in turkeys in the 1980s. Mammalian astroviruses are very resistant to inactivation. In these studies, we examined the stability of partially purified PEMS-associated astrovirus to inactivation with heat, laboratory disinfectants, and commercial disinfectants used in commercial turkey houses in an embryonated egg model system. Similar to mammalian astroviruses, the PEMS-associated astrovirus is resistant to inactivation by heat, acidification, detergent treatment, and treatment with phenolic, quaternary ammonium, or benzalkonium chloride-based products. Only treatment with formaldehyde, beta-propriolactone, or the peroxymonosulfate-based product Virkon S completely inactivated the astrovirus in the embryo model. These studies provide an alternate means to potentially control at least one virus associated with PEMS through the use of specific disinfectants.

Tamoxifen and estrogen effects on TGF-beta formation: role of thrombospondin-1, alphavbeta3, and integrin-associated protein.

We have found that the enhanced activation of latent TGF-beta by human breast carcinoma cell lines either treated with tamoxifen or deprived of estrogen is dependent upon thrombospondin-1 (TSP-1) since activation was blocked by anti-TSP-1 antibodies or by a TSP antagonist peptide. However, TGF-beta formation upon tamoxifen exposure to estrogen withdrawal is associated with decreased levels of soluble TSP-1. A concomitant increase in the expression of the TSP-1 receptors alphavbeta3 and integrin-associated protein (IAP) occurs under these conditions, and antibodies to TSP-1 or to these receptors inhibit increased TGF-beta formation. Therefore, increased cell surface associated TSP-1 enhances latent TGF-beta activation.

Influenza virus (A/HK/156/97) hemagglutinin expressed by an alphavirus replicon system protects chickens against lethal infection with Hong Kong-origin H5N1 viruses.

Venezuelan equine encephalitis virus replicon particles (VRP) containing the gene expressing hemagglutinin (HA) from the human Hong Kong Influenza A isolate (A/HK/156/97) were evaluated as vaccines in chicken embryos and young chicks. Expressed HA was readily detected in bird-tissue staining with anti-H5 HA antibody and in chicken cells infected with the replicon preparations following immunoprecipitation with monoclonal antibody. Birds challenged with a dose of the lethal parent virus were protected to different extents depending on the age of the bird. In ovo and 1-day-old inoculated animals that received no boost with the VRP were partially protected; birds 2 weeks of age were completely protected with a single dose of VRP.

Newcastle disease virus phosphoprotein gene analysis and transcriptional editing in avian cells.

Nucleotide sequence was determined for the phosphoprotein (P) gene from 23 Newcastle disease virus (NDV) isolates representing all defined pathotypes with different chronological and geographic origins. Sequence variation, with synonymous substitutions dominating, occurred throughout the P gene. An exception was a conserved central region containing the transcriptional editing site. Four G nucleotide additions were detected in NDV P gene mRNA potentially creating alternative open reading frames. However, only one in-frame stop codon exists with a single G addition among all isolates that would allow for a potential V protein. A second potential stop codon does not exist in the P gene consensus sequence among all isolates with more than one G nucleotide addition at the editing site. This precludes a possible W protein in these isolates. A second potential alternative in-frame start site exists among all isolates that could encode a predicted X protein for NDV. Comparison of the P gene editing sites among the Paramyxovirinae and predicted P gene usage demonstrates that NDV more closely resembles the respiroviruses and morbilliviruses. Phylogenetic analysis of P gene sequences among NDV isolates demonstrates there are two clades of these viruses. One group includes viruses isolated in the US prior to 1970, while a second cluster includes virulent viruses circulating worldwide.

Development of an RT-PCR diagnostic test for an avian astrovirus.

Astroviruses are small round viruses that cause enteric disease in the young of several species. Detection and diagnosis of astrovirus infection in non-human hosts relies heavily on electron microscopy and fluorescent antibody tests. Recently, our laboratory isolated and sequenced an avian astrovirus from poult enteritis mortality syndrome affected turkeys. These studies describe the development of RT-PCR methods, which specifically detect regions of the viral capsid and polymerase genes, and demonstrate their use in detecting astrovirus infection in commercial turkey flocks.

Identifying agent(s) associated with poult enteritis mortality syndrome: importance of the thymus.

Poult enteritis mortality syndrome (PEMS), a highly infectious disease of young turkeys, causes serious financial losses to the turkey industry. Clinically, PEMS is defined by mortality profiles, diarrhea, growth depression, and immunosuppression. Although many viruses, bacteria, and parasites are found in PEMS-infected birds, the inciting agent remains unknown. Experimentally, PEMS can be reproduced by exposing naïve poults to the intestinal contents from infected birds. Previous reports suggest that extraintestinal tissues fail to reproduce the disease. Histopathologic examination of tissues from PEMS-infected poults suggested that the thymus exhibited the earliest signs of pathology. On the basis of these observations, we hypothesized that the thymus harbors an agent(s) involved in PEMS. In these studies, naïve turkey poults were orally inoculated with a bacteria-free filtrate composed of either the intestines and feces or the thymus from PEMS-infected birds and were monitored for clinical signs of PEMS. Poults exposed to a filtrate composed solely of the thymus from PEMS-infected birds exhibited diarrhea, growth depression, mortality, pathology, and, most importantly, immunosuppression similar to poults exposed to the intestinal filtrate. The results of this study suggest that the thymus of infected birds harbors the agent(s) that can reproduce a PEMS-like disease in turkey poults.

Continued circulation in China of highly pathogenic avian influenza viruses encoding the hemagglutinin gene associated with the 1997 H5N1 outbreak in poultry and humans.

Since the outbreak in humans of an H5N1 avian influenza virus in Hong Kong in 1997, poultry entering the live-bird markets of Hong Kong have been closely monitored for infection with avian influenza. In March 1999, this monitoring system detected geese that were serologically positive for H5N1 avian influenza virus, but the birds were marketed before they could be sampled for virus. However, viral isolates were obtained by swabbing the cages that housed the geese. These samples, known collectively as A/Environment/Hong Kong/437/99 (A/Env/HK/437/99), contained four viral isolates, which were compared to the 1997 H5N1 Hong Kong isolates. Analysis of A/Env/HK/437/99 viruses revealed that the four isolates are nearly identical genetically and are most closely related to A/Goose/Guangdong/1/96. These isolates and the 1997 H5N1 Hong Kong viruses encode common hemagglutinin (H5) genes that have identical hemagglutinin cleavage sites. Thus, the pathogenicity of the A/Env/HK/437/99 viruses was compared in chickens and in mice to evaluate the potential for disease outbreaks in poultry and humans. The A/Env/HK/437/99 isolates were highly pathogenic in chickens but caused a longer mean death time and had altered cell tropism compared to A/Hong Kong/156/97 (A/HK/156/97). Like A/HK/156/97, the A/Env/HK/437/99 viruses replicated in mice and remained localized to the respiratory tract. However, the A/Env/HK/437/99 isolates caused only mild pathological lesions in these tissues and no clinical signs of disease or death. As a measure of the immune response to these viruses, transforming growth factor beta levels were determined in the serum of infected mice and showed elevated levels for the A/Env/HK/437/99 viruses compared to the A/HK/156/97 viruses. This study is the first to characterize the A/Env/HK/437/99 viruses in both avian and mammalian species, evaluating the H5 gene from the 1997 Hong Kong H5N1 isolates in a different genetic background. Our findings reveal that at least one of the avian influenza virus genes encoded by the 1997 H5N1 Hong Kong viruses continues to circulate in mainland China and that this gene is important for pathogenesis in chickens but is not the sole determinant of pathogenicity in mice. There is evidence that H9N2 viruses, which have internal genes in common with the 1997 H5N1 Hong Kong isolates, are still circulating in Hong Kong and China as well, providing a heterogeneous gene pool for viral reassortment. The implications of these findings for the potential for human disease are discussed.

Molecular characterization of an avian astrovirus.

Astroviruses are known to cause enteric disease in several animal species, including turkeys. However, only human astroviruses have been well characterized at the nucleotide level. Herein we report the nucleotide sequence, genomic organization, and predicted amino acid sequence of a turkey astrovirus isolated from poults with an emerging enteric disease.

Immunology of avian influenza virus: a review.

Avian influenza virus can cause serious disease in a wide variety of birds and mammals, but its natural host range is in wild ducks, gulls, and shorebirds. Infections in poultry can be inapparent or cause respiratory disease, decreases in production, or a rapidly fatal systemic disease known as highly pathogenic avian influenza (HPAI). For the protection of poultry, neutralizing antibody to the hemagglutinin and neuraminidase proteins provide the primary protection against disease. A variety of vaccines elicit neutralizing antibody, including killed whole virus vaccines and fowl-pox recombinant vaccines. Antigenic drift of influenza viruses appears to be less important in causing vaccine failures in poultry as compared to humans. The cytotoxic T lymphocyte response can reduce viral shedding in mildly pathogenic avian influenza viruses, but provides questionable protection against HPAI. Influenza viruses can directly affect the immune response of infected birds, and the role of the Mx gene, interferons, and other cytokines in protection from disease remains unknown.

Distinct pathogenesis of hong kong-origin H5N1 viruses in mice compared to that of other highly pathogenic H5 avian influenza viruses.

In 1997, an outbreak of virulent H5N1 avian influenza virus occurred in poultry in Hong Kong (HK) and was linked to a direct transmission to humans. The factors associated with transmission of avian influenza virus to mammals are not fully understood, and the potential risk of other highly virulent avian influenza A viruses infecting and causing disease in mammals is not known. In this study, two avian and one human HK-origin H5N1 virus along with four additional highly pathogenic H5 avian influenza viruses were analyzed for their pathogenicity in 6- to 8-week-old BALB/c mice. Both the avian and human HK H5 influenza virus isolates caused severe disease in mice, characterized by induced hypothermia, clinical signs, rapid weight loss, and 75 to 100% mortality by 6 to 8 days postinfection. Three of the non-HK-origin isolates caused no detectable clinical signs. One isolate, A/tk/England/91 (H5N1), induced measurable disease, and all but one of the animals recovered. Infections resulted in mild to severe lesions in both the upper and lower respiratory tracts. Most consistently, the viruses caused necrosis in respiratory epithelium of the nasal cavity, trachea, bronchi, and bronchioles with accompanying inflammation. The most severe and widespread lesions were observed in the lungs of HK avian influenza virus-infected mice, while no lesions or only mild lesions were evident with A/ck/Scotland/59 (H5N1) and A/ck/Queretaro/95 (H5N2). The A/ck/Italy/97 (H5N2) and the A/tk/England/91 (H5N1) viruses exhibited intermediate pathogenicity, producing mild to moderate respiratory tract lesions. In addition, infection by the different isolates could be further distinguished by the mouse immune response. The non-HK-origin isolates all induced production of increased levels of active transforming growth factor beta following infection, while the HK-origin isolates did not.