Search

Summary Previous phylogenies, built using a subset of genomic loci, split Pseudomonas syringae pv. pisi into two well-supported clades and implied convergence in host range for these lineages. The analysis of phenotypic and genotypic data within the context of this phylogenetic relationship implied further convergence at the level of virulence gene loss and acquisition. We generate draft genome assemblies for two additional P. syringae strains, isolated from diseased pea plants, and demonstrate incongruence between phylogenies created from a subset of the data compared with the whole genomes. Our whole-genome analysis demonstrates that strains classified as pv. pisi actually form a coherent monophyletic clade, so that apparent convergence is actually the product of shared ancestry. We use this example to urge caution when making evolutionary inferences across closely related strains of P. syringae.

Rooting phylogenies is critical for understanding evolution, yet the importance, intricacies and difficulties of rooting are often overlooked. For rooting, polymorphic characters among the group of interest (ingroup) must be compared to those of a relative (outgroup) that diverged before the last common ancestor (LCA) of the ingroup. Problems arise if an outgroup does not exist, is unknown, or is so distant that few characters are shared, in which case duplicated genes originating before the LCA can be used as proxy outgroups to root diverse phylogenies. Here, we describe a genome-wide expansion of this technique that can be used to solve problems at the other end of the evolutionary scale: where ingroup individuals are all very closely related to each other, but the next closest relative is very distant. We used shared orthologous single nucleotide polymorphisms (SNPs) from 10 whole genome sequences of Coxiella burnetii, the causative agent of Q fever in humans, to create a robust, but unrooted phylogeny. To maximize the number of characters informative about the rooting, we searched entire genomes for polymorphic duplicated regions where orthologs of each paralog could be identified so that the paralogs could be used to root the tree. Recent radiations, such as those of emerging pathogens, often pose rooting challenges due to a lack of ingroup variation and large genomic differences with known outgroups. Using a phylogenomic approach, we created a robust, rooted phylogeny for C. burnetii. [Coxiella burnetii; paralog SNPs; pathogen evolution; phylogeny; recent radiation; root; rooting using duplicated genes.].

The bacterium Burkholderia pseudomallei causes melioidosis, a rare but serious illness that can be fatal if untreated or misdiagnosed. Species-specific PCR assays provide a technically simple method for differentiating B. pseudomallei from near-neighbor species. However, substantial genetic diversity and high levels of recombination within this species reduce the likelihood that molecular signatures will differentiate all B. pseudomallei from other Burkholderiaceae. Currently available molecular assays for B. pseudomallei detection lack rigorous validation across large in silico datasets and isolate collections to test for specificity, and none have been subjected to stringent quality control criteria (accuracy, precision, selectivity, limit of quantitation (LoQ), limit of detection (LoD), linearity, ruggedness and robustness) to determine their suitability for environmental, clinical or forensic investigations. In this study, we developed two novel B. pseudomallei specific assays, 122018 and 266152, using a dual-probe approach to differentiate B. pseudomallei from B. thailandensis, B. oklahomensis and B. thailandensis-like species; other species failed to amplify. Species specificity was validated across a large DNA panel (>2,300 samples) comprising Burkholderia spp. and non-Burkholderia bacterial and fungal species of clinical and environmental relevance. Comparison of assay specificity to two previously published B. pseudomallei-specific assays, BurkDiff and TTS1, demonstrated comparable performance of all assays, providing between 99.7 and 100% specificity against our isolate panel. Last, we subjected 122018 and 266152 to rigorous quality control analyses, thus providing quantitative limits of assay performance. Using B. pseudomallei as a model, our study provides a framework for comprehensive quantitative validation of molecular assays and provides additional, highly validated B. pseudomallei assays for the scientific research community.

Heroin may have been accidentally contaminated by an animal-derived source along a major drug trafficking route., In December 2009, two unusual cases of anthrax were diagnosed in heroin users in Scotland. A subsequent anthrax outbreak in heroin users emerged throughout Scotland and expanded into England and Germany, sparking concern of nefarious introduction of anthrax spores into the heroin supply. To better understand the outbreak origin, we used established genetic signatures that provided insights about strain origin. Next, we sequenced the whole genome of a representative Bacillus anthracis strain from a heroin user (Ba4599), developed Ba4599-specific single-nucleotide polymorphism assays, and genotyped all available material from other heroin users with anthrax. Of 34 case-patients with B. anthracis–positive PCR results, all shared the Ba4599 single-nucleotide polymorphism genotype. Phylogeographic analysis demonstrated that Ba4599 was closely related to strains from Turkey and not to previously identified isolates from Scotland or Afghanistan, the presumed origin of the heroin. Our results suggest accidental contamination along the drug trafficking route through a cutting agent or animal hides used to smuggle heroin into Europe.

Francisella tularensis subsp. holarctica isolates from Austria, Germany, Hungary, Italy, and Romania were placed into an existing phylogeographic framework. Isolates from Italy were assigned to phylogenetic group B.FTNF002-00; the other isolates, to group B.13. Most F. tularensis subsp. holarctica isolates from Europe belong to these 2 geographically segregated groups.

Next-generation sequencing enables use of whole-genome sequence typing (WGST) as a viable and discriminatory tool for genotyping and molecular epidemiologic analysis. We used WGST to confi rm the linkage of a cluster of Coccidioides immitis isolates from 3 patients who received organ transplants from a single donor who later had positive test results for coccidioidomycosis. Isolates from the 3 patients were nearly genetically identical (a total of 3 singlenucleotide polymorphisms identifi ed among them), thereby demonstrating direct descent of the 3 isolates from an original isolate. We used WGST to demonstrate the genotypic relatedness of C. immitis isolates that were also epidemiologically linked. Thus, WGST offers unique benefi ts to public health for investigation of clusters considered to be linked to a single source.

