Insertion Polymorphisms of Mobile Genetic Elements in Sexual and Asexual Populations of Daphnia pulex

Trans-Specific Polymorphisms Between Cryptic Daphnia Species Affect Fitness and Behavior

Shared polymorphisms, loci with identical alleles across species, are of unique interest in evolutionary biology as they may represent cases of selection maintaining ancient genetic variation post-speciation, or contemporary selection promoting convergent evolution. In this study, we investigate the abundance of shared polymorphism between two members of the Daphnia pulex species complex. We test whether the presence of shared mutations is consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting or convergent evolution. We analyzed over 2,000 genomes from six taxa in the D. pulex species group and examined the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that North American and European D. pulex diverged over 10 million years ago, yet retained tens of thousands of shared polymorphisms. We suggest that the number of shared polymorphisms between North American and European D. pulex cannot be fully explained by hybridization or incomplete lineage sorting alone. We show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting convergent and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild.

Global gene expression profiling under nitrogen stress identifies key genes involved in nitrogen stress adaptation in maize (Zea mays L.)

Maize is a heavy consumer of fertilizer nitrogen (N) which not only results in the high cost of cultivation but may also lead to environmental pollution. Therefore, there is a need to develop N-use efficient genotypes, a prerequisite for which is a greater understanding of N-deficiency stress adaptation. In this study, comparative transcriptome analysis was performed using leaf and root tissues from contrasting inbred lines, viz., DMI 56 (tolerant to N stress) and DMI 81 (susceptible to N stress) to delineate the differentially expressed genes (DEGs) under low-N stress. The contrasting lines were grown hydroponically in modified Hoagland solution having either sufficient- or deficient-N, followed by high-throughput RNA-sequencing. A total of 8 sequencing libraries were prepared and 88–97% of the sequenced raw reads were mapped to the reference B73 maize genome. Genes with a p value ≤ 0.05 and fold change of ≥ 2.0 or ≤ − 2 were considered as DEGs in various combinations performed between susceptible and tolerant genotypes. DEGs were further classified into different functional categories and pathways according to their putative functions. Gene Ontology based annotation of these DEGs identified three different functional categories: biological processes, molecular function, and cellular component. The KEGG and Mapman based analysis revealed that most of the DEGs fall into various metabolic pathways, biosynthesis of secondary metabolites, signal transduction, amino acid metabolism, N-assimilation and metabolism, and starch metabolism. Some of the key genes involved in N uptake (high-affinity nitrate transporter 2.2 and 2.5), N assimilation and metabolism (glutamine synthetase, asparagine synthetase), redox homeostasis (SOD, POX), and transcription factors (MYB36, AP2-EREBP) were found to be highly expressed in the tolerant genotype compared to susceptible one. The candidate genes identified in the present study might be playing a pivotal role in low-N stress adaptation in maize and hence could be useful in augmenting further research on N metabolism and development of N-deficiency tolerant maize cultivars.

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.

Genomic Features of Parthenogenetic Animals

Evolution without sex is predicted to impact genomes in numerous ways. Case studies of individual parthenogenetic animals have reported peculiar genomic features that were suggested to be caused by their mode of reproduction, including high heterozygosity, a high abundance of horizontally acquired genes, a low transposable element load, or the presence of palindromes. We systematically characterized these genomic features in published genomes of 26 parthenogenetic animals representing at least 18 independent transitions to asexuality. Surprisingly, not a single feature was systematically replicated across a majority of these transitions, suggesting that previously reported patterns were lineage-specific rather than illustrating the general consequences of parthenogenesis. We found that only parthenogens of hybrid origin were characterized by high heterozygosity levels. Parthenogens that were not of hybrid origin appeared to be largely homozygous, independent of the cellular mechanism underlying parthenogenesis. Overall, despite the importance of recombination rate variation for the evolution of sexual animal genomes, the genome-wide absence of recombination does not appear to have had the dramatic effects which are expected from classical theoretical models. The reasons for this are probably a combination of lineage-specific patterns, the impact of the origin of parthenogenesis, and a survivorship bias of parthenogenetic lineages.

A Survey of Transposon Landscapes in the Putative Ancient Asexual Ostracod Darwinula stevensoni

How asexual reproduction shapes transposable element (TE) content and diversity in eukaryotic genomes remains debated. We performed an initial survey of TE load and diversity in the putative ancient asexual ostracod Darwinula stevensoni. We examined long contiguous stretches of DNA in clones from a genomic fosmid library, totaling about 2.5 Mb, and supplemented these data with results on TE abundance and diversity from an Illumina draft genome. In contrast to other TE studies in putatively ancient asexuals, which revealed relatively low TE content, we found that at least 19% of the fosmid dataset and 26% of the genome assembly corresponded to known transposons. We observed a high diversity of transposon families, including LINE, gypsy, PLE, mariner/Tc, hAT, CMC, Sola2, Ginger, Merlin, Harbinger, MITEs and helitrons, with the prevalence of DNA transposons. The predominantly low levels of sequence diversity indicate that many TEs are or have recently been active. In the fosmid data, no correlation was found between telomeric repeats and non-LTR retrotransposons, which are present near telomeres in other taxa. Most TEs in the fosmid data were located outside of introns and almost none were found in exons. We also report an N-terminal Myb/SANT-like DNA-binding domain in site-specific R4/Dong non-LTR retrotransposons. Although initial results on transposable loads need to be verified with high quality draft genomes, this study provides important first insights into TE dynamics in putative ancient asexual ostracods.

