Mitogenome phylogeographic analysis of a planktonic crustacean

Genome-Wide Allele Frequency Changes Reveal That Dynamic Metapopulations Evolve Differently

Two important characteristics of metapopulations are extinction-(re)colonization dynamics and gene flow between subpopulations. These processes can cause strong shifts in genome-wide allele frequencies that are generally not observed in "classical" (large, stable, and panmictic) populations. Subpopulations founded by one or a few individuals, the so-called propagule model, are initially expected to show intermediate allele frequencies at polymorphic sites until natural selection and genetic drift drive allele frequencies toward a mutation-selection-drift equilibrium characterized by a negative exponential-like distribution of the site frequency spectrum. We followed changes in site frequency spectrum distribution in a natural metapopulation of the cyclically parthenogenetic pond-dwelling microcrustacean Daphnia magna using biannual pool-seq samples collected over a 5-yr period from 118 ponds occupied by subpopulations of known age. As expected under the propagule model, site frequency spectra in newly founded subpopulations trended toward intermediate allele frequencies and shifted toward right-skewed distributions as the populations aged. Immigration and subsequent hybrid vigor altered this dynamic. We show that the analysis of site frequency spectrum dynamics is a powerful approach to understand evolution in metapopulations. It allowed us to disentangle evolutionary processes occurring in a natural metapopulation, where many subpopulations evolve in parallel. Thereby, stochastic processes like founder and immigration events lead to a pattern of subpopulation divergence, while genetic drift leads to converging site frequency spectrum distributions in the persisting subpopulations. The observed processes are well explained by the propagule model and highlight that metapopulations evolve differently from classical populations.

Long-term balancing selection for pathogen resistance maintains trans-species polymorphisms in a planktonic crustacean

Balancing selection is an evolutionary process that maintains genetic polymorphisms at selected loci and strongly reduces the likelihood of allele fixation. When allelic polymorphisms that predate speciation events are maintained independently in the resulting lineages, a pattern of trans-species polymorphisms may occur. Trans-species polymorphisms have been identified for loci related to mating systems and the MHC, but they are generally rare. Trans-species polymorphisms in disease loci are believed to be a consequence of long-term host-parasite coevolution by balancing selection, the so-called Red Queen dynamics. Here we scan the genomes of three crustaceans with a divergence of over 15 million years and identify 11 genes containing identical-by-descent trans-species polymorphisms with the same polymorphisms in all three species. Four of these genes display molecular footprints of balancing selection and have a function related to immunity. Three of them are located in or close to loci involved in resistance to a virulent bacterial pathogen, Pasteuria, with which the Daphnia host is known to coevolve. This provides rare evidence of trans-species polymorphisms for loci known to be functionally relevant in interactions with a widespread and highly specific parasite. These findings support the theory that specific antagonistic coevolution is able to maintain genetic diversity over millions of years.

Short lifespan is one's fate, long lifespan is one's achievement: lessons from Daphnia

Studies of longevity rely on baseline life expectancy of reference genotypes measured in standardized conditions. Variation among labs, protocols, and genotypes makes longevity intervention studies difficult to compare. Furthermore, extending lifespan under suboptimal conditions or that of a short-lived genotype may be of a lesser theoretical and translational value than extending the maximal possible lifespan. Daphnia is becoming a model organism of choice for longevity research complementing data obtained on traditional models. In this study, we report longevity of several genotypes of a long-lived species D. magna under a variety of protocols, aiming to document the highest lifespan, factors reducing it, and parameters that change with age and correlate with longevity. Combining longevity data from 25 experiments across two labs, we report a strong intraspecific variation, moderate effects of group size and medium composition, and strong genotype-by-environment interactions with respect to food level. Specifically, short-lived genotypes show no caloric restriction (CR) effect, while long-lived ones expand their lifespan even further under CR. We find that the CR non-responsive clones show little correlation between longevity and two measures of lipid peroxidation. In contrast, the long-lived, CR-responsive clones show a positive correlation between longevity and lipid hydroperoxide abundance, and a negative correlation with MDA concentration. This indicates differences among genotypes in age-related accumulation and detoxification of LPO products and their effects on longevity. Our observations support the hypothesis that a long lifespan can be affected by CR and levels of oxidative damage, while genetically determined short lifespan remains short regardless.

