Origin of intra-individual variation in PCR-amplified mitochondrial cytochrome oxidase I of Thrips tabaci (Thysanoptera: Thripidae): mitochondrial heteroplasmy or nuclear integration?

Thrips dynamics in Allium crops: Unraveling the role of reproductive mode and weather variables in Thrips tabaci population development

Thrips tabaci is the main thrips species affecting onion and related species. It is a cryptic species complex comprising three phylogenetic groups characterized by different reproductive modes (thelytoky or arrhenotoky) and host plant specialization. Thrips tabaci populations vary widely in genetic diversity, raising questions about the factor(s) that drive this diversity. We investigated the genetic diversity, reproductive mode, and heteroplasmy frequency in T. tabaci populations from different Allium spp fields in six locations in the Netherlands in 2021 and at two locations from the North and South of Uruguay over three years to unravel how the diversity is structured in the populations and if changes can be associated with weather variables. The thrips populations from each location studied were characterized by sequencing 33 individuals per sampling using the cytochrome oxidase subunit I (COI) gene. The reproductive mode was determined using specific primers and a phylogenetic analysis. Daily data for the weather variables was obtained from agrometeorological experimental stations located in the Uruguayan sampled crop fields. The diversity, reproductive mode, and heteroplasmy of T. tabaci populations in onion crops exhibited significant year-to-year variation depending on the location. Changes in the relative frequency of each reproductive mode in populations are associated with weather variables: precipitation, temperature, number of frosts, and relative humidity. Heteroplasmy frequency was associated with the same weather variables except temperature. In Uruguay and the Netherlands, T. tabaci thelytokous haplotype H1 was the most common, showing different heteroplasmy levels within and between the populations. In the field populations, a high frequency of heteroplasmic H1 individuals was associated with low precipitation, while all H1 individuals were also associated with high temperature and high relative humidity. In Uruguay, heteroplasmy was associated with arrhenotokous T. tabaci haplotypes, which were highly common in the North, pointing at specific adaptations leading to a faster population build-up. All this information may be instrumental for designing more precise integrated pest management techniques in both conventional and organic production.

Limited variation in microbial communities across populations of Macrosteles leafhoppers (Hemiptera: Cicadellidae)

Microbial symbionts play crucial roles in insect biology, yet their diversity, distribution, and temporal dynamics across host populations remain poorly understood. In this study, we investigated the spatio-temporal distribution of bacterial symbionts within the widely distributed and economically significant leafhopper genus Macrosteles, with a focus on Macrosteles laevis. Using host and symbiont marker gene amplicon sequencing, we explored the intricate relationships between these insects and their microbial partners. Our analysis of the cytochrome oxidase subunit I (COI) gene data revealed several intriguing findings. First, there was no strong genetic differentiation across M. laevis populations, suggesting gene flow among them. Second, we observed significant levels of heteroplasmy, indicating the presence of multiple mitochondrial haplotypes within individuals. Third, parasitoid infections were prevalent, highlighting the complex ecological interactions involving leafhoppers. The 16S rRNA data confirmed the universal presence of ancient nutritional endosymbionts-Sulcia and Nasuia-in M. laevis. Additionally, we found a high prevalence of Arsenophonus, another common symbiont. Interestingly, unlike most previously studied species, M. laevis exhibited only occasional cases of infection with known facultative endosymbionts and other bacteria. Notably, there was no significant variation in symbiont prevalence across different populations or among sampling years within the same population. Comparatively, facultative endosymbionts such as Rickettsia, Wolbachia, Cardinium and Lariskella were more common in other Macrosteles species. These findings underscore the importance of considering both host and symbiont dynamics when studying microbial associations. By simultaneously characterizing host and symbiont marker gene amplicons in large insect collections, we gain valuable insights into the intricate interplay between insects and their microbial partners. Understanding these dynamics contributes to our broader comprehension of host-microbe interactions in natural ecosystems.

Evolutionary origin and distribution of amino acid mutations associated with resistance to sodium channel modulators in onion thrips, Thrips tabaci

In onion thrips Thrips tabaci, reduced sensitivity of the sodium channel caused by several sodium channel mutations have been correlated with pyrethroid resistance. For this study, using mitochondrial cytochrome c oxidase subunit I gene sequences, we examined the phylogenetic relation among a total of 52 thelytokous and arrhenotokous strains with different genotypes of the sodium channel mutations. Then, we used flow cytometry to estimate their ploidy. Results showed that the strains are divisible into three groups: diploid thelytoky, triploid thelytoky, and diploid arrhenotoky. Using 23 whole genome resequencing data obtained from 20 strains out of 52, we examined their genetic relation further using principal component analysis, admixture analysis, and a fixation index. Results showed that diploid and triploid thelytokous groups are further classifiable into two based on the sodium channel mutations harbored by the respective group members (strains). The greatest genetic divergence was observed between thelytokous and arrhenotokous groups with a pair of T929I and K1774N. Nevertheless, they shared a genomic region with virtually no polymorphism around the sodium channel gene loci, suggesting a hard selective sweep. Based on these findings, we discuss the evolutionary origin and distribution of the sodium channel mutations in T. tabaci.

Frontiers Approaches to the Diagnosis of Thrips (Thysanoptera): How Effective Are the Molecular and Electronic Detection Platforms?

Simple Summary

Thrips are important agricultural and forest pests. They cause damage by sucking plant sap and transmitting several plant viruses. Correct identification is the key for epidemiological studies and formulating appropriate management strategies. The application of molecular and electronic detection platforms has improved the morphological character-based diagnosis of thrips species. This article reviews research on molecular and automated identification of thrips species and discusses future research strategies for rapid and high throughput thrips diagnosis.

