Ultrafast Evolution and Loss of CRISPRs Following a Host Shift in a Novel Wildlife Pathogen, Mycoplasma gallisepticum

Measureable rates of genome evolution are well documented in human pathogens but are less well understood in bacterial pathogens in the wild, particularly during and after host switches. Mycoplasma gallisepticum (MG) is a pathogenic bacterium that has evolved predominantly in poultry and recently jumped to wild house finches (Carpodacus mexicanus), a common North American songbird. For the first time we characterize the genome and measure rates of genome evolution in House Finch isolates of MG, as well as in poultry outgroups. Using whole-genome sequences of 12 House Finch isolates across a 13-year serial sample and an additional four newly sequenced poultry strains, we estimate a nucleotide diversity in House Finch isolates of only ∼2% of ancestral poultry strains and a nucleotide substitution rate of 0.8−1.2×10−5 per site per year both in poultry and in House Finches, an exceptionally fast rate rivaling some of the highest estimates reported thus far for bacteria. We also found high diversity and complete turnover of CRISPR arrays in poultry MG strains prior to the switch to the House Finch host, but after the invasion of House Finches there is progressive loss of CRISPR repeat diversity, and recruitment of novel CRISPR repeats ceases. Recent (2007) House Finch MG strains retain only ∼50% of the CRISPR repertoire founding (1994–95) strains and have lost the CRISPR–associated genes required for CRISPR function. Our results suggest that genome evolution in bacterial pathogens of wild birds can be extremely rapid and in this case is accompanied by apparent functional loss of CRISPRs

Stringent_2010_Masked_Val

Most stringent SNP alignmen

Stringent_Moderate_v2_2010_Masked_Val

Working SNP calling alignment

Moderate_2010_Masked_Val

Alignment of SNPs at least stringent SNP calling criteria as described in the text

Evolution of the CRISPR locus in <i>Mycoplasma gallisepticum</i> isolates collected from House Finches, chickens, and turkeys.

<p>Numbers by each strain indicate the number of repeats in each CRISPR array. The ancestral 71-repeat CRISPR array of the chicken MG strain is shown in simplified form at bottom. Diagnostic CRISPR repeats for House Finch MG isolates are indicated in repeat-specific patterns. The black ovals signify the cluster of four CRISP-associated (CAS) genes, which are deleted in the 2007 strains. The tree at left is broadly consistent with the tree based on SNPs (<a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1002511#pgen-1002511-g003" target="_blank">Figure 3</a>) but emphasizes strain clusters indicated by rare genomic changes and CRISPR deletions; it was constructed as described in <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1002511#pgen.1002511.s021" target="_blank">Text S3</a>.</p

MG_Rate_Estimate

One of many BEAST files formed with the working alignment that can produce the data

95% highest posterior density intervals on the estimated substitution rate.

<p>A) for House Finch <i>Mycoplasma</i> strains derived from 34 analyses using the different data and model combinations described in <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1002511#pgen.1002511.s020" target="_blank">Text S2</a>. The middle circle of each bar is the estimated mean; top and bottom circles are the upper and lower 95% bounds of each highest posterior density (HPDs). b) Root-to-tip graph of sampling date of House Finch <i>Mycoplasma</i> strains versus divergence from the closest sequence in the putative source population TK_2001. A simple regression gives an estimated substitution rate of 1.45×10⁻⁵, consistent with estimates from BEAST. See <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1002511#pgen.1002511.s020" target="_blank">Text S2</a> and <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1002511#pgen.1002511.s025" target="_blank">Text S7</a> for further information.</p

Phylogeny of <i>Mycoplasma gallisepticum</i> isolates collected at time points 1994–2007 following a host shift from poultry to House Finches.

<p>The basic topology and branch lengths of the tree come from the output for the BEAST analysis made while estimating evolutionary rates. From this tree we collapsed branches with less than 0.6 posterior probability or if there were no phylogenetically informative SNPs supporting that branch. Several strains are shown as polytomies because their genomic histories are shaped by recombination. Within the House Finch MG clade, branch lengths are proportional to time. Major genomic events are indicated on appropriate branches. The numbers of diagnostic SNPs indicated on various branches are minima. The numbers of CRISPR changes shown are only those that can be constructed with reasonable support (<a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1002511#pgen-1002511-g005" target="_blank">Figure 5</a>); one possible reconstruction is presented.</p