Incomplete sexual isolation in sympatry between subspecies of the butterfly Danaus chrysippus (L.) and the creation of a hybrid zone

Genome assembly of Danaus chrysippus and comparison with the Monarch Danaus plexippus

Milkweed butterflies in the genus Danaus are studied in a diverse range of research fields including the neurobiology of migration, biochemistry of plant detoxification, host–parasite interactions, evolution of sex chromosomes, and speciation. We have assembled a nearly chromosomal genome for Danaus chrysippus (known as the African Monarch, African Queen, and Plain Tiger) using long-read sequencing data. This species is of particular interest for the study of genome structural change and its consequences for evolution. Comparison with the genome of the North American Monarch Danaus plexippus reveals generally strong synteny but highlights 3 inversion differences. The 3 chromosomes involved were previously found to carry peaks of intraspecific differentiation in D. chrysippus in Africa, suggesting that these inversions may be polymorphic and associated with local adaptation. The D. chrysippus genome is over 40% larger than that of D. plexippus, and nearly all of the additional ∼100 Megabases of DNA comprises repeats. Future comparative genomic studies within this genus will shed light on the evolution of genome architecture.

Neo Sex Chromosomes, Colour Polymorphism and Male-Killing in the African Queen Butterfly, Danaus chrysippus (L.)

Danaus chrysippus (L.), one of the world's commonest butterflies, has an extensive range throughout the Old-World tropics. In Africa it is divided into four geographical subspecies which overlap and hybridise freely in the East African Rift: Here alone a male-killing (MK) endosymbiont, Spiroplasma ixodetis, has invaded, causing female-biased populations to predominate. In ssp. chrysippus, inside the Rift only, an autosome carrying a colour locus has fused with the W chromosome to create a neo-W chromosome. A total of 40-100% of Rift females are neo-W and carry Spiroplasma, thus transmitting a linked, matrilineal neo-W, MK complex. As neo-W females have no sons, half the mother's genes are lost in each generation. Paradoxically, although neo-W females have no close male relatives and are thereby forced to outbreed, MK restricts gene flow between subspecies and may thus promote speciation. The neo-W chromosome originated in the Nairobi region around 2.2 k years ago and subsequently spread throughout the Rift contact zone in some 26 k generations, possibly assisted by not having any competing brothers. Our work on the neo-W chromosome, the spread of Spiroplasma and possible speciation is ongoing.

Molecular data suggest population expansion and high level of gene flow in the Plain Tiger (Danaus chrysippus; Nymphalidae: Danainae)

In the present study, we sequenced the individuals of the Danaus chrysippus from the different altitudinal ranges in the western Himalayan state of Uttarakhand, India and compared with other global published data across its distribution range using the cytochrome oxidase c subunit-I (COI). Among the sequences generated in this study, we observed total six haplotypes with intra-species sequence divergence of 0.001-0.009. Whereas the combined data generated a total of 24 haplotypes. The genetic diversity and neutrality test indices suggested overall population expansion of the species. This has also been supported by the MJ Network as it shows the star-like topology and formation of one core haplotype with maximum frequency distribution to the multiple locations. Whereas, the phylogenetic tree shows mixing of the haplotype from the different locations to the same clade. Haplotype arrangement in both phylogenetic tree and MJ Network is not clustered as per their geographic affinity, which suggests incomplete lineage shorting and recent population expansion or colonizations to the new area.

Molecular evidence of hybridization in sympatric populations of the Enantia jethys complex (Lepidoptera: Pieridae)