Abstract

Thrips are insect pests of economically important agricultural, horticultural, and forest crops. They cause damage by sucking plant sap and by transmitting several tospoviruses, ilarviruses, carmoviruses, sobemoviruses, and machlomoviruses. Accurate and timely identification is the key to successful management of thrips species. However, their small size, cryptic nature, presence of color and reproductive morphs, and intraspecies genetic variability make the identification of thrips species challenging. The use of molecular and electronic detection platforms has made thrips identification rapid, precise, sensitive, high throughput, and independent of developmental stages. Multi-locus phylogeny based on mitochondrial, nuclear, and other markers has resolved ambiguities in morphologically indistinguishable thrips species. Microsatellite, RFLP, RAPD, AFLP, and CAPS markers have helped to explain population structure, gene flow, and intraspecies heterogeneity. Recent techniques such as LAMP and RPA have been employed for sensitive and on-site identification of thrips. Artificial neural networks and high throughput diagnostics facilitate automated identification. This review also discusses the potential of pyrosequencing, microarrays, high throughput sequencing, and electronic sensors in delimiting thrips species.

Mitochondrial heteroplasmy and pseudogenes in the freshwater prawn, Macrobrachium amazonicum (Heller, 1862): DNA barcoding and phylogeographic implications

The mitochondrial cytochrome oxidase c subunit 1 (COI) gene has been widely used in phylogenetic studies of crustaceans and analyses in population genetics. As COI studies have become more popular, there has been an increase in the number of reports of the presence of nuclear insertions of mitochondrial DNA (Numts) and mitochondrial heteroplasmy. Here, we provide evidence of both types of event in the COI sequences of Macrobrachium amazonicum, an economically important freshwater prawn, which is widespread in South America. Heteroplasmy and Numts were confirmed by different methods of DNA extraction (genomic, mitochondrial, and nuclear-enriched DNA), cloning, and sequencing, and were observed in 11 of the 14 populations sampled, primarily in the Amazon region. We discuss how the occurrence of these events affects the interpretation of the genetic relationships among the M. amazonicum populations, and we recommend caution when using COI for genetic inferences in prawns of the genus Macrobrachium, and in particular that any analysis should include nuclear markers.

DNA Barcoding: A Reliable Method for the Identification of Thrips Species (Thysanoptera, Thripidae) Collected on Sticky Traps in Onion Fields

Simple Summary

Thrips species (Insecta, Thysanoptera) identification using traditional approach is not an easy. In the present study, DNA barcoding was used to support the thrips species characterization of a wide collection sampled in onion fields. Our findings confirmed the selected method as a simple and accurate approach identifying major thrips species, characterizing successfully nearly 86% of the samples collected in nine main species. The results here reported underlined the role of genetic markers as a valuable and useful method for species identification, especially when the morphological approach is unsure or even impossible.

Abstract

Several thrips species (Insecta, Thysanoptera) are globally known as important crop pests and vectors of viral diseases, but their identification is difficult because of their small body size and inconspicuous morphological differences. Sequencing variation in the mitochondrial cytochrome c oxidase I (COI) region has been proven to be useful for the identification of species of many groups of insect pests. Here, DNA barcoding has been used to identify thrips species collected with the use of sticky traps placed in an open onion field. A total of 238 thrips specimens were analyzed, 151 of which could be identified to species and 27 to genera belonging to the family Thripidae. Fifty-one specimens could not be assigned to any genus, with the closest BLAST match in the GenBank queries being below 98%, whilst six specimens were not recognized as Thysanoptera. The results indicate that, although there are a few pest thrips species not yet barcoded, most of the species that may cause damage to crops in Europe are represented in GenBank and other databases, enabling correct identification. Additionally, DNA barcoding can be considered a valuable alternative to the classic morphology method for identification of major thrips species.

Novel Mitochondrial DNA Lineage Found among Ochlerotatus communis (De Geer, 1776) of the Nordic-Baltic Region

The Ochlerotatus (Oc.) communis complex consist of three Northern American species as well as a common Holarctic mosquito (Diptera: Culicidae) Oc. communis (De Geer, 1776). These sister species exhibit important ecological differences and are capable of transmitting various pathogens, but cannot always be differentiated by morphological traits. To investigate the Oc. communis complex in Europe, we compared three molecular markers (COI, ND5 and ITS2) from 54 Estonian mosquitoes as well as two COI marker sequences from Sweden. These sequences were subjected to phylogenetic analysis and screened for Wolbachia Hertig and Wolbach symbionts. Within and between groups, distances were calculated for each marker to better understand the relationships among individuals. Results demonstrate that a group of samples, extracted from adult female mosquitoes matching the morphology of Oc. communis, show a marked difference from the main species when comparing the mitochondrial markers COI and ND5. However, there is no variance between the same specimens when considering the nuclear ITS2. We conclude that Oc. communis encompasses two distinct mitochondrial DNA lineages in the Nordic-Baltic region. Further research is needed to investigate the origin and extent of these genetic differences.

Mitochondrial DNA intra-individual variation in a bumblebee species: A challenge for evolutionary studies and molecular identification

