NGS read classification using AI.

Clinical metagenomics is a powerful diagnostic tool, as it offers an open view into all DNA in a patient's sample. This allows the detection of pathogens that would slip through the cracks of classical specific assays. However, due to this unspecific nature of metagenomic sequencing, a huge amount of unspecific data is generated during the sequencing itself and the diagnosis only takes place at the data analysis stage where relevant sequences are filtered out. Typically, this is done by comparison to reference databases. While this approach has been optimized over the past years and works well to detect pathogens that are represented in the used databases, a common challenge in analysing a metagenomic patient sample arises when no pathogen sequences are found: How to determine whether truly no evidence of a pathogen is present in the data or whether the pathogen's genome is simply absent from the database and the sequences in the dataset could thus not be classified? Here, we present a novel approach to this problem of detecting novel pathogens in metagenomic datasets by classifying the (segments of) proteins encoded by the sequences in the datasets. We train a neural network on the sequences of coding sequences, labeled by taxonomic domain, and use this neural network to predict the taxonomic classification of sequences that can not be classified by comparison to a reference database, thus facilitating the detection of potential novel pathogens.

Analysis of DNA Sequence Classification Using CNN and Hybrid Models.

In a general computational context for biomedical data analysis, DNA sequence classification is a crucial challenge. Several machine learning techniques have used to complete this task in recent years successfully. Identification and classification of viruses are essential to avoid an outbreak like COVID-19. Regardless, the feature selection process remains the most challenging aspect of the issue. The most commonly used representations worsen the case of high dimensionality, and sequences lack explicit features. It also helps in detecting the effect of viruses and drug design. In recent days, deep learning (DL) models can automatically extract the features from the input. In this work, we employed CNN, CNN-LSTM, and CNN-Bidirectional LSTM architectures using Label and K-mer encoding for DNA sequence classification. The models are evaluated on different classification metrics. From the experimental results, the CNN and CNN-Bidirectional LSTM with K-mer encoding offers high accuracy with 93.16% and 93.13%, respectively, on testing data.

NanoString Technology for Human Papillomavirus Typing.

High-throughput HPV typing assays with increased automation, faster turnaround and type-specific digital readout would facilitate studies monitoring the impact of HPV vaccination. We evaluated the NanoString nCounter® platform for detection and digital readout of 48 HPV types in a single reaction. NanoString (NS) used proprietary software to design CodeSets: type-specific probe pairs targeting 48 HPV types and the globin gene. We tested residual DNA extracts from epidemiologic specimens and defined samples (HPV plasmids at 10 to 104 copies/reaction) directly (No-PCR) as well as after L1 consensus PCR of 45 (PCR-45) or 15 cycles (PCR-15). Assay and interpretation followed NS recommendations. We evaluated analytic performance by comparing NanoString results for types included in prior assays: Roche Linear Array (LA) or HPV TypeSeq assay. No-PCR results on 40 samples showed good type-specific agreement with LA (k = 0.621) but sensitivity was 65% with lower limit of detection (LOD) at 104 plasmid copies. PCR-45 results showed almost perfect type-specific agreement with LA (k = 0.862), 82% sensitivity and LOD at 10 copies. PCR-15 results on 75 samples showed substantial type-specific agreement with LA (k = 0.796, 92% sensitivity) and TypeSeq (k = 0.777, 87% sensitivity), and LOD at 10 copies of plasmids. This proof-of-principle study demonstrates the efficacy of the NS platform with HPV CodeSet for type-specific detection using a low number of PCR cycles (PCR-15). Studies are in progress to evaluate assay reproducibility and analytic validation with a larger number of samples.

Determining a sampling regime for PCR detection of respiratory tract viral infection at coronial post-mortem examinations.

