The utility of the neglected mitochondrial control region for evolutionary studies in lepidoptera (insecta)

Mitotyping of random bred cats and pure breed cats (Turkish Angora and Turkish Van) using non-repetitive mitochondrial DNA control region

The Fertile Crescent appears to be the most plausible region where the domestication of cats commenced through a mutually beneficial relationship between wild cats and early agrarian societies. These domesticated cats then journeyed across the globe mirroring the paths of human migration. An examination of mitochondrial DNA (mtDNA) control region-based mitotyping suggested that a significant majority, exceeding 80%, of globally sampled random-bred and pure-bred cats could be categorized into 12 predominant mitotypes. However, the extent of mitotype diversity within random-bred cats from regions proximate to the Fertile Crescent remains inadequately explored. In light of this we aimed to investigate the mitotype diversity in random bred cats sampled from various regions across Turkey. Additionally, we sought to establish a comparison with the mitotype profiles of locally recognized pure breeds, namely the Turkish Angora and Turkish Van. To unravel their evolutionary narratives, we engaged in comprehensive population genetics analyses at both the individual and mitotype-based levels. Our study encompassed a sample size of 240 specimens, forming the basis for both mitotyping and population genetics scrutiny. Our analysis yielded the identification of nine 'universal' mitotypes (A-J), alongside an 'outlier' mitotype group I. Notably mitotypes A and D emerged as particularly prevalent in contrast to the lesser occurrence mitotypes C, G, and H. With the realm of random bred cats the structure of haplotypes exhibited remarkable diversity presenting distinctions from Turkish Angora and Van breeds. Nucleotide diversity was higher compared to previous reports from Turkey and was one of the highest among reported world cat population estimates. Intriguingly, our investigations did not unveil any pronounced instances of strong selection, population expansions or contractions within any specific population or mitotype. To conclude, our study represents a pioneering effort in uncovering the mitotype profiles and haplotype structures inherent to both random-bred and pure breed cats in Turkey. This endeavor not only broadens our understanding of the feline genetic landscape within the region but also lays the foundation for future inquiries into the evolutionary trajectories and genetic legacies of these feline populations.

pgHMA: Application of the heteroduplex mobility assay analysis in phylogenetics and population genetics

The heteroduplex mobility assay (HMA) has proven to be a robust tool for the detection of genetic variation. Here, we describe a simple and rapid application of the HMA by microfluidic capillary electrophoresis, for phylogenetics and population genetic analyses (pgHMA). We show how commonly applied techniques in phylogenetics and population genetics have equivalents with pgHMA: phylogenetic reconstruction with bootstrapping, skyline plots, and mismatch distribution analysis. We assess the performance and accuracy of pgHMA by comparing the results obtained against those obtained using standard methods of analyses applied to sequencing data. The resulting comparisons demonstrate that: (a) there is a significant linear relationship (R²  = .992) between heteroduplex mobility and genetic distance, (b) phylogenetic trees obtained by HMA and nucleotide sequences present nearly identical topologies, (c) clades with high pgHMA parametric bootstrap support also have high bootstrap support on nucleotide phylogenies, (d) skyline plots estimated from the UPGMA trees of HMA and Bayesian trees of nucleotide data reveal similar trends, especially for the median trend estimate of effective population size, and (e) optimized mismatch distributions of HMA are closely fitted to the mismatch distributions of nucleotide sequences. In summary, pgHMA is an easily-applied method for approximating phylogenetic diversity and population trends.

Mitochondrial genome sequencing and phylogeny of Haemagogus albomaculatus, Haemagogus leucocelaenus, Haemagogus spegazzinii, and Haemagogus tropicalis (Diptera: Culicidae)

The genus Haemagogus (Diptera: Culicidae) comprises species of great epidemiological relevance, involved in transmission cycles of the Yellow fever virus and other arboviruses in South America. So far, only Haemagogus janthinomys has complete mitochondrial sequences available. Given the unavailability of information related to aspects of the evolutionary biology and molecular taxonomy of this genus, we report here, the first sequencing of the mitogenomes of Haemagogus albomaculatus, Haemagogus leucocelaenus, Haemagogus spegazzinii, and Haemagogus tropicalis. The mitogenomes showed an average length of 15,038 bp, average AT content of 79.3%, positive AT-skews, negative GC-skews, and comprised 37 functional subunits (13 PCGs, 22 tRNA, and 02 rRNA). The PCGs showed ATN as start codon, TAA as stop codon, and signs of purifying selection. The tRNAs had the typical leaf clover structure, except tRNA^Ser1. Phylogenetic analyzes of Bayesian inference and Maximum Likelihood, based on concatenated sequences from all 13 PCGs, produced identical topologies and strongly supported the monophyletic relationship between the Haemagogus and Conopostegus subgenera, and corroborated with the known taxonomic classification of the evaluated taxa, based on external morphological aspects. The information produced on the mitogenomes of the Haemagogus species evaluated here may be useful in carrying out future taxonomic and evolutionary studies of the genus.

First mitochondrial genome from Yponomeutidae (Lepidoptera, Yponomeutoidea) and the phylogenetic analysis for Lepidoptera

The complete mitochondrial genome (mitogenome) of Yponomeuta montanatus is sequenced and compared with other published yponomeutoid mitogenomes. The mitogenome is circular, 15,349 bp long, and includes the typical metazoan mitochondrial genes (13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes) and an A + T-rich region. All 13 protein-coding genes use a typical start codon ATN, the one exception being cox1, which uses CGA across yponomeutoid mitogenomes. Comparative analyses further show that the secondary structures of tRNAs are conserved, including loss of the Dihydorouidine (DHU) arm in trnS1 (AGN), but remarkable nucleotide variation has occurred mainly in the DHU arms and pseudouridine (TψC) loops. A + T-rich regions exhibit substantial length variation among yponomeutoid mitogenomes, and conserved sequence blocks are recognized but some of them are not present in all species. Multiple phylogenetic analyses confirm the position of Y. montanatus in Yponomeutoidea. However, the superfamily-level relationships in the Macroheterocera clade in Lepidoptera recovered herein show considerable difference with that recovered in previous mitogenomic studies, raising the necessity of extensive phylogenetic investigation when more mitogenomes become available for this clade.

