Development of a new barcode-based, multiplex-PCR, next-generation-sequencing assay and data processing and analytical pipeline for multiplicity of infection detection of Plasmodium falciparum.

BACKGROUND: Simultaneous infection with multiple malaria parasite strains is common in high transmission areas. Quantifying the number of strains per host, or the multiplicity of infection (MOI), provides additional parasite indices for assessing transmission levels but it is challenging to measure accurately with current tools. This paper presents new laboratory and analytical methods for estimating the MOI of Plasmodium falciparum. METHODS: Based on 24 single nucleotide polymorphisms (SNPs) previously identified as stable, unlinked targets across 12 of the 14 chromosomes within P. falciparum genome, three multiplex PCRs of short target regions and subsequent next generation sequencing (NGS) of the amplicons were developed. A bioinformatics pipeline including B4Screening pathway removed spurious amplicons to ensure consistent frequency calls at each SNP location, compiled amplicons by SNP site diversity, and performed algorithmic haplotype and strain reconstruction. The pipeline was validated by 108 samples generated from cultured-laboratory strain mixtures in different proportions and concentrations, with and without pre-amplification, and using whole blood and dried blood spots (DBS). The pipeline was applied to 273 smear-positive samples from surveys conducted in western Kenya, then providing results into StrainRecon Thresholding for Infection Multiplicity (STIM), a novel MOI estimator. RESULTS: The 24 barcode SNPs were successfully identified uniformly across the 12 chromosomes of P. falciparum in a sample using the pipeline. Pre-amplification and parasite concentration, while non-linearly associated with SNP read depth, did not influence the SNP frequency calls. Based on consistent SNP frequency calls at targeted locations, the algorithmic strain reconstruction for each laboratory-mixed sample had 98.5% accuracy in dominant strains. STIM detected up to 5 strains in field samples from western Kenya and showed declining MOI over time (q < 0.02), from 4.32 strains per infected person in 1996 to 4.01, 3.56 and 3.35 in 2001, 2007 and 2012, and a reduction in the proportion of samples with 5 strains from 57% in 1996 to 18% in 2012. CONCLUSION: The combined approach of new multiplex PCRs and NGS, the unique bioinformatics pipeline and STIM could identify 24 barcode SNPs of P. falciparum correctly and consistently. The methodology could be applied to field samples to reliably measure temporal changes in MOI.

Malignant Transformation of Hymenolepis nana in a Human Host.

Neoplasms occur naturally in invertebrates but are not known to develop in tapeworms. We observed nests of monomorphic, undifferentiated cells in samples from lymph-node and lung biopsies in a man infected with the human immunodeficiency virus (HIV). The morphologic features and invasive behavior of the cells were characteristic of cancer, but their small size suggested a nonhuman origin. A polymerase-chain-reaction (PCR) assay targeting eukaryotes identified Hymenolepis nana DNA. Although the cells were unrecognizable as tapeworm tissue, immunohistochemical staining and probe hybridization labeled the cells in situ. Comparative deep sequencing identified H. nana structural genomic variants that are compatible with mutations described in cancer. Invasion of human tissue by abnormal, proliferating, genetically altered tapeworm cells is a novel disease mechanism that links infection and cancer.

Spike gene deletion quasispecies in serum of patient with acute MERS-CoV infection.

The spike glycoprotein of the Middle East respiratory coronavirus (MERS-CoV) facilitates receptor binding and cell entry. During investigation of a multi-facility outbreak of MERS-CoV in Taif, Saudi Arabia, we identified a mixed population of wild-type and variant sequences with a large 530 nucleotide deletion in the spike gene from the serum of one patient. The out of frame deletion predicted loss of most of the S2 subunit of the spike protein leaving the S1 subunit with an intact receptor binding domain. This finding documents human infection with a novel genetic variant of MERS-CoV present as a quasispecies. J. Med. Virol. 89:542-545, 2017. © 2016 Wiley Periodicals, Inc.

Evolution of mitosome metabolism and invasion-related proteins in Cryptosporidium.

