Streamlined Whole-Genome Sequencing of Mumps Virus for High-Resolution Outbreak Analysis.

Since 2015, the United States has experienced a resurgence in the number of mumps cases and outbreaks in fully vaccinated populations. These outbreaks have occurred predominantly in close-quarter settings, such as camps, colleges, and detention centers. Phylogenetic analysis of 758 mumps-positive samples from outbreaks across the United States identified 743 (98%) as genotype G based on sequence analysis of the mumps small hydrophobic (SH) gene. Additionally, SH sequences in the genotype G samples showed almost no sequence diversity, with 675 (91%) of them having identical sequences or only one nucleotide difference. This uniformity of circulating genotype and strain created complications for epidemiologic investigations and necessitated the development of a system for rapidly generating mumps whole-genome sequences for more detailed analysis. In this study, we report a novel and streamlined assay for whole-genome sequencing (WGS) of mumps virus genotype G. The WGS procedure successfully generated 318 high-quality WGS sequences on nucleic acid from genotype G-positive respiratory samples collected during several mumps outbreaks in the United States between 2016 and 2019. Sequencing was performed by a rapid and highly sensitive custom Ion AmpliSeq mumps genotype G panel, with sample preparation performed on an Ion Chef and sequencing on an Ion S5. The WGS data generated by the AmpliSeq panel provided enhanced genomic resolution for epidemiological outbreak investigations. Translation and protein sequence analysis also identified several potentially important epitope changes in the circulating mumps genotype G strains compared to the Jeryl-Lynn strain (JL5) used in vaccines in the United States, which could explain the current level of vaccine escapes.

Whole-Genome Sequencing of SARS-CoV-2: Assessment of the Ion Torrent AmpliSeq Panel and Comparison with the Illumina MiSeq ARTIC Protocol.

Fast and effective methods are needed for sequencing of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome to track genetic mutations and to identify new and emerging variants during the ongoing pandemic. The objectives were to assess the performance of the SARS-CoV-2 AmpliSeq research panel and S5 plug-in analysis tools for whole-genome sequencing analysis of SARS-CoV-2 and to compare the results with those obtained with the MiSeq-based ARTIC analysis pipeline, using metrics such as depth, coverage, and concordance of single-nucleotide variant (SNV) calls. A total of 191 clinical specimens and a single cultured isolate were extracted and sequenced with AmpliSeq technology and analysis tools. Of the 191 clinical specimens, 83 (with threshold cycle [CT] values of 15.58 to 32.54) were also sequenced using an Illumina MiSeq-based method with the ARTIC analysis pipeline, for direct comparison. A total of 176 of the 191 clinical specimens sequenced on the S5XL system and prepared using the SARS-CoV-2 research panel had nearly complete coverage (>98%) of the viral genome, with an average depth of 5,031×. Similar coverage levels (>98%) were observed for 81/83 primary specimens that were sequenced with both methods tested. The sample with the lowest viral load (CT value of 32.54) achieved 89% coverage using the MiSeq method and failed to sequence with the AmpliSeq method. Consensus sequences produced by each method were identical for 81/82 samples in areas of equal coverage, with a single difference present in one sample. The AmpliSeq approach is as effective as the Illumina-based method using ARTIC v3 amplification for sequencing SARS-CoV-2 directly from patient specimens across a range of viral loads (CT values of 15.56 to 32.54 [median, 22.18]). The AmpliSeq workflow is very easily automated with the Ion Chef and S5 instruments and requires less training and experience with next-generation sequencing sample preparation than the Illumina workflow.

IL-13-regulated Macrophage Polarization during Granuloma Formation in an In Vitro Human Sarcoidosis Model.

