Shu-Li Yang

Human papillomavirus-16 infection in advanced oral cavity cancer patients is related to an increased risk of distant metastases and poor survival.

Background Human papillomavirus (HPV) is an oncogenic virus causing oropharyngeal cancers and resulting in a favorable outcome after the treatment. The role of HPV in oral cavity squamous cell carcinoma (OSCC) remains ambiguous. Objective This study aimed to examine the effect of HPV infection on disease control among patients with OSCC following radical surgery with radiation-based adjuvant therapy. Patients and Method We prospectively followed 173 patients with advanced OSCC (96% were stage III/IV) who had undergone radical surgery and adjuvant therapy between 2004 and 2006. They were followed between surgery and death or up to 60 months. Surgical specimens were examined using a PCR-based HPV blot test. The primary endpoints were the risk of relapse and the time to relapse; the secondary endpoints were disease-free survival, disease-specific survival, and overall survival. Results The prevalence of HPV-positive OSCC was 22%; HPV-16 (9%) and HPV-18 (7%) were the genotypes most commonly encountered. Solitary HPV-16 infection was a poor predictor of 5-year distant metastases (hazard ratio, 3.4; 95% confidence interval, 1.4–8.0; P = 0.005), disease-free survival (P = 0.037), disease-specific survival (P = 0.006), and overall survival (P = 0.010), whereas HPV-18 infection had no impact on 5-year outcomes. The rate of 5-year distant metastases was significantly higher in the HPV-16 or level IV/V metastasis group compared with both the extracapsular spread or tumor depth ≥11-mm group and patients without risk factors (P<0.001). Conclusions HPV infections in advanced OSCC patients are not uncommon and clinically relevant. Compared with HPV-16-negative advanced OSCC patients, those with a single HPV-16 infection are at higher risk of distant metastases and poor survival despite undergoing radiation-based adjuvant therapy and require a more aggressive adjuvant treatment and a more thorough follow-up.

Molecular evolution and the global reemergence of enterovirus D68 by genome-wide analysis

Abstract Human enterovirus D68 (EV-D68) was first reported in the United States in 1962; thereafter, a few cases were reported from 1970 to 2005, but 2 outbreaks occurred in the Philippines (2008) and the United States (2014). However, little is known regarding the molecular evolution of this globally reemerging virus due to a lack of whole-genome sequences and analyses. Here, all publically available sequences including 147 full and 1248 partial genomes from GenBank were collected and compared at the clade and subclade level; 11 whole genomes isolated in Taiwan (TW) in 2014 were also added to the database. Phylogenetic trees were constructed to identify a new subclade, B3, and represent clade circulations among strains. Nucleotide sequence identities of the VP1 gene were 94% to 95% based on a comparison of subclade B3 to B1 and B2 and 87% to 91% when comparing A, C, and D. The patterns of clade circulation need to be clarified to improve global monitoring of EV-D68, even though this virus showed lower diversity among clades compared with the common enterovirus EV-71. Notably, severe cases isolated from Taiwan and China in 2014 were found in subclade B3. One severe case from Taiwan occurred in a female patient with underlying angioimmunoblastic T-cell lymphoma, from whom a bronchoalveolar lavage specimen was obtained. Although host factors play a key role in disease severity, we cannot exclude the possibility that EV-D68 may trigger clinical symptoms or death. To further investigate the genetic diversity of EV-D68, we reported 34 amino acid (aa) polymorphisms identified by comparing subclade B3 to B1 and B2. Clade D strains had a 1-aa deletion and a 2-aa insertion in the VP1 gene, and 1 of our TW/2014 strains had a shorter deletion in the 5′ untranslated region than a previously reported deletion. In summary, a new subclade, genetic indels, and polymorphisms in global strains were discovered elucidating evolutionary and epidemiological trends of EV-D68, and 11 genomes were added to the database. Virus variants may contribute to disease severity and clinical manifestations, and further studies are needed to investigate the associations between genetic diversity and clinical outcomes.

