Male-biased mutation rates occur in a diverse array of organisms. The ratio of male-to-female mutation rate may have major ramifications for evolution across the genome, and for sex-linked genes in particular. In ZW species, the Z chromosome is carried by males two-thirds of the time, leading to the prediction that male-biased mutation rates will have a disproportionate effect on the evolution of Z-linked genes relative to autosomes and the W chromosome. Colubroid snakes (including colubrids, elapids, and viperids) have ZW sex determination, yet male-biased mutation rates have not been well studied in this group. Here we analyze a population genomic dataset from rattlesnakes to quantify genetic variation within and genetic divergence between species. We use a new method for unbiased estimation of population genetic summary statistics to compare variation between the Z chromosome and autosomes and to calculate net nucleotide differentiation between species. We find evidence for a 2.03-fold greater mutation rate in male rattlesnakes relative to females, corresponding to an average μZ/μA ratio of 1.1. Our results from snakes are quantitatively similar to birds, suggesting that male-biased mutation rates may be a common feature across vertebrate lineages with ZW sex determination.

Rising global temperatures threaten to disrupt population sex ratios, which can in turn cause mate shortages, reduce population growth and adaptive potential, and increase extinction risk, particularly when ratios are male biased. Sex ratio distortion can then have cascading effects across other species and even ecosystems. Our understanding of the problem is limited by how often studies measure temperature effects in both sexes. To address this, the current review surveyed 194 published studies of heat tolerance, finding that the majority did not even mention the sex of the individuals used, with <10% reporting results for males and females separately. Although the data are incomplete, this review assessed phylogenetic patterns of thermally induced sex ratio bias for 3 different mechanisms: sex-biased heat tolerance, temperature-dependent sex determination (TSD), and temperature-induced sex reversal. For sex-biased heat tolerance, documented examples span a large taxonomic range including arthropods, chordates, protists, and plants. Here, superior heat tolerance is more common in females than males, but the direction of tolerance appears to be phylogenetically fluid, perhaps due to the large number of contributing factors. For TSD, well-documented examples are limited to reptiles, where high temperature usually favors females, and fishes, where high temperature consistently favors males. For temperature-induced sex reversal, unambiguous cases are again limited to vertebrates, and high temperature usually favors males in fishes and amphibians, with mixed effects in reptiles. There is urgent need for further work on the full taxonomic extent of temperature-induced sex ratio distortion, including joint effects of the multiple contributing mechanisms.

Bison are an icon of the American West and an ecologically, commercially, and culturally important species. Despite numbering in the hundreds of thousands today, conservation concerns remain for the species, including the impact on genetic diversity of a severe bottleneck around the turn of the 20th century and genetic introgression from domestic cattle. Genetic diversity and admixture are best evaluated at genome-wide scale, for which a high-quality reference is necessary. Here, we use trio binning of long reads from a bison-Simmental cattle (Bos taurus taurus) male F1 hybrid to sequence and assemble the genome of the American plains bison (Bison bison bison). The male haplotype genome is chromosome-scale, with a total length of 2.65 Gb across 775 scaffolds (839 contigs) and a scaffold N50 of 87.8 Mb. Our bison genome is ~13× more contiguous overall and ~3400× more contiguous at the contig level than the current bison reference genome. The bison genome sequence presented here (ARS-UCSC_bison1.0) will enable new research into the evolutionary history of this iconic megafauna species and provide a new tool for the management of bison populations in federal and commercial herds.

Erratum in J Hered. 2021 Nov 1;112(6):565.

Erratum in J Hered. 2021 May 10;:

Sun bear populations are fragmented and at risk from habitat loss and exploitation for body parts. These threats are made worse by significant gaps in knowledge of sun bear population genetic diversity, population connectivity, and taxonomically significant management units. Using a complete sun bear mitochondrial genome, we developed a set of mitochondrial markers to assess haplotype variation and the evolutionary history of sun bears from Peninsular (West) Malaysia and Sabah (East Malaysia). Genetic samples from 28 sun bears from Peninsular Malaysia, 36 from Sabah, and 18 from Thailand were amplified with primers targeting a 1800 bp region of the mitochondrial genome including the complete mitochondrial control region and adjacent genes. Sequences were analyzed using phylogenetic methods. We identified 51 mitochondrial haplotypes among 82 sun bears. Phylogenetic and network analyses provided strong support for a deep split between Malaysian sun bears and sun bears in East Thailand and Yunnan province in China. The Malaysian lineage was further subdivided into two clades: Peninsular Malaysian and Malaysian Borneo (Sabah). Sun bears from Thailand occurred in both Sabah and Peninsular Malaysian clades. Our study supports recent findings that sun bears from Sundaland form a distinct clade from those in China and Indochina with Thailand possessing lineages from the three clades. Importantly, we demonstrate a more recent and clear genetic delineation between sun bears from the Malay Peninsula and Sabah indicating historical barriers to gene flow within the Sundaic region.