The mechanisms underlying abnormal granuloma formation in patients with sarcoidosis are complex and remain poorly understood. A novel in vitro human granuloma model was used to determine the molecular mechanisms of granuloma genesis in patients with sarcoidosis in response to putative disease-causing mycobacterial antigens. Peripheral blood mononuclear cells (PBMCs) from patients with active sarcoidosis and from normal, disease-free control subjects were incubated for 7 days with purified protein derivative-coated polystyrene beads. Molecular responses, as reflected by differential expression of genes, extracellular cytokine patterns, and cell surface receptor expression, were analyzed. Unbiased systems biology approaches were used to identify signaling pathways engaged during granuloma formation. Model findings were compared with human lung and mediastinal lymph node gene expression profiles. Compared with identically treated PBMCs of control subjects (n = 5), purified protein derivative-treated sarcoidosis PBMCs (n = 6) were distinguished by the formation of cellular aggregates resembling granulomas. Ingenuity Pathway Analysis of differential expression gene patterns identified molecular pathways that are primarily regulated by IL-13, which promotes alternatively activated (M2) macrophage polarization. M2 polarization was further demonstrated by immunohistochemistry performed on the in vitro sarcoidosis granuloma-like structures. IL-13-regulated gene pathways were confirmed in human sarcoidosis lung and mediastinal lymph node tissues. The in vitro human sarcoidosis granuloma model provides novel insights into early granuloma formation, particularly IL-13 regulation of molecular networks that regulate M2 macrophage polarization. M2 macrophages are predisposed to aggregation and multinucleated giant cell formation, which are characteristic features of sarcoidosis granulomas. Clinical trial registered with www.clinicaltrials.gov (NCT01857401).

Isocitrate dehydrogenase 1 mutations in melanoma frequently co-occur with NRAS mutations.

AIMS: Isocitrate dehydrogenase 1 (IDH1) is a metabolic enzyme that converts isocitrate to α-ketoglutarate. IDH1 mutations are associated with the accumulation of the oncometabolite D-2-hydroxyglutarate, which acts as an epigenetic modifier, and the development of multiple malignancies. METHODS AND RESULTS: From May 2013 to June 2017, 252 melanoma samples from 214 patients with advanced or distant metastatic disease were tested for somatic mutations with the 50-gene AmpliSeq version 2 Cancer Hotspot Panel. Two hundred and twenty-six samples were sequenced successfully from 206 patients with 26 samples being characterised as quantity not sufficient. Melanomas from 10 separate patients (4.9%) were positive for IDH1 R132C (nine) or R132S (one). In six cases, the tumours also had a co-existing NRAS mutation (p.Q61R, Q61L and Q61K in two patients each) (P = 0.0044), whereas three patients had BRAF non-V600E mutations (V600K, V600G and V600R). Two cases had a TP53 variant, two cases an ATM variant, one a CDKN2A variant and one had an APC variant. The patients' ages ranged from 45 to 82 years (mean = 65.3, median = 65 years) and three of 10 patients were female (M:F ratio = 2:3). Three patients were stage 3 and seven were stage 4. Two are deceased, five are alive with stable disease (four on pembrolizumab) and three have no evidence of disease. CONCLUSION: IDH mutations may define a unique subset of melanoma patients who are eligible for IDH1 targeted therapies or combined therapies, such as MEK inhibitors when there is co-existing NRAS mutations, or immunotherapy.

Comparative tolerance of two massively parallel sequencing systems to common PCR inhibitors.

