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.

Genomic Characterization and Copy Number Variation of Bacillus anthracis Plasmids pXO1 and pXO2 in a Historical Collection of 412 Strains.

Bacillus anthracis plasmids pXO1 and pXO2 carry the main virulence factors responsible for anthrax. However, the extent of copy number variation within the species and how the plasmids are related to pXO1/pXO2-like plasmids in other species of the Bacillus cereus sensu lato group remain unclear. To gain new insights into these issues, we sequenced 412 B. anthracis strains representing the total phylogenetic and ecological diversity of the species. Our results revealed that B. anthracis genomes carried, on average, 3.86 and 2.29 copies of pXO1 and pXO2, respectively, and also revealed a positive linear correlation between the copy numbers of pXO1 and pXO2. No correlation between the plasmid copy number and the phylogenetic relatedness of the strains was observed. However, genomes of strains isolated from animal tissues generally maintained a higher plasmid copy number than genomes of strains from environmental sources (P < 0.05 [Welch two-sample t test]). Comparisons against B. cereus genomes carrying complete or partial pXO1-like and pXO2-like plasmids showed that the plasmid-based phylogeny recapitulated that of the main chromosome, indicating limited plasmid horizontal transfer between or within these species. Comparisons of gene content revealed a closed pXO1 and pXO2 pangenome; e.g., plasmids encode <8 unique genes, on average, and a single large fragment deletion of pXO1 in one B. anthracis strain (2000031682) was detected. Collectively, our results provide a more complete view of the genomic diversity of B. anthracis plasmids, their copy number variation, and the virulence potential of other Bacillus species carrying pXO1/pXO2-like plasmids. IMPORTANCE Bacillus anthracis microorganisms are of historical and epidemiological importance and are among the most homogenous bacterial groups known, even though the B. anthracis genome is rich in mobile elements. Mobile elements can trigger the diversification of lineages; therefore, characterizing the extent of genomic variation in a large collection of strains is critical for a complete understanding of the diversity and evolution of the species. Here, we sequenced a large collection of B. anthracis strains (>400) that were recovered from human, animal, and environmental sources around the world. Our results confirmed the remarkable stability of gene content and synteny of the anthrax plasmids and revealed no signal of plasmid exchange between B. anthracis and pathogenic B. cereus isolates but rather predominantly vertical descent. These findings advance our understanding of the biology and pathogenomic evolution of B. anthracis and its plasmids.

Draft Genome Sequence of Bacillus cereus LA2007, a Human-Pathogenic Isolate Harboring Anthrax-Like Plasmids.

We present the genome sequence of Bacillus cereus LA2007, a strain isolated in 2007 from a fatal pneumonia case in Louisiana. Sequence-based genome analysis revealed that LA2007 carries a plasmid highly similar to Bacillus anthracis pXO1, including the genes responsible for the production and regulation of anthrax toxin.

Genome Sequences of Two Brucella suis Strains Isolated from the Same Patient, 8 Years Apart.

Brucella suis is a Gram-negative, facultative intracellular pathogen that has pigs as its preferred host, but it can also infect humans. Here, we report the draft genome sequences of two B. suis strains that were isolated from the same patient, 8 years apart.

Genome Sequences of Three Brucella canis Strains Isolated from Humans and a Dog.

Brucella canis is a facultative intracellular pathogen that preferentially infects members of the Canidae family. Here, we report the genome sequencing of two Brucella canis strains isolated from humans and one isolated from a dog host.

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.

Serro 2 Virus Highlights the Fundamental Genomic and Biological Features of a Natural Vaccinia Virus Infecting Humans.

