The complete mitochondrial genomes of the Fenton's wood white, Leptidea morsei, and the lemon emigrant, Catopsilia pomona

Complete mitochondrial genome of Colias sifanica Grum-Grshimailo, 1891 (Lepidoptera: Pieridae) from montane region in northwestern China

Colias sifanica Grum-Grshimailo, 1891, is a typical montane butterfly species which occurs on the Qinghai-Tibet Plateau (QTP) and adjacent regions in China. In this study, the complete mitochondrial genome of this species was assembled from data generated by next-generation sequencing. The mitogenome was 15,151 bp in length and comprised 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes and a control region. The base composition of the C. sifanica mitogenome was 39.7% A, 41.3% T, 11.3% C and 7.7% G, significantly AT biased as commonly found in other Pieridae mitogenomes. Phylogenetic analyses based on all PCGs using both the maximum likelihood and Bayesian inference methods indicated that C. sifanica is closely related to C. fieldii with high support values, and the phylogenetic relationship of (Dercas + ((Gandaca + Gonepteryx) + (Phoebis + (Anteos + (Catopsilia + (Zerene + Colias))))))) was shown for the tribe Coliadini. Though both the mitogenomic gene order and overall base composition were found to be conserved, different Ka/Ks ratios for several mitogenomic PCGs were detected between Colias and other species in the tribe Coliadini.

Characterization and Molecular Phylogenetic Analysis of Subfamily Erebinae (Lepidoptera: Noctuoidea: Erebidae) Using Five Complete Mitochondrial Genomes