Human remains can be severely affected by the environment, and the DNA may be damaged, degraded, and/or inhibited. In this study, a DNA sample (at 1 ng DNA target input in triplicate) was spiked with five concentrations of five inhibitors (humic acid, melanin, hematin, collagen, and calcium) and sequenced with both the HID-Ion AmpliSeq™ Library Kit and ID panel on the Ion PGM™ System and the ForenSeq DNA Signature Prep Kit on the MiSeq FGx™. The objective of this study was to compare the baseline tolerance of the two sequencing chemistries and platforms to common inhibitors encountered in human remains recovered from missing person cases. The two chemistries generally were comparable but not always susceptible to the same inhibitors or at the same capacity. The HID-Ion AmpliSeq™ Library Kit and ID panel and the ForenSeq DNA Signature Prep Kit both were susceptible to humic acid, melanin, and collagen; however, the ForenSeq kit showed greater inhibition to melanin and collagen than the AmpliSeq™ kit. In contrast, the ForenSeq kit was resistant to the effects of hematin and calcium, whereas the AmpliSeq™ kit was highly inhibited by hematin. Short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs) showed the same trend among inhibitors when using the ForenSeq kit. Generally, locus read depth, heterozygote allele balance, and the numbers of alleles typed were inversely correlated with increasing inhibitor concentration. The larger STR loci were affected more so by the presence of inhibitors compared to smaller STR amplicons and SNP loci. Additionally, it does not appear that sequence noise is affected by the inhibitors. The noise percentage, however, does increase as the inhibitor concentration increases, due to the decrease in locus read depth and not likely because of chemistry effects.

A Novel In Vitro Human Granuloma Model of Sarcoidosis and Latent Tuberculosis Infection.

Many aspects of pathogenic granuloma formation are poorly understood, requiring new relevant laboratory models that represent the complexity (genetics and diversity) of human disease. To address this need, we developed an in vitro model of granuloma formation using human peripheral blood mononuclear cells (PBMCs) derived from patients with active sarcoidosis, latent tuberculosis (TB) infection (LTBI), or normal healthy control subjects. PBMCs were incubated for 7 days with uncoated polystyrene beads or beads coated with purified protein derivative (PPD) or human serum albumin. In response to PPD-coated beads, PBMCs from donors with sarcoidosis and LTBI formed robust multicellular aggregates resembling granulomas, displaying a typical T-helper cell type 1 immune response, as assessed by cytokine analyses. In contrast, minimal PBMC aggregation occurred when control PBMCs were incubated with PPD-coated beads, whereas the response to uncoated beads was negligible in all groups. Sarcoidosis PBMCs responded to human serum albumin-coated beads with modest cellular aggregation and inflammatory cytokine release. Whereas the granuloma-like aggregates formed in response to PPD-coated beads were similar for sarcoidosis and LTBI, molecular profiles differed significantly. mRNA expression patterns revealed distinct pathways engaged in early granuloma formation in sarcoidosis and LTBI, and they resemble molecular patterns reported in diseased human tissues. This novel in vitro human granuloma model is proposed as a tool to investigate mechanisms of early granuloma formation and for preclinical drug discovery research of human granulomatous disorders. Clinical trial registered with www.clinicaltrials.gov (NCT01857401).

Novel missense mutation in the bZIP transcription factor, MAF, associated with congenital cataract, developmental delay, seizures and hearing loss (Aymé-Gripp syndrome).

BACKGROUND: Cataract is a major cause of severe visual impairment in childhood. The purpose of this study was to determine the genetic cause of syndromic congenital cataract in an Australian mother and son. METHOD: Fifty-one genes associated with congenital cataract were sequenced in the proband using a custom Ampliseq library on the Ion Torrent Personal Genome Machine (PGM). Reads were aligned against the human genome (hg19) and variants were annotated. Variants were prioritised for validation by Sanger sequencing if they were novel, rare or previously reported to be associated with paediatric cataract and were predicted to be protein changing. Variants were assessed for segregation with the phenotype in the affected mother. RESULT: A novel likely pathogenic variant was identified in the transactivation domain of the MAF gene (c.176C > G, p.(Pro59Arg)) in the proband and his affected mother., but was absent in 326 unrelated controls and absent from public variant databases. CONCLUSION: The MAF variant is the likely cause of the congenital cataract, Asperger syndrome, seizures, hearing loss and facial characteristics in the proband, providinga diagnosis of Aymé-Gripp syndrome for the family.

Targeted next generation sequencing identifies novel NOTCH3 gene mutations in CADASIL diagnostics patients.