Mitochondrial DNA (mtDNA) regions have been widely used as molecular markers in evolutionary studies and species identification. However, the presence of heteroplasmy and NUMTs may represent obstacles. Heteroplasmy is a state where an organism has different mitochondrial haplotypes. NUMTs are nuclear pseudogenes originating from mtDNA sequences transferred to nuclear DNA. Evidences of heteroplasmy were already verified in the bumblebee Bombus morio in an earlier study. The present work investigated in more detail the presence of intra-individual haplotypes variation in this species. Heteroplasmy was detected in individuals from all the ten sampled locations, with an average of six heteroplasmic haplotypes per individual. In addition, some of these heteroplasmic haplotypes were shared among individuals from different locations, suggesting the existence of stable heteroplasmy in B. morio. These results demonstrated that heteroplasmy is likely to affect inferences based on mtDNA analysis, especially in phylogenetic, phylogeographic and population genetics studies. In addition, NUMTs were also detected. These sequences showed divergence of 2.7% to 12% in relation to the mitochondrial haplotypes. These levels of divergence could mislead conclusions in evolutionary studies and affect species identification through DNA barcoding.

A simple molecular identification method of the Thrips tabaci (Thysanoptera: Thripidae) cryptic species complex

The onion thrips (Thrips tabaci Lindeman, 1889) is a key pest of a wide range of crops because of its ecological attributes such as polyphagy, high reproduction rate, ability to transmit tospoviruses and resistance to insecticides. Recent studies revealed that T. tabaci is a cryptic species complex and it has three lineages (leek-associated arrhenotokous L1-biotype, leek-associated thelytokous L2-biotype and tobacco-associated arrhenotokous T-biotype), however, the adults remain indistinguishable. T. tabaci individuals were collected from different locations of Hungary to create laboratory colonies from each biotypes. Mitochondrial COI (mtCOI) region was sequenced from morphologically identified individuals. After sequence analysis SNPs were identified and used for CAPS marker development, which were suitable for distinguishing the three T. tabaci lineages. Genetic analysis of the T. tabaci species complex based on mtCOI gene confirmed the three well-known biotypes (L1, L2, T) and a new biotype because the new molecular evidence presented in this study suggests T-biotype of T. tabaci forming two distinct (sub)clades (T1 and T2). This genetic finding indicates that the genetic variability of T. tabaci populations is still not fully mapped. We validated our developed marker on thrips individuals from our thrips colonies. The results demonstrated that the new marker effectively identifies the different T. tabaci biotypes. We believe that our reliable genotyping method will be useful in further studies focusing on T. tabaci biotypes and in pest management by scanning the composition of sympatric T. tabaci populations.

Seasonally varying marine influences on the coastal ecosystem detected through molecular gut analysis

Terrestrial predators on marine shores benefit from the inflow of organisms and matter from the marine ecosystem, often causing very high predator densities and indirectly affecting the abundance of other prey species on shores. This indirect effect may be particularly strong if predators shift diets between seasons. We therefore quantified the seasonal variation in diet of two wolf spider species that dominate the shoreline predator community, using molecular gut content analyses with general primers to detect the full prey range. Across the season, spider diets changed, with predominantly terrestrial prey from May until July and predominantly marine prey (mainly chironomids) from August until October. This pattern coincided with a change in the spider age and size structure, and prey abundance data and resource selection analyses suggest that the higher consumption of chironomids during autumn is due to an ontogenetic diet shift rather than to variation in prey abundance. The analyses suggested that small dipterans with a weak flight capacity, such as Chironomidae, Sphaeroceridae, Scatopsidae and Ephydridae, were overrepresented in the gut of small juvenile spiders during autumn, whereas larger, more robust prey, such as Lepidoptera, Anthomyidae and Dolichopodidae, were overrepresented in the diet of adult spiders during spring. The effect of the inflow may be that the survival and growth of juvenile spiders is higher in areas with high chironomid abundances, leading to higher densities of adult spiders and higher predation rates on the terrestrial prey next spring.

The analysis of genetic variation in the mitochondrial genome and its application for the identification of Papilio species

Mitochondrial DNA (mtDNA) markers are ideal for evolutionary studies, including phylogeography, population genetics, phylogeny, etc. However, different mitochondrial genes always own different evolutionary rate. In this study, we analysed the genetic variation across the 16 complete mtDNA from 13 species in the genus Papilio and recognized the best DNA barcoding for Papilio species. The mitochondrial gene arrangement for each species shares a similar order, similar to the typical Papilionidae species, which indicated the relatively conservative state of gene arrangement in Papilio. The sliding window of genetic diversity showed that there was a significant difference in the genetic diversity of each gene in the mitochondrial genome of Papilio. The relatively mean clock rate of the ND1 was broadly lower than the other genes in mitochondrial genome of Papilio; while the ATP8 owns the largest values of mean clock rate. Those results suggested that the rate of evolution of each gene is not balanced and all mitochondrial genes except ND1 and ATP8 could act as barcoding for the identification of Papilio species. The phylogenetic analyses of complete mtDNA data for 13 Papilio species divided them into five major branches, which keep the same topological structure with previous studies.

Unexpected High Intragenomic Variation in Two of Three Major Pest Thrips Species Does Not Affect Ribosomal Internal Transcribed Spacer 2 (ITS2) Utility for Thrips Identification