BACKGROUND: The switch from photosynthetic or predatory to parasitic life strategies by apicomplexans is accompanied with a reductive evolution of genomes and losses of metabolic capabilities. Cryptosporidium is an extreme example of reductive evolution among apicomplexans, with losses of both the mitosome genome and many metabolic pathways. Previous observations on reductive evolution were largely based on comparative studies of various groups of apicomplexans. In this study, we sequenced two divergent Cryptosporidium species and conducted a comparative genomic analysis to infer the reductive evolution of metabolic pathways and differential evolution of invasion-related proteins within the Cryptosporidium lineage. RESULTS: In energy metabolism, Cryptosporidium species differ from each other mostly in mitosome metabolic pathways. Compared with C. parvum and C. hominis, C. andersoni possesses more aerobic metabolism and a conventional electron transport chain, whereas C. ubiquitum has further reductions in ubiquinone and polyisprenoid biosynthesis and has lost both the conventional and alternative electron transport systems. For invasion-associated proteins, similar to C. hominis, a reduction in the number of genes encoding secreted MEDLE and insulinase-like proteins in the subtelomeric regions of chromosomes 5 and 6 was also observed in C. ubiquitum and C. andersoni, whereas mucin-type glycoproteins are highly divergent between the gastric C. andersoni and intestinal Cryptosporidium species. CONCLUSIONS: Results of the study suggest that rapidly evolving mitosome metabolism and secreted invasion-related proteins could be involved in tissue tropism and host specificity in Cryptosporidium spp. The finding of progressive reduction in mitosome metabolism among Cryptosporidium species improves our knowledge of organelle evolution within apicomplexans.

Comparative genomics reveals Cyclospora cayetanensis possesses coccidia-like metabolism and invasion components but unique surface antigens.

BACKGROUND: Cyclospora cayetanensis is an apicomplexan that causes diarrhea in humans. The investigation of foodborne outbreaks of cyclosporiasis has been hampered by a lack of genetic data and poor understanding of pathogen biology. In this study we sequenced the genome of C. cayetanensis and inferred its metabolism and invasion components based on comparative genomic analysis. RESULTS: The genome organization, metabolic capabilities and potential invasion mechanism of C. cayetanensis are very similar to those of Eimeria tenella. Propanoyl-CoA degradation, GPI anchor biosynthesis, and N-glycosylation are some apparent metabolic differences between C. cayetanensis and E. tenella. Unlike Eimeria spp., there are no active LTR-retrotransposons identified in C. cayetanensis. The similar repertoire of host cell invasion-related proteins possessed by all coccidia suggests that C. cayetanensis has an invasion process similar to the one in T. gondii and E. tenella. However, the significant reduction in the number of identifiable rhoptry protein kinases, phosphatases and serine protease inhibitors indicates that monoxenous coccidia, especially C. cayetanensis, have limited capabilities or use a different system to regulate host cell nuclear activities. C. cayetanensis does not possess any cluster of genes encoding the TA4-type SAG surface antigens seen in E. tenella, and may use a different family of surface antigens in initial host cell interactions. CONCLUSIONS: Our findings indicate that C. cayetanensis possesses coccidia-like metabolism and invasion components but unique surface antigens. Amino acid metabolism and post-translation modifications of proteins are some major differences between C. cayetanensis and other apicomplexans. The whole genome sequence data of C. cayetanensis improve our understanding of the biology and evolution of this major foodborne pathogen and facilitate the development of intervention measures and advanced diagnostic tools.

Multilocus Sequence Typing Tool for Cyclospora cayetanensis.

Because the lack of typing tools for Cyclospora cayetanensis has hampered outbreak investigations, we sequenced its genome and developed a genotyping tool. We observed 2 to 10 geographically segregated sequence types at each of 5 selected loci. This new tool could be useful for case linkage and infection/contamination source tracking.

Determination of Evolutionary Relationships of Outbreak-Associated Listeria monocytogenes Strains of Serotypes 1/2a and 1/2b by Whole-Genome Sequencing.

