Genome-wide evidence for efficient positive and purifying selection in Capsella grandiflora, a plant species with a large effective population size

Genetic Causes and Genomic Consequences of Breakdown of Distyly in Linum trigynum

Distyly is an iconic floral polymorphism governed by a supergene, which promotes efficient pollen transfer and outcrossing through reciprocal differences in the position of sexual organs in flowers, often coupled with heteromorphic self-incompatibility. Distyly has evolved convergently in multiple flowering plant lineages, but has also broken down repeatedly, often resulting in homostylous, self-compatible populations with elevated rates of self-fertilization. Here, we aimed to study the genetic causes and genomic consequences of the shift to homostyly in Linum trigynum, which is closely related to distylous Linum tenue. Building on a high-quality genome assembly, we show that L. trigynum harbors a genomic region homologous to the dominant haplotype of the distyly supergene conferring long stamens and short styles in L. tenue, suggesting that loss of distyly first occurred in a short-styled individual. In contrast to homostylous Primula and Fagopyrum, L. trigynum harbors no fixed loss-of-function mutations in coding sequences of S-linked distyly candidate genes. Instead, floral gene expression analyses and controlled crosses suggest that mutations downregulating the S-linked LtWDR-44 candidate gene for male self-incompatibility and/or anther height could underlie homostyly and self-compatibility in L. trigynum. Population genomic analyses of 224 whole-genome sequences further demonstrate that L. trigynum is highly self-fertilizing, exhibits significantly lower genetic diversity genome-wide, and is experiencing relaxed purifying selection and less frequent positive selection on nonsynonymous mutations relative to L. tenue. Our analyses shed light on the loss of distyly in L. trigynum, and advance our understanding of a common evolutionary transition in flowering plants.

Impact of pre-copulatory sexual cannibalism on genetic diversity and efficacy of selection

Factors that increase reproductive variance among individuals act to reduce effective population size (Ne), which accelerates the loss of genetic diversity and decreases the efficacy of purifying selection. These factors include sexual cannibalism, offspring investment and mating system. Pre-copulatory sexual cannibalism, where the female consumes the male prior to mating, exacerbates this effect. We performed comparative transcriptomics in two spider species, the cannibalistic Trechaleoides biocellata and the non-cannibalistic T. keyserlingi, to generate genomic evidence to support these predictions. First, we estimated heterozygosity and found that genetic diversity is relatively lower in the cannibalistic species. Second, we calculated dN/dS ratios as a measure of purifying selection; a higher dN/dS ratio indicated relaxed purifying selection in the cannibalistic species. These results are consistent with the hypothesis that sexual cannibalism impacts operational sex ratio and demographic processes, which interact with evolutionary forces to shape the genetic structure of populations. However, other factors such as the mating system and life-history traits contribute to shaping Ne. Comparative analyses across multiple contrasting species pairs would be required to disentangle these effects. Our study highlights that extreme behaviours such as pre-copulatory cannibalism may have profound eco-evolutionary effects.

Evolutionary radiation of the Eurasian Pinus species under pervasive gene flow

Evolutionary radiation, a pivotal aspect of macroevolution, offers valuable insights into evolutionary processes. The genus Pinus is the largest genus in conifers withc . $$ c. $$ 90% of the extant species emerged in the Miocene, which signifies a case of rapid diversification. Despite this remarkable history, our understanding of the mechanisms driving radiation within this expansive genus has remained limited. Using exome capture sequencing and a fossil-calibrated phylogeny, we investigated the divergence history, niche diversification, and introgression among 13 closely related Eurasian species spanning climate zones from the tropics to the boreal Arctic. We detected complex introgression among lineages in subsection Pinus at all stages of the phylogeny. Despite this widespread gene exchange, each species maintained its genetic identity and showed clear niche differentiation. Demographic analysis unveiled distinct population histories among these species, which further influenced the nucleotide diversity and efficacy of purifying and positive selection in each species. Our findings suggest that radiation in the Eurasian pines was likely fueled by interspecific recombination and further reinforced by their adaptation to distinct environments. Our study highlights the constraints and opportunities for evolutionary change, and the expectations of future adaptation in response to environmental changes in different lineages.

Forty Years of Inferential Methods in the Journals of the Society for Molecular Biology and Evolution

We are launching a series to celebrate the 40th anniversary of the first issue of Molecular Biology and Evolution. In 2024, we will publish virtual issues containing selected papers published in the Society for Molecular Biology and Evolution journals, Molecular Biology and Evolution and Genome Biology and Evolution. Each virtual issue will be accompanied by a perspective that highlights the historic and contemporary contributions of our journals to a specific topic in molecular evolution. This perspective, the first in the series, presents an account of the broad array of methods that have been published in the Society for Molecular Biology and Evolution journals, including methods to infer phylogenies, to test hypotheses in a phylogenetic framework, and to infer population genetic processes. We also mention many of the software implementations that make methods tractable for empiricists. In short, the Society for Molecular Biology and Evolution community has much to celebrate after four decades of publishing high-quality science including numerous important inferential methods.

Quantifying adaptive evolution and the effects of natural selection across the Norway spruce genome

Detecting natural selection is one of the major goals of evolutionary genomics. Here, we sequenced the whole genome of 25 Picea abies individuals and quantified the amount of selection across the genome. Using an estimate of the distribution of fitness effects, we showed that both negative selection and the rate of positively selected substitutions are very limited in coding regions. We found a positive correlation between the rate of adaptive substitutions and recombination rate and a negative correlation between the rate of adaptive substitutions and gene density, suggesting a widespread influence from Hill-Robertson interference on the efficiency of protein adaptation in P. abies. Finally, the distinct population statistics between genomic regions under either positive or balancing selection with that under neutral regions indicated the impact of natural selection on the genomic architecture of Norway spruce. Further gene ontology enrichment analysis for genes located in regions identified as undergoing either positive or long-term balancing selection also highlighted the specific molecular functions and biological processes that appear to be targets of selection in Norway spruce.

Inference of the distribution of fitness effects of mutations is affected by single nucleotide polymorphism filtering methods, sample size and population structure

The distribution of fitness effects (DFE) of new mutations has been of interest to evolutionary biologists since the concept of mutations arose. Modern population genomic data enable us to quantify the DFE empirically, but few studies have examined how data processing, sample size and cryptic population structure might affect the accuracy of DFE inference. We used simulated and empirical data (from Arabidopsis lyrata) to show the effects of missing data filtering, sample size, number of single nucleotide polymorphisms (SNPs) and population structure on the accuracy and variance of DFE estimates. Our analyses focus on three filtering methods-downsampling, imputation and subsampling-with sample sizes of 4-100 individuals. We show that (1) the choice of missing-data treatment directly affects the estimated DFE, with downsampling performing better than imputation and subsampling; (2) the estimated DFE is less reliable in small samples (<8 individuals), and becomes unpredictable with too few SNPs (<5000, the sum of 0- and 4-fold SNPs); and (3) population structure may skew the inferred DFE towards more strongly deleterious mutations. We suggest that future studies should consider downsampling for small data sets, and use samples larger than 4 (ideally larger than 8) individuals, with more than 5000 SNPs in order to improve the robustness of DFE inference and enable comparative analyses.

Influence of the mutation load on the genomic composition of hybrids between outcrossing and self-fertilizing species

