The Association Between Breeding System and Transposable Element Dynamics in Daphnia Pulex

Comparative genomics of tadpole shrimps (Crustacea, Branchiopoda, Notostraca): Dynamic genome evolution against the backdrop of morphological stasis

This analysis presents five genome assemblies of four Notostraca taxa. Notostraca origin dates to the Permian/Upper Devonian and the extant forms show a striking morphological similarity to fossil taxa. The comparison of sequenced genomes with other Branchiopoda genomes shows that, despite the morphological stasis, Notostraca share a dynamic genome evolution with high turnover for gene families' expansion/contraction and a transposable elements content comparable to other branchiopods. While Notostraca substitutions rate appears similar or lower in comparison to other branchiopods, a subset of genes shows a faster evolutionary pace, highlighting the difficulty of generalizing about genomic stasis versus dynamism. Moreover, we found that the variation of Triops cancriformis transposable elements content appeared linked to reproductive strategies, in line with theoretical expectations. Overall, besides providing new genomic resources for the study of these organisms, which appear relevant for their ecology and evolution, we also confirmed the decoupling of morphological and molecular evolution.

Engines of change: Transposable element mutation rates are high and variable within Daphnia magna

Transposable elements (TEs) represent a major portion of most eukaryotic genomes, yet little is known about their mutation rates or how their activity is shaped by other evolutionary forces. Here, we compare short- and long-term patterns of genome-wide mutation accumulation (MA) of TEs among 9 genotypes from three populations of Daphnia magna from across a latitudinal gradient. While the overall proportion of the genome comprised of TEs is highly similar among genotypes from Finland, Germany, and Israel, populations are distinguishable based on patterns of insertion site polymorphism. Our direct rate estimates indicate TE movement is highly variable (net rates ranging from -11.98 to 12.79 x 10-5 per copy per generation among genotypes), differing both among populations and TE families. Although gains outnumber losses when selection is minimized, both types of events appear to be highly deleterious based on their low frequency in control lines where propagation is not limited to random, single-progeny descent. With rate estimates 4 orders of magnitude higher than base substitutions, TEs clearly represent a highly mutagenic force in the genome. Quantifying patterns of intra- and interspecific variation in TE mobility with and without selection provides insight into a powerful mechanism generating genetic variation in the genome.

SINE Retrotransposon variation drives Ecotypic disparity in natural populations of Coilia nasus

Background

SINEs are a type of nonautonomous retrotransposon that can transpose from one site to be integrated elsewhere in an organism genome. SINE insertion can give rise to genetic variants and regulate gene expression, allowing organisms to acquire new adaptive capacity. Studies on this subject have focused on the impacts of SINEs on genes. However, ecological disparities in fish have not yet been explained by SINEs.

Results

New SINEs were isolated from Coilia nasus, which has two ecotypes—migratory and resident—that differ in their spawning and migration behaviors. The SINEs possess two structures that resemble a tRNA gene and a LINE retrotransposon tail. Comparison of olfactory tissue transcriptomes, intact SINE transcript copies were detected in only the migratory fish at the initial retrotransposition stage. The SINE DNA copy numbers were higher in the resident type than in the migratory type, while the frequency of SINE insertion was higher in the migratory type than in the resident type. Furthermore, SINE insertions can lead to new repeats of short DNA fragments in the genome, along with target site duplications. SINEs in the resident type have undergone excision via a mechanism in which predicted cleavage sites are formed by mutations, resulting in gaps that are then filled by microsatellites via microhomology-induced replication.

Conclusions

Notably, SINEs in the resident type have undergone strong natural selection, causing genomic heteroplasmy and driving ecological diversity of C. nasus. Our results reveal possible evolutionary mechanisms underlying the ecological diversity at the interface between SINE mobilization and organism defense.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

Proteomic Profile of Daphnia pulex using Data-Independent Acquisition Mass Spectrometry and Ion Mobility Separation

