Genes mirror geography inDaphnia magna

Copy number and sequence variation in rDNA of Daphnia pulex from natural populations: insights from whole-genome sequencing

Ribosomal DNA (rDNA) has a vital role in ribosome biogenesis as it contains the genes that encode ribosomal RNA (rRNA) separated by intergenic spacers (IGSs). The rRNA genes occur in hundreds to tens of thousands of copies per haploid genome in eukaryotes and are generally highly conserved with low variation within species. Due to the repetitive nature and large size of rDNA arrays, detecting intraindividual variation can be difficult. In this study, we use whole-genome sequences of 169 Daphnia pulex individuals from 10 natural populations to measure the copy number and sequence variation in rDNA. This revealed that variation in rDNA copy number between individuals spans an order of magnitude. We further observed a substantial level of sequence variation within individual genomes. As expected, single-nucleotide polymorphisms occurred in regions of lower functional constraint such as the IGS and expansion segments of the rRNA genes. The presence of strong linkage disequilibrium among variants facilitated identification of haplotypes within each population. Although there was evidence of recombination among haplotypes from different populations, it is insufficient to eliminate linkage disequilibrium within populations. Estimating copy number and haplotype diversity within individuals revealed that the level of intraindividual sequence variation is not strongly correlated with copy number. The observed patterns of variation highlight a complex evolutionary history of rDNA in D. pulex. Future research should explore the functional implications of rDNA copy number and sequence variation on organismal phenotypes.

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.

Strong philopatry in an estuarine-dependent fish

Understanding fish movement is critical in determining the spatial scales in which to appropriately manage wild populations. Genetic markers provide a natural tagging approach to assess the degree of gene flow and population connectivity across a species distribution. We investigated the genetic structure of black bream Acanthopagrus butcheri across its entire distribution range in Australia, as well as regional scale gene flow across south-eastern Australia by undertaking a comprehensive analysis of the populations in estuaries across the region. We applied genome-wide sampling of single-nucleotide polymorphism (SNP) markers generated from restriction site-associated DNA sequencing. Genetic structure and potential gene flow was assessed using principal component analyses and admixture analyses (STRUCTURE). Using 33,493 SNPs, we detected broad scale genetic structuring, with limited gene flow among regional clusters (i.e. Western Australia, South Australia and western Victoria; and eastern Victoria, Tasmania and New South Wales). This is likely the result of unsuitable habitats, strong ocean currents (e.g. the Leeuwin Current and the East Australian Current), large water bodies (e.g. Bass Strait) and known biogeographical provinces across the continent. Local-scale genetic structuring was also identified across the south-eastern Australian estuaries sampled, reflecting that the coexistence of both migratory and resident individuals within populations (i.e. partial migration), and the movement of fish into coastal waters, still results in strong philopatry across the region. Instances of movement among estuaries at this spatial scale were primarily found between adjacent estuaries and were likely attributed to lone migrants utilising inshore coastal currents for movement beyond nearby habitats. Targeting SNP markers in A. butcheri at this continental scale highlighted how neither spatial proximity of estuaries nor black bream's ability to move into coastal waters reflects increased gene flow. Overall, our findings highlight the importance of location-specific management.

The limits of stress-tolerance for zooplankton resting stages in freshwater ponds

In seasonal environments, many organisms evolve strategies such as diapause to survive stressful periods. Understanding the link between habitat stability and diapause strategy can help predict a population's survival in a changing world. Indeed, resting stages may be an important way freshwater organisms can survive periods of drought or freezing, and as the frequency and extent of drought or freezing vary strongly among habitats and are predicted to change with climate change, it raises questions about how organisms cope with, and survive, environmental stress. Using Daphnia magna as a model system, we tested the ability of resting stages from different populations to cope with stress during diapause. The combination of elevated temperatures and wet conditions during diapause shows to prevent hatching altogether. In contrast, hatching is relatively higher after a dry and warm diapause, but declines with rising temperatures, while time to hatch increases. Resting stages produced by populations from summer-dry habitats perform slightly, but consistently, better at higher temperatures and dryness, supporting the local adaptation hypothesis. A higher trehalose content in resting eggs from summer-dry habitat might explain such pattern. Considering that temperatures and summer droughts are projected to increase in upcoming years, it is fundamental to know how resting stages resist stressful conditions so as to predict and protect the ecological functioning of freshwater ecosystems.