Hybridization is a natural process whereby two diverging evolutionary lineages reproduce and create offspring of mixed ancestry. Differences in mating systems (e.g., self-fertilization and outcrossing) are expected to affect the direction and extent of hybridization and introgression in hybrid zones. Among other factors, selfers and outcrossers are expected to differ in their mutation loads. This has been studied both theoretically and empirically; however, conflicting predictions have been made on the effects mutation loads of parental species with different mating systems can have on the genomic composition of hybrids. Here, we develop a multi-locus, selective model to study how the different mutation load built up in selfers and outcrossers as a result of selective interference and homozygosity impact the long-term genetic composition of hybrid populations. Notably, our results emphasize that genes from the parental population with lesser mutation load get rapidly overrepresented in hybrid genomes, regardless of the hybrids own mating system. When recombination tends to be more important than mutation, outcrossers' genomes tend to be of higher quality and prevail. When recombination rates are low, however, selfers' genomes may reach higher quality than outcrossers' genomes and prevail in the hybrids. Taken together, these results provide concrete insights into one of the multiple factors influencing hybrid genome ancestry and introgression patterns in hybrid zones containing species with different mating systems.

Whole Transcriptome Sequencing Unveils the Genomic Determinants of Putative Somaclonal Variation in Mint (Mentha L.)

Mint (Mentha L., Lamiaceae) is a strongly scented herb of the family Lamiaceae that is grown mostly by clonal propagation, making it a valuable species for the study of somaclonal variation and its phenotypic consequences. The recent introduction of a few species of mint in South America, followed by a presumably rampant propagation, make this region particularly ideal for studying the extent of somaclonal genetic diversity. Hence, the objective of this work was to offer a preliminary characterization of somaclonal genetically coding diversity of the mint in the northern Andes in order to address the question of whether somaclonal variants may have emerged despite relatively recent introductions in a region where mint is not native. A total of 29 clonally propagated specimens, collected in mint export farms in the province of Antioquia, a major region for mint production in the northwest Andes of Colombia, were genotyped using RNA sequencing (RNA-Seq). SNP calling was carried out from the leaves’ transcriptome profiles of each plant by combining the GATK4 and TRINITY protocols, obtaining a total of 2033 loci across 912 transcripts with a minimum read depth of 20X and 4% of missing data. Unsupervised machine learning algorithms considered the K-means, AGNES and UPGMA approaches, all of which suggested three genetic clusters for M. spicata and a unique cluster for M. × piperita. The results indicate that at least two different origins of M. spicata reached the eastern region of the Antioquia province, clonally propagated in the locality ever since for local consumption and export. One of these ancestries had more population structure, possibly due to environmental or anthropological pressures that intervened in the fragmentation of this genetic group or to a higher somaclonal mutation rate. This work offers a first step into the study of the accumulation and transmission of presumably quasi-neutral somatic mutations at coding regions in an herbaceous clonally propagated scented species such as mint, likely favored by an expected population expansion after its Andean introduction. These ad hoc hypotheses warrant further study as part of future research.

Globally Relaxed Selection and Local Adaptation in Boechera stricta

The strength of selection varies among populations and across the genome, but the determinants of efficacy of selection remain unclear. In this study, we used whole-genome sequencing data from 467 Boechera stricta accessions to quantify the strength of selection and characterize the pattern of local adaptation. We found low genetic diversity on 0-fold degenerate sites and conserved non-coding sites, indicating functional constraints on these regions. The estimated distribution of fitness effects and the proportion of fixed substitutions suggest relaxed negative and positive selection in B. stricta. Among the four population groups, the NOR and WES groups have smaller effective population size (Ne), higher proportions of effectively neutral sites, and lower rates of adaptive evolution compared with UTA and COL groups, reflecting the effect of Ne on the efficacy of natural selection. We also found weaker selection on GC-biased sites compared with GC-conservative (unbiased) sites, suggested that GC-biased gene conversion has affected the strength of selection in B. stricta. We found mixed evidence for the role of the recombination rate on the efficacy of selection. The positive and negative selection was stronger in high-recombination regions compared with low-recombination regions in COL but not in other groups. By scanning the genome, we found different subsets of selected genes suggesting differential adaptation among B. stricta groups. These results show that differences in effective population size, nucleotide composition, and recombination rate are important determinants of the efficacy of selection. This study enriches our understanding of the roles of natural selection and local adaptation in shaping genomic variation.

OUP accepted manuscript

It is known that methods to estimate the rate of adaptive evolution, which are based on the McDonald–Kreitman test, can be biased by changes in effective population size. Here, we demonstrate theoretically that changes in population size can also generate an artifactual correlation between the rate of adaptive evolution and any factor that is correlated to the strength of selection acting against deleterious mutations. In this context, we have investigated whether several site-level factors influence the rate of adaptive evolution in the divergence of humans and chimpanzees, two species that have been inferred to have undergone population size contraction since they diverged. We find that the rate of adaptive evolution, relative to the rate of mutation, is higher for more exposed amino acids, lower for amino acid pairs that are more dissimilar in terms of their polarity, volume, and lower for amino acid pairs that are subject to stronger purifying selection, as measured by the ratio of the numbers of nonsynonymous to synonymous polymorphisms (p_N/p_S). All of these correlations are opposite to the artifactual correlations expected under contracting population size. We therefore conclude that these correlations are genuine.

Selection on Accessible Chromatin Regions in Capsella grandiflora

Accurate estimates of genome-wide rates and fitness effects of new mutations are essential for an improved understanding of molecular evolutionary processes. Although eukaryotic genomes generally contain a large noncoding fraction, functional noncoding regions and fitness effects of mutations in such regions are still incompletely characterized. A promising approach to characterize functional noncoding regions relies on identifying accessible chromatin regions (ACRs) tightly associated with regulatory DNA. Here, we applied this approach to identify and estimate selection on ACRs in Capsella grandiflora, a crucifer species ideal for population genomic quantification of selection due to its favorable population demography. We describe a population-wide ACR distribution based on ATAC-seq data for leaf samples of 16 individuals from a natural population. We use population genomic methods to estimate fitness effects and proportions of positively selected fixations (α) in ACRs and find that intergenic ACRs harbor a considerable fraction of weakly deleterious new mutations, as well as a significantly higher proportion of strongly deleterious mutations than comparable inaccessible intergenic regions. ACRs are enriched for expression quantitative trait loci (eQTL) and depleted of transposable element insertions, as expected if intergenic ACRs are under selection because they harbor regulatory regions. By integrating empirical identification of intergenic ACRs with analyses of eQTL and population genomic analyses of selection, we demonstrate that intergenic regulatory regions are an important source of nearly neutral mutations. These results improve our understanding of selection on noncoding regions and the role of nearly neutral mutations for evolutionary processes in outcrossing Brassicaceae species.

Mutation accumulation in inbreeding populations under evolution of the selfing rate