In nature, closely related species commonly display divergent mating behaviors, suggesting a central role for such traits in the origin of species. Elucidating the genetic basis of divergence in these traits is necessary to understand the evolutionary process leading to reproductive barriers and speciation. The rapidly speciating Hawaiian crickets of the genus Laupala provides an ideal system for dissecting the genetic basis of mating behavior divergence. In Laupala, closely related species differ markedly in male song pulse rate and female preference for pulse rate. These behaviors play an important role in determining mating patterns. Previous studies identified a genetic architecture consisting of numerous small to moderate effect loci causing interspecific differences in pulse rate and preference, including colocalizing pulse rate and preference QTL on linkage group one (LG1). To further interrogate these QTL, we conduct a fine mapping study using high-density SNP linkage maps. With improved statistical power and map resolution, we provide robust evidence for genetic coupling between song and preference, along with two additional pulse rate QTL on LG1, revealing a more resolved picture of the genetic architecture underlying mating behavior divergence. Our sequence-based genetic map, along with dramatically narrowed QTL confidence intervals, allowed us to annotate genes within the QTL regions and identify several exciting candidate genes underlying variation in pulse rate and preference divergence. Such knowledge suggests potential molecular mechanisms underlying the evolution of behavioral barriers.

In addition to including one of the most popular companion animals, species from the cat family Felidae serve as a powerful system for genetic analysis of inherited and infectious disease, as well as for the study of phenotypic evolution and speciation. Previous diploid-based genome assemblies for the domestic cat have served as the primary reference for genomic studies within the cat family. However, these versions suffered from poor resolution of complex and highly repetitive regions, with substantial amounts of unplaced sequence that is polymorphic or copy number variable. We sequenced the genome of a female F1 Bengal hybrid cat, the offspring of a domestic cat (Felis catus) x Asian leopard cat (Prionailurus bengalensis) cross, with PacBio long sequence reads and used Illumina sequence reads from the parents to phase >99.9% of the reads into the 2 species' haplotypes. De novo assembly of the phased reads produced highly continuous haploid genome assemblies for the domestic cat and Asian leopard cat, with contig N50 statistics exceeding 83 Mb for both genomes. Whole-genome alignments reveal the Felis and Prionailurus genomes are colinear, and the cytogenetic differences between the homologous F1 and E4 chromosomes represent a case of centromere repositioning in the absence of a chromosomal inversion. Both assemblies offer significant improvements over the previous domestic cat reference genome, with a 100% increase in contiguity and the capture of the vast majority of chromosome arms in 1 or 2 large contigs. We further demonstrated that comparably accurate F1 haplotype phasing can be achieved with members of the same species when one or both parents of the trio are not available. These novel genome resources will empower studies of feline precision medicine, adaptation, and speciation.

Selection acts upon genes linked together on chromosomes. This physical connection reduces the efficiency by which selection can act because, in the absence of sex, alleles must rise and fall together in frequency with the genome in which they are found. This selective interference underlies such phenomena as clonal interference and Muller's Ratchet and is broadly termed Hill-Robertson interference. In this review, I examine the potential for selective interference to account for the evolution and maintenance of sex, discussing the positive and negative evidence from both theoretical and empirical studies, and highlight the gaps that remain.

Organisms across the tree of life have complex life cycles that include both sexual and asexual reproduction or that are obligately asexual. These organisms include ecologically dominant species that structure many terrestrial and marine ecosystems, as well as many pathogens, pests, and invasive species. We must consider both the evolution and maintenance of these various reproductive modes and how these modes shape the genetic diversity, adaptive evolution, and ability to persist in the species that exhibit them. Thus, having a common framework is a key aspect of understanding the biodiversity that shapes our planet. In the 2019 AGA President's Symposium, Sex and Asex: The genetics of complex life cycles, researchers investigating a wide range of taxonomic models and using a variety of modes of investigation coalesced around a common theme-understanding not only how such complex life cycles may evolve, but how they are shaped by the evolutionary and ecological forces around them. In this introduction to the Special Issue from the symposium, we give an overview of some of the key ideas and areas of investigation (a common clonal lexicon, we might say) and introduce the breadth of work submitted by symposium participants.

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.