BACKGROUND: Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is a monogenic, hereditary, small vessel disease of the brain causing stroke and vascular dementia in adults. CADASIL has previously been shown to be caused by varying mutations in the NOTCH3 gene. The disorder is often misdiagnosed due to its significant clinical heterogeneic manifestation with familial hemiplegic migraine and several ataxia disorders as well as the location of the currently identified causative mutations. The aim of this study was to develop a new, comprehensive and efficient single assay strategy for complete molecular diagnosis of NOTCH3 mutations through the use of a custom next-generation sequencing (NGS) panel for improved routine clinical molecular diagnostic testing. RESULTS: Our custom NGS panel identified nine genetic variants in NOTCH3 (p.D139V, p.C183R, p.R332C, p.Y465C, p.C597W, p.R607H, p.E813E, p.C977G and p.Y1106C). Six mutations were stereotypical CADASIL mutations leading to an odd number of cysteine residues in one of the 34 NOTCH3 gene epidermal growth factor (EGF)-like repeats, including three new typical cysteine mutations identified in exon 11 (p.C597W; c.1791C>G); exon 18 (p.C977G; c.2929T>G) and exon 20 (p.Y1106C; c.3317A>G). Interestingly, a novel missense mutation in the CACNA1A gene was also identified in one CADASIL patient. All variants identified (novel and known) were further investigated using in silico bioinformatic analyses and confirmed through Sanger sequencing. CONCLUSIONS: NGS provides an improved and effective methodology for the diagnosis of CADASIL. The NGS approach reduced time and cost for comprehensive genetic diagnosis, placing genetic diagnostic testing within reach of more patients.

Profiling cancer gene mutations in clinical formalin-fixed, paraffin-embedded colorectal tumor specimens using targeted next-generation sequencing.

PURPOSE: The success of precision oncology relies on accurate and sensitive molecular profiling. The Ion AmpliSeq Cancer Panel, a targeted enrichment method for next-generation sequencing (NGS) using the Ion Torrent platform, provides a fast, easy, and cost-effective sequencing workflow for detecting genomic "hotspot" regions that are frequently mutated in human cancer genes. Most recently, the U.K. has launched the AmpliSeq sequencing test in its National Health Service. This study aimed to evaluate the clinical application of the AmpliSeq methodology. METHODS: We used 10 ng of genomic DNA from formalin-fixed, paraffin-embedded human colorectal cancer (CRC) tumor specimens to sequence 46 cancer genes using the AmpliSeq platform. In a validation study, we developed an orthogonal NGS-based resequencing approach (SimpliSeq) to assess the AmpliSeq variant calls. RESULTS: Validated mutational analyses revealed that AmpliSeq was effective in profiling gene mutations, and that the method correctly pinpointed "true-positive" gene mutations with variant frequency >5% and demonstrated high-level molecular heterogeneity in CRC. However, AmpliSeq enrichment and NGS also produced several recurrent "false-positive" calls in clinically druggable oncogenes such as PIK3CA. CONCLUSION: AmpliSeq provided highly sensitive and quantitative mutation detection for most of the genes on its cancer panel using limited DNA quantities from formalin-fixed, paraffin-embedded samples. For those genes with recurrent "false-positive" variant calls, caution should be used in data interpretation, and orthogonal verification of mutations is recommended for clinical decision making.

Alternative genome sequencing approaches of SARS-CoV-2 using Ion AmpliSeq Technology.