Increasing rates of low-risk human papillomavirus infections in patients with oral cavity squamous cell carcinoma: Association with clinical outcomes

BACKGROUND Although human papillomavirus (HPV) infections have been causally linked to oral cavity squamous cell carcinoma (OSCC), the potential role of low-risk HPV (LR-HPV) types in the pathogenesis of this malignancy remains unclear. OBJECTIVES We sought to investigate the distribution of HPV genotypes and their prognostic significance in OSCC patients treated by radical surgery, either with or without adjuvant therapy. STUDY DESIGN We studied two non-overlapping OSCC cohorts for the periods 2005-2006 (2005 cohort, n = 204) and 2010-2011 (2010 cohort, n = 206). Paraffin-embedded tissue blocks were collected, and the HPV genotype was determined using PCR plus HPV blot tests. The primary study endpoint was the prevalence of HPV genotypes. The secondary endpoints were the 2-year therapeutic outcomes. RESULTS The overall prevalence of HPV infections did not differ significantly in the two study cohorts. However, the prevalence of LR-HPV was significantly higher in the 2010 cohort than in the 2005 cohort (p = 0.002). The overall prevalence of HPV infections was not significantly associated with the 2-year outcomes. However, multivariate analysis demonstrated that LR-HPV infection was a predictor of poor 2-year disease-free survival (p = 0.036, hazard ratio [HR] = 3.1), disease-specific survival (p = 0.014, HR = 3.8), and overall survival (p = 0.016, HR = 3.2) in the subgroups of OSCC patients with poor differentiation, pN2 lymph node metastases, or extracapsular spread (n = 150). CONCLUSIONS LR-HPV infections may have an important role in determining the clinical outcomes of certain OSCC patients bearing specific risk factors.

Genomic Signatures for Avian H7N9 Viruses Adapting to Humans

An avian influenza A H7N9 virus emerged in March 2013 and caused a remarkable number of human fatalities. Genome variability in these viruses may provide insights into host adaptability. We scanned over 140 genomes of the H7N9 viruses isolated from humans and identified 104 positions that exhibited seven or more amino acid substitutions. Approximately half of these substitutions were identified in the influenza ribonucleoprotein (RNP) complex. Although PB2 627K of the avian virus promotes replication in humans, 45 of the 147 investigated PB2 sequences retained the E signature at this position, which is an avian characteristic. We discovered 10 PB2 substitutions that covaried with K627E. An RNP activity assay showed that Q591K, D701N, and M535L restored the polymerase activity in human cells when 627K transformed to an avian-like E. Genomic analysis of the human-isolated avian influenza virus is crucial in assessing genome variability, because relationships between position-specific variations can be observed and explored. In this study, we observed alternative positions that can potentially compensate for PB2 627K, a well-known marker for cross-species infection. An RNP assay suggested Q591K, D701N, and M535L as potential markers for an H7N9 virus capable of infecting humans.

Human papillomavirus 16/18 E7 viral loads predict distant metastasis in oral cavity squamous cell carcinoma

BACKGROUND Human papillomaviruses (HPV) seem to be related to distant metastasis (DM) in advanced oral cavity squamous cell carcinoma (OSCC) patients. OBJECTIVES This study aimed to investigate whether high-risk HPV viral load may predict DM among OSCC patients and stratify patients for risk-adapted treatment. STUDY DESIGN Viral loads of E7 oncogenes for HPV 16/18 were measured by quantitative PCR tests in paraffin-embedded lesional specimens from 312 OSCC of which the HPV genotypes had been determined previously. Multivariable Cox regression analysis was used to identify the independent prognostic factors for 5-year DM and C statistics were further computed. RESULTS By multivariable analysis, high HPV 16 E7 viral load (≥15.0 copies/genome); high HPV 18 E7 viral load (≥15.0 copies/genome); pathological N2 status (pN2); tumor depth ≥11 mm; extracapsular spread (ECS); and level IV/V metastases were independent risk factors for DM. We further identified three prognostic groups. In the high-risk group (level IV/V metastases or high HPV 16/18 E7 viral load plus pN2, tumor depth ≥11 mm, or ECS), the 5-year distant metastasis rate was 74%. In the intermediate-risk group (high HPV 16/18 E7 viral load, pN2, tumor depth ≥11 mm, or ECS), the 5-year DM rate was 17%. Finally, the 5-year DM rate was 1% in the low-risk group (no risk factors). The value of the C statistics was 0.78. CONCLUSIONS Among OSCC patients, high HPV 16/18 E7 viral load identifies a small subgroup of patients at high-risk of 5-year DM and suggest the need for more intensive treatments and follow-up strategies.

