Precision Tracing of Household Dengue Spread Using Inter- and Intra-Host Viral Variation Data, Kamphaeng Phet, Thailand.

Dengue control approaches are best informed by granular spatial epidemiology of these viruses, yet reconstruction of inter- and intra-household transmissions is limited when analyzing case count, serologic, or genomic consensus sequence data. To determine viral spread on a finer spatial scale, we extended phylogenomic discrete trait analyses to reconstructions of house-to-house transmissions within a prospective cluster study in Kamphaeng Phet, Thailand. For additional resolution and transmission confirmation, we mapped dengue intra-host single nucleotide variants on the taxa of these time-scaled phylogenies. This approach confirmed 19 household transmissions and revealed that dengue disperses an average of 70 m per day between households in these communities. We describe an evolutionary biology framework for the resolution of dengue transmissions that cannot be differentiated based on epidemiologic and consensus genome data alone. This framework can be used as a public health tool to inform control approaches and enable precise tracing of dengue transmissions.

Next Generation Sequencing and Bioinformatics Methodologies for Infectious Disease Research and Public Health: Approaches, Applications, and Considerations for Development of Laboratory Capacity.

Next generation sequencing (NGS) combined with bioinformatics has successfully been used in a vast array of analyses for infectious disease research of public health relevance. For instance, NGS and bioinformatics approaches have been used to identify outbreak origins, track transmissions, investigate epidemic dynamics, determine etiological agents of a disease, and discover novel human pathogens. However, implementation of high-quality NGS and bioinformatics in research and public health laboratories can be challenging. These challenges mainly include the choice of the sequencing platform and the sequencing approach, the choice of bioinformatics methodologies, access to the appropriate computation and information technology infrastructure, and recruiting and retaining personnel with the specialized skills and experience in this field. In this review, we summarize the most common NGS and bioinformatics workflows in the context of infectious disease genomic surveillance and pathogen discovery, and highlight the main challenges and considerations for setting up an NGS and bioinformatics-focused infectious disease research public health laboratory. We describe the most commonly used sequencing platforms and review their strengths and weaknesses. We review sequencing approaches that have been used for various pathogens and study questions, as well as the most common difficulties associated with these approaches that should be considered when implementing in a public health or research setting. In addition, we provide a review of some common bioinformatics tools and procedures used for pathogen discovery and genome assembly, along with the most common challenges and solutions. Finally, we summarize the bioinformatics of advanced viral, bacterial, and parasite pathogen characterization, including types of study questions that can be answered when utilizing NGS and bioinformatics.

Dengue diversity across spatial and temporal scales: Local structure and the effect of host population size.

A fundamental mystery for dengue and other infectious pathogens is how observed patterns of cases relate to actual chains of individual transmission events. These pathways are intimately tied to the mechanisms by which strains interact and compete across spatial scales. Phylogeographic methods have been used to characterize pathogen dispersal at global and regional scales but have yielded few insights into the local spatiotemporal structure of endemic transmission. Using geolocated genotype (800 cases) and serotype (17,291 cases) data, we show that in Bangkok, Thailand, 60% of dengue cases living <200 meters apart come from the same transmission chain, as opposed to 3% of cases separated by 1 to 5 kilometers. At distances <200 meters from a case (encompassing an average of 1300 people in Bangkok), the effective number of chains is 1.7. This number rises by a factor of 7 for each 10-fold increase in the population of the "enclosed" region. This trend is observed regardless of whether population density or area increases, though increases in density over 7000 people per square kilometer do not lead to additional chains. Within Thailand these chains quickly mix, and by the next dengue season viral lineages are no longer highly spatially structured within the country. In contrast, viral flow to neighboring countries is limited. These findings are consistent with local, density-dependent transmission and implicate densely populated communities as key sources of viral diversity, with home location the focal point of transmission. These findings have important implications for targeted vector control and active surveillance.

Frequency of influenza H3N2 intra-subtype reassortment: attributes and implications of reassortant spread.