Daphnia pulex is a keystone species for aquatic habitats and an ecological/evolution model organism. Although significant progress has been made on characterizing its genome, the D. pulex proteome remains largely uncharacterized partially due to abnormally high protein degradation during homogenization and emphasis on genomic analysis. In this study, various sample preparation and mass spectrometry acquisition methods are performed for the purpose of improving D. pulex proteome exploration. Benefits for employing both in-gel and in-solution methods of trypsin digestion are observed. Furthermore, acquisition methods employing ion mobility separation greatly increase peptide identification and more than doubled the proteome coverage. Bioinformatic analysis suggests that mitochondrial and hydrolytic activities are enriched in D. pulex compared to closely related invertebrates or Homo sapiens. Also, novel D. pulex proteins possessing putative genome modifying functional domains are identified. Data are available via ProteomeXchange with identifier PXD008455.

Obligatory parthenogenesis and TE load: Bacillus stick insects and the R2 non-LTR retrotransposon

Transposable elements (TEs) are selfish genetic elements whose self-replication is contrasted by the host genome. In this context, host reproductive strategies are predicted to impact on both TEs load and activity. The presence and insertion distribution of the non-LTR retrotransposon R2 was here studied in populations of the strictly bisexual Bacillus grandii maretimi and of the obligatory parthenogenetic Bacillus atticus atticus. Furthermore, data were also obtained from the offspring of selected B. a. atticus females. At the population level, the gonochoric B. g. maretimi showed a significantly higher R2 load than the obligatory parthenogenetic B. a. atticus. The comparison with bisexual and unisexual Bacillus rossius populations showed that their values were higher than those recorded for B. a. atticus and similar, or even higher, than those of B. g. maretimi. Consistently, an R2 load reduction is scored in B. a. atticus offspring even if with a great variance. On the whole, data here produced indicate that in the obligatory unisexual B. a. atticus R2 is active and that mechanisms of molecular turnover are effective. Furthermore, progeny analyses show that, at variance of the facultative parthenogenetic B. rossius, the R2 activity is held at a lower rate. Modeling parental-offspring inheritance, suggests that in B. a. atticus recombination plays a major role in eliminating insertions rather than selection, as previously suggested for unisexual B. rossius progeny, even if in both cases a high variance is observed. In addition to this, mechanisms of R2 silencing or chances of clonal selection cannot be ruled out.

Distribution of the DNA transposon family, Pokey in the Daphnia pulex species complex

BACKGROUND

The Pokey family of DNA transposons consists of two putatively autonomous groups, PokeyA and PokeyB, and two groups of Miniature Inverted-repeat Transposable Elements (MITEs), mPok1 and mPok2. This TE family is unusual as it inserts into a specific site in ribosomal (r)DNA, as well as other locations in Daphnia genomes. The goals of this study were to determine the distribution of the Pokey family in lineages of the Daphnia pulex species complex, and to test the hypothesis that unusally high PokeyA number in some isolates of Daphnia pulicaria is the result of recent transposition. To do this, we estimated the haploid number of Pokey, mPok, and rRNA genes in 45 isolates from five Daphnia lineages using quantitative PCR. We also cloned and sequenced partial copies of PokeyA from four isolates of D. pulicaria.

RESULTS

Haploid PokeyA and PokeyB number is generally less than 20 and tends to be higher outside rDNA in four lineages. Conversely, the number of both groups is much higher outside rDNA (~120) in D. arenata, and PokeyB is also somewhat higher inside rDNA. mPok1 was only detected in D. arenata. mPok2 occurs both outside (~30) and inside rDNA (~6) in D. arenata, but was rare (≤2) outside rDNA in the other four lineages. There is no correlation between Pokey and rRNA gene number (mean = 240 across lineages) in any lineage. Variation among cloned partial PokeyA sequences is significantly higher in isolates with high number compared to isolates with an average number.

CONCLUSIONS

The high Pokey number outside rDNA in D. arenata and inside rDNA in some D. pulicaria isolates is consistent with a recent increase in transposition rate. The D. pulicaria increase may have been triggered by insertion of PokeyA into a region of transcriptionally active rDNA. The expansion in D. arenata (thought to be of hybrid origin) may be a consequence of release from epigenetic repression following hybridization. Previous work found D. obtusa to be very different from the D. pulex complex; mean PokeyA is higher in rDNA (~75), rDNA array size is nearly twice as large (415), and the two are positively correlated. The predominance of Pokey in only one location could be explained by purifying selection against ectopic recombination between elements inside and outside rDNA.