Hybridization events are frequently demonstrated in natural butterfly populations. One interesting butterfly complex species is the Enantia jethys complex that has been studied for over a century; many debates exist regarding the species composition of this complex. Currently, three species that live sympatrically in the Gulf slope of Mexico (Enantia jethys, E. mazai, and E. albania) are recognized in this complex (based on morphological and molecular studies). Where these species live in sympatry, some cases of interspecific mating have been observed, suggesting hybridization events. Considering this, we employed a multilocus approach (analyses of mitochondrial and nuclear sequences: COI, RpS5, and Wg; and nuclear dominant markers: inter-simple sequence repeat (ISSRs) to study hybridization in sympatric populations from Veracruz, Mexico. Genetic diversity parameters were determined for all molecular markers, and species identification was assessed by different methods such as analyses of molecular variance (AMOVA), clustering, principal coordinate analysis (PCoA), gene flow, and PhiPT parameters. ISSR molecular markers were used for a more profound study of hybridization process. Although species of the Enantia jethys complex have a low dispersal capacity, we observed high genetic diversity, probably reflecting a high density of individuals locally. ISSR markers provided evidence of a contemporary hybridization process, detecting a high number of hybrids (from 17% to 53%) with significant differences in genetic diversity. Furthermore, a directional pattern of hybridization was observed from E. albania to other species. Phylogenetic study through DNA sequencing confirmed the existence of three clades corresponding to the three species previously recognized by morphological and molecular studies. This study underlines the importance of assessing hybridization in evolutionary studies, by tracing the lineage separation process that leads to the origin of new species. Our research demonstrates that hybridization processes have a high occurrence in natural populations.

The evolution of cardenolide-resistant forms of Na⁺,K⁺ -ATPase in Danainae butterflies

Cardenolides are a class of plant secondary compounds that inhibit the proper functioning of the Na(+) , K(+) -ATPase enzyme in susceptible animals. Nonetheless, many insect species are able to sequester cardenolides for their own defence. These include butterflies in the subfamily Danainae (Family: Nymphalidae) such as the monarch (Danaus plexippus). Previous studies demonstrated that monarchs harbour an asparagine (N) to histidine (H) substitution (N122H) in the α subunit of Na(+) , K(+) -ATPase (ATPα) that reduces this enzyme's sensitivity to cardenolides. More recently, it has been suggested that at ATPα position 111, monarchs may also harbour a leucine (L)/glutamine (Q) polymorphism. This later amino acid could also contribute to cardenolide insensitivity. However, here we find that incorrect annotation of the initially reported DNA sequence for ATPα has led to several erroneous conclusions. Using a population genetic and phylogenetic analysis of monarchs and their close relatives, we show that an ancient Q111L substitution occurred prior to the radiation of all Danainae, followed by a second substitution at the same site to valine (V), which arose before the diversification of the Danaus genus. In contrast, N122H appears to be a recent substitution specific to monarchs. Surprisingly, examination of a broader insect phylogeny reveals that the same progression of amino acid substitutions (Q111L → L111V + N122H) has also occurred in Chyrsochus beetles (Family: Chrysomelidae, Subfamily: Eumolpinae) that feed on cardenolide-containing host plants. The parallel pattern of amino acid substitution in these two distantly related lineages is consistent with an adaptive role for these substitutions in reducing cardenolide sensitivity and suggests that their temporal order may be limited by epistatic interactions.

Genetic differentiation and natural hybridization between the Sardinian endemic Maniola nurag and the European Maniola jurtina

The Mediterranean island of Sardinia is known for its multitude of unique genetic lineages. We view one of them in a larger phylogeographic context. The endemic Sardinian Meadow Brown butterfly, Maniola nurag, is restricted to the mountainous areas of the island, whereas its widespread close relative, Maniola jurtina, also occurs on the coast. At intermediate altitudes the species' distributions overlap. There, a number of individuals exhibit phenotypic characteristics intermediate between the two species. We examined patterns of intra- and interpopulation variation in 10 M. nurag populations from Sardinia and 16 M. jurtina populations from Sardinia and continental Europe, as well as 17 intermediate individuals, sampled in 1999-2002, by means of allozyme markers, combining it with a morphometric analysis based on 18 wing-characters of 52 males. At the 15 loci studied (aldolase, aat-1, aat-2, g6pdh, gpd, idh-1, idh-2, mdh-1, mdh-2, mpi, me, leu-ala, pgi, pgm, and 6pgdh), 76 different alleles were detected, 63 of which were shared by M. nurag and M. jurtina. None of the loci was found to be alternatively fixed between the two species. In that respect, this study testifies to the difficulties that may arise when trying to identify hybrids from genotypic data. Levels of genetic variation in island populations (M. jurtina: H(O) = 0.137-0.189; M. nurag: H(O) = 0.141-0.270) were comparable to those of mainland M. jurtina (H(O) = 0.141-0.236). A Bayesian admixture analysis supported the hypothesis of mixed (hybrid) ancestry of individuals occurring at intermediate altitudes. Similarly, neighbour-joining and unweighted pair-group method with arithmetic averaging (UPGMA) analyses, as well as morphometrics hinted at the existence of a Maniola-hybrid zone in Sardinia at intermediate altitudes. We discuss the results in the light of the phylogeography of other Sardinian taxa with the aim to reach a general understanding of the biogeographic history of this island's endemic species.

DISCORDANT DIVERGENCE TIMES AMONG Z-CHROMOSOME REGIONS BETWEEN TWO ECOLOGICALLY DISTINCT SWALLOWTAIL BUTTERFLY SPECIES

We investigate multilocus patterns of differentiation between parental populations of two swallowtail butterfly species that differ at a number of ecologically important sex-linked traits. Using a new coalescent-based approach, we show that there is significant heterogeneity in estimated divergence times among five Z-linked markers, rejecting a purely allopatric speciation model. We infer that the Z chromosome is a mosaic of regions that differ in the extent of historical gene flow, potentially due to isolating barriers that prevent the introgression of species-specific traits that result in hybrid incompatibilities. Surprisingly, a candidate region for a strong barrier to introgression, Ldh, does not show a significantly deeper divergence time than other markers on the Z chromosome. Our approach can be used to test alternative models of speciation and can potentially assign chronological order to the appearance of factors contributing to reproductive isolation between species.