Deleterious mutations and the evolution of sexual reproduction

Evolvability: filling the explanatory gap between adaptedness and the long-term mathematical conception of fitness

The new foundation for the propensity interpretation of fitness (PIF), developed by Pence and Ramsey (Br J Philos Sci 64:851-881, 2013), describes fitness as a probability distribution that encompasses all possible daughter populations to which the organism may give rise, including daughter populations in which traits might change and the possible environments that members of the daughter populations might encounter. This long-term definition of fitness is general enough to avoid counterexamples faced by previous mathematical conceptions of PIF. However, there seem to be downsides to its generality: the ecological role of fitness involves describing the degree of adaptedness between an organism and the specific environment it inhabits. When all possible changes in traits and all possible environments that a daughter population may encounter are included in the concept, it becomes difficult to see how fitness can fulfill this role. In this paper, we argue that this is a feature of Pence and Ramsey's view rather than a bug: long-term fitness accommodates evolvability considerations, which concern the role that variation plays in evolutionary processes. Building on the foundations, we show that Pence and Ramsey's fitness-F-can be partitioned into fourths: adaptedness, robustness of adaptedness, and two facets of evolvability. Conceptualizing these last three components forces us to consider the role played by grains of description of both organisms and the environment when thinking about long-term fitness. They track the possibility that there could be a change in type in a daughter population as a way of responding to environmental challenges, or that the type persists in the face of novel environments. We argue that these components are just as salient as adaptedness for long-term fitness. Together, this decomposition of F provides a more accurate picture of the factors involved in long-term evolutionary success.

Understanding the Genetic Basis of Variation in Meiotic Recombination: Past, Present, and Future

Meiotic recombination is a fundamental feature of sexually reproducing species. It is often required for proper chromosome segregation and plays important role in adaptation and the maintenance of genetic diversity. The molecular mechanisms of recombination are remarkably conserved across eukaryotes, yet meiotic genes and proteins show substantial variation in their sequence and function, even between closely related species. Furthermore, the rate and distribution of recombination shows a huge diversity within and between chromosomes, individuals, sexes, populations, and species. This variation has implications for many molecular and evolutionary processes, yet how and why this diversity has evolved is not well understood. A key step in understanding trait evolution is to determine its genetic basis-that is, the number, effect sizes, and distribution of loci underpinning variation. In this perspective, I discuss past and current knowledge on the genetic basis of variation in recombination rate and distribution, explore its evolutionary implications, and present open questions for future research.

Limits to species' range: the tension between local and global adaptation

We know that heritable variation is abundant, and that selection causes all but the smallest populations to rapidly shift beyond their original trait distribution. So then, what limits the range of a species? There are physical constraints and also population genetic limits to the effectiveness of selection, ultimately set by population size. Global adaptation, where the same genotype is favoured over the whole range, is most efficient when based on a multitude of weakly selected alleles and is effective even when local demes are small, provided that there is some gene flow. In contrast, local adaptation is sensitive to gene flow and may require alleles with substantial effect. How can populations combine the advantages of large effective size with the ability to specialise into local niches? To what extent does reproductive isolation help resolve this tension? I address these questions using eco-evolutionary models of polygenic adaptation, contrasting discrete demes with continuousspace.

Apomixis and the paradox of sex in plants

BACKGROUND

The predominance of sex in eukaryotes, despite the high costs of meiosis and mating, remains an evolutionary enigma. Many theories have been proposed, none of them being conclusive on its own, and they are, in part, not well applicable to land plants. Sexual reproduction is obligate in embryophytes for the great majority of species.

SCOPE

This review compares the main forms of sexual and asexual reproduction in ferns and angiosperms, based on the generation cycling of sporophyte and gametophyte (leaving vegetative propagation aside). The benefits of sexual reproduction for maintenance of genomic integrity in comparison to asexuality are discussed in the light of developmental, evolutionary, genetic and phylogenetic studies.

CONCLUSIONS

Asexual reproduction represents modifications of the sexual pathway, with various forms of facultative sexuality. For sexual land plants, meiosis provides direct DNA repair mechanisms for oxidative damage in reproductive tissues. The ploidy alternations of meiosis-syngamy cycles and prolonged multicellular stages in the haploid phase in the gametophytes provide a high efficiency of purifying selection against recessive deleterious mutations. Asexual lineages might buffer effects of such mutations via polyploidy and can purge the mutational load via facultative sexuality. The role of organelle-nuclear genome compatibility for maintenance of genome integrity is not well understood. In plants in general, the costs of mating are low because of predominant hermaphroditism. Phylogenetic patterns in the archaeplastid clade suggest that high frequencies of sexuality in land plants are concomitant with a stepwise increase of intrinsic and extrinsic stress factors. Furthermore, expansion of genome size in land plants would increase the potential mutational load. Sexual reproduction appears to be essential for keeping long-term genomic integrity, and only rare combinations of extrinsic and intrinsic factors allow for shifts to asexuality.

Population genomics and epigenomics of Spirodela polyrhiza provide insights into the evolution of facultative asexuality

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

Diversity and Convergence of Sex-Determination Mechanisms in Teleost Fish

Sexual reproduction is prevalent across diverse taxa. However, sex-determination mechanisms are so diverse that even closely related species often differ in sex-determination systems. Teleost fish is a taxonomic group with frequent turnovers of sex-determining mechanisms and thus provides us with great opportunities to investigate the molecular and evolutionary mechanisms underlying the turnover of sex-determining systems. Here, we compile recent studies on the diversity of sex-determination mechanisms in fish. We demonstrate that genes in the TGF-β signaling pathway are frequently used for master sex-determining (MSD) genes. MSD genes arise via two main mechanisms, duplication-and-transposition and allelic mutations, with a few exceptions. We also demonstrate that temperature influences sex determination in many fish species, even those with sex chromosomes, with higher temperatures inducing differentiation into males in most cases. Finally, we review theoretical models for the turnover of sex-determining mechanisms and discuss what questions remain elusive.

Mysterious Bamboo flowering phenomenon: A literature review and new perspectives

Bamboo, a globally distributed non-timber forest resource, plays a critical role in local ecosystems and economies. Despite its significance, the understanding of bamboo's long and unpredictable flowering cycles remains limited. Our bibliometric analysis of bamboo flowering-related literature from the Web of Science database reveals an initial focus on regeneration studies, with a recent trend shifting towards microscopic and molecular perspectives. Furthermore, our narrative review emphasizes the importance of considering factors such as the proportion of flowering culms and the duration of flowering in classifying bamboo flowering phenomena. While numerous studies have endorsed the predator saturation hypothesis as a suitable explanation for the synchronicity of bamboo flowering, no existing theory explains bamboo's prolonged flowering cycles. We propose a new natural selection hypothesis as a potential explanation for these extraordinary cycles, underscoring the need for further research in this area. Despite the substantial volume of data accumulated on bamboo flowering, these resources have not been fully exploited in recent research. Future studies would benefit from more comprehensive data collection methods, encompassing field observations, satellite remote sensing data, and omics data. The convergence of traditional ecological studies with molecular techniques may pave the way for significant advancements in bamboo flowering research.

Phased genome assemblies reveal haplotype-specific genetic load in the critically endangered Chinese Bahaba (Teleostei, Sciaenidae)

