The Accumulation of Deleterious Mutations as a Consequence of Domestication and Improvement in Sunflowers and Other Compositae Crops

A foundational large language model for edible plant genomes

Significant progress has been made in the field of plant genomics, as demonstrated by the increased use of high-throughput methodologies that enable the characterization of multiple genome-wide molecular phenotypes. These findings have provided valuable insights into plant traits and their underlying genetic mechanisms, particularly in model plant species. Nonetheless, effectively leveraging them to make accurate predictions represents a critical step in crop genomic improvement. We present AgroNT, a foundational large language model trained on genomes from 48 plant species with a predominant focus on crop species. We show that AgroNT can obtain state-of-the-art predictions for regulatory annotations, promoter/terminator strength, tissue-specific gene expression, and prioritize functional variants. We conduct a large-scale in silico saturation mutagenesis analysis on cassava to evaluate the regulatory impact of over 10 million mutations and provide their predicted effects as a resource for variant characterization. Finally, we propose the use of the diverse datasets compiled here as the Plants Genomic Benchmark (PGB), providing a comprehensive benchmark for deep learning-based methods in plant genomic research. The pre-trained AgroNT model is publicly available on HuggingFace at https://huggingface.co/InstaDeepAI/agro-nucleotide-transformer-1b  for future research purposes.

Genomic Reconstruction of the Successful Establishment of a Feralized Bovine Population on the Subantarctic Island of Amsterdam

The feral cattle of the subantarctic island of Amsterdam provide an outstanding case study of a large mammalian population that was established by a handful of founders and thrived within a few generations in a seemingly inhospitable environment. Here, we investigated the genetic history and composition of this population using genotyping and sequencing data. Our inference showed an intense but brief founding bottleneck around the late 19th century and revealed contributions from European taurine and Indian Ocean Zebu in the founder ancestry. Comparative analysis of whole-genome sequences further revealed a moderate reduction in genetic diversity despite high levels of inbreeding. The brief and intense bottleneck was associated with high levels of drift, a flattening of the site frequency spectrum and a slight relaxation of purifying selection on mildly deleterious variants. Unlike some populations that have experienced prolonged reductions in effective population size, we did not observe any significant purging of highly deleterious variants. Interestingly, the population's success in the harsh environment can be attributed to preadaptation from their European taurine ancestry, suggesting no strong bioclimatic challenge, and also contradicting evidence for insular dwarfism. Genome scan for footprints of selection uncovered a majority of candidate genes related to nervous system function, likely reflecting rapid feralization driven by behavioral changes and complex social restructuring. The Amsterdam Island cattle offers valuable insights into rapid population establishment, feralization, and genetic adaptation in challenging environments. It also sheds light on the unique genetic legacies of feral populations, raising ethical questions according to conservation efforts.

When birds of a feather flock together: Severe genomic erosion and the implications for genetic rescue in an endangered island passerine

The Seychelles magpie-robin's (SMR) five island populations exhibit some of the lowest recorded levels of genetic diversity among endangered birds, and high levels of inbreeding. These populations collapsed during the 20th century, and the species was listed as Critically Endangered in the IUCN Red List in 1994. An assisted translocation-for-recovery program initiated in the 1990s increased the number of mature individuals, resulting in its downlisting to Endangered in 2005. Here, we explore the temporal genomic erosion of the SMR based on a dataset of 201 re-sequenced whole genomes that span the past ~150 years. Our sample set includes individuals that predate the bottleneck by up to 100 years, as well as individuals from contemporary populations established during the species recovery program. Despite the SMR's recent demographic recovery, our data reveal a marked increase in both the genetic load and realized load in the extant populations when compared to the historical samples. Conservation management may have reduced the intensity of selection by increasing juvenile survival and relaxing intraspecific competition between individuals, resulting in the accumulation of loss-of-function mutations (i.e. severely deleterious variants) in the rapidly recovering population. In addition, we found a 3-fold decrease in genetic diversity between temporal samples. While the low genetic diversity in modern populations may limit the species' adaptability to future environmental changes, future conservation efforts (including IUCN assessments) may also need to assess the threats posed by their high genetic load. Our computer simulations highlight the value of translocations for genetic rescue and show how this could halt genomic erosion in threatened species such as the SMR.

Patterns of the Predicted Mutation Burden in 19,778 Domesticated Barley Accessions Conserved Ex Situ

Long-term conservation of more than 7 million plant germplasm accessions in 1750 genebanks worldwide is a challenging mission. The extent of deleterious mutations present in conserved germplasm and the genetic risk associated with accumulative mutations are largely unknown. This study took advantage of published barley genomic data to predict sample-wise mutation burdens for 19,778 domesticated barley (Hordeum vulgare L.) accessions conserved ex situ. It was found that the conserved germplasm harbored 407 deleterious mutations and 337 (or 82%) identified deleterious alleles were present in 20 (or 0.1%) or fewer barley accessions. Analysis of the predicted mutation burdens revealed significant differences in mutation burden for several groups of barley germplasm (landrace > cultivar (or higher burden estimate in landrace than in cultivar); winter barley > spring barley; six-rowed barley > two-rowed barley; and 1000-accession core collection > non-core germplasm). Significant differences in burden estimate were also found among seven major geographical regions. The sample-wise predicted mutation burdens were positively correlated with the estimates of sample average pairwise genetic difference. These findings are significant for barley germplasm management and utilization and for a better understanding of the genetic risk in conserved plant germplasm.

Deleterious phenotypes in wild Arabidopsis arenosa populations are common and linked to runs of homozygosity

In this study, we aimed to systematically assess the frequency at which potentially deleterious phenotypes appear in natural populations of the outcrossing model plant Arabidopsis arenosa, and to establish their underlying genetics. For this purpose, we collected seeds from wild A. arenosa populations and screened over 2,500 plants for unusual phenotypes in the greenhouse. We repeatedly found plants with obvious phenotypic defects, such as small stature and necrotic or chlorotic leaves, among first-generation progeny of wild A. arenosa plants. Such abnormal plants were present in about 10% of maternal sibships, with multiple plants with similar phenotypes in each of these sibships, pointing to a genetic basis of the observed defects. A combination of transcriptome profiling, linkage mapping and genome-wide runs of homozygosity patterns using a newly assembled reference genome indicated a range of underlying genetic architectures associated with phenotypic abnormalities. This included evidence for homozygosity of certain genomic regions, consistent with alleles that are identical by descent being responsible for these defects. Our observations suggest that deleterious alleles with different genetic architectures are segregating at appreciable frequencies in wild A. arenosa populations.

Porous borders at the wild-crop interface promote weed adaptation in Southeast Asia

High reproductive compatibility between crops and their wild relatives can provide benefits for crop breeding but also poses risks for agricultural weed evolution. Weedy rice is a feral relative of rice that infests paddies and causes severe crop losses worldwide. In regions of tropical Asia where the wild progenitor of rice occurs, weedy rice could be influenced by hybridization with the wild species. Genomic analysis of this phenomenon has been very limited. Here we use whole genome sequence analyses of 217 wild, weedy and cultivated rice samples to show that wild rice hybridization has contributed substantially to the evolution of Southeast Asian weedy rice, with some strains acquiring weed-adaptive traits through introgression from the wild progenitor. Our study highlights how adaptive introgression from wild species can contribute to agricultural weed evolution, and it provides a case study of parallel evolution of weediness in independently-evolved strains of a weedy crop relative.

Deleterious and Adaptive Mutations in Plant Germplasm Conserved Ex Situ