Advanced prostate cancer can progress to systemic metastatic tumors, which are generally androgen insensitive and ultimately lethal. Here, we report a comprehensive genomic survey for somatic events in systemic metastatic prostate tumors using both high-resolution copy number analysis and targeted mutational survey of 3508 exons from 577 cancer-related genes using next generation sequencing. Focal homozygous deletions were detected at 8p22, 10q23.31, 13q13.1, 13q14.11, and 13q14.12. Key genes mapping within these deleted regions include PTEN, BRCA2, C13ORF15, and SIAH3. Focal high-level amplifications were detected at 5p13.2-p12, 14q21.1, 7q22.1, and Xq12. Key amplified genes mapping within these regions include SKP2, FOXA1, and AR. Furthermore, targeted mutational analysis of normal-tumor pairs has identified somatic mutations in genes known to be associated with prostate cancer including AR and TP53, but has also revealed novel somatic point mutations in genes including MTOR, BRCA2, ARHGEF12, and CHD5. Finally, in one patient where multiple independent metastatic tumors were available, we show common and divergent somatic alterations that occur at both the copy number and point mutation level, supporting a model for a common clonal progenitor with metastatic tumor-specific divergence. Our study represents a deep genomic analysis of advanced metastatic prostate tumors and has revealed candidate somatic alterations, possibly contributing to lethal prostate cancer.

Francisella tularensis contains several highly pathogenic subspecies, including Francisella tularensis subsp. holarctica , whose distribution is circumpolar in the northern hemisphere. The phylogeography of these subspecies and their subclades was examined using whole-genome single nucleotide polymorphism (SNP) analysis, high-density microarray SNP genotyping, and real-time-PCR-based canonical SNP (canSNP) assays. Almost 30,000 SNPs were identified among 13 whole genomes for phylogenetic analysis. We selected 1,655 SNPs to genotype 95 isolates on a high-density microarray platform. Finally, 23 clade- and subclade-specific canSNPs were identified and used to genotype 496 isolates to establish global geographic genetic patterns. We confirm previous findings concerning the four subspecies and two Francisella tularensis subsp. tularensis subpopulations and identify additional structure within these groups. We identify 11 subclades within F. tularensis subsp. holarctica , including a new, genetically distinct subclade that appears intermediate between Japanese F. tularensis subsp. holarctica isolates and the common F. tularensis subsp. holarctica isolates associated with the radiation event (the B radiation) wherein this subspecies spread throughout the northern hemisphere. Phylogenetic analyses suggest a North American origin for this B-radiation clade and multiple dispersal events between North America and Eurasia. These findings indicate a complex transmission history for F. tularensis subsp. holarctica .

We identified nine naturally-occurring human single nucleotide polymorphisms (SNPs) in the alpha(1a)-adrenoceptor (alpha(1a)AR) coding region, seven of which result in amino acid change. Utilizing rat-1 fibroblasts stably expressing wild type alpha(1a)AR or each SNP at both high and low levels, we investigated the effect of these SNPs on receptor function. Compared with wild type, two SNPs (R166K, V311I) cause a decrease in binding affinity for agonists norepinephrine, epinephrine, and phenylephrine, and also shift the dose-response curve for norepinephrine stimulation of inositol phosphate (IP) production to the right (reduced potency) without altering maximal IP activity. In addition, SNP V311I and I200S display altered antagonist binding. Interestingly, a receptor with SNP G247R (located in the third intracellular loop) displays increased maximal receptor IP activity and stimulates cell growth. The increased receptor signaling for alpha(1a)AR G247R is not mediated by altered ligand binding or a deficiency in agonist-mediated desensitization, but appears to be related to enhanced receptor-G protein coupling. In conclusion, four naturally-occurring human alpha(1a)AR SNPs induce altered receptor pharmacology and/or biological activity. This finding has potentially important implications in many areas of medicine and can be used to guide alpha(1a)AR SNP choice for future clinical studies.

Duplications have long been postulated to be an important mechanism by which genomes evolve. Interspecies genomic comparisons are one method by which the origin and molecular mechanism of duplications can be inferred. By comparative mapping in human, mouse, and rat, we previously found evidence for a recent chromosome-fission event that occurred in the mouse lineage. Cytogenetic mapping revealed that the genomic segments flanking the fission site appeared to be duplicated, with copies residing near the centromere of multiple mouse chromosomes. Here we report the mapping and sequencing of the regions of mouse chromosomes 5 and 6 involved in this chromosome-fission event as well as the results of comparative sequence analysis with the orthologous human and rat genomic regions. Our data indicate that the duplications associated with mouse chromosomes 5 and 6 are recent and that the resulting duplicated segments share significant sequence similarity with a series of regions near the centromeres of the mouse chromosomes previously identified by cytogenetic mapping. We also identified pericentromeric duplicated segments shared between mouse chromosomes 5 and 1. Finally, novel mouse satellite sequences as well as putative chimeric transcripts were found to be associated with the duplicated segments. Together, these findings demonstrate that pericentromeric duplications are not restricted to primates and may be a common mechanism for genome evolution in mammals.[Supplemental material is available online at www.genome.org.]

Williams syndrome is a complex developmental disorder that results from the heterozygous deletion of a ∼1.6-Mb segment of human chromosome 7q11.23. These deletions are mediated by large (∼300 kb) duplicated blocks of DNA of near-identical sequence. Previously, we showed that the orthologous region of the mouse genome is devoid of such duplicated segments. Here, we extend our studies to include the generation of ∼3.3 Mb of genomic sequence from the mouse Williams syndrome region, of which just over 1.4 Mb is finished to high accuracy. Comparative analyses of the mouse and human sequences within and immediately flanking the interval commonly deleted in Williams syndrome have facilitated the identification of nine previously unreported genes, provided detailed sequence-based information regarding 30 genes residing in the region, and revealed a number of potentially interesting conserved noncoding sequences. Finally, to facilitate comparative sequence analysis, we implemented several enhancements to the program PipMaker, including the addition of links from annotated features within a generated percent-identity plot to specific records in public databases. Taken together, the results reported here provide an important comparative sequence resource that should catalyze additional studies of Williams syndrome, including those that aim to characterize genes within the commonly deleted interval and to develop mouse models of the disorder.[The sequence data described in this paper have been submitted to GenBank under accession nos. AF267747, AF289666,AF289667, AF289664, AF289665, AC091250, AC079938, AC084109, AC024607,AC074359, AC024608, AC083858, AC083948, AC084162, AC087420, AC083890,AC080158, AC084402, AC083889, AC083857, and AC079872.]

Sequence analyses and subtyping of Bacillus anthracis strains from Georgia reveal a single distinct lineage (Aust94) that is ecologically established. Phylogeographic analysis and comparisons to a global collection reveals a clade that is mostly restricted to Georgia. Within this clade, many groups are found around the country, however at least one subclade is only found in the eastern part. This pattern suggests that dispersal into and out of Georgia has been rare and despite historical dispersion within the country, for at least for one lineage, current spread is limited. © 2014 Khmaladze et al.