The complete mitochondrial genome of water flea Ceriodaphnia dubia (Crustacea: Cladocera) NIES strain

Water flea Ceriodaphnia dubia has been widely used for risk assessments of chemicals and environmental contamination. In this study, the complete mitochondrial genome (mitogenome) of this species NIES strain was determined using short-read high throughput and long-read sequencing technologies. The mitogenome of C. dubia was 15,170 bp in length and consisted of 13 protein-coding genes (PCGs), 2 ribosomal RNAs (rRNAs), and 22 transfer RNAs (tRNAs). The gene order was identical to the pattern conserved across crustaceans. The complete mitogenome of the NIES strain will serve as genetical reference in ecological risk assessments in Japan, as well as resources for future phylogenetical studies using cladocerans.

Population genetic analysis of the microsporidium Ordospora colligata reveals the role of natural selection and phylogeography on its extremely compact and reduced genome

The determinants of variation in a species' genome-wide nucleotide diversity include historical, environmental, and stochastic aspects. This diversity can inform us about the species' past and present evolutionary dynamics. In parasites, the mode of transmission and the interactions with the host might supersede the effects of these aspects in shaping parasite genomic diversity. We used genomic samples from 10 populations of the microsporidian parasite Ordospora colligata to investigate present genomic diversity and how it was shaped by evolutionary processes, specifically, the role of phylogeography, co-phylogeography (with the host), natural selection, and transmission mode. Although very closely related microsporidia cause diseases in humans, O. colligata is specific to the freshwater crustacean Daphnia magna and has one of the smallest known eukaryotic genomes. We found an overlapping phylogeography between O. colligata and its host highlighting the long-term, intimate relationship between them. The observed geographic distribution reflects previous findings that O. colligata exhibits adaptations to colder habitats, which differentiates it from other microsporidian gut parasites of D. magna predominantly found in warmer areas. The co-phylogeography allowed us to calibrate the O. colligata phylogeny and thus estimate its mutation rate. We identified several genetic regions under potential selection. Our whole-genome study provides insights into the evolution of one of the most reduced eukaryotic genomes and shows how different processes shape genomic diversity of an obligate parasite.

Climate Warming Since the Holocene Accelerates West-East Communication for the Eurasian Temperate Water Strider Species Aquarius paludum

Holocene climate warming has dramatically altered biological diversity and distributions. Recent human-induced emissions of greenhouse gases will exacerbate global warming and thus induce threats to cold-adapted taxa. However, the impacts of this major climate change on transcontinental temperate species are still poorly understood. Here, we generated extensive genomic datasets for a water strider, Aquarius paludum, which was sampled across its entire distribution in Eurasia and used these datasets in combination with ecological niche modeling (ENM) to elucidate the influence of the Holocene and future climate warming on its population structure and demographic history. We found that A. paludum consisted of two phylogeographic lineages that diverged in the middle Pleistocene, which resulted in a “west–east component” genetic pattern that was probably triggered by Central Asia-Mongoxin aridification and Pleistocene glaciations. The diverged western and eastern lineages had a second contact in the Holocene, which shaped a temporary hybrid zone located at the boundary of the arid–semiarid regions of China. Future predictions detected a potentially novel northern corridor to connect the western and eastern populations, indicating west–east gene flow would possibly continue to intensify under future warming climate conditions. Further integrating phylogeographic and ENM analyses of multiple Eurasian temperate taxa based on published studies reinforced our findings on the “west–east component” genetic pattern and the predicted future northern corridor for A. paludum. Our study provided a detailed paradigm from a phylogeographic perspective of how transcontinental temperate species differ from cold-adapted taxa in their response to climate warming.

