Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: the effects of inherited symbionts

The Complete Mitochondrial Genome of the Luciocyprinus langsoni (Cypriniformes: Cyprinidae): Characterization, Phylogeny, and Genetic Diversity Analysis

Background: Luciocyprinus langsoni is a species belonging to the Cyprinidae family. The objective of this study is to gain a comprehensive understanding of its evolutionary history and genetic characteristics. Methods: The complete mitochondrial genome of L. langsoni was determined using overlapping PCR. A phylogenetic analysis was conducted based on 13 protein-coding genes from 48 species. A population genetic diversity analysis using the COI gene and a selection analysis of 13 protein-coding genes were also performed. Results: The mitogenome is 16,586 base pairs long and consists of 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and two control regions. It has a high adenine-thymine (A + T) content. The phylogenetic analysis confirms the placement of L. langsoni within the subfamily Cyprininae. The population genetic diversity analysis reveals low variability in the Hechi Longjiang population. The selection analysis shows that all 13 protein-coding genes have evolved under purifying selection with Ka/Ks ratios below 1. Conclusions: These results enhance our understanding of L. langsoni's evolutionary history and lay a genetic foundation for future studies in population genetics and phylogenetics.

Phylogenetics and species delimitation of the recluse spider, Loxosceles rufescens (Araneae: Sicariidae) populations invading Bangkok, Thailand

The Mediterranean recluse spider, Loxosceles rufescens, has been discovered for the first time inhabiting human dwellings in Bangkok, Thailand. Expeditions across 39 localities revealed five establishments with L. rufescens populations. The highest density was recorded in a storage house on Yaowarat Road, located in the heart of Bangkok's Chinatown, where 315 individuals were found, including adults, juveniles, and spiderlings. This medically significant spider's presence in such a densely populated urban area raises concerns about potential envenomation risks. Thirteen specimens of L. rufescens were extracted for DNA and sequenced for molecular phylogenetic analyses. COI and ITS2 markers were used to investigate relationships within L. rufescens and across available Loxosceles species sequences. Results indicate COI is superior for resolving species-level genetic clusters compared to ITS2. Surprisingly, L. rufescens individuals from the same house were found in significantly distant COI lineages, suggesting mtDNA may not be suitable for studying intra-specific phylogeography in this case. Species delimitation methods ABGD and ASAP demonstrated promising results for both COI and ITS2, while bPTP and GMYC tended to overestimate species numbers. ITS2 exhibited high sequence similarity in L. rufescens, suggesting potential utility as a barcoding marker for identification of this globally distributed species. Genetic distance analyses revealed a potential barcoding gap (K2P) of 8-9 % for COI and <2 % for ITS2 in Loxosceles. This study contributes valuable sequence data for the medically important genus Loxosceles and highlights the need for integrative approaches in understanding its evolution and spread. The findings have important implications for pest management strategies and public health in urban environments.

Intraspecific diversity of Myrmecophilus acervorum (Orthoptera: Myrmecophilidae) indicating an ongoing cryptic speciation

Myrmecophilus acervorum, previously considered a parthenogenetic species widely-distributed in Europe, has been observed to have both sexes in populations inhabiting the central part of the distribution range. Specimens from those heterosexual populations have been found being infected with Wolbachia. New mitochondrial data (COI and 16S markers) revealed the well-supported differentiation of M. acervorum populations inhabiting western Polesie (Poland) and southern Europe. In turn, analyses of EF1α marker support the hypothesis on the unfinished lineage sorting at the nuclear DNA level. Interestingly, we found that parthenogenetic populations inhabiting western Polesie are infected with Wolbachia belonging to supergroup A, while endosymbionts occurring in sexual populations of M. acervorum observed in Romania belong to supergroup B. Furthermore, new and potentially diagnostic characteristics in the external structures of the eyes of M. acervorum were identified. The surface of ommatidia in specimens occurring in southern Europe was smooth. In contrast, the ommatidia surface of individuals collected in Poland was visibly sculptured. To sum up, the significant genetic variability found in the present case, and the differentiating morphological character, are almost certainly effects of cryptic species being present within M. acervorum. This is indicative of ongoing speciation within the populations of this insect, and of simultaneous unfinished lineage sorting at the nuclear DNA level.

Spatially Varying Wolbachia Frequencies Reveal the Invasion Origin of an Agricultural Pest Recently Introduced From Europe to North America

The introduction of non-native species across the world represents a major global challenge. Retracing invasion origin is an important first step in understanding the invasion process, often requiring detailed sampling within the native range. Insect species frequently host Wolbachia, a widespread endosymbiotic bacterium that manipulates host reproduction to increase infected female fitness. Here, we draw on the spatial variation in infection frequencies of an actively spreading Wolbachia strain wCer2 to investigate the invasion origin of the European cherry fruit fly, Rhagoletis cerasi. This pest of cherries was introduced from Europe to North America within the last decade. First, we screen the introduced fly population for the presence of Wolbachia. The introduced populations lack the wCer2 strain and the strongly associated mitochondrial haplotype, suggesting strain absence due to founder effects with invading individuals originating from wCer2-uninfected native population(s). To narrow down geographic regions of invasion origin, we perform spatial interpolation of the wCer2 infection frequency across the native range and predict the infection frequency in unsampled regions. For this, we use an extensive dataset of R. cerasi infection covering 238 populations across Europe over 25 years, complemented with 14 additional populations analyzed for this study. We find that R. cerasi was unlikely introduced from wCer2-infected populations in Central and Western Europe. We propose wCer2-uninfected populations from Eastern Europe and the Mediterranean region as the most likely candidates for the invasion origin. This work utilizes Wolbachia as an indirect instrument to provide insights into the invasion source of R. cerasi in North America, revealing yet another application for this multifaceted heritable endosymbiont. Given the prevalence of biological invasions, rapidly uncovering invasion origins gives fundamental insights into how invasive species adapt to new environments.

Patterns of molecular evolution in a parthenogenic terrestrial isopod (Trichoniscus pusillus)

The "paradox of sex" refers to the question of why sexual reproduction is maintained in the wild, despite how costly it is compared to asexual reproduction. Because of these costs, one might expect nature to select for asexual reproduction, yet sex seems to be continually selected for. Multiple hypotheses have been proposed to explain this incongruence, including the niche differentiation hypothesis, the Red Queen hypothesis, and accumulation of harmful mutations in asexual species due to inefficient purifying selection. This study focuses on the accumulation of mutations in two terrestrial isopods, Trichoniscus pusillus, which has sexual diploid and parthenogenic triploid forms, and Hyloniscus riparius, an obligately sexual relative. We surveyed sex ratios of both species in an upstate New York population and obtained RNA-seq data from wild-caught individuals of both species to examine within- and between-species patterns of molecular evolution in protein-coding genes. The sex ratio and RNA-seq data together provide strong evidence that this T. pusillus population is entirely asexual and triploid, while the H. riparius population is sexual and diploid. Although all the wild T. pusillus individuals used for sequencing shared identical genotypes at nearly all SNPs, supporting a clonal origin, heterozygosity and SNP density were much higher in T. pusillus than in the sexually reproducing H. riparius. This observation suggests this parthenogenic lineage may have arisen via mating between two divergent diploid lineages. Between-species sequence comparisons showed no evidence of ineffective purifying selection in the asexual T. pusillus lineage, as measured by the ratio of nonsynonymous to synonymous substitutions (dN/dS ratios). Likewise, there was no difference between T. pusillus and H. riparius in the ratios of nonsynonymous to synonymous SNPs overall (pN/pS). However, pN/pS ratios in T. pusillus were significantly higher when considering only SNPs that may have arisen via recent mutation after the transition to parthenogenesis. Thus, these recent SNPs are consistent with the hypothesis that purifying selection is less effective against new mutations in asexual lineages, but only over long time scales. This system provides a useful model for future studies on the evolutionary tradeoffs between sexual and asexual reproduction in nature.

Delimiting species, revealing cryptic diversity, and population divergence in Qinghai-Tibet Plateau weevils through DNA barcoding

The Leptomias group represents one of the most diverse taxonomic group of weevils in the Qinghai-Tibet Plateau and its adjacent areas. Despite the potential of hidden diversity, relatively few comprehensive studies have been conducted on species diversity in this taxonomic group. In this study, we performed DNA barcoding analysis for species of the Leptomias group using a comprehensive DNA barcode dataset that included 476 sequences representing 54 morphospecies. Within the dataset, our laboratory contributed 474 sequences, and 390 sequences were newly generated for this study. The average Kimura 2-parameter distances among morphospecies and genera were 0.76% and 19.15%, respectively. In 94.4% of the species, the minimum interspecific distances exceeded the maximum intraspecific distances, indicating the presence of barcode gaps in most species of Leptomias group. The application of Automatic Barcode Gap Discovery, Assemble Species by Automatic Partitioning, Barcode Index Number, Bayesian Poisson tree processes, jMOTU, and Neighbor-joining tree methods revealed 45, 45, 63, 54, and 55 distinct clusters representing single species, respectively. Additionally, a total of four morphospecies, Leptomias kangmarensis, L. midlineatus, L. siahus, and L. sp.9RL, were found to be assigned to multiple subclade each, indicating the geographical divergences and the presence of cryptic diversity. Our findings of this study demonstrate that Qinghai-Tibet Plateau exhibits a higher species diversity of the Leptomias group, and it is imperative to investigate cryptic species within certain morphospecies using integrative taxonomic approaches in future studies. Moreover, the construction of a DNA barcode reference library presented herein establishes a robust foundational dataset to support forthcoming research on weevil taxonomy, phylogenetics, ecology, and evolution.