A thorough understanding of SARS-CoV-2 genetic features is compulsory to track the ongoing pandemic across multiple geographical locations of the world. Thermo Fisher Scientific USA has developed the Ion AmpliSeq SARS-CoV-2 Research Panel for the targeted sequencing of SARS-CoV-2 complete genome with high coverage and lower error rate. In this study an alternative approach of complete genome sequencing has been validated using different commercial sequencing kits to sequence the SARS-CoV-2. Amplification of cDNA with the SARS-CoV-2 primer pool was performed separately using two different master mixes: 2X environmental master mix (EM) and Platinum™ PCR SuperMix High Fidelity master mix (PM) instead of 5X Ion AmpliSeq™ HiFi Mix whereas NEBNext® Fast DNA Library Prep Set for Ion Torrent™ kit was used as an alternative to Ion AmpliSeq Library Kit Plus for other reagents. This study demonstrated a successful procedure to sequence the SARS-CoV-2 whole genome with average ∼2351 depth and 98.1% of total the reads aligned against the reference sequence (SARS-CoV-2, isolate Wuhan-Hu-1, complete genome). Although genome coverage varied, complete genomes were retrieved for both reagent sets with a reduced cost. This study proposed an alternative approach of high throughput sequencing using Ion torrent technology for the sequencing of SARS-CoV-2 in developing countries where sequencing facilities are low. This blended sequencing technique also offers a low cost protocol in developing countries like Bangladesh.

Gene expression changes in therapeutic ultrasound-treated venous leg ulcers.

Introduction: Low-frequency, low-intensity ultrasound has been previously shown to promote healing of chronic wounds in humans, but mechanisms behind these effects are poorly understood. The purpose of this study was to evaluate gene expression differences in debrided human venous ulcer tissue from patients treated with low-frequency (20 kHz), low-intensity (100 mW/cm2) ultrasound compared to a sham treatment in an effort to better understand the potential biological mechanisms. Methods: Debrided venous ulcer tissue was collected from 32 subjects one week after sham treatment or low-frequency, low-intensity ultrasound treatment. Of these samples, 7 samples (3 ultrasound treated and 4 sham treated) yielded sufficient quality total RNA for analysis by ultra-high multiplexed PCR (Ampliseq) and expression of more than 24,000 genes was analyzed. 477 genes were found to be significantly differentially expressed between the ultrasound and sham groups using cut-off values of p < 0.05 and fold change of 2. Results and Discussion: The top differentially expressed genes included those involved in regulation of cell metabolism, proliferation, and immune cell signaling. Gene set enrichment analysis identified 20 significantly enriched gene sets from upregulated genes and 4 significantly enriched gene sets from downregulated genes. Most of the enriched gene sets from upregulated genes were related to cell-cell signaling pathways. The most significantly enriched gene set from downregulated genes was the inflammatory response gene set. These findings show that therapeutic ultrasound influences cellular behavior in chronic wounds as early as 1 week after application. Considering the well-known role of chronic inflammation in impairing wound healing in chronic wounds, these results suggest that a downregulation of inflammatory genes is a possible biological mechanism of ultrasound-mediated venous chronic wound healing. Such increased understanding may ultimately lead to the enhancement of ultrasound devices to accelerate chronic wound healing and increase patient quality of life.

Rapid Identification of Lineage and Drug Resistance in Clinical Samples of Mycobacterium tuberculosis.

BACKGROUND: Mycobacterium tuberculosis is a slow-growing bacterium, which could delay its diagnosis and, therefore, promote the spread of the disease. Whole-genome sequencing allows us to obtain the complete drug-resistance profile of the strain; however, bacterial cultivation of clinical samples, along with complex processing, is required. METHODS: In this work, we explore AmpliSeq, an amplicon-based enrichment method for preparing libraries for targeted next-generation sequencing, to identify lineage and drug resistance directly from clinical samples. RESULTS: In our study, 111 clinical samples were tested. The lineage was identified in 100% of the culture-derived samples (52/52), in 95% of the smear (BK)-positive clinical samples (38/40) and in 42.1% of the BK-negative clinical samples (8/19). The drug-resistance profile was accurately identified in all but 11 samples, in which some phenotypic and genotypic discrepancies were found. In this respect, our panels were not exact in the detection of streptomycin resistance for isolates derived from clinical samples, as an extremely high number of SNPs in the rrs and rrl genes were detected due to cross-contamination. CONCLUSION: This technique has demonstrated high sensitivity in obtaining the drug-resistance profile of the isolates, as even those samples with DNA concentrations below the detection limit of Qubit produced a result. AmpliSeq technology is cheaper than whole-genome sequencing, easy to perform by laboratory technicians and applicable to any microorganism using the Ion Torrent platform.