Complete mitochondrial genome of Ovalona pulchella (Branchiopoda, Cladocera) as the first representative in the family Chydoridae: Gene rearrangements and phylogenetic analysis of Cladocera

Chydoridae are phytophilic-benthic microcrustaceans that make up a significant proportion of species diversity and play an important role in the littoral zone of freshwater ecosystems worldwide. Here, we provide the complete mitochondrial genome of Ovalona pulchella (King, 1853), determined by next-generation sequencing. The entire mitochondrial genome is 15,362 bp in length; this is the first sequenced mitochondrial genome in the family Chydoridae. The base composition and codon usage were typical of Cladocera species. The mitochondrial gene arrangement (37 genes) was not consistent with that of other Branchiopoda. Both maximum likelihood and Bayesian analyses supported each suborder and family of Branchiopoda as monophyletic groups. The relationships among the families were as follows: [(Leptestheriidae + Limnadiidae) + (Sididae + (Bosminidae + (Chydoridae + Daphniidae)))] + Triopsidae. The newly sequenced O. pulchella was most closely related to the family Daphniidae. The complete mitochondrial genome of O. pulchella also provides valuable molecular information for further analysis of the phylogeny of the Chydoridae and the taxonomic status of the Branchiopoda.

Population-genomic analysis identifies a low rate of global adaptive fixation in the proteins of the cyclical parthenogen Daphnia magna

Daphnia are well-established ecological and evolutionary models, and the interaction between D. magna and its microparasites is widely considered a paragon of the host-parasite coevolutionary process. Like other well-studied arthropods such as Drosophila melanogaster and Anopheles gambiae, D. magna is a small, widespread, and abundant species that is therefore expected to display a large long-term population size and high rates of adaptive protein evolution. However, unlike these other species, D. magna is cyclically asexual and lives in a highly structured environment (ponds and lakes) with moderate levels of dispersal, both of which are predicted to impact upon long-term effective population size and adaptive protein evolution. To investigate patterns of adaptive protein fixation, we produced the complete coding genomes of 36 D. magna clones sampled from across the European range (Western Palaearctic), along with draft sequences for the close relatives D. similis and D. lumholtzi, used as outgroups. We analyzed genome-wide patterns of adaptive fixation, with a particular focus on genes that have an a priori expectation of high rates, such as those likely to mediate immune responses, RNA interference against viruses and transposable elements, and those with a strongly male-biased expression pattern. We find that, as expected, D. magna displays high levels of diversity and that this is highly structured among populations. However, compared with Drosophila, we find that D. magna proteins appear to have a high proportion of weakly deleterious variants and do not show evidence of pervasive adaptive fixation across its entire range. This is true of the genome as a whole, and also of putative ‘arms race’ genes that often show elevated levels of adaptive substitution in other species. In addition to the likely impact of extensive, and previously documented, local adaptation, we speculate that these findings may reflect reduced efficacy of selection associated with cyclical asexual reproduction.

Complete mitochondrial genomes of three fairy shrimps from snowmelt pools in Japan

Background

Fairy shrimps belong to order Anostraca, class Branchiopoda, subphylum Crustacea, and phylum Arthropoda. Three fairy shrimp species (Eubranchipus uchidai, E. asanumai, and E. hatanakai) that inhabit snowmelt pools are currently known in Japan. Whole mitochondrial genomes are useful genetic information for conducting phylogenetic analyses. Mitochondrial genome sequences for Branchiopoda members are gradually being collated.

Results

Six whole mitochondrial genomes from the three Eubranchipus species are presented here. Eubranchipus species share the anostracan pattern of gene arrangement in their mitochondrial genomes. The mitochondrial genomes of the Eubranchipus species have a higher GC content than those of other anostracans. Accelerated substitution rates in the lineage of Eubranchipus species were observed.

Conclusion