The complete mitochondrial genome of a medicinal insect, Hydrillodes repugnalis (Lepidoptera: Noctuoidea: Erebidae), and related phylogenetic analysis

The complete mitochondrial genome (mitogenome) of an important medicinal insect Hydrillodes repugnalis (Lepidoptera: Noctuoidea) was sequenced and analyzed. The mitogenome is circular with 15,570 bp long, and shows typical gene content and arrangement. Nucleotide composition is highly biased toward A + T nucleotides (81.1%). All protein-coding genes (PCGs) initiate with canonical start codon ATN, except for cox1 being CGA. The typical stop codon TAA is used for most PCGs, while the nad4l uses the TAG, and cox1 and cox2 use incomplete termination codon T. All tRNAs have a typical clover-leaf structure, except for trnS1 (AGN) lacking the dihydrouridine arm. Comparative mitogenome analysis showed that the motif "ATGATAA" between atp8 and atp6, and the motif "ATACTAA" between trnS2 and nad1 are commonly present in noctuoid mitogenomes. In A + T-rich regions, the motif "ATAGA" and subsequent poly-T structure, the motif "ATTTA" and followed macrosatellite (AT)_n element and an "A"-rich 3' end upstream of the trnM gene can be recognized across noctuoid mitogenomes. Phylogenetic analyses showed that H. repugnalis is nested into the Erebidae clade, consistently being sister to the Aganainae. Within Noctuoidea, the (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))) was consistently recovered firstly based on multiple mitochondrial datasets.

Complete mitochondrial genome of the soybean leaffolder, Omiodes indicata (Lepidoptera: Pyraloidea: Crambidae), and phylogenetic analysis for Pyraloidea

Complete mitochondrial genome (mitogenome) of the Omiodes indicata was sequenced and characterized. The circular mitogenome is 15,367bp long, including 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A+T-rich region. Nucleotide composition is highly biased toward A+T nucleotides (81.6%). All 13 PCGs initiate with canonical start codon (ATN), except for cox1 that initiates with CGA. All tRNAs have a typical clover-leaf structure, except for trnS1 (AGN) in which the base pairs of the dihydrouridine (DHU) arm are reduced. In O. indicata, the motifs "ATGATAA" and "ATACTAA" between atp8 and atp6, trnS2 and nad1, respectively, and the motifs "ATAG" and "ATTTA" in the A+T-rich region can be identified. Comparative phylogenetic analyses based on four datasets show that the dataset including all coding positions of 13 PCGs exhibit the highest informativeness in resolving higher phylogeny of Pyraloidea. Bayesian inference (BI) and maximum likelihood (ML) analyses yield generally well-supported phylogenetic relationships among the eleven pyraloid subfamilies involved. However, the relationships among the five subfamilies (Acentropinae, Crambinae, Glaphyriinae, Schoenobiinae and Scopariinae) in ML analysis are ambiguous, which might be resolved by ample sampling in future mitogenome-based phylogenetic studies.

Mitochondrial genome sequencing reveals potential origins of the scabies mite Sarcoptes scabiei infesting two iconic Australian marsupials

BACKGROUND

Debilitating skin infestations caused by the mite, Sarcoptes scabiei, have a profound impact on human and animal health globally. In Australia, this impact is evident across different segments of Australian society, with a growing recognition that it can contribute to rapid declines of native Australian marsupials. Cross-host transmission has been suggested to play a significant role in the epidemiology and origin of mite infestations in different species but a chronic lack of genetic resources has made further inferences difficult. To investigate the origins and molecular epidemiology of S. scabiei in Australian wildlife, we sequenced the mitochondrial genomes of S. scabiei from diseased wombats (Vombatus ursinus) and koalas (Phascolarctos cinereus) spanning New South Wales, Victoria and Tasmania, and compared them with the recently sequenced mitochondrial genome sequences of S. scabiei from humans.

RESULTS

We found unique S. scabiei haplotypes among individual wombat and koala hosts with high sequence similarity (99.1% - 100%). Phylogenetic analysis of near full-length mitochondrial genomes revealed three clades of S. scabiei (one human and two marsupial), with no apparent geographic or host species pattern, suggestive of multiple introductions. The availability of additional mitochondrial gene sequences also enabled a re-evaluation of a range of putative molecular markers of S. scabiei, revealing that cox1 is the most informative gene for molecular epidemiological investigations. Utilising this gene target, we provide additional evidence to support cross-host transmission between different animal hosts.

CONCLUSIONS

Our results suggest a history of parasite invasion through colonisation of Australia from hosts across the globe and the potential for cross-host transmission being a common feature of the epidemiology of this neglected pathogen. If this is the case, comparable patterns may exist elsewhere in the 'New World'. This work provides a basis for expanded molecular studies into mange epidemiology in humans and animals in Australia and other geographic regions.

Transcriptome sequencing and phylogenetic analysis of four species of luminescent beetles