Death due to respiratory infection is commonly encountered at autopsy. With only one opportunity to obtain samples for identification of a causative agent, it is important to ensure that sampling regimes are optimized to provide the greatest detection, without the expense and redundancy that can arise from over-sampling. This study was performed retrospectively using data from Coronial autopsies over the period 2012-2019 from which swabs from the nasopharyngeal region, trachea and lung parenchyma, in addition to samples of lung tissue, had been submitted for multiplex PCR detection of respiratory pathogens. From 97 cases with all four samples, there were 24 with at least one positive result for viral infection. Some cases had multiple positive results and a total of 27 respiratory tract viruses were identified, of which rhinovirus, influenza A virus and respiratory syncytial virus were the most common. Seventeen of the 27 viral infections (63%) were identified in all four samples. However, in nearly all cases (96%) the nasopharyngeal swab detected the infective agent when the multiplex PCR panel had detected infection in any of the four sample types. A nasopharyngeal swab is considered to be an optimal sample for detection of respiratory tract viral infection. As the samples analyzed were acquired before the appearance of the COVID-19 virus, the applicability of this finding for COVID-19 screening is not established.

Synergistically enhanced colorimetric molecular detection using smart cup: a case for instrument-free HPV-associated cancer screening.

Rationale: Early and accurate detection of disease is crucial for its prevention, identification, and treatment. However, most of disease diagnostics is still limited in clinical laboratories due to the need of complicated instruments and professional personnel. Herein, we reported a smartphone-based synergistically enhanced colorimetric method for molecular diagnostics in our point of care (POC) smart cup platform. Methods: A disposable microfluidic chip was developed for colorimetric loop-mediated isothermal amplification (LAMP) detection of multiple HPV DNA in our POC smart cup platform. The colorimetric detection takes advantage of synergistic effect of PPi4- and H+ ions, two byproducts of LAMP reaction. Color signal of LAMP assay was recorded and analyzed by our custom Android app (dubbed "Hue Analyzer"). Results: Our method not only significantly improves colorimetric readout, but also provides a 10-fold increase in detection sensitivity. It has been successfully applied for HPV-associated cancer screening with spiked saliva and clinical swab samples. Conclusion: The proposed POC diagnostic platform is completely compatible with other nucleic acid biomarkers and has great potential for personalized health monitoring and disease prevention.

Incidence and Clearance of Anal High-risk Human Papillomavirus Infections and Their Determinants Over 5 Years Among Human Immunodeficiency Virus-negative Men Who Have Sex With Men.

BACKGROUND: We aimed to assess the incidence and clearance of anal high-risk human papillomavirus (hrHPV) infections and determinants thereof among human immunodeficiency virus (HIV)-negative men who have sex with men (MSM) over a period of up to 5 years. METHODS: From 2010 to 2015, HIV-negative MSM were followed every 6 months. Anal self-swabs were collected at inclusion and every 6 months thereafter, and were HPV genotyped using the SPF10-PCR DEIA/LiPA25-system-v1. Incidence rates (IRs) and clearance rates (CRs) of incident anal hrHPV infections were assessed by hrHPV type (types 16, 18, 31, 33, 45, 52, and 58). Determinants of transitions between uninfected and infected states were assessed by hrHPV type using a time-homogenous multi-state Markov model. RESULTS: This study included 713 HIV-negative MSM, with a median age of 37 years (interquartile range [IQR] 31-43) and a median number of study visits of 6 (IQR 2-7). The IRs of anal infections had a median of 5.2 per 100 person-years (range: 2.2-7.9) across types, with HPV16 having the highest IR. The CRs of incident anal hrHPV infections had a median of 53.7 per 100 person-years (range: 33.4-65.3) across types, with HPV16 having the lowest CR. Having had over 100 lifetime sex partners was significantly associated with incident anal hrHPV infections in multivariable analyses. CONCLUSIONS: The high incidence and low clearance rates of anal HPV16 infection, compared to other hrHPV types, is consistent with HPV16 being implicated in the large majority of anal cancer cases.

Discovery and characterization of novel trans-spliced products of human polyoma JC virus late transcripts from PML patients.