Evolutionary Genomics of Sister Species Differing in Effective Population Sizes and Recombination Rates

Studies of closely related species with known ecological differences provide exceptional opportunities for understanding the genetic mechanisms of evolution. In this study, we compared population-genomics data between Daphnia pulex and Daphnia pulicaria, two reproductively compatible sister species experiencing ecological speciation, the first largely confined to intermittent ponds and the second to permanent lakes in the same geographic region. Daphnia pulicaria has lower genome-wide nucleotide diversity, a smaller effective population size, a higher incidence of private alleles, and a substantially more linkage disequilibrium than D. pulex. Positively selected genes in D. pulicaria are enriched in potentially aging-related categories such as cellular homeostasis, which may explain the extended life span in D. pulicaria. We also found that opsin-related genes, which may mediate photoperiodic responses, are under different selection pressures in these two species. Genes involved in mitochondrial functions, ribosomes, and responses to environmental stimuli are found to be under positive selection in both species. Additionally, we found that the two species have similar average evolutionary rates at the DNA-sequence level, although approximately 160 genes have significantly different rates in the two lineages. Our results provide insights into the physiological traits that differ within this regionally sympatric sister-species pair that occupies unique microhabitats.

Near chromosome-level genome assembly of the microsporidium Hamiltosporidium tvaerminnensis

Microsporidia are intracellular parasitic fungi whose genomes rank among the smallest of all known eukaryotes. A number of outstanding questions remain concerning the evolution of their large-scale variation in genome architecture, responsible for genome size variation of more than an order of magnitude. This genome report presents the first near-chromosomal assembly of a large-genome microsporidium, Hamiltosporidium tvaerminnensis. Combined Oxford Nanopore, Pacific Biosciences (PacBio), and Illumina sequencing led to a genome assembly of 17 contigs, 11 of which represent complete chromosomes. Our assembly is 21.64 Mb in length, has an N50 of 1.44 Mb, and consists of 39.56% interspersed repeats. We introduce a novel approach in microsporidia, PacBio Iso-Seq, as part of a larger annotation pipeline for obtaining high-quality annotations of 3,573 protein-coding genes. Based on direct evidence from the full-length Iso-Seq transcripts, we present evidence for alternative polyadenylation and variation in splicing efficiency, which are potential regulation mechanisms for gene expression in microsporidia. The generated high-quality genome assembly is a necessary resource for comparative genomics that will help elucidate the evolution of genome architecture in response to intracellular parasitism.

Out of the 'host' box: extreme off-host conditions alter the infectivity and virulence of a parasitic bacterium

Disease agents play an important role in the ecology and life history of wild and cultivated populations and communities. While most studies focus on the adaptation of parasites to their hosts, the adaptation of free-living parasite stages to their external (off-host) environment may tell us a lot about the factors that shape the distribution of parasites. Pasteuria ramosa is an endoparasitic bacterium of the water flea Daphnia with a wide geographical distribution. Its transmission stages rest outside of the host and thus experience varying environmental regimes. We examined the life history of P. ramosa populations from four environmental conditions (i.e. groups of habitats): the factorial combinations of summer-dry water bodies or not, and winter-freeze water bodies or not. Our goal was to examine how the combination of winter temperature and summer dryness affects the parasite's ability to attach to its host and to infect it. We subjected samples of the four groups of habitats to temperatures of 20, 33, 46 and 60°C in dry and wet conditions, and exposed a susceptible clone of Daphnia magna to the treated spores. We found that spores which had undergone desiccation endured higher temperatures better than spores kept wet, both regarding attachment and subsequent infection. Furthermore, spores treated with heightened temperatures were much less infective and virulent. Even under high temperatures (60°C), exposed spores from all populations were able to attach to the host cuticle, albeit they were unable to establish infection. Our work highlights the sensitivity of a host-free resting stage of a bacterial parasite to the external environment. Long heatwaves and harsh summers, which are becoming more frequent owing to recent climate changes, may therefore pose a problem for parasite survival. This article is part of the theme issue 'Infectious disease ecology and evolution in a changing world'.

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.