The evolution of bioluminescence has prompted scientific attention to illuminate phylogenetic relationships of luminescent beetles. However, genomic resources are virtually lacking in rhagophthalmids (Rhagophthalmidae) and their related firefly beetles lampyrids (Lampyridae). Here, we employed the Illumina Hiseq 2000 platform and sequenced the whole-body transcriptomes of the four luminescent beetles: one rhagophthalmid (Rhagophthalmus sp.) and three fireflies (Asymmetricata circumdata, Aquatica ficta, and Pyrocoelia pectoralis). We obtained 55.4, 43.4, 38.6, and 36.7 million clean reads for the four species, respectively. All reads were assembled into contigs from which unigenes were derived. All unigenes were annotated by publicly available databases, and a total of 4325 orthologous genes were identified. Using multiple phylogenetic approaches, our transcriptome data confirmed the distinctiveness of Rhagophthalmidae from Lampyridae, which was also supported by our mitogenome analysis using three newly determined mitogenome sequences and 12 previously published ones. Together, this study is the first report of whole transcriptome sequencing data in Rhagophthalmidae and Lampyridae species, representing a valuable genomic resource for studying the origin and evolution of some remarkable traits in these beetles such as bioluminescence. Moreover, our transcriptome and mitogenome data provide useful phylogenetic information that could be of importance in future studies of phylogenetic inference.

Extra-Mediterranean glacial refugia in a Mediterranean faunal element: the phylogeography of the chalk-hill blue Polyommatus coridon (Lepidoptera, Lycaenidae)

Most warm-adapted Central European species are thought to have survived ice ages exclusively in Mediterranean refugia. During recent years, this point of view has been questioned. Therefore, we tested the hypothesis that extra-Mediterranean refugia also played a role in warm-adapted insect species and selected the chalk-hill blue, Polyommatus coridon. We sequenced two mitochondrial loci (COI, CR) in 150 individuals from 30 populations covering nearly the complete range. Minimum spanning networks and other statistical analyses concordantly revealed four genetic lineages with strong phylogeographic signal: a western group in Italy, France and western/central Germany, an eastern lineage in the Balkan Peninsula, the Carpathian Basin and eastern Central Europe, an Alpine group with populations in the Alps and southern Germany and a Pyrenean group. Our results are generally consistent with previous analyses for P. coridon based on allozymes and DNA sequences, but provide additional insights. We propose that these four lineages have evolved during allopatry in different glacial refugia, two in typical Mediterranean refugia (Apennines and Balkan Peninsulas), but two in extra-Mediterranean areas south of the Alps and Pyrenees. This supports survival of warm-adapted organisms in these regions in close geographic proximity to the refugia of high mountain species.

Chromosomal and mitochondrial diversity in Melitaea didyma complex (Lepidoptera, Nymphalidae): eleven deeply diverged DNA barcode groups in one non-monophyletic species?

It is generally accepted that cases of species' polyphyly in COI trees arising as a result of deep intraspecific divergence are negligible, and the detected cases reflect misidentifications or/and methodological errors. Here we studied the problem of species' non-monophyly through chromosomal and molecular analysis of butterfly taxa close to Melitaea didyma (Esper, 1779) (Lepidoptera, Nymphalidae). We found absence or low interspecific chromosome number variation and presence of intraspecific variation, therefore we conclude that in this group, chromosome numbers have relatively low value as taxonomic markers. Despite low karyotype variability, the group was found to have unexpectedly high mitochondrial haplotype diversity. These haplotypes were clustered in 23 highly diverged haplogroups. Twelve of these haplogroups are associated with nine traditionally recognized and morphologically distinct species Melitaea chitralensis Moore, 1901, Melitaea deserticola Oberthür, 1909, Melitaea didymoides Eversmann, 1847, Melitaea gina Higgins, 1941, Melitaea interrupta Colenati, 1846, Melitaea latonigena Eversmann, 1847, Melitaea mixta Evans, 1912, Melitaea saxatilis Christoph, 1873 and Melitaea sutschana Staudinger, 1892. The rest of the haplogroups (11 lineages) belong to a well-known west-palaearctic species Melitaea didyma. The last species is particularly unusual in the haplotypes we obtained. First, it is clearly polyphyletic with respect to COI gene. Second, the differentiation in COI gene between these mostly allopatric (but in few cases sympatric) eleven lineages is extremely high (up to 7.4%), i.e. much deeper than the "standard" DNA barcode species threshold (2.7-3%). This level of divergence normally could correspond not even to different species, but to different genera. Despite this divergence, the bearers of these haplogroups were found to be morphologically indistinguishable and, most importantly, to share absolutely the same ecological niches, i.e. demonstrating the pattern which is hardly compatible with hypothesis of multiple cryptic species. Most likely such a profound irregularity in barcodes is caused by reasons other than speciation and represents an extraordinary example of intra-species barcode variability. Given the deep level of genetic differentiation between the lineages, we assume that there was a long period (up to 5.0 My) of allopatric differentiation when the lineages were separated by geographic or/and ecological barriers and evolved in late Pliocene and Pleistocene refugia of north Africa, the Iberian and Balkan Peninsulas, the Middle East and Central Asia. We discuss the refugia-within-refugia concept as a mechanism explaining the presence of additional diverged minor haplogroups within the areas of the major haplogroups. We also provide the first record of Melitaea gina in Azerbaijan and the record of Melitaea didyma turkestanica as a new taxon for Russia and Europe.

Phylogeography of the Spanish Moon Moth Graellsia isabellae (Lepidoptera, Saturniidae)

BACKGROUND

Geographic and demographic factors as well as specialisation to a new host-plant may lead to host-associated differentiation in plant-feeding insects. We explored the phylogeography of a protected moth, Graellsia isabellae, and its two recognised host-plant species (Pinus sylvestris and P. nigra) in order to seek for any concordance useful to disentangle the evolutionary history of this iconic lepidopteran.

RESULTS

DNA variation in one mitochondrial marker and nine nuclear microsatellite loci revealed a strong phylogeographic pattern across 28 populations of G. isabellae studied in Spain and France comprising six groups mostly distributed along different mountain ranges. Reanalysis of a previously published chloroplast microsatellite dataset revealed a three and two-group structure for Spanish P. sylvestris and P. nigra, respectively. Overall, the population groupings of this protected moth did not match the ones of P. sylvestris and P. nigra.