We used whole-genome sequencing to determine evolutionary relationships among 20 outbreak-associated clinical isolates of Listeria monocytogenes serotypes 1/2a and 1/2b. Isolates from 6 of 11 outbreaks fell outside the clonal groups or "epidemic clones" that have been previously associated with outbreaks, suggesting that epidemic potential may be widespread in L. monocytogenes and is not limited to the recognized epidemic clones. Pairwise comparisons between epidemiologically related isolates within clonal complexes showed that genome-level variation differed by 2 orders of magnitude between different comparisons, and the distribution of point mutations (core versus accessory genome) also varied. In addition, genetic divergence between one closely related pair of isolates from a single outbreak was driven primarily by changes in phage regions. The evolutionary analysis showed that the changes could be attributed to horizontal gene transfer; members of the diverse bacterial community found in the production facility could have served as the source of novel genetic material at some point in the production chain. The results raise the question of how to best utilize information contained within the accessory genome in outbreak investigations. The full magnitude and complexity of genetic changes revealed by genome sequencing could not be discerned from traditional subtyping methods, and the results demonstrate the challenges of interpreting genetic variation among isolates recovered from a single outbreak. Epidemiological information remains critical for proper interpretation of nucleotide and structural diversity among isolates recovered during outbreaks and will remain so until we understand more about how various population histories influence genetic variation.

Comparative genomic analysis reveals occurrence of genetic recombination in virulent Cryptosporidium hominis subtypes and telomeric gene duplications in Cryptosporidium parvum.

BACKGROUND: Cryptosporidium hominis is a dominant species for human cryptosporidiosis. Within the species, IbA10G2 is the most virulent subtype responsible for all C. hominis-associated outbreaks in Europe and Australia, and is a dominant outbreak subtype in the United States. In recent yearsIaA28R4 is becoming a major new subtype in the United States. In this study, we sequenced the genomes of two field specimens from each of the two subtypes and conducted a comparative genomic analysis of the obtained sequences with those from the only fully sequenced Cryptosporidium parvum genome. RESULTS: Altogether, 8.59-9.05 Mb of Cryptosporidium sequences in 45-767 assembled contigs were obtained from the four specimens, representing 94.36-99.47% coverage of the expected genome. These genomes had complete synteny in gene organization and 96.86-97.0% and 99.72-99.83% nucleotide sequence similarities to the published genomes of C. parvum and C. hominis, respectively. Several major insertions and deletions were seen between C. hominis and C. parvum genomes, involving mostly members of multicopy gene families near telomeres. The four C. hominis genomes were highly similar to each other and divergent from the reference IaA25R3 genome in some highly polymorphic regions. Major sequence differences among the four specimens sequenced in this study were in the 5' and 3' ends of chromosome 6 and the gp60 region, largely the result of genetic recombination. CONCLUSIONS: The sequence similarity among specimens of the two dominant outbreak subtypes and genetic recombination in chromosome 6, especially around the putative virulence determinant gp60 region, suggest that genetic recombination plays a potential role in the emergence of hyper-transmissible C. hominis subtypes. The high sequence conservation between C. parvum and C. hominis genomes and significant differences in copy numbers of MEDLE family secreted proteins and insulinase-like proteases indicate that telomeric gene duplications could potentially contribute to host expansion in C. parvum.

Chlamydia psittaci comparative genomics reveals intraspecies variations in the putative outer membrane and type III secretion system genes.

Chlamydia psittaci is an obligate intracellular bacterium that can cause significant disease among a broad range of hosts. In humans, this organism may cause psittacosis, a respiratory disease that can spread to involve multiple organs, and in rare untreated cases may be fatal. There are ten known genotypes based on sequencing the major outer-membrane protein gene, ompA, of C. psittaci. Each genotype has overlapping host preferences and virulence characteristics. Recent studies have compared C. psittaci among other members of the Chlamydiaceae family and showed that this species frequently switches hosts and has undergone multiple genomic rearrangements. In this study, we sequenced five genomes of C. psittaci strains representing four genotypes, A, B, D and E. Due to the known association of the type III secretion system (T3SS) and polymorphic outer-membrane proteins (Pmps) with host tropism and virulence potential, we performed a comparative analysis of these elements among these five strains along with a representative genome from each of the remaining six genotypes previously sequenced. We found significant genetic variation in the Pmps and tbl3SS genes that may partially explain differences noted in C. psittaci host infection and disease.