Conserving more than 7 million plant germplasm accessions in 1,750 genebanks worldwide raises the hope of securing the food supply for humanity for future generations. However, there is a genetic cost for such long-term germplasm conservation, which has been largely unaccounted for before. We investigated the extent and variation of deleterious and adaptive mutations in 490 individual plants representing barley, wheat, oat, soybean, maize, rapa, and sunflower collections in a seed genebank using RNA-Seq technology. These collections were found to have a range of deleterious mutations detected from 125 (maize) to 83,695 (oat) with a mean of 13,537 and of the averaged sample-wise mutation burden per deleterious locus from 0.069 to 0.357 with a mean of 0.200. Soybean and sunflower collections showed that accessions acquired earlier had increased mutation burdens. The germplasm with more years of storage in several collections carried more deleterious and fewer adaptive mutations. The samples with more cycles of germplasm regeneration revealed fewer deleterious and more adaptive mutations. These findings are significant for understanding mutational dynamics and genetic cost in conserved germplasm and have implications for long-term germplasm management and conservation.

Temporal collections to study invasion biology

Biological invasions represent an extraordinary opportunity to study evolution. This is because accidental or deliberate species introductions have taken place for centuries across large geographical scales, frequently prompting rapid evolutionary transitions in invasive populations. Until recently, however, the utility of invasions as evolutionary experiments has been hampered by limited information on the makeup of populations that were part of earlier invasion stages. Now, developments in ancient and historical DNA technologies, as well as the quickening pace of digitization for millions of specimens that are housed in herbaria and museums globally, promise to help overcome this obstacle. In this review, we first introduce the types of temporal data that can be used to study invasions, highlighting the timescale captured by each approach and their respective limitations. We then discuss how ancient and historical specimens as well as data available from prior invasion studies can be used to answer questions on mechanisms of (mal)adaptation, rates of evolution, or community-level changes during invasions. By bridging the gap between contemporary and historical invasive populations, temporal data can help us connect pattern to process in invasion science. These data will become increasingly important if invasions are to achieve their full potential as experiments of evolution in nature.

Crop adaptation to climate change: An evolutionary perspective

The disciplines of evolutionary biology and plant and animal breeding have been intertwined throughout their development, with responses to artificial selection yielding insights into the action of natural selection and evolutionary biology providing statistical and conceptual guidance for modern breeding. Here we offer an evolutionary perspective on a grand challenge of the 21st century: feeding humanity in the face of climate change. We first highlight promising strategies currently under way to adapt crops to current and future climate change. These include methods to match crop varieties with current and predicted environments and to optimize breeding goals, management practices, and crop microbiomes to enhance yield and sustainable production. We also describe the promise of crop wild relatives and recent technological innovations such as speed breeding, genomic selection, and genome editing for improving environmental resilience of existing crop varieties or for developing new crops. Next, we discuss how methods and theory from evolutionary biology can enhance these existing strategies and suggest novel approaches. We focus initially on methods for reconstructing the evolutionary history of crops and their pests and symbionts, because such historical information provides an overall framework for crop-improvement efforts. We then describe how evolutionary approaches can be used to detect and mitigate the accumulation of deleterious mutations in crop genomes, identify alleles and mutations that underlie adaptation (and maladaptation) to agricultural environments, mitigate evolutionary trade-offs, and improve critical proteins. Continuing feedback between the evolution and crop biology communities will ensure optimal design of strategies for adapting crops to climate change.

Genomic signatures of bottleneck and founder effects in dingoes

Dingoes arrived in Australia during the mid-Holocene and are the top-order terrestrial predator on the continent. Although dingoes subsequently spread across the continent, the initial founding population(s) could have been small. We investigated this hypothesis by sequencing the whole genomes of three dingoes and also obtaining the genome data from nine additional dingoes and 56 canines, including wolves, village dogs and breed dogs, and examined the signatures of bottlenecks and founder effects. We found that the nucleotide diversity of dingoes was low, 36% less than highly inbred breed dogs and 3.3 times lower than wolves. The number of runs of homozygosity (RoH) segments in dingoes was 1.6-4.7 times higher than in other canines. While examining deleterious mutational load, we observed that dingoes carried elevated ratios of nonsynonymous-to-synonymous diversities, significantly higher numbers of homozygous deleterious Single Nucleotide Variants (SNVs), and increased numbers of loss of function SNVs, compared to breed dogs, village dogs, and wolves. Our findings can be explained by bottlenecks and founder effects during the establishment of dingoes in mainland Australia. These findings highlight the need for conservation-based management of dingoes and the need for wildlife managers to be cognisant of these findings when considering the use of lethal control measures across the landscape.

Unraveling Prevalence and Effects of Deleterious Mutations in Maize Elite Lines across Decades of Modern Breeding

Future breeding is likely to involve the detection and removal of deleterious alleles, which are mutations that negatively affect crop fitness. However, little is known about the prevalence of such mutations and their effects on phenotypic traits in the context of modern crop breeding. To address this, we examined the number and frequency of deleterious mutations in 350 elite maize inbred lines developed over the past few decades in China and the United States. Our findings reveal an accumulation of weakly deleterious mutations and a decrease in strongly deleterious mutations, indicating the dominant effects of genetic drift and purifying selection for the two types of mutations, respectively. We also discovered that slightly deleterious mutations, when at lower frequencies, were more likely to be heterozygous in the developed hybrids. This is consistent with complementation as a potential explanation for heterosis. Subsequently, we found that deleterious mutations accounted for more of the variation in phenotypic traits than nondeleterious mutations with matched minor allele frequencies, especially for traits related to leaf angle and flowering time. Moreover, we detected fewer deleterious mutations in the promoter and gene body regions of differentially expressed genes across breeding eras than in nondifferentially expressed genes. Overall, our results provide a comprehensive assessment of the prevalence and impact of deleterious mutations in modern maize breeding and establish a useful baseline for future maize improvement efforts.

Evolutionary dynamics and adaptive benefits of deleterious mutations in crop gene pools

Mutations with deleterious consequences in nature may be conditionally deleterious in crop plants. That is, while some genetic variants may reduce fitness under wild conditions and be subject to purifying selection, they can be under positive selection in domesticates. Such deleterious alleles can be plant breeding targets, particularly for complex traits. The difficulty of distinguishing favorable from unfavorable variants reduces the power of selection, while favorable trait variation and heterosis may be attributable to deleterious alleles. Here, we review the roles of deleterious mutations in crop breeding and discuss how they can be used as a new avenue for crop improvement with emerging genomic tools, including HapMaps and pangenome analysis, aiding the identification, removal, or exploitation of deleterious mutations.

Deleterious mutation load in the admixed mice population

Deleterious mutation loads are known to correlate negatively with effective population size (Ne). Due to this reason, previous studies observed a higher proportion of harmful mutations in small populations than that in large populations. However, the mutational load in an admixed population that derived from introgression between individuals from two populations with vastly different Ne is not known. We investigated this using the whole genome data from two subspecies of the mouse (Mus musculus castaneus and Mus musculus musculus) with significantly different Ne. We used the ratio of diversities at nonsynonymous and synonymous sites (dN/dS) to measure the harmful mutation load. Our results showed that this ratio observed for the admixed population was intermediate between those of the parental populations. The dN/dS ratio of the hybrid population was significantly higher than that of M. m. castaneus but lower than that of M. m. musculus. Our analysis revealed a significant positive correlation between the proportion of M. m. musculus ancestry in admixed individuals and their dN/dS ratio. This suggests that the admixed individuals with high proportions of M. m. musculus ancestry have large dN/dS ratios. We also used the proportion of deleterious nonsynonymous SNVs as a proxy for deleterious mutation load, which also produced similar results. The observed results were in concordance with those expected by theory. We also show a shift in the distribution of fitness effects of nonsynonymous SNVs in the admixed genomes compared to the parental populations. These findings suggest that the deleterious mutation load of the admixed population is determined by the proportion of the ancestries of the subspecies. Therefore, it is important to consider the status and the level of genetic admixture of the populations whilst estimating the mutation loads.