Although the human neurotropic polyomavirus, JC virus (JCV), was isolated almost a half century ago, understanding the molecular mechanisms governing its biology remains highly elusive. JCV infects oligodendrocytes and astrocytes in the central nervous system (CNS) and causes a rare fatal brain disease known as progressive multifocal leukoencephalopathy (PML) in immunocompromised individuals including AIDS. It has a small circular DNA genome (∼5 kb) and generates two primary transcripts from its early and late coding regions, producing several predicted alternatively spliced products mainly by cis-splicing. Here, we report the discovery and characterization of two novel open reading frames (ORF1 and ORF2) associated with JCV late transcripts, generated by an unusual splicing process called trans-splicing. These ORFs result from (i) the trans-splicing of two different lengths of the 5'-short coding region of VP1 between the coding regions of agnoprotein and VP2 after replacing the intron located between these two coding regions and (ii) frame-shifts occurring within the VP2 coding sequences terminated by a stop codon. ORF1 and ORF2 are capable of encoding 58 and 72 aa long proteins respectively and are expressed in infected cells and PML patients. Each ORF protein shares a common coding region with VP1 and has a unique coding sequence of their own. When the expression of the unique coding regions of ORFs is blocked by a stop codon insertion in the viral background, the mutant virus replicates less efficiently when compared to wild-type, suggesting that the newly discovered ORFs play critical roles in the JCV life cycle.

Assessment of the Usefulness of Multiplex Real-Time PCR Tests in the Diagnostic and Therapeutic Process of Pneumonia in Hospitalized Children: A Single-Center Experience.

The aim of the study was assessment of the usefulness of multiplex real-time PCR tests in the diagnostic and therapeutic process in children hospitalized due to pneumonia and burdened with comorbidities. Methods. The study group included 97 children hospitalized due to pneumonia at the Karol Jonscher Teaching Hospital in Poznan, in whom multiplex real-time PCR tests (FTD respiratory pathogens 33; fast-track diagnostics) were used. Results. Positive test results of the test were achieved in 74 patients (76.3%). The average age in the group was 56 months. Viruses were detected in 61 samples (82% of all positive results); bacterial factors were found in 29 samples (39% of all positive results). The presence of comorbidities was established in 90 children (92.78%). On the basis of the obtained results, 5 groups of patients were established: viral etiology of infection, 34 patients; bacterial etiology, 7 patients; mixed etiology, 23 patients; pneumocystis, 9 patients; and no etiology diagnosed, 24 patients. Conclusions. Our analysis demonstrated that the participation of viruses in causing severe lung infections is significant in children with comorbidities. Multiplex real-time PCR tests proved to be more useful in establishing the etiology of pneumonia in hospitalized children than the traditional microbiological examinations.

Bonafide, type-specific human papillomavirus persistence among HIV-positive pregnant women: predictive value for cytological abnormalities, a longitudinal cohort study.

This study investigated the rate of human papillomavirus (HPV) persistence, associated risk factors, and predictors of cytological alteration outcomes in a cohort of human immunodeficiency virus-infected pregnant women over an 18-month period. HPV was typed through L1 gene sequencing in cervical smears collected during gestation and at 12 months after delivery. Outcomes were defined as nonpersistence (clearance of the HPV in the 2nd sample), re-infection (detection of different types of HPV in the 2 samples), and type-specific HPV persistence (the same HPV type found in both samples). An unfavourable cytological outcome was considered when the second exam showed progression to squamous intraepithelial lesion or high squamous intraepithelial lesion. Ninety patients were studied. HPV DNA persistence occurred in 50% of the cases composed of type-specific persistence (30%) or re-infection (20%). A low CD4+T-cell count at entry was a risk factor for type-specific, re-infection, or HPV DNA persistence. The odds ratio (OR) was almost three times higher in the type-specific group when compared with the re-infection group (OR = 2.8; 95% confidence interval: 0.43-22.79). Our findings show that bonafide (type-specific) HPV persistence is a stronger predictor for the development of cytological abnormalities, highlighting the need for HPV typing as opposed to HPV DNA testing in the clinical setting.

Bonafide, type-specific human papillomavirus persistence among HIV-positive pregnant women: predictive value for cytological abnormalities, a longitudinal cohort study