CONCLUSIONS

There was no evidence of host-associated differentiation between populations using P. sylvestris and the ones inhabiting P. nigra. The two major mitochondrial clades of G. isabellae likely diverged before the Last Glacial Maximum and geographically separated the species into a "southern" (Central and Southern Iberian clusters) and a "northern" lineage (Eastern Iberian, Pyrenean and French Alpine clusters). The Eastern Iberian System, where this insect uses both host-plants, harboured the highest level of genetic diversity. Such a group independently colonised the West and East parts of the Pyrenees. Our results point to a native origin for the French populations occurring in the Alps, genetically related to the Eastern Iberian and Pyrenean sites. The Central Iberian group derived from Southern Iberian ancestors. Secondary contacts were inferred between the Southern/Central Iberian populations and Eastern Iberian cluster as well as between the two Pyrenean ones. The mito-nuclear discordance observed with regard to the Eastern Iberian cluster is congruent with a secondary contact after the evolution of mito-nuclear incompatibilities in geographically isolated areas.

DNA barcoding reveals twelve lineages with properties of phylogenetic and biological species within Melitaea didyma sensu lato (Lepidoptera, Nymphalidae)

The complex of butterfly taxa close to Melitaea didyma includes the traditionally recognized species Melitaea didyma, Melitaea didymoides and Melitaea sutschana, the taxa that were recognized as species only relatively recently (Melitaea latonigena, Melitaea interrupta, Melitaea chitralensis and Melitaea mixta) as well as numerous described subspecies and forms with unclear taxonomic status. Here analysis of mitochondrial DNA barcodes is used to demonstrate that this complex is monophyletic group consisting of at least 12 major haplogroups strongly differentiated with respect to the gene COI. Six of these haplogroups are shown to correspond to six of the above-mentioned species (Melitaea didymoides, Melitaea sutschana, Melitaea latonigena, Melitaea interrupta, Melitaea chitralensis and Melitaea mixta). It is hypothesized that each of the remaining six haplogroups also represents a distinct species (Melitaea mauretanica, Melitaea occidentalis, Melitaea didyma, Melitaea neera, Melitaea liliputana and Melitaea turkestanica), since merging these haplogroups would result in a polyphyletic assemblage and the genetic distances between them are comparable with those found between the other six previously recognized species.

Three in One--Multiple Faunal Elements within an Endangered European Butterfly Species

Ice ages within Europe forced many species to retreat to refugia, of which three major biogeographic basic types can be distinguished: "Mediterranean", "Continental" and "Alpine / Arctic" species. However, this classification often fails to explain the complex phylogeography of European species with a wide range of latitudinal and altitudinal distribution. Hence, we tested for the possibility that all three mentioned faunal elements are represented within one species. Our data was obtained by scoring 1,307 Euphydryas aurinia individuals (46 European locations) for 17 allozyme loci, and sequencing a subset of 492 individuals (21 sites) for a 626 base pairs COI fragment. Genetic diversity indices, F statistics, hierarchical analyses of molecular variance, individual-based clustering, and networks were used to explore the phylogeographic patterns. The COI fragment represented 18 haplotypes showing a strong geographic structure. All but one allozyme loci analysed were polymorphic with a mean F_ST of 0.20, supporting a pronounced among population structure. Interpretation of both genetic marker systems, using several analytical tools, calls for the recognition of twelve genetic groups. These analyses consistently distinguished different groups in Iberia (2), Italy, Provence, Alps (3), Slovenia, Carpathian Basin, the lowlands of West and Central Europe as well as Estonia, often with considerable additional substructures. The genetic data strongly support the hypothesis that E. aurinia survived the last glaciation in Mediterranean, extra-Mediterranean and perialpine refugia. It is thus a rare example of a model organism that combines attributes of faunal elements from all three of these sources. The observed differences between allozymes and mtDNA most likely result from recent introgression of mtDNA into nuclear allozyme groups. Our results indicate discrepancies with the morphologically-based subspecies models, underlining the need to revise the current taxonomy.

Sequence and analysis of mtDNA control region in Cnaphalacrocis medinalis Guenée and Marasmia patnalis Bradley (Lepidoptera: Pyralidae)

Rice leaffolders Cnaphalocrocis medinalis and Marasmia patnalis (Lepidoptera: Pyralidae) are important rice pests and morphologically related with same damage methods. The control region of mitochondrial DNA (mtDNA) in two rice leaffolders C. medinalis and M. patnalis were sequenced and analyzed. The control regions of the two rice leaffolders are located between 12S ribosomal RNA and tRNA_Met, and exhibit conserved structural elements. We identified the poly T stretches, (TA)n block, and stem-loop structure. There are no long tandem repeats found in mtDNA control region in the two rice leaffolders studied. Nevertheless, we did not find the GA-rich block in mtDNA control regions of the two rice leaffolders. The molecular-based phylogenies support the traditional morphologically based view of relationships of Pyralidae within the Ditrysia. The addition of the C. medinalis and M. patnalis mtDNA control regions to the literature promotes the understanding of the molecular evolution of Pyralidae in Ditrysia.

Genetic identification of endangered North African ungulates using noninvasive sampling