While haplotype-specific genetic load shapes the evolutionary trajectory of natural and captive populations, mixed-haplotype assembly and genotyping hindered its characterization in diploids. Herein, we produced two phased genome assemblies of the critically endangered fish Chinese Bahaba (Bahaba taipingensis, Sciaenidae, Teleostei) and resequenced 20 whole genomes to quantify population genetic load at a haplotype level. We identified frame-shifting variants as the most deleterious type, followed by mutations in the 5'-UTR, 3'-UTR and missense mutations at conserved amino acids. Phased haplotypes revealed gene deletions and high-impact deleterious variants. We estimated ~1.12% of genes missing or interrupted per haplotype, with a significant overlap of disrupted genes (30.35%) between haplotype sets. Relative proportions of deleterious variant categories differed significantly between haplotypes. Simulations suggested that purifying selection struggled to purge slightly deleterious genetic load in captive breeding compared to genotyping interventions, and that higher inter-haplotypic variance of genetic load predicted more efficient purging by artificial selection. Combining the knowledge of haplotype-resolved genetic load with predictive modelling will be immensely useful for understanding the evolution of deleterious variants and guiding conservation planning.

Unzipped chromosome-level genomes reveal allopolyploid nematode origin pattern as unreduced gamete hybridization

The formation and consequences of polyploidization in animals with clonal reproduction remain largely unknown. Clade I root-knot nematodes (RKNs), characterized by parthenogenesis and allopolyploidy, show a widespread geographical distribution and extensive agricultural destruction. Here, we generated 4 unzipped polyploid RKN genomes and identified a putative novel alternative telomeric element. Then we reconstructed 4 chromosome-level assemblies and resolved their genome structures as AAB for triploid and AABB for tetraploid. The phylogeny of subgenomes revealed polyploid RKN origin patterns as hybridization between haploid and unreduced gametes. We also observed extensive chromosomal fusions and homologous gene expression decrease after polyploidization, which might offset the disadvantages of clonal reproduction and increase fitness in polyploid RKNs. Our results reveal a rare pathway of polyploidization in parthenogenic polyploid animals and provide a large number of high-precision genetic resources that could be used for RKN prevention and control.

Association of a single amino acid replacement with dorsal pigmentation in a lizard from the Qinghai-Tibetan Plateau

Reptiles can evolve adaptive colors in different environments, but relatively little is known about the genetic mechanisms. Here, we identified the MC1R gene and its association with intraspecific color variation in the lizard Phrynocephalus erythrurus. Analysis of the MC1R sequence in 143 individuals from dark South Qiangtang Plateau (SQP) and light North Qiangtang plateau (NQP) populations, revealed two amino acid sites that showed significant differences in frequency between two areas. One SNP, corresponding to Glu183Lys residue, was found to be a highly significant outlier and differentially fixed for SQP and NQP populations. This residue is located in an extracellular area in the second small extracellular loop within the secondary structure of MC1R, which represents an "attachment pocket" part of the 3D structure. Cytological expression of MC1R alleles with the Glu183Lys replacement showed a 39 % increase in intracellular agonist-induced cyclic AMP levels and a 23.18 % greater cell surface expression of MC1R protein in the SQP relative to the NQP allele. Further in silico 3D modeling and in vitro binding experiments indicated a higher MC1R-α-MSH binding for the SQP allele, and elevated melanin synthesis. We provide an overview of how a single amino acid replacement leads to fundamental changes in MC1R function, and hence shapes variation in dorsal pigmentation in lizards from different environments.

Mutations, sex, and genetic diversity: New arguments for partially directed evolution

A review of the theories of the existence of sexes, genetic diversity, and the distribution of mutations among organisms shows that all these concepts are not a product of random evolution and cannot be explained within the framework of Darwinism. Most mutations are the result of the genome acting on itself. This is an organized process that is implemented very differently in different species, in different places in the genome. Because of the fact that it is not random, this process must be directed and regulated, albeit with complex and not fully understood laws. This means that an additional reason must be included in order to model such mutations during evolution. The assumption of directionality must not only be explicitly included in evolutionary theory but must also occupy a central place in it. In this study an updated model of partially directed evolution is constructed, which is capable of qualitatively explaining the indicated features of evolution. Experiments are described that can confirm or disprove the proposed model.

Uncovering the Male Presence in Parthenogenetic Marchalina hellenica (Hemiptera: Marchalinidae): Insights into Its mtDNA Divergence and Reproduction Strategy

Marchalina hellenica (Hemiptera: Marchalinidae), an endemic species in Greece and Turkey, is a major contributor to the annual honey production in its native range. However, in the areas that it invades, lacking natural enemies, it has detrimental effects on pine trees and potentially contributes to tree mortality. Although it was originally reported as thelytokous, males were later reported in Turkey and on several of the islands of Greece. To further disambiguate the exact parthenogenetic reproduction strategy of M. hellenica, we studied the emergence pattern of male individuals in Greece for two consecutive years (2021 and 2022). Furthermore, we examined the genetic variation among 15 geographically distant populations of M. hellenica in Greece using a mitochondrial DNA marker and compared the results with data from Turkey. The findings of this study document the existence of an additional M. hellenica population in its native range that repeatedly produces males, apart from the areas of Greece and Turkey in which they were initially reported, suggesting that males play a major, so far unknown role in the reproduction of this species. The populations in Greece and Turkey exhibited a strong genetic affinity, while human-aided dispersal seems to have obscured the genetic pattern acquired.

What Has Genomics Taught An Evolutionary Biologist?

Genomics, an interdisciplinary field of biology on the structure, function, and evolution of genomes, has revolutionized many subdisciplines of life sciences, including my field of evolutionary biology, by supplying huge data, bringing high-throughput technologies, and offering a new approach to biology. In this review, I describe what I have learned from genomics and highlight the fundamental knowledge and mechanistic insights gained. I focus on three broad topics that are central to evolutionary biology and beyond-variation, interaction, and selection-and use primarily my own research and study subjects as examples. In the next decade or two, I expect that the most important contributions of genomics to evolutionary biology will be to provide genome sequences of nearly all known species on Earth, facilitate high-throughput phenotyping of natural variants and systematically constructed mutants for mapping genotype-phenotype-fitness landscapes, and assist the determination of causality in evolutionary processes using experimental evolution.

Conflicting effects of recombination on the evolvability and robustness in neutrally evolving populations

Understanding the benefits and costs of recombination under different scenarios of evolutionary adaptation remains an open problem for theoretical and experimental research. In this study, we focus on finite populations evolving on neutral networks comprising viable and unfit genotypes. We provide a comprehensive overview of the effects of recombination by jointly considering different measures of evolvability and mutational robustness over a broad parameter range, such that many evolutionary regimes are covered. We find that several of these measures vary non-monotonically with the rates of mutation and recombination. Moreover, the presence of unfit genotypes that introduce inhomogeneities in the network of viable states qualitatively alters the effects of recombination. We conclude that conflicting trends induced by recombination can be explained by an emerging trade-off between evolvability on the one hand, and mutational robustness on the other. Finally, we discuss how different implementations of the recombination scheme in theoretical models can affect the observed dependence on recombination rate through a coupling between recombination and genetic drift.

Accumulation and maintenance of information in evolution