This study investigated the rate of human papillomavirus (HPV) persistence, associated risk factors, and predictors of cytological alteration outcomes in a cohort of human immunodeficiency virus-infected pregnant women over an 18-month period. HPV was typed through L1 gene sequencing in cervical smears collected during gestation and at 12 months after delivery. Outcomes were defined as nonpersistence (clearance of the HPV in the 2nd sample), re-infection (detection of different types of HPV in the 2 samples), and type-specific HPV persistence (the same HPV type found in both samples). An unfavourable cytological outcome was considered when the second exam showed progression to squamous intraepithelial lesion or high squamous intraepithelial lesion. Ninety patients were studied. HPV DNA persistence occurred in 50% of the cases composed of type-specific persistence (30%) or re-infection (20%). A low CD4+T-cell count at entry was a risk factor for type-specific, re-infection, or HPV DNA persistence. The odds ratio (OR) was almost three times higher in the type-specific group when compared with the re-infection group (OR = 2.8; 95% confidence interval: 0.43-22.79). Our findings show that bonafide (type-specific) HPV persistence is a stronger predictor for the development of cytological abnormalities, highlighting the need for HPV typing as opposed to HPV DNA testing in the clinical setting.

Significance of "Not Detected but Amplified" Results by Real-Time PCR Method for HPV DNA Detection.

Human papillomavirus (HPV) infection is an important etiologic factor in cervical carcinogenesis. Various HPV DNA detection methods have been evaluated for clinicopathological level. For the specimens with normal cytological finding, discrepancies among the detection methods were frequently found and adequate interpretation can be difficult. 6,322 clinical specimens were submitted and evaluated for real-time PCR and Hybrid Capture 2 (HC2). 573 positive or "Not Detected but Amplified" (NDBA) specimens by real-time PCR were additionally tested using genetic analyzer. For the reliability of real-time PCR, 325 retests were performed. Optimal cut-off cycle threshold (CT ) value was evaluated also. 78.7% of submitted specimens showed normal or nonspecific cytological finding. The distributions of HPV types by real-time PCR were not different between positive and NDBA cases. For positive cases by fragment analysis, concordance rates with real-time PCR and HC2 were 94.2% and 84.2%. In NDBA cases, fragment analysis and real-time PCR showed identical results in 77.0% and HC2 revealed 27.6% of concordance with fragment analysis. Optimal cut-off CT value was different for HPV types. NDBA results in real-time PCR should be regarded as equivocal, not negative. The adjustment of cut-off CT value for HPV types will be helpful for the appropriate result interpretation.

A Population Genetics Perspective on Geminivirus Infection.

UNLABELLED: The selective accumulation of both DNA components of a bipartite geminivirus, Abutilon mosaic virus, was recorded during early systemic infection of Nicotiana benthamiana plants. Purified nuclei were diagnosed for viral DNA using hybridization specific for DNA A or DNA B to detect these individual genome components either alone or both simultaneously by dual-color staining. Although this virus needs both components for symptomatic infection, DNA A alone was transported to upper leaves, where it was imported into phloem nuclei and replicated autonomously. The coinfection with DNA A and DNA B revealed an independent spread of both molecules, which resulted in a stochastic distribution of DNA A- and DNA A/B-infected nuclei. A population genetics evaluation of the respective frequencies was compared to a model computation. This elucidated a surprisingly simple relationship between the initial frequencies of the viral DNA components and the number of susceptible cells during the course of early systemic infection. IMPORTANCE: For bipartite begomoviruses, DNA B-independent long-distance spread of DNA A has been described before, but it has never been shown whether viral DNA A alone invades nuclei of systemic tissues and replicates therein. This is demonstrated now for the first time. During infection with DNA A and DNA B, a similar solitary spread of DNA A can be recognized at early stages. We describe a population genetics model of how the hit probabilities of DNA A and DNA B for susceptible cells determine the relative frequencies of either genome component during the course of infection.

Identification of several clades of novel single-stranded circular DNA viruses with conserved stem-loop structures in pig feces.