DNA barcoding reveals cryptic species in the sea slater Ligia italica (Crustacea, Isopoda) from Tunisia

Barcoding studies have provided significant insights into phylogenetic relationships among species belonging to the genus Ligia (Crustacea, Isopoda). Herein the diversity of the Italian sea slater Ligia italica from Tunisia is studied for the first time. Samples were collected from 18 localities in Tunisia, and the analysis included previously published sequences from Italy and Greece available in GenBank. Bayesian and Maximum Likelihood phylogenetic analyses were carried out using a fragment of the mitochondrial COI gene. Putative cryptic species were explored using the 'barcode gap' approach in the software ASAP. A genetic landscape shape analysis was carried out using the program Alleles in Space. The analyses revealed highly divergent and well-supported clades of L. italica dispersed across Tunisia (Clades A1 and A2), Greece (Clade B) and Italy (Clades C1 and C2). High genetic dissimilarity among clades suggested that L. italica constitute a cryptic species complex. Divergence among different L. italica lineages (Clades A, B and C) occurred around 7-4.5 Ma. The detected high genetic distances among clades did not result from atypical mitochondrial DNAs or intracellular infection by Wolbachia bacteria. The complex history of the Mediterranean Sea appears to have played a significant role in shaping the phylogeographic pattern of Ligia italica. Additional morphological and molecular studies are needed to confirm the existence of cryptic species in Ligia italica in Mediterranean.

Ectoparasite and bacterial population genetics and community structure indicate extent of bat movement across an island chain

Few studies have examined the genetic population structure of vector-borne microparasites in wildlife, making it unclear how much these systems can reveal about the movement of their associated hosts. This study examined the complex host–vector–microbe interactions in a system of bats, wingless ectoparasitic bat flies (Nycteribiidae), vector-borne microparasitic bacteria (Bartonella) and bacterial endosymbionts of flies (Enterobacterales) across an island chain in the Gulf of Guinea, West Africa. Limited population structure was found in bat flies and Enterobacterales symbionts compared to that of their hosts. Significant isolation by distance was observed in the dissimilarity of Bartonella communities detected in flies from sampled populations of Eidolon helvum bats. These patterns indicate that, while genetic dispersal of bats between islands is limited, some non-reproductive movements may lead to the dispersal of ectoparasites and associated microbes. This study deepens our knowledge of the phylogeography of African fruit bats, their ectoparasites and associated bacteria. The results presented could inform models of pathogen transmission in these bat populations and increase our theoretical understanding of community ecology in host–microbe systems.

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

Genomic signatures of hybridization between Ixodes ricinus and Ixodes persulcatus in natural populations

Identifying hybridization between common pathogen vectors is essential due to the major public health implications through risks associated with hybrid's enhanced pathogen transmission potential. The hard-ticks Ixodes ricinus and Ixodes persulcatus are the two most common vectors of tick-borne pathogens that affect human and animal health in Europe. Ixodes ricinus is a known native species in Finland with a well-known distribution, whereas I. persulcatus has expanded in range and abundance over the past 60 years, and currently it appears the most common tick species in certain areas in Finland. Here we used double-digest restriction site-associated DNA (ddRAD) sequencing on 186 ticks (morphologically identified as 92 I. ricinus, and 94 I. persulcatus) collected across Finland to investigate whether RAD generated single nucleotide polymorphisms (SNPs) can discriminate tick species and identify potential hybridization events. Two different clustering methods were used to assign specific species based on how they clustered and identified hybrids among them. We were able to discriminate between the two tick species and identified 11 putative hybrids with admixed genomic proportions ranging from approximately 24 to 76 percent. Four of these hybrids were morphologically identified as I. ricinus while the remaining seven were identified as I. persulcatus. Our results thus indicate that RAD SNPs are robust in identifying both species of the ticks as well as putative hybrids. These results further suggest ongoing hybridization between I. ricinus and I. persulcatus in their natural populations in Finland. The unique ability of RAD markers to discriminate between tick species and hybrids adds a useful aspect to tick evolutionary studies. Our findings align with previous studies and suggest a shared evolutionary history between the species, with instances of individuals possessing a considerable proportion of the other species' genome. This study is a significant step in understanding the formation of hybridization zones due to range expansion potentially associated with climate change.

Large-Scale Cytochrome C Oxidase Subunit I Gene Data Analysis for the Development of a Multiplex Polymerase Chain Reaction Test Capable of Identifying Biting Midge Vector Species and Haplotypes (Diptera: Ceratopogonidae) of the Culicoides Subgenus Avaritia Fox, 1955

The emergence of culicoid-transmitted bluetongue and Schmallenberg viruses in several European countries demonstrated the ability of indigenous biting midge species to transmit pathogens. Entomologic research programs identified members of the Obsoletus Group (Culicoides subgenus Avaritia) as keyplayers in disease epidemiology in Europe. However, morphological identification of potential vectors is challenging due to the recent discovery of new genetic variants (haplotypes) of C. obsoletus sensu stricto (s.s.), forming distinct clades. In this study, 4422 GenBank entries of the mitochondrial cytochrome c oxidase subunit I (COI) gene of subgenus Avaritia members of the genus Culicoides were analyzed to develop a conventional multiplex PCR, capable of detecting all vector species and clades of the Western Palearctic in this subgenus. Numerous GenBank entries incorrectly assigned to a species were identified, analyzed and reassigned. The results suggest that the three C. obsoletus clades represent independent species, whereas C. montanus should rather be regarded as a genetic variant of C. obsoletus s.s. Based on these findings, specific primers were designed and validated with DNA material from field-caught biting midges which achieved very high diagnostic sensitivity (100%) when compared to an established reference PCR (82.6%).

Conservation Mitonuclear Replacement: Facilitated mitochondrial adaptation for a changing world

Most species will not be able to migrate fast enough to cope with climate change, nor evolve quickly enough with current levels of genetic variation. Exacerbating the problem are anthropogenic influences on adaptive potential, including the prevention of gene flow through habitat fragmentation and the erosion of genetic diversity in small, bottlenecked populations. Facilitated adaptation, or assisted evolution, offers a way to augment adaptive genetic variation via artificial selection, induced hybridization, or genetic engineering. One key source of genetic variation, particularly for climatic adaptation, are the core metabolic genes encoded by the mitochondrial genome. These genes influence environmental tolerance to heat, drought, and hypoxia, but must interact intimately and co-evolve with a suite of important nuclear genes. These coadapted mitonuclear genes form some of the important reproductive barriers between species. Mitochondrial genomes can and do introgress between species in an adaptive manner, and they may co-introgress with nuclear genes important for maintaining mitonuclear compatibility. Managers should consider the relevance of mitonuclear genetic variability in conservation decision-making, including as a tool for facilitating adaptation. I propose a novel technique dubbed Conservation Mitonuclear Replacement (CmNR), which entails replacing the core metabolic machinery of a threatened species-the mitochondrial genome and key nuclear loci-with those from a closely related species or a divergent population, which may be better-adapted to climatic changes or carry a lower genetic load. The most feasible route to CmNR is to combine CRISPR-based nuclear genetic editing with mitochondrial replacement and assisted reproductive technologies. This method preserves much of an organism's phenotype and could allow populations to persist in the wild when no other suitable conservation options exist. The technique could be particularly important on mountaintops, where rising temperatures threaten an alarming number of species with almost certain extinction in the next century.

Facilitating taxonomy and phylogenetics: An informative and cost-effective protocol integrating long amplicon PCRs and third-generation sequencing

Phylogenetic inference has become a standard technique in integrative taxonomy and systematics, as well as in biogeography and ecology. DNA barcodes are often used for phylogenetic inference, despite being strongly limited due to their low number of informative sites. Also, because current DNA barcodes are based on a fraction of a single, fast-evolving gene, they are highly unsuitable for resolving deeper phylogenetic relationships due to saturation. In recent years, methods that analyse hundreds and thousands of loci at once have improved the resolution of the Tree of Life, but these methods require resources, experience and molecular laboratories that most taxonomists do not have. This paper introduces a PCR-based protocol that produces long amplicons of both slow- and fast-evolving unlinked mitochondrial and nuclear gene regions, which can be sequenced by the affordable and portable ONT MinION platform with low infrastructure or funding requirements. As a proof of concept, we inferred a phylogeny of a sample of 63 spider species from 20 families using our proposed protocol. The results were overall consistent with the results from approaches based on hundreds and thousands of loci, while requiring just a fraction of the cost and labour of such approaches, making our protocol accessible to taxonomists worldwide.