Selection accumulates information in the genome-it guides stochastically evolving populations toward states (genotype frequencies) that would be unlikely under neutrality. This can be quantified as the Kullback-Leibler (KL) divergence between the actual distribution of genotype frequencies and the corresponding neutral distribution. First, we show that this population-level information sets an upper bound on the information at the level of genotype and phenotype, limiting how precisely they can be specified by selection. Next, we study how the accumulation and maintenance of information is limited by the cost of selection, measured as the genetic load or the relative fitness variance, both of which we connect to the control-theoretic KL cost of control. The information accumulation rate is upper bounded by the population size times the cost of selection. This bound is very general, and applies across models (Wright-Fisher, Moran, diffusion) and to arbitrary forms of selection, mutation, and recombination. Finally, the cost of maintaining information depends on how it is encoded: Specifying a single allele out of two is expensive, but one bit encoded among many weakly specified loci (as in a polygenic trait) is cheap.

The origination events of gametic sexual reproduction and anisogamy

The evolution of gametic sex (meiosis and fertilization) and subsequent transition from isogamy (fusion between two equal-sized gametes) to anisogamy (dimorphism into eggs and sperm, namely, females and males) is one of the largest enigmas of evolutionary biology. Meiosis entails genome-dilution cost and anisogamy entails male-production cost. Despite much progress has been made for the maintenance mechanisms of sex, its origination events under such “twofold cost of sex” are still unsolved. Here, we posit two hypothetical scenarios as follows: the “Seesaw Effect” hypothesizes that automictic selfing between isogametes effectively purged deleterious mutations from an organism’s lineage and simultaneously fixed the sex-controlling allele and all other loci (no genome-dilution cost raised). The high relatedness among homoeologous cell colonies led to multicellularization. The “inflated isogamy” hypothesizes that multicellularity increased the reproductive investment of both mates, resulting in excessively large isogametes. This redundancy induced cheating of one sex (evolving to male) to reduce gamete size. However, the other sex (evolving to female) allowed this cheat because her cost did not change. Therefore, anisogamy originated as a kind of commensalism but turned into beneficial for females because it solved the gamete limitation problem inherent to isogamy. Thus, smooth transition to anisogamy had been attained.

Insect herbivores drive sex allocation in angiosperm flowers

Why sex has evolved and is maintained is an open question in evolutionary biology. The Red Queen hypothesis predicts that host lineages subjected to more intense parasite pressure should invest more in sexual reproduction to continuously create novel defences against their rapidly evolving natural enemies. In this comparative study across the angiosperms, we show that hermaphrodite plant species associated with higher species richness of insect herbivores evolved flowers with higher biomass allocation towards the male sex, an indication of their greater outcrossing effort. This pattern remained robust after controlling for key vegetative, reproductive and biogeographical traits, suggesting that long-term herbivory pressure is a key factor driving the selfing-outcrossing gradient of higher plants. Although flower evolution is frequently associated with mutualistic pollinators, our findings support the Red Queen hypothesis and suggest that insect herbivores drive the sexual strategies of flowering plants and their genetic diversity.

The advantage of sex: Reinserting fluctuating selection in the pluralist approach

The advantage of sex, and its fixation in some clades and species all over the eukaryote tree of life, is considered an evolutionary enigma, especially regarding its assumed two-fold cost. Several likely hypotheses have been proposed such as (1) a better response to the negative frequency-dependent selection imposed by the “Red Queen” hypothesis; (2) the competition between siblings induced by the Tangled Bank hypothesis; (3) the existence of genetic and of (4) ecological factors that can diminish the cost of sex to less than the standard assumed two-fold; and (5) a better maintenance of genetic diversity and its resulting phenotypic variation, providing a selective advantage in randomly fluctuating environments. While these hypotheses have mostly been studied separately, they can also act simultaneously. This was advocated by several studies which presented a pluralist point of view. Only three among the five causes cited above were considered yet in such a framework: the Red Queen hypothesis, the Tangled Bank and the genetic factors lowering the cost of sex. We thus simulated the evolution of a finite mutating population undergoing negative frequency-dependent selection on phenotypes and a two-fold (or less) cost of sexuality, experiencing randomly fluctuating selection along generations. The individuals inherited their reproductive modes, either clonal or sexual. We found that exclusive sexuality begins to fix in populations exposed to environmental variation that exceeds the width of one ecological niche (twice the standard deviation of a Gaussian response to environment). This threshold was lowered by increasing negative frequency-dependent selection and when reducing the two-fold cost of sex. It contributes advocating that the different processes involved in a short-term advantage of sex and recombination can act in combination to favor the fixation of sexual reproduction in populations.

Mitochondrial Genomes of Two Asexual Trichogramma (Hymenoptera: Trichogrammatidae) Strains and Comparison with Their Sexual Relatives

Simple Summary

Sexual reproduction is dominant in animals, while asexual lineages are rare and evolutionarily short-lived. However, sexual reproduction has substantial costs, such as male production, inputs to courtship and mating, increased risk of predator exposure, and sexually transmitted diseases. A large body of theories has been proposed to explain the paradox of sex. One favored explanation is that asexuals are more likely to accumulate a greater number of deleterious mutations, known as Muller’s ratchet. Trichogramma is a genus of egg parasitoid wasps and is widely used as a biological control agent for agricultural and forest pests. With asexual lineages in at least 16 species, Trichogramma provides an excellent model to investigate the causes and consequences of asexual reproduction. In this study, we sequenced and assembled two asexual Trichogramma mitogenomes, representing two divergent origins of asexual reproduction. The asexual T. pretiosum is induced by the endosymbiont Wolbachia, while T. cacoeciae presumably originates from interspecific hybridization. To test Muller’s ratchet hypothesis, we compared these two asexual mitogenomes with their sexual relatives and found no association between asexual reproduction and mutation accumulation. This study provides a basis for further investigation into mitochondrial evolution and asexual reproduction in Trichogramma.

Abstract

Despite its substantial costs, sexual reproduction dominates in animals. One popular explanation for the paradox of sex is that asexual reproduction is more likely to accumulate deleterious mutations than sexual reproduction. To test this hypothesis, we compared the mitogenomes of two asexual wasp strains, Trichogramma cacoeciae and T. pretiosum, to their sexual relatives. These two asexual strains represent two different transition mechanisms in Trichogramma from sexual to asexual reproduction. Asexual T. pretiosum is induced by Wolbachia, while T. cacoeciae presumably originated from interspecific hybridization. We sequenced and assembled complete mitochondrial genomes of asexual T. cacoeciae and T. pretiosum. Compared to four sexual relatives, we found no evidence of higher mutation accumulation in asexual Trichogramma mitogenomes than in their sexual relatives. We also did not detect any relaxed selection in asexual Trichogramma mitogenomes. In contrast, the intensified selection was detected in Nad1 and Nad4 of the asexual T. pretiosum mitogenome, suggesting more purifying selection. In summary, no higher mitochondrial mutation accumulation was detected in these two asexual Trichogramma strains. This study provides a basis for further investigating mitochondrial evolution and asexual reproduction in Trichogramma.

The impact of genetic modifiers on variation in germline mutation rates within and among human populations