Diversity and evolution of Ty1-copia retroelements within Chalcidoidea by reverse transcriptase domain analysis

Ty1-copia retrotransposons are widespread and diverse in insects. Some features of their hosts, such as mating and genetic systems, are predicted to influence the spread of selfish genetic elements like Ty1-copia. Using part of the reverse transcriptase gene as a reference, we experimentally surveyed Ty1-copia elements in eight species of fig wasps (Hymenoptera: Chalcidoidea), and performed an in silico analysis of six available genomes of chalcid wasps. Contrary to initial expectations that selfish elements such as Ty1-copia would be purged from the genomes of these species because of inbreeding and haplodiploidy, almost all of these wasps harbour an abundance of diverse Ty1-copia elements. Phylogenetic analyses suggest that the families of Ty1-copia elements found in these species have had a long association with their chalcid hosts. These results suggest an evolutionary scenario in which there was ancestral polymorphism followed by some taxa-specific events including stochastic loss and further diversification. Furthermore, estimating natural selection within the internal and terminal portions of the Ty1-copia phylogenies demonstrated that the elements are under strong evolutionary constraints for their long-term survival, but evolve like pseudogenes in the short term, accompanied by the rise and fall of parasitic elements in the history of wasp lineage.

Copy number of the transposon, Pokey, in rDNA is positively correlated with rDNA copy number in Daphnia obtuse [corrected]

Pokey is a class II DNA transposon that inserts into 28S ribosomal RNA (rRNA) genes and other genomic regions of species in the subgenus, Daphnia. Two divergent lineages, PokeyA and PokeyB have been identified. Recombination between misaligned rRNA genes changes their number and the number of Pokey elements. We used quantitative PCR (qPCR) to estimate rRNA gene and Pokey number in isolates from natural populations of Daphnia obtusa, and in clonally-propagated mutation accumulation lines (MAL) initiated from a single D. obtusa female. The change in direction and magnitude of Pokey and rRNA gene number did not show a consistent pattern across ∼ 87 generations in the MAL; however, Pokey and rRNA gene number changed in concert. PokeyA and 28S gene number were positively correlated in the isolates from both natural populations and the MAL. PokeyB number was much lower than PokeyA in both MAL and natural population isolates, and showed no correlation with 28S gene number. Preliminary analysis did not detect PokeyB outside rDNA in any isolates and detected only 0 to 4 copies of PokeyA outside rDNA indicating that Pokey may be primarily an rDNA element in D. obtusa. The recombination rate in this species is high and the average size of the rDNA locus is about twice as large as that in other Daphnia species such as D. pulicaria and D. pulex, which may have facilitated expansion of PokeyA to much higher numbers in D. obtusa rDNA than these other species.

Evolutionary dynamics of R2 retroelement and insertion inheritance in the genome of bisexual and parthenogenetic Bacillus rossius populations (Insecta Phasmida)

Theoretical and empirical studies have shown differential management of transposable elements in organisms with different reproductive strategies. To investigate this issue, we analysed the R2 retroelement structure and variability in parthenogenetic and bisexual populations of Bacillus rossius stick insects, as well as insertions inheritance in the offspring of parthenogenetic isolates and of crosses. The B. rossius genome hosts a functional (R2Br(fun) ) and a degenerate (R2Br(deg) ) element, their presence correlating with neither reproductive strategies nor population distribution. The median-joining network method indicated that R2Br(fun) duplicates through a multiple source model, while R2Br(deg) is apparently still duplicating via a master gene model. Offspring analyses showed that unisexual and bisexual offspring have a similar number of R2Br-occupied sites. Multiple or recent shifts from gonochoric to parthenogenetic reproduction may explain the observed data. Moreover, insertion frequency spectra show that higher-frequency insertions in unisexual offspring significantly outnumber those in bisexual offspring. This suggests that unisexual offspring eliminate insertions with lower efficiency. A comparison with simulated insertion frequencies shows that inherited insertions in unisexual and bisexual offspring are significantly different from the expectation. On the whole, different mechanisms of R2 elimination in unisexual vs bisexual offspring and a complex interplay between recombination effectiveness, natural selection and time can explain the observed data.