Molecular Surveillance of Malaria Using the PF AmpliSeq Custom Assay for Plasmodium falciparum Parasites from Dried Blood Spot DNA Isolates from Peru.

Malaria molecular surveillance has great potential to support national malaria control programs (NMCPs), informing policy for its control and elimination. Here, we present a new three-day workflow for targeted resequencing of markers in 13 resistance-associated genes, histidine rich protein 2 and 3 (hrp2&3), a country (Peru)-specific 28 SNP-barcode for population genetic analysis, and apical membrane antigen 1 (ama1), using Illumina short-read sequencing technology. The assay applies a multiplex PCR approach to amplify all genomic regions of interest in a rapid and easily standardizable procedure and allows simultaneous amplification of a high number of targets at once, therefore having great potential for implementation into routine surveillance practice by NMCPs. The assay can be performed on routinely collected filter paper blood spots and can be easily adapted to different regions to investigate either regional trends or in-country epidemiological changes.

A custom hepatitis A virus assay for whole-genome sequencing.

Since 2016, the United States has experienced a resurgence in the number of hepatitis A virus (HAV) cases and outbreaks. These outbreaks have been sustained by person-to-person transmission with cases occurring predominantly in high-risk populations including intravenous drug users, individuals experiencing homelessness, and men who have sex with men. To investigate HAV transmission, a molecular-surveillance system consisting of real-time RT-PCR (rRT-PCR) for detection, and a conventional RT-PCR assay for genotyping of HAV, was established in New York State (NYS) in 2019. Since then, a total of 271 HAV-positive serum samples collected from cases across NYS between 2019 and 2021 were identified by rRT-PCR. To rapidly and efficiently generate HAV whole-genome sequences, a custom AmpliSeq™ panel was designed in collaboration with Thermo Fisher. To streamline the process, sample preparation was performed on an Ion Chef and sequencing on an Ion S5XL. Of the 271 HAV-positive samples, the whole-genome sequencing (WGS) assay successfully generated 134 near-complete, high-quality HAV sequences. Phylogenetic analysis of the VP1-2A region identified 216 IB, 48 IA, and 2 IIIA genotypes, while 5 were unable to be typed due to poor sequence in this key region. The HAV whole-genome sequencing approach provided a more efficient and streamlined approach for genotyping HAV compared to previous methods and resulted in phylogenetic trees with enhanced resolution compared to the HAV VP1-2A region alone.

Single nucleotide variants in Pseudomonas aeruginosa populations from sputum correlate with baseline lung function and predict disease progression in individuals with cystic fibrosis.

The severity and progression of lung disease are highly variable across individuals with cystic fibrosis (CF) and are imperfectly predicted by mutations in the human gene CFTR, lung microbiome variation or other clinical factors. The opportunistic pathogen Pseudomonas aeruginosa (Pa) dominates airway infections in most CF adults. Here we hypothesized that within-host genetic variation of Pa populations would be associated with lung disease severity. To quantify Pa genetic variation within CF sputum samples, we used deep amplicon sequencing (AmpliSeq) of 209 Pa genes previously associated with pathogenesis or adaptation to the CF lung. We trained machine learning models using Pa single nucleotide variants (SNVs), microbiome diversity data and clinical factors to classify lung disease severity at the time of sputum sampling, and to predict lung function decline after 5 years in a cohort of 54 adult CF patients with chronic Pa infection. Models using Pa SNVs alone classified lung disease severity with good sensitivity and specificity (area under the receiver operating characteristic curve: AUROC=0.87). Models were less predictive of lung function decline after 5 years (AUROC=0.74) but still significantly better than random. The addition of clinical data, but not sputum microbiome diversity data, yielded only modest improvements in classifying baseline lung function (AUROC=0.92) and predicting lung function decline (AUROC=0.79), suggesting that Pa AmpliSeq data account for most of the predictive value. Our work provides a proof of principle that Pa genetic variation in sputum tracks lung disease severity, moderately predicts lung function decline and could serve as a disease biomarker among CF patients with chronic Pa infections.