Brucella species include important zoonotic pathogens that have a substantial impact on both agriculture and human health throughout the world. Brucellae are thought of as “stealth pathogens” that escape recognition by the host innate immune response, modulate the acquired immune response, and evade intracellular destruction. We analyzed the genome sequences of members of the family Brucellaceae to assess its evolutionary history from likely free-living soil-based progenitors into highly successful intracellular pathogens. Phylogenetic analysis split the genus into two groups: recently identified and early-dividing “atypical” strains and a highly conserved “classical” core clade containing the major pathogenic species. Lateral gene transfer events brought unique genomic regions into Brucella that differentiated them from Ochrobactrum and allowed the stepwise acquisition of virulence factors that include a type IV secretion system, a perosamine-based O antigen, and systems for sequestering metal ions that are absent in progenitors. Subsequent radiation within the core Brucella resulted in lineages that appear to have evolved within their preferred mammalian hosts, restricting their virulence to become stealth pathogens capable of causing long-term chronic infections.

Prairie dogs (Cynomys spp.) are highly susceptible to Yersinia pestis, with ≥99% mortality reported from multiple studies of plague epizootics. A colony of Gunnison's prairie dogs (Cynomys gunnisoni) in the Aubrey Valley (AV) of northern Arizona appears to have survived several regional epizootics of plague, whereas nearby colonies have been severely affected by Y. pestis. To examine potential mechanisms accounting for survival in the AV colony, we conducted a laboratory Y. pestis challenge experiment on 60 wild-caught prairie dogs from AV and from a nearby, large colony with frequent past outbreaks of plague, Espee (n = 30 per colony). Test animals were challenged subcutaneously with the fully virulent Y. pestis strain CO92 at three doses: 50, 5,000, and 50,000 colony-forming units (cfu); this range is lethal in black-tailed prairie dogs (Cynomys ludovicianus). Contrary to our expectations, only 40% of the animals died. Although mortality trended higher in the Espee colony (50%) compared with AV (30%), the differences among infectious doses were not statistically significant. Only 39% of the survivors developed moderate to high antibody levels to Y. pestis, indicating that mechanisms other than humoral immunity are important in resistance to plague. The ratio of neutrophils to lymphocytes was not correlated with plague survival in this study. However, several immune proteins with roles in innate immunity (VCAM-1, CXCL-1, and vWF) were upregulated during plague infection and warrant further inquiry into their role for protection against this disease. These results suggest plague resistance exists in wild populations of the Gunnison's prairie dog and provide important directions for future studies.

Burkholderia pseudomallei causes the potentially fatal disease melioidosis. It is generally accepted that B. pseudomallei is a noncommensal bacterium and that any culture-positive clinical specimen denotes disease requiring treatment. Over a 23-year study of melioidosis cases in Darwin, Australia, just one patient from 707 survivors has developed persistent asymptomatic B. pseudomallei carriage. To better understand the mechanisms behind this unique scenario, we performed whole-genome analysis of two strains isolated 139 months apart. During this period, B. pseudomallei underwent several adaptive changes. Of 23 point mutations, 78% were nonsynonymous and 43% were predicted to be deleterious to gene function, demonstrating a strong propensity for positive selection. Notably, a nonsense mutation inactivated the universal stress response sigma factor RpoS, with pleiotropic implications. The genome underwent substantial reduction, with four deletions in chromosome 2 resulting in the loss of 221 genes. The deleted loci included genes involved in secondary metabolism, environmental survival, and pathogenesis. Of 14 indels, 11 occurred in coding regions and 9 resulted in frameshift mutations that dramatically affected predicted gene products. Disproportionately, four indels affected lipopolysaccharide biosynthesis and modification. Finally, we identified a frameshift mutation in both P314 isolates within wcbR, an important component of the capsular polysaccharide I locus, suggesting virulence attenuation early in infection. Our study illustrates a unique clinical case that contrasts a high-consequence infectious agent with a long-term commensal infection and provides further insights into bacterial evolution within the human host., IMPORTANCE Some bacterial pathogens establish long-term infections that are difficult or impossible to eradicate with current treatments. Rapid advances in genome sequencing technologies provide a powerful tool for understanding bacterial persistence within the human host. Burkholderia pseudomallei is considered a highly pathogenic bacterium because infection is commonly fatal. Here, we document within-host evolution of B. pseudomallei in a unique case of human infection with ongoing chronic carriage. Genomic comparison of isolates obtained 139 months (11.5 years) apart showed a strong signal of adaptation within the human host, including inactivation of virulence and immunogenic factors, and deletion of pathways involved in environmental survival. Two global regulatory genes were mutated in the 139-month isolate, indicating extensive regulatory changes favoring bacterial persistence. Our study provides insights into B. pseudomallei pathogenesis and, more broadly, identifies parallel evolutionary mechanisms that underlie chronic persistence of all bacterial pathogens.

A cluster of human plague cases occurred in the seaport city of Mahajanga, Madagascar, from 1991 to 1999 following 62 years with no evidence of plague, which offered insights into plague pathogen dynamics in an urban environment. We analyzed a set of 44 Mahajanga isolates from this 9-year outbreak, as well as an additional 218 Malagasy isolates from the highland foci. We sequenced the genomes of four Mahajanga strains, performed whole-genome sequence single-nucleotide polymorphism (SNP) discovery on those strains, screened the discovered SNPs, and performed a high-resolution 43-locus multilocus variable-number tandem-repeat analysis of the isolate panel. Twenty-two new SNPs were identified and defined a new phylogenetic lineage among the Malagasy isolates. Phylogeographic analysis suggests that the Mahajanga lineage likely originated in the Ambositra district in the highlands, spread throughout the northern central highlands, and was then introduced into and became transiently established in Mahajanga. Although multiple transfers between the central highlands and Mahajanga occurred, there was a locally differentiating and dominant subpopulation that was primarily responsible for the 1991-to-1999 Mahajanga outbreaks. Phylotemporal analysis of this Mahajanga subpopulation revealed a cycling pattern of diversity generation and loss that occurred during and after each outbreak. This pattern is consistent with severe interseasonal genetic bottlenecks along with large seasonal population expansions. The ultimate extinction of plague pathogens in Mahajanga suggests that, in this environment, the plague pathogen niche is tenuous at best. However, the temporary large pathogen population expansion provides the means for plague pathogens to disperse and become ecologically established in more suitable nonurban environments., IMPORTANCE Maritime spread of plague led to the global dissemination of this disease and affected the course of human history. Multiple historical plague waves resulted in massive human mortalities in three classical plague pandemics: Justinian (6th and 7th centuries), Middle Ages (14th to 17th centuries), and third (mid-1800s to the present). Key to these events was the pathogen’s entry into new lands by “plague ships” via seaport cities. Although initial disease outbreaks in ports were common, they were almost never sustained for long and plague pathogens survived only if they could become established in ecologically suitable habitats. Although plague pathogens’ ability to invade port cities has been essential for intercontinental spread, these regions have not proven to be a suitable long-term niche. The disease dynamics in port cities such as Mahajanga are thus critical to plague pathogen amplification and dispersal into new suitable ecological niches for the observed global long-term maintenance of plague pathogens.

