Search

Genomic regions that play a role in parasite defense are often found to be highly variable, with the major histocompatibility complex serving as an iconic example. Single nucleotide polymorphisms may represent only a small portion of this variability, with Indel polymorphisms and copy number variation further contributing. In extreme cases, haplotypes may no longer be recognized as orthologous. Understanding the evolution of such highly divergent regions is challenging because the most extreme variation is not visible using reference-assisted genomic approaches. Here we analyze the case of the Pasteuria Resistance Complex in the crustacean Daphnia magna , a defense complex in the host against the common and virulent bacterium Pasteuria ramosa. Two haplotypes of this region have been previously described, with parts of it being nonhomologous, and the region has been shown to be under balancing selection. Using pan-genome analysis and tree reconciliation methods to explore the evolution of the Pasteuria Resistance Complex and its characteristics within and between species of Daphnia and other Cladoceran species, our analysis revealed a remarkable diversity in this region even among host species, with many nonhomologous hyper-divergent haplotypes. The Pasteuria Resistance Complex is characterized by extensive duplication and losses of Fucosyltransferase (FuT) and Galactosyltransferase (GalT) genes that are believed to play a role in parasite defense. The Pasteuria Resistance Complex region can be traced back to common ancestors over 250 million years. The unique combination of an ancient resistance complex and a dynamic, hyper-divergent genomic environment presents a fascinating opportunity to investigate the role of such regions in the evolution and long-term maintenance of resistance polymorphisms. Our findings offer valuable insights into the evolutionary forces shaping disease resistance and adaptation, not only in the genus Daphnia , but potentially across the entire Cladocera class. [ABSTRACT FROM AUTHOR]

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

We have sequenced, assembled, and analyzed the nuclear and mitochondrial genomes and transcriptomes of Potamopyrgus estuarinus and Potamopyrgus kaitunuparaoa , two prosobranch snail species native to New Zealand that together span the continuum from estuary to freshwater. These two species are the closest known relatives of the freshwater species Potamopyrgus antipodarum— a model for studying the evolution of sex, host–parasite coevolution, and biological invasiveness—and thus provide key evolutionary context for understanding its unusual biology. The P. estuarinus and P. kaitunuparaoa genomes are very similar in size and overall gene content. Comparative analyses of genome content indicate that these two species harbor a near-identical set of genes involved in meiosis and sperm functions, including seven genes with meiosis-specific functions. These results are consistent with obligate sexual reproduction in these two species and provide a framework for future analyses of P. antipodarum— a species comprising both obligately sexual and obligately asexual lineages, each separately derived from a sexual ancestor. Genome-wide multigene phylogenetic analyses indicate that P. kaitunuparaoa is likely the closest relative to P. antipodarum. We nevertheless show that there has been considerable introgression between P. estuarinus and P. kaitunuparaoa. That introgression does not extend to the mitochondrial genome, which appears to serve as a barrier to hybridization between P. estuarinus and P. kaitunuparaoa. Nuclear-encoded genes whose products function in joint mitochondrial-nuclear enzyme complexes exhibit similar patterns of nonintrogression, indicating that incompatibilities between the mitochondrial and the nuclear genome may have prevented more extensive gene flow between these two species. [ABSTRACT FROM AUTHOR]

The angiosperm genus Silene has been the subject of extensive study in the field of ecology and evolution, but the availability of high-quality reference genome sequences has been limited for this group. Here, we report a chromosome-level assembly for the genome of Silene conica based on Pacific Bioscience HiFi, Hi-C, and Bionano technologies. The assembly produced 10 scaffolds (1 per chromosome) with a total length of 862 Mb and only ∼1% gap content. These results confirm previous observations that S. conica and its relatives have a reduced base chromosome number relative to the genus's ancestral state of 12. Silene conica has an exceptionally large mitochondrial genome (>11 Mb), predominantly consisting of sequence of unknown origins. Analysis of shared sequence content suggests that it is unlikely that transfer of nuclear DNA is the primary driver of this mitochondrial genome expansion. More generally, this assembly should provide a valuable resource for future genomic studies in Silene , including comparative analyses with related species that recently evolved sex chromosomes. [ABSTRACT FROM AUTHOR]