The mitochondrial cytochrome oxidase I gene (mtCO1) and the ribosomal internal transcribed spacer 2 region (ITS2) are among the most widely used molecular markers for insect taxonomic characterization. Three economically important species of thrips, Scirtothripsdorsalis, Thripspalmi, and Frankliniellaoccidentalis were selected to examine the extent of intragenomic variation within these two marker regions in the family Thripidae, and determine if this variation would affect the utility of markers in thrips molecular diagnostics. For each species, intragenomic (within individual) variation and intergenomic (among individuals) variation was assessed by cloning and sequencing PCR-amplified copies. Intergenomic variation was generally higher than intragenomic variation except in cases where intergenomic variation was very low, as in mtCO1 from S.dorsalis and F.occidentalis. Intragenomic variation was detected in both markers in all three of the thrips species, however, 2-3 times more intragenomic variation was observed for ITS2 than mtCO1 in both S.dorsalis and T.palmi. Furthermore, levels of intragenomic variation were low for both of the genes in F.occidentalis. In all of the three thrips species, no sex-based clustering of haplotypes was observed in either marker. Unexpected high intragenomic variation in ITS2 for two of three thrips species did not interfere with thrips diagnostics. However, caution should be taken in applying ITS2 to certain studies of S.dorsalis and T.palmi when high levels of intragenomic variation could be problematic or confounding. In such studies, mtCO1 may be a preferable marker. Possible reasons for discrepancies in intragenomic variation among genomic regions are discussed.

Close congruence between Barcode Index Numbers (bins) and species boundaries in the Erebidae (Lepidoptera: Noctuoidea) of the Iberian Peninsula

The DNA barcode reference library for Lepidoptera holds much promise as a tool for taxonomic research and for providing the reliable identifications needed for conservation assessment programs. We gathered sequences for the barcode region of the mitochondrial cytochrome c oxidase subunit I gene from 160 of the 176 nominal species of Erebidae moths (Insecta: Lepidoptera) known from the Iberian Peninsula. These results arise from a research project which constructing a DNA barcode library for the insect species of Spain. New records for 271 specimens (122 species) are coupled with preexisting data for 38 species from the Iberian fauna. Mean interspecific distance was 12.1%, while the mean nearest neighbour divergence was 6.4%. All 160 species possessed diagnostic barcode sequences, but one pair of congeneric taxa (Eublemma rosea and Eublemma rietzi) were assigned to the same BIN. As well, intraspecific sequence divergences higher than 1.5% were detected in four species which likely represent species complexes. This study reinforces the effectiveness of DNA barcoding as a tool for monitoring biodiversity in particular geographical areas and the strong correspondence between sequence clusters delineated by BINs and species recognized through detailed taxonomic analysis.

Heteroplasmy due to coexistence of mtCOI haplotypes from different lineages of the Thrips tabaci cryptic species group

Heteroplasmy is the existence of multiple mitochondrial DNA haplotypes within the cell. Although the number of reports of heteroplasmy is increasing for arthropods, the occurrence, number of variants, and origins are not well studied. In this research, the occurrence of heteroplasmy was investigated in Thrips tabaci, a putative species complex whose lineages can be distinguished by their mitochondrial DNA haplotypes. The results from this study showed that heteroplasmy was due to the occurrence of mitochondrial cytochrome oxydase I (mtCOI) haplotypes from two different T. tabaci lineages. An assay using flow cytometry and quantitative real-time PCR was then used to quantify the per cell copy number of the two mtCOI haplotypes present in individuals exhibiting heteroplasmy from nine geographically distant populations in India. All of the T. tabaci individuals in this study were found to exhibit heteroplasmy, and in every individual the per cell copy number of mtCOI from lineage 3 comprised 75-98% of the haplotypes detected and was variable among individuals tested. There was no evidence to suggest that the presense of lineage-specific haplotypes was due to nuclear introgression; however, further studies are needed to investigate nuclear introgression and paternal leakage during rare interbreeding between individuals from lineages 2 and 3.

Detection and Localization of Wolbachia in Thrips palmi Karny (Thysanoptera: Thripidae)

Thrips palmi Karny is a globally distributed polyphagous agricultural pest. It causes huge economic loss by its biological behaviors like feeding, reproduction and transmission of tospoviruses. Since T. palmi shows close morphological similarities with other thrips species, we employed mitochondrial cytochrome oxidase 1 (mtCO1) gene as a molecular marker. BLAST analysis of this sequence helped us to identify the collected specimen as T. palmi. We observed the female to male ratio of about 3:1 from collected samples and suspected the presence of Wolbachia. The presence of Wolbachia was detected by PCR using genus specific primers of 16S rRNA gene. Further confirmation of Wolbachia strain was achieved by conducting PCR amplification of three ubiquitous genes ftsZ, gatB and groEL. A phylogenetic tree was constructed with concatenated sequences of ftsZ and gatB gene to assign supergroup to Wolbachia. Finally, we localized Wolbachia in abdominal region of the insect using fluorescent in situ hybridization with the help of confocal microscope. Our result confirmed the presence of Wolbachia supergroup B strain for the first time in T. palmi.

Population structure of the Indonesian giant tiger shrimp Penaeus monodon: a window into evolutionary similarities between paralogous mitochondrial DNA sequences and their genomes

Here we used both microsatellites and mtCR (mitochondrial DNA control region) sequences as genetic markers to examine the genetic diversity and population structure of Penaeus monodon shrimp from six Indonesian regions. The microsatellite data showed that shrimp from the Indian and the Pacific Ocean were genetically distinct from each other. It has been reported previously that P. monodon mtCR sequences from the Indo-Pacific group into two major paralogous clades of unclear origin. Here we show that the population structure inferred from mtCR sequences matches the microsatellite-based population structure for one of these clades. This is consistent with the notion that this mtCR clade shares evolutionary history with nuclear DNA and may thus represent nuclear mitochondrial pseudogenes (Numts).

Phylogenetic incongruence inferred with two mitochondrial genes in Mepraia spp. and Triatoma eratyrusiformis(Hemiptera, Reduviidae)