Because thyroid hormones (THs) are conserved modulators of development and physiology, identification of compounds adversely affecting TH signaling is critical to human and wildlife health. Anurans are an established model for studying disruption of TH signaling because metamorphosis is dependent upon the thyroid system. In order to strengthen this model and identify new gene transcript biomarkers for TH disruption, we performed DNA microarray analysis of Xenopus laevis tadpole tail transcriptomes following treatment with triiodothyronine (T 3 ). Comparison of these results with previous studies in frogs and mammals identified 36 gene transcripts that were TH-sensitive across clades. We then tested molecular biomarkers for sensitivity to disruption by exposure to wastewater effluent (WWE). X. laevis tadpoles, exposed to WWE from embryo through metamorphosis, exhibited an increased developmental rate compared to controls. Cultured tadpole tails showed dramatic increases in levels of four TH-sensitive gene transcripts (thyroid hormone receptor β (TRβ), deiodinase type II (DIO2), and corticotropin releasing hormone binding protein (CRHBP), fibroblast activation protein α (FAPα)) when exposed to T 3 and WWE extracts. TRβ, DIO2, and CRHBP were identified as TH sensitive in other studies, while FAPα mRNA transcripts were highly TH sensitive in our array. The results validate the array and demonstrate TH-disrupting activity by WWE. Our findings demonstrate the usefulness of cross-clade analysis for identification of gene transcripts that provide sensitivity to endocrine disruption. Further, the results suggest that development is disrupted by exposure to complex mixes of compounds found in WWE possibly through interference with TH signaling.

Melioidosis is caused by Burkholderia pseudomallei, a Gram-negative bacillus, primarily found in soils in Southeast Asia and northern Australia. A recent case of melioidosis in non-endemic Arizona was determined to be the result of locally acquired infection, as the patient had no travel history to endemic regions and no previous history of disease. Diagnosis of the case was confirmed through multiple microbiologic and molecular techniques. To enhance the epidemiological analysis, we conducted several molecular genotyping procedures, including multi-locus sequence typing, SNP-profiling, and whole genome sequence typing. Each technique has different molecular epidemiologic advantages, all of which provided evidence that the infecting strain was most similar to those found in Southeast Asia, possibly originating in, or around, Malaysia. Advancements in new typing technologies provide genotyping resolution not previously available to public health investigators, allowing for more accurate source identification., Author Summary Melioidosis is a bacterial disease caused by percutaneous inoculation, aspiration or ingestion of the soil bacteria Burkholderia pseudomallei. Melioidosis is primarily found in Southeast Asia and Northern Australia, and, to a lesser degree, nearby regions. A recent case of melioidosis in Southwestern United States (Southern Arizona) prompted a detailed epidemiological and molecular investigation to discover the source of infection. The authors describe the use of multiple genomic analysis tools to aid in this investigation. The results of these analyses uniformly identified Southeast Asia as the source of the strain that infected the patient, however the epidemiological investigation had determined the patient had no international travel or known exposures to Southeast Asian products. New cutting edge technologies, such as next generation sequencing, are quickly being adapted into epidemiologic investigations, particularly for cases and outbreaks of unknown origin, although older, mature technologies with larger existing databases will still be needed for appropriate comparative analyses.

Summary Genomic structural variants (SVs) are abundant in humans, differing from other variation classes in extent, origin, and functional impact. Despite progress in SV characterization, the nucleotide resolution architecture of most SVs remains unknown. We constructed a map of unbalanced SVs (i.e., copy number variants) based on whole genome DNA sequencing data from 185 human genomes, integrating evidence from complementary SV discovery approaches with extensive experimental validations. Our map encompassed 22,025 deletions and 6,000 additional SVs, including insertions and tandem duplications. Most SVs (53%) were mapped to nucleotide resolution, which facilitated analyzing their origin and functional impact. We examined numerous whole and partial gene deletions with a genotyping approach and observed a depletion of gene disruptions amongst high frequency deletions. Furthermore, we observed differences in the size spectra of SVs originating from distinct formation mechanisms, and constructed a map constructed a map of SV hotspots formed by common mechanisms. Our analytical framework and SV map serves as a resource for sequencing-based association studies.