Genomic diversity and signals of selection processes in wild and farm-reared red-legged partridges (Alectoris rufa)

The genetic dynamics of wild populations with releases of farm-reared reinforcements are very complex. These releases can endanger wild populations through genetic swamping or by displacing them. We assessed the genomic differences between wild and farm-reared red-legged partridges (Alectoris rufa) and described differential selection signals between both populations. We sequenced the whole genome of 30 wild and 30 farm-reared partridges. Both partridges had similar nucleotide diversity (π). Farm-reared partridges had a more negative Tajima's D and more and longer regions of extended haplotype homozygosity than wild partridges. We observed higher inbreeding coefficients (F_IS and F_ROH) in wild partridges. Selective sweeps (Rsb) were enriched with genes that contribute to the reproductive, skin and feather colouring, and behavioural differences between wild and farm-reared partridges. The analysis of genomic diversity should inform future decisions for the preservation of wild populations.

Deleterious Variation in Natural Populations and Implications for Conservation Genetics

Deleterious mutations decrease reproductive fitness and are ubiquitous in genomes. Given that many organisms face ongoing threats of extinction, there is interest in elucidating the impact of deleterious variation on extinction risk and optimizing management strategies accounting for such mutations. Quantifying deleterious variation and understanding the effects of population history on deleterious variation are complex endeavors because we do not know the strength of selection acting on each mutation. Further, the effect of demographic history on deleterious mutations depends on the strength of selection against the mutation and the degree of dominance. Here we clarify how deleterious variation can be quantified and studied in natural populations. We then discuss how different demographic factors, such as small population size, nonequilibrium population size changes, inbreeding, and gene flow, affect deleterious variation. Lastly, we provide guidance on studying deleterious variation in nonmodel populations of conservation concern.

Crop germplasm: Current challenges, physiological-molecular perspective, and advance strategies towards development of climate-resilient crops

Germplasm is a long-term resource management mission and investment for civilization. An estimated ∼7.4 million accessions are held in 1750 plant germplasm centres around the world; yet, only 2% of these assets have been utilized as plant genetic resources (PGRs). According to recent studies, the current food yield trajectory will be insufficient to feed the world's population in 2050. Additionally, possible negative effects in terms of crop failure because of climate change are already being experienced across the world. Therefore, it is necessary to reconciliation of research advancement and innovation of practices for further exploration of the potential of crop germplasm especially for the complex traits associated with yield such as water- and nitrogen use efficiency. In this review, we tried to address current challenges, research gaps, physiological and molecular aspects of two broad spectrum complex traits such as water- and nitrogen-use efficiency, and advanced integrated strategies that could provide a platform for combined stress management for climate-smart crop development. Additionally, recent development in technologies that are directly related to germplasm characterization was highlighted for further molecular utilization towards the development of elite varieties.

Unintended Consequences of Plant Domestication

Human selection on wild populations mostly favored a common set of plant traits during domestication. This process of direct selection also altered other independent traits that were not directly perceived or desired during crop domestication and improvement. A deeper knowledge of the inadvertent and undesirable phenotypic effects and their underlying genetic causes can help design strategies to mitigate their effects and improve genetic gain in crop plants. We review different factors explaining the negative consequences of plant domestication at the phenotypic and genomic levels. We further describe the genetic causes of undesirable effects that originate from the selection of favorable alleles during plant domestication. In addition, we propose strategies that could be useful in attenuating such effects for crop improvement. With novel -omics and genome-editing tools, it is relatively approachable to understand and manipulate the genetic and biochemical mechanisms responsible for the undesirable phenotypes in domesticated plants.

Deleterious mutations and the rare allele burden on rice gene expression

Deleterious genetic variation is maintained in populations at low frequencies. Under a model of stabilizing selection, rare (and presumably deleterious) genetic variants are associated with increase or decrease in gene expression from some intermediate optimum. We investigate this phenomenon in a population of largely Oryza sativa ssp. indica rice landraces under normal unstressed wet and stressful drought field conditions. We include single nucleotide polymorphisms, insertion/deletion mutations, and structural variants in our analysis and find a stronger association between rare variants and gene expression outliers under the stress condition. We also show an association of the strength of this rare variant effect with linkage, gene expression levels, network connectivity, local recombination rate, and fitness consequence scores, consistent with the stabilizing selection model of gene expression.

Future of Bacterial Disease Management in Crop Production

Bacterial diseases are a constant threat to crop production globally. Current management strategies rely on an array of tactics, including improved cultural practices; application of bactericides, plant activators, and biocontrol agents; and use of resistant varieties when available. However, effective management remains a challenge, as the longevity of deployed tactics is threatened by constantly changing bacterial populations. Increased scrutiny of the impact of pesticides on human and environmental health underscores the need for alternative solutions that are durable, sustainable, accessible to farmers, and environmentally friendly. In this review, we discuss the strengths and shortcomings of existing practices and dissect recent advances that may shape the future of bacterial disease management. We conclude that disease resistance through genome modification may be the most effective arsenal against bacterial diseases. Nonetheless, more research is necessary for developing novel bacterial disease management tactics to meet the food demand of a growing global population.

Give and Take: Effects of Genetic Admixture on Mutation Load in Endangered Florida Panthers

Genetic admixture is a biological event inherent to genetic rescue programs aimed at the long-term conservation of endangered wildlife. Although the success of such programs can be measured by the increase in genetic diversity and fitness of subsequent admixed individuals, predictions supporting admixture costs to fitness due to the introduction of novel deleterious alleles are necessary. Here, we analyzed nonsynonymous variation from conserved genes to quantify and compare levels of mutation load (i.e., proportion of deleterious alleles and genotypes carrying these alleles) among endangered Florida panthers and non-endangered Texas pumas. Specifically, we used canonical (i.e., non-admixed) Florida panthers, Texas pumas, and F1 (canonical Florida x Texas) panthers dating from a genetic rescue program and Everglades National Park panthers with Central American ancestry resulting from an earlier admixture event. We found neither genetic drift nor selection significantly reduced overall proportions of deleterious alleles in the severely bottlenecked canonical Florida panthers. Nevertheless, the deleterious alleles identified were distributed into a disproportionately high number of homozygous genotypes due to close inbreeding in this group. Conversely, admixed Florida panthers (either with Texas or Central American ancestry) presented reduced levels of homozygous genotypes carrying deleterious alleles but increased levels of heterozygous genotypes carrying these variants relative to canonical Florida panthers. Although admixture is likely to alleviate the load of standing deleterious variation present in homozygous genotypes, our results suggest introduced novel deleterious alleles (temporarily present in heterozygous state) in genetically rescued populations could potentially be expressed in subsequent generations if their effective sizes remain small.

The genomic signature of wild‐to‐crop introgression during the domestication of scarlet runner bean ( Phaseolus coccineus L.)