Vaccinia virus (VACV) has been implicated in infections of dairy cattle and humans, and outbreaks have substantially impacted local economies and public health in Brazil. During a 2005 outbreak, a VACV strain designated Serro 2 virus (S2V) was collected from a 30-year old male milker. Our aim was to phenotypically and genetically characterize this VACV Brazilian isolate. S2V produced small round plaques without associated comets when grown in BSC40 cells. Furthermore, S2V was less virulent than the prototype strain VACV-Western Reserve (WR) in a murine model of intradermal infection, producing a tiny lesion with virtually no surrounding inflammation. The genome of S2V was sequenced by primer walking. The coding region spans 184,572 bp and contains 211 predicted genes. Mutations in envelope genes specifically associated with small plaque phenotypes were not found in S2V; however, other alterations in amino acid sequences within these genes were identified. In addition, some immunomodulatory genes were truncated in S2V. Phylogenetic analysis using immune regulatory-related genes, besides the hemagglutinin gene, segregated the Brazilian viruses into two clusters, grouping the S2V into Brazilian VACV group 1. S2V is the first naturally-circulating human-associated VACV, with a low passage history, to be extensively genetically and phenotypically characterized.

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.

Complete Genome Sequences for Three Chromosomes of the Burkholderia stabilis Type Strain (ATCC BAA-67).

We report here the complete annotated genome sequence of the Burkholderia stabilis type strain ATCC BAA-67. There were three circular chromosomes with a combined size of 8,527,947 bp and G+C composition of 66.4%. These characteristics closely resemble the genomes of other sequenced members of the Burkholderia cepacia complex.

Dynamics of genome change among Legionella species.

Legionella species inhabit freshwater and soil ecosystems where they parasitize protozoa. L. pneumonphila (LP) serogroup-1 (Lp1) is the major cause of Legionnaires' Disease (LD), a life-threatening pulmonary infection that can spread systemically. The increased global frequency of LD caused by Lp and non-Lp species underscores the need to expand our knowledge of evolutionary forces underlying disease pathogenesis. Whole genome analyses of 43 strains, including all known Lp serogroups 1-17 and 17 emergent LD-causing Legionella species (of which 33 were sequenced in this study) in addition to 10 publicly available genomes, resolved the strains into four phylogenetic clades along host virulence demarcations. Clade-specific genes were distinct for genetic exchange and signal-transduction, indicating adaptation to specific cellular and/or environmental niches. CRISPR spacer comparisons hinted at larger pools of accessory DNA sequences in Lp than predicted by the pan-genome analyses. While recombination within Lp was frequent and has been reported previously, population structure analysis identified surprisingly few DNA admixture events between species. In summary, diverse Legionella LD-causing species share a conserved core-genome, are genetically isolated from each other, and selectively acquire genes with potential for enhanced virulence.

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.

Three Genome Sequences of Legionella pneumophila subsp. pascullei Associated with Colonization of a Health Care Facility.

Here, we report the complete genome sequences of three Legionella pneumophila subsp. pascullei strains (including both serogroup 1 and 5 strains) that were found in the same health care facility in 1982 and 2012.

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.

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.

Identification and Characterization of Microsatellite Markers Derived from the Whole Genome Analysis of Taenia solium.

BACKGROUND: Infections with Taenia solium are the most common cause of adult acquired seizures worldwide, and are the leading cause of epilepsy in developing countries. A better understanding of the genetic diversity of T. solium will improve parasite diagnostics and transmission pathways in endemic areas thereby facilitating the design of future control measures and interventions. Microsatellite markers are useful genome features, which enable strain typing and identification in complex pathogen genomes. Here we describe microsatellite identification and characterization in T. solium, providing information that will assist in global efforts to control this important pathogen. METHODS: For genome sequencing, T. solium cysts and proglottids were collected from Huancayo and Puno in Peru, respectively. Using next generation sequencing (NGS) and de novo assembly, we assembled two draft genomes and one hybrid genome. Microsatellite sequences were identified and 36 of them were selected for further analysis. Twenty T. solium isolates were collected from Tumbes in the northern region, and twenty from Puno in the southern region of Peru. The size-polymorphism of the selected microsatellites was determined with multi-capillary electrophoresis. We analyzed the association between microsatellite polymorphism and the geographic origin of the samples. RESULTS: The predicted size of the hybrid (proglottid genome combined with cyst genome) T. solium genome was 111 MB with a GC content of 42.54%. A total of 7,979 contigs (>1,000 nt) were obtained. We identified 9,129 microsatellites in the Puno-proglottid genome and 9,936 in the Huancayo-cyst genome, with 5 or more repeats, ranging from mono- to hexa-nucleotide. Seven microsatellites were polymorphic and 29 were monomorphic within the analyzed isolates. T. solium tapeworms were classified into two genetic groups that correlated with the North/South geographic origin of the parasites. CONCLUSIONS/SIGNIFICANCE: The availability of draft genomes for T. solium represents a significant step towards the understanding the biology of the parasite. We report here a set of T. solium polymorphic microsatellite markers that appear promising for genetic epidemiology studies.

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.