Species Identification of Rehabilitated Critically Endangered Orangutans Through DNA Forensic: Implication for Conservation

Rehabilitating and releasing orangutans back into the wild is one of the conservation strategies being pursued to conserve orangutans. However, the species determination between Sumatran, Tapanuli, and Bornean orangutans is essential for reintroduction to avoid outbreeding depression, which could lead to DNA hybridisation and increase the probability of recessive characters. Here, we reported on an investigation of three orangutans in which DNA forensic techniques were used to identify the species before release and reintroduction to their habitat. By applying DNA forensic, the orangutan was successfully confirmed with high probabilities (100%) by identifying two orangutan species, Pongo abelii and Pongo pygmaeus wurmbii. Based on ambiguous morphology, we found the possibility of orangutan species being misidentified in rehabilitation. This case report demonstrates the importance of molecular diagnostics to identify the orangutan species. We also provide workflow recommendations from genetic aspect for rehabilitated orangutans. These recommendations will enable decision-makers to consider genetics when assessing future management decisions, which will help ensure that the orangutan species is effectively conserved.

Validating a Mitochondrial Sweep Accompanying the Rapid Spread of a Maternally Inherited Symbiont

Maternally inherited symbiotic bacteria that interfere with the reproduction of their hosts can contribute to selective sweeps of mitochondrial haplotypes through hitch-hiking or coordinate inheritance of cytoplasmic bacteria and host mitochondria. The sweep will be manifested by genetic variations of mitochondrial genomic DNA of symbiont-infected hosts relative to their uninfected counterparts. In particular, at the population level, infected specimens will show a reduced mitochondrial DNA polymorphism compared to that in the nuclear DNA. This may challenge the use of mitochondrial DNA sequences as neutral genetic markers, as the mitochondrial patterns will reflect the evolutionary history of parasitism, rather than the sole evolutionary history of the host. Here, I describe a detailed step-by-step procedure to infer the occurrence and timing of symbiont-induced mitochondrial sweeps in host species.

Two sympatric lineages of Australian Cnestus solidus share Ambrosiella symbionts but not Wolbachia

Sympatric lineages of inbreeding species provide an excellent opportunity to investigate species divergence patterns and processes. Many ambrosia beetle lineages (Curculionidae: Scolytinae) reproduce by predominant inbreeding through sib mating in nests excavated in woody plant parts wherein they cultivate symbiotic ambrosia fungi as their sole source of nutrition. The Xyleborini ambrosia beetle species Cnestus solidus and Cnestus pseudosolidus are sympatrically distributed across eastern Australia and have overlapping morphological variation. Using multilocus sequencing analysis of individuals collected from 19 sites spanning their sympatric distribution, we assessed their phylogenetic relationships, taxonomic status and microbial symbionts. We found no genetic differentiation between individuals morphologically identified as C. solidus and C. pseudosolidus confirming previous suggestions that C. pseudosolidus is synonymous to C. solidus. However, within C. solidus we unexpectedly discovered the sympatric coexistence of two morphologically indistinguishable but genetically distinct lineages with small nuclear yet large mitochondrial divergence. At all sites except one, individuals of both lineages carried the same primary fungal symbiont, a new Ambrosiella species, indicating that fungal symbiont differentiation may not be involved in lineage divergence. One strain of the maternally inherited bacterial endosymbiont Wolbachia was found at high prevalence in individuals of the more common lineage but not in the other, suggesting that it may influence host fitness. Our data suggest that the two Australian Cnestus lineages diverged allopatrically, and one lineage then acquired Wolbachia. Predominant inbreeding and Wolbachia infection may have reinforced reproductive barriers between these two lineages after their secondary contact contributing to their current sympatric distribution.

Phylogeography of Two Enigmatic Sulphur Butterflies, Colias mongola Alphéraky, 1897 and Colias tamerlana Staudinger, 1897 (Lepidoptera, Pieridae), with Relations to Wolbachia Infection

The genus Colias Fabricius, 1807 includes numerous taxa and forms with uncertain status and taxonomic position. Among such taxa are Colias mongola Alphéraky, 1897 and Colias tamerlana Staudinger, 1897, interpreted in the literature either as conspecific forms, as subspecies of different but morphologically somewhat similar Colias species or as distinct species-level taxa. Based on mitochondrial and nuclear DNA markers, we reconstructed a phylogeographic pattern of the taxa in question. We recover and include in our analysis DNA barcodes of the century-old type specimens, the lectotype of C. tamerlana deposited in the Natural History Museum (Museum für Naturkunde), Berlin, Germany (ZMHU) and the paralectotype of C. tamerlana and the lectotype of C. mongola deposited in the Zoological Institute, Russian Academy of Sciences, St. Petersburg, Russia (ZISP). Our analysis grouped all specimens within four (HP_I-HP_IV) deeply divergent but geographically poorly structured clades which did not support nonconspecifity of C. mongola-C. tamerlana. We also show that all studied females of the widely distributed haplogroup HP_II were infected with a single Wolbachia strain belonging to the supergroup B, while the males of this haplogroup, as well as all other investigated specimens of both sexes, were not infected. Our data highlight the relevance of large-scale sampling dataset analysis and the need for testing for Wolbachia infection to avoid erroneous phylogenetic reconstructions and species misidentification.

Disentangling Positive Selection from Relaxed Selection in Animal Mitochondrial Genomes

AbstractDisentangling different types of selection is a common goal in molecular evolution. Elevated dN/dS ratios (the ratio of nonsynonymous to synonymous substitution rates) in focal lineages are often interpreted as signs of positive selection. Paradoxically, relaxed purifying selection can also result in elevated dN/dS ratios, but tests to distinguish these two causes are seldomly implemented. Here, we reevaluated seven case studies describing elevated dN/dS ratios in animal mitochondrial DNA (mtDNA) and their accompanying hypotheses regarding selection. They included flightless lineages versus flighted lineages in birds, bats, and insects and physiological adaptations in snakes, two groups of electric fishes, and primates. We found that elevated dN/dS ratios were often not caused by the predicted mechanism, and we sometimes found strong support for the opposite mechanism. We discuss reasons why energetic hypotheses may be confounded by other selective forces acting on mtDNA and caution against overinterpreting singular molecular signals, including elevated dN/dS ratios.

New species based on the biological species concept within the complex of Lariophagus distinguendus (Hymenoptera, Chalcidoidea, Pteromalidae), a parasitoid of household pests

The pteromalid parasitoid Lariophagus distinguendus (Foerster) belongs to the Hymenoptera, a megadiverse insect order with high cryptic diversity. It attacks stored product pest beetles in human storage facilities. Recently, it has been shown to consist of two separate species. To further study its cryptic diversity, strains were collected to compare their relatedness using barcoding and nuclear genes. Nuclear genes identified two clusters which agree with the known two species, whereas the barcode fragment determined an additional third Clade. Total reproductive isolation (RI) according to the biological species concept (BSC) was investigated in crossing experiments within and between clusters using representative strains. Sexual isolation exists between all studied pairs, increasing from slight to strong with genetic distance. Postzygotic barriers mostly affected hybrid males, pointing to Haldane's rule. Hybrid females were only affected by unidirectional Spiroplasma-induced cytoplasmic incompatibility and behavioural sterility, each in one specific strain combination. RI was virtually absent between strains separated by up to 2.8% COI difference, but strong or complete in three pairs from one Clade each, separated by at least 7.2%. Apparently, each of these clusters represents one separate species according to the BSC, highlighting cryptic diversity in direct vicinity to humans. In addition, these results challenge the recent 'turbo-taxonomy' practice of using 2% COI differences to delimitate species, especially within parasitic Hymenoptera. The gradual increase in number and strength of reproductive barriers between strains with increasing genetic distance also sheds light on the emergence of barriers during the speciation process in L. distinguendus.

On the Origins of Phenotypic Parallelism in Benthic and Limnetic Stickleback

Rapid evolution of similar phenotypes in similar environments, giving rise to in situ parallel adaptation, is an important hallmark of ecological speciation. However, what appears to be in situ adaptation can also arise by dispersal of divergent lineages from elsewhere. We test whether two contrasting phenotypes repeatedly evolved in parallel, or have a single origin, in an archetypal example of ecological adaptive radiation: benthic-limnetic three-spined stickleback (Gasterosteus aculeatus) across species pair and solitary lakes in British Columbia. We identify two genomic clusters across freshwater populations, which differ in benthic-limnetic divergent phenotypic traits and separate benthic from limnetic individuals in species pair lakes. Phylogenetic reconstruction and niche evolution modeling both suggest a single evolutionary origin for each of these clusters. We detected strong phylogenetic signal in benthic-limnetic divergent traits, suggesting that they are ancestrally retained. Accounting for ancestral state retention, we identify local adaptation of body armor due to the presence of an intraguild predator, the sculpin (Cottus asper), and environmental effects of lake depth and pH on body size. Taken together, our results imply a predominant role for retention of ancestral characteristics in driving trait distribution, with further selection imposed on some traits by environmental factors.

