Ghost mtDNA haplotypes generated by fortuitous NUMTs can deeply disturb infra-specific genetic diversity and phylogeographic pattern

Evaluating species boundaries using coalescent delimitation in pine-killing Monochamus (Coleoptera: Cerambycidae) sawyer beetles

Plant-feeding beetle species are diverse and often individually highly variable. Accurate classifications can be difficult to establish yet are essential for study of evolutionary patterns and processes. Molecular data are key to further characterizing morphologically difficult groups and defining genus and species boundaries. Monochamus Dejean species are ecologically and economically significant, and in coniferous forests they vector the nematode that causes Pine Wilt Disease. This study uses nuclear and mitochondrial genes to test the monophyly and relationships of Monochamus and applies coalescent methods to further delimit the conifer-feeding species. Monochamus has also included approximately 120 Old World species associated with diverse angiosperm tree species. We sample from these additional morphologically diverse species to determine their placement in the Lamiini. Through supermatrix and coalescent methods, the higher-level relationships of Monochamus show that conifer-feeders are a monophyletic group that includes the type species and has split into Nearctic and Palearctic clades. Molecular dating indicates a single dispersal of conifer-feeders to North America over the second Bering Land Bridge circa 5.3 Ma. All other Monochamus sampled fall in different parts of the Lamiini tree. Small-bodied angiosperm-feeding Monochamus group with the monotypic genus Microgoes Casey. The African Monochamus subgenera sampled are distantly related to the conifer-feeding clade. The multispecies coalescent delimitation methods BPP and STACEY delimit 17 conifer-feeding Monochamus species for a total of 18 species, and supports the retention of all current species. An interrogation with nuclear gene allele phasing reveals that unphased data can be unreliable for accurate delimitations and divergence times. The delimited species are discussed with integrative evidence, highlighting real-world challenges in recognizing the completion of speciation trajectories.

Dispersion routes of the main vectors of human malaria in the Americas: Genetic evidence from the mitochondrial COI gene

BACKGROUND AND OBJECTIVES

In America, of the 44 species of Anopheles, nine are main vectors of malaria and, of these, genetic information exists for seven. Hence, this study sought to know the gene flow and diversity of the seven principal vectors of malaria at the Americas level.

METHODS

For the seven species and the sequences of the mitochondrial cytochrome c oxidase I (COI) gene obtained from the GenBank and Bold System, genetic analyzes of populations and genetic structure were performed and haplotype networks and phylogenetic trees were obtained.

RESULTS

For the seven species, 1440 sequences were analyzed and 519 haplotypes were detected. The Hd and π values were higher within a continental context than by countries. Neutrality tests indicated positive and negative values with most of these being significant (p < 0.05). Phylogenetic analyses for all the species recovered three clades with no geographic pattern among them.

INTERPRETATION & CONCLUSION

Studies suggest that native species of Anopheles from the Americas have greater haplotype diversity and low genetic differentiation due to the lack of physical barriers to impede gene flow among these populations. Moreover, all the species are interconnected by roadways. This scenario complicates the epidemiological picture of malaria in the Americas.

Do pseudogenes pose a problem for metabarcoding marine animal communities?

Because DNA metabarcoding typically employs sequence diversity among mitochondrial amplicons to estimate species composition, nuclear mitochondrial pseudogenes (NUMTs) can inflate diversity. This study quantifies the incidence and attributes of NUMTs derived from the 658 bp barcode region of cytochrome c oxidase I (COI) in 156 marine animal genomes. NUMTs were examined to ascertain if they could be recognized by their possession of indels or stop codons. In total, 309 NUMTs ≥ 150 bp were detected, with an average of 1.98 per species (range = 0-33) and a mean length of 391 bp ± 200 bp. Among this total, 75 (24.3 %) lacked indels or stop codons. NUMTs appear to pose the greatest interpretational risk when short (< 313 bp) amplicons are used, such as in eDNA studies, dietary analyses, or processed fish identification. Employing the standard amplicon length (313 bp) for marine metabarcoding, NUMTs could potentially inflate the OTU count by 21% above the true species count while also raising intraspecific variation at COI by 15%. However, when both amplicon length and position are considered, inflation in OTU counts and in barcode variation were just 9% and 10%, respectively, suggesting NUMTs will not seriously distort biodiversity assessments. There was a weak positive correlation between genome size and NUMT count but no variation among phyla or trophic groups. Until bioinformatic advances improve NUMT detection, the best defense involves targeting long amplicons and developing reference databases that include both mitochondrial sequences and their NUMT derivatives.