Impact of ploidy level on the distribution of Pokey element insertions in the Daphnia pulex complex

Background

Transposable elements (TEs) play a major role in genome evolution. Their capacity to move and/or multiply in the genome of their host may have profound impacts on phenotypes and dramatic consequences on genome structure. The population dynamics and distribution of TEs are influenced by their mode of transposition, the availability of niches in host genomes, and host population dynamics. Theories predict an increase in the number of TE insertions following hybridization or polyploidization. Evolution of TEs in hybrids and polyploids has mostly been studied in plants; few studies have examined the impacts of hybridization and/or polyploidization on TEs in animals. Hybrids and polyploids have arisen multiple times in the Daphnia pulex complex and are thought to reproduce by obligate parthenogenesis. Our study examines the effects of ploidy level on polymorphism and number of Pokey element insertions in diploid and polyploid hybrid isolates from the Daphnia pulex complex.

Results

The polymorphism of Pokey insertion sites did not depend solely on either the ploidy level or the genetic background of their host; therefore, it may be the result of interactions between these parameters and other parameters such as Pokey activity, selection and/or drift. No significant effect of ploidy level was found on the number of Pokey insertions using TE display and qPCR. However, the load of Pokey insertion sites and the number of unique insertion sites were slightly (but not significantly) higher in polyploids than in diploids.

Conclusions

These results suggest a lack of increase in the number of Pokey insertions following polyploidization but higher availability of Pokey insertion sites in polyploids than in diploids. Compared to previous TE display and qPCR results, the load of Pokey insertions in hybrid diploids was higher than in non-hybrid sexual and asexual diploids, which suggests an increase in the density of Pokey insertions following hybridization.

In and out of the rRNA genes: characterization of Pokey elements in the sequenced Daphnia genome

Background

Only a few transposable elements are known to exhibit site-specific insertion patterns, including the well-studied R-element retrotransposons that insert into specific sites within the multigene rDNA. The only known rDNA-specific DNA transposon, Pokey (superfamily: piggyBac) is found in the freshwater microcrustacean, Daphnia pulex. Here, we present a genome-wide analysis of Pokey based on the recently completed whole genome sequencing project for D. pulex.

Results

Phylogenetic analysis of Pokey elements recovered from the genome sequence revealed the presence of four lineages corresponding to two divergent autonomous families and two related lineages of non-autonomous miniature inverted repeat transposable elements (MITEs). The MITEs are also found at the same 28S rRNA gene insertion site as the Pokey elements, and appear to have arisen as deletion derivatives of autonomous elements. Several copies of the full-length Pokey elements may be capable of producing an active transposase. Surprisingly, both families of Pokey possess a series of 200 bp repeats upstream of the transposase that is derived from the rDNA intergenic spacer (IGS). The IGS sequences within the Pokey elements appear to be evolving in concert with the rDNA units. Finally, analysis of the insertion sites of Pokey elements outside of rDNA showed a target preference for sites similar to the specific sequence that is targeted within rDNA.

Conclusions

Based on the target site preference of Pokey elements and the concerted evolution of a segment of the element with the rDNA unit, we propose an evolutionary path by which the ancestors of Pokey elements have invaded the rDNA niche. We discuss how specificity for the rDNA unit may have evolved and how this specificity has played a role in the long-term survival of these elements in the subgenus Daphnia.

Genetic and epigenetic changes involving (retro)transposons in animal hybrids and polyploids

Transposable elements (TEs) are discrete genetic units that have the ability to change their location within chromosomal DNA, and constitute a major and rapidly evolving component of eukaryotic genomes. They can be subdivided into 2 distinct types: retrotransposons, which use an RNA intermediate for transposition, and DNA transposons, which move only as DNA. Rapid advances in genome sequencing significantly improved our understanding of TE roles in genome shaping and restructuring, and studies of transcriptomes and epigenomes shed light on the previously unknown molecular mechanisms underlying genetic and epigenetic TE controls. Knowledge of these control systems may be important for better understanding of reticulate evolution and speciation in the context of bringing different genomes together by hybridization and perturbing the established regulatory balance by ploidy changes. See also sister article focusing on plants by Bento et al. in this themed issue.

