Search

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.

Dengue virus (DENV) causes a profound burden of morbidity and mortality, and its global burden is rising due to the co-circulation of four divergent DENV serotypes in the ecological context of globalization, travel, climate change, urbanization, and expansion of the geographic range of the Ae.aegypti and Ae.albopictus vectors. Understanding DENV evolution offers valuable opportunities to enhance surveillance and response to DENV epidemics via advances in RNA virus sequencing, bioinformatics, phylogenetic and other computational biology methods. Here we provide a scoping overview of the evolution and molecular epidemiology of DENV and the range of ways that evolutionary analyses can be applied as a public health tool against this arboviral pathogen.

Next-generation sequencing technologies, exponential increases in the availability of virus genomic data, and ongoing advances in phylogenomic methods have made genomic epidemiology an increasingly powerful tool for public health response to a range of mosquito-borne virus outbreaks. In this review, we offer a brief primer on the scope and methods of phylogenomic analyses that can answer key epidemiological questions during mosquito-borne virus public health emergencies. We then focus on case examples of outbreaks, including those caused by dengue, Zika, yellow fever, West Nile, and chikungunya viruses, to demonstrate the utility of genomic epidemiology to support the prevention and control of mosquito-borne virus threats. We extend these case studies with operational perspectives on how to best incorporate genomic epidemiology into structured surveillance and response programs for mosquito-borne virus control. Many tools for genomic epidemiology already exist, but so do technical and nontechnical challenges to advancing their use. Frameworks to support the rapid sharing of multidimensional data and increased cross-sector partnerships, networks, and collaborations can support advancement on all scales, from research and development to implementation by public health agencies.

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.

SummaryThere is an increasing role for bioinformatic and phylogenetic analysis in tropical medicine research. However, scientists working in low- and middle-income regions may lack access to training opportunities in these methods. To help address this gap, a 5-day intensive bioinformatics workshop was offered in Lima, Peru. The syllabus is presented here for others who want to develop similar programs. To assess knowledge gained, a 20-point knowledge questionnaire was administered to participants (21 participants) before and after the workshop, covering topics on sequence quality control, alignment/formatting, database retrieval, models of evolution, sequence statistics, tree building, and results interpretation. Evolution/tree-building methods represented the lowest scoring domain at baseline and after the workshop. There was a considerable median gain in total knowledge scores (increase of 30%, p

While studying rickettsial infections in Peru, we detected Rickettsia asembonensis in fleas from domestic animals. We characterized 5 complete genomic regions (17kDa, gltA, ompA, ompB, and sca4) and conducted multilocus sequence typing and phylogenetic analyses. The molecular isolate from Peru is distinct from the original R. asembonensis strain from Kenya.

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 living200 meters apart come from the same transmission chain, as opposed to 3% of cases separated by 1 to 5 kilometers. At distances200 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.

Here, we present the complete genome sequences of two Zika virus (ZIKV) strains, Zika virus/ Homo sapiens -tc/THA/2014/SV0127-14 and Zika virus/ H. sapiens -tc/PHL/2012/CPC-0740, isolated from the blood of patients collected in Thailand, 2014, and the Philippines, 2012, respectively. Sequencing and phylogenetic analysis showed that both strains belong to the Asian lineage.

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 (

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. Electronic supplementary material The online version of this article (doi:10.1186/s12859-015-0840-5) contains supplementary material, which is available to authorized users.

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. © 2011 Elsevier B.V.

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.

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.

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.

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°C to 63.4°C (average 60°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, N 2 fixation in hot spring mats was attributed either to filamentous cyanobacteria (not present at >50°C in these mats) or to heterotrophic bacteria. To explore how energy-generating processes of the Synechococcus ecotypes track natural light and O 2 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 N 2 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.

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. Electronic supplementary material The online version of this article (doi:10.1186/s12915-016-0337-3) contains supplementary material, which is available to authorized users.