Microevolutionary variation in molar morphology of Onychomys leucogaster decoupled from genetic structure

In neutral models of quantitative trait evolution, both genetic and phenotypic divergence scale as random walks, producing a correlation between the two measures. However, complexity in the genotype-phenotype map may alter the correlation between genotypic and phenotypic divergence, even when both are evolving neutrally or nearly so. Understanding this correlation between phenotypic and genetic variation is critical for accurately interpreting the fossil record. This study compares the geographic structure and scaling of morphological variation of the shape of the first lower molar of 77 individuals of the northern grasshopper mouse Onychomys leucogaster to genome-wide SNP variation in the same sample. We found strong genetic structure but weak or absent morphological structure indicating that the scaling of each type of variation is decoupled from one another. Low P_ST values relative to F_ST values are consistent with a lack of morphological divergence in contrast to genetic divergence between groups. This lack of phenotypic structure and the presence of notable within-sample phenotypic variance are consistent with uniform selection or constraints on molar shape across a wide geographic and environmental range. Over time, this kind of decoupling may result in patterns of phenotypic stasis masking underlying genetic patterns.

Evolutionary Genomics of a Subdivided Species

The ways in which genetic variation is distributed within and among populations is a key determinant of the evolutionary features of a species. However, most comprehensive studies of these features have been restricted to studies of subdivision in settings known to have been driven by local adaptation, leaving our understanding of the natural dispersion of allelic variation less than ideal. Here, we present a geographic population-genomic analysis of 10 populations of the freshwater microcrustacean Daphnia pulex, an emerging model system in evolutionary genomics. These populations exhibit a pattern of moderate isolation-by-distance, with an average migration rate of 0.6 individuals per generation, and average effective population sizes of ∼650,000 individuals. Most populations contain numerous private alleles, and genomic scans highlight the presence of islands of excessively high population subdivision for more common alleles. A large fraction of such islands of population divergence likely reflect historical neutral changes, including rare stochastic migration and hybridization events. The data do point to local adaptive divergence, although the precise nature of the relevant variation is diffuse and cannot be associated with particular loci, despite the very large sample sizes involved in this study. In contrast, an analysis of between-species divergence highlights positive selection operating on a large set of genes with functions nearly nonoverlapping with those involved in local adaptation, in particular ribosome structure, mitochondrial bioenergetics, light reception and response, detoxification, and gene regulation. These results set the stage for using D. pulex as a model for understanding the relationship between molecular and cellular evolution in the context of natural environments.

Potential differences in chitin synthesis ability cause different sensitivities to diflubenzuron among three strains of Daphnia magna

Ecotoxicity testing of crustaceans using Daphnia magna has been implemented in the chemical management systems of various countries. While the chemical sensitivity of D. magna varies depending on genetically different clonal lineages, the strain used in ecotoxicity tests, including the acute immobilization test (OECD TG202), has not been specified. We hypothesized that comprehensive gene expression profiles could provide useful information on phenotypic differences among strains, including chemical sensitivity. To test this hypothesis, we performed mRNA sequencing on three different strains (NIES, England, and Clone 5) of D. magna under culture conditions. The resulting expression profile of the NIES strain was clearly different compared to the profiles of the other two strains. Gene ontology (GO) enrichment analysis suggested that chitin metabolism was significantly enriched in the NIES strain compared to that in the England strain. Consistent with the GO analysis, evidence of high levels of chitin metabolism in the NIES strain were observed across multiple levels of biological organization, such as expression of chitin synthase genes, chitin content, and chitinase activity, which suggested that the different strains would exhibit different sensitivities to chemicals used to inhibit chitin synthesis. We found that among all strains, the NIES strain was more tolerant to diflubenzuron, a chitin synthesis inhibitor, with a 14-fold difference in the 48 h-EC₅₀ value for the acute immobilization test compared to the England strain. The present study demonstrates that the differences among strains in chitin metabolism may lead to sensitivity difference to diflubenzuron, and serves as a case study of the usefulness of comprehensive gene expression profiles in finding sensitivity differences.

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.