The scarlet runner bean (Phaseolus coccineus) is one of the five domesticated Phaseolus species. It is cultivated in small‐scale agriculture in the highlands of Mesoamerica for its dry seeds and immature pods, and unlike the other domesticated beans, P. coccineus is an open‐pollinated legume. Contrasting with its close relative, the common bean, few studies focusing on its domestication history have been conducted. Demographic bottlenecks associated with domestication might reduce genetic diversity and facilitate the accumulation of deleterious mutations. Conversely, introgression from wild relatives could be a source of variation. Using Genotyping by Sequencing data (79,286 single‐nucleotide variants) from 237 cultivated and wild samples, we evaluated the demographic history of traditional varieties from different regions of Mexico and looked for evidence of introgression between sympatric wild and cultivated populations. Traditional varieties have high levels of diversity, even though there is evidence of a severe initial genetic bottleneck followed by a population expansion. Introgression from wild to domesticated populations was detected, which might contribute to the recovery of the genetic variation. Introgression has occurred at different times: constantly in the center of Mexico; recently in the North West; and anciently in the South. Several factors are acting together to increase and maintain genetic diversity in P. coccineus cultivars, such as demographic expansion and introgression. Wild relatives represent a valuable genetic resource and have played a key role in scarlet runner bean evolution via introgression into traditional varieties.

Rare and population-specific functional variation across pig lines

BACKGROUND

It is expected that functional, mainly missense and loss-of-function (LOF), and regulatory variants are responsible for most phenotypic differences between breeds and genetic lines of livestock species that have undergone diverse selection histories. However, there is still limited knowledge about the existing missense and LOF variation in commercial livestock populations, in particular regarding population-specific variation and how it can affect applications such as across-breed genomic prediction.

METHODS

We re-sequenced the whole genome of 7848 individuals from nine commercial pig lines (average sequencing coverage: 4.1×) and imputed whole-genome genotypes for 440,610 pedigree-related individuals. The called variants were categorized according to predicted functional annotation (from LOF to intergenic) and prevalence level (number of lines in which the variant segregated; from private to widespread). Variants in each category were examined in terms of their distribution along the genome, alternative allele frequency, per-site Wright's fixation index (F_ST), individual load, and association to production traits.

RESULTS

Of the 46 million called variants, 28% were private (called in only one line) and 21% were widespread (called in all nine lines). Genomic regions with a low recombination rate were enriched with private variants. Low-prevalence variants (called in one or a few lines only) were enriched for lower allele frequencies, lower F_ST, and putatively functional and regulatory roles (including LOF and deleterious missense variants). On average, individuals carried fewer private deleterious missense alleles than expected compared to alleles with other predicted consequences. Only a small subset of the low-prevalence variants had intermediate allele frequencies and explained small fractions of phenotypic variance (up to 3.2%) of production traits. The significant low-prevalence variants had higher per-site F_ST than the non-significant ones. These associated low-prevalence variants were tagged by other more widespread variants in high linkage disequilibrium, including intergenic variants.

CONCLUSIONS

Most low-prevalence variants have low minor allele frequencies and only a small subset of low-prevalence variants contributed detectable fractions of phenotypic variance of production traits. Accounting for low-prevalence variants is therefore unlikely to noticeably benefit across-breed analyses, such as the prediction of genomic breeding values in a population using reference populations of a different genetic background.

Mutation load in sunflower inversions is negatively correlated with inversion heterozygosity

Recombination is critical both for accelerating adaptation and purging deleterious mutations. Chromosomal inversions can act as recombination modifiers that suppress local recombination in heterozygotes and thus, under some conditions, are predicted to accumulate such mutations. In this study, we investigated patterns of recombination, transposable element abundance and coding sequence evolution across the genomes of 1,445 individuals from three sunflower species, as well as within nine inversions segregating within species. We also analyzed the effects of inversion genotypes on 87 phenotypic traits to test for overdominance. We found significant negative correlations of long terminal repeat retrotransposon abundance and deleterious mutations with recombination rates across the genome in all three species. However, we failed to detect an increase in these features in the inversions, except for a modest increase in the proportion of stop codon mutations in several very large or rare inversions. Consistent with this finding, there was little evidence of overdominance of inversions in phenotypes that may relate to fitness. On the other hand, significantly greater load was observed for inversions in populations polymorphic for a given inversion compared to populations monomorphic for one of the arrangements, suggesting that the local state of inversion polymorphism affects deleterious load. These seemingly contradictory results can be explained by the low frequency of inversion heterozygotes in wild sunflower populations, apparently due to divergent selection and associated geographic structure. Inversions contributing to local adaptation represent ideal recombination modifiers, acting to facilitate adaptive divergence with gene flow, while largely escaping the accumulation of deleterious mutations.

Comparative genomics and signatures of selection in North Atlantic eels

Comparative genomic approaches can identify putative private and shared signatures of selection. We performed a comparative genomic study of North Atlantic eels, European eel (Anguilla Anguilla) and American eel (A. rostrata). The two sister species are nearly undistinguishable at the phenotypic level and despite a wide non-overlapping continental distribution, they spawn in partial sympatry in the Sargasso Sea. Taking advantage of the newly assembled and annotated genome, we used genome wide RAD sequencing data of 359 individuals retrieved from Sequence Nucleotide Archive and state-of-the-art statistic tests to identify putative genomic signatures of selection in North Atlantic eels. First, using the F_ST and XP-EHH methods, we detected apparent islands of divergence on a total of 7 chromosomes, particularly on chromosomes 6 and 10. Gene ontology analyses suggested candidate genes mainly related to energy production, development and regulation, which could reflect strong selection on traits related to eel migration and larval duration time. Gene effect prediction using SNPeff showed a high number of SNPs in noncoding regions, pointing to a possible regulatory role. Second, using the iHS method we detected shared regions under selection on a total of 11 chromosomes. Several hypotheses might account for the detection of shared islands of selection in North Atlantic eels, including parallel evolution due to adaptation to similar environments and introgression. Future comparative genomic studies will be needed to further clarify the causes and consequences of introgression, including the directionality of these introgression events.

No disruption of rhizobial symbiosis during early stages of cowpea domestication

Modern agriculture intensely selects aboveground plant structures, while often neglecting belowground features, and evolutionary tradeoffs between these traits are predicted to disrupt host control over microbiota. Moreover, drift, inbreeding, and relaxed selection for symbiosis in crops might degrade plant mechanisms that support beneficial microbes. We studied the impact of domestication on the nitrogen fixing symbiosis between cowpea and root-nodulating Bradyrhizobium. We combined genome-wide analyses with a greenhouse inoculation study to investigate genomic diversity, heritability, and symbiosis trait variation among wild and early-domesticated cowpea genotypes. Cowpeas experienced modest decreases in genome-wide diversity during early domestication. Nonetheless, domesticated cowpeas responded efficiently to variation in symbiotic effectiveness, by forming more root nodules with nitrogen-fixing rhizobia and sanctioning non-fixing strains. Domesticated populations invested a larger proportion of host tissues into root nodules than wild cowpeas. Unlike soybean and wheat, cowpea showed no compelling evidence for degradation of symbiosis during domestication. Domesticated cowpeas experienced a less severe bottleneck than these crops and the low nutrient conditions in Africa where cowpea landraces were developed likely favored plant genotypes that gain substantial benefits from symbiosis. Breeders have largely neglected symbiosis traits, but artificial selection for improved plant responses to microbiota could increase plant performance and sustainability. This article is protected by copyright. All rights reserved.

Fitness Effects of Mutations: An Assessment of PROVEAN Predictions Using Mutation Accumulation Data