Genome studies of facultative sexual species, which can either reproduce sexually or asexually, are providing insight into the evolutionary consequences of mixed reproductive modes. It is currently unclear to what extent the evolutionary history of facultative sexuals' genomes can be approximated by the standard coalescent, and if a coalescent effective population size Ne exists. Here, I determine if and when these approximations can be made. When sex is frequent (occurring at a frequency much greater than 1/N per reproduction per generation, for N the actual population size), the underlying genealogy can be approximated by the standard coalescent, with a coalescent Ne≈N. When sex is very rare (at frequency much lower than 1/N), approximations for the pairwise coalescent time can be obtained, which is strongly influenced by the frequencies of sex and mitotic gene conversion, rather than N. However, these terms do not translate into a coalescent Ne. These results are used to discuss the best sampling strategies for investigating the evolutionary history of facultative sexual species.

Partially clonality is an incredibly common reproductive mode found across all the major eukaryotic lineages. Yet, population genetic theory is based on exclusive sexuality or exclusive asexuality, and partial clonality is often ignored. This is particularly true in haplodiplontic eukaryotes, including algae, ferns, mosses, and fungi, where somatic development occurs in both the haploid and diploid stages. Haplodiplontic life cycles are predicted to be correlated with asexuality, but tests of this prediction are rare. Moreover, there are unique consequences of having long-lived haploid and diploid stages in the same life cycle. For example, clonal processes uncouple the life cycle such that the repetition of the diploid stage via clonality leads to the loss of the haploid stage. Here, we surveyed the literature to find studies that had genotyped both haploid and diploid stages and recalculated population genetic summary metrics for seven red algae, one green alga, three brown algae, and three mosses. We compared these data to recent simulations that explicitly addressed the population genetic consequences of partial clonality in haplodiplontic life cycles. Not only was partial clonality found to act as a homogenizing force, but the combined effects of proportion of haploids, rate of clonality, and the relative strength of mutation versus genetic drift impacts the distributions of population genetic indices. We found remarkably similar patterns across commonly used population genetic metrics between our empirical and recent theoretical expectations. To facilitate future studies, we provide some recommendations for sampling and analyzing population genetic parameters for haplodiplontic taxa.

Sexual reproduction may pose myriad short-term costs to females. Despite these costs, sexual reproduction is near ubiquitous. Facultative parthenogenesis is theorized to mitigate some of the costs of sex, as individuals can participate in occasional sex to limit costs while obtaining many benefits. However, most theoretical models assume sexual reproduction is fixed following mating, with no possibility of clutches of mixed reproductive ontogeny. Therefore, we asked: if coercive males are present at high frequency in a population of facultative parthenogens, will their clutches be solely sexually produced, or will there be evidence of sexually and asexually-produced offspring? How will their offspring production compare to conspecifics in low-frequency male populations? We addressed our questions by collecting females and egg clutches of the facultatively parthenogenetic Opiliones species Leiobunum manubriatum and L. globosum. In L. manubriatum, females from populations with few males were not significantly more fecund than females from populations with higher male relative frequency, despite the potential release of the former from sexual conflict. We used 3 genotyping methods along with a custom set of DNA capture probes to reveal that offspring of L. manubriatum from these high male populations were primarily produced via asexual reproduction. This is surprising because sex ratios in these southern populations approach equality, increasing the probability for females to encounter mates and produce offspring sexually. We additionally found evidence for reproductive polymorphisms within populations. Rapid and accurate SNP genotyping data will continue to allow us to address broader evolutionary questions regarding the role of facultative reproductive modes in the maintenance of sex.