Macrogenetics reveals multifaceted influences of environmental variation on vertebrate population genetic diversity across the Americas

The broad scale distribution of population-specific genetic diversity (GDP ) across taxa remains understudied relative to species diversity gradients, despite its relevance for systematic conservation planning. We used nuclear DNA data collected from 3678 vertebrate populations across the Americas to assess the role of environmental and spatial variables in structuring the distribution of GDP , a key component of adaptive potential in the face of environmental change. We specifically assessed non-linear trends for a metric of GDP, expected heterozygosity (HE ), and found more evidence for spatial hotspots and cold spots in HE rather than a strict pattern with latitude. We also detected inconsistent relationships between HE and environmental variables, where only 11 of 30 environmental comparisons among taxa groups were statistically significant at the .05 level, and the shape of significant trends differed substantially across vertebrate groups. Only one of six taxonomic groups, freshwater fishes, consistently showed significant relationships between HE and most (four of five) environmental variables. The remaining groups had statistically significant relationships for either two (amphibians, reptiles), one (birds, mammals), or no variables (anadromous fishes). Our study highlights gaps in the theoretical foundation upon which macrogenetic predictions have been made thus far in the literature, as well as the nuances for assessing broad patterns in GDP among vertebrate groups. Overall, our results suggest a disconnect between patterns of species and genetic diversity, and underscores that large-scale factors affecting genetic diversity may not be the same factors as those shaping taxonomic diversity. Thus, careful spatial and taxonomic-specific considerations are needed for applying macrogenetics to conservation planning.

Multiple introgressions shape mitochondrial evolutionary history in Drosophila paulistorum and the Drosophila willistoni group

Hybridization and the consequent introgression of genomic elements is an important source of genetic diversity for biological lineages. This is particularly evident in young clades in which hybrid incompatibilities are still incomplete and mixing between species is more likely to occur. Drosophila paulistorum, a representative of the Neotropical Drosophila willistoni subgroup, is a classic model of incipient speciation. The species is divided into six semispecies that show varying degrees of pre- and post-mating incompatibility with each other. In the present study, we investigate the mitochondrial evolutionary history of D. paulistorum and the willistoni subgroup. For that, we perform phylogenetic and comparative analyses of the complete mitochondrial genomes and draft nuclear assemblies of 25 Drosophila lines of the willistoni and saltans species groups. Our results show that the mitochondria of D. paulistorum are polyphyletic and form two non-sister clades that we name α and β. Identification and analyses of nuclear mitochondrial insertions further reveal that the willistoni subgroup has an α-like mitochondrial ancestor and strongly suggest that both the α and β mitochondria of D. paulistorum were acquired through introgression from unknown fly lineages of the willistoni subgroup. We also uncover multiple mitochondrial introgressions across D. paulistorum semispecies and generate novel insight into the evolution of the species.

Molecular taxonomy of phlebotomine sand flies (Diptera, Psychodidae) with emphasis on DNA barcoding: A review

The taxonomy and systematics of sand flies (Diptera, Psychodidae, Phlebotominae) are one of the pillars of research aimed to identifying vector populations and the agents transmitted by these insects. Traditionally, the use of morphological traits has been the main line of evidence for the definition of species, but the use of DNA sequences is useful as an integrative approach for their delimitation. Here, we discuss the current status of the molecular taxonomy of sand flies, including their most sequenced molecular markers and the main results. Only about 37% of all sand fly species have been processed for any molecular marker and are publicly available in the NCBI GenBank or BOLD Systems databases. The genera Phlebotomus, Nyssomyia, Psathyromyia and Psychodopygus are well-sampled, accounting for more than 56% of their sequenced species. However, less than 34% of the species of Sergentomyia, Lutzomyia, Trichopygomyia and Trichophoromyia have been sampled, representing a major gap in the knowledge of these groups. The most sequenced molecular markers are those within mtDNA, especially the DNA barcoding fragment of the cytochrome c oxidase subunit I (coi) gene, which has shown promising results in detecting cryptic diversity within species. Few sequences of conserved genes have been generated, which hampers higher-level phylogenetic inferences. We argue that sand fly species should be sequenced for at least the coi DNA barcoding marker, but multiple markers with different mutation rates should be assessed, whenever possible, to generate multilocus analysis.

Ribosomal RNA (rRNA) sequences from 33 globally distributed mosquito species for improved metagenomics and species identification

Total RNA sequencing (RNA-seq) is an important tool in the study of mosquitoes and the RNA viruses they vector as it allows assessment of both host and viral RNA in specimens. However, there are two main constraints. First, as with many other species, abundant mosquito ribosomal RNA (rRNA) serves as the predominant template from which sequences are generated, meaning that the desired host and viral templates are sequenced far less. Second, mosquito specimens captured in the field must be correctly identified, in some cases to the sub-species level. Here, we generate mosquito rRNA datasets which will substantially mitigate both of these problems. We describe a strategy to assemble novel rRNA sequences from mosquito specimens and produce an unprecedented dataset of 234 full-length 28S and 18S rRNA sequences of 33 medically important species from countries with known histories of mosquito-borne virus circulation (Cambodia, the Central African Republic, Madagascar, and French Guiana). These sequences will allow both physical and computational removal of rRNA from specimens during RNA-seq protocols. We also assess the utility of rRNA sequences for molecular taxonomy and compare phylogenies constructed using rRNA sequences versus those created using the gold standard for molecular species identification of specimens-the mitochondrial cytochrome c oxidase I (COI) gene. We find that rRNA- and COI-derived phylogenetic trees are incongruent and that 28S and concatenated 28S+18S rRNA phylogenies reflect evolutionary relationships that are more aligned with contemporary mosquito systematics. This significant expansion to the current rRNA reference library for mosquitoes will improve mosquito RNA-seq metagenomics by permitting the optimization of species-specific rRNA depletion protocols for a broader range of species and streamlining species identification by rRNA sequence and phylogenetics.

Diversification of Amazonian spiny tree rats in genus Makalata (Rodentia, Echimyidae): Cryptic diversity, geographic structure and drivers of speciation

Amazonian mammal diversity is exceptionally high, yet new taxonomic discoveries continue to be made and many questions remain for understanding its diversification through time and space. Here we investigate the diversification of spiny rats in the genus Makalata, whose species are strongly associated with seasonally flooded forests, watercourses and flooded islands. We use a biogeographical approach based on a mitochondrial cytochrome b gene through divergence time estimation and reconstruction of ancestral areas and events. Our findings indicate an ancient origin of Makalata for the Guiana Shield and Eastern Amazonia as ancestral area. A first cladogenetic event led to a phylogeographic break into two broader clades of Makalata through dispersal, implying a pattern of western/Eastern Amazonian clades coinciding with the Purus Arch (middle Miocene). Most of subclades we infer originated between the late Pliocene to the early Pleistocene, with few recent exceptions in the early Pliocene through dispersal and vicariant events. The hypothesis of rivers as dispersal barriers is not corroborated for Makalata, as expected for mammalian species associated with seasonally flooded environments. We identify two key events for the expansion and diversification of Makalata species: the presence of geologically stable areas in the Guiana and Brazilian shields and the transition from lacustrine conditions in western Amazonia (Acre system) to a river system, with the establishment of the Amazon River transcontinental system and its tributaries. Our results are congruent with older geological scenarios for the Amazon basin formation (Miocene), but we do not discard the influence of recent dynamics on some speciation events and, mainly, on phylogeographic structuring processes.

From DNA barcodes to ecology: Meta-analysis of central European beetles reveal link with species ecology but also to data pattern and gaps

DNA barcoding has been used worldwide to identify biological specimens and to delimit species. It represents a cost-effective, fast, and efficient way to assess biodiversity with help of the public Barcode of Life Database (BOLD) accounting for more than 236,000 animal species and more than 10 million barcode sequences. Here, we performed a meta-analysis of available barcode data of central European Coleoptera to detect intraspecific genetic patterns among ecological groups in relation to geographic distance with the aim to investigate a possible link between infraspecific variation and species ecology. We collected information regarding feeding style, body size, as well as habitat and biotope preferences. Mantel tests and two variants of Procrustes analysis, both involving the Principal Coordinates Neighborhood Matrices (PCNM) approach, were applied on genetic and geographic distance matrices. However, significance levels were too low to further use the outcome for further trait investigation: these were in mean for all ecological guilds only 7.5, 9.4, or 15.6% for PCNM + PCA, NMDS + PCA, and Mantel test, respectively, or at best 28% for a single guild. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. Results suggest that increased numbers of species, sampling localities, and specimens for a chosen area of interest may give new insights to explore barcode data and species ecology for the scope of conservation on a larger scale. We performed a meta-analysis of available barcode data of central European beetles to detect intraspecific genetic patterns among ecological groups in relation to geographic distance, regarding feeding style, body size, as well as habitat and biotope preferences. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. However, significance levels were too low to further use the outcome for further trait investigation.