Mutation rates and spectra differ among human populations. Here, we examine whether this variation could be explained by evolution at mutation modifiers. To this end, we consider genetic modifier sites at which mutations, "mutator alleles", increase genome-wide mutation rates and model their evolution under purifying selection due to the additional deleterious mutations that they cause, genetic drift, and demographic processes. We solve the model analytically for a constant population size and characterize how evolution at modifier sites impacts variation in mutation rates within and among populations. We then use simulations to study the effects of modifier sites under a plausible demographic model for Africans and Europeans. When comparing populations that evolve independently, weakly selected modifier sites (2Nes ≈ 1), which evolve slowly, contribute the most to variation in mutation rates. In contrast, when populations recently split from a common ancestral population, strongly selected modifier sites (2Nes » 1), which evolve rapidly, contribute the most to variation between them. Moreover, a modest number of modifier sites (e.g., 10 per mutation type in the standard classification into 96 types) subject to moderate to strong selection (2Nes > 1) could account for the variation in mutation rates observed among human populations. If such modifier sites indeed underlie differences among populations, they should also cause variation in mutation rates within populations and their effects should be detectable in pedigree studies.

Genomic Analyses of Wild and Cultivated Bacanora Agave (Agave angustifolia var. pacifica) Reveal Inbreeding, Few Signs of Cultivation History and Shallow Population Structure

Due to the recent increase in demand for agave-based beverages, many wild agave populations have experienced rapid decline and fragmentation, whereas cultivated plants are now managed at monocultural plantations, in some cases involving clonal propagation. We examined the relative effect of migration, genetic drift, natural selection and human activities on the genetic repertoire of Agave angustifolia var. pacifica, an agave used for bacanora (an alcoholic spirit similar to tequila) production in northwestern Mexico. We sampled 34 wild and cultivated sites and used over eleven thousand genome-wide SNPs. We found shallow genetic structure among wild samples, although we detected differentiation between coastal and inland sites. Surprisingly, no differentiation was found between cultivated and wild populations. Moreover, we detected moderate inbreeding (F_IS ~ 0.13) and similar levels of genomic diversity in wild and cultivated agaves. Nevertheless, the cultivated plants had almost no private alleles and presented evidence of clonality. The overall low genetic structure in A. angustifolia var. pacifica is apparently the result of high dispersibility promoted by pollinators and the possibility of clonal reproduction. Incipient cultivation history and reliance on wild seeds and plants are probably responsible for the observed patterns of high genetic connectivity and considerable diversity in cultivated samples.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

This essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

Domestication reprogrammed the budding yeast life cycle

Domestication of plants and animals is the foundation for feeding the world human population but can profoundly alter the biology of the domesticated species. Here we investigated the effect of domestication on one of our prime model organisms, the yeast Saccharomyces cerevisiae, at a species-wide level. We tracked the capacity for sexual and asexual reproduction and the chronological life span across a global collection of 1,011 genome-sequenced yeast isolates and found a remarkable dichotomy between domesticated and wild strains. Domestication had systematically enhanced fermentative and reduced respiratory asexual growth, altered the tolerance to many stresses and abolished or impaired the sexual life cycle. The chronological life span remained largely unaffected by domestication and was instead dictated by clade-specific evolution. We traced the genetic origins of the yeast domestication syndrome using genome-wide association analysis and genetic engineering and disclosed causative effects of aneuploidy, gene presence/absence variations, copy number variations and single-nucleotide polymorphisms. Overall, we propose domestication to be the most dramatic event in budding yeast evolution, raising questions about how much domestication has distorted our understanding of the natural biology of this key model species.

Recombination resolves the cost of horizontal gene transfer in experimental populations of Helicobacter pylori

Horizontal gene transfer (HGT)—the transfer of DNA between lineages—is responsible for a large proportion of the genetic variation that contributes to evolution in microbial populations. While HGT can bring beneficial genetic innovation, the transfer of DNA from other species or strains can also have deleterious effects. In this study, we evolve populations of the bacteria Helicobacter pylori and use DNA sequencing to identify over 40,000 genetic variants transferred by HGT. We measure the cost of many of these and find that both strongly beneficial mutations and deleterious mutations are genetic variants transferred by natural transformation. Importantly, we also show how recombination that separates linked beneficial and deleterious mutations resolves the cost of HGT.

Horizontal gene transfer (HGT) is important for microbial evolution, yet we know little about the fitness effects and dynamics of horizontally transferred genetic variants. In this study, we evolve laboratory populations of Helicobacter pylori, which take up DNA from their environment by natural transformation, and measure the fitness effects of thousands of transferred genetic variants. We find that natural transformation increases the rate of adaptation but comes at the cost of significant genetic load. We show that this cost is circumvented by recombination, which increases the efficiency of selection by decoupling deleterious and beneficial genetic variants. Our results show that adaptation with HGT, pervasive in natural microbial populations, is shaped by a combination of selection, recombination, and genetic drift not accounted for in existing models of evolution.

Genetic load and extinction in peripheral populations: the roles of migration, drift and demographic stochasticity

We analyse how migration from a large mainland influences genetic load and population numbers on an island, in a scenario where fitness-affecting variants are unconditionally deleterious, and where numbers decline with increasing load. Our analysis shows that migration can have qualitatively different effects, depending on the total mutation target and fitness effects of deleterious variants. In particular, we find that populations exhibit a genetic Allee effect across a wide range of parameter combinations, when variants are partially recessive, cycling between low-load (large-population) and high-load (sink) states. Increased migration reduces load in the sink state (by increasing heterozygosity) but further inflates load in the large-population state (by hindering purging). We identify various critical parameter thresholds at which one or other stable state collapses, and discuss how these thresholds are influenced by the genetic versus demographic effects of migration. Our analysis is based on a 'semi-deterministic' analysis, which accounts for genetic drift but neglects demographic stochasticity. We also compare against simulations which account for both demographic stochasticity and drift. Our results clarify the importance of gene flow as a key determinant of extinction risk in peripheral populations, even in the absence of ecological gradients. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

Reproduction dynamics of planktonic microbial eukaryotes in the open ocean

Understanding the biology of reproduction of an organismal lineage is important for retracing key evolutionary processes, yet gaining detailed insights often poses major challenges. Planktonic Foraminifera are globally distributed marine microbial eukaryotes and important contributors to the global carbon cycle. They cannot routinely be cultured under laboratory conditions across generations, and thus details of their life cycle remain incomplete. The production of flagellated gametes has long been taken as an indication of exclusively sexual reproduction, but recent research suggests the existence of an additional asexual generation in the life cycle. To gain a better understanding of the reproductive biology of planktonic Foraminifera, we applied a dynamic, individual-based modelling approach with parameters based on laboratory and field observations to test if sexual reproduction is sufficient for maintaining viable populations. We show that temporal synchronization and potentially spatial concentration of gamete release seems inevitable for maintenance of the population under sexual reproduction. We hypothesize that sexual reproduction is likely beneficial during the adaptation to new environments, while population sustenance in stable environments can be ensured through asexual reproduction.

Identification of a major QTL and genome-wide epistatic interactions for single vs. paired spikelets in a maize-teosinte F2 population