Background Plague was introduced to Madagascar in 1898 and continues to be a significant human health problem. It exists mainly in the central highlands, but in the 1990s was reintroduced to the port city of Mahajanga, where it caused extensive human outbreaks. Despite its prevalence, the phylogeography and molecular epidemiology of Y. pestis in Madagascar has been difficult to study due to the great genetic similarity among isolates. We examine island-wide geographic-genetic patterns based upon whole-genome discovery of SNPs, SNP genotyping and hypervariable variable-number tandem repeat (VNTR) loci to gain insight into the maintenance and spread of Y. pestis in Madagascar. Methodology/Principal Findings We analyzed a set of 262 Malagasy isolates using a set of 56 SNPs and a 43-locus multi-locus VNTR analysis (MLVA) system. We then analyzed the geographic distribution of the subclades and identified patterns related to the maintenance and spread of plague in Madagascar. We find relatively high levels of VNTR diversity in addition to several SNP differences. We identify two major groups, Groups I and II, which are subsequently divided into 11 and 4 subclades, respectively. Y. pestis appears to be maintained in several geographically separate subpopulations. There is also evidence for multiple long distance transfers of Y. pestis, likely human mediated. Such transfers have resulted in the reintroduction and establishment of plague in the port city of Mahajanga, where there is evidence for multiple transfers both from and to the central highlands. Conclusions/Significance The maintenance and spread of Y. pestis in Madagascar is a dynamic and highly active process that relies on the natural cycle between the primary host, the black rat, and its flea vectors as well as human activity., Author Summary Plague, caused by the bacterium Yersinia pestis, has been a problem in Madagascar since it was introduced in 1898. It mainly affects the central highlands, but also has caused several large outbreaks in the port city of Mahajanga, after it was reintroduced there in the 1990s. Despite its prevalence, the genetic diversity and related geographic distribution of different genetic groups of Y. pestis in Madagascar has been difficult to study due to the great genetic similarity among isolates. We subtyped a set of Malagasy isolates and identified two major genetic groups that were subsequently divided into 11 and 4 subgroups, respectively. Y. pestis appears to be maintained in several geographically separate subpopulations. There is also evidence for multiple long distance transfers of Y. pestis, likely human mediated. Such transfers have resulted in the reintroduction and establishment of plague in the port city of Mahajanga where there is evidence for multiple transfers both from and to the central highlands. The maintenance and spread of Y. pestis in Madagascar is a dynamic and highly active process that relies on the natural cycle between the primary host, the black rat, and its flea vectors as well as human activity.

Coccidioidomycosis is an infection caused by Coccidioides immitis or C. posadasii. We developed a TaqMan real-time PCR assay that rapidly and accurately differentiates the species. This assay can be used as a tool to improve disease surveillance, increase understanding of the natural history of the infection, and assist in clinical differentiation studies.

Burkholderia mallei (Bm), the causative agent of the predominately equine disease glanders, is a genetically uniform species that is very closely related to the much more diverse species Burkholderia pseudomallei (Bp), an opportunistic human pathogen and the primary cause of melioidosis. To gain insight into the relative lack of genetic diversity within Bm, we performed whole-genome comparative analysis of seven Bm strains and contrasted these with eight Bp strains. The Bm core genome (shared by all seven strains) is smaller in size than that of Bp, but the inverse is true for the variable gene sets that are distributed across strains. Interestingly, the biological roles of the Bm variable gene sets are much more homogeneous than those of Bp. The Bm variable genes are found mostly in contiguous regions flanked by insertion sequence (IS) elements, which appear to mediate excision and subsequent elimination of groups of genes that are under reduced selection in the mammalian host. The analysis suggests that the Bm genome continues to evolve through random IS-mediated recombination events, and differences in gene content may contribute to differences in virulence observed among Bm strains. The results are consistent with the view that Bm recently evolved from a single strain of Bp upon introduction into an animal host followed by expansion of IS elements, prophage elimination, and genome rearrangements and reduction mediated by homologous recombination across IS elements.

Background Phylogeographic reconstruction of some bacterial populations is hindered by low diversity coupled with high levels of lateral gene transfer. A comparison of recombination levels and diversity at seven housekeeping genes for eleven bacterial species, most of which are commonly cited as having high levels of lateral gene transfer shows that the relative contributions of homologous recombination versus mutation for Burkholderia pseudomallei is over two times higher than for Streptococcus pneumoniae and is thus the highest value yet reported in bacteria. Despite the potential for homologous recombination to increase diversity, B. pseudomallei exhibits a relative lack of diversity at these loci. In these situations, whole genome genotyping of orthologous shared single nucleotide polymorphism loci, discovered using next generation sequencing technologies, can provide very large data sets capable of estimating core phylogenetic relationships. We compared and searched 43 whole genome sequences of B. pseudomallei and its closest relatives for single nucleotide polymorphisms in orthologous shared regions to use in phylogenetic reconstruction. Results Bayesian phylogenetic analyses of >14,000 single nucleotide polymorphisms yielded completely resolved trees for these 43 strains with high levels of statistical support. These results enable a better understanding of a separate analysis of population differentiation among >1,700 B. pseudomallei isolates as defined by sequence data from seven housekeeping genes. We analyzed this larger data set for population structure and allele sharing that can be attributed to lateral gene transfer. Our results suggest that despite an almost panmictic population, we can detect two distinct populations of B. pseudomallei that conform to biogeographic patterns found in many plant and animal species. That is, separation along Wallace's Line, a biogeographic boundary between Southeast Asia and Australia. Conclusion We describe an Australian origin for B. pseudomallei, characterized by a single introduction event into Southeast Asia during a recent glacial period, and variable levels of lateral gene transfer within populations. These patterns provide insights into mechanisms of genetic diversification in B. pseudomallei and its closest relatives, and provide a framework for integrating the traditionally separate fields of population genetics and phylogenetics for other bacterial species with high levels of lateral gene transfer.

Helianthus petiolaris andH. niveus are polytypic species which are morphologically distinct at the periphery of their ranges but intergrade in areas of sympatry.Helianthus niveus includes both annual and perennial members, whereasH. petiolaris is strictly annual. Chloroplast DNA and nuclear ribosomal DNA restriction site data were used to reconstruct the evolutionary history of populations of the two species. Cladistic analyses reveal the following: (1) neither species is monophyletic; (2) the annual habit is derived once in this complex; and (3) the region of morphological intergradation appears to be primary in origin. The significance of interbreeding versus common descent in defining species concepts is discussed in relation to the above cladistic analyses.

Background Burkholderia pseudomallei is the etiologic agent of melioidosis, a significant cause of morbidity and mortality where this infection is endemic. Genomic differences among strains of B. pseudomallei are predicted to be one of the major causes of the diverse clinical manifestations observed among patients with melioidosis. The purpose of this study was to examine the role of genomic islands (GIs) as sources of genomic diversity in this species. Results We found that genomic islands (GIs) vary greatly among B. pseudomallei strains. We identified 71 distinct GIs from the genome sequences of five reference strains of B. pseudomallei: K96243, 1710b, 1106a, MSHR668, and MSHR305. The genomic positions of these GIs are not random, as many of them are associated with tRNA gene loci. In particular, the 3' end sequences of tRNA genes are predicted to be involved in the integration of GIs. We propose the term "tRNA-mediated site-specific recombination" (tRNA-SSR) for this mechanism. In addition, we provide a GI nomenclature that is based upon integration hotspots identified here or previously described. Conclusion Our data suggest that acquisition of GIs is one of the major sources of genomic diversity within B. pseudomallei and the molecular mechanisms that facilitate horizontally-acquired GIs are common across multiple strains of B. pseudomallei. The differential presence of the 71 GIs across multiple strains demonstrates the importance of these mobile elements for shaping the genetic composition of individual strains and populations within this bacterial species.