Target-enriched sequencing enables accurate identification of bloodstream infections in whole blood.

Bloodstream infections are within the top ten causes of death globally, with a mortality rate of up to 70%. Gold standard blood culture testing is time-consuming, resulting in delayed, but accurate, treatment. Molecular methods, such as RT-qPCR, have limited targets in one run. We present a new Ampliseq detection system (ADS) combining target amplification and next-generation sequencing for accurate identification of bacteria, fungi, and antimicrobial resistance determinants directly from blood samples. In this study, we included removal of human genomic DNA during nucleic acid extraction, optimized the target sequence set and drug resistance genes, performed antimicrobial resistance profiling of clinical isolates, and evaluated mock specimens and clinical samples by ADS. ADS successfully identified pathogens at the species-level in 36 h, from nucleic acid extraction to results. Besides pathogen identification, ADS can also present drug resistance profiles. ADS enabled detection of all bacteria and accurate identification of 47 pathogens. In 20 spiked samples and 8 clinical specimens, ADS detected at least 92.81% of reads mapped to pathogens. ADS also showed consistency with the three culture-negative samples, and correctly identified pathogens in four of five culture-positive clinical blood specimens. This Ampliseq-based technology promises broad coverage and accurate pathogen identification, helping clinicians to accurately diagnose and treat bloodstream infections.

Whole-genome sequencing of Kaposi sarcoma-associated herpesvirus (KSHV/HHV8) reveals evidence for two African lineages.

Kaposi sarcoma (KS)-associated herpesvirus (KSHV/HHV-8) was first sequenced from the body cavity (BC) lymphoma cell line, BC-1, in 1996. Few other KSHV genomes have been reported. Our knowledge of sequence variation for this virus remains spotty. This study reports additional genomes from historical US patient samples and from African KS biopsies. It describes an assay that spans regions of the virus that cannot be covered by short read sequencing. These include the terminal repeats, the LANA repeats, and the origins of replication. A phylogenetic analysis, based on 107 genomes, identified three distinct clades; one containing isolates from USA/Europe/Japan collected in the 1990s and two of Sub-Saharan Africa isolates collected since 2010. This analysis indicates that the KSHV strains circulating today differ from the isolates collected at the height of the AIDS epidemic. This analysis helps experimental designs and potential vaccine studies.

Targeted next generation sequencing of Cyclospora cayetanensis mitochondrial genomes from seeded fresh produce and other seeded food samples.

Cyclospora cayetanensis, a coccidian apicomplexan parasite, causes large outbreaks of foodborne diarrheal disease globally. Tracking the source of C. cayetanensis oocyst contamination in food items is essential to reduce, even prevent outbreaks. We previously showed that a genotyping method based on mitochondrial single nucleotide polymorphism (SNP) profiles had discriminatory power in classifying C. cayetanensis clinical isolates. In food specimens, low level contamination by oocysts and difficulties in DNA extraction present significant challenges in genotyping method development. Here, we report the development of a highly sensitive, custom-designed, targeted sequencing method based on the Illumina AmpliSeq platform; our method was capable of consistently generating near-complete mitochondrial genome sequences of C. cayetanensis from foods with low levels of contamination. To simulate environmentally observed contamination levels in foods, we seeded various food matrices, such as fresh produce and prepared dishes, with known quantities of oocysts, and isolated genomic DNA from washed food samples. Using the Ampliseq Targeted Sequencing method, we obtained near-complete mitochondrial genome sequences of C. cayetanensis from food samples seeded with as low as five to ten oocysts and used the data in downstream analysis. The flexibility of the AmpliSeq platform could potentially allow for more genomic targets to be added to achieve higher discriminatory power. This level of sensitivity in capturing high resolution genome data from contaminated food samples is a critical milestone towards the potential development of a comprehensive genotyping method for C. cayetanensis.