Asexual reproduction reduces transposable element load in experimental yeast populations

Theory predicts that sexual reproduction can either facilitate or restrain transposable element (TE) accumulation by providing TEs with a means of spreading to all individuals in a population, versus facilitating TE load reduction via purifying selection. By quantifying genomic TE loads over time in experimental sexual and asexual Saccharomyces cerevisiae populations, we provide direct evidence that TE loads decrease rapidly under asexual reproduction. We show, using simulations, that this reduction may occur via evolution of TE activity, most likely via increased excision rates. Thus, sex is a major driver of genomic TE loads and at the root of the success of TEs.

The genetic information of most living organisms contains parasitic invaders known as transposable elements. These genetic sequences multiply by copying and pasting themselves through the genome, but this process can disrupt the activity of important genes and put the organism at risk.

How transposable elements proliferate in a population depends on the way organisms reproduce. If they simply clone themselves asexually, the selfish elements cannot spread between the different clones. If the organisms mate together their respective transposable elements get mixed, which helps the sequences to spread more easily and to potentially become more virulent. However, sexual reproduction also comes with mechanisms that keep transposable elements in check.

Bast, Jaron et al. took advantage of the fact that yeasts can reproduce with or without mating to explore whether sexual or asexual organisms are better at controlling the spread of transposable elements. The number of copies of transposable elements in the genomes of yeast grown sexually or asexually was assessed. The results showed that sexual populations kept constant numbers of selfish elements, while asexual organisms lost these genomic parasites over time. Simulations then revealed that this difference emerged because a defense gene that helps to delete transposable elements was spreading more quickly in the asexual group.

The work by Bast, Jaron et al. therefore suggests that sex is responsible for the evolutionary success of transposable elements, while asexual populations can discard these sequences over time. Sex therefore helps genetic parasites, somewhat similar to sexually transmitted diseases, to spread between individuals and remain virulent.

Genetic control of male production in Daphnia pulex

Daphnia normally reproduce by cyclical parthenogenesis, with offspring sex being determined by environmental cues. However, some females have lost the ability to produce males. Our results demonstrate that this loss of male-producing ability is controlled by a dominant allele at a single locus. We identified the locus by comparing whole-genome sequences of 67 nonmale-producing (NMP) and 100 male-producing (MP) clones from 5 Daphnia pulex populations, revealing 132 NMP-linked SNPs and 59 NMP-linked indels within a single 1.1-Mb nonrecombining region on chromosome I. These markers include 7 nonsynonymous mutations, all of which are located within one unannotated protein-coding gene (gene 8960). Within this single gene, all of the marker-linked NMP haplotypes from different populations form a monophyletic clade, suggesting a single origin of the NMP phenotype, with the NMP haplotype originating by introgression from a sister species, Daphnia pulicaria Methyl farnesoate (MF) is the innate juvenile hormone in daphnids, which induces the production of males and whose inhibition results in female-only production. Gene 8960 is sensitive to treatment by MF in MP clones, but such responsiveness is greatly reduced in NMP clones. Thus, we hypothesize that gene 8960 is located downstream of the MF-signaling pathway in D. pulex, with the NMP phenotype being caused by expression change of gene 8960.

A Maximum-Likelihood Approach to Estimating the Insertion Frequencies of Transposable Elements from Population Sequencing Data

Transposable elements (TEs) contribute to a large fraction of the expansion of many eukaryotic genomes due to the capability of TEs duplicating themselves through transposition. A first step to understanding the roles of TEs in a eukaryotic genome is to characterize the population-wide variation of TE insertions in the species. Here, we present a maximum-likelihood (ML) method for estimating allele frequencies and detecting selection on TE insertions in a diploid population, based on the genotypes at TE insertion sites detected in multiple individuals sampled from the population using paired-end (PE) sequencing reads. Tests of the method on simulated data show that it can accurately estimate the allele frequencies of TE insertions even when the PE sequencing is conducted at a relatively low coverage (=5X). The method can also detect TE insertions under strong selection, and the detection ability increases with sample size in a population, although a substantial fraction of actual TE insertions under selection may be undetected. Application of the ML method to genomic sequencing data collected from a natural Daphnia pulex population shows that, on the one hand, most (>90%) TE insertions present in the reference D. pulex genome are either fixed or nearly fixed (with allele frequencies >0.95); on the other hand, among the nonreference TE insertions (i.e., those detected in some individuals in the population but absent from the reference genome), the majority (>70%) are still at low frequencies (<0.1). Finally, we detected a substantial fraction (∼9%) of nonreference TE insertions under selection.

A New Reference Genome Assembly for the Microcrustacean Daphnia pulex

Comparing genomes of closely related genotypes from populations with distinct demographic histories can help reveal the impact of effective population size on genome evolution. For this purpose, we present a high quality genome assembly of Daphnia pulex (PA42), and compare this with the first sequenced genome of this species (TCO), which was derived from an isolate from a population with >90% reduction in nucleotide diversity. PA42 has numerous similarities to TCO at the gene level, with an average amino acid sequence identity of 98.8 and >60% of orthologous proteins identical. Nonetheless, there is a highly elevated number of genes in the TCO genome annotation, with ∼7000 excess genes appearing to be false positives. This view is supported by the high GC content, lack of introns, and short length of these suspicious gene annotations. Consistent with the view that reduced effective population size can facilitate the accumulation of slightly deleterious genomic features, we observe more proliferation of transposable elements (TEs) and a higher frequency of gained introns in the TCO genome.