Population Genetics of Overwintering Monarch Butterflies, Danaus plexippus (Linnaeus), from Central Mexico Inferred from Mitochondrial DNA and Microsatellite Markers

Population genetic variation and demographic history in Danaus plexippus (L.), from Mexico were assessed based on analyses of mitochondrial cytochrome c oxidase subunit I (COI; 658 bp) and subunit II (COII; 503 bp) gene segments and 7 microsatellite loci. The sample of 133 individuals included both migratory monarchs, mainly from 4 overwintering sites within the Monarch Butterfly Biosphere Reserve (MBBR) in central Mexico (states of Michoacán and México), and a nonmigratory population from Irapuato, Guanajuato. Haplotype (h) and nucleotide (π) diversities were relatively low, averaging 0.466 and 0.00073, respectively, for COI, and 0.629 and 0.00245 for COII. Analysis of molecular variance of the COI data set, which included additional GenBank sequences from a nonmigratory Costa Rican population, showed significant population structure between Mexican migratory monarchs and nonmigratory monarchs from both Mexico and Costa Rica, suggesting limited gene flow between the 2 behaviorally distinct groups. Interestingly, while the COI haplotype frequencies of the nonmigratory populations differed from the migratory, they were similar to each other, despite the great physical distance between them. Microsatellite analyses, however, suggested a lack of structure between the 2 groups, possibly owing to the number of significant deviations from Hardy–Weinberg equilibrium resulting from heterzoygote deficiencies found for most of the loci. Estimates of demographic history of the combined migratory MBBR monarch population, based on the mismatch distribution and Bayesian skyline analyses of the concatenated COI and COII data set (n = 89) suggested a population expansion dating to the late Pleistocene (~35000–40000 years before present) followed by a stable effective female population size (N(ef)) of about 6 million over the last 10000 years.