Fifty-percent majority-rule consensus phylogenetic tree of <i>Creatonotos</i> spp. recovered from Bayesian inference analysis of an alignment comprising 61 <i>COI</i> haplotypes of <i>Creatonotos</i> spp. and two haplotypes of the out-group taxa (<i>Arctia menetriesii</i> and <i>A</i>. <i>tundrana</i>).

<p>The island localities are in bold. Black numbers near branches are Bayesian posterior probabilities (BPP). The red number near each primary clade is the probability of each species-level MOTU based on the highest Bayesian supported solution of the PTP model.</p

DNA barcoding unravels contrasting evolutionary history of two widespread Asian tiger moth species during the Late Pleistocene

<div><p>Populations of widespread pest insects in tropical areas are characterized by a complex evolutionary history, with overlapping natural and human-mediated dispersal events, sudden expansions, and bottlenecks. Here, we provide biogeographic reconstructions for two widespread pest species in the tiger moth genus <i>Creatonotos</i> (Lepidoptera: Erebidae: Arctiinae) based on the mitochondrial <i>cytochrome c oxidase subunit I</i> (<i>COI</i>) gene. The Asian <i>Creatonotos transiens</i> reveals shallow genetic divergence between distant populations that does not support its current intraspecific systematics with several local subspecies. In contrast, the more widespread <i>Creatonotos gangis</i> comprises at least three divergent subclades corresponding to certain geographic areas, i.e. Australia, Arabia + South Asia and Southeast Asia. With respect to our approximate Bayesian computation (ABC) model, the expansion of <i>Creatonotos gangis</i> into Australia is placed in the Late Pleistocene (~65–63 ka). This dating coincide with an approximate time of the earliest human migration into the continent (~65–54 ka) and the period of intervisibility between Timor and Australia (~65–62 ka). Our findings highlight that the drying Sunda and Sahul shelf areas likely support successful migrations of Asian taxa into Australia during the Pleistocene. The phylogeographic patterns discovered in this study can be used to improve the effectiveness of integrated pest control programs that is a task of substantial practical importance to a broad range of agricultural stakeholders.</p></div

DNA barcoding unravels contrasting evolutionary history of two widespread Asian tiger moth species during the Late Pleistocene - Fig 2