BACKGROUND: Increasing evidence suggests that influenza reassortment not only contributes to the emergence of new human pandemics but also plays an important role in seasonal influenza epidemics, disease severity, evolution, and vaccine efficacy. We studied this process within 2091 H3N2 full genomes utilizing a combination of the latest reassortment detection tools and more conventional phylogenetic analyses. RESULTS: We found that the amount of H3N2 intra-subtype reassortment depended on the number of sampled genomes, occurred with a steady frequency of 3.35%, and was not affected by the geographical origins, evolutionary patterns, or previous reassortment history of the virus. We identified both single reassortant genomes and reassortant clades, each clade representing one reassortment event followed by successful spread of the reassorted variant in the human population. It was this spread that was mainly responsible for the observed high presence of H3N2 intra-subtype reassortant genomes. The successfully spread variants were generally sampled within one year of their formation, highlighting the risk of their rapid spread but also presenting an opportunity for their rapid detection. Simultaneous spread of several different reassortant lineages was observed, and despite their limited average lifetime, second and third generation reassortment was detected, as well as reassortment between viruses belonging to different vaccine-associated clades, likely displaying differing antigenic properties. Some of the spreading reassortants remained confined to certain geographical regions, while others, sharing common properties in amino acid positions of the HA, NA, and PB2 segments, were found throughout the world. CONCLUSIONS: Detailed surveillance of seasonal influenza reassortment patterns and variant properties may provide unique information needed for prediction of spread and construction of future influenza vaccines.

Human Sentinel Surveillance of Influenza and Other Respiratory Viral Pathogens in Border Areas of Western Cambodia.

Little is known about circulation of influenza and other respiratory viruses in remote populations along the Thai-Cambodia border in western Cambodia. We screened 586 outpatients (median age 5, range 1-77) presenting with influenza-like-illness (ILI) at 4 sentinel sites in western Cambodia between May 2010 and December 2012. Real-time reverse transcriptase (rRT) PCR for influenza was performed on combined nasal and throat specimens followed by viral culture, antigenic analysis, antiviral susceptibility testing and full genome sequencing for phylogenetic analysis. ILI-specimens negative for influenza were cultured, followed by rRT-PCR for enterovirus and rhinovirus (EV/RV) and EV71. Influenza was found in 168 cases (29%) and occurred almost exclusively in the rainy season from June to November. Isolated influenza strains had close antigenic and phylogenetic relationships, matching vaccine and circulating strains found elsewhere in Cambodia. Influenza vaccination coverage was low (<20%). Western Cambodian H1N1(2009) isolate genomes were more closely related to 10 earlier Cambodia isolates (94.4% genome conservation) than to 13 Thai isolates (75.9% genome conservation), despite sharing the majority of the amino acid changes with the Thai references. Most genes showed signatures of purifying selection. Viral culture detected only adenovirus (5.7%) and parainfluenza virus (3.8%), while non-polio enteroviruses (10.3%) were detected among 164 culture-negative samples including coxsackievirus A4, A6, A8, A9, A12, B3, B4 and echovirus E6 and E9 using nested RT-PCR methods. A single specimen of EV71 was found. Despite proximity to Thailand, influenza epidemiology of these western Cambodian isolates followed patterns observed elsewhere in Cambodia, continuing to support current vaccine and treatment recommendations from the Cambodian National Influenza Center. Amino acid mutations at non-epitope sites, particularly hemagglutinin genes, require further investigation in light of an increasingly important role of permissive mutations in influenza virus evolution. Further research about the burden of adenovirus and non-polio enteroviruses as etiologic agents in acute respiratory infections in Cambodia is also needed.

Pathosphere.org: pathogen detection and characterization through a web-based, open source informatics platform.

BACKGROUND: The detection of pathogens in complex sample backgrounds has been revolutionized by wide access to next-generation sequencing (NGS) platforms. However, analytical methods to support NGS platforms are not as uniformly available. Pathosphere (found at Pathosphere.org) is a cloud - based open - sourced community tool that allows for communication, collaboration and sharing of NGS analytical tools and data amongst scientists working in academia, industry and government. The architecture allows for users to upload data and run available bioinformatics pipelines without the need for onsite processing hardware or technical support. RESULTS: The pathogen detection capabilities hosted on Pathosphere were tested by analyzing pathogen-containing samples sequenced by NGS with both spiked human samples as well as human and zoonotic host backgrounds. Pathosphere analytical pipelines developed by Edgewood Chemical Biological Center (ECBC) identified spiked pathogens within a common sample analyzed by 454, Ion Torrent, and Illumina sequencing platforms. ECBC pipelines also correctly identified pathogens in human samples containing arenavirus in addition to animal samples containing flavivirus and coronavirus. These analytical methods were limited in the detection of sequences with limited homology to previous annotations within NCBI databases, such as parvovirus. Utilizing the pipeline-hosting adaptability of Pathosphere, the analytical suite was supplemented by analytical pipelines designed by the United States Army Medical Research Insititute of Infectious Diseases and Walter Reed Army Institute of Research (USAMRIID-WRAIR). These pipelines were implemented and detected parvovirus sequence in the sample that the ECBC iterative analysis previously failed to identify. CONCLUSIONS: By accurately detecting pathogens in a variety of samples, this work demonstrates the utility of Pathosphere and provides a platform for utilizing, modifying and creating pipelines for a variety of NGS technologies developed to detect pathogens in complex sample backgrounds. These results serve as an exhibition for the existing pipelines and web-based interface of Pathosphere as well as the plug-in adaptability that allows for integration of newer NGS analytical software as it becomes available.