Population connectivity, dispersal, and swimming behavior in Daphnia

The water flea Daphnia has the capacity to respond rapidly to environmental stressors, to disperse over large geographical scales, and to preserve its genetic material by forming egg banks in the sediment. Spatial and temporal distributions of D. magna have been extensively studied over the last decades using behavioral or genetic tools, although the correlation between the two has rarely been the focus. In the present study, we therefore investigated the population genetic structure and behavioral response to a lethal threat, ultraviolet radiation (UVR), among individuals from two different water bodies. Our results show two genetic populations with moderate gene flow, highly correlated with geographical location and with inheritable traits through generations. However, despite the strong genetic differences between populations, we show homogeneous refuge demand between populations when exposed to the lethal threat solar UVR.

Host-parasite co-evolution and its genomic signature

Studies in diverse biological systems have indicated that host-parasite co-evolution is responsible for the extraordinary genetic diversity seen in some genomic regions, such as major histocompatibility (MHC) genes in jawed vertebrates and resistance genes in plants. This diversity is believed to evolve under balancing selection on hosts by parasites. However, the mechanisms that link the genomic signatures in these regions to the underlying co-evolutionary process are only slowly emerging. We still lack a clear picture of the co-evolutionary concepts and of the genetic basis of the co-evolving phenotypic traits in the interacting antagonists. Emerging genomic tools that provide new options for identifying underlying genes will contribute to a fuller understanding of the co-evolutionary process.

Genetic diversity of a widespread annual killifish from coastal Tanzania

Background

African annual killifishes (Nothobranchius spp.) are adapted to seasonally desiccating habitats (ephemeral pools), surviving dry periods as dormant eggs. Given their peculiar life history, geographic aspects of their diversity uniquely combine patterns typical for freshwater taxa (river basin structure and elevation gradient) and terrestrial animals (rivers acting as major dispersal barriers). However, our current knowledge on fine-scale inter-specific and intra-specific genetic diversity of African annual fish is limited to a single, particularly dry region of their distribution (subtropical Mozambique). Using a widespread annual killifish from coastal Tanzania and Kenya, we tested whether the same pattern of genetic divergence pertains to a wet equatorial region in the centre of Nothobranchius distribution.

Results

In populations of Nothobranchius melanospilus species group across its range, we genotyped a part of mitochondrial cytochrome oxidase subunit 1 (COI) gene (83 individuals from 22 populations) and 10 nuclear microsatellite markers (251 individuals from 16 populations). We found five lineages with a clear phylogeographic structure but frequent secondary contact. Mitochondrial lineages were largely congruent with main population genetic clusters identified on microsatellite markers. In the upper Wami basin, populations are isolated as a putative Nothobranchius prognathus, but include also a population from a periphery of the middle Ruvu basin. Other four lineages (including putative Nothobranchius kwalensis) coexisted in secondary contact zones, but possessed clear spatial pattern. Main river channels did not form apparent barriers to dispersal. The most widespread lineage had strong signal of recent population expansion.

Conclusions

We conclude that dispersal of a Nothobranchius species from a wet part of the genus distribution (tropical lowland) is not constrained by main river channels and closely related lineages frequently coexist in secondary contact zones. We also demonstrate contemporary connection between the Ruvu and Rufiji river basins. Our data do not provide genetic support for existence of recently described cryptic species from N. melanospilus complex, but cannot resolve this issue.

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.

Fitness and Genomic Consequences of Chronic Exposure to Low Levels of Copper and Nickel in Daphnia pulex Mutation Accumulation Lines