BackgroundYersinia pestis, the causative agent of plague, is responsible for some of the greatest epidemic scourges of mankind. It is widespread in the western United States, although it has only been present there for just over 100 years. As a result, there has been very little time for diversity to accumulate in this region. Much of the diversity that has been detected among North American isolates is at loci that mutate too quickly to accurately reconstruct large-scale phylogenetic patterns. Slowly-evolving but stable markers such as SNPs could be useful for this purpose, but are difficult to identify due to the monomorphic nature of North American isolates.Methodology/principal findingsTo identify SNPs that are polymorphic among North American populations of Y. pestis, a gapped genome sequence of Y. pestis strain FV-1 was generated. Sequence comparison of FV-1 with another North American strain, CO92, identified 19 new SNP loci that differ among North American isolates.Conclusions/significanceThe 19 SNP loci identified in this study should facilitate additional studies of the genetic population structure of Y. pestis across North America.

The systematic comparison of genomic sequences from different organisms represents a central focus of contemporary genome analysis. Comparative analyses of vertebrate sequences can identify coding1,2,3,4,5,6 and conserved non-coding4,6,7 regions, including regulatory elements8,9,10, and provide insight into the forces that have rendered modern-day genomes6. As a complement to whole-genome sequencing efforts3,5,6, we are sequencing and comparing targeted genomic regions in multiple, evolutionarily diverse vertebrates. Here we report the generation and analysis of over 12 megabases (Mb) of sequence from 12 species, all derived from the genomic region orthologous to a segment of about 1.8 Mb on human chromosome 7 containing ten genes, including the gene mutated in cystic fibrosis. These sequences show conservation reflecting both functional constraints and the neutral mutational events that shaped this genomic region. In particular, we identify substantial numbers of conserved non-coding segments beyond those previously identified experimentally, most of which are not detectable by pair-wise sequence comparisons alone. Analysis of transposable element insertions highlights the variation in genome dynamics among these species and confirms the placement of rodents as a sister group to the primates.

Background The zoonosis brucellosis causes economically significant reproductive problems in livestock and potentially debilitating disease of humans. Although the causative agent, organisms from the genus Brucella, can be differentiated into a number of species based on phenotypic characteristics, there are also significant differences in genotype that are concordant with individual species. This paper describes the development of a five target multiplex assay to identify five terrestrial Brucella species using real-time polymerase chain reaction (PCR) and subsequent high resolution melt curve analysis. This technology offers a robust and cost effective alternative to previously described hydrolysis-probe Single Nucleotide Polymorphism (SNP)-based species defining assays. Results Through the use of Brucella whole genome sequencing five species defining SNPs were identified. Individual HRM assays were developed to these target these changes and, following optimisation of primer concentrations, it was possible to multiplex all five assays in a single tube. In a validation exercise using a panel of 135 Brucella strains of terrestrial and marine origin, it was possible to distinguish the five target species from the other species within this panel. Conclusion The HRM multiplex offers a number of diagnostic advantages over previously described SNP-based typing approaches. Further, and uniquely for HRM, the successful multiplexing of five assays in a single tube allowing differentiation of five Brucella species in the diagnostic laboratory in a cost-effective and timely manner is described. However there are possible limitations to using this platform on DNA extractions direct from clinical material. Electronic supplementary material The online version of this article (doi:10.1186/1756-0500-7-903) contains supplementary material, which is available to authorized users.

