Population genomics of a symbiont in the early stages of a pest invasion.

Invasive species often depend on microbial symbionts, but few studies have examined the evolutionary dynamics of symbionts during the early stages of an invasion. The insect Megacopta cribraria and its bacterial nutritional symbiont Candidatus Ishikawaella capsulata invaded the southeastern US in 2009. While M. cribraria was initially discovered on wild kudzu plants, it was found as a pest on soybeans within 1 year of infestation. Because prior research suggests Ishikawaella confers the pest status--that is, the ability to thrive on soybeans--in some Megacopta species, we performed a genomic study on Ishikawaella from US. Megacopta cribraria populations to understand the role of the symbiont in driving host plant preferences. We included Ishikawaella samples collected in the first days of the invasion in 2009 and from 23 locations across the insect's 2011 US range. The 0.75 Mb symbiont genome revealed only 47 fixed differences from the pest-conferring Ishikawaella in Japan, with only one amino acid change in a nutrition-provisioning gene. This similarity, along with a lack of fixed substitutions in the US symbiont population, indicates that Ishikawella likely arrived in the US capable of being a soybean pest. Analyses of allele frequency changes between 2009 and 2011 uncover signatures of both positive and negative selection and suggest that symbionts on soybeans and kudzu experience differential selection for genes related to nutrient provisioning. Our data reveal the evolutionary trajectory of an important insect-bacteria symbiosis in the early stages of an invasion, highlighting the role microbial symbionts may play in the spread of invasive species.

Contrasting population genetic patterns within the white-throated sparrow genome (Zonotrichia albicollis).

BACKGROUND: The level of nucleotide diversity observed across the genome is positively correlated with the local rate of recombination. Avian karyotypes are typified by large variation in chromosome size and the rate of recombination in birds has been shown to be negatively correlated with chromosome size. It has thus been predicted that nucleotide diversity is negatively correlated with chromosome size in aves. However, there is limited empirical evidence to support this prediction. RESULTS: Here we sequenced 27 autosomal and 12 sex chromosome-linked loci in the white-throated sparrow (Zonotrichia albicollis) to quantify and compare patterns of recombination, linkage disequilibrium (LD), and genetic diversity across the genome of this North American songbird. Genetic diversity on the autosomes varied up to 8-fold, with the lowest diversity observed on the macrochromosomes and the highest diversity on the microchromosomes. Genetic diversity on the sex chromosomes was reduced compared to the autosomes, the most extreme difference being a ~300-fold difference between the W chromosome and the microchromosomes. LD and population structure associated with a common inversion polymorphism (ZAL2/2m) in this species were found to be atypical compared to other macrochromosomes, and nucleotide diversity within this inversion on the two chromosome arrangements was more similar to that observed on the Z chromosome. CONCLUSIONS: A negative correlation between nucleotide diversity and autosome size was observed in the white-throated sparrow genome, as well as low levels of diversity on the sex chromosomes comparable to those reported in other birds. The population structure and extended LD associated with the ZAL2/2m chromosomal polymorphism are exceptional compared to the rest of the white-throated sparrow genome.

Evolution of a bitter taste receptor gene cluster in a New World sparrow.

Bitter taste perception likely evolved as a protective mechanism against the ingestion of harmful compounds in food. The evolution of the taste receptor type 2 (TAS2R) gene family, which encodes the chemoreceptors that are directly responsible for the detection of bitter compounds, has therefore been of considerable interest. Though TAS2R repertoires have been characterized for a number of species, to date the complement of TAS2Rs from just one bird, the chicken, which had a notably small number of TAS2Rs, has been established. Here, we used targeted mapping and genomic sequencing in the white-throated sparrow (Zonotrichia albicollis) and sample sequencing in other closely related birds to reconstruct the history of a TAS2R gene cluster physically linked to the break points of an evolutionary chromosomal rearrangement. In the white-throated sparrow, this TAS2R cluster encodes up to 18 functional bitter taste receptors and likely underwent a large expansion that predates and/or coincides with the radiation of the Emberizinae subfamily into the New World. In addition to signatures of gene birth-and-death evolution within this cluster, estimates of Ka/Ks for the songbird TAS2Rs were similar to those previously observed in mammals, including humans. Finally, comparison of the complete genomic sequence of the cluster from two common haplotypes in the white-throated sparrow revealed a number of nonsynonymous variants and differences in functional gene content within this species. These results suggest that interspecies and intraspecies genetic variability does exist in avian TAS2Rs and that these differences could contribute to variation in bitter taste perception in birds.

Comparison of two multiple-locus variable-number tandem-repeat analysis methods for molecular strain typing of human Brucella melitensis isolates from the Middle East.