Predicting fitness in natural populations is a major challenge in biology. It may be possible to leverage fast-accumulating genomic data sets to infer the fitness effects of mutant alleles, allowing evolutionary questions to be addressed in any organism. In this paper, we investigate the utility of one such tool, called PROVEAN. This program compares a query sequence with existing data to provide an alignment-based score for any protein variant, with scores categorized as neutral or deleterious based on a pre-set threshold. PROVEAN has been used widely in evolutionary studies, for example, to estimate mutation load in natural populations, but has not been formally tested as a predictor of aggregate mutational effects on fitness. Using three large published data sets on the genome sequences of laboratory mutation accumulation lines, we assessed how well PROVEAN predicted the actual fitness patterns observed, relative to other metrics. In most cases, we find that a simple count of the total number of mutant proteins is a better predictor of fitness than the number of proteins with variants scored as deleterious by PROVEAN. We also find that the sum of all mutant protein scores explains variation in fitness better than the number of mutant proteins in one of the data sets. We discuss the implications of these results for studies of populations in the wild.

Phylogenomic analysis points to a South American origin of Manihot and illuminates the primary gene pool of cassava

The genus Manihot, with around 120 known species, is native to a wide range of habitats and regions in the tropical and subtropical Americas. Its high species richness and recent diversification only c. 6 million years ago have significantly complicated previous phylogenetic analyses. Several basic elements of Manihot evolutionary history therefore remain unresolved. Here, we conduct a comprehensive phylogenomic analysis of Manihot, focusing on exhaustive sampling of South American taxa. We find that two recently described species from northeast Brazil's Atlantic Forest were the earliest to diverge, strongly suggesting a South American common ancestor of Manihot. Ancestral state reconstruction indicates early Manihot diversification in dry forests, with numerous independent episodes of new habitat colonization, including into savannas and rainforests within South America. We identify the closest wild relatives to Manihot esculenta, including the crop cassava, and we quantify extensive wild introgression into the cassava gene pool from at least five wild species, including Manihot glaziovii, a species used widely in breeding programs. Finally, we show that this wild-to-crop introgression substantially shapes the mutation load in cassava. Our findings provide a detailed case study for neotropical evolutionary history in a diverse and widespread group, and a robust phylogenomic framework for future Manihot and cassava research.

Deleterious protein-coding variants in diverse cattle breeds of the world

The domestication of wild animals has resulted in a reduction in effective population sizes, which can affect the deleterious mutation load of domesticated breeds. In addition, artificial selection contributes to the accumulation of deleterious mutations because of an increased rate of inbreeding among domesticated animals. Since founder population sizes and artificial selection differ between cattle breeds, their deleterious mutation load can vary. We investigated this question by using whole-genome data from 432 animals belonging to 54 worldwide cattle breeds. Our analysis revealed a negative correlation between genomic heterozygosity and nonsynonymous-to-silent diversity ratio, which suggests a higher proportion of single nucleotide variants (SNVs) affecting proteins in low-diversity breeds. Our results also showed that low-diversity breeds had a larger number of high-frequency (derived allele frequency (DAF) > 0.51) deleterious SNVs than high-diversity breeds. An opposite trend was observed for the low-frequency (DAF ≤ 0.51) deleterious SNVs. Overall, the number of high-frequency deleterious SNVs was larger in the genomes of taurine cattle breeds than of indicine breeds, whereas the number of low-frequency deleterious SNVs was larger in the genomes of indicine cattle than in those of taurine cattle. Furthermore, we observed significant variation in the counts of deleterious SNVs within taurine breeds. The variations in deleterious mutation load between taurine and indicine breeds could be attributed to the population sizes of the wild progenitors before domestication, whereas the variations observed within taurine breeds could be due to differences in inbreeding level, strength of artificial selection, and/or founding population size. Our findings imply that the incidence of genetic diseases can vary between cattle breeds.

Interplay between extreme drift and selection intensities favors the fixation of beneficial mutations in selfing maize populations

Population and quantitative genetic models provide useful approximations to predict long-term selection responses sustaining phenotypic shifts, and underlying multilocus adaptive dynamics. Valid across a broad range of parameters, their use for understanding the adaptive dynamics of small selfing populations undergoing strong selection intensity (thereafter High Drift-High selection regime, HDHS) remains to be explored. Saclay Divergent Selection Experiments (DSEs) on maize flowering time provide an interesting example of populations evolving under HDHS, with significant selection responses over 20 generations in two directions. We combined experimental data from Saclay DSEs, forward individual-based simulations, and theoretical predictions to dissect the evolutionary mechanisms at play in the observed selection responses. We asked two main questions: How do mutations arise, spread, and reach fixation in populations evolving under HDHS? How does the interplay between drift and selection influence observed phenotypic shifts? We showed that the long-lasting response to selection in small populations is due to the rapid fixation of mutations occurring during the generations of selection. Among fixed mutations, we also found a clear signal of enrichment for beneficial mutations revealing a limited cost of selection. Both environmental stochasticity and variation in selection coefficients likely contributed to exacerbate mutational effects, thereby facilitating selection grasp and fixation of small-effect mutations. Together our results highlight that despite a small number of polymorphic loci expected under HDHS, adaptive variation is continuously fueled by a vast mutational target. We discuss our results in the context of breeding and long-term survival of small selfing populations.

Genomic signatures of inbreeding and mutation load in a threatened rattlesnake

Theory predicts that threatened species living in small populations will experience high levels of inbreeding that will increase their genetic load, but recent work suggests that the impact of load may be minimized by purging resulting from long-term population bottlenecks. Empirical studies that examine this idea using genome-wide estimates of inbreeding and genetic load in threatened species are limited. Here we use individual genome resequencing data to compare levels of inbreeding, levels of genetic load (estimated as mutation load) and population history in threatened Eastern massasauga rattlesnakes (Sistrurus catenatus), which exist in small isolated populations, and closely related yet outbred Western massasauga rattlesnakes (Sistrurus tergeminus). In terms of inbreeding, S. catenatus genomes had a greater number of runs of homozygosity of varying sizes, indicating sustained inbreeding through repeated bottlenecks when compared to S. tergeminus. At the species level, outbred S. tergeminus had higher genome-wide levels of mutation load in the form of greater numbers of derived deleterious mutations compared to S. catenatus, presumably due to long-term purging of deleterious mutations in S. catenatus. In contrast, mutations that escaped species-level drift effects within S. catenatus populations were in general more frequent and more often found in homozygous genotypes than in S. tergeminus, suggesting a reduced efficiency of purifying selection in smaller S. catenatus populations for most mutations. Our results support an emerging idea that the historical demography of a threatened species has a significant impact on the type of genetic load present, which impacts implementation of conservation actions such as genetic rescue.

Harmful mutation load in the mitochondrial genomes of cattle breeds

Objective

Domestication of wild animals results in a reduction in the effective population size, and this could affect the deleterious mutation load of domesticated breeds. Furthermore, artificial selection will also contribute to the accumulation of deleterious mutations due to the increased rate of inbreeding among these animals. The process of domestication, founder population size, and artificial selection differ between cattle breeds, which could lead to a variation in their deleterious mutation loads. We investigated this using mitochondrial genome data from 364 animals belonging to 18 cattle breeds of the world.

Results

Our analysis revealed more than a fivefold difference in the deleterious mutation load among cattle breeds. We also observed a negative correlation between the breed age and the proportion of deleterious amino acid-changing polymorphisms. This suggests a proportionally higher deleterious SNPs in young breeds compared to older breeds. Our results highlight the magnitude of difference in the deleterious mutations present in the mitochondrial genomes of various breeds. The results of this study could be useful in predicting the rate of incidence of genetic diseases in different breeds.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-021-05664-y.

Ancestral sequence reconstruction - An underused approach to understand the evolution of gene function in plants?