North African ungulates include several threatened and emblematic species, yet are poorly studied mainly due to their remoteness and elusiveness. Noninvasive sampling provides a useful approach to obtain ecological and genetic information essential to guide conservation actions. The very first and most important step in conservation planning is to accurately identify species, and molecular genetics has been proved to be a useful tool. Several molecular genetics protocols are available for species identification, even for samples with poor quality DNA, such as faeces, hairs or bones. Most of these protocols use mitochondrial DNA for barcoding despite this marker being especially prone to problems, including mtDNA introgression, nuclear insert copies, high intraspecific diversity or heteroplasmy. In this work, we developed a molecular method based on polymorphisms in small fragments of the mitochondrial cytochrome b (cytb, mtDNA) and the nuclear kappa casein genes (KCAS, nDNA) for identifying endangered North African ungulates. These fragments revealed polymorphisms, including species-specific variation, which allowed species identification of nine ungulate species that co-occur in North Africa. The method was validated across more than 400 samples, including different types of noninvasive samples collected in the field. The simplicity, high reliability and relative low cost of the described method make it a promising tool to improve ecological studies of the North African ungulates and consequently, the implementation of more efficient management and conservation plans for these endangered ungulates.

Characterization of the complete mitochondrial genome of Diaphania pyloalis (Lepidoptera: Pyralididae)

The complete mitochondrial genome (mitogenome) of Diaphania pyloalis (Lepidoptera: Pyralididae) was determined to be 15,298 bp and has the typical gene organization of mitogenomes from lepidopteran insects. It consists of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A+T-rich region. The A+T content of this mitogenome is 80.83% and the AT skew is slightly positive. All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which is initiated by CGA. Only the cox2 gene has an incomplete stop codon consisting of just a T. All the tRNA genes display a typical clover-leaf structure of mitochondrial tRNA. The A+T-rich region of the mitogenome is 332 bp in length, including several common features found in lepidopteran mitogenomes. Phylogenetic analysis showed that the D. pyloalis is close to Pyralididae.

The complete mitogenome of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) and comparison with other lepidopteran insects

The complete mitochondrial genome (mitogenome) of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) was determined to be 15,653bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a A+T-rich region. It has the typical gene organization and order of mitogenomes from lepidopteran insects. The AT skew of this mitogenome was slightly positive and the nucleotide composition was also biased toward A+T nucleotides (81.31%). All PCGs were initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which was initiated by CGA. The cox1 and cox2 genes had incomplete stop codons consisting of just a T. All the tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The A+T-rich region of the mitogenome was 495bp in length and consisted of several features common to the lepidopteras. Phylogenetic analysis showed that the B. mori Dazao was close to Bombycidae.

Comparative mitogenomic analysis of damsel bugs representing three tribes in the family Nabidae (Insecta: Hemiptera)

BACKGROUND

Nabidae, a family of predatory heteropterans, includes two subfamilies and five tribes. We previously reported the complete mitogenome of Alloeorhynchus bakeri, a representative of the tribe Prostemmatini in the subfamily Prostemmatinae. To gain a better understanding of architecture and evolution of mitogenome in Nabidae, mitogenomes of five species representing two tribes (Gorpini and Nabini) in the subfamily Nabinae were sequenced, and a comparative mitogenomic analysis of three nabid tribes in two subfamilies was carried out.

METHODOLOGY/PRINCIPAL FINDINGS

Nabid mitogenomes share a similar nucleotide composition and base bias, except for the control region, where differences are observed at the subfamily level. In addition, the pattern of codon usage is influenced by the GC content and consistent with the standard invertebrate mitochondrial genetic code and the preference for A+T-rich codons. The comparison among orthologous protein-coding genes shows that different genes have been subject to different rates of molecular evolution correlated with the GC content. The stems and anticodon loops of tRNAs are extremely conserved, and the nucleotide substitutions are largely restricted to TψC and DHU loops and extra arms, with insertion-deletion polymorphisms. Comparative analysis shows similar rates of substitution between the two rRNAs. Long non-coding regions are observed in most Gorpini and Nabini mtDNAs in-between trnI-trnQ and/or trnS2-nad1. The lone exception, Nabis apicalis, however, has lost three tRNAs. Overall, phylogenetic analysis using mitogenomic data is consistent with phylogenies constructed mainly form morphological traits.

CONCLUSIONS/SIGNIFICANCE

This comparative mitogenomic analysis sheds light on the architecture and evolution of mitogenomes in the family Nabidae. Nucleotide diversity and mitogenomic traits are phylogenetically informative at subfamily level. Furthermore, inclusion of a broader range of samples representing various taxonomic levels is critical for the understanding of mitogenomic evolution in damsel bugs.

The complete mitochondrial genome of the wild silkworm moth, Actias selene

The complete mitochondrial genome (mitogenome) of Actias selene (Lepidoptera: Saturniidae) was determined to be 15,236 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. The arrangement of 13 PCGs was similar to that of other sequenced lepidopterans. The AT skew of the mitogenome of A. selene was slightly negative, indicating a higher number of T compared to A nucleotides. The nucleotide composition of the mitogenome of A. selene was also biased toward A+T nucleotides (78.91%). All PCGs were initiated by ATN codons, except for the gene encoding cytochrome c oxidase subunit 1 (cox1), which may be initiated by the TTAG, as observed in other lepidopterans. Three genes, including cox1, cox2, and nad5, had incomplete stop codons consisting of just a T. With an exception for trnS1(AGN), all the other tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The A+T-rich region of the mitogenome of A. selene was 339 bp in length, and contains several features common to the Lepidopteras, including non-repetitive sequences, a conserved structure combining the motif ATAGA and an 18-bp poly-T stretch and a poly-A element upstream of trnM gene. Phylogenetic analysis showed that A. selene was close to Saturniidae.

The complete mitochondrial genome of the fall webworm, Hyphantria cunea (Lepidoptera: Arctiidae)