Metagenomic analysis of fecal samples collected from diarrheal swine detected sequences encoding a replication initiator protein (Rep). The genomes of ten novel single-stranded DNA viruses were determined, and they exhibited a similar genome organization. The two putative open reading frames (ORFs) encoding Rep and the capsid protein are bidirectionally transcribed and separated by two intergenic regions. Stem-loop structure(s) typical of genomes that undergo the rolling-circle DNA replication mechanism were observed. Phylogenetically, these ten genomes are in a monophyletic clade with the previously described porcine stool-associated virus (PoSCV) but are divergent enough to be further classified into to six distinct virus clades.

Development of a virus detection and discovery pipeline using next generation sequencing.

Next generation sequencing (NGS) has revolutionized virus discovery. Notwithstanding, a vertical pipeline, from sample preparation to data analysis, has not been available to the plant virology community. We developed a degenerate oligonucleotide primed RT-PCR method with multiple barcodes for NGS, and constructed VirFind, a bioinformatics tool specifically for virus detection and discovery able to: (i) map and filter out host reads, (ii) deliver files of virus reads with taxonomic information and corresponding Blastn and Blastx reports, and (iii) perform conserved domain search for reads of unknown origin. The pipeline was used to process more than 30 samples resulting in the detection of all viruses known to infect the processed samples, the extension of the genomic sequences of others, and the discovery of several novel viruses. VirFind was tested by four external users with datasets from plants or insects, demonstrating its potential as a universal virus detection and discovery tool.

False-positive results in metagenomic virus discovery: a strong case for follow-up diagnosis.

A viral metagenomic approach using virion enrichment, random amplification and next-generation sequencing was used to investigate an undiagnosed cluster of dairy cattle presenting with high persistent fever, unresponsive to anti-microbial and anti-inflammatory treatment, diarrhoea and redness of nose and teat. Serum and whole blood samples were taken in the predicted hyperviraemic state of an animal that a few days later presented with these clinical signs. Bioinformatics analysis of the resulting data from the DNA virus identification workflow (a total of 32 757 sequences with average read length 335 bases) initially demonstrated the presence of parvovirus-like sequences in the tested blood sample. Thorough follow-up using specific real-time RT-PCR assays targeting the detected sequence fragments confirmed the presence of these sequences in the original sample as well as in a sample of an additional animal, but a contamination with an identical genetic signature in negative extraction controls was demonstrated. Further investigation using an alternative extraction method identified a contamination of the originally used Qiagen extraction columns with parvovirus-like nucleic acids or virus particles. Although we did not find any relevant virus that could be associated with the disease, these observations clearly illustrate the importance of using a proper control strategy and follow-up diagnostic tests in any viral metagenomic study.

Characterization of human papillomavirus type 154 and tissue tropism of gammapapillomaviruses.

The novel human papillomavirus type 154 (HPV154) was characterized from a wart on the crena ani of a three-year-old boy. It was previously designated as the putative HPV type FADI3 by sequencing of a subgenomic FAP amplicon. We obtained the complete genome by combined methods including rolling circle amplification (RCA), genome walking through an adapted method for detection of integrated papillomavirus sequences by ligation-mediated PCR (DIPS-PCR), long-range PCR, and finally by cloning of four overlapping amplicons. Phylogenetically, the HPV154 genome clustered together with members of the proposed species Gammapapillomavirus 11, and demonstrated the highest identity in L1 to HPV136 (68.6%). The HPV154 was detected in 3% (2/62) of forehead skin swabs from healthy children. In addition, the different detection sites of 62 gammapapillomaviruses were summarized in order to analyze their tissue tropism. Several of these HPV types have been detected from multiple sources such as skin, oral, nasal, and genital sites, suggesting that the gammapapillomaviruses are generalists with a broader tissue tropism than previously appreciated. The study expands current knowledge concerning genetic diversity and tropism among HPV types in the rapidly growing gammapapillomavirus genus.

Mutability dynamics of an emergent single stranded DNA virus in a naïve host.