Mitochondrial DNA (mtDNA) is widely used to clarify phylogenetic relationships among and within species, and to determine population structure. Due to the linked nature of mtDNA genes it is expected that different genes will show similar results. Phylogenetic incongruence using mtDNA genes may result from processes such as heteroplasmy, nuclear integration of mitochondrial genes, polymerase errors, contamination, and recombination. In this study we used sequences from two mitochondrial genes (cytochrome b and cytochrome oxidase subunit I) from the wild vectors of Chagas disease, Triatoma eratyrusiformis and Mepraia species to test for topological congruence. The results showed some cases of phylogenetic incongruence due to misplacement of four haplotypes of four individuals. We discuss the possible causes of such incongruence and suggest that the explanation is an intra-individual variation likely due to heteroplasmy. This phenomenon is an independent evidence of common ancestry between these taxa.

DNA barcoding to fishes: current status and future directions

DNA barcoding appears to be a promising approach for taxonomic identification, characterization, and discovery of newer species, facilitating biodiversity studies. It helps researchers to appreciate genetic and evolutionary associations by collection of molecular, morphological, and distributional data. Fish DNA barcoding, based on the sequencing of a uniform area of Cytochrome C Oxidase type I (COI) gene, has received significant interest as an accurate tool for species identification, authentication, and phylogenetic analysis. The aim of this review article was to investigate recent global status, approaches, and future direction of DNA barcoding in fisheries sectors. We have tried to highlight its possible impacts, complications, and validation issues at species levels for biodiversity analysis. Moreover, an effort has been put forward to understand issues related to various marker genes associated with barcode process as primer sequences and have concluded barcode promotion as an indispensable tool of molecular biology for the development of taxonomic support systems.

Extent and divergence of heteroplasmy of the DNA barcoding region in Anapodisma miramae (Orthoptera: Acrididae)

A partial sequence of the mitochondrial cytochrome oxidase subunit I (COI) gene is widely used as a molecular marker for species identification in animals, also termed a DNA barcode. However, the presence of more than one sequence type in a single individual, also known as heteroplasmy, is one of the shortcomings of barcode identification. In this study, we examined the extent and divergence of COI heteroplasmy, including nuclear-encoded mitochondrial pseudogenes (NUMTs), at the genomic-DNA level from 13 insect species including orthopteran Anapodisma miramae, and a long fragment of mitochondrial DNA and cDNA from A. miramae as templates. When multiple numbers of clones originated from genomic DNA were sequenced, heteroplasmy was prevalent in all species and NUMTs were observed in five species. Long fragment DNA (∼13.5 kb) also is a source of heteroplasmic amplification, but the divergent haplotypes and NUMTs obtained from genomic DNA were not detected in A. miramae. On the other hand, cDNA was relatively heteroplasmy-free. Consistently, one dominant haplotype was always obtained from the genomic DNA-origin clones in all species and also from the long fragment- and cDNA-origin clones in the two tested individuals of A. miramae. Furthermore, the dominant haplotype was identical in sequence, regardless of the DNA source in A. miramae. Thus, one possible solution to avoid the barcoding problem in relationship to heteroplasmy could be the acquisition of multiple numbers of barcoding sequences to determine a dominant haplotype that can be assigned as barcoding sequence for a given species.

Identification of species-specific nuclear insertions of mitochondrial DNA (numts) in gorillas and their potential as population genetic markers

The first hyper-variable region (HV1) of the mitochondrial control region (MCR) has been widely used as a molecular tool in population genetics, but inadvertent amplification of nuclear translocated copies of mitochondrial DNA (numts) in gorillas has compromised the use of mitochondrial DNA in population genetic studies. At least three putative classes (I, II, III) of gorilla-specific HV1 MCR numts have been uncovered over the past decade. However, the number, size and location of numt loci in gorillas and other apes are completely unknown. Furthermore, little work to date has assessed the utility of numts as candidate population genetic markers. In the present study, we screened Bacterial Artificial Chromosome (BAC) genomic libraries in the chimpanzee and gorilla to compare patterns of mitochondrial-wide insertion in both taxa. We conducted an intensive BLAST search for numts in the gorilla genome and compared the prevalence of numt loci originating from the MCR with other great ape taxa. Additional gorilla-specific MCR numts were retrieved either through BAC library screens or using an anchored-PCR (A-PCR) amplification using genomic DNA from five unrelated gorillas. Locus-specific primers were designed to identify numt insertional polymorphisms and evaluate their potential as population genetic markers. Mitochondrial-wide surveys of chimpanzee and gorilla BACs showed that the number of numts does not differ between these two taxa. However, MCR numts are more abundant in chimpanzees than in other great apes. We identified and mapped 67 putative gorilla-specific numts, including two that contain the entire HV1 domain, cluster with sequences from two numt classes (I, IIb) and will likely co-amplify with mitochondrial sequences using most published HV1 primers. However, phylogenetic analysis coupled with post-hoc analysis of mitochondrial variation can successfully differentiate nuclear sequences. Insertional polymorphisms were evident in three out of five numts examined, indicating their potential utility as molecular markers. Taken together, these findings demonstrate the potentially powerful insight that numts could make in uncovering population history in gorillas and other mammals.

Molecular systematics of the Simulium jenningsi species group (Diptera: Simuliidae), with three new fast-evolving nuclear genes for phylogenetic inference