Translocation of mitochondrial DNA into the nuclear genome blurs phylogeographic and conservation genetic studies in seabirds

Mitochondrial DNA (mtDNA) translocated into the nuclear genome (numt), when co-analysed with genuine mtDNA, could plague phylogeographic studies. To evaluate numt-related biases in population genetics parameters in birds, which are prone to accumulating numts, we targeted the mitochondrial mt-cytb gene. We looked at 13 populations of Audubon's shearwater (Puffinus lherminieri), including five mitochondrial lineages. mt-cytb homologue and paralogue (numt) sequences were determined by Sanger sequencing with and without prior exonuclease digestion of nuclear DNA. Numts formed monophyletic clades corresponding to three of the five mitochondrial lineages tested (the remaining two forming a paraphyletic group). Nineteen percent of numt alleles fell outside of their expected mitochondrial clade, a pattern consistent with multiple translocation events, incomplete lineage sorting (ILS), and/or introgression. When co-analysing mt-cytb paralogues and homologues, excluding individuals with ambiguities underestimates genetic diversity (4%) and differentiation (11%) among least-sampled populations. Removing ambiguous sites drops the proportion of inter-lineage genetic variance by 63%. While co-analysing numts with mitochondrial sequences can lead to severe bias and information loss in bird phylogeographic studies, the separate analysis of genuine mitochondrial loci and their nuclear paralogues can shed light on numt molecular evolution, as well as evolutionary processes such as ILS and introgression.

Barcoding pest species in a biodiversity hot-spot: the South African polyphagous broad-nosed weevils (Coleoptera, Curculionidae, Entiminae)

Polyphagous broad nosed weevils (Curculionidae: Entiminae) constitute a large and taxonomically challenging subfamily that contains economically significant agricultural pests worldwide. South Africa is a hot-spot for biodiversity and several species of indigenous and endemic genera of Entiminae have shifted on to cultivated plants, with some being phytosanitary pests. The sporadic pest status of many species (where the species has an occasional economic impact on the agricultural industry, but is not encountered often enough that is is readily recognisable by researchers and agricultural extension workers) and the presence of pest complexes and cryptic species represent an identification challenge to non-specialists. Furthermore, no comprehensive identification tools exist to identify immature stages that may be found in crops/soil. In this paper, a curated barcoding database with 70 COI sequences from 41 species (39 Entiminae, 2 Cyclominae) is initiated, to assist with the complexity of identification of species in this group.

DNA barcoding for bio-surveillance of emerging pests and species identification in Afrotropical Prioninae (Coleoptera, Cerambycidae)

DNA barcoding has been succesfully used for bio-surveillance of forest and agricultural pests in temperate areas, but has few applications in the tropics and particulary in Africa. Cacoscelesnewmannii (Coleoptera: Cerambycidae) is a Prioninae species that is locally causing extensive damage in commercially-grown sugarcane in the KwaZulu-Natal Province in South Africa. Due to the risk of spread of this species to the rest of southern Africa and to other sugarcane growing regions, clear and easy identification of this pest is critical for monitoring and for phytosanitary services. The genus Cacosceles Newman, 1838 includes four species, most being very similar in morphology. The damaging stage of the species is the larva, which is inherently difficult to distinguish morphologically from other Cerambycidae species. A tool for rapid and reliable identification of this species was needed by plant protection and quarantine agencies to monitor its potential abundance and spread. Here, we provide newly-generated barcodes for C.newmannii that can be used to reliably identify any life stage, even by non-trained taxonomists. In addition, we compiled a curated DNA barcoding reference library for 70 specimens of 20 named species of Afrotropical Prioninae to evaluate DNA barcoding as a valid tool to identify them. We also assessed the level of deeply conspecific mitochondrial lineages. Sequences were assigned to 42 different Barcode Index Numbers (BINs), 28 of which were new to BOLD. Out of the 20 named species barcoded, 11 (52.4%) had their own unique Barcode Index Number (BIN). Eight species (38.1%) showed multiple BINs with no morphological differentiation. Amongst them, C.newmannii showed two highly divergent genetic clusters which co-occur sympatrically, but further investigation is required to test whether they could represent new cryptic species.