In at least some unicellular organisms, mutation rates are temporarily raised upon exposure to environmental stress, potentially contributing to the evolutionary response to stress. Whether this is true for multicellular organisms, however, has received little attention. This study investigated the effects of chronic mild stress, in the form of low-level copper and nickel exposure, on mutational processes in Daphnia pulex using a combination of mutation accumulation, whole genome sequencing and life-history assays. After over 100 generations of mutation accumulation, we found no effects of metal exposure on the rates of single nucleotide mutations and of loss of heterozygosity events, the two mutation classes that occurred in sufficient numbers to allow statistical analysis. Similarly, rates of decline in fitness, as measured by intrinsic rate of population increase and of body size at first reproduction, were negligibly affected by metal exposure. We can reject the possibility that Daphnia were insufficiently stressed to invoke genetic responses as we have previously shown rates of large-scale deletions and duplications are elevated under metal exposure in this experiment. Overall, the mutation accumulation lines did not significantly depart from initial values for phenotypic traits measured, indicating the lineage used was broadly mutationally robust. Taken together, these results indicate that the mutagenic effects of chronic low-level exposure to these metals are restricted to certain mutation classes and that fitness consequences are likely minor and therefore unlikely to be relevant in determining the evolutionary responses of populations exposed to these stressors.

The Combined Effect of Temperature and Host Clonal Line on the Microbiota of a Planktonic Crustacean

Host-associated microbiota vary across host individuals and environmental conditions, but the relative importance of their genetic background versus their environment is difficult to disentangle. We sought to experimentally determine the factors shaping the microbiota of the planktonic Crustacean, Daphnia magna. We used clonal lines from a wide geographic distribution, which had been kept under standardized conditions for over 75 generations. Replicate populations were kept for three generations at 20 and 28 °C. The interaction of the host clonal line and environment (i.e., temperature) influenced microbiota community characteristics, including structure, the relative abundance of common microbial species, and the microbial richness and phylogenetic diversity. We did not find any correlation between host-associated microbiota and the geographic origin of the clones or their temperature tolerance. Our results highlight the prominent effects that host clonal lineage and its interaction with the environment has on host-associated microbiota composition.

Phylogeography of Daphnia magna Straus (Crustacea: Cladocera) in Northern Eurasia: Evidence for a deep longitudinal split between mitochondrial lineages

Species with a large geographic distributions present a challenge for phylogeographic studies due to logistic difficulties of obtaining adequate sampling. For instance, in most species with a Holarctic distribution, the majority of studies has concentrated on the European or North American part of the distribution, with the Eastern Palearctic region being notably understudied. Here, we study the phylogeography of the freshwater cladoceran Daphnia magna Straus, 1820 (Crustacea: Cladocera), based on partial mitochondrial COI sequences and using specimens from populations spread longitudinally from westernmost Europe to easternmost Asia, with many samples from previously strongly understudied regions in Siberia and Eastern Asia. The results confirm the previously suspected deep split between Eastern and Western mitochondrial haplotype super-clades. We find a narrow contact zone between these two super-clades in the eastern part of Western Siberia, with proven co-occurrence in a single lake in the Novosibirsk region. However, at present there is no evidence suggesting that the two mitochondrial super-clades represent cryptic species. Rather, they may be explained by secondary contact after expansion from different refugia. Interestingly, Central Siberia has previously been found to be an important contact zone also in other cladoceran species, and may thus be a crucial area for understanding the Eurasian phylogeography of freshwater invertebrates. Together, our study provides an unprecedented complete, while still not global, picture of the phylogeography of this important model species.

QTL mapping of a natural genetic polymorphism for long-term parasite persistence in Daphnia populations

Knowing the determinants of the geographic ranges of parasites is important for understanding their evolutionary ecology, epidemiology and their potential to expand their range. Here we explore the determinants of geographic range in the peculiar case of a parasite species – the microsporidian Hamiltosporidium tvaerminnensis – that has a limited geographic distribution in a wide-spread host – Daphnia magna. We conducted a quantitative trait loci (QTLs) analysis with monoclonal F2D. magna populations originating from a cross between a susceptible northern European genotype and a resistant central European genotype. Contrary to our expectations, long-term persistence turned out to be a quantitative trait across the F2 offspring. Evidence for two QTLs, one epistatic interaction and for further minor QTL was found. This finding contrasts markedly with the previously described bimodal pattern for long-term parasite persistence in natural host genotypes across Europe and leaves open the question of how a quantitative genetic trait could determine the disjunct geographic distribution of the parasite across Europe.