It is theoretically established that self-fertilization can facilitate mutation accumulation, thus increasing extinction risk. However, in previous studies, selfing rates are often set as fixed parameters, but in natural systems, evolution of selfing rates and deleterious mutations may mutually affect each other. I carried out simulations to investigate the dynamics of selfing rates and mutation accumulation, by allowing deleterious mutations to coevolve with alleles that modify the selfing rate (selfing modifiers). I found that selfing rates will often fluctuate over time, due to successive invasion of alleles that increase selfing and outcrossing. Since mutation fixation is mainly caused by Muller's ratchet, its rate is sensitive to the change of the selfing rate mutations will accumulate in a punctuated pattern. The dynamics are influenced by several factors, such as recombination and the selfing rate effects of selfing modifier loci. Also, such temporal variation produces variation of selfing rates and mutation accumulation rates between multiple conspecific populations, which can increase the average fitness across populations. As factors, such as the genomic mutation rate of deleterious mutations, can simultaneously influence the selfing rate and mutation fixation, effects of these factors on mutation accumulation rates can be complicated and non-monotonic.

The Genomic Selfing Syndrome Accompanies the Evolutionary Breakdown of Heterostyly

The evolutionary transition from outcrossing to selfing can have important genomic consequences. Decreased effective population size and the reduced efficacy of selection are predicted to play an important role in the molecular evolution of the genomes of selfing species. We investigated evidence for molecular signatures of the genomic selfing syndrome using 66 species of Primula including distylous (outcrossing) and derived homostylous (selfing) taxa. We complemented our comparative analysis with a microevolutionary study of P. chungensis, which is polymorphic for mating system and consists of both distylous and homostylous populations. We generated chloroplast and nuclear genomic data sets for distylous, homostylous, and distylous–homostylous species and identified patterns of nonsynonymous to synonymous divergence (d_N/d_S) and polymorphism (π_N/π_S) in species or lineages with contrasting mating systems. Our analysis of coding sequence divergence and polymorphism detected strongly reduced genetic diversity and heterozygosity, decreased efficacy of purifying selection, purging of large-effect deleterious mutations, and lower rates of adaptive evolution in samples from homostylous compared with distylous populations, consistent with theoretical expectations of the genomic selfing syndrome. Our results demonstrate that self-fertilization is a major driver of molecular evolutionary processes with genomic signatures of selfing evident in both old and relatively young homostylous populations.

The scales of coevolution: comparative phylogeography and genetic demography of a locally adapted venomous predator and its prey

Coevolutionary theory predicts that differences in the genetic demography of interacting species can influence patterns of local adaptation by affecting the potential of local populations to respond to selection. We conducted a comparative phylogeographical study of venomous rattlesnakes and their venom-resistant ground squirrel prey across California, and assessed how effective population size (Ne) estimates correspond with a previously documented pattern of rattlesnake local adaptation. Using RAD sequencing markers, we detected lineage relationships among both the rattlesnakes (Crotalus oreganus ssp.) and ground squirrels (Otospermophilus sp.) that are incongruent with previous phylogenetic hypotheses. Both rattlesnakes and squirrels share a deep divergence at the Sacramento–San Joaquin River Delta. At this broad phylogeographical scale, we found that the locally adapted rattlesnakes had higher Ne than squirrels. At the population scale, snakes also had larger Ne accompanied by larger values of several metrics of population genetic diversity. However, the specific magnitude of local adaptation of venom activity to ground squirrel venom resistance was not significantly correlated with local differences in Ne or other diversity statistics between predator and prey populations, suggesting that other factors in the geographic mosaic of coevolution contribute to the specific local-scale outcomes of this interaction. These results suggest an evolutionary mechanism that may explain some (but clearly not all) of rattlesnake local adaptation in this coevolutionary interaction – larger population sizes raise the adaptive potential of rattlesnakes compared to ground squirrels.

Glucosinolate profile and Myrosinase gene expression are modulated upon Plasmodiophora brassicae infection in cabbage

Clubroot is a devastating disease of Brassicaceae caused by the biotrophic protist Plasmodiophora brassicae. The progression of clubroot disease is modulated by the glucosinolate (GSL) profile of the host plant. GSL is hydrolysed by the enzyme myrosinase upon cell disruption and gives rise to metabolites like isothiocyanate, nitriles, thiocyanates, epithionitriles and oxazolidines. Some of these metabolites play important roles in the plant's defence mechanism. We identified 13 Myrosinase (Myro) and 28 Myrosinase-Binding Protein-like (MBP) genes from Brassica oleracea L. using a comparative genomics approach and characterised them through in silico analyses. We compared the expression patterns of these genes in a clubroot-susceptible line and a resistant line following inoculation with P. brassicae. Two BolMyro and 12 BolMBP genes were highly expressed in the susceptible line, whereas only one BolMyro and five BolMBP genes were highly expressed in the resistant line. Principal component analysis confirmed that specific GSL profiles and gene expression were modulated due to pathogen infection. Plants with higher levels of neoglucobrassicin, glucobrassicin and methooxyglucobrassicin produced disease symptoms and formed galls, whereas, plants with higher levels of sinigrin, hydroxyglucobrassicin and progoitrin produced less symptoms with almost no galls. Our results provide insights into the roles of Myro and MBP genes in GSL hydrolysis during P. brassicae infection, which will help for developing clubroot resistant cabbage lines.

Purifying Selection in Corvids Is Less Efficient on Islands

Theory predicts that deleterious mutations accumulate more readily in small populations. As a consequence, mutation load is expected to be elevated in species where life-history strategies and geographic or historical contingencies reduce the number of reproducing individuals. Yet, few studies have empirically tested this prediction using genome-wide data in a comparative framework. We collected whole-genome sequencing data for 147 individuals across seven crow species (Corvus spp.). For each species, we estimated the distribution of fitness effects of deleterious mutations and compared it with proxies of the effective population size N_e. Island species with comparatively smaller geographic range sizes had a significantly increased mutation load. These results support the view that small populations have an elevated risk of mutational meltdown, which may contribute to the higher extinction rates observed in island species.

Biased Gene Conversion Constrains Adaptation in Arabidopsis thaliana

Reduction of fitness due to deleterious mutations imposes a limit to adaptive evolution. By characterizing features that influence this genetic load we may better understand constraints on responses to both natural and human-mediated selection. Here, using whole-genome, transcriptome, and methylome data from >600 Arabidopsis thaliana individuals, we set out to identify important features influencing selective constraint. Our analyses reveal that multiple factors underlie the accumulation of maladaptive mutations, including gene expression level, gene network connectivity, and gene-body methylation. We then focus on a feature with major effect, nucleotide composition. The ancestral vs. derived status of segregating alleles suggests that GC-biased gene conversion, a recombination-associated process that increases the frequency of G and C nucleotides regardless of their fitness effects, shapes sequence patterns in A. thaliana Through estimation of mutational effects, we present evidence that biased gene conversion hinders the purging of deleterious mutations and contributes to a genome-wide signal of decreased efficacy of selection. By comparing these results to two outcrossing relatives, Arabidopsis lyrata and Capsella grandiflora, we find that protein evolution in A. thaliana is as strongly affected by biased gene conversion as in the outcrossing species. Last, we perform simulations to show that natural levels of outcrossing in A. thaliana are sufficient to facilitate biased gene conversion despite increased homozygosity due to selfing. Together, our results show that even predominantly selfing taxa are susceptible to biased gene conversion, suggesting that it may constitute an important constraint to adaptation among plant species.

The Evolutionary Forces Shaping Cis- and Trans-Regulation of Gene Expression within a Population of Outcrossing Plants