This study is the first to obtain whole mitochondrial genomes for Far Eastern Eubranchipus species. We show that the nucleotide sequences of cytochrome oxidase subunit I and the 16S ribosomal RNA of E. asanumai presented in a previous study were nuclear mitochondrial DNA segments. Higher GC contents and accelerated substitution rates are specific characteristics of the mitochondrial genomes of Far Eastern Eubranchipus. The results will be useful for further investigations of the evolution of Anostraca as well as Branchiopoda.

Mitogenomics of Cladocera (Branchiopoda): Marked gene order rearrangements and independent predation roots

Cladocera (Crustacea: Branchiopoda) is a key group of invertebrates. Despite a long history of phylogenetic research, relationships within this group remain disputed. We here provide new insights based on 15 new mitochondrial genomes obtained from high-throughput sequencing (HTS) and 40 mitogenomes extracted from published HTS datasets. Together with 25 mitogenomes from GenBank, we generated a matrix of 80 mitogenomes, 44 of them belonging to Cladocera. We also obtained a matrix with 168 nuclear orthologous genes to further assess the phylogenetic result from mitogenomes based on published data and one new HTS data ofLeptodora. Maximum likelihood and Bayesian phylogenetic analyses recovered all Branchiopoda orders as monophyletic and supported a sister-group relationship between Anomopoda and Onychopoda, making the taxon Gymnomera paraphyletic and supporting an independent origin of predatory Haplopoda and Onychopoda. The nuclear phylogeny and topological tests also support Gymnomera as paraphyletic, and the nuclear phylogeny strongly supports a sister-group relationship between Ctenopoda and Haplopoda. We provide a fossil-calibrated time tree, congruent with a Carboniferous origin for Cladocera and a subsequent diversification of the crown group of Anomopoda, Onychopoda, and Ctenopoda, at least in the Triassic. Despite their long evolutionary history, non-Cladoceran Branchiopoda exhibited high mitogenome structural stability. On the other hand, 21 out of 24 gene rearrangements occurred within the relatively younger Cladocera. We found the differential base compositional skewness patterns between Daphnia s.s. and Ctenodaphnia, which might be related to the divergence between these taxa. We also provide evidence to support the recent finding that Spinicaudata possesses mitogenomes with inversed compositional skewness without gene rearrangement. Such a pattern has only been reported in Spinicaudata.

Balancing selection for pathogen resistance reveals an intercontinental signature of Red Queen coevolution

The link between long-term host–parasite coevolution and genetic diversity is key to understanding genetic epidemiology and the evolution of resistance. The model of Red Queen host–parasite coevolution posits that high genetic diversity is maintained when rare host resistance variants have a selective advantage, which is believed to be the mechanistic basis for the extraordinarily high levels of diversity at disease-related genes such as the major histocompatibility complex in jawed vertebrates and R-genes in plants. The parasites that drive long-term coevolution are, however, often elusive. Here we present evidence for long-term balancing selection at the phenotypic (variation in resistance) and genomic (resistance locus) level in a particular host–parasite system: the planktonic crustacean Daphnia magna and the bacterium Pasteuria ramosa. The host shows widespread polymorphisms for pathogen resistance regardless of geographic distance, even though there is a clear genome-wide pattern of isolation by distance at other sites. In the genomic region of a previously identified resistance supergene, we observed consistent molecular signals of balancing selection, including higher genetic diversity, older coalescence times, and lower differentiation between populations, which set this region apart from the rest of the genome. We propose that specific long-term coevolution by negative-frequency-dependent selection drives this elevated diversity at the host's resistance loci on an intercontinental scale and provide an example of a direct link between the host’s resistance to a virulent pathogen and the large-scale diversity of its underlying genes.

Characterization and Phylogenetic Analysis of the Mitochondrial Genome Sequence of Nisia fuliginosa (Hemiptera: Fulgoroidea: Meenoplidae)