Maize ear carries paired spikelets, whereas the ear of its wild ancestor, teosinte, bears single spikelets. However, little is known about the genetic basis of the processes of transformation of single spikelets in teosinte ear to paired spikelets in maize ear. In this study, a two-ranked, paired-spikelets primitive maize and a two-ranked, single-spikelet teosinte were utilized to develop an F2 population, and quantitative trait locus (loci) (QTL) mapping for single vs. paired spikelets (PEDS) was performed. One major QTL (qPEDS3.1) for PEDS located on chromosome 3S was identified in the 162 F2 plants using the inclusive composite interval mapping of additive (ICIM-ADD) module, explaining 23.79% of the phenotypic variance. Out of the 409 F2 plants, 43 plants with PEDS = 0% and 43 plants with PEDS > 20% were selected for selective genotyping, and the QTL (qPEDS3.1) was detected again. Moreover, the QTL (qPEDS3.1) was validated in three environments, which explained 31.05%, 38.94%, and 23.16% of the phenotypic variance, respectively. In addition, 50 epistatic QTLs were detected in the 162 F2 plants using the two-locus epistatic QTL (ICIM-EPI) module; they were distributed on all 10 chromosomes and explained 94.40% of the total phenotypic variance. The results contribute to a better understanding of the genetic basis of domestication of paired spikelets and provide a genetic resource for future map-based cloning; in addition, the systematic dissection of epistatic interactions underlies a theoretical framework for overcoming epistatic effects on QTL fine mapping.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01276-x.

A PARTHENOGENESIS allele from apomictic dandelion can induce egg cell division without fertilization in lettuce

Apomixis, the clonal formation of seeds, is a rare yet widely distributed trait in flowering plants. We have isolated the PARTHENOGENESIS (PAR) gene from apomictic dandelion that triggers embryo development in unfertilized egg cells. PAR encodes a K2-2 zinc finger, EAR-domain protein. Unlike the recessive sexual alleles, the dominant PAR allele is expressed in egg cells and has a miniature inverted-repeat transposable element (MITE) transposon insertion in the promoter. The MITE-containing promoter can invoke a homologous gene from sexual lettuce to complement dandelion LOSS OF PARTHENOGENESIS mutants. A similar MITE is also present in the promoter of the PAR gene in apomictic forms of hawkweed, suggesting a case of parallel evolution. Heterologous expression of dandelion PAR in lettuce egg cells induced haploid embryo-like structures in the absence of fertilization. Sexual PAR alleles are expressed in pollen, suggesting that the gene product releases a block on embryogenesis after fertilization in sexual species while in apomictic species PAR expression triggers embryogenesis in the absence of fertilization.

Genomic Signature of Sexual Reproduction in the Bdelloid Rotifer Macrotrachella quadricornifera

Bdelloid rotifers, common freshwater invertebrates of ancient origin and worldwide distribution have long been thought to be entirely asexual, being the principal exception to the view that in eukaryotes the loss of sex leads to early extinction. That bdelloids are facultatively sexual is shown by a study of allele sharing within a group of closely related bdelloids of the species Macrotrachella quadricornifera, supporting the view that sexual reproduction is essential for long-term success in all eukaryotes.

Intraspecific Variation in the Rates of Mutations Causing Structural Variation in Daphnia magna

Mutations that cause structural variation are important sources of genetic variation upon which other evolutionary forces can act, however, they are difficult to observe and therefore few direct estimates of their rate and spectrum are available. Understanding mutation rate evolution, however, requires adding to the limited number of species for which direct estimates are available, quantifying levels of intraspecific variation in mutation rates, and assessing whether rate estimates co-vary across types of mutation. Here, we report structural variation-causing mutation rates (svcMRs) for six categories of mutations (short insertions and deletions, long deletions and duplications, and deletions and duplications at copy number variable sites) from nine genotypes of Daphnia magna collected from three populations in Finland, Germany, and Israel using a mutation accumulation approach. Based on whole-genome sequence data and validated using simulations, we find svcMRs are high (two orders of magnitude higher than base substitution mutation rates measured in the same lineages), highly variable among populations, and uncorrelated across categories of mutation. Furthermore, to assess the impact of scvMRs on the genome, we calculated rates while adjusting for the lengths of events and ran simulations to determine if the mutations occur in genic regions more or less frequently than expected by chance. Our results pose a challenge to most prevailing theories aimed at explaining the evolution of the mutation rate, underscoring the importance of obtaining additional mutation rate estimates in more genotypes, for more types of mutation, in more species, in order to improve our future understanding of mutation rates, their variation, and their evolution.

Precise measurement of the fitness effects of spontaneous mutations by droplet digital PCR in Burkholderia cenocepacia

Understanding how mutations affect survivability is a key component to knowing how organisms and complex traits evolve. However, most mutations have a minor effect on fitness and these effects are difficult to resolve using traditional molecular techniques. Therefore, there is a dire need for more accurate and precise fitness measurements methods. Here, we measured the fitness effects in Burkholderia cenocepacia HI2424 mutation accumulation (MA) lines using droplet-digital polymerase chain reaction (ddPCR). Overall, the fitness measurements from ddPCR-MA are correlated positively with fitness measurements derived from traditional phenotypic marker assays (r = 0.297, P = 0.05), but showed some differences. First, ddPCR had significantly lower measurement variance in fitness (F = 3.78, P < 2.6 × 10-13) in control experiments. Second, the mean fitness from ddPCR-MA measurements were significantly lower than phenotypic marker assays (-0.0041 vs -0.0071, P = 0.006). Consistent with phenotypic marker assays, ddPCR-MA measurements observed multiple (27/43) lineages that significantly deviated from mean fitness, suggesting that a majority of the mutations are neutral or slightly deleterious and intermixed with a few mutations that have extremely large effects. Of these mutations, we found a significant excess of mutations within DNA excinuclease and Lys R transcriptional regulators that have extreme deleterious and beneficial effects, indicating that modifications to transcription and replication may have a strong effect on organismal fitness. This study demonstrates the power of ddPCR as a ubiquitous method for high-throughput fitness measurements in both DNA- and RNA-based organisms regardless of cell type or physiology.

Long-term exhaustion of the inbreeding load in Drosophila melanogaster

Inbreeding depression, the decline in fitness of inbred individuals, is a ubiquitous phenomenon of great relevance in evolutionary biology and in the fields of animal and plant breeding and conservation. Inbreeding depression is due to the expression of recessive deleterious alleles that are concealed in heterozygous state in noninbred individuals, the so-called inbreeding load. Genetic purging reduces inbreeding depression by removing these alleles when expressed in homozygosis due to inbreeding. It is generally thought that fast inbreeding (such as that generated by full-sib mating lines) removes only highly deleterious recessive alleles, while slow inbreeding can also remove mildly deleterious ones. However, a question remains regarding which proportion of the inbreeding load can be removed by purging under slow inbreeding in moderately large populations. We report results of two long-term slow inbreeding Drosophila experiments (125-234 generations), each using a large population and a number of derived lines with effective sizes about 1000 and 50, respectively. The inbreeding load was virtually exhausted after more than one hundred generations in large populations and between a few tens and over one hundred generations in the lines. This result is not expected from genetic drift alone, and is in agreement with the theoretical purging predictions. Computer simulations suggest that these results are consistent with a model of relatively few deleterious mutations of large homozygous effects and partially recessive gene action.

Adaptive Divergence of Meiotic Recombination Rate in Ecological Speciation

Theories predict that directional selection during adaptation to a novel habitat results in elevated meiotic recombination rate. Yet the lack of population-level recombination rate data leaves this hypothesis untested in natural populations. Here, we examine the population-level recombination rate variation in two incipient ecological species, the microcrustacean Daphnia pulex (an ephemeral-pond species) and Daphnia pulicaria (a permanent-lake species). The divergence of D. pulicaria from D. pulex involved habitat shifts from pond to lake habitats as well as strong local adaptation due to directional selection. Using a novel single-sperm genotyping approach, we estimated the male-specific recombination rate of two linkage groups in multiple populations of each species in common garden experiments and identified a significantly elevated recombination rate in D. pulicaria. Most importantly, population genetic analyses show that the divergence in recombination rate between these two species is most likely due to divergent selection in distinct ecological habitats rather than neutral evolution.