Understanding the persistence of genetic variation within populations has long been a goal of evolutionary biology. One promising route toward achieving this goal is using population genetic approaches to describe how selection acts on the loci associated with trait variation. Gene expression provides a model trait for addressing the challenge of the maintenance of variation because it can be measured genome-wide without information about how gene expression affects traits. Previous work has shown that loci affecting the expression of nearby genes (local or cis-eQTLs) are under negative selection, but we lack a clear understanding of the selective forces acting on variants that affect the expression of genes in trans. Here, we identify loci that affect gene expression in trans using genomic and transcriptomic data from one population of the obligately outcrossing plant, Capsella grandiflora. The allele frequencies of trans-eQTLs are consistent with stronger negative selection acting on trans-eQTLs than cis-eQTLs, and stronger negative selection acting on trans-eQTLs associated with the expression of multiple genes. However, despite this general pattern, we still observe the presence of a trans-eQTL at intermediate frequency that affects the expression of a large number of genes in the same coexpression module. Overall, our work highlights the different selective pressures shaping variation in cis- and trans-regulation.

Transposable Elements Are Important Contributors to Standing Variation in Gene Expression in Capsella Grandiflora

Transposable elements (TEs) make up a significant portion of eukaryotic genomes and are important drivers of genome evolution. However, the extent to which TEs affect gene expression variation on a genome-wide scale in comparison with other types of variants is still unclear. We characterized TE insertion polymorphisms and their association with gene expression in 124 whole-genome sequences from a single population of Capsella grandiflora, and contrasted this with the effects of single nucleotide polymorphisms (SNPs). Population frequency of insertions was negatively correlated with distance to genes, as well as density of conserved noncoding elements, suggesting that the negative effects of TEs on gene regulation are important in limiting their abundance. Rare TE variants strongly influence gene expression variation, predominantly through downregulation. In contrast, rare SNPs contribute equally to up- and down-regulation, but have a weaker individual effect than TEs. An expression quantitative trait loci (eQTL) analysis shows that a greater proportion of common TEs are eQTLs as opposed to common SNPs, and a third of the genes with TE eQTLs do not have SNP eQTLs. In contrast with rare TE insertions, common insertions are more likely to increase expression, consistent with recent models of cis-regulatory evolution favoring enhancer alleles. Taken together, these results imply that TEs are a significant contributor to gene expression variation and are individually more likely than rare SNPs to cause extreme changes in gene expression.

The Impact of Protein Architecture on Adaptive Evolution

Adaptive mutations play an important role in molecular evolution. However, the frequency and nature of these mutations at the intramolecular level are poorly understood. To address this, we analyzed the impact of protein architecture on the rate of adaptive substitutions, aiming to understand how protein biophysics influences fitness and adaptation. Using Drosophila melanogaster and Arabidopsis thaliana population genomics data, we fitted models of distribution of fitness effects and estimated the rate of adaptive amino-acid substitutions both at the protein and amino-acid residue level. We performed a comprehensive analysis covering genome, gene, and protein structure, by exploring a multitude of factors with a plausible impact on the rate of adaptive evolution, such as intron number, protein length, secondary structure, relative solvent accessibility, intrinsic protein disorder, chaperone affinity, gene expression, protein function, and protein-protein interactions. We found that the relative solvent accessibility is a major determinant of adaptive evolution, with most adaptive mutations occurring at the surface of proteins. Moreover, we observe that the rate of adaptive substitutions differs between protein functional classes, with genes encoding for protein biosynthesis and degradation signaling exhibiting the fastest rates of protein adaptation. Overall, our results suggest that adaptive evolution in proteins is mainly driven by intermolecular interactions, with host-pathogen coevolution likely playing a major role.

Population Genomics of Transitions to Selfing in Brassicaceae Model Systems

Many plants harbor complex mechanisms that promote outcrossing and efficient pollen transfer. These include floral adaptations as well as genetic mechanisms, such as molecular self-incompatibility (SI) systems. The maintenance of such systems over long evolutionary timescales suggests that outcrossing is favorable over a broad range of conditions. Conversely, SI has repeatedly been lost, often in association with transitions to self-fertilization (selfing). This transition is favored when the short-term advantages of selfing outweigh the costs, primarily inbreeding depression. The transition to selfing is expected to have major effects on population genetic variation and adaptive potential, as well as on genome evolution. In the Brassicaceae, many studies on the population genetic, gene regulatory, and genomic effects of selfing have centered on the model plant Arabidopsis thaliana and the crucifer genus Capsella. The accumulation of population genomics datasets have allowed detailed investigation of where, when and how the transition to selfing occurred. Future studies will take advantage of the development of population genetics theory on the impact of selfing, especially regarding positive selection. Furthermore, investigation of systems including recent transitions to selfing, mixed mating populations and/or multiple independent replicates of the same transition will facilitate dissecting the effects of mating system variation from processes driven by demography.

Variation of the adaptive substitution rate between species and within genomes

The importance of adaptive mutations in molecular evolution is extensively debated. Recent developments in population genomics allow inferring rates of adaptive mutations by fitting a distribution of fitness effects to the observed patterns of polymorphism and divergence at sites under selection and sites assumed to evolve neutrally. Here, we summarize the current state-of-the-art of these methods and review the factors that affect the molecular rate of adaptation. Several studies have reported extensive cross-species variation in the proportion of adaptive amino-acid substitutions (α) and predicted that species with larger effective population sizes undergo less genetic drift and higher rates of adaptation. Disentangling the rates of positive and negative selection, however, revealed that mutations with deleterious effects are the main driver of this population size effect and that adaptive substitution rates vary comparatively little across species. Conversely, rates of adaptive substitution have been documented to vary substantially within genomes. On a genome-wide scale, gene density, recombination and mutation rate were observed to play a role in shaping molecular rates of adaptation, as predicted under models of linked selection. At the gene level, it has been reported that the gene functional category and the macromolecular structure substantially impact the rate of adaptive mutations. Here, we deliver a comprehensive review of methods used to infer the molecular adaptive rate, the potential drivers of adaptive evolution and how positive selection shapes molecular evolution within genes, across genes within species and between species.

Population genomics of the facultatively asexual duckweed Spirodela polyrhiza

Clonal propagation allows some plant species to achieve massive population sizes quickly but also reduces the evolutionary independence of different sites in the genome. We examine genome-wide genetic diversity in Spirodela polyrhiza, a duckweed that reproduces primarily asexually. We find that this geographically widespread and numerically abundant species has very low levels of genetic diversity. Diversity at nonsynonymous sites relative to synonymous sites is high, suggesting that purifying selection is weak. A potential explanation for this observation is that a very low frequency of sex renders selection ineffective. However, there is a pronounced decay in linkage disequilibrium over 40 kb, suggesting that though sex may be rare at the individual level it is not too infrequent at the population level. In addition, neutral diversity is affected by the physical proximity of selected sites, which would be unexpected if sex was exceedingly rare at the population level. The amount of genetic mixing as assessed by the decay in linkage disequilibrium is not dissimilar from selfing species such as Arabidopsis thaliana, yet selection appears to be much less effective in duckweed. We discuss alternative explanations for the signature of weak purifying selection.

Impact of demography on linked selection in two outcrossing Brassicaceae species