Selection Constrains High Rates of Tandem Repetitive DNA Mutation in Daphnia pulex

A long-standing evolutionary puzzle is that all eukaryotic genomes contain large amounts of tandemly-repeated DNA whose sequence motifs and abundance vary greatly among even closely related species. To elucidate the evolutionary forces governing tandem repeat dynamics, quantification of the rates and patterns of mutations in repeat copy number and tests of its selective neutrality are necessary. Here, we used whole-genome sequences of 28 mutation accumulation (MA) lines of Daphnia pulex, in addition to six isolates from a non-MA population originating from the same progenitor, to both estimate mutation rates of abundances of repeat sequences and evaluate the selective regime acting upon them. We found that mutation rates of individual repeats were both high and highly variable, ranging from additions/deletions of 0.29-105 copies per generation (reflecting changes of 0.12-0.80% per generation). Our results also provide evidence that new repeat sequences are often formed from existing ones. The non-MA population isolates showed a signal of either purifying or stabilizing selection, with 33% lower variation in repeat copy number on average than the MA lines, although the level of selective constraint was not evenly distributed across all repeats. The changes between many pairs of repeats were correlated, and the pattern of correlations was significantly different between the MA lines and the non-MA population. Our study demonstrates that tandem repeats can experience extremely rapid evolution in copy number, which can lead to high levels of divergence in genome-wide repeat composition between closely related species.

Digital gene expression analysis with sample multiplexing and PCR duplicate detection: A straightforward protocol

Tag-Seq is a high-throughput approach used for discovering SNPs and characterizing gene expression. In comparison to RNA-Seq, Tag-Seq eases data processing and allows detection of rare mRNA species using only one tag per transcript molecule. However, reduced library complexity raises the issue of PCR duplicates, which distort gene expression levels. Here we present a novel Tag-Seq protocol that uses the least biased methods for RNA library preparation combined with a novel approach for joint PCR template and sample labeling. In our protocol, input RNA is fragmented by hydrolysis, and poly(A)-bearing RNAs are selected and directly ligated to mixed DNA-RNA P5 adapters. The P5 adapters contain i5 barcodes composed of sample-specific (moderately) degenerate base regions (mDBRs), which later allow detection of PCR duplicates. The P7 adapter is attached via reverse transcription with individual i7 barcodes added during the amplification step. The resulting libraries can be sequenced on an Illumina sequencer. After sample demultiplexing and PCR duplicate removal with a free software tool we designed, the data are ready for downstream analysis. Our protocol was tested on RNA samples from predator-induced and control Daphnia microcrustaceans.

Phylogeographic and population genetic analyses reveal multiple species of Boa and independent origins of insular dwarfism

Boa is a Neotropical genus of snakes historically recognized as monotypic despite its expansive distribution. The distinct morphological traits and color patterns exhibited by these snakes, together with the wide diversity of ecosystems they inhabit, collectively suggest that the genus may represent multiple species. Morphological variation within Boa also includes instances of dwarfism observed in multiple offshore island populations. Despite this substantial diversity, the systematics of the genus Boa has received little attention until very recently. In this study we examined the genetic structure and phylogenetic relationships of Boa populations using mitochondrial sequences and genome-wide SNP data obtained from RADseq. We analyzed these data at multiple geographic scales using a combination of phylogenetic inference (including coalescent-based species delimitation) and population genetic analyses. We identified extensive population structure across the range of the genus Boa and multiple lines of evidence for three widely-distributed clades roughly corresponding with the three primary land masses of the Western Hemisphere. We also find both mitochondrial and nuclear support for independent origins and parallel evolution of dwarfism on offshore island clusters in Belize and Cayos Cochinos Menor, Honduras.