Application of a High-Throughput Targeted Sequence AmpliSeq Procedure to Assess the Presence and Variants of Virulence Genes in Salmonella.

We have developed a targeted, amplicon-based next-generation sequencing method to detect and analyze 227 virulence genes (VG) of Salmonella (AmpliSeqSalm_227VG) for assessing the pathogenicity potential of Salmonella. The procedure was developed using 80 reference genomes representing 75 epidemiologically-relevant serovars associated with human salmonellosis. We applied the AmpliSeqSalm_227VG assay to (a) 35 previously characterized field strains of Salmonella consisting of serovars commonly incriminated in foodborne illnesses and (b) 34 Salmonella strains with undisclosed serological or virulence attributes, and were able to divide Salmonella VGs into two groups: core VGs and variable VGs. The commonest serovars causing foodborne illnesses such as Enteritidis, Typhimurium, Heidelberg and Newport had a high number of VGs (217-227). In contrast, serovars of subspecies not commonly associated with human illnesses, such as houtenae, arizonae and salame, tended to have fewer VGs (177-195). Variable VGs were not only infrequent but, when present, displayed considerable sequence variation: safC, sseL, sseD, sseE, ssaK and stdB showed the highest variation and were linked to strain pathogenicity. In a chicken infection model, VGs belonging to rfb and sse operons showed differences and were linked with pathogenicity. The high-throughput, targeted NGS-based AmpliSeqSalm_227VG procedure provided previously unknown information about variation in select virulence genes that can now be applied to a much larger population of Salmonella for evaluating pathogenicity of various serovars of Salmonella and for risk assessment of foodborne salmonellosis.

Testing the Ion AmpliSeq™ HID Y-SNP Research Panel v1 for performance and resolution in admixed South Americans of haplogroup Q.

Y haplogroups, defined by Y-SNPs, allow the reconstruction of the human Y chromosome genealogy, which is important for population, evolutionary and forensic genetics. In this study, Y-SNPs were typed and haplogroups inferred with the MPS Ion AmpliSeq™ HID Y-SNP Research Panel v1, as a high-throughput approach. Firstly, the performance of the panel was evaluated with different DNA input amounts, reagent volumes and cycle numbers. DNA-inputs from 0.5 to 1 ng generated the most balanced read depth. Combined with full reagent and 19 cycles, this offered the highest number of amplicons with a sequencing read depth of at least 20 reads. Secondly, the sub-haplogroups of 182 admixed South Americans and Greenlanders belonging to haplogroup Q were inferred and tested for potential improvement in resolution. Most samples were assigned to lineage Q-M3 with some samples assigned to lineages upstream (Q-M346, L56, L57; Q-L331, L53; Q-L54; Q-CTS11969, CTS11970) or parallel (Q-L330, L334; Q-Z780/M971) to Q-M3. Only one sample was assigned to a downstream lineage (Q-Z35615, Z35616). Most individuals of haplogroup Q with NAM ancestry could neither be distinguished from each other, nor from half of the Greenlandic samples. Typing additional, known SNPs within lineage Q-M3, Z19483 and SA05, increased the resolution of predicted haplogroups. The search for novel variants in the sequenced regions allowed the detection of 42 variants and the subdivision of lineage Q-M3 into new subclades. The variants found in six of these subclades were exclusive to certain South American countries. In light of the limited differentiation of haplogroup Q samples, the additional information on known or novel SNPs disclosed in this study when using MPS Ion AmpliSeq™ HID Y-SNP Research Panel v1 should be included in the Yleaf software, to increase the differentiation of lineage Q-M3.