Evolution of a transposon in Daphnia hybrid genomes

Background

Transposable elements play a major role in genome evolution. Their capacity to move and/or multiply in the genome of their host may have profound impacts on phenotypes, and may have dramatic consequences on genome structure. Hybrid and polyploid clones have arisen multiple times in the Daphnia pulex complex and are thought to reproduce by obligate parthenogenesis. Our study examines the evolution of a DNA transposable element named Pokey in the D. pulex complex.

Results

Portions of Pokey elements inserted in the 28S rRNA genes from various Daphnia hybrids (diploids and polyploids) were sequenced and compared to sequences from a previous study to understand the evolutionary history of the elements. Pokey sequences show a complex phylogenetic pattern. We found evidence of recombination events in numerous Pokey alleles from diploid and polyploid hybrids and also from non-hybrid diploids. The recombination rate in Pokey elements is comparable to recombination rates previously estimated for 28S rRNA genes in the congener, Daphnia obtusa. Some recombinant Pokey alleles were encountered in Daphnia isolates from multiple locations and habitats.

Conclusions

Phylogenetic and recombination analyses showed that recombination is a major force that shapes Pokey evolution. Based on Pokey phylogenies, reticulation has played and still plays an important role in shaping the diversity of the D. pulex complex. Horizontal transfer of Pokey seems to be rare and hybrids often possess Pokey elements derived from recombination among alleles encountered in the putative parental species. The insertion of Pokey in hotspots of recombination may have important impacts on the diversity and fitness of this transposable element.

Genome evolution in outcrossing versus selfing versus asexual species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes, and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. However, while asexual species essentially conform to theoretical predictions, the situation is less simple in selfing species. We discuss the possible reasons to potentially explain this paradox. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

Transposon proliferation in an asexual parasitoid

The widespread occurrence of sex is one of the most elusive problems in evolutionary biology. Theory predicts that asexual lineages can be driven to extinction by uncontrolled proliferation of vertically transmitted transposable elements (TEs), which accumulate because of the inefficiency of purifying selection in the absence of sex and recombination. To test this prediction, we compared genome-wide TE load between a sexual lineage of the parasitoid wasp Leptopilina clavipes and a lineage of the same species that is rendered asexual by Wolbachia-induced parthenogenesis. We obtained draft genome sequences at 15-20× coverage of both the sexual and the asexual lineages using next-generation sequencing. We identified transposons of most major classes in both lineages. Quantification of TE abundance using coverage depth showed that copy numbers in the asexual lineage exceeded those in the sexual lineage for DNA transposons, but not LTR and LINE-like elements. However, one or a small number of gypsy-like LTR elements exhibited a fourfold higher coverage in the asexual lineage. Quantitative PCR showed that high loads of this gypsy-like TE were characteristic for 11 genetically distinct asexual wasp lineages when compared to sexual lineages. We found no evidence for an overall increase in copy number for all TE types in asexuals as predicted by theory. Instead, we suggest that the expansions of specific TEs are best explained as side effects of (epi)genetic manipulations of the host genome by Wolbachia. Asexuality is achieved in a myriad of ways in nature, many of which could similarly result in TE proliferation.

Copy number variation of ribosomal DNA and Pokey transposons in natural populations of Daphnia

BACKGROUND

Despite their ubiquity and high diversity in eukaryotic genomes, DNA transposons are rarely encountered in ribosomal DNA (rDNA). In contrast, R-elements, a diverse group of non-LTR retrotransposons, specifically target rDNA. Pokey is a DNA transposon that targets a specific rDNA site, but also occurs in many other genomic locations, unlike R-elements. However, unlike most DNA transposons, Pokey has been a stable component of Daphnia genomes for over 100 million years. Here we use qPCR to estimate the number of 18S and 28S ribosomal RNA genes and Pokey elements in rDNA (rPokey), as well as other genomic locations (gPokey) in two species of Daphnia. Our goals are to estimate the correlation between (1) the number of 18S and 28S rRNA genes, (2) the number of 28S genes and rPokey, and (3) the number of rPokey and gPokey. In addition, we ask whether Pokey number and distribution in both genomic compartments are affected by differences in life history between D. pulex and D. pulicaria.