A molecular phylogeny was inferred for the 22 nominal species of black flies in the Simulium jenningsi species group, which includes major pests of humans and livestock in North America. Females are structurally monomorphic, presenting a problem for identification of the pests. For each species, we sequenced approximately two kilobases from the mitochondrial genome (ND2, Cox I, proximal one-half of Cox II) and about six kilobases from the nuclear genome (ca. 2 kilobases each from 3 rapidly evolving nuclear genes: big zinc finger [BZF], "5-intron gene" [5intG], and elongation complex protein 1 [ECP1]) and analyzed them phylogenetically using maximum likelihood and Bayesian methods. The three nuclear loci have not previously been used in phylogenetic studies. The mitochondrial region recovered 6 group members as monophyletic. BZF, 5intG, and ECP1 sequences each permitted identification of 13 species and recovered the S. fibrinflatum and S. taxodium subgroups. Simulium aranti Stone and Snoddy and S. luggeri Nicholson and Mickel were consistently recovered at the base of the group. Simulium ozarkense Moulton and Adler, S. dixiense Stone and Snoddy, S. krebsorum Moulton and Adler, and S. haysi Stone and Snoddy branched off before two well-supported sister groups of the remaining species. This remainder consisted of species occupying slow, sandy lowland streams-S. definitum Moulton and Adler, S. jonesi Stone and Snoddy, and the S. taxodium subgroup (S. taxodium Snoddy and Beshear, S. chlorum Moulton and Adler, S. confusum Moulton and Adler, and S. lakei Snoddy)-as sister to two clades of species inhabiting swift, rocky upland streams-the S. fibrinflatum subgroup (S. fibrinflatum Twinn, S. notiale Stone and Snoddy, and S. snowi Stone and Snoddy) and a clade comprised of S. anchistinum Moulton and Adler, S. jenningsi Malloch, and S. nyssa Stone and Snoddy, plus species having cocoons without anterolateral apertures (S. infenestrum Moulton and Adler, S. podostemi Snoddy, S. penobscotense Snoddy and Bauer, and S. remissum Moulton and Adler). Simulium snowi Stone and Snoddy is here considered a synonym of S. notiale Stone and Snoddy. Trees inferred from BZF and 5intG were largely concordant with those from ECP1, but slightly less resolved. Combining mitochondrial and nuclear data sets did not greatly improve the performance of the ECP1 data set alone. We, therefore, propose ECP1 as the gold standard for identification of members of the S. jenningsi group. Maximum likelihood analysis of combined sequences from all three nuclear genes, with three morphological constraints imposed, yielded a tree proposed as the best hypothesis of relationships among group members, based on all available data.

Genetic differentiation of the mitochondrial cytochrome oxidase C subunit I gene in genus Paramecium (Protista, Ciliophora)

Background

The mitochondrial cytochrome c oxidase subunit I (COI) gene is being used increasingly for evaluating inter- and intra-specific genetic diversity of ciliated protists. However, very few studies focus on assessing genetic divergence of the COI gene within individuals and how its presence might affect species identification and population structure analyses.

Methodology/Principal findings

We evaluated the genetic variation of the COI gene in five Paramecium species for a total of 147 clones derived from 21 individuals and 7 populations. We identified a total of 90 haplotypes with several individuals carrying more than one haplotype. Parsimony network and phylogenetic tree analyses revealed that intra-individual diversity had no effect in species identification and only a minor effect on population structure.

Conclusions

Our results suggest that the COI gene is a suitable marker for resolving inter- and intra-specific relationships of Paramecium spp.

Frequency matrix approach demonstrates high sequence quality in avian BARCODEs and highlights cryptic pseudogenes

The accuracy of DNA barcode databases is critical for research and practical applications. Here we apply a frequency matrix to assess sequencing errors in a very large set of avian BARCODEs. Using 11,000 sequences from 2,700 bird species, we show most avian cytochrome c oxidase I (COI) nucleotide and amino acid sequences vary within a narrow range. Except for third codon positions, nearly all (96%) sites were highly conserved or limited to two nucleotides or two amino acids. A large number of positions had very low frequency variants present in single individuals of a species; these were strongly concentrated at the ends of the barcode segment, consistent with sequencing error. In addition, a small fraction (0.1%) of BARCODEs had multiple very low frequency variants shared among individuals of a species; these were found to represent overlooked cryptic pseudogenes lacking stop codons. The calculated upper limit of sequencing error was 8 × 10(-5) errors/nucleotide, which was relatively high for direct Sanger sequencing of amplified DNA, but unlikely to compromise species identification. Our results confirm the high quality of the avian BARCODE database and demonstrate significant quality improvement in avian COI records deposited in GenBank over the past decade. This approach has potential application for genetic database quality control, discovery of cryptic pseudogenes, and studies of low-level genetic variation.

Mitochondrial heteroplasmy and DNA barcoding in Hawaiian Hylaeus (Nesoprosopis) bees (Hymenoptera: Colletidae)

BACKGROUND

The past several years have seen a flurry of papers seeking to clarify the utility and limits of DNA barcoding, particularly in areas such as species discovery and paralogy due to nuclear pseudogenes. Heteroplasmy, the coexistence of multiple mitochondrial haplotypes in a single organism, has been cited as a potentially serious problem for DNA barcoding but its effect on identification accuracy has not been tested. In addition, few studies of barcoding have tested a large group of closely-related species with a well-established morphological taxonomy. In this study we examine both of these issues, by densely sampling the Hawaiian Hylaeus bee radiation.

RESULTS

Individuals from 21 of the 49 a priori morphologically-defined species exhibited coding sequence heteroplasmy at levels of 1-6% or more. All homoplasmic species were successfully identified by COI using standard methods of analysis, but only 71% of heteroplasmic species. The success rate in identifying heteroplasmic species was increased to 86% by treating polymorphisms as character states rather than ambiguities. Nuclear pseudogenes (numts) were also present in four species, and were distinguishable from heteroplasmic sequences by patterns of nucleotide and amino acid change.