The complete mitochondrial genome (mitogenome) of the fall webworm, Hyphantria cunea (Lepidoptera: Arctiidae) was determined. The genome is a circular molecule 15 481 bp long. It presents a typical gene organization and order for completely sequenced lepidopteran mitogenomes, but differs from the insect ancestral type for the placement of tRNA(Met). The nucleotide composition of the genome is also highly A + T biased, accounting for 80.38%, with a slightly positive AT skewness (0.010), indicating the occurrence of more As than Ts, as found in the Noctuoidea species. All protein-coding genes (PCGs) are initiated by ATN codons, except for COI, which is tentatively designated by the CGA codon as observed in other lepidopterans. Four of 13 PCGs harbor the incomplete termination codon, T or TA. All tRNAs have a typical clover-leaf structure of mitochondrial tRNAs, except for tRNA(Ser)(AGN), the DHU arm of which could not form a stable stem-loop structure. The intergenic spacer sequence between tRNA(Ser)(AGN) and ND1 also contains the ATACTAA motif, which is conserved across the Lepidoptera order. The H. cunea A+T-rich region of 357 bp is comprised of non-repetitive sequences, but harbors several features common to the Lepidoptera insects, including the motif ATAGA followed by an 18 bp poly-T stretch, a microsatellite-like (AT)(8) element preceded by the ATTTA motif, an 11 bp poly-A present immediately upstream tRNA(Met). The phylogenetic analyses support the view that the H. cunea is closerly related to the Lymantria dispar than Ochrogaster lunifer, and support the hypothesis that Noctuoidea (H. cunea, L. dispar, and O. lunifer) and Geometroidea (Phthonandria atrilineata) are monophyletic. However, in the phylogenetic trees based on mitogenome sequences among the lepidopteran superfamilies, Papillonoidea (Artogeia melete, Acraea issoria, and Coreana raphaelis) joined basally within the monophyly of Lepidoptera, which is different to the traditional classification.

The complete nucleotide sequence of the mitochondrial genome of the cabbage butterfly, Artogeia melete (Lepidoptera: Pieridae)

The complete mitochondrial genome (mitogenome) of Artogeia melete was determined as being composed of 15,140 bp, including 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and one control region. The gene order of A. melete mitogenome is typical of Lepidoptera and differs from the insect ancestral type in the location of trnM. The A. melete mitogenome has a total of 119 bp of intergenic spacer sequences spread over 10 regions, ranging in sizes between 1 and 48 bp. The nucleotide composition of the A. melete mitogenome is also biased toward A 1 T nucleotides (79.77%), which is higher than that of Ochrogaster lunifer (77.84%), but lower than nine other lepidopterans sequenced. The PCGs have typical mitochondrial start codons, except for cox1, which contains the unusual CGA. The cox1, cox2, nad2, and nad5 genes of the A. melete mitogenome have incomplete stop codons (T). The A. melete A 1 T-rich region contains some conserved structures that are similar to those found in other lepidopteran mitogenomes, including a structure combining the motif ATAGA, a 19-bp poly(T) stretch, a microsatellite (AT)n element, and a 9-bp poly(A) upstream trnM. The A. melete mitogenome contains a duplicated 36-bp repeat element, which consists of a 26- bp core sequence flanked by 10-bp perfectly inverted repeats.

The complete mitochondrial genome of the bag-shelter moth Ochrogaster lunifer (Lepidoptera, Notodontidae)

BACKGROUND

Knowledge of animal mitochondrial genomes is very important to understand their molecular evolution as well as for phylogenetic and population genetic studies. The Lepidoptera encompasses more than 160,000 described species and is one of the largest insect orders. To date only nine lepidopteran mitochondrial DNAs have been fully and two others partly sequenced. Furthermore the taxon sampling is very scant. Thus advance of lepidopteran mitogenomics deeply requires new genomes derived from a broad taxon sampling. In present work we describe the mitochondrial genome of the moth Ochrogaster lunifer.

RESULTS

The mitochondrial genome of O. lunifer is a circular molecule 15593 bp long. It includes the entire set of 37 genes usually present in animal mitochondrial genomes. It contains also 7 intergenic spacers. The gene order of the newly sequenced genome is that typical for Lepidoptera and differs from the insect ancestral type for the placement of trnM. The 77.84% A+T content of its alpha strand is the lowest among known lepidopteran genomes. The mitochondrial genome of O. lunifer exhibits one of the most marked C-skew among available insect Pterygota genomes. The protein-coding genes have typical mitochondrial start codons except for cox1 that present an unusual CGA. The O. lunifer genome exhibits the less biased synonymous codon usage among lepidopterans. Comparative genomics analysis study identified atp6, cox1, cox2 as cox3, cob, nad1, nad2, nad4, and nad5 as potential markers for population genetics/phylogenetics studies. A peculiar feature of O. lunifer mitochondrial genome it that the intergenic spacers are mostly made by repetitive sequences.

CONCLUSION

The mitochondrial genome of O. lunifer is the first representative of superfamily Noctuoidea that account for about 40% of all described Lepidoptera. New genome shares many features with other known lepidopteran genomes. It differs however for its low A+T content and marked C-skew. Compared to other lepidopteran genomes it is less biased in synonymous codon usage. Comparative evolutionary analysis of lepidopteran mitochondrial genomes allowed the identification of previously neglected coding genes as potential phylogenetic markers. Presence of repetitive elements in intergenic spacers of O. lunifer genome supports the role of DNA slippage as possible mechanism to produce spacers during replication.

Structure and evolution of the mitochondrial control region of the pollen beetleMeligethes thalassophilus(Coleoptera: Nitidulidae)

The organization of the mitochondrial DNA (mtDNA) control region (CR) of the pollen beetle Meligethes thalassophilus is described. This mtDNA CR represents the longest sequenced for beetles so far, since the entire nucleotide sequence ranges from ~5000 to ~5500 bp. The CR of M. thalassophilus is organized in three distinct domains: a conserved domain near the tRNAIlegene, a variable domain flanking the 12S rRNA gene, and a relatively large central tandem array made up of a variable number of ~170 bp repeats that is responsible for the intraspecific length variation observed. Like other CRs found in insects, the M. thalassophilus CR contains two long homopolymeric runs that may be involved in mtDNA replication. Furthermore, conserved stem-and-loop structures in the repetitive domain were identified and their possible role in generating length variation is examined. Intraspecific comparison of the tandem repeat elements of M. thalassophilus suggests mechanisms of concerted evolution leading to homogenization of the repetitive region. The utility of such an array of tandem repeats as a genetic marker for assessing population-level variability and evolutionary relationships among populations is discussed. Finally, the technical difficulties found in isolating the mtDNA CR in beetles are remarked upon.