Genetic similarities between Cyclospora cayetanensis and cecum-infecting avian Eimeria spp. in apicoplast and mitochondrial genomes.

BACKGROUND: Cyclospora cayetanensis is an important cause for diarrhea in children in developing countries and foodborne outbreaks of cyclosporiasis in industrialized nations. To improve understanding of the basic biology of Cyclospora spp. and development of molecular diagnostic tools and therapeutics, we sequenced the complete apicoplast and mitochondrial genomes of C. cayetanensis. METHODS: The genome of one Chinese C. cayetanensis isolate was sequenced using Roche 454 and Illumina technologies. The assembled genomes of the apicoplast and mitochondrion were retrieved, annotated, and compared with reference genomes for other apicomplexans to infer genome organizations and phylogenetic relationships. Sequence variations in the mitochondrial genome were identified by comparison of two C. cayetanensis nucleotide sequences from this study and a recent publication. RESULTS: The apicoplast and mitochondrial genomes of C. cayetanensis are 34,155 and 6,229 bp in size and code for 65 and 5 genes, respectively. Comparative genomic analysis showed high similarities between C. cayetanensis and Eimeria tenella in both genomes; they have 85.6% and 90.4% nucleotide sequence similarities, respectively, and complete synteny in gene organization. Phylogenetic analysis of the genomic sequences confirmed the genetic similarities between cecum-infecting avian Eimeria spp. and C. cayetanensis. Like in other coccidia, both genomes of C. cayetanensis are transcribed bi-directionally. The apicoplast genome is circular, codes for the complete machinery for protein biosynthesis, and contains two inverted repeats that differ slightly in LSU rRNA gene sequences. In contrast, the mitochondrial genome has a linear concatemer or circular mapping topology. Eight single-nucleotide and one 7-bp multiple-nucleotide variants were detected between the mitochondrial genomes of C. cayetanensis from this and recent studies. CONCLUSIONS: The apicoplast and mitochondrial genomes of C. cayetanensis are highly similar to those of cecum-infecting avian Eimeria spp. in both genome organization and sequences. The availability of sequence data beyond rRNA and heat shock protein genes could facilitate studies of C. cayetanensis biology and development of genotyping tools for investigations of cyclosporiasis outbreaks.

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.

Insights into the environmental reservoir of pathogenic Vibrio parahaemolyticus using comparative genomics.

Vibrio parahaemolyticus is an aquatic halophilic bacterium that occupies estuarine and coastal marine environments, and is a leading cause of seafood-borne food poisoning cases. To investigate the environmental reservoir and potential gene flow that occurs among V. parahaemolyticus isolates, the virulence-associated gene content and genome diversity of a collection of 133 V. parahaemolyticus isolates were analyzed. Phylogenetic analysis of housekeeping genes, and pulsed-field gel electrophoresis, demonstrated that there is genetic similarity among V. parahaemolyticus clinical and environmental isolates. Whole-genome sequencing and comparative analysis of six representative V. parahaemolyticus isolates was used to identify genes that are unique to the clinical and environmental isolates examined. Comparative genomics demonstrated an O3:K6 environmental isolate, AF91, which was cultured from sediment collected in Florida in 2006, has significant genomic similarity to the post-1995 O3:K6 isolates. However, AF91 lacks the majority of the virulence-associated genes and genomic islands associated with these highly virulent post-1995 O3:K6 genomes. These findings demonstrate that although they do not contain most of the known virulence-associated regions, some V. parahaemolyticus environmental isolates exhibit significant genetic similarity to clinical isolates. This highlights the dynamic nature of the V. parahaemolyticus genome allowing them to transition between aquatic and host-pathogen states.

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.