Understanding the genomic basis of infectious disease is a fundamental objective in co-evolutionary theory with relevance to healthcare, agriculture, and epidemiology. Models of host-parasite co-evolution often assume that infection requires specific combinations of host and parasite genotypes. Co-evolving host and parasite loci are, therefore, expected to show associations that reflect an underlying infection/resistance allele matrix, yet little evidence for such genome-to-genome interactions has been observed among natural populations. We conducted a study to search for this genomic signature across 258 linked host (Daphnia magna) and parasite (Pasteuria ramosa) genomes. Our results show a clear signal of genomic association between multiple epistatically interacting loci in the host genome, and a family of genes encoding for collagen-like protein in the parasite genome. These findings are supported by laboratory-based infection trials, which show strong correspondence between phenotype and genotype at the identified loci. Our study provides clear genomic evidence of antagonistic co-evolution among wild populations. [ABSTRACT FROM AUTHOR]

Recent advances in long-read sequencing technology have allowed for single-molecule sequencing of entire mitochondrial genomes, opening the door for direct investigation of the mitochondrial genome architecture and recombination. We used PacBio sequencing to reassemble mitochondrial genomes from two species of New Zealand freshwater snails, Potamopyrgus antipodarum and Potamopyrgus estuarinus. These assemblies revealed a ∼1.7 kb structure within the mitochondrial genomes of both species that was previously undetected by an assembly of short reads and likely corresponding to a large noncoding region commonly present in the mitochondrial genomes. The overall architecture of these Potamopyrgus mitochondrial genomes is reminiscent of the chloroplast genomes of land plants, harboring a large single-copy (LSC) region and a small single-copy (SSC) region separated by a pair of inverted repeats (IRa and IRb). Individual sequencing reads that spanned across the Potamopyrgus IRa-SSC-IRb structure revealed the occurrence of a "flip-flop" recombination. We also detected evidence for two distinct IR haplotypes and recombination between them in wild-caught P. estuarinus , as well as extensive intermolecular recombination between single-nucleotide polymorphisms in the LSC region. The chloroplast-like architecture and repeat-mediated mitochondrial recombination we describe here raise fundamental questions regarding the origins and commonness of inverted repeats in cytoplasmic genomes and their role in mitochondrial genome evolution. [ABSTRACT FROM AUTHOR]

The dynamics of extinction and (re)colonization in habitat patches are characterizing features of dynamic metapopulations, causing them to evolve differently than large, stable populations. The propagule model, which assumes genetic bottlenecks during colonization, posits that newly founded subpopulations have low genetic diversity and are genetically highly differentiated from each other. Immigration may then increase diversity and decrease differentiation between subpopulations. Thus, older and/or less isolated subpopulations are expected to have higher genetic diversity and less genetic differentiation. We tested this theory using whole-genome pool-sequencing to characterize nucleotide diversity and differentiation in 60 subpopulations of a natural metapopulation of the cyclical parthenogen Daphnia magna. For comparison, we characterized diversity in a single, large, and stable D. magna population. We found reduced (synonymous) genomic diversity, a proxy for effective population size, weak purifying selection, and low rates of adaptive evolution in the metapopulation compared with the large, stable population. These differences suggest that genetic bottlenecks during colonization reduce effective population sizes, which leads to strong genetic drift and reduced selection efficacy in the metapopulation. Consistent with the propagule model, we found lower diversity and increased differentiation in younger and also in more isolated subpopulations. Our study sheds light on the genomic consequences of extinction–(re)colonization dynamics to an unprecedented degree, giving strong support for the propagule model. We demonstrate that the metapopulation evolves differently from a large, stable population and that evolution is largely driven by genetic drift. [ABSTRACT FROM AUTHOR]