We explored characterization of the mitochondrial genome (mitogenome or mtGenome) and phylogenetic analysis between 32 Fulgoroid species by sequencing and analyzing the mitogenome of Nisia fuliginosa Yang and Hu, 1985 (Hemiptera: Fulgoroidea: Meenoplidae), thereby making it the first determined mitogenome from the family Meenoplidae. The mitogenome was found to be 15,754 bp in length and contained 13 protein-coding genes (PCGs), 22 tRNA genes, two ribosomal RNA genes (rRNAs), and a control region. All PCGs started with typical ATN codons, except for nad1, which used GTG as the start codon. Canonical TAA termination codons were found in 10 PCGs and the remaining three genes (cox2, nad6, and nad1) had incomplete stop codons T. All tRNAs could fold into typical cloverleaf secondary structures, with the exception of trnC, trnV, and trnS1. Additionally, we compared the AT and GC skews of 13 PCGs of 32 Fulgoroidea mitogenomes, on the L-strand, the AT and GC skews were negative and positive, respectively. However, on the H-strand, the AT skew could be positive or negative and the GC skew was always negative. Phylogenetic results showed that the eight families of Fulgoroidea were divided into two large groups. Delphacidae formed a monophyletic group sister to a clade comprising Meenoplidae and other six families (Fulgoridae, Ricaniidae, Flatidae, Issidae, Caliscelidae, and Achilidae). Meenoplidae was located near the clade of Delphacidae, and Fulgoridae was located near the clade of Meenoplidae. Furthermore, Caliscelidae, Issidae, Ricaniidae, and Flatidae are closely related and they collectively formed a sister group to Achilidae.

Does geographic variation in thermal tolerance in Daphnia represent trade-offs or conditional neutrality?

Geographic variation in thermal tolerance in Daphnia seems to represent genetic load at the loci specifically responsible for heat tolerance resulting from conditional neutrality. We see no evidence of trade-offs between fitness-related traits at 25 °C vs. 10 °C or between two algal diets across Daphnia magna clones from a variety of locations representing the opposite ends of the distribution of long-term heat tolerance. Likewise, we found no evidence of within-environment trade-offs between heat tolerance and fitness-related traits in any of the environments. Neither short-term and long-term heat tolerance shows any consistent relationship with lipid fluorescence polarization and lipid peroxidation across clones or environments. Pervasive positive correlations between fitness-related traits indicate differences in genetic load rather than trade-off based local adaptation or thermal specialization. For heat tolerance such differences may be caused by either relaxation of stabilizing selection due to lower exposure to high temperature extremes, i.e., conditional neutrality, or by small effective population size followed by the recent range expansion.

Global phylogeography and invasion history of the spotted lanternfly revealed by mitochondrial phylogenomics

Biological invasion has been a serious global threat due to increasing international trade and population movements. Tracking the source and route of invasive species and evaluating the genetic differences in their native regions have great significance for the effective monitoring and management, and further resolving the invasive mechanism. The spotted lanternfly Lycorma delicatula is native to China and invaded South Korea, Japan, and the United States during the last decade, causing severe damages to the fruits and timber industries. However, its global phylogeographic pattern and invasion history are not clearly understood. We applied high-throughput sequencing to obtain 392 whole mitochondrial genome sequences from four countries to ascertain the origin, dispersal, and invasion history of the spotted lanternfly. Phylogenomic analyses revealed that the spotted lanternfly originated from southwestern China, diverged into six phylogeographic lineages, and experienced northward expansion across the Yangtze River in the late Pleistocene. South Korea populations were derived from multiple invasions from eastern China and Japan with two different genetic sources of northwestern (Loess Plateau) and eastern (East Plain) lineages in China, whereas the each of Japan and the United States had only one. The United States populations originated through single invasive event from South Korea, which served as a bridgehead of invasion. The environmental conditions, especially the distribution of host Ailanthus trees, and adaptability possibly account for the rapid spread of the spotted lanternfly in the native and introduced regions.