Nephropathic cystinosis is a rare autosomal recessive, lysosomal storage disease with an incidence estimated at 1 per 100,000–200,000 live births (see http://www.ncbi.nlm.nih.gov/omim; OMIM 219800). The classic disorder is characterized clinically by renal tubular Fanconi syndrome in the first year of life, growth retardation in childhood, renal glomerular failure at ∼10 years of age, hypothyroidism, and a variety of other complications, including photophobia and corneal crystal formation (Gahl 1986; Gahl et al. 1995). After renal transplantation, cystine accumulation continues in nonrenal organs, frequently causing a distal vacuolar myopathy (Charnas et al. 1994), swallowing difficulty (Sonies et al. 1990), or retinal dysfunction (Kaiser-Kupfer et al. 1986), and occasionally causing diabetes mellitus (Fivush et al. 1987), pancreatic exocrine insufficiency (Fivush et al. 1988), or neurological deterioration (Ehrich et al. 1979; Fink et al. 1989). These complications arise because defective lysosomal transport of the disulfide cystine (Gahl et al. 1982a) causes this amino acid to accumulate within the lysosomes of many different cell types, which then triggers cystine crystal formation (Gahl et al. 1982b). The cystine transporter is the first of many lysosomal membrane carriers to be characterized biochemically (Thoene 1992), and cystinosis is the most common of a group of lysosomal transport disorders (Gahl et al. 1995). The gene altered in patients with cystinosis (CTNS) was recently identified by a positional cloning strategy (Town et al. 1998). CTNS is a 12-exon gene that is transcribed into a ∼2.6-kb mRNA. The encoded protein, named cystinosin, consists of a predicted 367 amino acids, appears to be an integral membrane protein, and most likely functions as a cystine transporter. A number of cystinosis-causing CTNS mutations have now been reported (Shotelersuk et al. 1998a; Town et al. 1998). The most prevalent mutation reported to date is a large (>55-kb) deletion, with 33%–44% of affected patients being homozygous for this deletion (Town et al. 1998; Anikster et al. 1999). In addition, at least 11 other smaller disease-causing deletions have been reported (Shotelersuk et al. 1998a; Forestier et al. 1999), suggesting that this genomic region may be prone to rearrangement. We sought to establish the long-range organization of the segment of chromosome 17p13 harboring CTNS and to determine the sequence of this clinically important gene and its surrounding DNA. Here we report the assembly of a detailed bacterial artificial chromosome (BAC)-based physical map encompassing CTNS. In addition, two BAC clones spanning >200 kb were sequenced to high accuracy, providing insight into the molecular architecture of the CTNS gene and the genomic segment commonly deleted in cystinosis patients.

The identification of the cystic fibrosis transmembrane conductance regulator gene ( CFTR ) in 1989 represents a landmark accomplishment in human genetics. Since that time, there have been numerous advances in elucidating the function of the encoded protein and the physiological basis of cystic fibrosis. However, numerous areas of cystic fibrosis biology require additional investigation, some of which would be facilitated by information about the long-range sequence context of the CFTR gene. For example, the latter might provide clues about the sequence elements responsible for the temporal and spatial regulation of CFTR expression. We thus sought to establish the sequence of the chromosomal segments encompassing the human CFTR and mouse Cftr genes, with the hope of identifying conserved regions of biologic interest by sequence comparison. Bacterial clone-based physical maps of the relevant human and mouse genomic regions were constructed, and minimally overlapping sets of clones were selected and sequenced, eventually yielding ≈1.6 Mb and ≈358 kb of contiguous human and mouse sequence, respectively. These efforts have produced the complete sequence of the ≈189-kb and ≈152-kb segments containing the human CFTR and mouse Cftr genes, respectively, as well as significant amounts of flanking DNA. Analyses of the resulting data provide insights about the organization of the CFTR / Cftr genes and potential sequence elements regulating their expression. Furthermore, the generated sequence reveals the precise architecture of genes residing near CFTR / Cftr , including one known gene ( WNT2 / Wnt2 ) and two previously unknown genes that immediately flank CFTR / Cftr .

Background: Brucellosis is a worldwide disease of mammals caused by Alphaproteobacteria in the genus Brucella. The genus is genetically monomorphic, requiring extensive genotyping to differentiate isolates. We utilized two different genotyping strategies to characterize isolates. First, we developed a microarray-based assay based on 1000 single nucleotide polymorphisms (SNPs) that were identified from whole genome comparisons of two B. abortus isolates, one B. melitensis, and one B. suis. We then genotyped a diverse collection of 85 Brucella strains at these SNP loci and generated a phylogenetic tree of relationships. Second, we developed a selective primer-extension assay system using capillary electrophoresis that targeted 17 high value SNPs across 8 major branches of the phylogeny and determined their genotypes in a large collection (n=340) of diverse isolates. Results: Our 1000 SNP microarray readily distinguished B. abortus, B. melitensis, and B. suis, differentiating B. melitensis and B. suis into two clades each. Brucella abortus was divided into four major clades. Our capillary-based SNP genotyping confirmed all major branches from the microarray assay and assigned all samples to defined lineages. Isolates from these lineages and closely related isolates, among the most commonly encountered lineages worldwide, can now be quickly and easily identified and genetically characterized. Conclusions: We have identified clade-specific SNPs in Brucella that can be used for rapid assignment into major groups below the species level in the three main Brucella species. Our assays represent SNP genotyping approaches that can reliably determine the evolutionary relationships of bacterial isolates without the need for whole genome sequencing of all isolates. Background Brucella are Gram-negative bacteria and the causative agent of brucellosis in domesticated animals, wildlife, and humans. Although the bacteria exhibit relatively strong host preference, separating the various Brucella species has proven extremely difficult due to minimal genetic differentiation [1]. Based on DNA-DNA hybridization, brucellae have previously been proposed to comprise a single species, with a series of biovars [2]. However, phylogenetic approaches explicitly incorporating host preference and virulence have upheld the six classical Brucella species: B. abortus (bovine), B. melitensis (caprine and ovine), B. suis (porcine), B. canis (canine), B. neotomae (desert woodrat), and B. ovis

Castor bean is an important oil-producing plant in the Euphorbiaceae family. Its high-quality oil contains up to 90% of the unusual fatty acid ricinoleate, which has many industrial and medical applications. Castor bean seeds also contain ricin, a highly toxic Type 2 ribosome-inactivating protein, which has gained relevance in recent years due to biosafety concerns. In order to gain knowledge on global genetic diversity in castor bean and to ultimately help the development of breeding and forensic tools, we carried out an extensive chloroplast sequence diversity analysis. Taking advantage of the recently published genome sequence of castor bean, we assembled the chloroplast and mitochondrion genomes extracting selected reads from the available whole genome shotgun reads. Using the chloroplast reference genome we used the methylation filtration technique to readily obtain draft genome sequences of 7 geographically and genetically diverse castor bean accessions. These sequence data were used to identify single nucleotide polymorphism markers and phylogenetic analysis resulted in the identification of two major clades that were not apparent in previous population genetic studies using genetic markers derived from nuclear DNA. Two distinct sub-clades could be defined within each major clade and large-scale genotyping of castor bean populations worldwide confirmed previously observed low levels of genetic diversity and showed a broad geographic distribution of each sub-clade.

The phytochemical/ecogeographic database (PHYTOCHEMECO) of the National Germplasm Resources Laboratory (NGRL), U.S. Department of Agriculture, is a unique blend of phytochemical, taxonomic, ecological, geographic and climatic aspects, and yields information and germplasm databases in a UNIX minicomputer environment. The phytochemical database contains information about 16,332 chemical compounds found in 1,164 plants of economic importance, including quantitative data where available, as well as localization and specific activities (1,517 activities of 4,329 chemicals) of the compounds.1,2 The taxonomic, ecological and yield databases contain plant names of 8,771 taxa, growing locations of 6,242 taxa, and crop yields of 239 taxa, respectively. The geographic and climatic3 databases contain specific information about more than 18,000 world-wide locations, including latitude, longitude, elevation, soil type, rainfall, temperature and lifezones. The Germplasm Resources Information Network database (GRIN) lists accessions of plant clonal and seed material available for study through the National Plant Germplasm System.

Background Burkholderia pseudomallei is intrinsically resistant to aminoglycosides and macrolides, mostly due to AmrAB-OprA efflux pump expression. We investigated the molecular mechanisms of aminoglycoside susceptibility exhibited by Thai strains 708a, 2188a, and 3799a. Methodology/Principal Findings qRT-PCR revealed absence of amrB transcripts in 708a and greatly reduced levels in 2188a and 3799a. Serial passage on increasing gentamicin concentrations yielded 2188a and 3799a mutants that became simultaneously resistant to other aminoglycosides and macrolides, whereas such mutants could not be obtained with 708a. Transcript analysis showed that the resistance of the 2188a and 3799a mutants was due to upregulation of amrAB-oprA expression by unknown mechanism(s). Use of a PCR walking strategy revealed that the amrAB-oprA operon was missing in 708a and that this loss was associated with deletion of more than 70 kb of genetic material. Rescue of the amrAB-oprB region from a 708a fosmid library and sequencing showed the presence of a large chromosome 1 deletion (131 kb and 141 kb compared to strains K96243 and 1710b, respectively). This deletion not only removed the amrAB-oprA operon, but also the entire gene clusters for malleobactin and cobalamin synthesis. Other genes deleted included the anaerobic arginine deiminase pathway, putative type 1 fimbriae and secreted chitinase. Whole genome sequencing and PCR analysis confirmed absence of these genes from 708a. Despite missing several putative virulence genes, 708a was fully virulent in a murine melioidosis model. Conclusions/Significance Strain 708a may be a natural candidate for genetic manipulation experiments that use Select Agent compliant antibiotics for selection and validates the use of laboratory-constructed Δ(amrAB-oprA) mutants in such experiments., Author Summary Burkholderia pseudomallei is the etiologic agent of melioidosis, an emerging tropical disease. Because of low infectious dose, broad-host-range infectivity, intrinsic antibiotic resistance and historic precedent as a bioweapon, B. pseudomallei was listed in the United States as a Select Agent and Priority Pathogen of biodefense concern by the US Centers for Disease Control and Prevention and the National Institute of Allergy and Infectious Diseases. The mechanisms governing antibiotic resistance and/or susceptibility and virulence in this bacterium are not well understood. Most clinical and environmental B. pseudomallei isolates are highly resistant to aminoglycosides, but susceptible variants do exist. The results of our studies with three such variants from Thailand reveal that lack of expression or deletion of an efflux pump is responsible for this susceptibility. The large deletion present in one strain not only removes an efflux pump but also several putative virulence genes, including an entire siderophore gene cluster. Despite this deletion, the strain is fully virulent in an acute mouse melioidosis model. In summary, our findings shed light on mechanisms of antibiotic resistance and pathogenesis. They also validate the previously advocated use of laboratory-constructed, aminoglycoside susceptible efflux pump mutants in genetic manipulation experiments.

Francisella tularensis is the causative agent of tularemia, which is a highly lethal disease from nature and potentially from a biological weapon. This species contains four recognized subspecies including the North American endemic F. tularensis subsp. tularensis (type A), whose genetic diversity is correlated with its geographic distribution including a major population subdivision referred to as A.I and A.II. The biological significance of the A.I - A.II genetic differentiation is unknown, though there are suggestive ecological and epidemiological correlations. In order to understand the differentiation at the genomic level, we have determined the complete sequence of an A.II strain (WY96-3418) and compared it to the genome of Schu S4 from the A.I population. We find that this A.II genome is 1,898,476 bp in size with 1,820 genes, 1,303 of which code for proteins. While extensive genomic variation exists between "WY96" and Schu S4, there is only one whole gene difference. This one gene difference is a hypothetical protein of unknown function. In contrast, there are numerous SNPs (3,367), small indels (1,015), IS element differences (7) and large chromosomal rearrangements (31), including both inversions and translocations. The rearrangement borders are frequently associated with IS elements, which would facilitate intragenomic recombination events. The pathogenicity island duplicated regions (DR1 and DR2) are essentially identical in WY96 but vary relative to Schu S4 at 60 nucleotide positions. Other potential virulence-associated genes (231) varied at 559 nucleotide positions, including 357 non-synonymous changes. Molecular clock estimates for the divergence time between A.I and A.II genomes for different chromosomal regions ranged from 866 to 2131 years before present. This paper is the first complete genomic characterization of a member of the A.II clade of Francisella tularensis subsp. tularensis.

Chloroplast DNA (cpDNA) and isozyme data were used to reconstruct the evolu- tionary history of the 21 taxa comprising Helianthus sect Helianthus. Low levels of cpDNA and isozyme divergence suggest a recent origin, perhaps within the last 1-2 million years, with subsequent rapid evolution and diversification throughout North America. This, combined with evolutionary phe- nomena such as phylogenetic sorting and introgression, makes phylogenetic reconstruction difficult. A strict consensus tree derived from Wagner parsimony analysis of the isozyme data set showed almost complete lack of resolution. The three clades resolved were, however, largely consistent with previous phylogenetic hypotheses. The most parsimonious cpDNA-based Wagner tree was poorly resolved, but showed that sect. Helianthus as presently circumscribed is monophyletic. Sur- prisingly, the cpDNA-based phylogeny does not resemble any of the previous phylogenetic hy- potheses for this group based on evidence from morphology, crossability, sesquiterpene lactone chemistry, and inferred chromosomal end arrangements. In general, species geographically most proximal are most closely related in terms of cpDNA, maternal inheritance of which is demonstrated herein. All four polytypic species in sect. Helianthus are polyphyletic in terms of cpDNA. Morpho- logical classification and cpDNA genotype were discordant for populations or species of H. annuus, H. debilis subsp. cucumerifolius, H. neglectus, H. petiolaris subsp. fallax, H. petiolaris subsp. petiolaris, and H. anomalus. Cytoplasmic introgression appears to be the most parsimonius explanation for these discrepancies. The edaphic differentiation and generally allopatric distribution of most taxa in sect. Helianthus are consistent with an allopatric mode of speciation. Molecular and chromosomal evi- dence, however, suggest that allopatric, parapatric, and quantum speciation all operate. In addition, at least one geographic race may have been derived through introgression and three species appear to be derived via selection or random fixation of recombinant genotypes following interspecific hybridization. We also tested the hypothesis that levels of genetic polymorphism will be lower in geographically restricted species than in more widespre-ad species. Eight of the 11 narrow endemics analyzed had lower levels of genetic diversity than any of their more widespread congeners.

Two-dimensional polyacrylamide gel (2-D PAGE) electrophoresis was appraised as an experimental technique for assessing systematic relationships among higher plants and to determine at which level in the taxonomic hierarchy this technique is most generally applicable. 2-D PAGE was performed on denatured extracts of mature leaves from 25 species representing five families of the order Centrospermae (Caryophyllales, Chenopodiales) in the Angiospermae as well as Welwitschia mirabilis (Gymnospermae). Cluster analysis of a 256 spot binary-coded data set derived from the computer-encoded spot patterns of the 25 species successfully separated taxa from the individual to the familial levels of the taxonomic hierarchy in accordance with traditional taxonomic delineations of the taxa.