An optimized methodology for whole genome sequencing of RNA respiratory viruses from nasopharyngeal aspirates.

Over the last decade, the number of viral genome sequences deposited in available databases has grown exponentially. However, sequencing methodology vary widely and many published works have relied on viral enrichment by viral culture or nucleic acid amplification with specific primers rather than through unbiased techniques such as metagenomics. The genome of RNA viruses is highly variable and these enrichment methodologies may be difficult to achieve or may bias the results. In order to obtain genomic sequences of human respiratory syncytial virus (HRSV) from positive nasopharyngeal aspirates diverse methodologies were evaluated and compared. A total of 29 nearly complete and complete viral genomes were obtained. The best performance was achieved with a DNase I treatment to the RNA directly extracted from the nasopharyngeal aspirate (NPA), sequence-independent single-primer amplification (SISPA) and library preparation performed with Nextera XT DNA Library Prep Kit with manual normalization. An average of 633,789 and 1,674,845 filtered reads per library were obtained with MiSeq and NextSeq 500 platforms, respectively. The higher output of NextSeq 500 was accompanied by the increasing of duplicated reads percentage generated during SISPA (from an average of 1.5% duplicated viral reads in MiSeq to an average of 74% in NextSeq 500). HRSV genome recovery was not affected by the presence or absence of duplicated reads but the computational demand during the analysis was increased. Considering that only samples with viral load ≥ E+06 copies/ml NPA were tested, no correlation between sample viral loads and number of total filtered reads was observed, nor with the mapped viral reads. The HRSV genomes showed a mean coverage of 98.46% with the best methodology. In addition, genomes of human metapneumovirus (HMPV), human rhinovirus (HRV) and human parainfluenza virus types 1-3 (HPIV1-3) were also obtained with the selected optimal methodology.

Dual-strain genital herpes simplex virus type 2 (HSV-2) infection in the US, Peru, and 8 countries in sub-Saharan Africa: A nested cross-sectional viral genotyping study.

BACKGROUND: Quantitative estimation of the extent to which the immune system's protective effect against one herpes simplex virus type 2 (HSV-2) infection protects against infection with additional HSV-2 strains is important for understanding the potential for HSV-2 vaccine development. Using viral genotyping, we estimated the prevalence of HSV-2 dual-strain infection and identified risk factors. METHODS AND FINDINGS: People with and without HIV infection participating in HSV-2 natural history studies (University of Washington Virology Research Clinic) and HIV prevention trials (HIV Prevention Trials Network 039 and Partners in Prevention HSV/HIV Transmission Study) in the US, Africa, and Peru with 2 genital specimens each containing ≥105 copies herpes simplex virus DNA/ml collected a median of 5 months apart (IQR: 2-11 months) were included. It is unlikely that 2 strains would be detected in the same sample simultaneously; therefore, 2 samples were required to detect dual-strain infection. We identified 85 HSV-2 SNPs that, in aggregate, could determine whether paired HSV-2 strains were the same or different with >90% probability. These SNPs were then used to create a customized high-throughput array-based genotyping assay. Participants were considered to be infected with more than 1 strain of HSV-2 if their samples differed by ≥5 SNPs between the paired samples, and dual-strain infection was confirmed using high-throughput sequencing (HTS). We genotyped pairs of genital specimens from 459 people; 213 (46%) were men, the median age was 34 years (IQR: 27-44), and 130 (28%) were HIV seropositive. Overall, 272 (59%) people were from the US, 59 (13%) were from Peru, and 128 (28%) were from 8 countries in Africa. Of the 459 people, 18 (3.9%) met the criteria for dual-strain infection. HTS and phylogenetic analysis of paired specimens confirmed shedding of 2 distinct HSV-2 strains collected at different times in 17 pairs, giving an estimated dual-strain infection prevalence of 3.7% (95% CI = 2.0%-5.4%). Paired samples with dual-strain infection differed by a median of 274 SNPs in the UL_US region (range 129-413). Matching our observed dual-strain infection frequency to simulated data of varying prevalences and allowing only 2 samples per person, we inferred the true prevalence of dual-strain infection to be 7%. In multivariable analysis, controlling for HIV status and continent of origin, people from Africa had a higher risk for dual-strain infection (risk ratio [RR] = 9.20, 95% CI = 2.05-41.32), as did people who were HIV seropositive (RR = 4.06, 95% CI = 1.42-11.56). CONCLUSIONS: HSV-2 dual-strain infection was detected in 3.7% of paired samples from individual participants, and was more frequent among people with HIV infection. Simulations suggest that the true prevalence of dual-strain infection is 7%. Our data indicate that naturally occurring immunity to HSV-2 may be protective against infection with a second strain. This study is limited by the inability to determine the timing of acquisition of the second strain.

