We report genomes of two species of cactophilic Drosophila: Drosophila arizonae and D. navojoa. These two are the closest relatives of D. mojavensis, forming the D. mojavensis cluster. D. mojavensis and D. arizonae diverged from D. navojoa ∼5.8 Mya, while the split between D. arizonae and D. mojavensis is more recent, at 1.5 Mya. Together the three genomes provide opportunities to examine genomic changes associated with speciation and host shifts in this ecologically defined group of flies. The three species are also separated by fixed inversion differences in three of their six chromosomes. While the levels of nucleotide divergence in the colinear chromosomes are significantly lower than in the inverted chromosomes, consistent with a past role of the inversions in preventing gene flow, the patterns differ among the inverted chromosomes when the locations of nucleotides inside or outside of the inversions are considered. For Muller element E, there is greater divergence external to the inversion breakpoints. For Muller A, the divergence is slightly higher inside the inversions, while for Muller B, the breakpoints and hence the difference in substitutions in relation to the inversions could not be determined. The differences among the inverted chromosomes, especially once the breakpoints are clearly established, could aid in dating the origins of the inversions.

A detailed knowledge of gene function in the monarch butterfly is still lacking. Here we generate a genome assembly from a Mexican nonmigratory population and used RNA-seq data from 14 biological samples for gene annotation and to construct an atlas portraying the breadth of gene expression during most of the monarch life cycle. Two thirds of the genes show expression changes, with long noncoding RNAs being particularly finely regulated during adulthood, and male-biased expression being four times more common than female-biased. The two portions of the monarch heterochromosome Z, one ancestral to the Lepidoptera and the other resulting from a chromosomal fusion, display distinct association with sex-biased expression, reflecting sample-dependent incompleteness or absence of dosage compensation in the ancestral but not the novel portion of the Z. This study presents extended genomic and transcriptomic resources that will facilitate a better understanding of the monarch’s adaptation to a changing environment.

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 (Nef) of about 6 million over the last 10000 years.

Only two parasite interactions are known for Drosophila to date: Allantonematid nematodes associated with mycophagous Drosophilids and the ectoparasitic mite Macrocheles subbadius with the Sonoran Desert endemic Drosophila nigrospiracula. Unlike the nematode-Drosophila association, breadth of mite parasitism on Drosophila species is unknown. As M. subbadius is a generalist, parasitism of additional Drosophilids is expected. We determined the extent and distribution of mite parasitism in nature Drosophilids collected in Mexico and southern California. Thirteen additional species of Drosophilids were infested. Interestingly, 10 belong to the repleta species group of the subgenus Drosophila, despite the fact that the majority of flies collected were of the subgenus Sophophora. In all cases but 2, the associated mites were M. subbadius. Drosophila hexastigma was found to have not only M. subbadius, but another Mesostigmatid mite, Paragarmania bakeri, as well. One D. hydei was also found to have a mite from genus Lasioseius attached. In both choice and no-choice experiments, mites were more attracted to repleta group species than to Sophophoran. The extent of mite parasitism clearly is much broader than previously reported and suggests a host bias mediated either by mite preference and/or some mechanism of resistance in particular Drosophilid lineages.

No abstract