RESULTS

We found differences in 18S and 28S gene number within isolates that are too large to be explained by experimental variation. In general, Pokey number within isolates is modest (< 20), and most are gPokey. There is no correlation between the number of rRNA genes and rPokey, or between rPokey and gPokey. However, we identified three isolates with unusually high numbers of both rPokey and gPokey, which we infer is a consequence of recent transposition. We also detected other rDNA insertions (rInserts) that could be degraded Pokey elements, R- elements or the divergent PokeyB lineage recently detected in the Daphnia genome sequence. Unlike rPokey, rInserts are positively correlated with rRNA genes, suggesting that they are amplified by the same mechanisms that amplify rDNA units even though rPokey is not. Overall, Pokey frequency and distribution are similar in D. pulex and D. pulicaria suggesting that differences in life history have no impact on Pokey.

CONCLUSIONS

The possibility that many rDNA units do not contain a copy of both 18S and 28S genes suggests that rDNA is much more complicated than once thought, and warrants further study. In addition, the lack of correlation between rPokey, gPokey and rDNA unit numbers suggests that Pokey transposition rate is generally very low, and that recombination, in combination with natural selection, eliminates rPokey much faster than gPokey. Our results suggest that further research to determine the mechanisms by which Pokey has escaped complete inactivation by its host (the usual fate of DNA transposons), would provide important insights into transposon biology.

LTR retroelements in the genome of Daphnia pulex

BACKGROUND

Long terminal repeat (LTR) retroelements represent a successful group of transposable elements (TEs) that have played an important role in shaping the structure of many eukaryotic genomes. Here, we present a genome-wide analysis of LTR retroelements in Daphnia pulex, a cyclical parthenogen and the first crustacean for which the whole genomic sequence is available. In addition, we analyze transcriptional data and perform transposon display assays of lab-reared lineages and natural isolates to identify potential influences on TE mobility and differences in LTR retroelements loads among individuals reproducing with and without sex.

RESULTS

We conducted a comprehensive de novo search for LTR retroelements and identified 333 intact LTR retroelements representing 142 families in the D. pulex genome. While nearly half of the identified LTR retroelements belong to the gypsy group, we also found copia (95), BEL/Pao (66) and DIRS (19) retroelements. Phylogenetic analysis of reverse transcriptase sequences showed that LTR retroelements in the D. pulex genome form many lineages distinct from known families, suggesting that the majority are novel. Our investigation of transcriptional activity of LTR retroelements using tiling array data obtained from three different experimental conditions found that 71 LTR retroelements are actively transcribed. Transposon display assays of mutation-accumulation lines showed evidence for putative somatic insertions for two DIRS retroelement families. Losses of presumably heterozygous insertions were observed in lineages in which selfing occurred, but never in asexuals, highlighting the potential impact of reproductive mode on TE abundance and distribution over time. The same two families were also assayed across natural isolates (both cyclical parthenogens and obligate asexuals) and there were more retroelements in populations capable of reproducing sexually for one of the two families assayed.

CONCLUSIONS

Given the importance of LTR retroelements activity in the evolution of other genomes, this comprehensive survey provides insight into the potential impact of LTR retroelements on the genome of D. pulex, a cyclically parthenogenetic microcrustacean that has served as an ecological model for over a century.

DNA transposons and the role of recombination in mutation accumulation in Daphnia pulex

Background

We identify DNA transposons from the completed draft genome sequence of Daphnia pulex, a cyclically parthenogenetic, aquatic microcrustacean of the class Branchiopoda. In addition, we experimentally quantify the abundance of six DNA transposon families in mutation-accumulation lines in which sex is either promoted or prohibited in order to better understand the role of recombination in transposon proliferation.

Results