Isolation and enrichment of Cryptosporidium DNA and verification of DNA purity for whole-genome sequencing.

Whole-genome sequencing of Cryptosporidium spp. is hampered by difficulties in obtaining sufficient, highly pure genomic DNA from clinical specimens. In this study, we developed procedures for the isolation and enrichment of Cryptosporidium genomic DNA from fecal specimens and verification of DNA purity for whole-genome sequencing. The isolation and enrichment of genomic DNA were achieved by a combination of three oocyst purification steps and whole-genome amplification (WGA) of DNA from purified oocysts. Quantitative PCR (qPCR) analysis of WGA products was used as an initial quality assessment of amplified genomic DNA. The purity of WGA products was assessed by Sanger sequencing of cloned products. Next-generation sequencing tools were used in final evaluations of genome coverage and of the extent of contamination. Altogether, 24 fecal specimens of Cryptosporidium parvum, C. hominis, C. andersoni, C. ubiquitum, C. tyzzeri, and Cryptosporidium chipmunk genotype I were processed with the procedures. As expected, WGA products with low (<16.0) threshold cycle (CT) values yielded mostly Cryptosporidium sequences in Sanger sequencing. The cloning-sequencing analysis, however, showed significant contamination in 5 WGA products (proportion of positive colonies derived from Cryptosporidium genomic DNA, ≤25%). Following this strategy, 20 WGA products from six Cryptosporidium species or genotypes with low (mostly <14.0) CT values were submitted to whole-genome sequencing, generating sequence data covering 94.5% to 99.7% of Cryptosporidium genomes, with mostly minor contamination from bacterial, fungal, and host DNA. These results suggest that the described strategy can be used effectively for the isolation and enrichment of Cryptosporidium DNA from fecal specimens for whole-genome sequencing.

A distinct lineage of influenza A virus from bats.

Influenza A virus reservoirs in animals have provided novel genetic elements leading to the emergence of global pandemics in humans. Most influenza A viruses circulate in waterfowl, but those that infect mammalian hosts are thought to pose the greatest risk for zoonotic spread to humans and the generation of pandemic or panzootic viruses. We have identified an influenza A virus from little yellow-shouldered bats captured at two locations in Guatemala. It is significantly divergent from known influenza A viruses. The HA of the bat virus was estimated to have diverged at roughly the same time as the known subtypes of HA and was designated as H17. The neuraminidase (NA) gene is highly divergent from all known influenza NAs, and the internal genes from the bat virus diverged from those of known influenza A viruses before the estimated divergence of the known influenza A internal gene lineages. Attempts to propagate this virus in cell cultures and chicken embryos were unsuccessful, suggesting distinct requirements compared with known influenza viruses. Despite its divergence from known influenza A viruses, the bat virus is compatible for genetic exchange with human influenza viruses in human cells, suggesting the potential capability for reassortment and contributions to new pandemic or panzootic influenza A viruses.

Novel paramyxovirus associated with severe acute febrile disease, South Sudan and Uganda, 2012.

In 2012, a female wildlife biologist experienced fever, malaise, headache, generalized myalgia and arthralgia, neck stiffness, and a sore throat shortly after returning to the United States from a 6-week field expedition to South Sudan and Uganda. She was hospitalized, after which a maculopapular rash developed and became confluent. When the patient was discharged from the hospital on day 14, arthralgia and myalgia had improved, oropharynx ulcerations had healed, the rash had resolved without desquamation, and blood counts and hepatic enzyme levels were returning to reference levels. After several known suspect pathogens were ruled out as the cause of her illness, deep sequencing and metagenomics analysis revealed a novel paramyxovirus related to rubula-like viruses isolated from fruit bats.