Recombination Does Not Hinder Formation or Detection of Ecological Species of Synechococcus Inhabiting a Hot Spring Cyanobacterial Mat.

Recent studies of bacterial speciation have claimed to support the biological species concept-that reduced recombination is required for bacterial populations to diverge into species. This conclusion has been reached from the discovery that ecologically distinct clades show lower rates of recombination than that which occurs among closest relatives. However, these previous studies did not attempt to determine whether the more-rapidly recombining close relatives within the clades studied may also have diversified ecologically, without benefit of sexual isolation. Here we have measured the impact of recombination on ecological diversification within and between two ecologically distinct clades (A and B') of Synechococcus in a hot spring microbial mat in Yellowstone National Park, using a cultivation-free, multi-locus approach. Bacterial artificial chromosome (BAC) libraries were constructed from mat samples collected at 60°C and 65°C. Analysis of multiple linked loci near Synechococcus 16S rRNA genes showed little evidence of recombination between the A and B' lineages, but a record of recombination was apparent within each lineage. Recombination and mutation rates within each lineage were of similar magnitude, but recombination had a somewhat greater impact on sequence diversity than mutation, as also seen in many other bacteria and archaea. Despite recombination within the A and B' lineages, there was evidence of ecological diversification within each lineage. The algorithm Ecotype Simulation identified sequence clusters consistent with ecologically distinct populations (ecotypes), and several hypothesized ecotypes were distinct in their habitat associations and in their adaptations to different microenvironments. We conclude that sexual isolation is more likely to follow ecological divergence than to precede it. Thus, an ecology-based model of speciation appears more appropriate than the biological species concept for bacterial and archaeal diversification.

Characterization based on the 56-Kda type-specific antigen gene of Orientia tsutsugamushi genotypes isolated from Leptotrombidium mites and the rodent host post-infection.

Abstract. Characterization of the 56-kDa type-specific antigen (TSA) genes of Orientia tsutsugamushi (OT) from three naturally infected, laboratory-reared mite colonies comprising three species (Leptotrombidium deliense [Ld], Leptotrombidium imphalum [Li], and Leptotrombidium chiangraiensis [Lc]) has revealed the presence of single and coexisting OT genotypes found in individual chiggers. The Karp genotype was found in all of the chiggers examined, whereas Gilliam and UT302 genotypes were only observed in combination with the Karp genotype. From analysis of these OT genotypes after transmission from chiggers to mice it was determined that with the Lc and Li mites, the OT genotype composition in the rodent spleens post-infection had not changed and therefore resembled that observed in the feeding chiggers. However, only the Karp genotype was found in rodents after feeding by Ld chiggers carrying Karp and Gilliam genotypes. The current findings reveal a complex association among the host, pathogen, and vector.

Hepatitis E virus is prevalent in the pig population of Lao People's Democratic Republic and evidence exists for homogeneity with Chinese Genotype 4 human isolates.

The objective of this study was to determine the prevalence and genotypic range of Hepatitis E virus (HEV) in the pig population of northern Lao People's Democratic Republic (PDR). We collected 181 faecal samples from indigenous-breed pigs ≤ 6 months of age and the faeces was stored in RNA stabilisation buffer due to cold-chain and transport limitations. Twenty-one (11.6%) pigs had detectable HEV RNA and 43.5% of village pig herds were infected. Based on a 240 base pair-nucleotide sequence flanking the junction of open reading frames 1, 2 and 3 (ORF1, ORF2 and ORF3) the isolates were phylogenetically classified within genotype 4. Phylogenetic analyses revealed distinct genetic groupings of the Lao HEV isolates and two groups clustered with human and pig HEV isolates from China. This was the first study to demonstrate genotype 4 HEV in Lao PDR and indicates pigs are a potential reservoir for human HEV infection.

Influence of molecular resolution on sequence-based discovery of ecological diversity among Synechococcus populations in an alkaline siliceous hot spring microbial mat.

Previous research has shown that sequences of 16S rRNA genes and 16S-23S rRNA internal transcribed spacer regions may not have enough genetic resolution to define all ecologically distinct Synechococcus populations (ecotypes) inhabiting alkaline, siliceous hot spring microbial mats. To achieve higher molecular resolution, we studied sequence variation in three protein-encoding loci sampled by PCR from 60°C and 65°C sites in the Mushroom Spring mat (Yellowstone National Park, WY). Sequences were analyzed using the ecotype simulation (ES) and AdaptML algorithms to identify putative ecotypes. Between 4 and 14 times more putative ecotypes were predicted from variation in protein-encoding locus sequences than from variation in 16S rRNA and 16S-23S rRNA internal transcribed spacer sequences. The number of putative ecotypes predicted depended on the number of sequences sampled and the molecular resolution of the locus. Chao estimates of diversity indicated that few rare ecotypes were missed. Many ecotypes hypothesized by sequence analyses were different in their habitat specificities, suggesting different adaptations to temperature or other parameters that vary along the flow channel.

Population level functional diversity in a microbial community revealed by comparative genomic and metagenomic analyses.

In microbial mat communities of Yellowstone hot springs, ribosomal RNA (rRNA) sequence diversity patterns indicate the presence of closely related bacterial populations along environmental gradients of temperature and light. To identify the functional bases for adaptation, we sequenced the genomes of two cyanobacterial (Synechococcus OS-A and OS-B') isolates representing ecologically distinct populations that dominate at different temperatures and are major primary producers in the mat. There was a marked lack of conserved large-scale gene order between the two Synechococcus genomes, indicative of extensive genomic rearrangements. Comparative genomic analyses showed that the isolates shared a large fraction of their gene content at high identity, yet, differences in phosphate and nitrogen utilization pathways indicated that they have adapted differentially to nutrient fluxes, possibly by the acquisition of genes by lateral gene transfer or their loss in certain populations. Comparisons of the Synechococcus genomes to metagenomic sequences derived from mats where these Synechococcus stains were originally isolated, revealed new facets of microbial diversity. First, Synechococcus populations at the lower temperature regions of the mat showed greater sequence diversity than those at high temperatures, consistent with a greater number of ecologically distinct populations at the lower temperature. Second, we found evidence of a specialized population that is apparently very closely related to Synechococcus OS-B', but contains genes that function in the uptake of reduced ferrous iron. In situ expression studies demonstrated that these genes are differentially expressed over the diel cycle, with highest expression when the mats are anoxic and iron may be in the reduced state. Genomic information from these mat-specific isolates and metagenomic information can be coupled to detect naturally occurring populations that are associated with different functionalities, not always represented by isolates, but which may nevertheless be important for niche partitioning and the establishment of microbial community structure.

In situ analysis of nitrogen fixation and metabolic switching in unicellular thermophilic cyanobacteria inhabiting hot spring microbial mats.

Genome sequences of two Synechococcus ecotypes inhabiting the Octopus Spring microbial mat in Yellowstone National Park revealed the presence of all genes required for nitrogenase biosynthesis. We demonstrate that nif genes of the Synechococcus ecotypes are expressed in situ in a region of the mat that varies in temperature from 53.5 degrees C to 63.4 degrees C (average 60 degrees C); transcripts are only detected at the end of the day when the mat becomes anoxic. Nitrogenase activity in mat samples was also detected in the evening. Hitherto, N2 fixation in hot spring mats was attributed either to filamentous cyanobacteria (not present at >50 degrees C in these mats) or to heterotrophic bacteria. To explore how energy-generating processes of the Synechococcus ecotypes track natural light and O2 conditions, we evaluated accumulation of transcripts encoding proteins involved in photosynthesis, respiration, and fermentation. Transcripts from photosynthesis (cpcF, cpcE, psaB, and psbB) and respiration (coxA and cydA) genes declined in the evening. In contrast, transcripts encoding enzymes that may participate in fermentation fell into two categories; some (ldh, pdhB, ald, and ackA) decreased in the evening, whereas others (pflB, pflA, adhE, and acs) increased at the end of the day and remained high into the night. Energy required for N2 fixation during the night may be derived from fermentation pathways that become prominent as the mat becomes anoxic. In a broader context, our data suggest that there are critical regulatory switches in situ that are linked to the diel cycle and that these switches alter many metabolic processes within the microbial mat.