<p>Phylogeography of <i>Creatonotos</i> spp. (<b>A</b>) Median-joining network of <i>COI</i> sequences (see <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0194200#pone.0194200.s003" target="_blank">S1 Table</a> for details). Photos (male specimens): <i>C</i>. <i>gangis</i> [Indonesia, Flores Island, voucher no. Sph0595] and <i>C</i>. <i>transiens</i> [Thailand, near Tham Lod Cave, voucher no. Sph0624] by Vitaly M. Spitsyn and <i>C</i>. <i>leucanioides</i> [Tanzania] by Roy Goff (with his permission; <a href="https://www.africanmoths.com" target="_blank">www.africanmoths.com</a>). (<b>B</b>) Map of approximate collection localities of the specimens in accordance with the respective countries (color squares). Small squares indicate island sites. The map was created using ESRI ArcGIS 10 software (<a href="https://www.esri.com/arcgis" target="_blank">www.esri.com/arcgis</a>); the base of the map was created with ESRI Data and Maps. (<b>C</b>) Frequency histogram of the barcoding gap analysis.</p

Simplified summary of expansion routes inferred across populations of <i>Creatonotos gangis</i>, and examples of male specimens and habitat of Southeast Asian population.

<p>(<b>A</b>) Map of expansion routes. Color circles indicate an approximate range of each population: Southeast Asian population (blue), Arabian–South Asian population (green), and Australian population (orange). Red arrows indicate the directions of expansion routes. The numbers near arrows show the mean age of putative expansion events (in thousands of years) obtained from the ABC model. The map was created using ESRI ArcGIS 10 software (<a href="https://www.esri.com/arcgis" target="_blank">www.esri.com/arcgis</a>); the base of the map was created with ESRI Data and Maps. (<b>B</b>) Male specimen, Maehongson, Thailand. (<b>C</b>) Male specimen, Kachin, Myanmar. (<b>D</b>) Paddy field, a typical habitat of the species, Thanh Hoa Province, Vietnam. (Photos: Vitaly M. Spitsyn).</p

Mismatch distributions of <i>Creatonotos</i> spp. samples based on the mitochondrial <i>COI</i> gene.

<p>Solid black lines indicate observed distribution, and solid red lines represent simulated distribution under a spatial expansion model. Dashed lines represent lower and upper confidence intervals (<i>p</i> = 0.01). (<b>A</b>) <i>C</i>. <i>gangis</i>, Eurasia (<i>N</i> = 40 sequences; Raggedness <i>P</i> = 0.342; Model (SSD) <i>P</i> = 0.317). (<b>B</b>) <i>C</i>. <i>gangis</i>, Arabia and South Asia (<i>N</i> = 20 sequences; Raggedness <i>P</i> = 0.264; Model (SSD) <i>P</i> = 0.273). (<b>C</b>) <i>C</i>. <i>gangis</i>, mainland Southeast Asia (<i>N</i> = 20 sequences; Raggedness <i>P</i> = 0.089; Model (SSD) <i>P</i> = 0.331). (<b>D</b>) <i>C</i>. <i>gangis</i>, Australia (<i>N</i> = 15 sequences; Raggedness <i>P</i> = 0.910; Model (SSD) <i>P</i> = 0.736). (<b>E</b>) <i>C</i>. <i>gangis</i>, the entire range (<i>N</i> = 57 sequences; Raggedness <i>P</i> = 0.715; Model (SSD) <i>P</i> = 0.557). (<b>F</b>) <i>C</i>. <i>transiens</i>, Eurasia (<i>N</i> = 39 sequences; Raggedness <i>P</i> = 0.671; Model (SSD) <i>P</i> = 0.596).</p

ABC modeling of origin of the <i>Creatonotos gangis</i> populations.

<p>(<b>A</b>) Biogeographic scenarios that were tested under an ABC framework using <i>COI</i> gene sequences. Southeast Asian population: samples from Myanmar, Vietnam, Thailand, and South China, Arabian–South Asian population: samples from Oman, Pakistan, India, and Nepal, and Australian population. Effective population size: <i>N</i>₁ –Southeast Asian population; <i>N</i>₂ –Arabian–South Asian population; <i>N</i>₃ –Australian population; <i>N</i>_2b and <i>N</i>_3b –hypothetical founder population for Arabian–South Asian and Australian populations, respectively. Time intervals: <i>t</i>_1a and <i>t</i>₁ –time of the primary split between populations <i>Pop1</i> and <i>Pop2</i> under scenarios 1-CG (before the mid-Pleistocene) and 2-CG (since the mid-Pleistocene), respectively; <i>t</i>_2a and <i>t</i>₂ –time of split between Southeast Asian and Australian populations under scenarios 1-CG (before the Late Pleistocene) and 2-CG (during the Late Pleistocene), respectively; <i>t</i>-<i>db</i>–period of low effective population size <i>N</i>_2b and <i>N</i>_3b since colonization of the Arabian–South Asian Region and Australia. Prior settings are presented in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0194200#pone.0194200.s007" target="_blank">S5 Table</a>. Evaluating the confidence in scenario choice using the direct (<b>B</b>) and linear regression (<b>C</b>) approaches.</p

Fifty-percent majority-rule consensus phylogenetic tree of <i>Creatonotos</i> spp. recovered from Bayesian inference analysis of an alignment comprising 61 <i>COI</i> haplotypes of <i>Creatonotos</i> spp. and two haplotypes of the out-group taxa (<i>Arctia menetriesii</i> and <i>A</i>. <i>tundrana</i>).

<p>The island localities are in bold. Black numbers near branches are Bayesian posterior probabilities (BPP). The red number near each primary clade is the probability of each species-level MOTU based on the highest Bayesian supported solution of the PTP model.</p

Shells of <i>Margaritifera dahurica</i> (Middendorff, 1850).

<p>A—Lectotype (ZISP: no. 7a). B—specimen from Komarovka River, Razdolnaya River basin, Primorye. Photos by I.V. Vikhrev. Scale bar—2 cm.</p

General plan of soft tissue of three Far Eastern margaritiferid species.

<p>A—<i>Margaritifera dahurica</i> (Middendorff, 1850). B—<i>M</i>. <i>laevis</i> (Haas, 1910). C—<i>M</i>. <i>middendorffi</i> (Rosén, 1926). Scale bar—1 cm.</p

Typical habitats of the Eastern Asian margaritiferid species.

<p>A—<i>Margaritifera dahurica</i> (Middendorff, 1850): the Ilystaya River, Primorsky Kray. B—<i>M</i>. <i>middendorffi</i> (Rosén, 1926): the Nachilova River, Kamchatka. C—<i>M</i>. <i>middendorffi</i> (Rosén, 1926) & <i>M</i>. <i>laevis</i> (Haas, 1910): the Golovnina River, Kunashir Island. Photos by Y.V. Bespalaya, Y.S. Kolosova & I.V. Vikhrev.</p