Genomic characterization of selfing in the cyclic parthenogen Daphnia magna

Inbreeding refers to the fusion of related individuals' gametes, with self-fertilization (selfing) being an extreme form of inbreeding-involving gametes produced by the same individual. Selfing is expected to reduce heterozygosity by an average of 50% in one generation; however, little is known about the empirical variation on a genome level surrounding this figure and the factors that affect variation. We selfed genotypes of the cyclic parthenogen Daphnia magna and analysed whole genomes of mothers and selfed offspring, observing the predicted 50% heterozygosity reduction on average. We also saw substantial variation around this value and significant differences among mother-offspring pairs. Crossover analysis confirmed the known trend of recombination occurring more often towards the telomeres. This effect was shown, through simulations, to increase the variance of heterozygosity reduction compared to when a uniform distribution of crossovers was used. Similarly, we simulated inbred line production after several generations of selfing and we observed higher variance in achieved homozygosity when we consider a higher recombination rate towards the telomeres. Our empirical and simulation study highlights that the expected mean values of heterozygosity reduction show remarkable variation, which can help understand, for example, differences among inbred individuals.

Microsporidia with Vertical Transmission Were Likely Shaped by Nonadaptive Processes

Microsporidia have the leanest genomes among eukaryotes, and their physiological and genomic simplicity has been attributed to their intracellular, obligate parasitic life-style. However, not all microsporidia genomes are small or lean, with the largest dwarfing the smallest ones by at least an order of magnitude. To better understand the evolutionary mechanisms behind this genomic diversification, we explore here two clades of microsporidia with distinct life histories, Ordospora and Hamiltosporidium, parasitizing the same host species, Daphnia magna. Based on seven newly assembled genomes, we show that mixed-mode transmission (the combination of horizontal and vertical transmission), which occurs in Hamiltosporidium, is found to be associated with larger and AT-biased genomes, more genes, and longer intergenic regions, as compared with the exclusively horizontally transmitted Ordospora. Furthermore, the Hamiltosporidium genome assemblies contain a variety of repetitive elements and long segmental duplications. We show that there is an excess of nonsynonymous substitutions in the microsporidia with mixed-mode transmission, which cannot be solely attributed to the lack of recombination, suggesting that bursts of genome size in these microsporidia result primarily from genetic drift. Overall, these findings suggest that the switch from a horizontal-only to a mixed mode of transmission likely produces population bottlenecks in Hamiltosporidium species, therefore reducing the effectiveness of natural selection, and allowing their genomic features to be largely shaped by nonadaptive processes.

Optimized and affordable high-throughput sequencing workflow for preserved and nonpreserved small zooplankton specimens

Genomic analysis of hundreds of individuals is increasingly becoming standard in evolutionary and ecological research. Individual-based sequencing generates large amounts of valuable data from experimental and field studies, while using preserved samples is an invaluable resource for studying biodiversity in remote areas or across time. Yet, small-bodied individuals or specimens from collections are often of limited use for genomic analyses due to a lack of suitable extraction and library preparation protocols for preserved or small amounts of tissues. Currently, high-throughput sequencing in zooplankton is mostly restricted to clonal species, that can be maintained in live cultures to obtain sufficient amounts of tissue, or relies on a whole-genome amplification step that comes with several biases and high costs. Here, we present a workflow for high-throughput sequencing of single small individuals omitting the need for prior whole-genome amplification or live cultures. We establish and demonstrate this method using 27 species of the genus Daphnia, aquatic keystone organisms, and validate it with small-bodied ostracods. Our workflow is applicable to both live and preserved samples at low costs per sample. We first show that a silica-column based DNA extraction method resulted in the highest DNA yields for nonpreserved samples while a precipitation-based technique gave the highest yield for ethanol-preserved samples and provided the longest DNA fragments. We then successfully performed short-read whole genome sequencing from single Daphnia specimens and ostracods. Moreover, we assembled a draft reference genome from a single Daphnia individual (>50× coverage) highlighting the value of the workflow for non-model organisms.