A Next-Generation Sequencing Data Analysis Pipeline for Detecting Unknown Pathogens from Mixed Clinical Samples and Revealing Their Genetic Diversity

Forty-two cytopathic effect (CPE)-positive isolates were collected from 2008 to 2012. All isolates could not be identified for known viral pathogens by routine diagnostic assays. They were pooled into 8 groups of 5-6 isolates to reduce the sequencing cost. Next-generation sequencing (NGS) was conducted for each group of mixed samples, and the proposed data analysis pipeline was used to identify viral pathogens in these mixed samples. Polymerase chain reaction (PCR) or enzyme-linked immunosorbent assay (ELISA) was individually conducted for each of these 42 isolates depending on the predicted viral types in each group. Two isolates remained unknown after these tests. Moreover, iteration mapping was implemented for each of these 2 isolates, and predicted human parechovirus (HPeV) in both. In summary, our NGS pipeline detected the following viruses among the 42 isolates: 29 human rhinoviruses (HRVs), 10 HPeVs, 1 human adenovirus (HAdV), 1 echovirus and 1 rotavirus. We then focused on the 10 identified Taiwanese HPeVs because of their reported clinical significance over HRVs. Their genomes were assembled and their genetic diversity was explored. One novel 6-bp deletion was found in one HPeV-1 virus. In terms of nucleotide heterogeneity, 64 genetic variants were detected from these HPeVs using the mapped NGS reads. Most importantly, a recombination event was found between our HPeV-3 and a known HPeV-4 strain in the database. Similar event was detected in the other HPeV-3 strains in the same clade of the phylogenetic tree. These findings demonstrated that the proposed NGS data analysis pipeline identified unknown viruses from the mixed clinical samples, revealed their genetic identity and variants, and characterized their genetic features in terms of viral evolution.

Influenza A/B virus detection and influenza A virus subtyping with emphasis on the novel H7N9 virus by using multiplex real-time RT-PCR

Infections of the novel avian influenza A H7N9 virus cause severe respiratory diseases and death. In this study, to develop highly sensitive methods for differentially detecting the H7N9 virus, multiplex and singular real-time reverse transcription polymerase chain reaction (RT-PCR) assays were established and examined by targeting the H7 and N9 genes of the H7N9 virus. Furthermore, an additional multiplex assay combining previous real time RT-PCR designs was established to subtype the pandemic H1N1, H3, and H5 influenza viruses. Applying the proposed assay system to analyze 100 clinical specimens collected from respiratory infection cases identified influenza A viruses (pandemic H1N1 and H3) in 23 samples. It has been demonstrated that other common respiratory viruses will not be detected by using this platform.

Incidence of and factors associated with false positives in laboratory diagnosis of norovirus infection by amplification of the RNA-dependent RNA polymerase gene.