Diversification, selective sweep, and body size in the invasive Palearctic alfalfa weevil infected with Wolbachia

The alfalfa weevil Hypera postica, native to the Western Palearctic, is an invasive legume pest with two divergent mitochondrial clades in its invading regions, the Western clade and the Eastern/Egyptian clade. However, knowledge regarding the native populations is limited. The Western clade is infected with the endosymbiotic bacteria Wolbachia that cause cytoplasmic incompatibility in host weevils. Our aim was to elucidate the spatial genetic structure of this insect and the effect of Wolbachia on its population diversity. We analyzed two mitochondrial and two nuclear genes of the weevil from its native ranges. The Western clade was distributed in western/central Europe, whereas the Eastern/Egyptian clade was distributed from the Mediterranean basin to central Asia. Intermediate mitotypes were found from the Balkans to central Asia. Most Western clade individuals in western Europe were infected with an identical Wolbachia strain. Mitochondrial genetic diversity of the infected individuals was minimal. The infected clades demonstrated a higher nonsynonymous/synonymous substitution rate ratio than the uninfected clades, suggesting a higher fixation of nonsynonymous mutations due to a selective sweep by Wolbachia. Trans-Mediterranean and within-European dispersal routes were supported. We suggest that the ancestral populations diversified by geographic isolation due to glaciations and that the diversity was reduced in the west by a recent Wolbachia-driven sweep(s). The intermediate clade exhibited a body size and host plant that differed from the other clades. Pros and cons of the possible use of infected-clade males to control uninfected populations are discussed.

Genetics of wild and mass-reared populations of a generalist aphid parasitoid and improvement of biological control

Due to their ability to parasitize various insect species, generalist parasitoids are widely used as biological control agents. They can be mass-reared and released in agroecosystems to control several pest species in various crops. However, the existence of genetic differentiation among populations of generalist parasitoid species is increasingly recognized and this can be associated with an adaptation to local conditions or to a reduced range of host species. Moreover, constraints of mass-rearing conditions can alter genetic variation within parasitoid populations released. These features could be associated with a reduced efficiency of the control of targeted pest species. Here, we focused on strawberry greenhouses where the control of aphids with the generalist parasitoid Aphidius ervi appears to be inefficient. We investigated whether this inefficiency may have both genetic and ecological bases comparing wild and commercial populations of A. ervi. We used two complementary genetic approaches: one based on the mitochondrial marker COI and one based on microsatellite markers. COI analysis showed a genetic differentiation within the A. ervi species, but the structure was neither associated with the commercial/wild status nor with host species factors. On the other hand, using microsatellite markers, we showed a genetic differentiation between commercial and wild A. ervi populations associated with a loss of genetic diversity within the mass-reared populations. Our ecological genetics study may potentially explain the weak efficiency of biological control of aphids in protected strawberry crops and enable to provide some insights to improve biological control.

Multilocus phylogeography of the world populations of Elaeidobius kamerunicus (Coleoptera, Curculionidae), pollinator of the palm Elaeis guineensis

Elaeidobius kamerunicus Faust (Coleoptera, Curculionidae) is one of the specific pollinators on inflorescences of the African oil palm Elaeis guineensis Jacquin. This derelomine weevil is native to tropical Africa. During the late 20th century, it was introduced into all tropical regions where E. guineensis is grown, in order to improve its pollination and fruit set. Despite an overall success, a decline in pollination efficiency has been documented in several regions. In this study, we reconstructed a multilocus phylogeography of the world populations of E. kamerunicus, in order to explore its genetic diversity in its native and introduced ranges. Our results showed that African populations of E. kamerunicus are forming two differentiated mitochondrial clusters in West and central Africa, forming a contact zone along the Cameroon Volcanic Line. The existence of this sharp contact zone along this weak altitudinal barrier suggests that other parameters, such as climate, may be driving the distribution of populations. A differential genetic structure between mitochondrial and nuclear genes, and the strong level of genetic structure of the mitochondrial gene, also suggest sex-biased dispersal in this species, with males dispersing more than females. The genetic structure inferred from Asian and South American populations suggests that they originate from populations of both western and central tropical Africa and that a bottleneck has probably been experienced by these populations.

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.