Although phylogeny estimation is notoriously difficult in radiations that occurred several hundred million years ago, phylogenomic approaches offer new ways to examine relationships among ancient lineages and evaluate hypotheses that are key to evolutionary biology. Here, we reconstruct the deep-rooted relationships of one of the oldest living arthropod clades, the branchiopod crustaceans, using a kaleidoscopic approach. We use concatenation and coalescent tree-building methods to analyze a large multigene data set at the nucleotide and amino acid level and examine gene tree versus species tree discordance. We unequivocally resolve long-debated relationships among extant orders of the Cladocera, the waterfleas, an ecologically relevant zooplankton group in global aquatic and marine ecosystems that is famous for its model systems in ecology and evolution. To build the data set, we assembled eight de novo genomes of key taxa including representatives of all extant cladoceran orders and suborders. Our phylogenetic analysis focused on a BUSCO-based set of 823 conserved single-copy orthologs shared among 23 representative taxa spanning all living branchiopod orders, including 11 cladoceran families. Our analysis supports the monophyly of the Cladocera and reveals remarkable homoplasy in their body plans. We found large phylogenetic distances between lineages with similar ecological specializations, indicating independent evolution in major body plans, such as in the pelagic predatory orders Haplopoda and Onychopoda (the "Gymnomera"). In addition, we assessed rapid cladogenesis by estimating relative timings of divergence in major lineages using reliable fossil-calibrated priors on eight nodes in the branchiopod tree, suggesting a Paleozoic origin around 325 Ma for the cladoceran ancestor and an ancient rapid radiation around 252 Ma at the Perm/Triassic boundary. These findings raise new questions about the roles of homoplasy and rapid radiation in the diversification of the cladocerans and help examine trait evolution from a genomic perspective in a functionally well understood, ancient arthropod group. [Cladocera; Daphnia ; evolution; homoplasy; molecular clock; phylogenomics; systematics; waterfleas.] [ABSTRACT FROM AUTHOR]

Intracellular transfers of mitochondrial DNA continue to shape nuclear genomes. Chromosome 2 of the model plant Arabidopsis thaliana contains one of the largest known nuclear insertions of mitochondrial DNA (numts). Estimated at over 600 kb in size, this numt is larger than the entire Arabidopsis mitochondrial genome. The primary Arabidopsis nuclear reference genome contains less than half of the numt because of its structural complexity and repetitiveness. Recent data sets generated with improved long-read sequencing technologies (PacBio HiFi) provide an opportunity to finally determine the accurate sequence and structure of this numt. We performed a de novo assembly using sequencing data from recent initiatives to span the Arabidopsis centromeres, producing a gap-free sequence of the Chromosome 2 numt, which is 641 kb in length and has 99.933% nucleotide sequence identity with the actual mitochondrial genome. The numt assembly is consistent with the repetitive structure previously predicted from fiber-based fluorescent in situ hybridization. Nanopore sequencing data indicate that the numt has high levels of cytosine methylation, helping to explain its biased spectrum of nucleotide sequence divergence and supporting previous inferences that it is transcriptionally inactive. The original numt insertion appears to have involved multiple mitochondrial DNA copies with alternative structures that subsequently underwent an additional duplication event within the nuclear genome. This work provides insights into numt evolution, addresses one of the last unresolved regions of the Arabidopsis reference genome, and represents a resource for distinguishing between highly similar numt and mitochondrial sequences in studies of transcription, epigenetic modifications, and de novo mutations. [ABSTRACT FROM AUTHOR]

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

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. [ABSTRACT FROM AUTHOR]

Understanding how diversity is maintained in natural populations is a major goal of evolutionary biology. In coevolving hosts and parasites, negative frequency‐dependent selection is one mechanism predicted to maintain genetic variation. While much is known about host diversity, parasite diversity remains understudied in coevolutionary research. Here, we survey natural diversity in a bacterial parasite by characterizing infection phenotypes for over 50 isolates in relation to 12 genotypes of their host, Daphnia magna. We find striking phenotypic variation among parasite isolates, and we discover the parasite can infect its host through at least five different attachment sites. Variation in attachment success at each site is explained to varying degrees by host and parasite genotypes. A spatial correlation analysis showed that infectivity of different isolates does not correlate with geographic distance, meaning isolates from widespread populations are equally able to infect the host. Overall, our results reveal that infection phenotypes of this parasite are highly diverse. Our results are consistent with the prediction that under Red Queen coevolutionary dynamics both the host and the parasite should show high genetic diversity for traits of functional importance in their interactions. [ABSTRACT FROM AUTHOR]