Genetic diversity is shaped by mutation, genetic drift, gene flow, recombination, and selection. The dynamics and interactions of these forces shape genetic diversity across different parts of the genome, between populations and species. Here, we have studied the effects of linked selection on nucleotide diversity in outcrossing populations of two Brassicaceae species, Arabidopsis lyrata and Capsella grandiflora, with contrasting demographic history. In agreement with previous estimates, we found evidence for a modest population size expansion thousands of generations ago, as well as efficient purifying selection in C. grandiflora. In contrast, the A. lyrata population exhibited evidence for very recent strong population size decline and weaker efficacy of purifying selection. Using multiple regression analyses with recombination rate and other genomic covariates as explanatory variables, we can explain 47% of the variance in neutral diversity in the C. grandiflora population, while in the A. lyrata population, only 11% of the variance was explained by the model. Recombination rate had a significant positive effect on neutral diversity in both species, suggesting that selection at linked sites has an effect on patterns of neutral variation. In line with this finding, we also found reduced neutral diversity in the vicinity of genes in the C. grandiflora population. However, in A. lyrata no such reduction in diversity was evident, a finding that is consistent with expectations of the impact of a recent bottleneck on patterns of neutral diversity near genes. This study thus empirically demonstrates how differences in demographic history modulate the impact of selection at linked sites in natural populations.

Expansion history and environmental suitability shape effective population size in a plant invasion

The margins of an expanding range are predicted to be challenging environments for adaptation. Marginal populations should often experience low effective population sizes (N_e ) where genetic drift is high due to demographic expansion and/or census population size is low due to unfavourable environmental conditions. Nevertheless, invasive species demonstrate increasing evidence of rapid evolution and potential adaptation to novel environments encountered during colonization, calling into question whether significant reductions in N_e are realized during range expansions in nature. Here we report one of the first empirical tests of the joint effects of expansion dynamics and environment on effective population size variation during invasive range expansion. We estimate contemporary values of N_e using rates of linkage disequilibrium among genome-wide markers within introduced populations of the highly invasive plant Centaurea solstitialis (yellow starthistle) in North America (California, USA), and within native Eurasian populations. As predicted, we find that N_e within the invaded range is positively correlated with both expansion history (time since founding) and habitat quality (abiotic climate). History and climate had independent additive effects with similar effect sizes, indicating an important role for both factors in this invasion. These results support theoretical expectations for the population genetics of range expansion, though whether these processes can ultimately arrest the spread of an invasive species remains an unanswered question.

Deleterious Mutation Burden and Its Association with Complex Traits in Sorghum (Sorghum bicolor)

Sorghum (Sorghum bicolor (L.) Moench) is a major staple food cereal for millions of people worldwide. Valluru et al. identify putative deleterious mutations among ∼5.5M segregating variants of 229 diverse sorghum...

Sorghum (Sorghum bicolor L.) is a major food cereal for millions of people worldwide. The sorghum genome, like other species, accumulates deleterious mutations, likely impacting its fitness. The lack of recombination, drift, and the coupling with favorable loci impede the removal of deleterious mutations from the genome by selection. To study how deleterious variants impact phenotypes, we identified putative deleterious mutations among ∼5.5 M segregating variants of 229 diverse biomass sorghum lines. We provide the whole-genome estimate of the deleterious burden in sorghum, showing that ∼33% of nonsynonymous substitutions are putatively deleterious. The pattern of mutation burden varies appreciably among racial groups. Across racial groups, the mutation burden correlated negatively with biomass, plant height, specific leaf area (SLA), and tissue starch content (TSC), suggesting that deleterious burden decreases trait fitness. Putatively deleterious variants explain roughly one-half of the genetic variance. However, there is only moderate improvement in total heritable variance explained for biomass (7.6%) and plant height (average of 3.1% across all stages). There is no advantage in total heritable variance for SLA and TSC. The contribution of putatively deleterious variants to phenotypic diversity therefore appears to be dependent on the genetic architecture of traits. Overall, these results suggest that incorporating putatively deleterious variants into genomic models slightly improves prediction accuracy because of extensive linkage. Knowledge of deleterious variants could be leveraged for sorghum breeding through either genome editing and/or conventional breeding that focuses on the selection of progeny with fewer deleterious alleles.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

The Role of Small RNA-Based Epigenetic Silencing for Purifying Selection on Transposable Elements in Capsella grandiflora

To avoid negative effects of transposable element (TE) proliferation, plants epigenetically silence TEs using a number of mechanisms, including RNA-directed DNA methylation. These epigenetic modifications can extend outside the boundaries of TE insertions and lead to silencing of nearby genes, resulting in a trade-off between TE silencing and interference with nearby gene regulation. Therefore, purifying selection is expected to remove silenced TE insertions near genes more efficiently and prevent their accumulation within a population. To explore how effects of TE silencing on gene regulation shapes purifying selection on TEs, we analyzed whole genome sequencing data from 166 individuals of a large population of the outcrossing species Capsella grandiflora. We found that most TEs are rare, and in chromosome arms, silenced TEs are exposed to stronger purifying selection than those that are not silenced by 24-nucleotide small RNAs, especially with increasing proximity to genes. An age-of-allele test of neutrality on a subset of TEs supports our inference of purifying selection on silenced TEs, suggesting that our results are robust to varying transposition rates. Our results provide new insights into the processes affecting the accumulation of TEs in an outcrossing species and support the view that epigenetic silencing of TEs results in a trade-off between preventing TE proliferation and interference with nearby gene regulation. We also suggest that in the centromeric and pericentromeric regions, the negative aspects of epigenetic TE silencing are missing.

Patterns of polymorphism and selection in the subgenomes of the allopolyploid Arabidopsis kamchatica

Genome duplication is widespread in wild and crop plants. However, little is known about genome-wide selection in polyploids due to the complexity of duplicated genomes. In polyploids, the patterns of purifying selection and adaptive substitutions may be affected by masking owing to duplicated genes or homeologs as well as effective population size. Here, we resequence 25 accessions of the allotetraploid Arabidopsis kamchatica, which is derived from the diploid species A. halleri and A. lyrata. We observe a reduction in purifying selection compared with the parental species. Interestingly, proportions of adaptive non-synonymous substitutions are significantly positive in contrast to most plant species. A recurrent pattern observed in both frequency and divergence–diversity neutrality tests is that the genome-wide distributions of both subgenomes are similar, but the correlation between homeologous pairs is low. This may increase the opportunity of different evolutionary trajectories such as in the HMA4 gene involved in heavy metal hyperaccumulation.

Despite the prevalence of genome duplication in plants, little is known about the evolutionary patterns of entire subgenomes. Here the authors resequence allopolyploid Arabidopsis kamchatica genome to estimate diversity, linkage disequilibrium and strengths of both positive and purifying selection.

Competitive ability of Capsella species with different mating systems and ploidy levels

BACKGROUND AND AIMS

Capsella is a model genus for studying the transition from outcrossing to selfing, with or without change in ploidy levels. The genomic consequences and changes in reproductive traits (selfing syndrome) associated with these shifts have been studied in depth. However, potential ecological divergence among species of the genus has not been determined. Among ecological traits, competitive ability could be relevant for selfing evolution, as selfing has been shown to be statistically associated with reduced competitiveness in a recent meta-analysis.

METHODS

We assessed the effect of competition on three Capsella species differing in their mating system and ploidy level. We used an experimental design where fitness related traits were measured in focal individuals with and without competitors.

KEY RESULTS

The diploid selfer (C. rubella) was most sensitive to competition, whereas the tetraploid selfer (C. bursa-pastoris) performed the best, with the diploid outcrosser (C. grandiflora) being intermediate.

CONCLUSIONS