Whilst substantial research effort has been placed on understanding the interactions of plant proteins with their molecular partners, relatively few studies in plants - by contrast to work in other organisms - address how these interactions evolve. It is thought that ancestral proteins were more promiscuous than modern proteins and that specificity often evolved following gene duplication and subsequent functional refining. However, ancestral protein resurrection studies have found that some modern proteins have evolved de novo from ancestors lacking those functions. Intriguingly, the new interactions evolved as a consequence of just a few mutations and, as such, acquisition of new functions appears to be neither difficult nor rare, however, only a few of them are incorporated into biological processes before they are lost to subsequent mutations. Here, we detail the approach of ancestral sequence reconstruction (ASR), providing a primer to reconstruct the sequence of an ancestral gene. We will present case studies from a range of different eukaryotes before discussing the few instances where ancestral reconstructions have been used in plants. As ASR is used to dig into the remote evolutionary past, we will also present some alternative genetic approaches to investigate molecular evolution on shorter timescales. We argue that the study of plant secondary metabolism is particularly well suited for ancestral reconstruction studies. Indeed, its ancient evolutionary roots and highly diverse landscape provide an ideal context in which to address the focal issue around the emergence of evolutionary novelties and how this affects the chemical diversification of plant metabolism.

Ancient genomes reveal early Andean farmers selected common beans while preserving diversity

All crops are the product of a domestication process that started less than 12,000 years ago from one or more wild populations^1,2. Farmers selected desirable phenotypic traits (such as improved energy accumulation, palatability of seeds and reduced natural shattering³) while leading domesticated populations through several more or less gradual demographic contractions^2,4. As a consequence, the erosion of wild genetic variation⁵ is typical of modern cultivars, making them highly susceptible to pathogens, pests and environmental change^6,7. The loss of genetic diversity hampers further crop improvement programmes to increase food production in a changing world, posing serious threats to food security^8,9. Using both ancient and modern seeds, we analysed the temporal dynamics of genetic variation and selection during the domestication process of the common bean (Phaseolus vulgaris) in the southern Andes. Here, we show that most domestic traits were selected for before 2,500 years ago, with no or only minor loss of whole-genome heterozygosity. In fact, most of the changes at coding genes and linked regions that differentiate wild and domestic genomes are already present in the ancient genomes analysed here, and all ancient domestic genomes dated between 600 and 2,500 years ago are highly variable (at least as variable as modern genomes from the wild). Single seeds from modern cultivars show reduced variation when compared with ancient seeds, indicating that intensive selection within cultivars in the past few centuries probably partitioned ancestral variation within different genetically homogenous cultivars. When cultivars from different Andean regions are pooled, the genomic variation of the pool is higher than that observed in the pool of ancient seeds from north and central western Argentina. Considering that most desirable phenotypic traits are probably controlled by multiple polymorphic genes¹⁰, a plausible explanation of this decoupling of selection and genetic erosion is that early farmers applied a relatively weak selection pressure² by using many phenotypically similar but genetically diverse individuals as parents. Our results imply that selection strategies during the past few centuries, as compared with earlier times, more intensively reduced genetic variation within cultivars and produced further improvements by focusing on a few plants carrying the traits of interest, at the cost of marked genetic erosion within Andean landraces.

Comparative analysis of chloroplast genomes of seven perennial Helianthus species

Sequencing and a comparative analysis of the complete chloroplast genomes of seven perennial Helianthus species were carried out. The chloroplast genomes have a typical quadripartite structure, including large and small single regions and a pair of inverted repeats. Genome sizes were between 151,152 bp and 151,289 bp. The genome of H. grosseserratus was the smallest, while that of H. microcephalus was the largest. The size variation of the chloroplast genomes is substantially due to the change in the length of simple sequence repeats (SSRs) in non-coding regions. An analysis of these SSRs revealed 35 polymorphic loci (average PIC value > 0.5) that can be used to examine ecological and evolutionary processes in wild Helianthus species. Eight divergence hotspots, including five intergenic regions (petN-psbM, clpP intron, rps3-rpl16, ndhD-ccsA, and ndhF-rpl32) and three gene regions (rbcL, ycf1, and ndhF) were also identified in Helianthus chloroplast genomes. The evolutionary selection pressure analysis revealed a strong purifying selection. Only the rbcL gene experienced efficiency of positive selection at the annual/perennial transitions. The inverted repeat (IR)/single copy (SC) boundaries were identical in all of these (Helianthus) species. In general, the comparison of the genomes revealed low levels of sequence variability (Pi = 0.00051). This indicates that the chloroplast genomes of the studied perennial species of Helianthus, in addition to purifying selection, are closely related and have a recent divergence time.

The patterns of deleterious mutations during the domestication of soybean

Globally, soybean is a major protein and oil crop. Enhancing our understanding of the soybean domestication and improvement process helps boost genomics-assisted breeding efforts. Here we present a genome-wide variation map of 10.6 million single-nucleotide polymorphisms and 1.4 million indels for 781 soybean individuals which includes 418 domesticated (Glycine max), 345 wild (Glycine soja), and 18 natural hybrid (G. max/G. soja) accessions. We describe the enhanced detection of 183 domestication-selective sweeps and the patterns of putative deleterious mutations during domestication and improvement. This predominantly selfing species shows 7.1% reduction of overall deleterious mutations in domesticated soybean relative to wild soybean and a further 1.4% reduction from landrace to improved accessions. The detected domestication-selective sweeps also show reduced levels of deleterious alleles. Importantly, genotype imputation with this resource increases the mapping resolution of genome-wide association studies for seed protein and oil traits in a soybean diversity panel.

Harnessing Knowledge from Maize and Rice Domestication for New Crop Breeding

Crop domestication has fundamentally altered the course of human history, causing a shift from hunter-gatherer to agricultural societies and stimulating the rise of modern civilization. A greater understanding of crop domestication would provide a theoretical basis for how we could improve current crops and develop new crops to deal with environmental challenges in a sustainable manner. Here, we provide a comprehensive summary of the similarities and differences in the domestication processes of maize and rice, two major staple food crops that feed the world. We propose that maize and rice might have evolved distinct genetic solutions toward domestication. Maize and rice domestication appears to be associated with distinct regulatory and evolutionary mechanisms. Rice domestication tended to select de novo, loss-of-function, coding variation, while maize domestication more frequently favored standing, gain-of-function, regulatory variation. At the gene network level, distinct genetic paths were used to acquire convergent phenotypes in maize and rice domestication, during which different central genes were utilized, orthologous genes played different evolutionary roles, and unique genes or regulatory modules were acquired for establishing new traits. Finally, we discuss how the knowledge gained from past domestication processes, together with emerging technologies, could be exploited to improve modern crop breeding and domesticate new crops to meet increasing human demands.

Evidence for the Accumulation of Nonsynonymous Mutations and Favorable Pleiotropic Alleles During Wheat Breeding

Plant breeding leads to the genetic improvement of target traits by selecting a small number of genotypes from among typically large numbers of candidate genotypes after careful evaluation. In this study, we first investigated how mutations at conserved nucleotide sites normally viewed as deleterious, such as nonsynonymous sites, accumulated in a wheat, Triticum aestivum, breeding lineage. By comparing a 150 year old ancestral and modern cultivar, we found recent nucleotide polymorphisms altered amino acids and occurred within conserved genes at frequencies expected in the absence of purifying selection. Mutations that are deleterious in other contexts likely had very small or no effects on target traits within the breeding lineage. Second, we investigated if breeders selected alleles with favorable effects on some traits and unfavorable effects on others and used different alleles to compensate for the latter. An analysis of a segregating population derived from the ancestral and modern parents provided one example of this phenomenon. The recent cultivar contains the Rht-B1b green revolution semi-dwarfing allele and compensatory alleles that reduce its negative effects. However, improvements in traits other than plant height were due to pleiotropic loci with favorable effects on traits and to favorable loci with no detectable pleiotropic effects. Wheat breeding appears to tolerate mutations at conserved nucleotide sites and to only select for alleles with both favorable and unfavorable effects on traits in exceptional situations.

