Mangalathu S. Rajeevan

Most human papillomavirus (HPV) surveillance studies target 30–50 of the more than 200 known types. We applied our recently described enriched whole-genome sequencing (eWGS) assay to demonstrate the impact of detecting all known and novel HPV types in male genital samples (n = 50). HPV was detected in nearly all (82%) samples, (mean number of types/samples 13.6; range 1–85), and nearly all HPV-positive samples included types in multiple genera (88%). A total of 560 HPV detections (237 unique HPV types: 46 alpha, 55 beta, 135 gamma, and 1 mu types) were made. The most frequently detected HPV types were alpha (HPV90, 43, and 74), beta (HPV115, 195, and 120), and gamma (HPV134, mSD2, and HPV50). High-risk alpha types (HPV16, 18, 31, 39, 52, and 58) were not common. A novel gamma type was identified (now officially HPV229) along with 90 unclassified types. This pilot study demonstrates the utility of the eWGS assay for broad-spectrum type detection and suggests a significantly higher type diversity in males compared to females that warrants further study.

Human papillomavirus (HPV) causes cervical cancer among women and is associated with other anogenital cancers in men and women. Prophylactic particulate vaccines that are affordable, self-administered and efficacious could improve uptake of HPV vaccines world-wide. The goal of this research is to develop a microparticulate HPV16 vaccine for transdermal administration using AdminPatch^® and assess its immunogenicity in a pre-clinical mouse model. HPV16 microparticles were prepared using a biocompatible polymer and characterized in terms of size, zeta potential, encapsulation efficiency and microparticle yield. Scanning and transmission electron microscopy were conducted to confirm particle image and to visualize the conformation of HPV16 vaccine particles released from microparticle formulation. In vivo studies performed to evaluate the potential of the microparticulate vaccine initiated a robust and sustained immune response. HPV16 IgG antibodies were significantly elevated in the microparticle group compared to antigen solutions administered by the transdermal route. Results show significant expansion of CD4+, CD45R, CD27 and CD62L cell populations in the vaccinated mice group, indicating the high efficacy of the microparticulate vaccine when administered via transdermal route. The findings of this study call attention to the use of minimally invasive, pain-free routes to deliver vaccine.

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 10⁴ 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 10⁴ 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.

We recently developed a test based on the Agilent SureSelect target enrichment system capturing genomic fragments from 191 human papillomaviruses (HPV) types for Illumina sequencing. This enriched whole genome sequencing (eWGS) assay provides an approach to identify all HPV types in a sample. Here we present a machine learning algorithm that calls HPV types based on the eWGS output. The algorithm based on the support vector machine (SVM) technique was trained on eWGS data from 122 control samples with known HPV types. The new algorithm demonstrated good performance in HPV type detection for designed samples with 25 or greater HPV plasmid copies per sample. We compared the results of HPV typing made by the new algorithm for 261 residual epidemiologic samples with the results of the typing delivered by the standard HPV Linear Array (LA). The agreement between methods (97.4%) was substantial (kappa = 0.783). However, the new algorithm identified additionally 428 instances of HPV types not detectable by the LA assay by design. Overall, we have demonstrated that the bioinformatics pipeline is an accurate tool for calling HPV types by analyzing data generated by eWGS processing of DNA fragments extracted from control and epidemiological samples.

We recently described a method for unbiased detection of all known human papillomaviruses (HPV) types with the potential for the determination of their variant and integration from the resulting whole genome sequence data. Considering the complex workflow for target-enriched next generation sequencing (NGS), we focused on the reproducibility and limit of detection (LOD) of this new universal HPV typing assay in this study. We evaluated the reproducibility and LOD for HPV genotyping based on our recently published method that used RNA-baits targeting whole genomes of 191 HPV types, Agilent SureSelect protocol for target enrichment and Illumina HiSeq 2500 for sequencing (eWGS, enriched whole genome sequencing). Two libraries, prepared from pooled plasmids representing 9 vaccine HPV types at varying input (1–625 copies/reaction), were sequenced twice giving four replicates for evaluating reproducibility and LOD. eWGS showed high correlation in the number of reads mapped to HPV reference genomes between the two flow-cell lanes within (R2 = 1) and between experiments (R2 = 0.99). The number of mapped reads was positively correlated to copy number (β = 13.9, p < 0.0001). The limit of blank (LOB) could be calculated based on mapped reads to HPV types not included in each sample. HPV genotyping was reproducible for all 9 types at 625 copies using multiple cut-off criteria but LOD was 25 copies based on number of reads above LOB even when multiple types were present. eWGS showed no bias for HPV genotyping under single or multiple infection (p = 0.16–0.99). The universal eWGS method for HPV genotyping has sensitivity, competitive with widely used consensus PCR methods with reduced type competition, and with the potential for determination of variant and integration status. The protocol used in this study, using defined samples varying in complexity and copy number, analyzed in replicate and duplicate assays, is applicable to most WGS methods.