Highly conserved intragenic HSV-2 sequences: Results from next-generation sequencing of HSV-2 UL and US regions from genital swabs collected from 3 continents.

INTRODUCTION: Understanding the variability in circulating herpes simplex virus type 2 (HSV-2) genomic sequences is critical to the development of HSV-2 vaccines. METHODS: Genital lesion swabs containing ≥ 107log10 copies HSV DNA collected from Africa, the USA, and South America underwent next-generation sequencing, followed by K-mer based filtering and de novo genomic assembly. Sites of heterogeneity within coding regions in unique long and unique short (UL_US) regions were identified. Phylogenetic trees were created using maximum likelihood reconstruction. RESULTS: Among 46 samples from 38 persons, 1468 intragenic base-pair substitutions were identified. The maximum nucleotide distance between strains for concatenated UL_US segments was 0.4%. Phylogeny did not reveal geographic clustering. The most variable proteins had non-synonymous mutations in < 3% of amino acids. CONCLUSIONS: Unenriched HSV-2 DNA can undergo next-generation sequencing to identify intragenic variability. The use of clinical swabs for sequencing expands the information that can be gathered directly from these specimens.

Distinct Zika Virus Lineage in Salvador, Bahia, Brazil.

Sequencing of isolates from patients in Bahia, Brazil, where most Zika virus cases in Brazil have been reported, resulted in 11 whole and partial Zika virus genomes. Phylogenetic analyses revealed a well-supported Bahia-specific Zika virus lineage, which indicates sustained Zika virus circulation in Salvador, Bahia's capital city, since mid-2014.

DNA microarray for detection of gastrointestinal viruses.

Gastroenteritis is a clinical illness of humans and other animals that is characterized by vomiting and diarrhea and caused by a variety of pathogens, including viruses. An increasing number of viral species have been associated with gastroenteritis or have been found in stool samples as new molecular tools have been developed. In this work, a DNA microarray capable in theory of parallel detection of more than 100 viral species was developed and tested. Initial validation was done with 10 different virus species, and an additional 5 species were validated using clinical samples. Detection limits of 1 × 10(3) virus particles of Human adenovirus C (HAdV), Human astrovirus (HAstV), and group A Rotavirus (RV-A) were established. Furthermore, when exogenous RNA was added, the limit for RV-A detection decreased by one log. In a small group of clinical samples from children with gastroenteritis (n = 76), the microarray detected at least one viral species in 92% of the samples. Single infection was identified in 63 samples (83%), and coinfection with more than one virus was identified in 7 samples (9%). The most abundant virus species were RV-A (58%), followed by Anellovirus (15.8%), HAstV (6.6%), HAdV (5.3%), Norwalk virus (6.6%), Human enterovirus (HEV) (9.2%), Human parechovirus (1.3%), Sapporo virus (1.3%), and Human bocavirus (1.3%). To further test the specificity and sensitivity of the microarray, the results were verified by reverse transcription-PCR (RT-PCR) detection of 5 gastrointestinal viruses. The RT-PCR assay detected a virus in 59 samples (78%). The microarray showed good performance for detection of RV-A, HAstV, and calicivirus, while the sensitivity for HAdV and HEV was low. Furthermore, some discrepancies in detection of mixed infections were observed and were addressed by reverse transcription-quantitative PCR (RT-qPCR) of the viruses involved. It was observed that differences in the amount of genetic material favored the detection of the most abundant virus. The microarray described in this work should help in understanding the etiology of gastroenteritis in humans and animals.

Discovery of a novel polyomavirus in acute diarrheal samples from children.

Polyomaviruses are small circular DNA viruses associated with chronic infections and tumors in both human and animal hosts. Using an unbiased deep sequencing approach, we identified a novel, highly divergent polyomavirus, provisionally named MX polyomavirus (MXPyV), in stool samples from children. The ∼5.0 kB viral genome exhibits little overall homology (<46% amino acid identity) to known polyomaviruses, and, due to phylogenetic variation among its individual proteins, cannot be placed in any existing taxonomic group. PCR-based screening detected MXPyV in 28 of 834 (3.4%) fecal samples collected from California, Mexico, and Chile, and 1 of 136 (0.74%) of respiratory samples from Mexico, but not in blood or urine samples from immunocompromised patients. By quantitative PCR, the measured titers of MXPyV in human stool at 10% (weight/volume) were as high as 15,075 copies. No association was found between the presence of MXPyV and diarrhea, although girls were more likely to shed MXPyV in the stool than boys (p=0.012). In one child, viral shedding was observed in two stools obtained 91 days apart, raising the possibility of chronic infection by MXPyV. A multiple sequence alignment revealed that MXPyV is a closely related variant of the recently reported MWPyV and HPyV10 polyomaviruses. Further studies will be important to determine the association, if any, of MXPyV with disease in humans.