The complete mitochondrial genome of the tobacco hornworm, Manduca sexta, (Insecta: Lepidoptera: Sphingidae), and an examination of mitochondrial gene variability within butterflies and moths

The entire mitochondrial genome of the tobacco hornworm, Manduca sexta (Lepidoptera: Spinghidae) was sequenced -- a circular molecular 15516 bp in size. The arrangement of the protein coding genes (PCGs) was the same as that found in the ancestral insect, however Manduca possessed the derived tRNA arrangement of CR-M-I-Q which has been found in all Lepidoptera sequenced to date. Additionally, Manduca, like all lepidopteran mt genomes, has numerous large intergenic spacer regions and microsatellite-like repeat regions. Nucleotide composition is highly A+T biased, and the lepidopterans have the second most biased nucleotide composition of the insect orders after Hymenoptera. Secondary structural features of the PCGs identified in other Lepidoptera were present but highly modified by the presence of microsatellite-like repeat regions which may significantly alter their function in the post-transcriptional modification of pre-mRNAs. Secondary structure models of the ribosomal RNA genes of Manduca are presented and are similar to those proposed for other insect orders. Conserved regions were identified within non-translated spacer regions which correspond to sites for the origin and termination of replication and transcription. Comparisons of gene variability across the order suggest that the mitochondrial genes most frequently used in phylogenetic analysis of the Lepidoptera, cox1 and cox2, are amongst the least variable genes in the genome and phylogenetic resolution could be improved by using alternative, higher variability genes such as nad2, nad3, nad4 and nad5.

Delimiting species boundaries and the conservation genetics of the endangered maritime ringlet butterfly (Coenonympha nipisiquit McDunnough)

Species delimitation is a difficult problem that has implications across organismal biology, yet no single method has proved wholly satisfactory. We tested the utility of combining species-delimitation methods based on phylogeny and gene flow statistics using two parapatric members of the Coenonympha tullia group as an example: the endangered maritime ringlet butterfly (Coenonympha nipisiquit McDunnough) and the common inornate ringlet butterfly (Coenonympha inornata Edwards). We reconstructed the phylogeny of the nearctic C. tullia-group taxa from mitochondrial DNA (mtDNA) sequences (cytochrome oxidase I and mitochondrial control region) to explore the ancestry of the C. nipisiquit lineage within the group. We investigated the extent of gene flow between the two taxa with F-statistics using 587 nuclear amplified fragment length polymorphism markers, accounting for the effect of potential scoring 'collisions' where a marker may represent more than one DNA fragment. Combining species-delimitation methods was especially effective because it uncovered both historical and recent evolutionary patterns. Phylogenetic analysis of mtDNA revealed the early divergence of C. nipisiquit from other C. tullia-group taxa, including the morphologically similar C. inornata. F-statistics and gene-by-gene introgression profiles demonstrated clear isolation between the two taxa and revealed strong population structure within C. nipisiquit. C. nipisiquit is the first taxon in the nearctic C. tullia group showing strong evidence of genetic isolation. The methods we used are relatively inexpensive and can be widely used to delimit taxonomic boundaries near the species level, both generally and in particular for taxa that may be targets of conservation efforts.

Does the DNA barcoding gap exist? - a case study in blue butterflies (Lepidoptera: Lycaenidae)

BACKGROUND

DNA barcoding, i.e. the use of a 648 bp section of the mitochondrial gene cytochrome c oxidase I, has recently been promoted as useful for the rapid identification and discovery of species. Its success is dependent either on the strength of the claim that interspecific variation exceeds intraspecific variation by one order of magnitude, thus establishing a "barcoding gap", or on the reciprocal monophyly of species.

RESULTS

We present an analysis of intra- and interspecific variation in the butterfly family Lycaenidae which includes a well-sampled clade (genus Agrodiaetus) with a peculiar characteristic: most of its members are karyologically differentiated from each other which facilitates the recognition of species as reproductively isolated units even in allopatric populations. The analysis shows that there is an 18% overlap in the range of intra- and interspecific COI sequence divergence due to low interspecific divergence between many closely related species. In a Neighbour-Joining tree profile approach which does not depend on a barcoding gap, but on comprehensive sampling of taxa and the reciprocal monophyly of species, at least 16% of specimens with conspecific sequences in the profile were misidentified. This is due to paraphyly or polyphyly of conspecific DNA sequences probably caused by incomplete lineage sorting.

CONCLUSION

Our results indicate that the "barcoding gap" is an artifact of insufficient sampling across taxa. Although DNA barcodes can help to identify and distinguish species, we advocate using them in combination with other data, since otherwise there would be a high probability that sequences are misidentified. Although high differences in DNA sequences can help to identify cryptic species, a high percentage of well-differentiated species has similar or even identical COI sequences and would be overlooked in an isolated DNA barcoding approach.

Short‐term impact of disturbance on genetic diversity and structure of Indonesian populations of the butterflyDrupadia thedain East Kalimantan