Prevalence of sequence types among clinical and environmental isolates of Legionella pneumophila serogroup 1 in the United States from 1982 to 2012.

Since the establishment of sequence-based typing as the gold standard for DNA-based typing of Legionella pneumophila, the Legionella laboratory at the Centers for Disease Control and Prevention (CDC) has conducted routine sequence-based typing (SBT) analysis of all incoming L. pneumophila serogroup 1 (Lp1) isolates to identify potential links between cases and to better understand genetic diversity and clonal expansion among L. pneumophila bacteria. Retrospective genotyping of Lp1 isolates from sporadic cases and Legionnaires' disease (LD) outbreaks deposited into the CDC reference collection since 1982 has been completed. For this study, we compared the distribution of sequence types (STs) among Lp1 isolates implicated in 26 outbreaks in the United States, 571 clinical isolates from sporadic cases of LD in the United States, and 149 environmental isolates with no known association with LD. The Lp1 isolates under study had been deposited into our collection between 1982 and 2012. We identified 17 outbreak-associated STs, 153 sporadic STs, and 49 environmental STs. We observed that Lp1 STs from outbreaks and sporadic cases are more similar to each other than either group is to environmental STs. The most frequent ST for both sporadic and environmental isolates was ST1, accounting for 25% and 49% of the total number of isolates, respectively. The STs shared by both outbreak-associated and sporadic Lp1 included ST1, ST35, ST36, ST37, and ST222. The STs most commonly found in sporadic and outbreak-associated Lp1 populations may have an increased ability to cause disease and thus may require special attention when detected.

Whole-genome analysis of Exserohilum rostratum from an outbreak of fungal meningitis and other infections.

Exserohilum rostratum was the cause of most cases of fungal meningitis and other infections associated with the injection of contaminated methylprednisolone acetate produced by the New England Compounding Center (NECC). Until this outbreak, very few human cases of Exserohilum infection had been reported, and very little was known about this dematiaceous fungus, which usually infects plants. Here, we report using whole-genome sequencing (WGS) for the detection of single nucleotide polymorphisms (SNPs) and phylogenetic analysis to investigate the molecular origin of the outbreak using 22 isolates of E. rostratum retrieved from 19 case patients with meningitis or epidural/spinal abscesses, 6 isolates from contaminated NECC vials, and 7 isolates unrelated to the outbreak. Our analysis indicates that all 28 isolates associated with the outbreak had nearly identical genomes of 33.8 Mb. A total of 8 SNPs were detected among the outbreak genomes, with no more than 2 SNPs separating any 2 of the 28 genomes. The outbreak genomes were separated from the next most closely related control strain by ∼136,000 SNPs. We also observed significant genomic variability among strains unrelated to the outbreak, which may suggest the possibility of cryptic speciation in E. rostratum.

Novel epidemic clones of Listeria monocytogenes, United States, 2011.

We identified a novel serotype 1/2a outbreak strain and 2 novel epidemic clones of Listeria monocytogenes while investigating a foodborne outbreak of listeriosis associated with consumption of cantaloupe during 2011 in the United States. Comparative analyses of strains worldwide are essential to identification of novel outbreak strains and epidemic clones.

Phylogenetic and ecologic perspectives of a monkeypox outbreak, southern Sudan, 2005.

Identification of human monkeypox cases during 2005 in southern Sudan (now South Sudan) raised several questions about the natural history of monkeypox virus (MPXV) in Africa. The outbreak area, characterized by seasonally dry riverine grasslands, is not identified as environmentally suitable for MPXV transmission. We examined possible origins of this outbreak by performing phylogenetic analysis of genome sequences of MPXV isolates from the outbreak in Sudan and from differing localities. We also compared the environmental suitability of study localities for monkeypox transmission. Phylogenetically, the viruses isolated from Sudan outbreak specimens belong to a clade identified in the Congo Basin. This finding, added to the political instability of the area during the time of the outbreak, supports the hypothesis of importation by infected animals or humans entering Sudan from the Congo Basin, and person-to-person transmission of virus, rather than transmission of indigenous virus from infected animals to humans.