How does asexual reproduction influence genome evolution? Although is it clear that genomic structural variation is common and important in natural populations, we know very little about how one of the most fundamental of eukaryotic traits—mode of genomic inheritance—influences genome structure. We address this question with the New Zealand freshwater snail Potamopyrgus antipodarum , which features multiple separately derived obligately asexual lineages that coexist and compete with otherwise similar sexual lineages. We used whole-genome sequencing reads from a diverse set of sexual and asexual individuals to analyze genomic abundance of a critically important gene family, rDNA (the genes encoding rRNAs), that is notable for dynamic and variable copy number. Our genomic survey of rDNA in P. antipodarum revealed two striking results. First, the core histone and 5S rRNA genes occur between tandem copies of the 18S–5.8S–28S gene cluster, a unique architecture for these crucial gene families. Second, asexual P. antipodarum harbor dramatically more rDNA–histone copies than sexuals, which we validated through molecular and cytogenetic analysis. The repeated expansion of this genomic region in asexual P. antipodarum lineages following distinct transitions to asexuality represents a dramatic genome structural change associated with asexual reproduction—with potential functional consequences related to the loss of sexual reproduction. [ABSTRACT FROM AUTHOR]

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

Knowledge of the genetic architecture of pathogen infectivity and host resistance is essential for a mechanistic understanding of coevolutionary processes, yet the genetic basis of these interacting traits remains unknown for most host–pathogen systems. We used a comparative genomic approach to explore the genetic basis of infectivity in Pasteuria ramosa , a Gram-positive bacterial pathogen of planktonic crustaceans that has been established as a model for studies of Red Queen host–pathogen coevolution. We sequenced the genomes of a geographically, phenotypically, and genetically diverse collection of P. ramosa strains and performed a genome-wide association study to identify genetic correlates of infection phenotype. We found multiple polymorphisms within a single gene, Pcl7 , that correlate perfectly with one common and widespread infection phenotype. We then confirmed this perfect association via Sanger sequencing in a large and diverse sample set of P. ramosa clones. Pcl7 codes for a collagen-like protein, a class of adhesion proteins known or suspected to be involved in the infection mechanisms of a number of important bacterial pathogens. Consistent with expectations under Red Queen coevolution, sequence variation of Pcl7 shows evidence of balancing selection, including extraordinarily high diversity and absence of geographic structure. Based on structural homology with a collagen-like protein of Bacillus anthracis , we propose a hypothesis for the structure of Pcl7 and the physical location of the phenotype-associated polymorphisms. Our results offer strong evidence for a gene governing infectivity and provide a molecular basis for further study of Red Queen dynamics in this model host–pathogen system. [ABSTRACT FROM AUTHOR]

The rate and spectrum of spontaneous mutations are critical parameters in basic and applied biology because they dictate the pace and character of genetic variation introduced into populations, which is a prerequisite for evolution. We use a mutation–accumulation approach to estimate mutation parameters from whole-genome sequence data from multiple genotypes from multiple populations of Daphnia magna , an ecological and evolutionary model system. We report extremely high base substitution mutation rates (⁠ µ - n , bs = 8.96 × 10−9/bp/generation [95% CI: 6.66–11.97 × 10−9/bp/generation] in the nuclear genome and µ - m , bs = 8.7 × 10−7/bp/generation [95% CI: 4.40–15.12 × 10−7/bp/generation] in the mtDNA), the highest of any eukaryote examined using this approach. Levels of intraspecific variation based on the range of estimates from the nine genotypes collected from three populations (Finland, Germany, and Israel) span 1 and 3 orders of magnitude, respectively, resulting in up to a ∼300-fold difference in rates among genomic partitions within the same lineage. In contrast, mutation spectra exhibit very consistent patterns across genotypes and populations, suggesting the mechanisms underlying the mutational process may be similar, even when the rates at which they occur differ. We discuss the implications of high levels of intraspecific variation in rates, the importance of estimating gene conversion rates using a mutation–accumulation approach, and the interacting factors influencing the evolution of mutation parameters. Our findings deepen our knowledge about mutation and provide both challenges to and support for current theories aimed at explaining the evolution of the mutation rate, as a trait, across taxa. [ABSTRACT FROM AUTHOR]