Why Do Some Vertebrates Have Microchromosomes?

With more than 70,000 living species, vertebrates have a huge impact on the field of biology and research, including karyotype evolution. One prominent aspect of many vertebrate karyotypes is the enigmatic occurrence of tiny and often cytogenetically indistinguishable microchromosomes, which possess distinctive features compared to macrochromosomes. Why certain vertebrate species carry these microchromosomes in some lineages while others do not, and how they evolve remain open questions. New studies have shown that microchromosomes exhibit certain unique characteristics of genome structure and organization, such as high gene densities, low heterochromatin levels, and high rates of recombination. Our review focuses on recent concepts to expand current knowledge on the dynamic nature of karyotype evolution in vertebrates, raising important questions regarding the evolutionary origins and ramifications of microchromosomes. We introduce the basic karyotypic features to clarify the size, shape, and morphology of macro- and microchromosomes and report their distribution across different lineages. Finally, we characterize the mechanisms of different evolutionary forces underlying the origin and evolution of microchromosomes.

Multicellularity and sex helped shape the Tree of Life

Across the Tree of Life, there are dramatic differences in species numbers among groups. However, the factors that explain the differences among the deepest branches have remained unknown. We tested whether multicellularity and sexual reproduction might explain these patterns, since the most species-rich groups share these traits. We found that groups with multicellularity and sexual reproduction have accelerated rates of species proliferation (diversification), and that multicellularity has a stronger effect than sexual reproduction. Patterns of species richness among clades are then strongly related to these differences in diversification rates. Taken together, these results help explain patterns of biodiversity among groups of organisms at the very broadest scales. They may also help explain the mysterious preponderance of sexual reproduction among species (the 'paradox of sex') by showing that organisms with sexual reproduction proliferate more rapidly.

TopDomain: Exhaustive Protein Domain Boundary Metaprediction Combining Multisource Information and Deep Learning

Protein domains are independent, functional, and stable structural units of proteins. Accurate protein domain boundary prediction plays an important role in understanding protein structure and evolution, as well as for protein structure prediction. Current domain boundary prediction methods differ in terms of boundary definition, methodology, and training databases resulting in disparate performance for different proteins. We developed TopDomain, an exhaustive metapredictor, that uses deep neural networks to combine multisource information from sequence- and homology-based features of over 50 primary predictors. For this purpose, we developed a new domain boundary data set termed the TopDomain data set, in which the true annotations are informed by SCOPe annotations, structural domain parsers, human inspection, and deep learning. We benchmark TopDomain against 2484 targets with 3354 boundaries from the TopDomain test set and achieve F1 scores of 78.4% and 73.8% for multidomain boundary prediction within ±20 residues and ±10 residues of the true boundary, respectively. When examined on targets from CASP11-13 competitions, TopDomain achieves F1 scores of 47.5% and 42.8% for multidomain proteins. TopDomain significantly outperforms 15 widely used, state-of-the-art ab initio and homology-based domain boundary predictors. Finally, we implemented TopDomain_TMC, which accurately predicts whether domain parsing is necessary for the target protein.

Negative linkage disequilibrium between amino acid changing variants reveals interference among deleterious mutations in the human genome

Evolutionary forces like Hill-Robertson interference and negative epistasis can lead to deleterious mutations being found on distinct haplotypes. However, the extent to which these forces depend on the selection and dominance coefficients of deleterious mutations and shape genome-wide patterns of linkage disequilibrium (LD) in natural populations with complex demographic histories has not been tested. In this study, we first used forward-in-time simulations to predict how negative selection impacts LD. Under models where deleterious mutations have additive effects on fitness, deleterious variants less than 10 kb apart tend to be carried on different haplotypes relative to pairs of synonymous SNPs. In contrast, for recessive mutations, there is no consistent ordering of how selection coefficients affect LD decay, due to the complex interplay of different evolutionary effects. We then examined empirical data of modern humans from the 1000 Genomes Project. LD between derived alleles at nonsynonymous SNPs is lower compared to pairs of derived synonymous variants, suggesting that nonsynonymous derived alleles tend to occur on different haplotypes more than synonymous variants. This result holds when controlling for potential confounding factors by matching SNPs for frequency in the sample (allele count), physical distance, magnitude of background selection, and genetic distance between pairs of variants. Lastly, we introduce a new statistic HR(j) which allows us to detect interference using unphased genotypes. Application of this approach to high-coverage human genome sequences confirms our finding that nonsynonymous derived alleles tend to be located on different haplotypes more often than are synonymous derived alleles. Our findings suggest that interference may play a pervasive role in shaping patterns of LD between deleterious variants in the human genome, and consequently influences genome-wide patterns of LD.

Movements of transposable elements contribute to the genomic plasticity and species diversification in an asexually reproducing nematode pest

Despite reproducing without sexual recombination, Meloidogyne incognita is an adaptive and versatile phytoparasitic nematode. This species displays a global distribution, can parasitize a large range of plants, and can overcome plant resistance in a few generations. The mechanisms underlying this adaptability remain poorly known. At the whole-genome level, only a few single nucleotide variations have been observed across different geographical isolates with distinct ranges of compatible hosts. Exploring other factors possibly involved in genomic plasticity is thus important. Transposable elements (TEs), by their repetitive nature and mobility, can passively and actively impact the genome dynamics. This is particularly expected in polyploid hybrid genomes such as the one of M. incognita. Here, we have annotated the TE content of M. incognita, analyzed the statistical properties of this TE landscape, and used whole-genome pool-seq data to estimate the mobility of these TEs across twelve geographical isolates, presenting variations in ranges of compatible host plants. DNA transposons are more abundant than retrotransposons, and the high similarity of TE copies to their consensus sequences suggests they have been at least recently active. We have identified loci in the genome where the frequencies of presence of a TE showed substantial variations across the different isolates. Overall, variations in TE frequencies across isolates followed their phylogenetic divergence, suggesting TEs participate in the species diversification. Compared with the M. incognita reference genome, we detected isolate and lineage-specific de novo insertion of some TEs, including within genic regions or in the upstream regulatory regions. We validated by PCR the insertion of some of these TEs inside genic regions, confirming TE movements have possible functional impacts. Overall, we show DNA transposons can drive genomic plasticity in M. incognita and their role in genome evolution of other parthenogenetic animal deserves further investigation.

Sex as information processing: Optimality and evolution

The long-term growth rate of populations in varying environments quantifies the evolutionary value of processing the information that biological individuals inherit from their ancestors and acquire from their environment. Previous models were limited to asexual reproduction with inherited information coming from a single parent with no recombination. We present a general extension to sexual reproduction and an analytical solution for a particular but important case, the infinitesimal model of quantitative genetics which assumes traits to be normally distributed. We study with this model the conditions under which sexual reproduction is advantageous and can evolve in the context of autocorrelated or directionally varying environments, mutational biases, spatial heterogeneities, and phenotypic plasticity. Our results generalize and unify previous analyses. We also examine the proposal made by Geodakyan that the presence of two phenotypically distinct sexes permits an optimal adaptation to varying environments. We verify that conditions exists where sexual dimorphism is adaptive but find that its evolutionary value does not generally compensate for the twofold cost of males.