Quasispecies variants and recombination were studied longitudinally in an emergent outbreak of beak and feather disease virus (BFDV) infection in the orange-bellied parrot (Neophema chrysogaster). Detailed health monitoring and the small population size (<300 individuals) of this critically endangered bird provided an opportunity to longitudinally track viral replication and mutation events occurring in a circular, single-stranded DNA virus over a period of four years within a novel bottleneck population. Optimized PCR was used with different combinations of primers, primer walking, direct amplicon sequencing and sequencing of cloned amplicons to analyze BFDV genome variants. Analysis of complete viral genomes (n = 16) and Rep gene sequences (n = 35) revealed that the outbreak was associated with mutations in functionally important regions of the normally conserved Rep gene and immunogenic capsid (Cap) gene with a high evolutionary rate (3.41×10(-3) subs/site/year) approaching that for RNA viruses; simultaneously we observed significant evidence of recombination hotspots between two distinct progenitor genotypes within orange-bellied parrots indicating early cross-transmission of BFDV in the population. Multiple quasispecies variants were also demonstrated with at least 13 genotypic variants identified in four different individual birds, with one containing up to seven genetic variants. Preferential PCR amplification of variants was also detected. Our findings suggest that the high degree of genetic variation within the BFDV species as a whole is reflected in evolutionary dynamics within individually infected birds as quasispecies variation, particularly when BFDV jumps from one host species to another.

Characterization of a novel enterovirus serotype and an enterovirus EV-B93 isolated from acute flaccid paralysis patients.

Non-polio enteroviruses (NPEVs) are among the most common viruses infecting humans worldwide. Most of these infections are asymptomatic but few can lead to systemic and neurological disorders like Acute Flaccid Paralysis (AFP). Acute Flaccid Paralysis is a clinical syndrome and NPEVs have been isolated frequently from the patients suffering from AFP but little is known about their causal relationship. The objective of this study was to identify and characterize the NPEV serotypes recovered from 184 stool samples collected from AFP patients in Federally Administered Tribal Areas (FATA) in north-west of Pakistan. Overall, 44 (95.6 %) isolates were successfully typed through microneutralization assay as a member of enterovirus B species including echovirus (E)-2, E-3, E-4, E-6, E-7, E-11, E-13, E-14, E-21 and E-29 while two isolates (PAK NIH SP6545B and PAK NIH SP1202B) remained untypeable. The VP1 and capsid regions analysis characterized these viruses as EV-B93 and EV-B106. Phylogenetic analysis confirmed that PAK NIH isolates had high genetic diversity and represent distinct genotypes circulating in the country. Our findings highlight the role of NPEVs in AFP cases to be thoroughly investigated especially in high disease risk areas, with limited surveillance activities and health resources.

Distribution of hepatitis B virus subgenotype F2a in São Paulo, Brazil.

BACKGROUND: HBV genotype F is primarily found in indigenous populations from South America and is classified in four subgenotypes (F1 to F4). Subgenotype F2a is the most common in Brazil among genotype F cases. The aim of this study was to characterize HBV genotype F2a circulating in 16 patients from São Paulo, Brazil. Samples were collected between 2006 and 2012 and sent to Hospital Israelita Albert Einstein. A fragment of 1306 bp partially comprising HBsAg and DNA polymerase coding regions was amplified and sequenced. Viral sequences were genotyped by phylogenetic analysis using reference sequences from GenBank (n=198), including 80 classified as subgenotype F2a. Bayesian Markov chain Monte Carlo simulation implemented in BEAST v.1.5.4 was applied to obtain the best possible estimates using the model of nucleotide substitutions GTR+G+I. FINDINGS: It were identified three groups of sequences of subgenotype F2a: 1) 10 sequences from São Paulo state; 2) 3 sequences from Rio de Janeiro and one from São Paulo states; 3) 8 sequences from the West Amazon Basin. CONCLUSIONS: These results showing for the first time the distribution of F2a subgenotype in Brazil. The spreading and the dynamic of subgenotype F2a in Brazil requires the study of a higher number of samples from different regions as it is unfold in almost all Brazilian populations studied so far. We cannot infer with certainty the origin of these different groups due to the lack of available sequences. Nevertheless, our data suggest that the common origin of these groups probably occurred a long time ago.