The planktonic freshwater crustacean of the genus Daphnia are a model system for biomedical research and, in particular, invertebrate-parasite interactions. Up until now, no virus has been characterized for this system. Here we report the discovery of an iridovirus as the causative agent of White Fat Cell Disease (WFCD) in Daphnia. WFCD is a highly virulent disease of Daphnia that can easily be cultured under laboratory conditions. Although it has been studied from sites across Eurasia for more than 60 years, its causative agent had not been described, nor had an iridovirus been connected to WFCD before now. Here we find that an iridovirus--the Daphnia iridescent virus 1 (DIV-1)--is the causative agent of WFCD. DIV-1 has a genome sequence of about 288 kbp, with 39% G+C content and encodes 367 predicted open reading frames. DIV-1 clusters together with other invertebrate iridoviruses but has by far the largest genome among all sequenced iridoviruses. Comparative genomics reveal that DIV-1 has apparently recently lost a substantial number of unique genes but has also gained genes by horizontal gene transfer from its crustacean host. DIV-1 represents the first invertebrate iridovirus that encodes proteins to purportedly cap RNA, and it contains unique genes for a DnaJ-like protein, a membrane glycoprotein and protein of the immunoglobulin superfamily, which may mediate host-pathogen interactions and pathogenicity. Our findings end a 60-year search for the causative agent of WFCD and add to our knowledge of iridovirus genomics and invertebrate-virus interactions. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 fulllength 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. [ABSTRACT FROM AUTHOR]

Outbred laboratory mice (Mus musculus) are readily available and have high fecundity, making them a popular choice in biomedical research, especially toxicological and pharmacological applications. Direct high throughput genome sequencing (HTS) of these widely used research animals is an important genetic quality control measure that enhances research reproducibility. HTS data have been used to confirm the common origin of outbred stocks and to molecularly define distinct outbred populations. But these data have also revealed unexpected population structure and homozygosity in some populations; genetic features that emerge when outbred stocks are not properly maintained. We used exome sequencing to discover and interrogate protein-coding variation in a newly established population of Swiss-derived outbred stock (J:ARC) that is closely related to other, commonly used CD-1 outbred populations. We used these data to describe the genetic architecture of the J:ARC population including heterozygosity, minor allele frequency, LD decay, and we defined novel, protein-coding sequence variation. These data reveal the expected genetic architecture for a properly maintained outbred stock and provide a basis for the on-going genetic quality control. We also compared these data to proteincoding variation found in a multiparent outbred stock, the Diversity Outbred (J:DO). We found that the more recently derived, multiparent outbred stock has significantly higher interindividual variability, greater overall genetic variation, higher heterozygosity, and fewer novel variants than the Swiss-derived J:ARC stock. However, among the novel variants found in the J:DO stock, significantly more are predicted to be protein-damaging. The fact that individuals from this population can tolerate a higher load of potentially damaging variants highlights the buffering effects of allelic diversity and the differing selective pressures in these stocks. While both outbred stocks offer significant individual heterozygosity, our data provide a molecular basis for their intended applications, where the J:DO are best suited for studies requiring maximum, population-level genetic diversity and power for mapping, while the J:ARC are best suited as a general-purpose outbred stock with robust fecundity, relatively low allelic diversity, and less potential for extreme phenotypic variability. [ABSTRACT FROM AUTHOR]