Multiple invasions, Wolbachia and human-aided transport drive the genetic variability of Aedes albopictus in the Iberian Peninsula

The Asian tiger mosquito, Aedes albopictus, is one of the most invasive species in the world. Native to the tropical forests of Southeast Asia, over the past 30 years it has rapidly spread throughout tropical and temperate regions of the world. Its dramatic expansion has resulted in public health concerns as a consequence of its vector competence for at least 16 viruses. Previous studies showed that Ae. albopictus spread has been facilitated by human-mediated transportation, but much remains unknown about how this has affected its genetic attributes. Here we examined the factors that contributed to shaping the current genetic constitution of Ae. albopictus in the Iberian Peninsula, where the species was first found in 2004, by combining population genetics and Bayesian modelling. We found that both mitochondrial and nuclear DNA markers showed a lack of genetic structure and the presence of worldwide dominant haplotypes, suggesting regular introductions from abroad. Mitochondrial DNA showed little genetic diversity compared to nuclear DNA, likely explained by infection with maternally transmitted bacteria of the genus Wolbachia. Multilevel models revealed that greater mosquito fluxes (estimated from commuting patterns and tiger mosquito population distribution) and spatial proximity between sampling sites were associated with lower nuclear genetic distance, suggesting that rapid short- and medium-distance dispersal is facilitated by humans through vehicular traffic. This study highlights the significant role of human transportation in shaping the genetic attributes of Ae. albopictus and promoting regional gene flow, and underscores the need for a territorially integrated surveillance across scales of this disease-carrying mosquito.

Genome-wide SNPs reveal novel patterns of spatial genetic structure in Aedes albopictus (Diptera Culicidae) population in China

Introduction

Since the second half of the 20th century, Aedes albopictus, a vector for more than 20 arboviruses, has spread worldwide. Aedes albopictus is the main vector of infectious diseases transmitted by Aedes mosquitoes in China, and it has caused concerns regarding public health. A comprehensive understanding of the spatial genetic structure of this vector species at a genomic level is essential for effective vector control and the prevention of vector-borne diseases.

Methods

During 2016-2018, adult female Ae. albopictus mosquitoes were collected from eight different geographical locations across China. Restriction site-associated DNA sequencing (RAD-seq) was used for high-throughput identification of single nucleotide polymorphisms (SNPs) and genotyping of the Ae. albopictus population. The spatial genetic structure was analyzed and compared to those exhibited by mitochondrial cytochrome c oxidase subunit 1 (cox1) and microsatellites in the Ae. albopictus population.

Results

A total of 9,103 genome-wide SNP loci in 101 specimens and 32 haplotypes of cox1 in 231 specimens were identified in the samples from eight locations in China. Principal component analysis revealed that samples from Lingshui and Zhanjiang were more genetically different than those from the other locations. The SNPs provided a better resolution and stronger signals for novel spatial population genetic structures than those from the cox1 data and a set of previously genotyped microsatellites. The fixation indexes from the SNP dataset showed shallow but significant genetic differentiation in the population. The Mantel test indicated a positive correlation between genetic distance and geographical distance. However, the asymmetric gene flow was detected among the populations, and it was higher from south to north and west to east than in the opposite directions.

Conclusions

The genome-wide SNPs revealed seven gene pools and fine spatial genetic structure of the Ae. albopictus population in China. The RAD-seq approach has great potential to increase our understanding of the spatial dynamics of population spread and establishment, which will help us to design new strategies for controlling vectors and mosquito-borne diseases.

Genetic variations and microbiome of the poultry red mite Dermanyssus gallinae

The poultry red mite Dermanyssus gallinae poses a significant threat to the health of hens and poultry production. A comprehensive understanding of D. gallinae is necessary to develop sustainable and efficacious control methods. Here we examined 144 D. gallinae collected from 18 poultry farms throughout the Japanese Archipelago for their genetic variations based on mitochondrial cytochrome c oxidase subunit I (COI) sequences, and microbiome variations based on amplicon sequencing of the 16S ribosomal RNA gene. According to COI sequencing, the Japanese samples were categorized into three haplogroups, which did not reflect the geographical distribution. Microbiome analyses found that the major bacteria associated with D. gallinae were Bartonella, Cardinium, Wolbachia, and Tsukamurella, with Bartonella being most predominant. Among 144 individual mites, all possessed one of the two major types of Bartonella (Bartonella sp. A), while 140 mites possessed the other type (Bartonella sp. B). The presence of the two strains of Bartonella was also confirmed by a single copy gene, rpoB. The presence of Bartonella in laid eggs suggested transovarial vertical transmission. Given that obligate blood-feeding arthropods generally require a supply of B vitamins from symbiotic bacteria, Bartonella may play an important role in mite survival. Rickettsiella, a major symbiont in European D. gallinae populations, and suggested to be an important symbiont by genomic data, was rarely found in Japanese populations. Cardinium detected from D. gallinae fell into a major clade found widely in arthropods, whereas Wolbachia detected in Japanese D. gallinae appear to be a new lineage, located at the base of Wolbachia phylogeny. Of the mitochondrial phylogeny, infection patterns of Cardinium and Wolbachia were strongly correlated, possibly suggesting one or both of the symbionts induce reproductive manipulations and increase spread in the host populations.

Polymorphisms in the hypervariable control region of the mitochondrial DNA differentiate BPH populations

The brown planthopper (BPH; Nilaparvata lugens) is one of India's most destructive pests of rice. BPH, a monophagous migratory insect, reported from all major rice-growing ecosystems of the country, is capable of traversing large distances and causing massive crop loss. A crucial step for developing viable management strategies is understanding its population dynamics. Very few reliable markers are currently available to screen BPH populations for their diversity. In the current investigation, we developed a combinatorial approach using the polymorphism present within the mitochondrial Control Region of BPH and in the nuclear genome (genomic simple sequence repeats; gSSRs) to unravel the diversity present in BPH populations collected from various rice-growing regions of India. Using two specific primer pairs, the complete Control Region (1112 to 2612 bp) was PCR amplified as two overlapping fragments, cloned and sequenced from BPH individuals representing nine different populations. Results revealed extensive polymorphism within this region due to a variable number of tandem repeats. The three selected gSSR markers also exhibited population-specific amplification patterns. Overall genetic diversity between the nine populations was high (>5%). Further, in silico double-digestion of the consensus sequences of the Control Region, with HpyCH4IV and Tsp45I restriction enzymes, revealed unique restriction fragment length polymorphisms (digital-RFLPs; dRFLPs) that differentiated all the nine BPH populations. To the best of our knowledge, this is the first report of markers developed from the Control Region of the BPH mitogenome that can differentiate populations. Eventually, such reliable and rapid marker-based identification of BPH populations will pave the way for an efficient pest management strategy.

Extremely divergent COI sequences within an amphipod species complex: A possible role for endosymbionts?

Some heritable endosymbionts can affect host mtDNA evolution in various ways. Amphipods host diverse endosymbionts, but whether their mtDNA has been influenced by these endosymbionts has yet to be considered. Here, we investigated the role of endosymbionts (microsporidians and Rickettsia) in explaining highly divergent COI sequences in Paracalliope fluviatilis species complex, the most common freshwater amphipods in New Zealand. We first contrasted phylogeographic patterns using COI, ITS, and 28S sequences. While molecular species delimitation methods based on 28S sequences supported 3-4 potential species (N, C, SA, and SB) among freshwater lineages, COI sequences supported 17-27 putative species reflecting high inter-population divergence. The deep divergence between NC and S lineages (~20%; 28S) and the substitution saturation on the 3rd codon position of COI detected even within one lineage (SA) indicate a very high level of morphological stasis. Interestingly, individuals infected and uninfected by Rickettsia comprised divergent COI lineages in one of four populations tested, suggesting a potential influence of endosymbionts in mtDNA patterns. We propose several plausible explanations for divergent COI lineages, although they would need further testing with multiple lines of evidence. Lastly, due to common morphological stasis and the presence of endosymbionts, phylogeographic patterns of amphipods based on mtDNA should be interpreted with caution.

Genomes From Historic DNA Unveil Massive Hidden Extinction and Terminal Endangerment in a Tropical Asian Songbird Radiation