Progressive vaccinia: case description and laboratory-guided therapy with vaccinia immune globulin, ST-246, and CMX001.

Progressive vaccinia (PV) is a rare but potentially lethal complication that develops in smallpox vaccine recipients with severely impaired cellular immunity. We describe a patient with PV who required treatment with vaccinia immune globulin and who received 2 investigational agents, ST-246 and CMX001. We describe the various molecular, pharmacokinetic, and immunologic studies that provided guidance to escalate and then successfully discontinue therapy. Despite development of resistance to ST-246 during treatment, the patient had resolution of PV. This case demonstrates the need for continued development of novel anti-orthopoxvirus pharmaceuticals and the importance of both intensive and timely clinical and laboratory support in management of PV.

Genomic basis of a polyagglutinating isolate of Neisseria meningitidis.

Containment strategies for outbreaks of invasive Neisseria meningitidis disease are informed by serogroup assays that characterize the polysaccharide capsule. We sought to uncover the genomic basis of conflicting serogroup assay results for an isolate (M16917) from a patient with acute meningococcal disease. To this end, we characterized the complete genome sequence of the M16917 isolate and performed a variety of comparative sequence analyses against N. meningitidis reference genome sequences of known serogroups. Multilocus sequence typing and whole-genome sequence comparison revealed that M16917 is a member of the ST-11 sequence group, which is most often associated with serogroup C. However, sequence similarity comparisons and phylogenetic analysis showed that the serogroup diagnostic capsule polymerase gene (synD) of M16917 belongs to serogroup B. These results suggest that a capsule-switching event occurred based on homologous recombination at or around the capsule locus of M16917. Detailed analysis of this locus uncovered the locations of recombination breakpoints in the M16917 genome sequence, which led to the introduction of an ∼2-kb serogroup B sequence cassette into the serogroup C genomic background. Since there is no currently available vaccine for serogroup B strains of N. meningitidis, this kind capsule-switching event could have public health relevance as a vaccine escape mutant.

Analysis of a unique Clostridium botulinum strain from the Southern hemisphere producing a novel type E botulinum neurotoxin subtype.

BACKGROUND: Clostridium botulinum strains that produce botulinum neurotoxin type E (BoNT/E) are most commonly isolated from botulism cases, marine environments, and animals in regions of high latitude in the Northern hemisphere. A strain of C. botulinum type E (CDC66177) was isolated from soil in Chubut, Argentina. Previous studies showed that the amino acid sequences of BoNT/E produced by various strains differ by < 6% and that the type E neurotoxin gene cluster inserts into the rarA operon. RESULTS: Genetic and mass spectral analysis demonstrated that the BoNT/E produced by CDC66177 is a novel toxin subtype (E9). Toxin gene sequencing indicated that BoNT/E9 differed by nearly 11% at the amino acid level compared to BoNT/E1. Mass spectrometric analysis of BoNT/E9 revealed that its endopeptidase substrate cleavage site was identical to other BoNT/E subtypes. Further analysis of this strain demonstrated that its 16S rRNA sequence clustered with other Group II C. botulinum (producing BoNT types B, E, and F) strains. Genomic DNA isolated from strain CDC66177 hybridized with fewer probes using a Group II C. botulinum subtyping microarray compared to other type E strains examined. Whole genome shotgun sequencing of strain CDC66177 revealed that while the toxin gene cluster inserted into the rarA operon similar to other type E strains, its overall genome content shared greater similarity with a Group II C. botulinum type B strain (17B). CONCLUSIONS: These results expand our understanding of the global distribution of C. botulinum type E strains and suggest that the type E toxin gene cluster may be able to insert into C. botulinum strains with a more diverse genetic background than previously recognized.

Genome sequences for five strains of the emerging pathogen Haemophilus haemolyticus.

We report the first whole-genome sequences for five strains, two carried and three pathogenic, of the emerging pathogen Haemophilus haemolyticus. Preliminary analyses indicate that these genome sequences encode markers that distinguish H. haemolyticus from its closest Haemophilus relatives and provide clues to the identity of its virulence factors.