Lineage Divergence of Dendrolimus punctatus in Southern China Based on Mitochondrial Genome

In southern China, the masson pine caterpillar, Dendrolimus punctatus, has caused serious damage to the Pinus massoniana (Lamb.) pine forests. Here, the whole mitochondrial DNA (mtDNA) was employed to analyze the population evolution of D. punctatus and to understand the process underlying its current phylogenetic pattern. D. punctatus populations within its distribution range in China were categorized into five subgroups: central and eastern China (CEC), southwestern China (SWC), Yibin in Sichuan (SC), Baise in Guangxi (GX), and Luoding in Guangdong (GD), with a high level of haplotype diversity and nucleotide diversity among them. The genetic distances between subgroups are relatively large; however, the genetic distances between populations within the CEC subgroup were relatively small, suggesting that many populations were closely related in this subgroup. The mantel test showed that geographic distance had an important impact on the genetic distance of different geographic populations (r = 0.3633, P < 0.001). The neutrality tests, Bayesian skyline plot, and haplotype network showed that D. punctatus experienced a population expansion around 100,000 years ago. The divergence times of GX/SC, SWC, GD, and CEC were 0.347, 0.236, 0.200, and 0.110 million years ago, respectively. The SWC, CEC, and GD subgroups might have evolved from GX or SC subgroups. The population genetic structure of D. punctatus was closely related to its host tree species, geographic distance among populations, the weak flight capacity, and many eco-environment conditions.

Unexpected endemism in the Daphnia longispina complex (Crustacea: Cladocera) in Southern Siberia

The biological significance of regional cladoceran morphotypes in the montane regions of the central Palearctic remains poorly understood. In the Holarctic Daphnia longispina complex (Cladocera: Daphniidae), several variants, lineages and species have been proposed as endemic for Southern Siberia. Daphnia turbinata Sars, for example, named after its unusual head shape, is known only from Southern Siberia. Here we sequence DNA of Daphnia from three mitochondrial genes (12S rRNA, 16S rRNA, and NADH dehydrogenase subunit 2, ND2) from 57 localities in Russia and Mongolia (the majority being from Southern Siberia) and place them in evolutionary context with existing data. Our aim was to examine regional endemism of the Daphnia longispina complex in Southern Siberian; to improve the phylogenetic understanding with improved taxonomic and regional sampling, and to better understand the influence of Pleistocene glaciation on the biogeography of these lineages. At least three lineages showed genetic evidence for endemism in Southern Siberia. There was strong support for D. turbinata as a sister lineage to to D. longispina/D. dentifera. Another endemic, Siberian D. cf. longispina, is a sister group to the longispina group in general. Within D. longispina s. str. there was an endemic Siberian clade with a western range boundary near the Yenisei River Basin. Gene flow estimates among populations (based on F_ST values) were very low for clades of D. longispina on a regional (the original 12S dataset), and on a pan-Eurasian (the extended 12S dataset) scale. Negative values of Fu’s F_S and Tajima’s D tests prevailed for the species examined with significant values found for two D. longispina clades, D. dentifera, D. galeata and D. cristata. Our results support the notion that Southern Siberia is an important biogeographic region for cladocerans as it contained unexpected diversity of endemics (such as D. turbinata, D. cf. longispina and lineages of D. umbra and D. longsipina s.str.) and from being the geographic meeting place of expanding postglacial lineages from eastern and western refugia.

Temperature- versus precipitation-limitation shape local temperature tolerance in a Holarctic freshwater crustacean