Global distribution and evolutionary transitions of angiosperm sexual systems

Angiosperm sexual systems are fundamental to the evolution and distribution of plant diversity, yet spatiotemporal patterns in angiosperm sexual systems and their drivers remain poorly known. Using data on sexual systems and distributions of 68453 angiosperm species, we present the first global maps of sexual system frequencies and evaluate sexual system evolution during the Cenozoic. Frequencies of dioecy and monoecy increase with latitude, while hermaphrodites are more frequent in warm and arid regions. Transitions to dioecy from other states were higher than to hermaphroditism, but transitions away from dioecy increased since the Cenozoic, suggesting that dioecy is not an evolutionary end point. Transitions between hermaphroditism and dioecy increased, while transitions to monoecy decreased with paleo-temperature when paleo-temperature >0℃. Our study demonstrates the biogeography of angiosperm sexual systems from a macroecological perspective, and enhances our understanding of plant diversity patterns and their response to climate change.

Selection Maintains Protein Interactome Resilience in the Long-Term Evolution Experiment with Escherichia coli

Most cellular functions are carried out by a dynamic network of interacting proteins. An open question is whether the network properties of protein interactomes represent phenotypes under natural selection. One proposal is that protein interactomes have evolved to be resilient, such that they tend to maintain connectivity when proteins are removed from the network. This hypothesis predicts that interactome resilience should be maintained by natural selection during long-term experimental evolution. I tested this prediction by modeling the evolution of protein-protein interaction (PPI) networks in Lenski's long-term evolution experiment with Escherichia coli (LTEE). In this test, I removed proteins affected by nonsense, insertion, deletion, and transposon mutations in evolved LTEE strains, and measured the resilience of the resulting networks. I compared the rate of change of network resilience in each LTEE population to the rate of change of network resilience for corresponding randomized networks. The evolved PPI networks are significantly more resilient than networks in which random proteins have been deleted. Moreover, the evolved networks are generally more resilient than networks in which the random deletion of proteins was restricted to those disrupted in LTEE. These results suggest that evolution in the LTEE has favored PPI networks that are, on average, more resilient than expected from the genetic variation across the evolved strains. My findings therefore support the hypothesis that selection maintains protein interactome resilience over evolutionary time.

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.

Genome Number and Size Polymorphism in Zika Virus Infectious Units

The arthropod-borne Zika virus (ZIKV) infects humans and can cause severe neurological sequelae, particularly in fetuses infected in utero. How this virus has been able to spread across vast geological ranges and evolve in new host populations is not yet understood.

Zika virus (ZIKV; Flaviviridae, Flavivirus) is an arthropod-borne infection that can result in severe outcomes, particularly in fetuses infected in utero. It has been assumed that infection by ZIKV, as well as other viruses, is largely initiated by individual virus particles binding to and entering a cell. However, recent studies have demonstrated that multiple virus particles are frequently delivered to a cell simultaneously and that this collective particle delivery enhances infection. ZIKV is maintained in nature between Aedes aegypti mosquitos and vertebrate hosts, including humans. Human infection is initiated through the injection of a relatively small initial inoculum comprised of a genetically complex virus population. Since most mutations decrease virus fitness, collective particle transmission could benefit ZIKV and other arthropod-borne diseases by facilitating the maintenance of genetic complexity and adaptability during infection or through other mechanisms. Therefore, we utilized a barcoded ZIKV to quantify the number of virus genomes that initiate a plaque. We found that individual plaques contain a mean of 10 infecting viral genomes (range, 1 to 212). Few plaques contained more than two dominant genomes. To determine whether multigenome infectious units consist of collectively transmitting virions, infectious units of ZIKV were then separated mechanically by centrifugation, and heavier fractions were found to contain more genomes per plaque-forming unit, with larger diameters. Finally, larger/heavier infectious units reformed after removal. These data suggest that ZIKV populations consist of a variety of infectious unit sizes, likely mostly made up of aggregates, and only rarely begin with a single virus genome.

IMPORTANCE The arthropod-borne Zika virus (ZIKV) infects humans and can cause severe neurological sequelae, particularly in fetuses infected in utero. How this virus has been able to spread across vast geological ranges and evolve in new host populations is not yet understood. This research demonstrates a novel mechanism of ZIKV transmission through multigenome aggregates, providing insight into ZIKV evolution, immunologic evasion, and better future therapeutic design. This study shows that ZIKV plaques result from collections of genomes rather than individual genomes, increasing the potential for interactions between ZIKV genotypes.

Uniparental Genome Elimination in Australian Carp Gudgeons

Metazoans usually reproduce sexually, blending the unique identity of parental genomes for the next generation through functional crossing-over and recombination in meiosis. However, some metazoan lineages have evolved reproductive systems where offspring are either full (clonal) or partial (hemiclonal) genetic replicas. In the latter group, the process of uniparental genome elimination selectively eliminates either the maternal or paternal genome from germ cells, and only one parental genome is selected for transmission. Although fairly common in plants, hybridogenesis (i.e., clonal haploidization via chromosome elimination) remains a poorly understood process in animals. Here, we explore the proximal cytogenomic mechanisms of somatic and germ cell chromosomes in sexual and hybrid genotypes of Australian carp gudgeons (Hypseleotris) by tracing the fate of each set during mitosis (in somatic tissues) and meiosis (in gonads). Our comparative study of diploid hybrid and sexual individuals revealed visually functional gonads in male and female hybrid genotypes and generally high karyotype variability, although the number of chromosome arms remains constant. Our results delivered direct evidence for classic hybridogenesis as a reproductive mode in carp gudgeons. Two parental sets with integral structure in the hybrid soma (the F1 constitution) contrasted with uniparental chromosomal inheritance detected in gonads. The inheritance mode happens through premeiotic genome duplication of the parental genome to be transmitted, whereas the second parental genome is likely gradually eliminated already in juvenile individuals. The role of metacentric chromosomes in hybrid evolution is also discussed.

Competitive history shapes rapid evolution in a seasonal climate

Eco-evolutionary dynamics will play a critical role in determining species' fates as climatic conditions change. Unfortunately, we have little understanding of how rapid evolutionary responses to climate play out when species are embedded in the competitive communities that they inhabit in nature. We tested the effects of rapid evolution in response to interspecific competition on subsequent ecological and evolutionary trajectories in a seasonally changing climate using a field-based evolution experiment with Drosophila melanogaster Populations of D. melanogaster were either exposed, or not exposed, to interspecific competition with an invasive competitor, Zaprionus indianus, over the summer. We then quantified these populations' ecological trajectories (abundances) and evolutionary trajectories (heritable phenotypic change) when exposed to a cooling fall climate. We found that competition with Z. indianus in the summer affected the subsequent evolutionary trajectory of D. melanogaster populations in the fall, after all interspecific competition had ceased. Specifically, flies with a history of interspecific competition evolved under fall conditions to be larger and have lower cold fecundity and faster development than flies without a history of interspecific competition. Surprisingly, this divergent fall evolutionary trajectory occurred in the absence of any detectible effect of the summer competitive environment on phenotypic evolution over the summer or population dynamics in the fall. This study demonstrates that competitive interactions can leave a legacy that shapes evolutionary responses to climate even after competition has ceased, and more broadly, that evolution in response to one selective pressure can fundamentally alter evolution in response to subsequent agents of selection.

Phenotypic plasticity rather than genotype drives reproductive choices in Hydra populations