Studies of mRNA and miRNA expression profiling increasingly use stabilized whole blood. Commercial RNA extraction kits do not provide information about the simultaneous recovery of both mRNA and miRNA. This study evaluated yield, quality, integrity and representation of mRNA and miRNA from whole blood stabilized in Tempus tubes using three RNA extraction kits; two filter-based (Tempus and Norgen) and one bead-based (MagMax; manual vs. semi-automated, and with and without DNase treatment). All RNA extraction kits and methods resulted in similar yields of mRNA (total RNA yield, quality, integrity and representation) whereas there were differences in yields of miRNA. MagMax, either manual or semi-automated, with or without DNase treatment, yielded 1.6–2.2-fold more miRNA than Tempus and Norgen kits. In addition, MagMax and Norgen methods gave greater than 12-fold more and 3.3-fold less enrichment of specific miRNA targets, respectively, in comparison to Tempus extraction reagents. This study identified MagMax kit for simultaneous recovery of both mRNA and miRNA from whole blood collected in Tempus tubes.

Chronic fatigue syndrome (CFS), also known as myalgic encephalomyelitis (ME), is a severely debilitating condition of unknown etiology. The symptoms and risk factors of ME/CFS share features of accelerated aging implicated in several diseases. Using telomere length as a marker, this study was performed to test the hypothesis that ME/CFS is associated with accelerated aging. Participant (n = 639) data came from the follow-up time point of the Georgia CFS surveillance study. Using the 1994 CFS Research Case Definition with questionnaire-based subscale thresholds for fatigue, function, and symptoms, participants were classified into four illness groups: CFS if all criteria were met (n = 64), CFS-X if CFS with exclusionary conditions (n = 77), ISF (insufficient symptoms/fatigue) if only some criteria were met regardless of exclusionary conditions (n = 302), and NF (non-fatigued) if no criteria and no exclusionary conditions (n = 196). Relative telomere length (T/S ratio) was measured using DNA from whole blood and real-time PCR. General linear models were used to estimate the association of illness groups or T/S ratio with demographics, biological measures and covariates with significance set at p < 0.05. The mean T/S ratio differed significantly by illness group (p = 0.0017); the T/S ratios in CFS (0.90 ± 0.03) and ISF (0.94 ± 0.02) were each significantly lower than in NF (1.06 ± 0.04). Differences in T/S ratio by illness groups remained significant after adjustment for covariates of age, sex, body mass index, waist–hip ratio, post-exertional malaise and education attainment. Telomere length was shorter by 635, 254 and 424 base pairs in CFS, CFS-X and ISF, respectively, compared to NF. This shorter telomere length translates to roughly 10.1–20.5, 4.0–8.2 and 6.6–13.7 years of additional aging in CFS, CFS-X and ISF compared to NF respectively. Further, stratified analyses based on age and sex demonstrated that the association of ME/CFS with short telomeres is largely moderated by female subjects < 45 years old. This study found a significant association of ME/CFS with premature telomere attrition that is largely moderated by female subjects < 45 years old. Our results indicate that ME/CFS could be included in the list of conditions associated with accelerated aging. Further work is needed to evaluate the functional significance of accelerated aging in ME/CFS.

The current practice of using only a few strongly associated genetic markers in regression models results in generally low power in prediction or accounting for heritability of complex human traits. We illustrate here a Bayesian joint estimation of single nucleotide polymorphism (SNP) effects principle to improve prediction of phenotype status from pathway-focused sets of SNPs. Chronic fatigue syndrome (CFS), a complex disease of unknown etiology with no laboratory methods for diagnosis, was chosen to demonstrate the power of this Bayesian method. For CFS, such a genetic predictive model in combination with clinical evidence might lead to an earlier diagnosis than one based solely on clinical findings. One of our goals is to model disease status using Bayesian statistics which perform variable selection and parameter estimation simultaneously and which can induce the sparseness and smoothness of the SNP effects. Smoothness of the SNP effects is obtained by explicit modeling of the covariance structure of the SNP effects. The Bayesian model achieved perfect goodness of fit when tested within the sampled data. Tenfold cross-validation resulted in 80 % accuracy, one of the best so far for CFS in comparison to previous prediction models. Model reduction aspects were investigated in a computationally feasible manner. Additionally, genetic variation estimates provided by the model identified specific genetic markers for their biological role in the disease pathophysiology. This proof-of-principle study provides a powerful approach combining Bayesian methods, SNPs representing multiple pathways and rigorous case ascertainment for accurate genetic risk prediction modeling of complex diseases like CFS and other chronic diseases.

Whole blood is generally processed for plasma and peripheral blood mononuclear cells (PBMCs) from granulocytes/erythrocytes using gradient centrifugation of blood with Histopaue-Ficoll. After separation of plasma and PBMCs, the residual erythrocytes/granulocytes, a rich source of DNA, is often discarded along with the separation medium. In order to isolate DNA from the granulocytes, current methods require the removal of the separation medium and subsequent purification of granulocytes. This report provides a method for extracting DNA using the PAXgene Blood DNA kit from granulocytes without purifying them from the separation medium.