Quantifying the magnitude of the global extinction crisis is important but remains challenging, as many extinction events pass unnoticed owing to our limited taxonomic knowledge of the world's organisms. The increasing rarity of many taxa renders comprehensive sampling difficult, further compounding the problem. Vertebrate lineages such as birds, which are thought to be taxonomically well understood, are therefore used as indicator groups for mapping and quantifying global extinction. To test whether extinction patterns are adequately gauged in well-studied groups, we implemented ancient-DNA protocols and retrieved whole genomes from the historic DNA of museum specimens in a widely known songbird radiation of shamas (genus Copsychus) that is assumed to be of least conservation concern. We uncovered cryptic diversity and an unexpected degree of hidden extinction and terminal endangerment. Our analyses reveal that >40% of the phylogenetic diversity of this radiation is already either extinct in the wild or nearly so, including the two genomically most distinct members of this group (omissus and nigricauda), which have so far flown under the conservation radar as they have previously been considered subspecies. Comparing the genomes of modern samples with those from roughly a century ago, we also found a significant decrease in genetic diversity and a concomitant increase in homozygosity affecting various taxa, including small-island endemics that are extinct in the wild as well as subspecies that remain widespread across the continental scale. Our application of modern genomic approaches demonstrates elevated levels of allelic and taxonomic diversity loss in a songbird clade that has not been listed as globally threatened, highlighting the importance of ongoing reassessments of extinction incidence even across well-studied animal groups. Key words: extinction, introgression, white-rumped shama, conservation.

Population structure and demographic history of the gastropod Thaisella chocolata (Duclos, 1832) from the Southeast Pacific inferred from mitochondrial DNA analyses

The present-day population structure of a species reflects the combination of oceanographic currents, life-history traits, and historical events. However, little is known about the mechanisms that have shaped the gene lineage distribution of marine species inhabiting the Southeast Pacific. Here, we provide a comprehensive phylogeographical study of a species distributed along the Southeast Pacific coastal region by analyzing the endemic gastropod Thaisella chocolata (Duclos, 1832). Sequencing of mitochondrial cytochrome c oxidase subunit 1 (CO1) and 16S rRNA revealed strikingly high haplotypic nucleotide and genetic diversity but a lack of significant population differentiation within the survey area. In addition, a star-shaped phylogeny and significantly negative Tajima's D and Fu's Fs tests of neutrality suggested historical occurrence of rapid demographic expansion. Mismatch distributions and Bayesian inference analyses also confirmed T. chocolata to have undergone two ancestral demographic expansions. Calculations suggested that these expansions began in the lower and middle Pleistocene epoch, likely due to continental shelf development and climatic conditions. These findings could help establish a genetic baseline for T. chocolata as the first step toward sustainable spatial management of this species, as well as understand this species' response to future climate change.

Linking large-scale genetic structure of three Argynnini butterfly species to geography and environment

Understanding which factors and processes are associated with genetic differentiation within and among species remains a major goal in evolutionary biology. To explore differences and similarities in genetic structure and its association with geographical and climatic factors in sympatric sister species, we conducted a large‐scale (>32° latitude and >36° longitude) comparative phylogeographical study on three Argynnini butterfly species (Speyeria aglaja, Fabriciana adippe and F. niobe) that have similar life histories, but differ in ecological generalism and dispersal abilities. Analyses of nuclear (ddRAD‐sequencing derived SNP markers) and mitochondrial (COI sequences) data revealed differences between species in genetic structure and how genetic differentiation was associated with climatic factors (temperature, solar radiation, precipitation, wind speed). Geographical proximity accounted for much of the variation in nuclear and mitochondrial structure and evolutionary relationships in F. adippe and F. niobe, but only explained the pattern observed in the nuclear data in S. aglaja, for which mitonuclear discordance was documented. In all species, Iberian and Balkan individuals formed genetic clusters, suggesting isolation in glacial refugia and limited postglacial expansion. Solar radiation and precipitation were associated with the genetic structure on a regional scale in all species, but the specific combinations of environmental and geographical factors linked to variation within species were unique, pointing to species‐specific responses to common environments. Our findings show that the species share similar colonization histories, and that the same ecological factors, such as niche breadth and dispersal capacity, covary with genetic differentiation within these species to some extent, thereby highlighting the importance of comparative phylogeographical studies in sympatric sister species.

Phylogeography and Wolbachia Infections Reveal Postglacial Recolonization Routes of the Parthenogenetic Plant Louse Cacopsylla myrtilli (W. Wagner 1947), (Hemiptera, Psylloidea)

To reveal the phylogeographic pattern of the parthenogenetic psyllid Cacopsylla myrtilli (W. Wagner 1947) (Hemiptera, Psylloidea), we sequenced a 638 bp fragment of the mitochondrial COI gene from 962 individuals. These insects originated from 46 sampling sites, which cover a significant part of the northern Palearctic distribution range of the species. The sequence data revealed 40 haplotypes, with three main (H1, H2, and H3) and 37 derived ones. The main haplotypes H1 or H2 or both were present at all sampling sites. The star-like shape of the haplotype networks indicated recent population expansion. In most cases, the derived haplotypes were specific for each country, suggesting that the main haplotypes H1 and H2 are of refugial origin, and the derived haplotypes have emerged after the postglacial recolonization process. Based on the haplotype sequences, we suggest H3 to represent the ancestral haplotype from which H1 and H2 have evolved. We suggest that the main haplotype H3 together with its derived haplotypes represents bisexual C. myrtilli, which shows a limited distribution on both sides of the border between Finland and Russia in northern Fennoscandia. The genetic diversity was the highest in Sjoa in southern Norway and also high in the White Sea region in northwest Russia. Higher diversity in Sjoa was attributed to both earlier recolonizations compared to that of the White Sea region and the absence of Wolbachia infection. We suggest that these sites were colonized from different Pleistocene refugia, i.e., from western and eastern refugia, respectively. From the White Sea region, recolonization continued eastwards to Ural Mountains and westwards to Finland and further north to Kola Peninsula. From northern Finland, recolonization continued to Finnmark, Norway, and further to Sweden and finally reached a secondary contact zone with colonizers from Norway in Central Sweden. The Caucasus and Siberian/Manchurian refugial regions have played an important role in the origin of C. myrtilli populations in Siberia and the Russian Far East.

Inheritance through the cytoplasm

Most heritable information in eukaryotic cells is encoded in the nuclear genome, with inheritance patterns following classic Mendelian segregation. Genomes residing in the cytoplasm, however, prove to be a peculiar exception to this rule. Cytoplasmic genetic elements are generally maternally inherited, although there are several exceptions where these are paternally, biparentally or doubly-uniparentally inherited. In this review, we examine the diversity and peculiarities of cytoplasmically inherited genomes, and the broad evolutionary consequences that non-Mendelian inheritance brings. We first explore the origins of vertical transmission and uniparental inheritance, before detailing the vast diversity of cytoplasmic inheritance systems across Eukaryota. We then describe the evolution of genomic organisation across lineages, how this process has been shaped by interactions with the nuclear genome and population genetics dynamics. Finally, we discuss how both nuclear and cytoplasmic genomes have evolved to co-inhabit the same host cell via one of the longest symbiotic processes, and all the opportunities for intergenomic conflict that arise due to divergence in inheritance patterns. In sum, we cannot understand the evolution of eukaryotes without understanding hereditary symbiosis.

Weak population genetic structure in Eurasian spruce bark beetle over large regional scales in Sweden

The Eurasian spruce bark beetle, Ips typographus, is a major pest, capable of killing spruce forests during large population outbreaks. Recorded dispersal distances of individual beetles are typically within hundreds of meters or a few kilometers. However, the connectivity between populations at larger distances and longer time spans and how this is affected by the habitat is less studied, despite its importance for understanding at which distances local outbreaks may spread. Previous population genetic studies in I. typographus typically used low resolution markers. Here, we use genome-wide data to assess population structure and connectivity of I. typographus in Sweden. We used 152 individuals from 19 population samples, distributed over 830 km from Strömsund (63° 46' 8″ N) in the north to Nyteboda (56° 8' 50″ N) in the south, to capture processes at a large regional scale, and a transect sampling design adjacent to a recent outbreak to capture processes at a smaller scale (76 km). Using restriction site-associated DNA sequencing (RADseq) markers capturing 1409-1997 SNPs throughout the genome, we document a weak genetic structure over the large scale, potentially indicative of high connectivity with extensive gene flow. No differentiation was detected at the smaller scale. We find indications of isolation-by-distance both for relative (F ST) and absolute divergence (Dxy). The two northernmost populations are most differentiated from the remaining populations, and diverge in parallel to the southern populations for a set of outlier loci. In conclusion, the population structure of I. typographus in Sweden is weak, suggesting a high capacity to disperse and establish outbreak populations in new territories.

Past and recent farming degrades aquatic insect genetic diversity