Species with wide geographical distributions are often adapted locally to the prevailing temperatures. To understand how they respond to ongoing climatic change, we must appreciate the interplay between temperature, seasonality and the organism's life cycle. The temperature experienced by many organisms results from an often-overlooked combination of climate and phenology. Summer-active (high latitude) populations are expected to adapt to local summer temperatures, but this is not expected for populations that outlive the summer in their dormant stage (low latitude, precipitation-limited). We recorded reproduction and survival in genotypes from 123 Holarctic populations of Daphnia magna during a multi-generation thermal ramp experiment. Genotypes from summer-active populations showed a positive relationship between heat tolerance and local summer temperature, whereas winter-active populations did not. These findings are consistent with the hypothesis that D. magna adapts to the local temperatures the animals experience during their planktonic phase. We conclude that predicting local temperature adaptation, in particular in the light of climate change, needs to consider the phenology of geographically wide-ranging species.

Mitochondrial Genomics Reveals Shared Phylogeographic Patterns and Demographic History among Three Periodical Cicada Species Groups

The mass application of whole mitogenome (MG) sequencing has great potential for resolving complex phylogeographic patterns that cannot be resolved by partial mitogenomic sequences or nuclear markers. North American periodical cicadas (Magicicada) are well known for their periodical mass emergence at 17- and 13-year intervals in the north and south, respectively. Magicicada comprises three species groups, each containing one 17-year species and one or two 13-year species. Within each life cycle, single-aged cohorts, called broods, of periodical cicadas emerge in different years, and most broods contain members of all three species groups. There are 12 and three extant broods of 17- and 13-year cicadas, respectively. The phylogeographic relationships among the populations and broods within the species groups have not been clearly resolved. We analyzed 125 whole MG sequences from all broods and seven species within three species groups to ascertain the divergence history of the geographic and allochronic populations and their life cycles. Our mitogenomic phylogeny analysis clearly revealed that each of the three species groups had largely similar phylogeographic subdivisions (east, middle, and west) and demographic histories (rapid population expansion after the last glacial period). The mitogenomic phylogeny also partly resolved the brood diversification process, which could be explained by hypothetical temporary life cycle shifts, and showed that none of the 13- and 17-year species within the species groups was monophyletic, possibly due to gene flow between them. Our findings clearly reveal phylogeographic structures in the three Magicicada species groups, demonstrating the advantage of whole MG sequence data in phylogeographic studies.

Accessible molecular phylogenomics at no cost: obtaining 14 new mitogenomes for the ant subfamily Pseudomyrmecinae from public data

The advent of Next Generation Sequencing has reduced sequencing costs and increased genomic projects from a huge amount of organismal taxa, generating an unprecedented amount of genomic datasets publicly available. Often, only a tiny fraction of outstanding relevance of the genomic data produced by researchers is used in their works. This fact allows the data generated to be recycled in further projects worldwide. The assembly of complete mitogenomes is frequently overlooked though it is useful to understand evolutionary relationships among taxa, especially those presenting poor mtDNA sampling at the level of genera and families. This is exactly the case for ants (Hymenoptera:Formicidae) and more specifically for the subfamily Pseudomyrmecinae, a group of arboreal ants with several cases of convergent coevolution without any complete mitochondrial sequence available. In this work, we assembled, annotated and performed comparative genomics analyses of 14 new complete mitochondria from Pseudomyrmecinae species relying solely on public datasets available from the Sequence Read Archive (SRA). We used all complete mitogenomes available for ants to study the gene order conservation and also to generate two phylogenetic trees using both (i) concatenated set of 13 mitochondrial genes and (ii) the whole mitochondrial sequences. Even though the tree topologies diverged subtly from each other (and from previous studies), our results confirm several known relationships and generate new evidences for sister clade classification inside Pseudomyrmecinae clade. We also performed a synteny analysis for Formicidae and identified possible sites in which nucleotidic insertions happened in mitogenomes of pseudomyrmecine ants. Using a data mining/bioinformatics approach, the current work increased the number of complete mitochondrial genomes available for ants from 15 to 29, demonstrating the unique potential of public databases for mitogenomics studies. The wide applications of mitogenomes in research and presence of mitochondrial data in different public dataset types makes the "no budget mitogenomics" approach ideal for comprehensive molecular studies, especially for subsampled taxa.