In this study, the complete mitogenomes of Sympis rufibasis, Lacera noctilio, Oxyodes scrobiculata, Mocis undata, and Artena dotata were newly sequenced to bring up-to-date the database using the next-generation sequencing methods. The gene order of all sequenced mitogenomes was identical consisting of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a non-coding A+T-rich region, which were common to other Lepidopteran insects. All protein-coding genes (PCGs) initiated with a canonical ATN codon and ended with TAN or an incomplete stop codon, single T. The A+T-rich region of S. rufibasis, L. noctilio, O. scrobiculata, M. undata, and A. dotata are 406 bp, 462 bp, 372 bp, 410 bp, and 406 bp long, respectively, containing number of characteristics that are distinctive to Noctuoidea moths. We analyzed concatenated amino acid sequences of protein-coding genes not including rRNAs, using Maximum Likelihood and Bayesian Inference methods. The phylogenetic analyses indicated that the tribe relationships within Erebinae were reconstructed as (Sypnini+((Erebini 1+Poaphilini 1)+((Euclidiini+Catocalini+(Hypopyrini+Erebini 2))+((Hulodini+(Poaphilini 2+Ophiusini))))). Phylogenetic analyses supported and confirmed the monophyly of the subfamilies' relationships as follows: (Hypeninae+Lymantriinae)+((Scoliopterginae+((Calpinae+Erebinae)+((Herminiinae+Aganainae)+Arctiinae)))) within Erebidae.

Variation in the Mitochondrial Genome of the Chagas Disease Vector Triatoma infestans (Hemiptera: Reduviidae)

Chagas' disease is transmitted mainly by members of the subfamily Triatominae (Hemiptera: Reduviidae). Among them, Triatoma infestans (Klug) is the main vector of the disease in Southern Cone of Latin America. In order to contribute to knowledge of the genetic variation between triatomine vectors, in the present study, we analyzed the intraspecific and interspecific variations of the seven mitogenomes available from Triatominae. In addition, in order to examine their evolutionary relationships with others species of Reduviidae and to estimate the divergence time of the main lineages, we constructed phylogenetic trees including mitogenome sequences of 30 species from Reduviidae. Comparative analysis between mitochondrial DNA sequences from two specimens of T. infestans revealed a total of 54 variable sites. Triatoma infestans, Triatoma dimidiata (Latreille), Triatoma rubrofasciata (De Geer), Triatoma migrans (Breddin), Rhodnius pictipes (Stål), and Panstrongylus rufotuberculatus (Champion) present similar mitogenome organization and the length differences observed among these species are primarily caused by variations in control region (CR) and intergenic spacers (IGS). The relative synonymous codon usage values (RSCU) were similar in the six species of Triatominae, and in agreement with the observed in other insects, a biased use of A and C nucleotides in the majority strand was detected. The monophyly of five subfamilies was strongly supported (Phymatinae, Peiratinae, Triatominae, Stenopodainae, and Harpactorinae), while the sampled species of Reduviinae were grouped with one specie from the Salyavatinae subfamily. The oldest subfamily is Phymatinae at 100.3 Mya (99.6-102.2 Mya) and the youngest is Triatominae and Stenopodainae at 52.6 Mya (42.5-63.7 Mya). The estimated diversification time for the Triatominae subfamily agrees with the Andean uplift geological event. An analysis with more mitogenomes from more Triatominae species would be necessary to provide sufficient evidence to support this finding.

Rapid Increase in Genome Size as a Consequence of Transposable Element Hyperactivity in Wood-White (Leptidea) Butterflies

Characterizing and quantifying genome size variation among organisms and understanding if genome size evolves as a consequence of adaptive or stochastic processes have been long-standing goals in evolutionary biology. Here, we investigate genome size variation and association with transposable elements (TEs) across lepidopteran lineages using a novel genome assembly of the common wood-white (Leptidea sinapis) and population re-sequencing data from both L. sinapis and the closely related L. reali and L. juvernica together with 12 previously available lepidopteran genome assemblies. A phylogenetic analysis confirms established relationships among species, but identifies previously unknown intraspecific structure within Leptidea lineages. The genome assembly of L. sinapis is one of the largest of any lepidopteran taxon so far (643 Mb) and genome size is correlated with abundance of TEs, both in Lepidoptera in general and within Leptidea where L. juvernica from Kazakhstan has considerably larger genome size than any other Leptidea population. Specific TE subclasses have been active in different Lepidoptera lineages with a pronounced expansion of predominantly LINEs, DNA elements, and unclassified TEs in the Leptidea lineage after the split from other Pieridae. The rate of genome expansion in Leptidea in general has been in the range of four Mb/Million year (My), with an increase in a particular L. juvernica population to 72 Mb/My. The considerable differences in accumulation rates of specific TE classes in different lineages indicate that TE activity plays a major role in genome size evolution in butterflies and moths.

Speciation in Cloudless Sulphurs Gleaned from Complete Genomes

For 200 years, zoologists have relied on phenotypes to learn about the evolution of animals. A glance at the genotype, even through several gene markers, revolutionized our understanding of animal phylogeny. Recent advances in sequencing techniques allow researchers to study speciation mechanisms and the link between genotype and phenotype using complete genomes. We sequenced and assembled a complete genome of the Cloudless Sulphur (Phoebis sennae) from a single wild-caught specimen. This genome was used as reference to compare genomes of six specimens, three from the eastern populations (Oklahoma and north Texas), referred to as a subspeciesPhoebis sennae eubule, and three from the southwestern populations (south Texas) known as a subspeciesPhoebis sennae marcellina While the two subspecies differ only subtly in phenotype and mitochondrial DNA, comparison of their complete genomes revealed consistent and significant differences, which are more prominent than those between tiger swallowtailsPterourus canadensisandPterourus glaucus The two sulphur taxa differed in histone methylation regulators, chromatin-associated proteins, circadian clock, and early development proteins. Despite being well separated on the whole-genome level, the two taxa show introgression, with gene flow mainly fromP. s. marcellinatoP. s. eubule Functional analysis of introgressed genes reveals enrichment in transmembrane transporters. Many transporters are responsible for nutrient uptake, and their introgression may be of selective advantage for caterpillars to feed on more diverse food resources. Phylogenetically, complete genomes place family Pieridae away from Papilionidae, which is consistent with previous analyses based on several gene markers.

The complete mitochondrial genome of Papilio glaucus and its phylogenetic implications

Due to the intriguing morphology, lifecycle, and diversity of butterflies and moths, Lepidoptera are emerging as model organisms for the study of genetics, evolution and speciation. The progress of these studies relies on decoding Lepidoptera genomes, both nuclear and mitochondrial. Here we describe a protocol to obtain mitogenomes from Next Generation Sequencing reads performed for whole-genome sequencing and report the complete mitogenome of Papilio (Pterourus) glaucus. The circular mitogenome is 15,306 bp in length and rich in A and T. It contains 13 protein-coding genes (PCGs), 22 transfer-RNA-coding genes (tRNA), and 2 ribosomal-RNA-coding genes (rRNA), with a gene order typical for mitogenomes of Lepidoptera. We performed phylogenetic analyses based on PCG and RNA-coding genes or protein sequences using Bayesian Inference and Maximum Likelihood methods. The phylogenetic trees consistently show that among species with available mitogenomes Papilio glaucus is the closest to Papilio (Agehana) maraho from Asia.