Adaptive and maladaptive genetic diversity in small populations: Insights from the Brook Charr (Salvelinus fontinalis) case study

Investigating the relative importance of neutral versus selective processes governing the accumulation of genetic variants is a key goal in both evolutionary and conservation biology. This is particularly true in the context of small populations, where genetic drift can counteract the effect of selection. Using Brook Charr (Salvelinus fontinalis) from Québec, Canada, as a case study, we investigated the importance of demographic versus selective processes governing the accumulation of both adaptive and maladaptive mutations in closed versus open and connected populations to assess gene flow effect. This was achieved by using 14,779 high-quality filtered SNPs genotyped among 1,416 fish representing 50 populations from three life history types: lacustrine (closed populations), riverine and anadromous (connected populations). Using the PROVEAN algorithm, we observed a considerable accumulation of putative deleterious mutations across populations. The absence of correlation between the occurrence of putatively beneficial or deleterious mutations and local recombination rate supports the hypothesis that genetic drift might be the main driver of the accumulation of such variants. However, despite a lower genetic diversity observed in lacustrine than in riverine or anadromous populations, lacustrine populations do not exhibit more deleterious mutations than the two other history types, suggesting that the negative effect of genetic drift in lacustrine populations may be mitigated by that of relaxed purifying selection. Moreover, we also identified genomic regions associated with anadromy, as well as an overrepresentation of transposable elements associated with variation in environmental variables, thus supporting the importance of transposable elements in adaptation.

A draft genome of sweet cherry (Prunus avium L.) reveals genome-wide and local effects of domestication

Sweet cherry (Prunus avium L.) trees are both economically important fruit crops but also important components of natural forest ecosystems in Europe, Asia and Africa. Wild and domesticated trees currently coexist in the same geographic areas with important questions arising on their historical relationships. Little is known about the effects of the domestication process on the evolution of the sweet cherry genome. We assembled and annotated the genome of the cultivated variety "Big Star*" and assessed the genetic diversity among 97 sweet cherry accessions representing three different stages in the domestication and breeding process (wild trees, landraces and modern varieties). The genetic diversity analysis revealed significant genome-wide losses of variation among the three stages and supports a clear distinction between wild and domesticated trees, with only limited gene flow being detected between wild trees and domesticated landraces. We identified 11 domestication sweeps and five breeding sweeps covering, respectively, 11.0 and 2.4 Mb of the P. avium genome. A considerable fraction of the domestication sweeps overlaps with those detected in the related species, Prunus persica (peach), indicating that artificial selection during domestication may have acted independently on the same regions and genes in the two species. We detected 104 candidate genes in sweep regions involved in different processes, such as the determination of fruit texture, the regulation of flowering and fruit ripening and the resistance to pathogens. The signatures of selection identified will enable future evolutionary studies and provide a valuable resource for genetic improvement and conservation programs in sweet cherry.

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

Survival in the Tropics despite isolation, inbreeding and asexual reproduction: insights from the genome of the world's southernmost poplar (Populus ilicifolia)

Species are becoming extinct at unprecedented rates as a consequence of human activity. Hence it is important to understand the evolutionary dynamics of species with already small population sizes. Populus ilicifolia is a vulnerable poplar species that is isolated from other poplar species and is uniquely adapted to the Tropics. It has a very limited size, reproduces partly clonally and is therefore an excellent case study for conservation genomics. We present here the first annotated draft genome of P. ilicifolia, characterize genome-wide patterns of polymorphisms and compare those to other poplar species with larger natural ranges. P. ilicifolia experienced a more prolonged and severe decline of effective population size (N_e ) and signs of genetic erosion than any other poplar species with which it was compared. At present, the species has the lowest genome-wide genetic diversity, the highest abundance of long runs of homozygosity, high inbreeding levels as well as a high overall accumulation of deleterious variants. However, more effective purging of severely deleterious variants and adaptation to the Tropics may have contributed to its survival. Hence, in spite of its limited genetic variation, it is certainly worth pursuing the conservation efforts of this unique species.

Accumulation of Deleterious Mutations in Landlocked Threespine Stickleback Populations

Colonization of new habitats often reduces population sizes and may result in the accumulation of deleterious mutations by genetic drift. Compared with the genomic basis for adaptation to new environments, genome-wide analysis of deleterious mutations in isolated populations remains limited. In the present study, we investigated the accumulation of deleterious mutations in five endangered freshwater populations of threespine stickleback (Gasterosteus aculeatus) in the central part of the mainland of Japan. Using whole-genome resequencing data, we first conducted phylogenomic analysis and confirmed at least two independent freshwater colonization events in the central mainland from ancestral marine ecotypes. Next, analyses of single nucleotide polymorphisms showed a substantial reduction of heterozygosity in freshwater populations compared with marine populations. Reduction in heterozygosity was more apparent at the center of each chromosome than the peripheries and on X chromosomes compared with autosomes. Third, bioinformatic analysis of deleterious mutations showed increased accumulation of putatively deleterious mutations in the landlocked freshwater populations compared with marine populations. For the majority of populations examined, the frequencies of putatively deleterious mutations were higher on X chromosomes than on autosomes. The interpopulation comparison indicated that the majority of putatively deleterious mutations may have accumulated independently. Thus, whole-genome resequencing of endangered populations can help to estimate the accumulation of deleterious mutations and inform us of which populations are the most severely endangered. Furthermore, analysis of variation among chromosomes can give insights into whether any particular chromosomes are likely to accumulate deleterious mutations.

From Summary Statistics to Gene Trees: Methods for Inferring Positive Selection

Methods to detect signals of natural selection from genomic data have traditionally emphasized the use of simple summary statistics. Here, we review a new generation of methods that consider combinations of conventional summary statistics and/or richer features derived from inferred gene trees and ancestral recombination graphs (ARGs). We also review recent advances in methods for population genetic simulation and ARG reconstruction. Finally, we describe opportunities for future work on a variety of related topics, including the genetics of speciation, estimation of selection coefficients, and inference of selection on polygenic traits. Together, these emerging methods offer promising new directions in the study of natural selection.

Genome-wide analysis reveals molecular convergence underlying domestication in 7 bird and mammals

Background

In response to ecological niche of domestication, domesticated mammals and birds developed adaptively phenotypic homoplasy in behavior modifications like fearlessness, altered sociability, exploration and cognition, which partly or indirectly result in consequences for economic productivity. Such independent adaptations provide an excellent model to investigate molecular mechanisms and patterns of evolutionary convergence driven by artificial selection.

Results

First performing population genomic and brain transcriptional comparisons in 68 wild and domesticated chickens, we revealed evolutionary trajectories, genetic architectures and physiologic bases of adaptively behavioral alterations. To extensively decipher molecular convergence on behavioral changes thanks to domestication, we investigated selection signatures in hundreds of genomes and brain transcriptomes across chicken and 6 other domesticated mammals. Although no shared substitution was detected, a common enrichment of the adaptive mutations in regulatory sequences was observed, presenting significance to drive adaptations. Strong convergent pattern emerged at levels of gene, gene family, pathway and network. Genes implicated in neurotransmission, semaphorin, tectonic protein and modules regulating neuroplasticity were central focus of selection, supporting molecular repeatability of homoplastic behavior reshapes. Genes at nodal positions in trans-regulatory networks were preferably targeted. Consistent down-regulation of majority brain genes may be correlated with reduced brain size during domestication. Up-regulation of splicesome genes in chicken rather mammals highlights splicing as an efficient way to evolve since avian-specific genomic contraction of introns and intergenics. Genetic burden of domestication elicits a general hallmark. The commonly selected genes were relatively evolutionary conserved and associated with analogous neuropsychiatric disorders in human, revealing trade-off between adaption to life with human at the cost of neural changes affecting fitness in wild.