Recent declines in once-common species are triggering concern that an environmental crisis point has been reached. Yet, the lack of long abundance time series data for most species can make it difficult to attribute these changes to anthropogenic causes, and to separate them from normal cycles. Genetic diversity, on the other hand, is sensitive to past and recent environmental changes, and reflects a measure of a populations' potential to adapt to future stressors. Here, we consider whether patterns of genetic diversity among aquatic insects can be linked to historical and recent patterns of land use change. We collated mitochondrial cytochrome c oxidase subunit I (COI) variation for >700 aquatic insect species across the United States, where patterns of agricultural expansion and intensification have been documented since the 1800s. We found that genetic diversity was lowest in regions where cropland was historically (pre-1950) most extensive, suggesting a legacy of past environmental harm. Genetic diversity further declined where cropland has since expanded, even after accounting for climate and sampling effects. Notably though, genetic diversity also appeared to rebound where cropland has diminished. Our study suggests that genetic diversity at the community level can be a powerful tool to infer potential population declines and rebounds over longer time spans than is typically possible with ecological data. For the aquatic insects that we considered, patterns of land use many decades ago appear to have left long-lasting damage to genetic diversity that could threaten evolutionary responses to rapid global change.

Wolbachia Infection in Native Populations of the Invasive Tawny Crazy Ant Nylanderia fulva

Antagonistic interactions can affect population growth and dispersal of an invasive species. Wolbachia are intracellular endosymbiont bacteria that infect arthropod and nematode hosts and are able to manipulate reproduction, which in some cases leads to cocladogenesis. Moreover, the presence of the strictly maternally transferred Wolbachia in a population can indirectly induce selective sweeps on the hosts' mitochondria. Ants have a Wolbachia infection rate of about 34%, which makes phylogenetic studies using mitochondrial markers vulnerable of being confounded by the effect of the endosymbiont. Nylanderia fulva is an invasive ant native to South America, considered a pest in the United States. Its distribution and biology are poorly known in its native range, and the taxonomic identity of this and its closely related species, Nylanderia pubens, has only recently been understood with the aid of molecular phylogenies. Aiming at estimating robust phylogenetic relationships of N. fulva in its native range, we investigated the presence and pattern of Wolbachia infection in populations of N. fulva from Argentina, part of its native range, to account for its possible effect on the host population structure. Using the ftsZ gene, 30 nests of N. fulva and four from sympatric Nylanderia species were screened for the presence of Wolbachia. We sequenced the MLST genes, the highly variable gene wsp, as well as glyQ, a novel target gene for which new primers were designed. Phylogeny of the ants was estimated using mtDNA (COI). We found supergroup A Wolbachia strains infecting 73% of N. fulva nests and two nests of Nylanderia sp. 1. Wolbachia phylogenetic tree inferred with MLST genes is partially congruent with the host phylogeny topology, with the exception of a lineage of strains shared by ants from different N. fulva clades. Furthermore, by comparing with Wolbachia sequences infecting other ants, we found that the strains infecting different N. fulva clades are not monophyletic. Our findings suggest there are three recent independent horizontally transmitted Wolbachia infections in N. fulva, and we found no evidence of influence of Wolbachia in the host mtDNA based phylogeny.

Unveiling the Mycodrosophila projectans (Diptera, Drosophilidae) species complex: Insights into the evolution of three Neotropical cryptic and syntopic species

The Zygothrica genus group has been shown to be speciose, with a high number of cryptic species. DNA barcoding approaches have been a valuable tool to uncover cryptic diversity in this lineage, as recently suggested for the Neotropical Mycodrosophila projectans complex, which seems to comprise at least three different species. The aim of this study was to confirm the subdivision of the M. projectans complex while shedding some light on the patterns and processes related to its diversification. In this sense, the use of single and multi-locus datasets under phylogenetic, distance, coalescence, and diagnostic nucleotide approaches confirmed the presence of at least three species under the general morphotype previously described as M. projectans. Only a few subtle morphological differences were found for the three species in terms of aedeagus morphology and abdominal color patterns. Ecologically, sympatry and syntopy seem to be recurrent for these three cryptic species, which present widely overlapping niches, implying niche conservatism. This morphological and ecological similarity has persisted though cladogenesis within the complex, which dates back to the Miocene, providing an interesting example of morphological conservation despite ancient divergence. These results, in addition to contrasting patterns of past demographic fluctuations, allowed us to hypothesize patterns of allopatric or parapatric diversification with secondary contact in Southern Brazil. Nevertheless, genetic diversity was generally high within species, suggesting that migration may encompass an adaptive response to the restrictions imposed by the ephemerality of resources.

A population genetic analysis of the Critically Endangered Madagascar big-headed turtle, Erymnochelys madagascariensis across captive and wild populations

Erymnochelys madagascariensis is a Critically Endangered turtle endemic to Madagascar. Anthropogenic activity has depleted the wild population by 70% in the last century, and effective conservation management is essential to ensuring its persistence. Captive breeding was implemented to augment depleted populations in the southern part of Ankarafantsika National Park (ANP), when no genetic data were available for E. madagascariensis. It is unknown how much of the natural population's diversity is encapsulated in captivity. We used eight microsatellite loci and fragments of two mitochondrial genes to identify the genetic structure of E. madagascariensis in the wild. Captive bred turtles were compared with wild populations in order to assess the representativeness of this ex situ conservation strategy for ANP. Six microsatellite clusters, ten cytochrome b, and nine COI haplotypes were identified across wild populations, with high genetic divergence found between populations in two groups of watersheds. Captive bred individuals represent three out of six sampled microsatellite clusters found in the wild and just one mitochondrial haplotype, possibly due to genetic drift. To improve genetic representation, the strategy of frequent interchange between captive and wild breeders within ANP should be revitalised and, as originally planned, hatchlings or juveniles should not be released beyond ANP.

Universal Mitochondrial Multi-Locus Sequence Analysis (mtMLSA) to Characterise Populations of Unanticipated Plant Pest Biosecurity Detections

Simple Summary

Agricultural and environmental sustainability requires effective biosecurity responses that prevent the establishment or spread of exotic insect pests. Understanding where new detections may have come from or if recurrent detections are connected contributes to this. Suitable population genetic markers use relatively rapidly evolving gene regions which render the PCR method species-specific at best. Because resource limitations mean these are pre-emptively developed for the highest risk species, populations of other exotic pests are unable to be characterised at the time. Here we have developed a generic method that is useful across species within the same taxonomic Order, including where there is little or no prior knowledge of their gene sequences. Markers are formed by concomitant sequencing of four gene regions. Sequence concatenation was shown to retrieve higher resolution signatures than standard DNA barcoding. The method is encouragingly universal, as illustrated across species in ten fly and 11 moth superfamilies. Although as-yet untested in a biosecurity situation, this relatively low-tech, off-the-shelf method makes a proactive contribution to the toolbox of quarantine agencies at the time of detection without the need for impromptu species-specific research and development.

Abstract

Biosecurity responses to post-border exotic pest detections are more effective with knowledge of where the species may have originated from or if recurrent detections are connected. Population genetic markers for this are typically species-specific and not available in advance for any but the highest risk species, leaving other less anticipated species difficult to assess at the time. Here, new degenerate PCR primer sets are designed for within the Lepidoptera and Diptera for the 3′ COI, ND3, ND6, and 3′ plus 5′ 16S gene regions. These are shown to be universal at the ordinal level amongst species of 14 and 15 families across 10 and 11 dipteran and lepidopteran superfamilies, respectively. Sequencing the ND3 amplicons as an example of all the loci confirmed detection of population-level variation. This supported finding multiple population haplotypes from the publicly available sequences. Concatenation of the sequences also confirmed that higher population resolution is achieved than for the individual genes. Although as-yet untested in a biosecurity situation, this method is a relatively simple, off-the-shelf means to characterise populations. This makes a proactive contribution to the toolbox of quarantine agencies at the time of detection without the need for unprepared species-specific research and development.

Integrative taxonomy reveals cryptic diversity in North American Lasius ants, and an overlooked introduced species

Biological invasions are a grave threat to ecosystems. The black garden ant (Lasius niger) is a pest species in Europe. Current literature states that L. niger occupies a disjunct native distribution in the Holarctic, however, based on recent work, we re-evaluate this distribution. The native range of L. niger is reconsidered based on phylogenetic relationships (nine mitochondrial and nuclear markers, 5670 bp), DNA-barcoding (98 Holarctic specimens), morphometry (88 Holarctic specimens, 19 different measurements) and subjective assessment of phenotype. The potential spread of this species is estimated using ecological niche modeling. Lasius niger is more closely related to other Palearctic species than to the Nearctic ants known under this name. The latter are described as a distinct species, L. ponderosae sp. nov. However, DNA-barcoding discovered established populations of L. niger in metropolitan areas in Canada (Vancouver and Halifax). We describe a morphometrical method to delineate L. ponderosae sp. nov. and L. niger. MtDNA diversity and divergence is high within L. ponderosae sp. nov., but low within L. niger. More than 1,000,000 km² are suitable as a habitat for L. niger in North America. This case emphasizes the critical role of integrative taxonomy to detect cryptic species and identify potential biological invasions in their nascent stages.