Hidden biodiversity revealed by integrated morphology and genetic species delimitation of spring dwelling water mite species (Acari, Parasitengona: Hydrachnidia)

BACKGROUND

Water mites are among the most diverse organisms inhabiting freshwater habitats and are considered as substantial part of the species communities in springs. As parasites, Hydrachnidia influence other invertebrates and play an important role in aquatic ecosystems. In Europe, 137 species are known to appear solely in or near springheads. New species are described frequently, especially with the help of molecular species identification and delimitation methods. The aim of this study was to verify the mainly morphology-based taxonomic knowledge of spring-inhabiting water mites of central Europe and to build a genetic species identification library.

METHODS

We sampled 65 crenobiontic species across the central Alps and tested the suitability of mitochondrial (cox1) and nuclear (28S) markers for species delimitation and identification purposes. To investigate both markers, distance- and phylogeny-based approaches were applied. The presence of a barcoding gap was tested by using the automated barcoding gap discovery tool and intra- and interspecific genetic distances were investigated. Furthermore, we analyzed phylogenetic relationships between different taxonomic levels.

RESULTS

A high degree of hidden diversity was observed. Seven taxa, morphologically identified as Bandakia concreta Thor, 1913, Hygrobates norvegicus (Thor, 1897), Ljania bipapillata Thor, 1898, Partnunia steinmanni Walter, 1906, Wandesia racovitzai Gledhill, 1970, Wandesia thori Schechtel, 1912 and Zschokkea oblonga Koenike, 1892, showed high intraspecific cox1 distances and each consisted of more than one phylogenetic clade. A clear intraspecific threshold between 5.6-6.0% K2P distance is suitable for species identification purposes. The monophyly of Hydrachnidia and the main superfamilies is evident with different species clearly separated into distinct clades. cox1 separates water mite species but is unsuitable for resolving higher taxonomic levels.

CONCLUSIONS

Water mite species richness in springs is higher than has been suggested based on morphological species identification alone and further research is needed to evaluate the true diversity. The standard molecular species identification marker cox1 can be used to identify species but should be complemented by a nuclear marker, e.g. 28S, to resolve taxonomic relationships. Our results contribute to the taxonomical knowledge on spring inhabiting Hydrachnida, which is indispensable for the development and implementation of modern environment assessment methods, e.g. metabarcoding, in spring ecology.

Beyond Biodiversity: Can Environmental DNA (eDNA) Cut It as a Population Genetics Tool?

Population genetic data underpin many studies of behavioral, ecological, and evolutionary processes in wild populations and contribute to effective conservation management. However, collecting genetic samples can be challenging when working with endangered, invasive, or cryptic species. Environmental DNA (eDNA) offers a way to sample genetic material non-invasively without requiring visual observation. While eDNA has been trialed extensively as a biodiversity and biosecurity monitoring tool with a strong taxonomic focus, it has yet to be fully explored as a means for obtaining population genetic information. Here, we review current research that employs eDNA approaches for the study of populations. We outline challenges facing eDNA-based population genetic methodologies, and suggest avenues of research for future developments. We advocate that with further optimizations, this emergent field holds great potential as part of the population genetics toolkit.

A DNA barcode library for ground beetles of Germany: the genus Amara Bonelli, 1810 (Insecta, Coleoptera, Carabidae)

The genus Amara Bonelli, 1810 is a very speciose and taxonomically difficult genus of the Carabidae. The identification of many of the species is accomplished with considerable difficulty, in particular for females and immature stages. In this study the effectiveness of DNA barcoding, the most popular method for molecular species identification, was examined to discriminate various species of this genus from Central Europe. DNA barcodes from 690 individuals and 47 species were analysed, including sequences from previous studies and more than 350 newly generated DNA barcodes. Our analysis revealed unique BINs for 38 species (81%). Interspecific K2P distances below 2.2% were found for three species pairs and one species trio, including haplotype sharing between Amara alpina/Amara torrida and Amara communis/Amara convexior/Amara makolskii. This study represents another step in generating an extensive reference library of DNA barcodes for carabids, highly valuable bioindicators for characterizing disturbances in various habitats.

Mansonella ozzardi mitogenome and pseudogene characterisation provides new perspectives on filarial parasite systematics and CO-1 barcoding