Facultative clonality is associated with complex life cycles where sexual and asexual forms can be exposed to contrasting selection pressures. Facultatively clonal animals often have distinct developmental capabilities that depend on reproductive mode (e.g., negligible senescence and exceptional regeneration ability in asexual individuals, which are lacking in sexual individuals). Understanding how these differences in life history strategies evolved is hampered by limited knowledge of the population structure underlying sexual and asexual forms in nature. Here we studied genetic differentiation of coexisting sexual and asexual Hydra oligactis polyps, a freshwater cnidarian where reproductive mode-dependent life history patterns are observed. We collected asexual and sexual polyps from 13 Central European water bodies and used restriction-site associated DNA sequencing to infer population structure. We detected high relatedness among populations and signs that hydras might spread with resting eggs through zoochory. We found no genetic structure with respect to mode of reproduction (asexual vs. sexual). On the other hand, clear evidence was found for phenotypic plasticity in mode of reproduction, as polyps inferred to be clones differed in reproductive mode. Moreover, we detected two cases of apparent sex change (males and females found within the same clonal lineages) in this species with supposedly stable sexes. Our study describes population genetic structure in Hydra for the first time, highlights the role of phenotypic plasticity in generating patterns of life history variation, and contributes to understanding the evolution of reproductive mode-dependent life history variation in coexisting asexual and sexual forms.

Episodic evolution of coadapted sets of amino acid sites in mitochondrial proteins

The rate of evolution differs between protein sites and changes with time. However, the link between these two phenomena remains poorly understood. Here, we design a phylogenetic approach for distinguishing pairs of amino acid sites that evolve concordantly, i.e., such that substitutions at one site trigger subsequent substitutions at the other; and also pairs of sites that evolve discordantly, so that substitutions at one site impede subsequent substitutions at the other. We distinguish groups of amino acid sites that undergo coordinated evolution and evolve discordantly from other such groups. In mitochondrion-encoded proteins of metazoans and fungi, we show that concordantly evolving sites are clustered in protein structures. By analysing the phylogenetic patterns of substitutions at concordantly and discordantly evolving site pairs, we find that concordant evolution has two distinct causes: epistatic interactions between amino acid substitutions and episodes of selection independently affecting substitutions at different sites. The rate of substitutions at concordantly evolving groups of protein sites changes in the course of evolution, indicating episodes of selection limited to some of the lineages. The phylogenetic positions of these changes are consistent between proteins, suggesting common selective forces underlying them.

The mode and rate of evolution of a protein site depends on the effect of its mutations on protein fitness. The fitness effect of a mutation itself can change in the course of evolution for at least two reasons. First, it can be modulated by substitutions occurring at other sites, a phenomenon called epistasis. Second, changes in selection can be non-epistatic, affecting sites independently of one another. Here, we analyse substitutions accumulated by the evolving lineages of the five proteins encoded by the mitochondrial genomes of thousands of species of metazoans and fungi. We show that substitutions at different amino acid sites occur in a coordinated fashion, and this coordination is caused both by epistasis and by episodes of selection affecting groups of sites. We partition each protein into several groups of concordantly evolving sites such that evolution of sites from different groups is discordant, and show that the proteins encoded by the mitochondrial genome consist of coevolving structural blocks. Some of these blocks have a clear functional specialization, e.g. are associated with interfaces between proteins composing respiratory complexes. Together, our results reveal a previously unrecognized complexity in the causes of variation in evolutionary rates between protein sites.

Island songbirds as windows into evolution in small populations

Due to their limited ranges and inherent isolation, island species have long been recognized as crucial systems for tackling a range of evolutionary questions, including in the early study of speciation.^1,2 Such species have been less studied in the understanding of the evolutionary forces driving DNA sequence evolution. Island species usually have lower census population sizes (N) than continental species and, supposedly, lower effective population sizes (Ne). Given that both the rates of change caused by genetic drift and by selection are dependent upon Ne, island species are theoretically expected to exhibit (1) lower genetic diversity, (2) less effective natural selection against slightly deleterious mutations,^3,4 and (3) a lower rate of adaptive evolution.^5-8 Here, we have used a large set of newly sequenced and published whole-genome sequences of Passerida species (14 insular and 11 continental) to test these predictions. We confirm that island species exhibit lower census size and Ne, supporting the hypothesis that the smaller area available on islands constrains the upper bound of Ne. In the insular species, we find lower nucleotide diversity in coding regions, higher ratios of non-synonymous to synonymous polymorphisms, and lower adaptive substitution rates. Our results provide robust evidence that the lower Ne experienced by island species has affected both the ability of natural selection to efficiently remove weakly deleterious mutations and also the adaptive potential of island species, therefore providing considerable empirical support for the nearly neutral theory. We discuss the implications for both evolutionary and conservation biology.

Present Status and Future Management Strategies for Sugarcane Yellow Leaf Virus: A Major Constraint to the Global Sugarcane Production

Sugarcane yellow leaf virus (SCYLV) is a distinct member of the Polerovirus genus of the Luteoviridae family. SCYLV is the major limitation to sugarcane production worldwide and presently occurring in most of the sugarcane growing countries. SCYLV having high genetic diversity within the species and presently ten genotypes are known to occur based on the complete genome sequence information. SCYLV is present in almost all the states of India where sugarcane is grown. Virion comprises of 180 coat protein units and are 24-29 nm in diameter. The genome of SCYLV is a monopartite and comprised of single-stranded (ss) positive-sense (+) linear RNA of about 6 kb in size. Virus genome consists of six open reading frames (ORFs) that are expressed by sub-genomic RNAs. The SCYLV is phloem-limited and transmitted by sugarcane aphid Melanaphis sacchari in a circulative and non-propagative manner. The other aphid species namely, Ceratovacuna lanigera, Rhopalosiphum rufiabdominalis, and R. maidis also been reported to transmit the virus. The virus is not transmitted mechanically, therefore, its transmission by M. sacchari has been studied in different countries. SCYLV has a limited natural host range and mainly infect sugarcane (Sachharum hybrid), grain sorghum (Sorghum bicolor), and Columbus grass (Sorghum almum). Recent insights in the protein-protein interactions of Polerovirus through protein interaction reporter (PIR) technology enable us to understand viral encoded proteins during virus replication, assembly, plant defence mechanism, short and long-distance travel of the virus. This review presents the recent understandings on virus biology, diagnosis, genetic diversity, virus-vector and host-virus interactions and conventional and next generation management approaches.

Idiosyncratic epistasis creates universals in mutational effects and evolutionary trajectories

Patterns of epistasis and shapes of fitness landscapes are of wide interest because of their bearings on a number of evolutionary theories. The common phenomena of slowing fitness increases during adaptations and diminishing returns from beneficial mutations are believed to reflect a concave fitness landscape and a preponderance of negative epistasis. Paradoxically, fitness decreases tend to decelerate and harm from deleterious mutations shrinks during the accumulation of random mutations-patterns thought to indicate a convex fitness landscape and a predominance of positive epistasis. Current theories cannot resolve this apparent contradiction. Here, we show that the phenotypic effect of a mutation varies substantially depending on the specific genetic background and that this idiosyncrasy in epistasis creates all of the above trends without requiring a biased distribution of epistasis. The idiosyncratic epistasis theory explains the universalities in mutational effects and evolutionary trajectories as emerging from randomness due to biological complexity.