Relationship of stranded cetaceans in Thai territorial waters to global populations: Mitochondrial DNA diversity of Cuvier's beaked whale, Indo Pacific finless porpoise, pygmy sperm whale, and dwarf sperm whale

Cetaceans inhabit oceans throughout the world. Four specific odontocetes, namely Cuvier's beaked whale (Ziphius cavirostris), Indo Pacific finless porpoise (Neophocaena phocaenoides), pygmy sperm whale (Kogia breviceps), and dwarf sperm whale (Kogia sima), have occasionally been found stranded along Thailand's coastal waters (the Andaman Sea and the Gulf of Thailand). Although shared haplotypes of each species for many locations have been found, and some species have revealed genetic structure through haplotype networks, cetaceans in Thai waters have never been investigated in terms of comparing haplotypes to those that have existed before. Herein, we have illustrated the matrilineally phylogeographic relationships among worldwide populations through Bayesian Phylogenetic tree computations using Markov Chain Monte Carlo (MCMC) and Median-Joining Networks (MJNs). Unique haplotypes of the control region mitochondrial DNA of Thai odontocetes were found for all species. Moreover, a high degree of worldwide haplotype diversity (hd) above 0.8 among the four species was detected, while the lowest degree of nucleotide diversity (π) was observed in the Indo Pacific finless porpoise (1.12% ± 0.184%). An expansion of the effective female population size worldwide of three odontocete species was detected using Bayesian Skyline Reconstruction, but this did not include the Indo Pacific finless porpoise. Because Thai seas are located within the Indo Polynesian province, where this biodiversity hotspot exists, we speculate that these odontocetes may also inhabit specific habitats within the Malay Peninsula and Thailand's territorial waters. Therefore, closer attention and monitoring of these cetacean populations will be necessary for future conservation efforts.

Hybridization fuelled diversification in Spialia butterflies

The importance of hybridization and introgression is well documented in the evolution of plants but, in insects, their role is not fully understood. Given the fact that insects are the most diverse group of organisms, assessing the impact of reticulation events on their evolution may be key to comprehend the emergence of such remarkable diversity. Here, we used an insect model, the Spialia butterflies, to gather genomic evidence of hybridization as a promoter of novel diversity. By using double‐digest RADseq (ddRADseq), we explored the phylogenetic relationships between Spialia orbifer, S. rosae and S. sertorius, and documented two independent events of interspecific gene flow. Our data support that the Iberian endemism S. rosae probably received genetic material from S. orbifer in both mitochondrial and nuclear DNA, which could have contributed to a shift in the ecological preferences of S. rosae. We also show that admixture between S. sertorius and S. orbifer probably occurred in Italy. As a result, the admixed Sicilian populations of S. orbifer are differentiated from the rest of populations both genetically and morphologically, and display signatures of reproductive character displacement in the male genitalia. Additionally, our analyses indicated that genetic material from S. orbifer is present in S. sertorius along the Italian Peninsula. Our findings add to the view that hybridization is a pervasive phenomenon in nature and in butterflies in particular, with important consequences for evolution due to the emergence of novel phenotypes.

Quality over quantity: unraveling the contributions to cytoplasmic incompatibility caused by two coinfecting Cardinium symbionts

Cytoplasmic incompatibility (CI) is a common form of reproductive sabotage caused by maternally inherited bacterial symbionts of arthropods. CI is a two-step manipulation: first, the symbiont modifies sperm in male hosts which results in the death of fertilized, uninfected embryos. Second, when females are infected with a compatible strain, the symbiont reverses sperm modification in the fertilized egg, allowing offspring of infected females to survive and spread the symbiont to high frequencies in a population. Although CI plays a role in arthropod evolution, the mechanism of CI is unknown for many symbionts. Cardinium hertigii is a common CI-inducing symbiont of arthropods, including parasitoid wasps like Encarsia partenopea. This wasp harbors two Cardinium strains, cEina2 and cEina3, and exhibits strong CI. The strains infect wasps at different densities, with the cEina3 present at a lower density than cEina2, and it was previously not known which strain caused CI. By differentially curing wasps of cEina3, we found that this low-density symbiont is responsible for CI and modifies males during their pupal stage. cEina2 does not modify host reproduction and may spread by 'hitchhiking' with cEina3 CI or by conferring an unknown benefit. The cEina3 strain also shows a unique localization pattern in male reproductive tissues. Instead of infecting sperm like other CI-inducing symbionts, cEina3 cells are found in somatic cells at the testis base and around the seminal vesicle. This may allow the low-density cEina3 to efficiently modify host males and suggests that cEina3 uses a different modification strategy than sperm-infecting CI symbionts.

Genomics and ecological modelling clarify species integrity in a confusing group of butterflies

Recent advances in both genomics and ecological modelling present new, multidisciplinary opportunities for resolving species boundaries and understanding the mechanisms that maintain their integrity in regions of contact. Here, we use a combination of high-throughput DNA sequencing and ecological niche modelling to resolve species boundaries and niche divergence within the Speyeria atlantis-hesperis (Lepidoptera: Nymphalidae) complex, a confusing group of North American butterflies. This complex is notorious for its muddled species delimitations, morphological ambiguity, and extensive mito-nuclear discordance. Our admixture and multispecies coalescent-based analyses of single nucleotide polymorphisms identified substantial divergences between S. atlantis and S. hesperis in areas of contact, as well as between distinct northern and southern lineages within S. hesperis. Our results also provide evidence of past introgression relating to another species, S. zerene, which previous work has shown to be more distantly related to the S. atlantis-hesperis complex. We then used ecological models to predict habitat suitability for each of the three recovered genomic lineages in the S. atlantis-hesperis complex and assess their pairwise niche divergence. These analyses resolved that these three lineages are significantly diverged in their respective niches and are not separated by discontinuities in suitable habitat that might present barriers to gene flow. We therefore infer that ecologically-mediated selection resulting in disparate habitat associations is a principal mechanism reinforcing their genomic integrity. Overall, our results unambiguously support significant evolutionary and ecological divergence between the northern and southern lineages of S. hesperis, sufficient to recognize the southern evolutionary lineage as a distinct species, called S. nausicaa based on taxonomic priority.

Surprisingly high genetic divergence of the mitochondrial DNA barcode fragment (COI) within Central European woodlice species (Crustacea, Isopoda, Oniscidea)

DNA barcoding has become the most popular approach for species identification in recent years. As part of the German Barcode of Life project, the first DNA barcode library for terrestrial and freshwater isopods from Germany is presented. The analyzed barcode library included 38 terrestrial (78% of the documented species of Germany) and five freshwater (63%) species. A total of 513 new barcodes was generated and 518 DNA barcodes were analyzed. This analysis revealed surprisingly high intraspecific genetic distances for numerous species, with a maximum of 29.4% for Platyarthrus hoffmannseggii Brandt, 1833. The number of BINs per species ranged from one (32 species, 68%) to a maximum of six for Trachelipus rathkii (Brandt, 1833). In spite of such high intraspecific variability, interspecific distances with values between 12.6% and 29.8% allowed a valid species assignment of all analyzed isopods. The observed high intraspecific distances presumably result from phylogeographic events, Wolbachia infections, atypical mitochondrial DNAs, heteroplasmy, or various combinations of these factors. Our study represents the first step in generating an extensive reference library of DNA barcodes for terrestrial and freshwater isopods for future molecular biodiversity assessment studies.

Substantial rearrangements, single nucleotide frameshift deletion and low diversity in mitogenome of Wolbachia-infected strepsipteran endoparasitoid in comparison to its tephritid hosts

Insect mitogenome organisation is highly conserved, yet, some insects, especially with parasitic life cycles, have rearranged mitogenomes. Furthermore, intraspecific mitochondrial diversity can be reduced by fitness-affecting bacterial endosymbionts like Wolbachia due to their maternal coinheritance with mitochondria. We have sequenced mitogenomes of the Wolbachia-infected endoparasitoid Dipterophagus daci (Strepsiptera: Halictophagidae) and four of its 22 known tephritid fruit fly host species using total genomic extracts of parasitised flies collected across > 700 km in Australia. This halictophagid mitogenome revealed extensive rearrangements relative to the four fly mitogenomes which exhibited the ancestral insect mitogenome pattern. Compared to the only four available other strepsipteran mitogenomes, the D. daci mitogenome had additional transpositions of one rRNA and two tRNA genes, and a single nucleotide frameshift deletion in nad5 requiring translational frameshifting or, alternatively, resulting in a large protein truncation. Dipterophagus daci displays an almost completely endoparasitic life cycle when compared to Strepsiptera that have maintained the ancestral state of free-living adults. Our results support the hypothesis that the transition to extreme endoparasitism evolved together with increased levels of mitogenome changes. Furthermore, intraspecific mitogenome diversity was substantially smaller in D. daci than the parasitised flies suggesting Wolbachia reduced mitochondrial diversity because of a role in D. daci fitness.