We identified 55 families belonging to 10 of the known superfamilies of DNA transposons in the genome of D. pulex. DNA transposons constitute approximately 0.7% of the genome. We characterized each family and, in many cases, identified elements capable of activity in the genome. Based on assays of six putatively active element families in mutation-accumulation lines, we compared DNA transposon abundance in lines where sex was either promoted or prohibited. We find the major difference in abundance in sexuals relative to asexuals in lab-reared lines is explained by independent assortment of heterozygotes in lineages where sex has occurred.

Conclusions

Our examination of the duality of sex as a mechanism for both the spread and elimination of DNA transposons in the genome reveals that independent assortment of chromosomes leads to significant copy loss in lineages undergoing sex. Although this advantage may offset the so-called 'two fold cost of sex' in the short-term, if insertions become homozygous at specific loci due to recombination, the advantage of sex may be decreased over long time periods. Given these results, we discuss the potential effects of sex on the dynamics of DNA transposons in natural populations of D. pulex.

DNA transposon dynamics in populations of Daphnia pulex with and without sex

We investigate the role of recombination in transposable element (TE) proliferation in the cyclical parthenogen, Daphnia pulex. Recombination provides a mechanism by which the rate of both TE gain and loss can be accelerated, a duality that has long intrigued many biologists interested in the influence of sex on mutation accumulation. We compared TE loads among populations of D. pulex where sex occurs regularly (cyclical parthenogens or 'sexuals') with those in which the ability to reproduce sexually has been completely lost (obligate 'asexuals') for six different families of DNA transposons. Transposon display assays showed that sexuals have more TEs than asexuals, contrary to the expectations under Muller's ratchet but consistent with the idea that sex facilitates TE spread. Sexuals also exhibit higher insertion site polymorphism among lineages, as predicted because recombination accelerates rates of loss and gain. Asexuals, however, have proportionally more singletons (loci occupied in a single isolate), which differs from previous studies where selfing and outcrossing were used as a proxy for high and low recombination. Our multi-element survey reveals that the impact of sex on TE proliferation is consistent among different Class II TE families and we discuss the genomic consequences of different reproductive strategies over long time periods.

Transposable elements and factors influencing their success in eukaryotes

Recent advances in genome sequencing have led to a vast accumulation of transposable element data. Consideration of the genome sequencing projects in a phylogenetic context reveals that despite the hundreds of eukaryotic genomes that have been sequenced, a strong bias in sampling exists. There is a general under-representation of unicellular eukaryotes and a dearth of genome projects in many branches of the eukaryotic phylogeny. Among sequenced genomes, great variation in genome size exists, however, little difference in the total number of cellular genes is observed. For many eukaryotes, the remaining genomic space is extremely dynamic and predominantly composed of a menagerie of populations of transposable elements. Given the dynamic nature of the genomic niche filled by transposable elements, it is evident that these elements have played an important role in genome evolution. The contribution of transposable elements to genome architecture and to the advent of genetic novelty is likely to be dependent, at least in part, on the transposition mechanism, diversity, number, and rate of turnover of transposable elements in the genome at any given time. The focus of this review is the discussion of some of the forces that act to shape transposable element diversity within and between genomes.

The effect of transposon Pokey insertions on sequence variation in the 28S rRNA gene of Daphnia pulex

The goal of this study was to determine the impact of breeding system and the presence of the transposon Pokey on intraindividual variation in 28S rRNA genes. We PCR-amplified, cloned, and sequenced 1000 nucleotides downstream of the Pokey insertion site in genes with and without insertions from 10 obligately and 10 cyclically parthenogenetic isolates of Daphnia pulex. Variation among genes with Pokey insertions was higher than variation among genes without insertions in both cyclic and obligate isolates. Although the differences were not quite significant (p = 0.06 in both cases), the results suggest that Pokey insertions are likely to inhibit the homogenization of their host genes to some extent. We also observed that the complement of 28S rRNA alleles differed between genes with and without inserts in some isolates, suggesting that a particular inserted gene can persist for substantial periods of time and even spread within the rDNA array, despite the fact that insertions are deleterious. This apparently contradictory pattern can be explained if homogenization of rRNA genes occurs primarily by gene conversion, but copies with Pokey inserts can occasionally increase in frequency within arrays owing to unequal crossing over events that do not originate in the inserted genes themselves.