These results add to the detailed characterization of Capsella species and highlight the possible roles of ecological context and ploidy in the evolutionary trajectories of selfing species.

High rate of adaptive evolution in two widespread European pines

Comparing related organisms with differing ecological requirements and evolutionary histories can shed light on the mechanisms and drivers underlying genetic adaptation. Here, by examining a common set of hundreds of loci, we compare patterns of nucleotide diversity and molecular adaptation of two European conifers (Scots pine and maritime pine) living in contrasted environments and characterized by distinct population genetic structure (low and clinal in Scots pine, high and ecotypic in maritime pine) and demographic histories. We found higher nucleotide diversity in Scots pine than in maritime pine, whereas rates of new adaptive substitutions (ω_a ), as estimated from the distribution of fitness effects, were similar across species and among the highest found in plants. Sample size and population genetic structure did not appear to have resulted in significant bias in estimates of ω_a . Moreover, population contraction-expansion dynamics for each species did not affect differentially the rate of adaptive substitution in these two pines. Several methodological and biological factors may underlie the unusually high rate of adaptive evolution of Scots pine and maritime pine. By providing two new case studies with contrasting evolutionary histories, we contribute to disentangling the multiple factors potentially affecting adaptive evolution in natural plant populations.

Natural Selection Drives Rapid Functional Evolution of Young Drosophila Duplicate Genes

Gene duplication is thought to play a major role in phenotypic evolution. Yet the forces involved in the functional divergence of young duplicate genes remain unclear. Here, we use population-genetic inference to elucidate the role of natural selection in the functional evolution of young duplicate genes in Drosophila melanogaster. We find that negative selection acts on young duplicates with ancestral functions, and positive selection on those with novel functions, suggesting that natural selection may determine whether and how young duplicate genes are retained. Moreover, evidence of natural selection is strongest in protein-coding regions and 3' UTRs of young duplicates, indicating that selection may primarily target encoded proteins and regulatory sequences specific to 3' UTRs. Further analysis reveals that natural selection acts immediately after duplication and weakens over time, possibly explaining the observed bias toward the acquisition of new functions by young, rather than old, duplicate gene copies. Last, we find an enrichment of testis-related functions in young duplicates that underwent recent positive selection, but not in young duplicates that did not undergo recent positive selection, or in old duplicates that either did or did not undergo recent positive selection. Thus, our findings reveal that natural selection is a key player in the functional evolution of young duplicate genes, acts rapidly and in a region-specific manner, and may underlie the origin of novel testis-specific phenotypes in Drosophila.

Transcriptome Signatures of Selection, Drift, Introgression, and Gene Duplication in the Evolution of an Extremophile Endemic Plant

Plants on serpentine soils provide extreme examples of adaptation to environment, and thus offer excellent models for the study of evolution at the molecular and genomic level. Serpentine outcrops are derived from ultramafic rock and have extremely low levels of essential plant nutrients (e.g., N, P, K, and Ca), as well as toxic levels of heavy metals (e.g., Ni, Cr, and Co) and low moisture availability. These outcrops provide habitat to a number of endemic plant species, including the annual mustard Caulanthus amplexicaulis var. barbarae (Cab) (Brassicaceae). Its sister taxon, C. amplexicaulis var. amplexicaulis (Caa), is intolerant to serpentine soils. Here, we assembled and annotated comprehensive reference transcriptomes of both Caa and Cab for use in protein coding sequence comparisons. A set of 29,443 reciprocal best Blast hit (RBH) orthologs between Caa and Cab was compared with identify coding sequence variants, revealing a high genome-wide dN/dS ratio between the two taxa (mean = 0.346). We show that elevated dN/dS likely results from the composite effects of genetic drift, positive selection, and the relaxation of negative selection. Further, analysis of paralogs within each taxon revealed the signature of a period of elevated gene duplication (∼10 Ma) that is shared with other species of the tribe Thelypodieae, and may have played a role in the striking morphological and ecological diversity of this tribe. In addition, distribution of the synonymous substitution rate, dS, is strongly bimodal, indicating a history of reticulate evolution that may have contributed to serpentine adaptation.

Contrasting Rates of Molecular Evolution and Patterns of Selection among Gymnosperms and Flowering Plants

The majority of variation in rates of molecular evolution among seed plants remains both unexplored and unexplained. Although some attention has been given to flowering plants, reports of molecular evolutionary rates for their sister plant clade (gymnosperms) are scarce, and to our knowledge differences in molecular evolution among seed plant clades have never been tested in a phylogenetic framework. Angiosperms and gymnosperms differ in a number of features, of which contrasting reproductive biology, life spans, and population sizes are the most prominent. The highly conserved morphology of gymnosperms evidenced by similarity of extant species to fossil records and the high levels of macrosynteny at the genomic level have led scientists to believe that gymnosperms are slow-evolving plants, although some studies have offered contradictory results. Here, we used 31,968 nucleotide sites obtained from orthologous genes across a wide taxonomic sampling that includes representatives of most conifers, cycads, ginkgo, and many angiosperms with a sequenced genome. Our results suggest that angiosperms and gymnosperms differ considerably in their rates of molecular evolution per unit time, with gymnosperm rates being, on average, seven times lower than angiosperm species. Longer generation times and larger genome sizes are some of the factors explaining the slow rates of molecular evolution found in gymnosperms. In contrast to their slow rates of molecular evolution, gymnosperms possess higher substitution rate ratios than angiosperm taxa. Finally, our study suggests stronger and more efficient purifying and diversifying selection in gymnosperm than in angiosperm species, probably in relation to larger effective population sizes.

Inference of Distribution of Fitness Effects and Proportion of Adaptive Substitutions from Polymorphism Data

The distribution of fitness effects (DFE) encompasses the fraction of deleterious, neutral, and beneficial mutations. It conditions the evolutionary trajectory of populations, as well as the rate of adaptive molecular evolution (α). Inferring DFE and α from patterns of polymorphism, as given through the site frequency spectrum (SFS) and divergence data, has been a longstanding goal of evolutionary genetics. A widespread assumption shared by previous inference methods is that beneficial mutations only contribute negligibly to the polymorphism data. Hence, a DFE comprising only deleterious mutations tends to be estimated from SFS data, and α is then predicted by contrasting the SFS with divergence data from an outgroup. We develop a hierarchical probabilistic framework that extends previous methods to infer DFE and α from polymorphism data alone. We use extensive simulations to examine the performance of our method. While an outgroup is still needed to obtain an unfolded SFS, we show that both a DFE, comprising both deleterious and beneficial mutations, and α can be inferred without using divergence data. We also show that not accounting for the contribution of beneficial mutations to polymorphism data leads to substantially biased estimates of the DFE and α. We compare our framework with one of the most widely used inference methods available and apply it on a recently published chimpanzee exome data set.

The Evolutionary Basis of Translational Accuracy in Plants