Through analyses of diverse microeukaryotes, we have previously argued that eukaryotic genomes are dynamic systems that rely on epigenetic mechanisms to distinguish germline (i.e., DNA to be inherited) from soma (i.e., DNA that undergoes polyploidization, genome rearrangement, etc.), even in the context of a single nucleus. Here, we extend these arguments by including two well-documented observations: (1) eukaryotic genomes interact frequently with mobile genetic elements (MGEs) like viruses and transposable elements (TEs), creating genetic conflict, and (2) epigenetic mechanisms regulate MGEs. Synthesis of these ideas leads to the hypothesis that genetic conflict with MGEs contributed to the evolution of a dynamic eukaryotic genome in the last eukaryotic common ancestor (LECA), and may have contributed to eukaryogenesis (i.e., may have been a driver in the evolution of FECA, the first eukaryotic common ancestor). Sex (i.e., meiosis) may have evolved within the context of the development of germline-soma distinctions in LECA, as this process resets the germline genome by regulating/eliminating somatic (i.e., polyploid, rearranged) genetic material. Our synthesis of these ideas expands on hypotheses of the origin of eukaryotes by integrating the roles of MGEs and epigenetics.

A growing body of research suggests that many clonal populations maintain genetic diversity even without occasional sexual reproduction. The purpose of our study was to document variation in single-nucleotide polymorphism (SNP) diversity, DNA content, and pathogen susceptibility in clonal lineages of the New Zealand freshwater snail, Potamopyrgus antipodarum. We studied snails that were collected from multiple field sites around a single lake (Lake Alexandrina), as well as isofemale clonal lineages that had been isolated and maintained in the laboratory. We used the kompetitive allele specific PCR (KASP) method to genotype our samples at 46 nuclear SNP sites, and we used flow cytometry to estimate DNA content. We found high levels of SNP diversity, both in our field samples and in our clonal laboratory lines. We also found evidence of high variation in DNA content among clones, even among clones with identical genotypes across all SNP sites. Controlled pathogen exposures of the laboratory populations revealed variation in susceptibility among distinct clonal genotypes, which was independent of DNA content. Taken together, these results show high levels of diversity among asexual snails, especially for DNA content, and they suggest rapid genome evolution in asexuals.

Despite decades of research, the evolution of sex remains an enigma in evolutionary biology. Typically, research addresses the costs of sex and asexuality to characterize the circumstances favoring one reproductive mode. Surprisingly few studies address the influence of common traits that are, in many organisms, obligately correlated with asexuality, including hybridization and polyploidy. These characteristics have substantial impacts on traits under selection. In particular, the fitness consequences of hybridization (i.e., reduced fitness due to interspecific reproductive isolation) will influence the evolution of sex. This may comprise a cost of either sex or asexuality due to the link between hybridity and asexuality. We examined reproductive isolation in the formation of de novo hybrid lineages between 2 widespread species in the ecological model system Boechera. Seventeen percent of 664 crosses produced F1 fruits, and only 10% of these were viable, suggesting that postmating prezygotic and postzygotic barriers inhibit hybrid success in this system. The postmating prezygotic barrier was asymmetrical, with 110 of 115 total F1 fruits produced when Boechera stricta acted as maternal parent. This asymmetry was confirmed in wild-collected lineages, using a chloroplast phylogeny of wild-collected B. stricta, Boechera retrofracta, and hybrids. We next compared fitness of F2 hybrids and selfed parental B. stricta lines, finding that F2 fitness was reduced by substantial hybrid sterility. Multiple reproductively isolating barriers influence the formation and fitness of hybrid lineages in the wild, and the costs of hybridization likely have profound impacts on the evolution of sex in the natural environment.

Ecological and evolutionary processes differ depending on how genetic diversity is organized in space. For clonal organisms, the organization of both genetic and genotypic diversity can influence the fitness effects of competition, the mating system, and reproductive mode, which are key drivers of life cycle evolution. Understanding how individual reproductive behavior contributes to population genetic structure is essential for disentangling these forces, particularly in species with complex and plastic life cycles. The widespread sea anemone, Diadumene lineata, exhibits temperature-dependent fission, which contributes to predictable variation in clonal rate along the Atlantic coast of the United States, part of its non-native range. Because warmer conditions lead to higher rates of clonality, we expected to find lower genotypic and genetic diversity in lower versus higher latitude populations. We developed primers for 11 microsatellite loci and genotyped 207 anemones collected from 8 sites ranging from Florida to Massachusetts. We found clonal influence at all sites, and as predicted, the largest clones were found at lower latitude sites. We also found genetic signatures of sex in the parts of the range where gametogenesis is most common. Evidence of sex outside the native range is novel for this species and provides insights into the dynamics of this successful invader. Our findings also illustrate challenges that partially clonal taxa pose for eco-evolutionary studies, such as difficulty sampling statistically robust numbers of genets and interpretating common population genetic metrics. For example, we found high among-locus variation in FIS, which makes the meaning of mean multilocus FIS unclear.