We investigated the short-term impact of disturbance on genetic diversity and structure of the tropical butterfly Drupadia theda Felder (Lepidoptera: Lycaenidae). Populations were sampled from five landscapes in East Kalimantan (Borneo, Indonesia) which were differentially disturbed by selective logging and the 1997/1998 El Niño Southern Oscillation (ENSO)-induced drought and fires. Sampling occurred before (in 1997) and after the forest fires (in 1998, 1999, 2000, and 2004). Drupadia theda populations underwent serious population size reductions following the 1997/1998 ENSO event. For a total of 208 individuals, we sequenced a 509-bp segment of mtDNA containing the control region plus the 5' end of the 12S rDNA gene. Haplotype diversity in D. theda populations ranged from 0.468 to 0.953. Just after the 1997/1998 ENSO event, number of recorded individuals and genetic diversity were very low in D. theda populations sampled in the two severely burned areas and in a small pristine forest fragment that was surrounded by burned forest and thereby affected by drought. Interestingly, higher levels of genetic diversity were observed in logged forest compared to proximate pristine forest. After 1998, the genetic composition within the three ENSO-disturbed areas diverged. In the twice-burned forest, the genetic diversity in 1999 already approached pre-fire levels, while it remained nearly unchanged in proximate once-burned forest. Our data suggest that the 1997/1998 ENSO-induced drought and fires caused massive reductions in the genetic diversity of D. theda and that population recoveries were linked to their geographical position relative to patches of unburned forest (and thus to source populations).

Genetic Approaches for Studying Myiasis-causing Flies: Molecular Markers and Mitochondrial Genomics

"Myiasis-causing flies" is a generic term that includes species from numerous dipteran families, mainly Calliphoridae and Oestridae, of which blowflies, screwworm flies and botflies are among the most important. This group of flies is characterized by the ability of their larvae to develop in animal flesh. When the host is a live vertebrate, such parasitism by dipterous larvae is known as primary myiasis. Myiasis-causing flies can be classified as saprophagous (free-living species), facultative or obligate parasites. Many of these flies are of great medical and veterinary importance in Brazil because of their role as key livestock insect-pests and vectors of pathogens, in addition to being considered important legal evidence in forensic entomology. The characterization of myiasis-causing flies using molecular markers to study mtDNA (by RFLP) and nuclear DNA (by RAPD and microsatellite) has been used to identify the evolutionary mechanisms responsible for specific patterns of genetic variability. These approaches have been successfully used to analyze the population structures of the New World screwworm fly Cochliomyia hominivorax and the botfly Dermatobia hominis. In this review, various aspects of the organization, evolution and potential applications of the mitochondrial genome of myiasis-causing flies in Brazil, and the analysis of nuclear markers in genetic studies of populations, are discussed.

The importance of time scale and multiple refugia: Incipient speciation and admixture of lineages in the butterfly Erebia triaria (Nymphalidae)

Evidence of four different Quaternary glacial stages has been found in NW Iberia. The different magnitude of these episodes probably conditioned the distribution of currently montane species. We examined if the population history of NW Iberian Erebia triaria butterflies reflected such an influence of different glacial stages. We also investigated whether these populations were post-glacially colonised from a single refugial area or several sources. For this, we performed phylogenetic analyses and coalescence simulations on mitochondrial DNA sequences of individuals from five NW Iberian locations. We analysed three additional populations as reference, i.e., Central Spain, the Pyrenees, and the Alps. One of the NW Iberian populations, a subspecies endemic to a particular mountain range (Xistral), showed a high level of genetic divergence from all other populations, regardless of their geographic distance. Isolation after an ancient glacial stage, and followed by allopatric differentiation, may account for such high differentiation. The genetic pattern shown by the rest of the NW Iberian population samples consisted of two lineages, likely reflecting that the ancestors of these populations sought for refuge in at least two different areas during a subsequent glacial stage. We showed evidence of both temporal and spatial divisions in the phylogeography of these butterflies. The congruence of this pattern with the geological history suggests that the Iberian Peninsula hosted several refugial areas that differed in area and location depending on the glacial stage.

The history and geography of diversification within the butterfly genus Lycaeides in North America

The Lycaeides butterfly species complex in North America consists of two nominal, morphologically defined species. These butterflies are ecologically diverse and appear to be distributed as a geographically complex mosaic of locally differentiated populations that may be undergoing adaptive radiation. We asked whether patterns of molecular genetic variation within the species complex are congruent with currently recognized morphological species and whether the distribution of molecular variation is consistent with the hypothesis that Pleistocene climate changes contributed to the process of differentiation within the genus. Variation in the form of the genitalia from 726 males from 59 populations clearly distinguishes both species with only six populations containing morphologically intermediate or ambiguous individuals. However, partitioning of molecular variance in a 236 bp section of the mitochondrial AT-rich region from 628 individuals (57 populations) surveyed using single strand conformation polymorphism analysis (SSCP) indicates that only 26% of the total genetic variation is distributed along nominal species boundaries as defined by morphology. Instead, three phylogeographical groups were detected, represented by three major haplotype clades, which account for 90% of the total genetic variance. Pleistocene glaciations appear to have fostered divergence during glacial maxima, while post-glacial range expansions created opportunities for gene exchange and reticulation along suture zones between geographical groups. Data presented here allow us to make inferences about the history of the species complex. However, evidence of ancestral polymorphism and reticulation limit our ability to define species boundaries based on mitochondrial DNA sequence variation.

Testing biennialism in the butterfly Erebia palarica (Nymphalidae: Satyrinae) by mtDNA sequencing

Abstract The montane butterfly Erebia palarica has an unconfirmed life-cycle length. The genus Erebia shows variability in length of diapause depending on latitude and altitude, which may lead to allochronic differentiation. To test for biennialism and possible asynchronous differentiation, we sequenced the most variable fragments of the mitochondrial control region and the cytochrome oxidase I in forty individuals from the same population, from two consecutive years. No differentiation between cohorts could be detected and therefore the idea of a current annual cycle is supported. We conclude that a biennial cycle may never have been important in this population and we further discuss adaptive strategies this species may have used to confront colder climates.