Mistranslation errors compromise fitness by wasting resources on nonfunctional proteins. In order to reduce the cost of mistranslations, natural selection chooses the most accurately translated codons at sites that are particularly important for protein structure and function. We investigated the determinants underlying selection for translational accuracy in several species of plants belonging to three clades: Brassicaceae, Fabidae, and Poaceae. Although signatures of translational selection were found in genes from a wide range of species, the underlying factors varied in nature and intensity. Indeed, the degree of synonymous codon bias at evolutionarily conserved sites varied among plant clades while remaining uniform within each clade. This is unlikely to solely reflect the diversity of tRNA pools because there is little correlation between synonymous codon bias and tRNA abundance, so other factors must affect codon choice and translational accuracy in plant genes. Accordingly, synonymous codon choice at a given site was affected not only by the selection pressure at that site, but also its participation in protein domains or mRNA secondary structures. Although these effects were detected in all the species we analyzed, their impact on translation accuracy was distinct in evolutionarily distant plant clades. The domain effect was found to enhance translational accuracy in dicot and monocot genes with a high GC content, but to oppose the selection of more accurate codons in monocot genes with a low GC content.

Adaptive sequence evolution is driven by biotic stress in a pair of orchid species (Dactylorhiza) with distinct ecological optima

The orchid family is the largest in the angiosperms, but little is known about the molecular basis of the significant variation they exhibit. We investigate here the transcriptomic divergence between two European terrestrial orchids, Dactylorhiza incarnata and Dactylorhiza fuchsii, and integrate these results in the context of their distinct ecologies that we also document. Clear signals of lineage-specific adaptive evolution of protein-coding sequences are identified, notably targeting elements of biotic defence, including both physical and chemical adaptations in the context of divergent pools of pathogens and herbivores. In turn, a substantial regulatory divergence between the two species appears linked to adaptation/acclimation to abiotic conditions. Several of the pathways affected by differential expression are also targeted by deviating post-transcriptional regulation via sRNAs. Finally, D. incarnata appears to suffer from insufficient sRNA control over the activity of RNA-dependent DNA polymerase, resulting in increased activity of class I transposable elements and, over time, in larger genome size than that of D. fuchsii. The extensive molecular divergence between the two species suggests significant genomic and transcriptomic shock in their hybrids and offers insights into the difficulty of coexistence at the homoploid level. Altogether, biological response to selection, accumulated during the history of these orchids, appears governed by their microenvironmental context, in which biotic and abiotic pressures act synergistically to shape transcriptome structure, expression and regulation.

Aging asexual lineages and the evolutionary maintenance of sex

Finite populations of asexual and highly selfing species suffer from a reduced efficacy of selection. Such populations are thought to decline in fitness over time due to accumulating slightly deleterious mutations or failing to adapt to changing conditions. These within-population processes that lead nonrecombining species to extinction may help maintain sex and outcrossing through species level selection. Although inefficient selection is proposed to elevate extinction rates over time, previous models of species selection for sex assumed constant diversification rates. For sex to persist, classic models require that asexual species diversify at rates lower than sexual species; the validity of this requirement is questionable, both conceptually and empirically. We extend past models by allowing asexual lineages to decline in diversification rates as they age, that is nonrecombining lineages "senesce" in diversification rates. At equilibrium, senescing diversification rates maintain sex even when asexual lineages, at young ages, diversify faster than their sexual progenitors. In such cases, the age distribution of asexual lineages contains a peak at intermediate values rather than showing the exponential decline predicted by the classic model. Coexistence requires only that the average rate of diversification in asexuals be lower than that of sexuals.

Genomic analysis reveals major determinants of cis-regulatory variation in Capsella grandiflora

Understanding the causes of cis-regulatory variation is a long-standing aim in evolutionary biology. Although cis-regulatory variation has long been considered important for adaptation, we still have a limited understanding of the selective importance and genomic determinants of standing cis-regulatory variation. To address these questions, we studied the prevalence, genomic determinants, and selective forces shaping cis-regulatory variation in the outcrossing plant Capsella grandiflora We first identified a set of 1,010 genes with common cis-regulatory variation using analyses of allele-specific expression (ASE). Population genomic analyses of whole-genome sequences from 32 individuals showed that genes with common cis-regulatory variation (i) are under weaker purifying selection and (ii) undergo less frequent positive selection than other genes. We further identified genomic determinants of cis-regulatory variation. Gene body methylation (gbM) was a major factor constraining cis-regulatory variation, whereas presence of nearby transposable elements (TEs) and tissue specificity of expression increased the odds of ASE. Our results suggest that most common cis-regulatory variation in C. grandiflora is under weak purifying selection, and that gene-specific functional constraints are more important for the maintenance of cis-regulatory variation than genome-scale variation in the intensity of selection. Our results agree with previous findings that suggest TE silencing affects nearby gene expression, and provide evidence for a link between gbM and cis-regulatory constraint, possibly reflecting greater dosage sensitivity of body-methylated genes. Given the extensive conservation of gbM in flowering plants, this suggests that gbM could be an important predictor of cis-regulatory variation in a wide range of plant species.

Patterns of Substitution Rate Variation at Many Nuclear Loci in Two Species Trios in the Brassicaceae Partitioned with ANOVA

There are marked variations among loci and among lineages in rates of nucleotide substitution. The generation time hypothesis (GTH) is a neutral explanation for substitution rate heterogeneity that has genomewide application, predicting that species with shorter generation times accumulate DNA sequence substitutions faster than species with longer generation times do since faster genome replication provides more opportunities for mutations to occur and reach fixation by genetic drift. Relatively few studies have rigorously evaluated the GTH in plants, and there are numerous alternative hypotheses for plant substitution rate variation. One major challenge has been finding pairs of closely related plant species with contrasting generation times and appropriate outgroup taxa that all also have DNA sequence data for numerous loci. To test for causes of rate variation, we obtained sequence data for 256 genes for Arabidopsis thaliana, normally reproducing every year, and the biennial Arabidopsis lyrata with three closely related outgroup taxa (Brassica rapa, Capsella grandiflora, and Neslia paniculata) as well as the biennial Brassica oleracea and the annual B. rapa lineage with the outgroup N. paniculata. A sign test indicated that more loci than expected by chance have faster rates of substitution on the branch leading to the annual than to the perennial for one three-species trio but not another. Tajima's 1D and 2D tests, and a likelihood ratio test that incorporated saturation correction, rejected rate homogeneity for up to 26 genes (up to 14 genes when correcting for multiple tests), consistently showing faster rates for the annual lineage in the Arabidopsis species trio. ANOVA showed significant rate heterogeneity between the Arabidopsis and Brassica species trios (about 6 % of rate variation) and among loci (about 26-32 % of rate variation). The lineage-by-locus interaction which would be caused by locus- and lineage-specific natural selection explained about 13 % of substitution rate variation in one ANOVA model using substitution rates from genes partitioned into odd and even codons but was not a significant effect without partitioned genes. Annual/perennial lineage and species trio by annual/perennial lineage each explained about 1 % of substitution rate variation.

Limits to Adaptation in Partially Selfing Species

In outcrossing populations, "Haldane's sieve" states that recessive beneficial alleles are less likely to fix than dominant ones, because they are less exposed to selection when rare. In contrast, selfing organisms are not subject to Haldane's sieve and are more likely to fix recessive types than outcrossers, as selfing rapidly creates homozygotes, increasing overall selection acting on mutations. However, longer homozygous tracts in selfers also reduce the ability of recombination to create new genotypes. It is unclear how these two effects influence overall adaptation rates in partially selfing organisms. Here, we calculate the fixation probability of beneficial alleles if there is an existing selective sweep in the population. We consider both the potential loss of the second beneficial mutation if it has a weaker advantage than the first one, and the possible replacement of the initial allele if the second mutant is fitter. Overall, loss of weaker adaptive alleles during a first selective sweep has a larger impact on preventing fixation of both mutations in highly selfing organisms. Furthermore, the presence of linked mutations has two opposing effects on Haldane's sieve. First, recessive mutants are disproportionally likely to be lost in outcrossers, so it is likelier that dominant mutations will fix. Second, with elevated rates of adaptive mutation, selective interference annuls the advantage in selfing organisms of not suffering from Haldane's sieve; outcrossing organisms are more able to fix weak beneficial mutations of any dominance value. Overall, weakened recombination effects can greatly limit adaptation in selfing organisms.