Studying alternative forms of reproduction in natural populations is of fundamental importance for understanding the costs and benefits of sex. Mayflies are one of the few animal groups where sexual reproduction co-occurs with different types of parthenogenesis, providing ideal conditions for identifying benefits of sex in natural populations. Here, we establish a catalog of all known mayfly species capable of reproducing by parthenogenesis, as well as species unable to do so. Overall, 1.8% of the described species reproduce parthenogenetically, which is an order of magnitude higher than reported in other animal groups. This frequency even reaches 47.8% if estimates are based on the number of studied rather than described mayfly species, as reproductive modes have thus far been studied in only 17 out of 42 families. We find that sex is a more successful strategy than parthenogenesis (associated with a higher hatching success of eggs), with a trade-off between the hatching success of parthenogenetic and sexual eggs. This means that improving the capacity for parthenogenesis may come at a cost for sexual reproduction. Such a trade-off can help explain why facultative parthenogenesis is extremely rare among animals despite its potential to combine the benefits of sexual and parthenogenetic reproduction. We argue that parthenogenesis is frequently selected in mayflies in spite of this probable trade-off because their typically low dispersal ability and short and fragile adult life may frequently generate situations of mate limitation in females. Mayflies are currently clearly underappreciated for understanding the benefits of sex under natural conditions.

In plants, long-distance dispersal is both attenuated and directed by specific movement vectors, including animals, wind, and/or water. Hence, movement vectors partly shape metapopulation genetic patterns that are, however, also influenced by other life-history traits such as clonal growth. We studied the relationship between area, isolation, plant-species richness, reproduction, and dispersal mechanisms with genetic diversity and divergence in 4 widespread wetland plant-species in a total of 20 island-like kettle-hole habitats surrounded by an intensive agricultural landscape. Our results showed that genetic parameters reflect the reproduction strategies with the highest genetic diversity being observed in the non-clonal, outcrossing Oenanthe aquatica compared to the clonal Lycopus europaeus, Typha latifolia, and Phragmites australis. Lycopus showed a positive relationship between genetic diversity and kettle-hole area, but a negative relationship with the number of neighboring kettle holes (less isolation). Genetic diversity increased with plant-species richness in the clonal species Phragmites and Lycopus; while it decreased in the non-clonal Oenanthe. Finally, genetic divergence and, therefore, connectivity differed between alternative dispersal strategies, where wind-dispersed Typha and Phragmites had a higher gene flow between the analyzed kettle holes compared with the insect-pollinated, hydrochorous Lycopus and Oenanthe. Our study provides information on genetic patterns related to reproduction and dispersal mechanisms of 4 common wetland species contributing to the understanding of the functioning of plant metacommunities occurring in kettle holes embedded in agricultural landscapes.

Partial clonality is known to affect the genetic composition and evolutionary trajectory of diplontic (single, free-living diploid stage) populations. However, many partially clonal eukaryotes exhibit life cycles in which somatic development occurs in both haploid and diploid individuals (haplodiplontic life cycles). Here, we studied how haplodiplontic life cycles and partial clonality structurally constrain, as immutable parameters, the reshuffling of genetic diversity and its dynamics in populations over generations. We assessed the distribution of common population genetic indices at different proportions of haploids, rates of clonality, mutation rates, and sampling efforts. Our results showed that haplodiplontic life cycles alone in finite populations affect effective population sizes and the ranges of distributions of population genetic indices. With nonoverlapping generations, haplodiplonty allowed the evolution of 2 temporal genetic pools that may diverge in sympatry due to genetic drift under full sexuality and clonality. Partial clonality in these life cycles acted as a homogenizing force between those 2 pools. Moreover, the combined effects of proportion of haploids, rate of clonality, and the relative strength of mutation versus genetic drift impacts the distributions of population genetics indices, rendering it difficult to transpose and use knowledge accumulated from diplontic or haplontic species. Finally, we conclude by providing recommendations for sampling and analyzing the population genetics of partially clonal haplodiplontic taxa.