CONCLUSIONS

Heteroplasmy significantly decreased the reliability of species identification. In addition, the practical issue of dealing with large numbers of polymorphisms- and resulting increased time and labor required - makes the development of DNA barcode databases considerably more complex than has previously been suggested. The impact of heteroplasmy on the utility of DNA barcoding as a bulk specimen identification tool will depend upon its frequency across populations, which remains unknown. However, DNA barcoding is still likely to remain an important identification tool for those species that are difficult or impossible to identify through morphology, as is the case for the ecologically important solitary bee fauna.

Molecular identification of larval stages of Otiorhynchus (Coleoptera: Curculionidae) species based on polymerase chain reaction-restriction fragment length polymorphism analysis

A couple of different members of the coleopteran genus Otiorhynchus (Coleoptera: Curculionidae) are becoming increasingly important as pests of nursery and ornamental plants in global horticulture. Although adult weevils are morphologically distinguishable by skilled personnel, high potential for misidentification is given for cryptic larval stages. For developing and applying efficient pest management strategies the determination of the respective species is however a prerequisite, because each species may have a different phenology or a varying susceptibility to pesticides. Here, we report on the development of a diagnostic polymerase chain reaction (PCR)-restriction fragment length polymorphism (RFLP) method for differentiation among 16 Otiorhynchus and seven other weevil species independent of their developmental stage. An approximately 780-bp fragment of the mitochondrial cytochrome oxidase subunit II was amplified and subsequently digested with at most four restriction enzymes generating species-specific fragment patterns. The assay was validated on a total of 127 individuals and the obtained fragment patterns correctly identified 23 different weevil species. The PCR-RFLP method reported here is cost-effective, robust, and fast and could be used in the future by plant protection services for diagnostic purposes.

Characterization of nucleotide misincorporation patterns in the iceman's mitochondrial DNA

Background

The degradation of DNA represents one of the main issues in the genetic analysis of archeological specimens. In the recent years, a particular kind of post-mortem DNA modification giving rise to nucleotide misincorporation (“miscoding lesions”) has been the object of extensive investigations.

Methodology/Principal Findings

To improve our knowledge regarding the nature and incidence of ancient DNA nucleotide misincorporations, we have utilized 6,859 (629,975 bp) mitochondrial (mt) DNA sequences obtained from the 5,350–5,100-years-old, freeze-desiccated human mummy popularly known as the Tyrolean Iceman or Ötzi. To generate the sequences, we have applied a mixed PCR/pyrosequencing procedure allowing one to obtain a particularly high sequence coverage. As a control, we have produced further 8,982 (805,155 bp) mtDNA sequences from a contemporary specimen using the same system and starting from the same template copy number of the ancient sample. From the analysis of the nucleotide misincorporation rate in ancient, modern, and putative contaminant sequences, we observed that the rate of misincorporation is significantly lower in modern and putative contaminant sequence datasets than in ancient sequences. In contrast, type 2 transitions represent the vast majority (85%) of the observed nucleotide misincorporations in ancient sequences.

Conclusions/Significance

This study provides a further contribution to the knowledge of nucleotide misincorporation patterns in DNA sequences obtained from freeze-preserved archeological specimens. In the Iceman system, ancient sequences can be clearly distinguished from contaminants on the basis of nucleotide misincorporation rates. This observation confirms a previous identification of the ancient mummy sequences made on a purely phylogenetical basis. The present investigation provides further indication that the majority of ancient DNA damage is reflected by type 2 (cytosine→thymine/guanine→adenine) transitions and that type 1 transitions are essentially PCR artifacts.

Tissue segregation of mitochondrial haplotypes in heteroplasmic Hawaiian bees: implications for DNA barcoding

The issue of mitochondrial heteroplasmy has been cited as a theoretical problem for DNA barcoding but is only beginning to be examined in natural systems. We sequenced multiple DNA extractions from 20 individuals of four Hawaiian Hylaeus bee species known to be heteroplasmic. All species showed strong differences at polymorphic sites between abdominal and muscle tissue in most individuals, and only two individuals had no obvious segregation. Two specimens produced completely clean sequences from abdominal DNA. The fact that these differences are clearly visible by direct sequencing indicates that substantial intra-individual mtDNA diversity may be overlooked when DNA is taken from small tissue fragments. At the same time, differences in haplotype distribution among individuals may result in incorrect recognition of cryptic species. Because DNA barcoding studies typically use only a small fragment of an organism, they are particularly vulnerable to sequencing bias where heteroplasmy and haplotype segregation are present. It is important to anticipate this possibility prior to undertaking large-scale barcoding projects to reduce the likelihood of haplotype segregation confounding the results.

Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified

Nuclear mitochondrial pseudogenes (numts) are nonfunctional copies of mtDNA in the nucleus that have been found in major clades of eukaryotic organisms. They can be easily coamplified with orthologous mtDNA by using conserved universal primers; however, this is especially problematic for DNA barcoding, which attempts to characterize all living organisms by using a short fragment of the mitochondrial cytochrome c oxidase I (COI) gene. Here, we study the effect of numts on DNA barcoding based on phylogenetic and barcoding analyses of numt and mtDNA sequences in two divergent lineages of arthropods: grasshoppers and crayfish. Single individuals from both organisms have numts of the COI gene, many of which are highly divergent from orthologous mtDNA sequences, and DNA barcoding analysis incorrectly overestimates the number of unique species based on the standard metric of 3% sequence divergence. Removal of numts based on a careful examination of sequence characteristics, including indels, in-frame stop codons, and nucleotide composition, drastically reduces the incorrect inferences of the number of unique species, but even such rigorous quality control measures fail to identify certain numts. We also show that the distribution of numts is lineage-specific and the presence of numts cannot be known a priori. Whereas DNA barcoding strives for rapid and inexpensive generation of molecular species tags, we demonstrate that the presence of COI numts makes this goal difficult to achieve when numts are prevalent and can introduce serious ambiguity into DNA barcoding.