Genomic legacies of the progenitors and the evolutionary consequences of allopolyploidy

The formation of an allopolyploid species involves the merger of genomes with separate evolutionary histories and thereby different genomic legacies. Contrary to expectations from theory, genes from one are often lost preferentially in allopolyploids - there is biased fractionation. Here, we provide an overview of two ways in which the genomic legacies of the progenitors may impact the fate of duplicated genes in allopolyploids. Specifically, we discuss the role of homeolog expression biases in setting the stage for biased fractionation, and the evidence for transposable element silencing as a possible mechanism for homeolog expression biases. Finally, we highlight how differences between the progenitors with respect to accumulation of deleterious variation may affect trajectories of duplicate gene evolution in allopolyploids.

Phylogenetic analysis suggests that sociality is associated with reduced effectiveness of selection

The evolution of sociality in spiders is associated with female bias, reproductive skew and an inbreeding mating system, factors that cause a reduction in effective population size and increase effects of genetic drift. These factors act to decrease the effectiveness of selection, thereby increasing the fixation probability of deleterious mutations. Comparative studies of closely related species with contrasting social traits and mating systems provide the opportunity to test consequences of low effective population size on the effectiveness of selection empirically. We used phylogenetic analyses of three inbred social spider species and seven outcrossing subsocial species of the genus Stegodyphus, and compared dN/dS ratios and codon usage bias between social Inbreeding and subsocial outcrossing mating systems to assess the effectiveness of selection. The overall results do not differ significantly between the social inbreeding and outcrossing species, but suggest a tendency for lower codon usage bias and higher dN/dS ratios in the social inbreeding species compared with their outcrossing congeners. The differences in dN/dS ratio and codon usage bias between social and subsocial species are modest but consistent with theoretical expectations of reduced effectiveness of selection in species with relatively low effective population size. The modest differences are consistent with relatively recent evolution of social mating systems. Additionally, the short terminal branches and lack of speciation of the social lineages, together with low genetic diversity lend support for the transient state of permanent sociality in spiders.

Natural Selection and Recombination Rate Variation Shape Nucleotide Polymorphism Across the Genomes of Three Related Populus Species

A central aim of evolutionary genomics is to identify the relative roles that various evolutionary forces have played in generating and shaping genetic variation within and among species. Here we use whole-genome resequencing data to characterize and compare genome-wide patterns of nucleotide polymorphism, site frequency spectrum, and population-scaled recombination rates in three species of Populus: Populus tremula, P. tremuloides, and P. trichocarpa. We find that P. tremuloides has the highest level of genome-wide variation, skewed allele frequencies, and population-scaled recombination rates, whereas P. trichocarpa harbors the lowest. Our findings highlight multiple lines of evidence suggesting that natural selection, due to both purifying and positive selection, has widely shaped patterns of nucleotide polymorphism at linked neutral sites in all three species. Differences in effective population sizes and rates of recombination largely explain the disparate magnitudes and signatures of linked selection that we observe among species. The present work provides the first phylogenetic comparative study on a genome-wide scale in forest trees. This information will also improve our ability to understand how various evolutionary forces have interacted to influence genome evolution among related species.

Association mapping reveals the role of purifying selection in the maintenance of genomic variation in gene expression

The evolutionary forces that maintain genetic variation in quantitative traits within populations remain poorly understood. One hypothesis suggests that variation is under purifying selection, resulting in an excess of low-frequency variants and a negative correlation between minor allele frequency and selection coefficients. Here, we test these predictions using the genetic loci associated with total expression variation (eQTLs) and allele-specific expression variation (aseQTLs) mapped within a single population of the plant Capsella grandiflora. In addition to finding eQTLs and aseQTLs for a large fraction of genes, we show that alleles at these loci are rarer than expected and exhibit a negative correlation between phenotypic effect size and frequency. Overall, our results show that the distribution of frequencies and effect sizes of the loci responsible for local expression variation within a single outcrossing population are consistent with the effects of purifying selection.

Adaptive Evolution Is Substantially Impeded by Hill-Robertson Interference in Drosophila

Hill-Robertson interference (HRi) is expected to reduce the efficiency of natural selection when two or more linked selected sites do not segregate freely, but no attempt has been done so far to quantify the overall impact of HRi on the rate of adaptive evolution for any given genome. In this work, we estimate how much HRi impedes the rate of adaptive evolution in the coding genome of Drosophila melanogaster. We compiled a data set of 6,141 autosomal protein-coding genes from Drosophila, from which polymorphism levels in D. melanogaster and divergence out to D. yakuba were estimated. The rate of adaptive evolution was calculated using a derivative of the McDonald-Kreitman test that controls for slightly deleterious mutations. We find that the rate of adaptive amino acid substitution at a given position of the genome is positively correlated to both the rate of recombination and the mutation rate, and negatively correlated to the gene density of the region. These correlations are robust to controlling for each other, for synonymous codon bias and for gene functions related to immune response and testes. We show that HRi diminishes the rate of adaptive evolution by approximately 27%. Interestingly, genes with low mutation rates embedded in gene poor regions lose approximately 17% of their adaptive substitutions whereas genes with high mutation rates embedded in gene rich regions lose approximately 60%. We conclude that HRi hampers the rate of adaptive evolution in Drosophila and that the variation in recombination, mutation, and gene density along the genome affects the HRi effect.

Cis-Regulatory Changes Associated with a Recent Mating System Shift and Floral Adaptation in Capsella

The selfing syndrome constitutes a suite of floral and reproductive trait changes that have evolved repeatedly across many evolutionary lineages in response to the shift to selfing. Convergent evolution of the selfing syndrome suggests that these changes are adaptive, yet our understanding of the detailed molecular genetic basis of the selfing syndrome remains limited. Here, we investigate the role of cis-regulatory changes during the recent evolution of the selfing syndrome in Capsella rubella, which split from the outcrosser Capsella grandiflora less than 200 ka. We assess allele-specific expression (ASE) in leaves and flower buds at a total of 18,452 genes in three interspecific F1 C. grandiflora x C. rubella hybrids. Using a hierarchical Bayesian approach that accounts for technical variation using genomic reads, we find evidence for extensive cis-regulatory changes. On average, 44% of the assayed genes show evidence of ASE; however, only 6% show strong allelic expression biases. Flower buds, but not leaves, show an enrichment of cis-regulatory changes in genomic regions responsible for floral and reproductive trait divergence between C. rubella and C. grandiflora. We further detected an excess of heterozygous transposable element (TE) insertions near genes with ASE, and TE insertions targeted by uniquely mapping 24-nt small RNAs were associated with reduced expression of nearby genes. Our results suggest that cis-regulatory changes have been important during the recent adaptive floral evolution in Capsella and that differences in TE dynamics between selfing and outcrossing species could be important for rapid regulatory divergence in association with mating system shifts.