Conclusions

After a comprehensive investigation on genomic diversity and evolutionary trajectories in chickens, we revealed basis, pattern and evolutionary significance of molecular convergence in domesticated bird and mammals, highlighted the genetic basis of a compromise on utmost adaptation to the lives with human at the cost of high risk of neurophysiological changes affecting animals’ fitness in wild.

Purging of highly deleterious mutations through severe bottlenecks in Alpine ibex

Human activity has caused dramatic population declines in many wild species. The resulting bottlenecks have a profound impact on the genetic makeup of a species with unknown consequences for health. A key genetic factor for species survival is the evolution of deleterious mutation load, but how bottleneck strength and mutation load interact lacks empirical evidence. We analyze 60 complete genomes of six ibex species and the domestic goat. We show that historic bottlenecks rather than the current conservation status predict levels of genome-wide variation. By analyzing the exceptionally well-characterized population bottlenecks of the once nearly extinct Alpine ibex, we find genomic evidence of concurrent purging of highly deleterious mutations but accumulation of mildly deleterious mutations. This suggests that recolonization bottlenecks induced both relaxed selection and purging, thus reshaping the landscape of deleterious mutation load. Our findings highlight that even populations of ~1000 individuals can accumulate mildly deleterious mutations. Conservation efforts should focus on preventing population declines below such levels to ensure long-term survival of species.

Although there is extensive theory predicting the effects of population bottlenecks on mutation load, there is little empirical evidence from recent bottlenecks. Here, Grossen et al. compare the consequences of population bottlenecks in six ibex species for genome-wide variation and mutation load.

The type of bottleneck matters: Insights into the deleterious variation landscape of small managed populations

Predictions about the consequences of a small population size on genetic and deleterious variation are fundamental to population genetics. As small populations are more affected by genetic drift, purifying selection acting against deleterious alleles is predicted to be less efficient, therefore increasing the risk of inbreeding depression. However, the extent to which small populations are subjected to genetic drift depends on the nature and time frame in which the bottleneck occurs. Domesticated species are an excellent model to investigate the consequences of population bottlenecks on genetic and deleterious variation in small populations. This is because their history is dominated by known bottlenecks associated with domestication, breed formation and intense selective breeding. Here, we use whole-genome sequencing data from 97 chickens representing 39 traditional fancy breeds to directly examine the consequences of two types of bottlenecks for deleterious variation: the severe domestication bottleneck and the recent population decline accompanying breed formation. We find that recently bottlenecked populations have a higher proportion of deleterious variants relative to populations that have been kept at small population sizes since domestication. We also observe that long tracts of homozygous genotypes (runs of homozygosity) are proportionally more enriched in deleterious variants than the rest of the genome. This enrichment is particularly evident in recently bottlenecked populations, suggesting that homozygosity of these variants is likely to occur due to genetic drift and recent inbreeding. Our results indicate that the timing and nature of population bottlenecks can substantially shape the deleterious variation landscape in small populations.

Pod indehiscence is a domestication and aridity resilience trait in common bean

Plant domestication has strongly modified crop morphology and development. Nevertheless, many crops continue to display atavistic characteristics that were advantageous to their wild ancestors but are deleterious under cultivation, such as pod dehiscence (PD). Here, we provide the first comprehensive assessment of the inheritance of PD in the common bean (Phaseolus vulgaris), a major domesticated grain legume. Using three methods to evaluate the PD phenotype, we identified multiple, unlinked genetic regions controlling PD in a biparental population and two diversity panels. Subsequently, we assessed patterns of orthology among these loci and those controlling the trait in other species. Our results show that different genes were selected in each domestication and ecogeographic race. A chromosome Pv03 dirigent-like gene, involved in lignin biosynthesis, showed a base-pair substitution that is associated with decreased PD. This haplotype may underlie the expansion of Mesoamerican domesticates into northern Mexico, where arid conditions promote PD. The rise in frequency of the decreased-PD haplotype may be a consequence of the markedly different fitness landscape imposed by domestication. Environmental dependency and genetic redundancy can explain the maintenance of atavistic traits under domestication.

The Fate of Deleterious Variants in a Barley Genomic Prediction Population

Targeted identification and purging of deleterious genetic variants has been proposed as a novel approach to animal and plant breeding. This strategy is motivated, in part, by the observation that demographic events and strong selection associated with cultivated species pose a "cost of domestication." This includes an increase in the proportion of genetic variants that are likely to reduce fitness. Recent advances in DNA resequencing and sequence constraint-based approaches to predict the functional impact of a mutation permit the identification of putatively deleterious SNPs (dSNPs) on a genome-wide scale. Using exome capture resequencing of 21 barley lines, we identified 3855 dSNPs among 497,754 total SNPs. We generated whole-genome resequencing data of Hordeum murinum ssp. glaucum as a phylogenetic outgroup to polarize SNPs as ancestral vs. derived. We also observed a higher proportion of dSNPs per synonymous SNPs (sSNPs) in low-recombination regions of the genome. Using 5215 progeny from a genomic prediction experiment, we examined the fate of dSNPs over three breeding cycles. Adjusting for initial frequency, derived alleles at dSNPs reduced in frequency or were lost more often than other classes of SNPs. The highest-yielding lines in the experiment, as chosen by standard genomic prediction approaches, carried fewer homozygous dSNPs than randomly sampled lines from the same progeny cycle. In the final cycle of the experiment, progeny selected by genomic prediction had a mean of 5.6% fewer homozygous dSNPs relative to randomly chosen progeny from the same cycle.

The role of recombination on genome-wide patterns of local ancestry exemplified by supplemented brook charr populations

Assessing the immediate and long-term evolutionary consequences of human-mediated hybridization is of major concern for conservation biology. Several studies have documented how selection in interaction with recombination modulates introgression at a genome-wide scale, but few have considered the dynamics of this process within and among chromosomes. Here, we used an exploited freshwater fish, the brook charr (Salvelinus fontinalis), for which decades of stocking practices have resulted in admixture between wild populations and an introduced domestic strain, to assess both the temporal dynamics and local chromosomal variation in domestic ancestry. We provide a detailed picture of the domestic ancestry patterns across the genome using about 33,000 mapped single nucleotide polymorphisms genotyped in 611 individuals from 24 supplemented populations. For each lake, we distinguished early- and late-generation hybrids using information regarding admixture tracts. To assess the selective outcomes following admixture we then evaluated the relationship between recombination and admixture proportions at three different scales: the whole genome, chromosomes and within 2-Mb windows. This allowed us to detect a wide range of evolutionary mechanisms varying along the genome, as reflected by the finding of favoured or disfavoured introgression of domestic haplotypes. Among these, the main factor modulating local ancestry was probably the presence of deleterious recessive mutations in the wild populations, which can be efficiently hidden to selection in the presence of long admixture tracts. Overall, our results emphasize the relevance of taking into consideration local ancestry information to assess both the temporal and the chromosomal variation in local admixture ancestry toward better understanding post-hybridization evolutionary outcomes.