Knowledge of the genome-wide rate and spectrum of mutations is necessary to understand the origin of disease and the genetic variation driving all evolutionary processes. Here, we provide a genome-wide analysis of the rate and spectrum of mutations obtained in two Daphnia pulex genotypes via separate mutation-accumulation (MA) experiments. Unlike most MA studies that utilize haploid, homozygous, or self-fertilizing lines, D. pulex can be propagated ameiotically while maintaining a naturally heterozygous, diploid genome, allowing the capture of the full spectrum of genomic changes that arise in a heterozygous state. While base-substitution mutation rates are similar to those in other multicellular eukaryotes (about 4 × 10 −9 per site per generation), we find that the rates of large-scale (>100 kb) de novo copy-number variants (CNVs) are significantly elevated relative to those seen in previous MA studies. The heterozygosity maintained in this experiment allowed for estimates of gene-conversion processes. While most of the conversion tract lengths we report are similar to those generated by meiotic processes, we also find larger tract lengths that are indicative of mitotic processes. Comparison of MA lines to natural isolates reveals that a majority of large-scale CNVs in natural populations are removed by purifying selection. The mutations observed here share similarities with disease-causing, complex, large-scale CNVs, thereby demonstrating that MA studies in D. pulex serve as a system for studying the processes leading to such alterations.

Hybridization plays a potentially important role in the origin of obligate parthenogenesis (OP) in many organisms. However, it remains controversial whether hybridization directly triggers the transition from sexual reproduction to obligate asexuality or a hybrid genetic background enables asexual species to persist. Furthermore, we know little about the specific genetic elements from the divergent, yet still hybridizing lineages responsible for this transition and how these elements are further spread to create other OP lineages. In this study, we address these questions in Daphnia pulex, where cyclically parthenogenetic (CP) and OP lineages coexist. Ancestry estimates and whole-genome association mapping using 32 OP isolates suggest that a complex hybridization history between the parental species D. pulex and D. pulicaria is responsible for the introgression of a set of 647 D. pulicaria single nucleotide polymorphism alleles that show perfect association with OP. Crossing experiments using males of OP lineages and females of CP lineages strongly support a polygenic basis for OP. Single-sperm analyses show that although normal meiotic recombination occurs in the production of haploid sperm by males of OP lineages, a significant proportion of such sperm are polyploid, suggesting that the spread of asexual elements through these males (i.e., contagious asexuality) is much less efficient than previously envisioned. Although the current Daphnia genome annotation does not provide mechanistic insight into the nature of the asexuality-associated alleles, these alleles should be considered as candidates for future investigations on the genetic underpinnings of OP.

Background A mechanistic understanding of the spread of SARS-CoV-2 and diligent tracking of ongoing mutagenesis are of key importance to plan robust strategies for confining its transmission. Large numbers of available sequences and their dates of transmission provide an unprecedented opportunity to analyze evolutionary adaptation in novel ways. Addition of high-resolution structural information can reveal the functional basis of these processes at the molecular level. Integrated systems biology-directed analyses of these data layers afford valuable insights to build a global understanding of the COVID-19 pandemic. Results Here we identify globally distributed haplotypes from 15,789 SARS-CoV-2 genomes and model their success based on their duration, dispersal, and frequency in the host population. Our models identify mutations that are likely compensatory adaptive changes that allowed for rapid expansion of the virus. Functional predictions from structural analyses indicate that, contrary to previous reports, the Asp614Gly mutation in the spike glycoprotein (S) likely reduced transmission and the subsequent Pro323Leu mutation in the RNA-dependent RNA polymerase led to the precipitous spread of the virus. Our model also suggests that two mutations in the nsp13 helicase allowed for the adaptation of the virus to the Pacific Northwest of the USA. Finally, our explainable artificial intelligence algorithm identified a mutational hotspot in the sequence of S that also displays a signature of positive selection and may have implications for tissue or cell-specific expression of the virus. Conclusions These results provide valuable insights for the development of drugs and surveillance strategies to combat the current and future pandemics.

Many organisms survive fluctuating and extreme environmental conditions by manifesting multiple distinct phenotypes during adulthood by means of developmental processes that enable phenotypic plasticity. We report on the discovery of putative plasticity-enabling genes that are involved in transforming the gill of the euryhaline teleost fish, Fundulus heteroclitus, from its freshwater to its seawater gill-type, a process that alters both morphology and function. Gene expression that normally enables osmotic plasticity is inhibited by arsenic. Gene sets defined by antagonistic interactions between arsenic and salinity show reduced transcriptional variation among individual fish, suggesting unusually accurate and precise regulatory control of these genes, consistent with the hypothesis that they participate in a canalized developmental response. We observe that natural selection acts to preserve canalized gene expression in populations of killifish that are most tolerant to abrupt salinity change and that these populations show the least variability in their transcription of genes enabling plasticity of the gill. We found that genes participating in this highly canalized and conserved plasticity-enabling response had significantly fewer and less complex associations with transcriptional regulators than genes that respond only to arsenic or salinity. Collectively these findings, which are drawn from the relationships between environmental challenge, plasticity, and canalization among populations, suggest that the selective processes that facilitate phenotypic plasticity do so by targeting the regulatory networks that gives rise to the response. These findings also provide a generalized, conceptual framework of how genes might interact with the environment and evolve toward the development of plastic traits.

Copy number variation (CNV) of genes coding for certain enzymes has been shown to be responsible for adaptation of arthropods to anthropogenic toxins. Natural toxins produced by cyanobacteria in freshwater ecosystems, that is, protease inhibitors (PIs), have been demonstrated to increase in frequency over the last decades due to eutrophication and global warming. These PIs inhibit digestive proteases of Daphnia, the major herbivore of phytoplankton and cyanobacteria. The adjustment of isoforms, differences in gene expression, and activity of gut proteases determine tolerance to dietary PIs in single Daphnia genotypes. Here, we tested whether similar mechanisms are also responsible for differences in tolerance among Daphnia population. We developed a droplet digital PCR (ddPCR) method for the analysis of CNV of Daphnia proteases. We report that one Daphnia protease gene showed CNV between populations and that CNV correlates with chymotrypsin gene expression among populations. We showed that populations of Daphnia magna differ in tolerance to cyanobacterial PIs according to the cyanobacterial background of their lake of origin, which hints at local adaptation. The tolerance of the populations correlates with IC values of their chymotrypsins, which is probably due to a combined effect of CNV (translating into gene expression differences) and positive selection of tolerant protease isoforms. This is the first study using ddPCR to demonstrate CNV of a gene with ecologically relevant function, and the first report of differences in tolerance to cyanobacterial PIs among Daphnia populations in combination with the assessment of underlying molecular mechanisms.