Background Conventional reverse transcription-polymerase chain reaction (RT-PCR) amplification of the RNA-dependent RNA polymerase (RdRp) gene remains a used method for the rapid detection of norovirus (NV) in clinical laboratories. The incidence of and factors associated with false positives in this assay have not been previously evaluated. Methods/Principal Findings After an NV outbreak caused by the GII.4 Sydney strain in 2012, we reanalysed 250 stool samples positive for NV by RdRp gene detection. True positives were confirmed in 154 (61.6%) samples by successful amplification and sequencing confirmation of the viral protein 1 gene. Of the remaining 96 samples that underwent RT-PCR for the RdRp gene, 34 samples yielded PCR products of the expected length. However, the sequences of the amplicons belonged to the human genome, with 91–97% matched nucleotide sequences, indicating false positives. Multivariate analysis of the clinical features of the patients identified a positive stool culture for bacteria (adjusted odds ratio [aOR] 9.07, 95% adjusted confidence interval [aCI] 2.17–37.92, P = .003) and the use of parenteral antibiotics (aOR 5.55, 95% aCI 1.21–24.73, P = .027) as significant and independent factors associated with false positives. Conclusion Conventional RT-PCR targeting the RdRp gene of NV can lead to false positives in patients with bacterial enterocolitis by incidental amplification of DNA from a human source.

Molecular and serologic markers of HPV 16 infection are associated with local recurrence in patients with oral cavity squamous cell carcinoma

Human papillomavirus (HPV) infections predict mortality in Taiwanese patients with oral cavity squamous cell carcinoma (OCSCC). To address their prognostic significance for local recurrence (LR), in this retrospective cohort study we investigated different serologic and molecular markers of HPV 16 infection in 85 consecutive patients with primary OCSCC who received standard treatment and had their sera stored before treatment. Resected tumor specimens were examined with PCR-based assays for HPV 16 E6/E7 mRNA expression. Sera were tested with suspension arrays for the presence of HPV-specific antibodies using synthetic L1 and E6 peptides as well as a synthetic E7 protein. HPV 16 E6/E7 mRNA, anti-L1, anti-E6, and anti-E7 antibodies tested positive in 12%, 25%, 38%, and 41% of the study patients, respectively. Multivariate analysis identified pathological T3/T4, E6/E7 mRNA, and anti-E7 antibodies as independent risk factors for LR, whereas anti-E6 antibodies were an independent protective factor. In patients with ≥ 3 (high-risk group), 2 (intermediate-risk), and ≤ 1 (low-risk) independent risk factors (predictors), the 5-year LR rates were 75%, 42%, and 4%, respectively. Results were validated in an independent cohort. Together, our preliminary data indicate that HPV 16 infections as well as low and high serum levels of anti-E6 and anti-E7 antibodies, respectively, can serve as biomarkers of LR in patients with OCSCC, whereas the clinical usefulness of anti-HPV 16 antibodies for risk stratification of newly diagnosed cases deserves further scrutiny.

Erratum to: A metagenomics study for the identification of respiratory viruses in mixed clinical specimens: an application of the iterative mapping approach

Metagenomic approaches to detect viral genomes and variants in clinical samples have various challenges, including low viral titers and bacterial and human genome contamination. To address these limitations, we examined a next-generation sequencing (NGS) and iterative mapping approach for virus detection in clinical samples. We analyzed 40 clinical specimens from hospitalized children diagnosed with acute bronchiolitis, croup, or respiratory tract infections in which virus identification by viral culture or polymerase chain reaction (PCR) was unsuccessful. For our NGS data analysis pipeline, clinical samples were pooled into two NGS groups to reduce sequencing costs, and the depth and coverage of assembled contigs were effectively increased using an iterative mapping approach. PCR was individually performed for each specimen according to the NGS-predicted viral type. We successfully detected previously unidentified respiratory viruses in 26 of 40 specimens using our proposed NGS pipeline. Two dominant populations within the detected viruses were human rhinoviruses (HRVs; n = 14) and human coronavirus NL63 (n = 8), followed by human parainfluenza virus (HPIV), human parechovirus, influenza A virus, respiratory syncytial virus (RSV), and human metapneumovirus. This is the first study reporting the complete genome sequences of HRV-A101, HRV-C3, HPIV-4a, and RSV, as well as an analysis of their genetic variants, in Taiwan. These results demonstrate that this NGS pipeline allows to detect viruses which were not identified by routine diagnostic assays, directly from clinical samples.