Molecular differences in the mitochondrial cytochrome oxidase I (mtCOI) gene and development of a species-specific marker for onion thrips, Thrips tabaci Lindeman, and melon thrips, T. palmi Karny (Thysanoptera: Thripidae), vectors of tospoviruses (Bunyaviridae)

A quick and developmental-stage non-limiting method of the identification of vectors of tospoviruses, such as Thrips tabaci and T. palmi, is important in the study of vector transmission, insecticide resistance, biological control, etc. Morphological identification of these thrips vectors is often a stumbling block in the absence of a specialist and limited by polymorphism, sex, stage of development, etc. Molecular identification, on the other hand, is not hampered by the above factors and can easily be followed by a non-specialist with a little training. The mitochondrial cytochrome oxidase I (mtCOI) exhibits reliable inter-species variations as compared to the other markers. In this communication, we present the differences in the mtCOI partial sequence of morphologically identified specimens of T. tabaci and T. palmi collected from onion and watermelon, respectively. Species-specific markers, identified in this study, could successfully determine T. tabaci and T. palmi, which corroborated the morphological identification. Phylogenetic analyses showed that both T. tabaci and T. palmi formed different clades as compared to the other NCBI accessions. The implication of these variations in vector efficiency has to be investigated further. The result of this investigation is useful in the quick identification of T. tabaci and T. palmi, a critical factor in understanding the epidemiology of the tospoviruses, their management and also in quarantine.

A genomic perspective on the shortcomings of mitochondrial DNA for "barcoding" identification

Approximately 600-bp sequences of mitochondrial DNA (mtDNA) have been designated as "DNA barcodes" and have become one of the most contentious and animated issues in the application of genetic information to global biodiversity assessment and species identification. Advocates of DNA barcodes have received extensive attention and promotion in many popular and refereed scientific publications. However, we suggest that the utility of barcodes is suspect and vulnerable to technical challenges that are particularly pertinent to mtDNA. We review the natural history of mtDNA and discuss problems for barcoding which are particularly associated with mtDNA and inheritance, including reduced effective population size, maternal inheritance, recombination, inconsistent mutation rate, heteroplasmy, and compounding evolutionary processes. The aforementioned could significantly limit the application and utility of mtDNA barcoding efforts. Furthermore, global use of barcodes will require application and acceptance of a barcode-based species concept that has not been evaluated in the context of the extensive literature concerning species designation. Implementation of mtDNA barcodes in spite of technical and practical shortcomings we discuss may degrade the longstanding synthesis of genetic and organism-based research and will not advance studies ranging from genomic evolution to biodiversity assessment.

Genetic differentiation of geographically separated populations of the southern green stink bug Nezara viridula (Hemiptera: Pentatomidae)

Genetic variation in the southern green stink bug Nezara viridula (Linnaeus) from 11 geographically separated sampling locations (Slovenia, France, Greece, Italy, Madeira, Japan, Guadeloupe, Galapagos, California, Brazil and Botswana) was studied by sequencing 16S and 28S rDNA, cytochrome b and cytochrome c oxidase subunit I gene fragments and random amplified polymorphic DNA (RAPD) analysis. Sequencing revealed 11 distinct haplotypes clustering into lineages A, B and C. Lineage C was characteristic for a single analysed specimen from Botswana. Lineage B was detected in Japan, and it probably arose in Asia. Haplotypes of European and American specimens belonged to lineage A; specimens from France, Slovenia, Madeira and Brazil shared highly similar haplotypes (>99%) from subgroup A1, while all the specimens from Greece, California, Galapagos and Guadeloupe shared a haplotype from subgroup A2. RAPD data were more variable but consistent with mtDNA sequences, revealing the same clustering. They separated the Botswanian specimen from Japanese specimens and from a group of more closely related specimens from Europe and America. Sequence and RAPD results both support the African origin of N. viridula, followed by dispersal to Asia (lineage B) and, more recently, by expansion to Europe and America (lineage A). RAPD analysis revealed two highly supported subgroups in Japan, congruent with mtDNA lineages A2 and B, suggesting multiple colonization of Japan. Invariant sequences at the 28S rDNA combined with other results do not support the hypothesis that cryptic (sibling) species exist within the populations investigated in this study.

Quantitative assessment of heteroplasmy levels in Senecio vulgaris chloroplast DNA

Heteroplasmy in coding chloroplast DNA was only recently shown to occur and was so far not quantitatively assessed. We present a quantitative analysis of cpDNA heteroplasmy levels at a triazine-resistance determining site within and between individual Senecio vulgaris plants. Detectable levels of heteroplasmic haplotypes were observed in all tested plants. As expected, the levels of heteroplasmy vary greatly between plants. However, even within individual plants, the fraction of one haplotype may cover a range from below 1% to well over 90%. Our results suggest that heteroplasmy may be a common phenomenon in S. vulgaris. Possible consequences for molecular diagnostics of chloroplast encoded traits as well as evolutionary consequences of chloroplast heteroplasmy are discussed.