Despite the broad distribution of M. ozzardi in Latin America and the Caribbean, there is still very little DNA sequence data available to study this neglected parasite’s epidemiology. Mitochondrial DNA (mtDNA) sequences, especially the cytochrome oxidase (CO1) gene’s barcoding region, have been targeted successfully for filarial diagnostics and for epidemiological, ecological and evolutionary studies. MtDNA-based studies can, however, be compromised by unrecognised mitochondrial pseudogenes, such as Numts. Here, we have used shot-gun Illumina-HiSeq sequencing to recover the first complete Mansonella genus mitogenome and to identify several mitochondrial-origin pseudogenes. Mitogenome phylogenetic analysis placed M. ozzardi in the Onchocercidae “ONC5” clade and suggested that Mansonella parasites are more closely related to Wuchereria and Brugia genera parasites than they are to Loa genus parasites. DNA sequence alignments, BLAST searches and conceptual translations have been used to compliment phylogenetic analysis showing that M. ozzardi from the Amazon and Caribbean regions are near-identical and that previously reported Peruvian M. ozzardi CO1 reference sequences are probably of pseudogene origin. In addition to adding a much-needed resource to the Mansonella genus’s molecular tool-kit and providing evidence that some M. ozzardi CO1 sequence deposits are pseudogenes, our results suggest that all Neotropical M. ozzardi parasites are closely related.

Metabarcoding for the parallel identification of several hundred predators and their prey: Application to bat species diet analysis

Assessing diet variability is of main importance to better understand the biology of bats and design conservation strategies. Although the advent of metabarcoding has facilitated such analyses, this approach does not come without challenges. Biases may occur throughout the whole experiment, from fieldwork to biostatistics, resulting in the detection of false negatives, false positives or low taxonomic resolution. We detail a rigorous metabarcoding approach based on a short COI minibarcode and two-step PCR protocol enabling the "all at once" taxonomic identification of bats and their arthropod prey for several hundreds of samples. Our study includes faecal pellets collected in France from 357 bats representing 16 species, as well as insect mock communities that mimic bat meals of known composition, negative and positive controls. All samples were analysed using three replicates. We compare the efficiency of DNA extraction methods, and we evaluate the effectiveness of our protocol using identification success, taxonomic resolution, sensitivity and amplification biases. Our parallel identification strategy of predators and prey reduces the risk of mis-assigning prey to wrong predators and decreases the number of molecular steps. Controls and replicates enable to filter the data and limit the risk of false positives, hence guaranteeing high confidence results for both prey occurrence and bat species identification. We validate 551 COI variants from arthropod including 18 orders, 117 family, 282 genus and 290 species. Our method therefore provides a rapid, resolutive and cost-effective screening tool for addressing evolutionary ecological issues or developing "chirosurveillance" and conservation strategies.

Quantifying the Number of Independent Organelle DNA Insertions in Genome Evolution and Human Health

Fragments of organelle genomes are often found as insertions in nuclear DNA. These fragments of mitochondrial DNA (numts) and plastid DNA (nupts) are ubiquitous components of eukaryotic genomes. They are, however, often edited out during the genome assembly process, leading to systematic underestimation of their frequency. Numts and nupts, once inserted, can become further fragmented through subsequent insertion of mobile elements or other recombinational events that disrupt the continuity of the inserted sequence relative to the genuine organelle DNA copy. Because numts and nupts are typically identified through sequence comparison tools such as BLAST, disruption of insertions into smaller fragments can lead to systematic overestimation of numt and nupt frequencies. Accurate identification of numts and nupts is important, however, both for better understanding of their role during evolution, and for monitoring their increasingly evident role in human disease. Human populations are polymorphic for 141 numt loci, five numts are causal to genetic disease, and cancer genomic studies are revealing an abundance of numts associated with tumor progression. Here, we report investigation of salient parameters involved in obtaining accurate estimates of numt and nupt numbers in genome sequence data. Numts and nupts from 44 sequenced eukaryotic genomes reveal lineage-specific differences in the number, relative age and frequency of insertional events as well as lineage-specific dynamics of their postinsertional fragmentation. Our findings outline the main technical parameters influencing accurate identification and frequency estimation of numts in genomic studies pertinent to both evolution and human health.