Brucella species are highly monomorphic, with minimal genetic variation among species, hindering the development of reliable subtyping tools for epidemiologic and phylogenetic analyses. Our objective was to compare two distinct multiple-locus variable-number tandem-repeat analysis (MLVA) subtyping methods on a collection of 101 Brucella melitensis isolates from sporadic human cases of brucellosis in Egypt (n = 83), Qatar (n = 17), and Libya (n = 1). A gel-based MLVA technique, MLVA-15(IGM), was compared to an automated capillary electrophoresis-based method, MLVA-15(NAU), with each MLVA scheme examining a unique set of variable-number tandem repeats. Both the MLVA(IGM) and MLVA(NAU) methods were highly discriminatory, resolving 99 and 101 distinct genotypes, respectively, and were able to largely separate genotypes from Egypt and Qatar. The MLVA-15(NAU) scheme presented higher strain-to-strain diversity in our test population than that observed with the MLVA-15(IGM) assay. Both schemes were able to genetically correlate some strains originating from the same hospital or region within a country. In addition to comparing the genotyping abilities of these two schemes, we also compared the usability, limitations, and advantages of the two MLVA systems and their applications in the epidemiological genotyping of human B. melitensis strains.

Global genetic population structure of Bacillus anthracis.

Anthrax, caused by the bacterium Bacillus anthracis, is a disease of historical and current importance that is found throughout the world. The basis of its historical transmission is anecdotal and its true global population structure has remained largely cryptic. Seven diverse B. anthracis strains were whole-genome sequenced to identify rare single nucleotide polymorphisms (SNPs), followed by phylogenetic reconstruction of these characters onto an evolutionary model. This analysis identified SNPs that define the major clonal lineages within the species. These SNPs, in concert with 15 variable number tandem repeat (VNTR) markers, were used to subtype a collection of 1,033 B. anthracis isolates from 42 countries to create an extensive genotype data set. These analyses subdivided the isolates into three previously recognized major lineages (A, B, and C), with further subdivision into 12 clonal sub-lineages or sub-groups and, finally, 221 unique MLVA15 genotypes. This rare genomic variation was used to document the evolutionary progression of B. anthracis and to establish global patterns of diversity. Isolates in the A lineage are widely dispersed globally, whereas the B and C lineages occur on more restricted spatial scales. Molecular clock models based upon genome-wide synonymous substitutions indicate there was a massive radiation of the A lineage that occurred in the mid-Holocene (3,064-6,127 ybp). On more recent temporal scales, the global population structure of B. anthracis reflects colonial-era importation of specific genotypes from the Old World into the New World, as well as the repeated industrial importation of diverse genotypes into developed countries via spore-contaminated animal products. These findings indicate humans have played an important role in the evolution of anthrax by increasing the proliferation and dispersal of this now global disease. Finally, the value of global genotypic analysis for investigating bioterrorist-mediated outbreaks of anthrax is demonstrated.

Tandem repeat regions within the Burkholderia pseudomallei genome and their application for high resolution genotyping.

BACKGROUND: The facultative, intracellular bacterium Burkholderia pseudomallei is the causative agent of melioidosis, a serious infectious disease of humans and animals. We identified and categorized tandem repeat arrays and their distribution throughout the genome of B. pseudomallei strain K96243 in order to develop a genetic typing method for B. pseudomallei. We then screened 104 of the potentially polymorphic loci across a diverse panel of 31 isolates including B. pseudomallei, B. mallei and B. thailandensis in order to identify loci with varying degrees of polymorphism. A subset of these tandem repeat arrays were subsequently developed into a multiple-locus VNTR analysis to examine 66 B. pseudomallei and 21 B. mallei isolates from around the world, as well as 95 lineages from a serial transfer experiment encompassing ~18,000 generations. RESULTS: B. pseudomallei contains a preponderance of tandem repeat loci throughout its genome, many of which are duplicated elsewhere in the genome. The majority of these loci are composed of repeat motif lengths of 6 to 9 bp with 4 to 10 repeat units and are predominately located in intergenic regions of the genome. Across geographically diverse B. pseudomallei and B.mallei isolates, the 32 VNTR loci displayed between 7 and 28 alleles, with Nei's diversity values ranging from 0.47 and 0.94. Mutation rates for these loci are comparable (>10-5 per locus per generation) to that of the most diverse tandemly repeated regions found in other less diverse bacteria. CONCLUSION: The frequency, location and duplicate nature of tandemly repeated regions within the B. pseudomallei genome indicate that these tandem repeat regions may play a role in generating and maintaining adaptive genomic variation. Multiple-locus VNTR analysis revealed extensive diversity within the global isolate set containing B. pseudomallei and B. mallei, and it detected genotypic differences within clonal lineages of both species that were identical using previous typing methods. Given the health threat to humans and livestock and the potential for B. pseudomallei to be released intentionally, MLVA could prove to be an important tool for fine-scale epidemiological or forensic tracking of this increasingly important environmental pathogen.

Strain-specific single-nucleotide polymorphism assays for the Bacillus anthracis Ames strain.

Highly precise diagnostics and forensic assays can be developed through a combination of evolutionary analysis and the exhaustive examination of genomic sequences. In Bacillus anthracis, whole-genome sequencing efforts revealed ca. 3,500 single-nucleotide polymorphisms (SNPs) among eight different strains and evolutionary analysis provides the identification of canonical SNPs. We have previously shown that SNPs are highly evolutionarily stable, and the clonal nature of B. anthracis makes them ideal signatures for subtyping this pathogen. Here we identified SNPs that define the lineage of B. anthracis that contains the Ames strain, the strain used in the 2001 bioterrorist attacks in the United States. Sequencing and real-time PCR were used to validate these SNPs across B. anthracis strains, including (i) 88 globally and genetically diverse isolates; (ii) isolates that were shown to be genetic relatives of the Ames strain by multiple-locus variable number tandem repeat analysis (MLVA); and (iii) several different lab stocks of the Ames strain, including a clinical isolate from the 2001 letter attack. Six SNPs were found to be highly specific for the Ames strain; four on the chromosome, one on the pX01 plasmid, and one on the pX02 plasmid. All six SNPs differentiated the B. anthracis Ames strain from the 88 unique B. anthracis strains, while five of the six separated Ames from its close genetic relatives. The use of these SNPs coupled with real-time PCR allows specific and sensitive (<100 fg of template DNA) identification of the Ames strain. This evolutionary and genomics-based approach provides an effective means for the discovery of strain-specific SNPs in B. anthracis.

Mass spectrometry provides accurate characterization of two genetic marker types in Bacillus anthracis.

Epidemiological and forensic analyses of bioterrorism events involving Bacillus anthracis could be improved if both variable number tandem repeats (VNTRs) and single nucleotide polymorphisms (SNPs) could be combined on a single analysis platform. Here we present the use of electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS) to characterize 24 alleles from 6 VNTR loci and 11 alleles from 7 SNP loci in B. anthracis. The results obtained with ESI-FTICR-MS were consistent with independent results obtained from traditional approaches using electrophoretic detection of fluorescent products. However, ESI-FTICR-MS improves on the traditional approaches because it does not require fluorescent labeling of PCR products, minimizes post-PCR processing, obviates electrophoresis, and provides unambiguous base composition of both SNP and VNTR PCR products. In addition, ESI-FTICR-MS allows both marker types to be examined simultaneously and at a rate of approximately 1 sample per min. This technology represents a significant advance in our ability to rapidly characterize B. anthracis isolates using VNTR and SNP loci.

Phylogenetic discovery bias in Bacillus anthracis using single-nucleotide polymorphisms from whole-genome sequencing.

Phylogenetic reconstruction using molecular data is often subject to homoplasy, leading to inaccurate conclusions about phylogenetic relationships among operational taxonomic units. Compared with other molecular markers, single-nucleotide polymorphisms (SNPs) exhibit extremely low mutation rates, making them rare in recently emerged pathogens, but they are less prone to homoplasy and thus extremely valuable for phylogenetic analyses. Despite their phylogenetic potential, ascertainment bias occurs when SNP characters are discovered through biased taxonomic sampling; by using whole-genome comparisons of five diverse strains of Bacillus anthracis to facilitate SNP discovery, we show that only polymorphisms lying along the evolutionary pathway between reference strains will be observed. We illustrate this in theoretical and simulated data sets in which complex phylogenetic topologies are reduced to linear evolutionary models. Using a set of 990 SNP markers, we also show how divergent branches in our topologies collapse to single points but provide accurate information on internodal distances and points of origin for ancestral clades. These data allowed us to determine the ancestral root of B. anthracis, showing that it lies closer to a newly described "C" branch than to either of two previously described "A" or "B" branches. In addition, subclade rooting of the C branch revealed unequal evolutionary rates that seem to be correlated with ecological parameters and strain attributes. Our use of nonhomoplastic whole-genome SNP characters allows branch points and clade membership to be estimated with great precision, providing greater insight into epidemiological, ecological, and forensic questions.

Anthrax molecular epidemiology and forensics: using the appropriate marker for different evolutionary scales.

Precise identification of Bacillus anthracis isolates has aided forensic and epidemiological analyses of natural anthrax cases, bioterrorism acts and industrial scale accidents by state-sponsored bioweapons programs. Because there is little molecular variation among B. anthracis isolates, identifying and using rare variation is crucial for precise strain identification. We think that mutation is the primary diversifying force in a clonal, recently emerged pathogen, such as B. anthracis, since mutation rate is correlated with diversity on a per locus basis. While single nucleotide polymorphisms (SNPs) are rare, their detection is facilitated by whole genome discovery approaches. As highly stable phylogenetic markers, SNPs are useful for identifying long branches or key phylogenetic positions. Selection of single, diagnostic "Canonical SNPs" (canSNPs) for these phylogenetic positions allows for efficient and defining assays. We have taken a nested hierarchal strategy for subtyping B. anthracis, which is consistent with traditional diagnostics and applicable to a wide range of pathogens. Progressive hierarchical resolving assays using nucleic acids (PHRANA) uses a progression of diagnostic genomic loci that are initially highly stable but with low resolution and, ultimately, very unstable but with high resolution. This approach mitigates the need for data weighting and provides both a deeply rooted phylogenetic hypothesis and high resolution discrimination among closely related isolates.