Massive haplotypes underlie ecotypic differentiation in sunflowers

Genomic architecture of adaptive radiation and hybridization in Alpine whitefish

Adaptive radiations represent some of the most remarkable explosions of diversification across the tree of life. However, the constraints to rapid diversification and how they are sometimes overcome, particularly the relative roles of genetic architecture and hybridization, remain unclear. Here, we address these questions in the Alpine whitefish radiation, using a whole-genome dataset that includes multiple individuals of each of the 22 species belonging to six ecologically distinct ecomorph classes across several lake-systems. We reveal that repeated ecological and morphological diversification along a common environmental axis is associated with both genome-wide allele frequency shifts and a specific, larger effect, locus, associated with the gene edar. Additionally, we highlight the possible role of introgression between species from different lake-systems in facilitating the evolution and persistence of species with unique trait combinations and ecology. These results highlight the importance of both genome architecture and secondary contact with hybridization in fuelling adaptive radiation.

In this genomic study on Alpine whitefish radiations, the authors reveal details on the genetic architecture underlying the repeated eco-morphological diversification and the role of hybridization in the evolution of endemic whitefish species.

Association mapping of colour variation in a butterfly provides evidence that a supergene locks together a cluster of adaptive loci

Supergenes are genetic architectures associated with discrete and concerted variation in multiple traits. It has long been suggested that supergenes control these complex polymorphisms by suppressing recombination between sets of coadapted genes. However, because recombination suppression hinders the dissociation of the individual effects of genes within supergenes, there is still little evidence that supergenes evolve by tightening linkage between coadapted genes. Here, combining a landmark-free phenotyping algorithm with multivariate genome-wide association studies, we dissected the genetic basis of wing pattern variation in the butterfly Heliconius numata. We show that the supergene controlling the striking wing pattern polymorphism displayed by this species contains several independent loci associated with different features of wing patterns. The three chromosomal inversions of this supergene suppress recombination between these loci, supporting the hypothesis that they may have evolved because they captured beneficial combinations of alleles. Some of these loci are, however, associated with colour variations only in a subset of morphs where the phenotype is controlled by derived inversion forms, indicating that they were recruited after the formation of the inversions. Our study shows that supergenes and clusters of adaptive loci in general may form via the evolution of chromosomal rearrangements suppressing recombination between co-adapted loci but also via the subsequent recruitment of linked adaptive mutations. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Three problems in the genetics of speciation by selection

Speciation is the process by which barriers to gene flow evolve between populations. Although we now know that speciation is largely driven by natural selection, knowledge of the agents of selection and the genetic and genomic mechanisms that facilitate divergence is required for a satisfactory theory of speciation. In this essay, we highlight three advances/problems in our understanding of speciation that have arisen from studies of the genes and genomic regions that underlie the evolution of reproductive isolation. First, we describe how the identification of “speciation” genes makes it possible to identify the agents of selection causing the evolution of reproductive isolation, while also noting that the link between the genetics of phenotypic divergence and intrinsic postzygotic reproductive barriers remains tenuous. Second, we discuss the important role of recombination suppressors in facilitating speciation with gene flow, but point out that the means and timing by which reproductive barriers become associated with recombination cold spots remains uncertain. Third, we establish the importance of ancient genetic variation in speciation, although we argue that the focus of speciation studies on evolutionarily young groups may bias conclusions in favor of ancient variation relative to new mutations.

A Quantitative Genetic Study of Sclerotinia Head Rot Resistance Introgressed from the Wild Perennial Helianthus maximiliani into Cultivated Sunflower (Helianthus annuus L.)

Sclerotinia head rot (HR), caused by Sclerotinia sclerotiorum, is an economically important disease of sunflower with known detrimental effects on yield and quality in humid climates worldwide. The objective of this study was to gain insight into the genetic architecture of HR resistance from a sunflower line HR21 harboring HR resistance introgressed from the wild perennial Helianthus maximiliani. An F2 population derived from the cross of HA 234 (susceptible-line)/HR21 (resistant-line) was evaluated for HR resistance at two locations during 2019−2020. Highly significant genetic variations (p < 0.001) were observed for HR disease incidence (DI) and disease severity (DS) in both individual and combined analyses. Broad sense heritability (H2) estimates across environments for DI and DS were 0.51 and 0.62, respectively. A high-density genetic map of 1420.287 cM was constructed with 6315 SNP/InDel markers developed using genotype-by-sequencing technology. A total of 16 genomic regions on eight sunflower chromosomes, 1, 2, 10, 12, 13, 14, 16 and 17 were associated with HR resistance, each explaining between 3.97 to 16.67% of the phenotypic variance for HR resistance. Eleven of these QTL had resistance alleles from the HR21 parent. Molecular markers flanking the QTL will facilitate marker-assisted selection breeding for HR resistance in sunflower.

Rapid adaptive radiation of Darwin's finches depends on ancestral genetic modules

Recent adaptive radiations are models for investigating mechanisms contributing to the evolution of biodiversity. An unresolved question is the relative importance of new mutations, ancestral variants, and introgressive hybridization for phenotypic evolution and speciation. Here, we address this issue using Darwin's finches and investigate the genomic architecture underlying their phenotypic diversity. Admixture mapping for beak and body size in the small, medium, and large ground finches revealed 28 loci showing strong genetic differentiation. These loci represent ancestral haplotype blocks with origins predating speciation events during the Darwin's finch radiation. Genes expressed in the developing beak are overrepresented in these genomic regions. Ancestral haplotypes constitute genetic modules for selection and act as key determinants of the unusual phenotypic diversity of Darwin's finches. Such ancestral haplotype blocks can be critical for how species adapt to environmental variability and change.

Dispersed emergence and protracted domestication of polyploid wheat uncovered by mosaic ancestral haploblock inference

Major crops are all survivors of domestication bottlenecks. Studies have focused on the genetic loci related to the domestication syndrome, while the contribution of ancient haplotypes remains largely unknown. Here, an ancestral genomic haploblock dissection method is developed and applied to a resequencing dataset of 386 tetraploid/hexaploid wheat accessions, generating a pan-ancestry haploblock map. Together with cytoplastic evidences, we reveal that domesticated polyploid wheat emerged from the admixture of six founder wild emmer lineages, which contributed the foundation of ancestral mosaics. The key domestication-related loci, originated over a wide geographical range, were gradually pyramided through a protracted process. Diverse stable-inheritance ancestral haplotype groups of the chromosome central zone are identified, revealing the expanding routes of wheat and the trends of modern wheat breeding. Finally, an evolution model of polyploid wheat is proposed, highlighting the key role of wild-to-crop and interploidy introgression, that increased genomic diversity following bottlenecks introduced by domestication and polyploidization.

The contribution of ancient haplotypes to domestication is largely unknown. Here, the authors develop an ancestral genomic haploblock dissection method to generate a mosaic pan-ancestry genomic map and reveal that the domesticated polyploidy wheat emerged from the admixture of six founder wild emmer linages.

HAPPE: A Tool for Population Haplotype Analysis and Visualization in Editable Excel Tables

Haplotype identification, characterization and visualization are important for large-scale analysis and use in population genomics. Many tools have been developed to visualize haplotypes, but it is challenging to display both the pattern of haplotypes and the genotypes for each single SNP in the context of a large amount of genomic data. Here, we describe the tool HAPPE, which uses the agglomerative hierarchical clustering algorithm to characterize and visualize the genotypes and haplotypes in a phylogenetic context. The tool displays the plots by coloring the cells and/or their borders in Excel tables for any given gene and genomic region of interest. HAPPE facilitates informative displays wherein data in plots are easy to read and access. It allows parallel display of several lines of values, such as phylogenetic trees, P values of GWAS, the entry of genes or SNPs, and the sequencing depth at each position. These features are informative for the detection of insertion/deletions or copy number variations. Overall, HAPPE provides editable plots consisting of cells in Excel tables, which are user-friendly to non-programmers. This pipeline is coded in Python and is available at https://github.com/fengcong3/HAPPE.

Sheltering of deleterious mutations explains the stepwise extension of recombination suppression on sex chromosomes and other supergenes

Many organisms have sex chromosomes with large nonrecombining regions that have expanded stepwise, generating "evolutionary strata" of differentiation. The reasons for this remain poorly understood, but the principal hypotheses proposed to date are based on antagonistic selection due to differences between sexes. However, it has proved difficult to obtain empirical evidence of a role for sexually antagonistic selection in extending recombination suppression, and antagonistic selection has been shown to be unlikely to account for the evolutionary strata observed on fungal mating-type chromosomes. We show here, by mathematical modeling and stochastic simulation, that recombination suppression on sex chromosomes and around supergenes can expand under a wide range of parameter values simply because it shelters recessive deleterious mutations, which are ubiquitous in genomes. Permanently heterozygous alleles, such as the male-determining allele in XY systems, protect linked chromosomal inversions against the expression of their recessive mutation load, leading to the successive accumulation of inversions around these alleles without antagonistic selection. Similar results were obtained with models assuming recombination-suppressing mechanisms other than chromosomal inversions and for supergenes other than sex chromosomes, including those without XY-like asymmetry, such as fungal mating-type chromosomes. However, inversions capturing a permanently heterozygous allele were found to be less likely to spread when the mutation load segregating in populations was lower (e.g., under large effective population sizes or low mutation rates). This may explain why sex chromosomes remain homomorphic in some organisms but are highly divergent in others. Here, we model a simple and testable hypothesis explaining the stepwise extensions of recombination suppression on sex chromosomes, mating-type chromosomes, and supergenes in general.

The Relative Power of Structural Genomic Variation versus SNPs in Explaining the Quantitative Trait Growth in the Marine Teleost Chrysophrys auratus

Background: Genetic diversity provides the basic substrate for evolution. Genetic variation consists of changes ranging from single base pairs (single-nucleotide polymorphisms, or SNPs) to larger-scale structural variants, such as inversions, deletions, and duplications. SNPs have long been used as the general currency for investigations into how genetic diversity fuels evolution. However, structural variants can affect more base pairs in the genome than SNPs and can be responsible for adaptive phenotypes due to their impact on linkage and recombination. In this study, we investigate the first steps needed to explore the genetic basis of an economically important growth trait in the marine teleost finfish Chrysophrys auratus using both SNP and structural variant data. Specifically, we use feature selection methods in machine learning to explore the relative predictive power of both types of genetic variants in explaining growth and discuss the feature selection results of the evaluated methods. Methods: SNP and structural variant callers were used to generate catalogues of variant data from 32 individual fish at ages 1 and 3 years. Three feature selection algorithms (ReliefF, Chi-square, and a mutual-information-based method) were used to reduce the dataset by selecting the most informative features. Following this selection process, the subset of variants was used as features to classify fish into small, medium, or large size categories using KNN, naïve Bayes, random forest, and logistic regression. The top-scoring features in each feature selection method were subsequently mapped to annotated genomic regions in the zebrafish genome, and a permutation test was conducted to see if the number of mapped regions was greater than when random sampling was applied. Results: Without feature selection, the prediction accuracies ranged from 0 to 0.5 for both structural variants and SNPs. Following feature selection, the prediction accuracy increased only slightly to between 0 and 0.65 for structural variants and between 0 and 0.75 for SNPs. The highest prediction accuracy for the logistic regression was achieved for age 3 fish using SNPs, although generally predictions for age 1 and 3 fish were very similar (ranging from 0–0.65 for both SNPs and structural variants). The Chi-square feature selection of SNP data was the only method that had a significantly higher number of matches to annotated genomic regions of zebrafish than would be explained by chance alone. Conclusions: Predicting a complex polygenic trait such as growth using data collected from a low number of individuals remains challenging. While we demonstrate that both SNPs and structural variants provide important information to help understand the genetic basis of phenotypic traits such as fish growth, the full complexities that exist within a genome cannot be easily captured by classical machine learning techniques. When using high-dimensional data, feature selection shows some increase in the prediction accuracy of classification models and provides the potential to identify unknown genomic correlates with growth. Our results show that both SNPs and structural variants significantly impact growth, and we therefore recommend that researchers interested in the genotype–phenotype map should strive to go beyond SNPs and incorporate structural variants in their studies as well. We discuss how our machine learning models can be further expanded to serve as a test bed to inform evolutionary studies and the applied management of species.

Three-nucleotide periodicity of nucleotide diversity in a population enables the identification of open reading frames

Accurate prediction of open reading frames (ORFs) is important for studying and using genome sequences. Ribosomes move along mRNA strands with a step of three nucleotides and datasets carrying this information can be used to predict ORFs. The ribosome-protected footprints (RPFs) feature a significant 3-nt periodicity on mRNAs and are powerful in predicting translating ORFs, including small ORFs (sORFs), but the application of RPFs is limited because they are too short to be accurately mapped in complex genomes. In this study, we found a significant 3-nt periodicity in the datasets of populational genomic variants in coding sequences, in which the nucleotide diversity increases every three nucleotides. We suggest that this feature can be used to predict ORFs and develop the Python package ‘OrfPP’, which recovers ~83% of the annotated ORFs in the tested genomes on average, independent of the population sizes and the complexity of the genomes. The novel ORFs, including sORFs, identified from single-nucleotide polymorphisms are supported by protein mass spectrometry evidence comparable to that of the annotated ORFs. The application of OrfPP to tetraploid cotton and hexaploid wheat genomes successfully identified 76.17% and 87.43% of the annotated ORFs in the genomes, respectively, as well as 4704 sORFs, including 1182 upstream and 2110 downstream ORFs in cotton and 5025 sORFs, including 232 upstream and 234 downstream ORFs in wheat. Overall, we propose an alternative and supplementary approach for ORF prediction that can extend the studies of sORFs to more complex genomes.

Advancing designer crops for climate resilience through an integrated genomics approach

Climate change and exponential population growth are exposing an immediate need for developing future crops that are highly resilient and adaptable to changing environments to maintain global food security in the next decade. Rigorous selection from long domestication history has rendered cultivated crops genetically disadvantaged, raising concerns in their ability to adapt to these new challenges and limiting their usefulness in breeding programmes. As a result, future crop improvement efforts must rely on integrating various genomic strategies ranging from high-throughput sequencing to machine learning, in order to exploit germplasm diversity and overcome bottlenecks created by domestication, expansive multi-dimensional phenotypes, arduous breeding processes, complex traits and big data.

Implications of Large-Effect Loci for Conservation: A Review and Case Study with Pacific Salmon

The increasing feasibility of assembling large genomic datasets for non-model species presents both opportunities and challenges for applied conservation and management. A popular theme in recent studies is the search for large-effect loci that explain substantial portions of phenotypic variance for a key trait(s). If such loci can be linked to adaptations, 2 important questions arise: 1) Should information from these loci be used to reconfigure conservation units (CUs), even if this conflicts with overall patterns of genetic differentiation? 2) How should this information be used in viability assessments of populations and larger CUs? In this review, we address these questions in the context of recent studies of Chinook salmon and steelhead (anadromous form of rainbow trout) that show strong associations between adult migration timing and specific alleles in one small genomic region. Based on the polygenic paradigm (most traits are controlled by many genes of small effect) and genetic data available at the time showing that early-migrating populations are most closely related to nearby late-migrating populations, adult migration differences in Pacific salmon and steelhead were considered to reflect diversity within CUs rather than separate CUs. Recent data, however, suggest that specific alleles are required for early migration, and that these alleles are lost in populations where conditions do not support early-migrating phenotypes. Contrasting determinations under the US Endangered Species Act and the State of California's equivalent legislation illustrate the complexities of incorporating genomics data into CU configuration decisions. Regardless how CUs are defined, viability assessments should consider that 1) early-migrating phenotypes experience disproportionate risks across large geographic areas, so it becomes important to identify early-migrating populations that can serve as reliable sources for these valuable genetic resources; and 2) genetic architecture, especially the existence of large-effect loci, can affect evolutionary potential and adaptability.

Effects of the neo-X chromosome on genomic signatures of hybridization in Rumex hastatulus

Natural hybrid zones provide opportunities for studies of the evolution of reproductive isolation in wild populations. Although recent investigations have found that the formation of neo-sex chromosomes is associated with reproductive isolation, the mechanisms remain unclear in most cases. Here, we assess the contemporary structure of gene flow in the contact zone between largely allopatric cytotypes of the dioecious plant Rumex hastatulus, a species with evidence of sex chromosome turn-over. Males to the west of the Mississippi river, USA, have an X and a single Y chromosome, whereas populations to the east of the river have undergone a chromosomal rearrangement giving rise to a larger X and two Y chromosomes. Using reduced-representation sequencing, we provide evidence that hybrids form readily and survive multiple backcross generations in the field, demonstrating the potential for ongoing gene flow between the cytotypes. Cline analysis of each chromosome separately captured no signals of difference in cline shape between chromosomes. However, principal component regression revealed a significant increase in the contribution of individual SNPs to inter-cytotype differentiation on the neo-X chromosome, but no correlation with recombination rate. Cline analysis revealed that the only SNPs with significantly steeper clines than the genome average were located on the neo-X. Our data are consistent with a role for neo-sex chromosomes in reproductive isolation between R. hastatulus cytotypes. Our investigation highlights the importance of studying plant hybrid zones for understanding the evolution of sex chromosomes.

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.

Sex determination and sex chromosome evolution in land plants

Linnaeus's very first opus, written when he was 22 years old, dealt with the analogy that exists between plants and animals in how they 'propagate their species', and a revised version with a plate depicting the union of male and female Mercurialis annua plants became a foundational text on the sexuality of plants. The question how systems with separate males and females have evolved in sedentary organisms that appear ancestrally bisexual has fascinated biologists ever since. The phenomenon, termed dioecy, has important consequences for plant reproductive success and is of commercial interest since it affects seed quality and fruit production. This theme issue presents a series of articles that synthesize and challenge the current understanding of how plants achieve dioecy. The articles deal with a broad set of taxa, including Coccinia, Ginkgo, Mercurialis, Populus, Rumex and Silene, as well as overarching topics, such as the field's terminology, analogies with animal sex determination systems, evolutionary pathways to dioecy, dosage compensation, and the longevity of the two sexes. In this introduction, we focus on four topics, each addressed by several articles from different angles and with different conclusions. Our highlighting of unclear or controversial issues may help future studies to build on the current understanding and to ask new questions that will expand our knowledge of plant sexual systems. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Population differentiation and structural variation in the Manduca sexta genome across the United States

Many species that are extensively studied in the laboratory are less well characterized in their natural habitat, and laboratory strains represent only a small fraction of the variation in a species’ genome. Here we investigate genomic variation in 3 natural North American populations of an agricultural pest and a model insect for many scientific disciplines, the tobacco hornworm (Manduca sexta). We show that hornworms from Arizona, Kansas, and North Carolina are genetically distinct, with Arizona being particularly differentiated from the other 2 populations using Illumina whole-genome resequencing. Peaks of differentiation exist across the genome, but here, we focus in on the most striking regions. In particular, we identify 2 likely segregating inversions found in the Arizona population. One inversion on the Z chromosome may enhance adaptive evolution of the sex chromosome. The larger, 8 Mb inversion on chromosome 12 contains a pseudogene which may be involved in the exploitation of a novel hostplant in Arizona, but functional genetic assays will be required to support this hypothesis. Nevertheless, our results reveal undiscovered natural variation and provide useful genomic data for both pest management and evolutionary genetics of this insect species.

Meiotic Chromosomal Abnormality Detected in a Heterozygote of Elymus nutans

Elymus nutans is an allopolyploid with a genome constitution of StStYYHH (2n = 6x = 42). Highly frequent intergenomic translocations and chromosomal variations with repeat amplification and deletions in E. nutans have been identified in the previous studies. However, more complicated structural variations such as chromosomal inversions or intra-genomic translocations are still unknown in this species, so does the reason for the origin of the chromosomal variations. Heterozygotes with rearranged chromosomes always present irregular meiosis behaviors, which subsequently cause the secondary chromosome rearrangements. Investigation on the meiosis of heterozygotes, especially on the individual chromosome level, may provide the important clues to identify the more complicated chromosome structural variations in the populations, and clarify the origin of the chromosome variations. In this study, meiotic analysis was conducted on a heterozygote plant of Elymus nutans, which showed high intra- and inter-genome chromosomal variations, by sequential fluorescence in situ hybridization (FISH) and genomic in situ hybridization (GISH), with each chromosome clearly recognized. The results showed chromosomal abnormalities at every meiotic stage and abnormalities in frequency variations between different sub-genomes and different individual chromosomes. The abnormalities were revealed as univalent, fragment, rod, or Y shape bivalent in diakinesis; univalent and rod bivalent in metaphase I; lagged and segregated chromatid, bridge, fragment of the sister chromatid, fragment, bridge accompanied with fragment, and unequal segregated chromosome in anaphase I; bridge and lagged chromatid in ana-telophase II; and micronucleus at uninucleate stage. Generally, the St and H genomes harbor more abnormalities than the Y genome. Moreover, a paracentric inversion in 2St was exclusively determined, and another paracentric inversion in 6Y was tentatively identified. In addition, novel deletions were clearly detected in 3H, 4H, 1Y, and 3Y homologous chromosomes; in particular, de novo pericentric inversion in 3H was repeatedly identified in metaphase I. The study revealed the chromosomal inversions pre-existed in parents or populations, as well as de novo inversions and deletions originated in the meiosis of the heterozygote in E. nutans. Moreover, it indicated wide range of meiosis abnormalities on different stages and different chromosomes, and suggests that secondary rearrangements contribute much to the chromosome variations in E. nutans.

Expression complementation of gene presence/absence polymorphisms in hybrids contributes importantly to heterosis in sunflower

INTRODUCTION

Numerous crops have transitioned to hybrid seed production to increase yields and yield stability through heterosis. However, the molecular mechanisms underlying heterosis and its stability across environments are not yet fully understood.

OBJECTIVES

This study aimed to (1) elucidate the genetic and molecular mechanisms underlying heterosis in sunflower, and (2) determine how heterosis is maintained under different environments.

METHODS

Genome-wide association (GWA) analyses were employed to assess the effects of presence/absence variants (PAVs) and stop codons on 16 traits phenotyped in the sunflower association mapping population at three locations. To link the GWA results to transcriptomic variation, we sequenced the transcriptomes of two sunflower cultivars and their F₁ hybrid (INEDI) under both control and drought conditions and analyzed patterns of gene expression and alternative splicing.

RESULTS

Thousands of PAVs were found to affect phenotypic variation using a relaxed significance threshold, and at most such loci the "absence" allele reduced values of heterotic traits, but not those of non-heterotic traits. This pattern was strengthened for PAVs that showed expression complementation in INEDI. Stop codons were much rarer than PAVs and less likely to reduce heterotic trait values. Hybrid expression patterns were enriched for the GO category, sensitivity to stimulus, but all genotypes responded to drought similarily - by up-regulating water stress response pathways and down-regulating metabolic pathways. Changes in alternative splicing were strongly negatively correlated with expression variation, implying that alternative splicing in this system largely acts to reinforce expression responses.

CONCLUSION

Our results imply that complementation of expression of PAVs in hybrids is a major contributor to heterosis in sunflower, consistent with the dominance model of heterosis. This mechanism can account for yield stability across different environments. Moreover, given the much larger numbers of PAVs in plant vs. animal genomes, it also offers an explanation for the stronger heterotic responses seen in the former.

Chromosome evolution and the genetic basis of agronomically important traits in greater yam

The nutrient-rich tubers of the greater yam, Dioscorea alata L., provide food and income security for millions of people around the world. Despite its global importance, however, greater yam remains an orphan crop. Here, we address this resource gap by presenting a highly contiguous chromosome-scale genome assembly of D. alata combined with a dense genetic map derived from African breeding populations. The genome sequence reveals an ancient allotetraploidization in the Dioscorea lineage, followed by extensive genome-wide reorganization. Using the genomic tools, we find quantitative trait loci for resistance to anthracnose, a damaging fungal pathogen of yam, and several tuber quality traits. Genomic analysis of breeding lines reveals both extensive inbreeding as well as regions of extensive heterozygosity that may represent interspecific introgression during domestication. These tools and insights will enable yam breeders to unlock the potential of this staple crop and take full advantage of its adaptability to varied environments.

Using structural variants to understand the ecological and evolutionary dynamics of fungal plant pathogens

Deletions, duplications, insertions, inversions and translocations are commonly referred to as structural variants (SVs). Fungal plant pathogens have compact genomes, facilitating the generation of accurate maps of SVs for these species in recent studies. Structural variants have been found to constitute a significant proportion of the standing genetic variation in fungal plant pathogen populations, potentially leading to the generation of accessory genes, regions or chromosomes enriched in pathogenicity factors. Structural variants are involved in the rapid adaptation and ecological traits of pathogens, including host specialization and mating. Long-read sequencing techniques coupled with theoretical and experimental approaches have considerable potential for elucidating the phenotypic effects of SVs and deciphering the evolutionary and genomic mechanisms underlying the formation of SVs in fungal plant pathogens.

Population genomic signatures of the oriental fruit moth related to the Pleistocene climates

The Quaternary climatic oscillations are expected to have had strong impacts on the evolution of species. Although legacies of the Quaternary climates on population processes have been widely identified in diverse groups of species, adaptive genetic changes shaped during the Quaternary have been harder to decipher. Here, we assembled a chromosome-level genome of the oriental fruit moth and compared genomic variation among refugial and colonized populations of this species that diverged in the Pleistocene. High genomic diversity was maintained in refugial populations. Demographic analysis showed that the effective population size of refugial populations declined during the penultimate glacial maximum (PGM) but remained stable during the last glacial maximum (LGM), indicating a strong impact of the PGM rather than the LGM on this pest species. Genome scans identified one chromosomal inversion and a mutation of the circadian gene Clk on the neo-Z chromosome potentially related to the endemicity of a refugial population. In the colonized populations, genes in pathways of energy metabolism and wing development showed signatures of selection. These different genomic signatures of refugial and colonized populations point to multiple impacts of Quaternary climates on adaptation in an extant species.

The oriental fruit moth is a pest species native to East Asia with refugial and colonized populations throughout the region. Here, a chromosome-level assembly for the species is reported and used to identify genomic signatures related to Quaternary climate change.

Wild sunflower goes viral: citizen science and comparative genomics allow tracking the origin and establishment of invasive sunflower in the Levant

Globalization and intensified volume of trade and transport around the world are accelerating the rate of biological invasions. It is therefore increasingly important to understand the processes through which invasive species colonize new habitats, often to the detriment of native flora. The initial steps of an invasion are particularly critical, as the introduced species relies on limited genetic diversity to adapt to a new environment. However, our understanding of this critical stage of the invasion is currently limited. We used a citizen science approach and social media to survey the distribution of invasive sunflower in Israel. We then sampled and sequenced a representative collection and compared it with available genomic data sets of North American wild sunflower, landraces and cultivars. We show that invasive wild sunflower is rapidly establishing throughout Israel, probably from a single, recent introduction from Texas, while maintaining high genetic diversity through ongoing gene flow. Since its introduction, invasive sunflower has spread quickly to most regions, and differentiation was detected despite extensive gene flow between clusters. Our findings suggest that rapid spread followed by continuous gene flow between diverging populations can serve as an efficient mechanism for maintaining sufficient genetic diversity at the early stages of invasion, promoting rapid adaptation and establishment in the new territory.

Performance of African-ancestry-specific polygenic hazard score varies according to local ancestry in 8q24

BACKGROUND

We previously developed an African-ancestry-specific polygenic hazard score (PHS46+African) that substantially improved prostate cancer risk stratification in men with African ancestry. The model consists of 46 SNPs identified in Europeans and 3 SNPs from 8q24 shown to improve model performance in Africans. Herein, we used principal component (PC) analysis to uncover subpopulations of men with African ancestry for whom the utility of PHS46+African may differ.

MATERIALS AND METHODS

Genotypic data were obtained from the PRACTICAL consortium for 6253 men with African genetic ancestry. Genetic variation in a window spanning 3 African-specific 8q24 SNPs was estimated using 93 PCs. A Cox proportional hazards framework was used to identify the pair of PCs most strongly associated with the performance of PHS46+African. A calibration factor (CF) was formulated using Cox coefficients to quantify the extent to which the performance of PHS46+African varies with PC.

RESULTS

CF of PHS46+African was strongly associated with the first and twentieth PCs. Predicted CF ranged from 0.41 to 2.94, suggesting that PHS46+African may be up to 7 times more beneficial to some African men than others. The explained relative risk for PHS46+African varied from 3.6% to 9.9% for individuals with low and high CF values, respectively. By cross-referencing our data set with 1000 Genomes, we identified significant associations between continental and calibration groupings.

CONCLUSION

We identified PCs within 8q24 that were strongly associated with the performance of PHS46+African. Further research to improve the clinical utility of polygenic risk scores (or models) is needed to improve health outcomes for men of African ancestry.

The Build-Up of Population Genetic Divergence along the Speciation Continuum during a Recent Adaptive Radiation of Rhagoletis Flies

New species form through the evolution of genetic barriers to gene flow between previously interbreeding populations. The understanding of how speciation proceeds is hampered by our inability to follow cases of incipient speciation through time. Comparative approaches examining different diverging taxa may offer limited inferences, unless they fulfill criteria that make the comparisons relevant. Here, we test for those criteria in a recent adaptive radiation of the Rhagoletis pomonella species group (RPSG) hypothesized to have diverged in sympatry via adaptation to different host fruits. We use a large-scale population genetic survey of 1568 flies across 33 populations to: (1) detect on-going hybridization, (2) determine whether the RPSG is derived from the same proximate ancestor, and (3) examine patterns of clustering and differentiation among sympatric populations. We find that divergence of each in-group RPSG taxon is occurring under current gene flow, that the derived members are nested within the large pool of genetic variation present in hawthorn-infesting populations of R. pomonella, and that sympatric population pairs differ markedly in their degree of genotypic clustering and differentiation across loci. We conclude that the RPSG provides a particularly robust opportunity to make direct comparisons to test hypotheses about how ecological speciation proceeds despite on-going gene flow.

Approximate Bayesian computation untangles signatures of contemporary and historical hybridization between two endangered species

Contemporary gene flow, when resumed after a period of isolation, can have crucial consequences for endangered species, as it can both increase the supply of adaptive alleles and erode local adaptation. Determining the history of gene flow and thus the importance of contemporary hybridization, however, is notoriously difficult. Here, we focus on two endangered plant species, Arabis nemorensis and A. sagittata, which hybridize naturally in a sympatric population located on the banks of the Rhine. Using reduced genome sequencing, we determined the phylogeography of the two taxa but report only a unique sympatric population. Molecular variation in chloroplast DNA indicated that A. sagittata is the principal receiver of gene flow. Applying classical D-statistics and its derivatives to whole-genome data of 35 accessions, we detect gene flow not only in the sympatric population but also among allopatric populations. Using an Approximate Bayesian computation approach, we identify the model that best describes the history of gene flow between these taxa. This model shows that low levels of gene flow have persisted long after speciation. Around 10 000 years ago, gene flow stopped and a period of complete isolation began. Eventually, a hotspot of contemporary hybridization was formed in the unique sympatric population. Occasional sympatry may have helped protect these lineages from extinction in spite of their extremely low diversity.

Genomic population structure and local adaptation of the wild strawberry Fragaria nilgerrensis

The crop wild relative, Fragaria nilgerrensis, is adapted to a variety of diverse habitats across its native range in China. Thus, discoveries made in this species could serve useful to guide the development of new superior strawberry cultivars that are resilient to new or variable environments. However, the genetic diversity and genetic architecture of traits in this species underlying important adaptive traits remain poorly understood. Here, we used whole-genome resequencing data from 193 F. nilgerrensis individuals spanning the distribution range in China to investigate the genetic diversity, population structure and the genomic basis of local adaptation. We identified four genetic groups, with the western group located in Hengduan Mountains exhibited the highest genetic diversity. Redundancy analysis suggests that both environment and geographic variables shaped a significant proportion of genomic variation. Our analyses revealed that the environmental difference explains more of the observed genetic variation than geographic distance. This suggests that adaptation to distinct habitats, unique combination of abiotic factors, likely drove genetic differentiation. Lastly, by implementing selective sweeps scans and genome-environment association analysis throughout the genome, we identified the genetic variation associated with local adaptation and investigated the functions of putative candidate genes in F. nilgerrensis.

Genetic basis and dual adaptive role of floral pigmentation in sunflowers

Variation in floral displays, both between and within species, has been long known to be shaped by the mutualistic interactions that plants establish with their pollinators. However, increasing evidence suggests that abiotic selection pressures influence floral diversity as well. Here, we analyse the genetic and environmental factors that underlie patterns of floral pigmentation in wild sunflowers. While sunflower inflorescences appear invariably yellow to the human eye, they display extreme diversity for patterns of ultraviolet pigmentation, which are visible to most pollinators. We show that this diversity is largely controlled by cis-regulatory variation affecting a single MYB transcription factor, HaMYB111, through accumulation of ultraviolet (UV)-absorbing flavonol glycosides in ligules (the ‘petals’ of sunflower inflorescences). Different patterns of ultraviolet pigments in flowers are strongly correlated with pollinator preferences. Furthermore, variation for floral ultraviolet patterns is associated with environmental variables, especially relative humidity, across populations of wild sunflowers. Ligules with larger ultraviolet patterns, which are found in drier environments, show increased resistance to desiccation, suggesting a role in reducing water loss. The dual role of floral UV patterns in pollinator attraction and abiotic response reveals the complex adaptive balance underlying the evolution of floral traits.

Flowers are an important part of how many plants reproduce. Their distinctive colours, shapes and patterns attract specific pollinators, but they can also help to protect the plant from predators and environmental stresses.

Many flowers contain pigments that absorb ultraviolet (UV) light to display distinct UV patterns – although invisible to the human eye, most pollinators are able to see them. For example, when seen in UV, sunflowers feature a ‘bullseye’ with a dark centre surrounded by a reflective outer ring. The sizes and thicknesses of these rings vary a lot within and between flower species, and so far, it has been unclear what causes this variation and how it affects the plants.

To find out more, Todesco et al. studied the UV patterns in various wild sunflowers across North America by considering the ecology and molecular biology of different plants. This revealed great variation between the UV patterns of the different sunflower populations. Moreover, Todesco et al. found that a gene called HaMYB111 is responsible for the diverse UV patterns in the sunflowers. This gene controls how plants make chemicals called flavonols that absorb UV light.

Flavonols also help to protect plants from damage caused by droughts and extreme temperatures. Todesco et al. showed that plants with larger bullseyes had more flavonols, attracted more pollinators, and were better at conserving water. Accordingly, these plants were found in drier locations.

This study suggests that, at least in sunflowers, UV patterns help both to attract pollinators and to control water loss. These insights could help to improve pollination – and consequently yield – in cultivated plants, and to develop plants with better resistance to extreme weather. This work also highlights the importance of combining biology on small and large scales to understand complex processes, such as adaptation and evolution.

Signatures of natural selection in a foundation tree along Mediterranean climatic gradients

Temperature and precipitation regimes are rapidly changing, resulting in forest dieback and extinction events, particularly in Mediterranean‐type climates (MTC). Forest management that enhance forests’ resilience is urgently required, however adaptation to climates in heterogeneous landscapes with multiple selection pressures is complex. For widespread trees in MTC we hypothesized that: patterns of local adaptation are associated with climate; precipitation is a stronger factor of adaptation than temperature; functionally related genes show similar signatures of adaptation; and adaptive variants are independently sorting across the landscape. We sampled 28 populations across the geographic distribution of Eucalyptus marginata (jarrah), in South‐west Western Australia, and obtained 13,534 independent single nucleotide polymorphic (SNP) markers across the genome. Three genotype‐association analyses that employ different ways of correcting population structure were used to identify putatively adapted SNPs associated with independent climate variables. While overall levels of population differentiation were low (F_ST = 0.04), environmental association analyses found a total of 2336 unique SNPs associated with temperature and precipitation variables, with 1440 SNPs annotated to genic regions. Considerable allelic turnover was identified for SNPs associated with temperature seasonality and mean precipitation of the warmest quarter, suggesting that both temperature and precipitation are important factors in adaptation. SNPs with similar gene functions had analogous allelic turnover along climate gradients, while SNPs among temperature and precipitation variables had uncorrelated patterns of adaptation. These contrasting patterns provide evidence that there may be standing genomic variation adapted to current climate gradients, providing the basis for adaptive management strategies to bolster forest resilience in the future.

Predictability and parallelism in the contemporary evolution of hybrid genomes

Hybridization between species is widespread across the tree of life. As a result, many species, including our own, harbor regions of their genome derived from hybridization. Despite the recognition that this process is widespread, we understand little about how the genome stabilizes following hybridization, and whether the mechanisms driving this stabilization tend to be shared across species. Here, we dissect the drivers of variation in local ancestry across the genome in replicated hybridization events between two species pairs of swordtail fish: Xiphophorus birchmanni × X. cortezi and X. birchmanni × X. malinche. We find unexpectedly high levels of repeatability in local ancestry across the two types of hybrid populations. This repeatability is attributable in part to the fact that the recombination landscape and locations of functionally important elements play a major role in driving variation in local ancestry in both types of hybrid populations. Beyond these broad scale patterns, we identify dozens of regions of the genome where minor parent ancestry is unusually low or high across species pairs. Analysis of these regions points to shared sites under selection across species pairs, and in some cases, shared mechanisms of selection. We show that one such region is a previously unknown hybrid incompatibility that is shared across X. birchmanni × X. cortezi and X. birchmanni × X. malinche hybrid populations.

Wheat speciation and adaptation: perspectives from reticulate evolution

Reticulate evolution through the interchanging of genetic components across organisms can impact significantly on the fitness and adaptation of species. Bread wheat (Triticum aestivum subsp. aestivum) is one of the most important crops in the world. Allopolyploid speciation, frequent hybridization, extensive introgression, and occasional horizontal gene transfer (HGT) have been shaping a typical paradigm of reticulate evolution in bread wheat and its wild relatives, which is likely to have a substantial influence on phenotypic traits and environmental adaptability of bread wheat. In this review, we outlined the evolutionary history of bread wheat and its wild relatives with a highlight on the interspecific hybridization events, demonstrating the reticulate relationship between species/subspecies in the genera Triticum and Aegilops. Furthermore, we discussed the genetic mechanisms and evolutionary significance underlying the introgression of bread wheat and its wild relatives. An in-depth understanding of the evolutionary process of Triticum species should be beneficial to future genetic study and breeding of bread wheat.

An accurate assignment test for extremely low-coverage whole-genome sequence data

Genomic assignment tests can provide important diagnostic biological characteristics, such as population of origin or ecotype. Yet, assignment tests often rely on moderate- to high-coverage sequence data that can be difficult to obtain for fields such as molecular ecology and ancient DNA. We have developed a novel approach that efficiently assigns biologically relevant information (i.e., population identity or structural variants such as inversions) in extremely low-coverage sequence data. First, we generate databases from existing reference data using a subset of diagnostic single nucleotide polymorphisms (SNPs) associated with a biological characteristic. Low-coverage alignment files are subsequently compared to these databases to ascertain allelic state, yielding a joint probability for each association. To assess the efficacy of this approach, we assigned haplotypes and population identity in Heliconius butterflies, Atlantic herring, and Atlantic cod using chromosomal inversion sites and whole-genome data. We scored both modern and ancient specimens, including the first whole-genome sequence data recovered from ancient Atlantic herring bones. The method accurately assigns biological characteristics, including population membership, using extremely low-coverage data (as low as 0.0001x) based on genome-wide SNPs. This approach will therefore increase the number of samples in evolutionary, ecological and archaeological research for which relevant biological information can be obtained.

Origins and evolution of extreme life span in Pacific Ocean rockfishes

Pacific Ocean rockfishes (genus Sebastes) exhibit extreme variation in life span, with some species being among the most long-lived extant vertebrates. We de novo assembled the genomes of 88 rockfish species and from these identified repeated signatures of positive selection in DNA repair pathways in long-lived taxa and 137 longevity-associated genes with direct effects on life span through insulin signaling and with pleiotropic effects through size and environmental adaptations. A genome-wide screen of structural variation reveals copy number expansions in the immune modulatory butyrophilin gene family in long-lived species. The evolution of different rockfish life histories is coupled to genetic diversity and reshapes the mutational spectrum driving segregating CpG→TpG variants in long-lived species. These analyses highlight the genetic innovations that underlie life history trait adaptations and, in turn, how they shape genomic diversity.

Current status of structural variation studies in plants

Structural variations (SVs) including gene presence/absence variations and copy number variations are a common feature of genomes in plants and, together with single nucleotide polymorphisms and epigenetic differences, are responsible for the heritable phenotypic diversity observed within and between species. Understanding the contribution of SVs to plant phenotypic variation is important for plant breeders to assist in producing improved varieties. The low resolution of early genetic technologies and inefficient methods have previously limited our understanding of SVs in plants. However, with the rapid expansion in genomic technologies, it is possible to assess SVs with an ever-greater resolution and accuracy. Here, we review the current status of SV studies in plants, examine the roles that SVs play in phenotypic traits, compare current technologies and assess future challenges for SV studies.

Genetic control of arbuscular mycorrhizal colonization by Rhizophagus intraradices in Helianthus annuus (L.)

Plant symbiosis with arbuscular mycorrhizal (AM) fungi provides many benefits, including increased nutrient uptake, drought tolerance, and belowground pathogen resistance. To develop a better understanding of the genetic architecture of mycorrhizal symbiosis, we conducted a genome-wide association study (GWAS) of this plant-fungal interaction in cultivated sunflower. A diversity panel of cultivated sunflower (Helianthus annuus L.) was phenotyped for root colonization under inoculation with the AM fungus Rhizophagus intraradices. Using a mixed linear model approach with a high-density genetic map, we identified genomic regions that are likely associated with R. intraradices colonization in sunflower. Additionally, we used a set of twelve diverse lines to assess the effect that inoculation with R. intraradices has on dried shoot biomass and macronutrient uptake. Colonization among lines in the mapping panel ranged from 0-70% and was not correlated with mycorrhizal growth response, shoot phosphorus response, or shoot potassium response among the Core 12 lines. Association mapping yielded three single-nucleotide polymorphisms (SNPs) that were significantly associated with R. intraradices colonization. This is the first study to use GWAS to identify genomic regions associated with AM colonization in an Asterid eudicot species. Three genes of interest identified from the regions containing these SNPs are likely related to plant defense.

Biodiversity, resilience and the stability of evolutionary systems

Various macro-evolutionary phenomena, such as long-term stability punctuated by bursts of evolution, are difficult to explain via the micro-evolutionary process of weak selection acting steadily on individual mutations. In contrast, bursts of change are expected if evolutionary systems are complex and balanced, with occasional disruption of balance. Such disruption represents the collapse of resilience, akin to the snapping of an elastic band. It can be driven by external factors, or by self-propagating feedback loops internal to a system. Thus, evolutionary resilience could help explain how evolution generates broader patterns of biodiversity. We outline evidence and tests for this hypothesis, which emphasizes the processes balancing evolution, as urged fifty years ago in ecological genetics and via modern results in a range of systems.

Allen Orr and the genetics of adaptation

Over most of the 20th century, evolutionary biologists predominantly subscribed to a strong form of "micro-mutationism," in which adaptive phenotypic divergence arises from allele frequency changes at many loci, each with a small effect on the phenotype. To be sure, there were well-known examples of large-effect alleles contributing to adaptation, yet such cases were generally regarded as atypical and unrepresentative of evolutionary change in general. In 1998, Allen Orr published a landmark theoretical paper in Evolution, which showed that both small- and large-effect mutations are likely to contribute to "adaptive walks" of a population to an optimum. Coupled with a growing set of empirical examples of large-effect alleles contributing to divergence (e.g., from QTL studies), Orr's paper provided a mathematical formalism that converted many evolutionary biologists from micro-mutationism to a more pluralistic perspective on the genetic basis of evolutionary change. We revisit the theoretical insights emerging from Orr's paper within the historical context leading up to 1998, and track the influence of this paper on the field of evolutionary biology through an examination of its citations over the last two decades and an analysis of the extensive body of theoretical and empirical research that Orr's pioneering paper inspired.

Genomic signatures of parallel alpine adaptation in recently evolved flightless insects

Natural selection along elevational gradients has potential to drive predictable adaptations across distinct lineages, but the extent of such repeated evolution remains poorly studied for many widespread alpine taxa. We present parallel genomic analyses of two recently evolved flightless alpine insect lineages to test for molecular signatures of repeated alpine adaptation. Specifically, we compare low-elevation vs. alpine stonefly ecotypes from parallel stream populations in which flightless upland ecotypes have been independently derived. We map 67,922 polymorphic genetic markers, generated across 176 Zelandoperla fenestrata specimens from two independent alpine stream populations in New Zealand's Rock and Pillar Range, to a newly developed plecopteran reference genome. Genome-wide scans revealed 31 regions with outlier single nucleotide polymorphisms (SNPs) differentiating lowland vs. alpine ecotypes in Lug Creek, and 37 regions with outliers differentiating ecotypes in Six Mile Creek. Of these regions, 13% (8/60) yielded outlier SNPs across both within-stream ecotype comparisons, implying comparable genomic shifts contribute to this repeated alpine adaptation. Candidate genes closely linked to repeated outlier regions include several with documented roles in insect wing-development (e.g., dishevelled), suggesting that they may contribute to repeated alpine wing reduction. Additional candidate genes have been shown to influence insect fecundity (e.g., ovo) and lifespan (e.g., Mrp4), implying that they might contribute to life history differentiation between upland and lowland ecotypes. Additional outlier genes have potential roles in the evolution of reproductive isolation among ecotypes (hedgehog and Desaturase 1). These results demonstrate how replicated outlier tests across independent lineages can potentially contribute to the discovery of genes underpinning repeated adaptation.

The evolutionary genomics of species' responses to climate change

Climate change is a threat to biodiversity. One way that this threat manifests is through pronounced shifts in the geographical range of species over time. To predict these shifts, researchers have primarily used species distribution models. However, these models are based on assumptions of niche conservatism and do not consider evolutionary processes, potentially limiting their accuracy and value. To incorporate evolution into the prediction of species' responses to climate change, researchers have turned to landscape genomic data and examined information about local genetic adaptation using climate models. Although this is an important advancement, this approach currently does not include other evolutionary processes-such as gene flow, population dispersal and genomic load-that are critical for predicting the fate of species across the landscape. Here, we briefly review the current practices for the use of species distribution models and for incorporating local adaptation. We next discuss the rationale and theory for considering additional processes, reviewing how they can be incorporated into studies of species' responses to climate change. We summarize with a conceptual framework of how manifold layers of information can be combined to predict the potential response of specific populations to climate change. We illustrate all of the topics using an exemplar dataset and provide the source code as potential tutorials. This Perspective is intended to be a step towards a more comprehensive integration of population genomics with climate change science.

Aberrant RNA splicing due to genetic incompatibilities in sunflower hybrids

Genome-scale studies have revealed divergent mRNA splicing patterns between closely related species or populations. However, it is unclear whether splicing differentiation is a simple byproduct of population divergence, or whether it also acts as a mechanism for reproductive isolation. We examined mRNA splicing in wild × domesticated sunflower hybrids and observed 45 novel splice forms that were not found in the wild or domesticated parents, in addition to 16 high-expression parental splice forms that were absent in one or more hybrids. We identify loci associated with variation in the levels of these splice forms, finding that many aberrant transcripts were regulated by multiple alleles with nonadditive interactions. We identified particular spliceosome components that were associated with 21 aberrant isoforms, more than half of which were located in or near regulatory QTL. These incompatibilities often resulted in alteration in the protein-coding regions of the novel transcripts in the form of frameshifts and truncations. By associating the splice variation in these genes with size and growth rate measurements, we found that the cumulative expression of all aberrant transcripts was correlated with a significant reduction in growth rate. Our results lead us to propose a model where divergent splicing regulatory loci could act as incompatibility loci that contribute to the evolution of reproductive isolation.

Genetic trade-offs underlie divergent life history strategies for local adaptation in white clover

Local adaptation is common in plants, yet characterization of its underlying genetic basis is rare in herbaceous perennials. Moreover, while many plant species exhibit intraspecific chemical defence polymorphisms, their importance for local adaptation remains poorly understood. We examined the genetic architecture of local adaptation in a perennial, obligately-outcrossing herbaceous legume, white clover (Trifolium repens). This widespread species displays a well-studied chemical defence polymorphism for cyanogenesis (HCN release following tissue damage) and has evolved climate-associated cyanogenesis clines throughout its range. Two biparental F₂  mapping populations, derived from three parents collected in environments spanning the U.S. latitudinal species range (Duluth, MN, St. Louis, MO and Gainesville, FL), were grown in triplicate for two years in reciprocal common garden experiments in the parental environments (6,012 total plants). Vegetative growth and reproductive fitness traits displayed trade-offs across reciprocal environments, indicating local adaptation. Genetic mapping of fitness traits revealed a genetic architecture characterized by allelic trade-offs between environments, with 100% and 80% of fitness QTL in the two mapping populations showing significant QTL×E interactions, consistent with antagonistic pleiotropy. Across the genome there were three hotspots of QTL colocalization. Unexpectedly, we found little evidence that the cyanogenesis polymorphism contributes to local adaptation. Instead, divergent life history strategies in reciprocal environments were major fitness determinants: selection favoured early investment in flowering at the cost of multiyear survival in the southernmost site versus delayed flowering and multiyear persistence in the northern environments. Our findings demonstrate that multilocus genetic trade-offs contribute to contrasting life history characteristics that allow for local adaptation in this outcrossing herbaceous perennial.

Introgression among North American wild grapes (Vitis) fuels biotic and abiotic adaptation

BACKGROUND

Introgressive hybridization can reassort genetic variants into beneficial combinations, permitting adaptation to new ecological niches. To evaluate evolutionary patterns and dynamics that contribute to introgression, we investigate six wild Vitis species that are native to the Southwestern United States and useful for breeding grapevine (V. vinifera) rootstocks.

RESULTS

By creating a reference genome assembly from one wild species, V. arizonica, and by resequencing 130 accessions, we focus on identifying putatively introgressed regions (pIRs) between species. We find six species pairs with signals of introgression between them, comprising up to ~ 8% of the extant genome for some pairs. The pIRs tend to be gene poor, located in regions of high recombination and enriched for genes implicated in disease resistance functions. To assess potential pIR function, we explore SNP associations to bioclimatic variables and to bacterial levels after infection with the causative agent of Pierce's disease (Xylella fastidiosa). pIRs are enriched for SNPs associated with both climate and bacterial levels, suggesting that introgression is driven by adaptation to biotic and abiotic stressors.

CONCLUSIONS

Altogether, this study yields insights into the genomic extent of introgression, potential pressures that shape adaptive introgression, and the evolutionary history of economically important wild relatives of a critical crop.

Genomic structural variants constrain and facilitate adaptation in natural populations of Theobroma cacao, the chocolate tree

Genomic structural variants (SVs) can play important roles in adaptation and speciation. Yet the overall fitness effects of SVs are poorly understood, partly because accurate population-level identification of SVs requires multiple high-quality genome assemblies. Here, we use 31 chromosome-scale, haplotype-resolved genome assemblies of Theobroma cacao-an outcrossing, long-lived tree species that is the source of chocolate-to investigate the fitness consequences of SVs in natural populations. Among the 31 accessions, we find over 160,000 SVs, which together cover eight times more of the genome than single-nucleotide polymorphisms and short indels (125 versus 15 Mb). Our results indicate that a vast majority of these SVs are deleterious: they segregate at low frequencies and are depleted from functional regions of the genome. We show that SVs influence gene expression, which likely impairs gene function and contributes to the detrimental effects of SVs. We also provide empirical support for a theoretical prediction that SVs, particularly inversions, increase genetic load through the accumulation of deleterious nucleotide variants as a result of suppressed recombination. Despite the overall detrimental effects, we identify individual SVs bearing signatures of local adaptation, several of which are associated with genes differentially expressed between populations. Genes involved in pathogen resistance are strongly enriched among these candidates, highlighting the contribution of SVs to this important local adaptation trait. Beyond revealing empirical evidence for the evolutionary importance of SVs, these 31 de novo assemblies provide a valuable resource for genetic and breeding studies in T cacao.

Threespine Stickleback: A Model System For Evolutionary Genomics

The repeated adaptation of oceanic threespine sticklebacks to fresh water has made it a premier organism to study parallel evolution. These small fish have multiple distinct ecotypes that display a wide range of diverse phenotypic traits. Ecotypes are easily crossed in the laboratory, and families are large and develop quickly enough for quantitative trait locus analyses, positioning the threespine stickleback as a versatile model organism to address a wide range of biological questions. Extensive genomic resources, including linkage maps, a high-quality reference genome, and developmental genetics tools have led to insights into the genomic basis of adaptation and the identification of genomic changes controlling traits in vertebrates. Recently, threespine sticklebacks have been used as a model system to identify the genomic basis of highly complex traits, such as behavior and host-microbiome and host-parasite interactions. We review the latest findings and new avenues of research that have led the threespine stickleback to be considered a supermodel of evolutionary genomics.

The Insights into Mitochondrial Genomes of Sunflowers

The significant difference in the mtDNA size and structure with simultaneous slow evolving genes makes the mitochondrial genome paradoxical among all three DNA carriers in the plant cell. Such features make mitochondrial genome investigations of particular interest. The genus Helianthus is a diverse taxonomic group, including at least two economically valuable species-common sunflower (H. annuus) and Jerusalem artichoke (H. tuberosus). The successful investigation of the sunflower nuclear genome provided insights into some genomics aspects and significantly intensified sunflower genetic studies. However, the investigations of organelles' genetic information in Helianthus, especially devoted to mitochondrial genomics, are presented by limited studies. Using NGS sequencing, we assembled the complete mitochondrial genomes for H. occidentalis (281,175 bp) and H. tuberosus (281,287 bp) in the current investigation. Besides the master circle chromosome, in the case of H. tuberosus, the 1361 bp circular plasmid was identified. The mitochondrial gene content was found to be identical for both sunflower species, counting 32 protein-coding genes, 3 rRNA, 23 tRNA genes, and 18 ORFs. The comparative analysis between perennial sunflowers revealed common and polymorphic SSR and SNPs. Comparison of perennial sunflowers with H. annuus allowed us to establish similar rearrangements in mitogenomes, which have possibly been inherited from a common ancestor after the divergence of annual and perennial sunflower species. It is notable that H. occidentalis and H. tuberosus mitogenomes are much more similar to H. strumosus than H. grosseserratus.

Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.

Expanding the conservation genomics toolbox: Incorporating structural variants to enhance genomic studies for species of conservation concern

Structural variants (SVs) are large rearrangements (>50 bp) within the genome that impact gene function and the content and structure of chromosomes. As a result, SVs are a significant source of functional genomic variation, that is, variation at genomic regions underpinning phenotype differences, that can have large effects on individual and population fitness. While there are increasing opportunities to investigate functional genomic variation in threatened species via single nucleotide polymorphism (SNP) data sets, SVs remain understudied despite their potential influence on fitness traits of conservation interest. In this future-focused Opinion, we contend that characterizing SVs offers the conservation genomics community an exciting opportunity to complement SNP-based approaches to enhance species recovery. We also leverage the existing literature-predominantly in human health, agriculture and ecoevolutionary biology-to identify approaches for readily characterizing SVs and consider how integrating these into the conservation genomics toolbox may transform the way we manage some of the world's most threatened species.

Arabidopsis-Based Dual-Layered Biological Network Analysis Elucidates Fully Modulated Pathways Related to Sugarcane Resistance on Biotrophic Pathogen Infection

We assembled a dual-layered biological network to study the roles of resistance gene analogs (RGAs) in the resistance of sugarcane to infection by the biotrophic fungus causing smut disease. Based on sugarcane-Arabidopsis orthology, the modeling used metabolic and protein-protein interaction (PPI) data from Arabidopsis thaliana (from Kyoto Encyclopedia of Genes and Genomes (KEGG) and BioGRID databases) and plant resistance curated knowledge for Viridiplantae obtained through text mining of the UniProt/SwissProt database. With the network, we integrated functional annotations and transcriptome data from two sugarcane genotypes that differ significantly in resistance to smut and applied a series of analyses to compare the transcriptomes and understand both signal perception and transduction in plant resistance. We show that the smut-resistant sugarcane has a larger arsenal of RGAs encompassing transcriptionally modulated subnetworks with other resistance elements, reaching hub proteins of primary metabolism. This approach may benefit molecular breeders in search of markers associated with quantitative resistance to diseases in non-model systems.

Complex evolution of novel red floral color in Petunia

Red flower color has arisen multiple times and is generally associated with hummingbird pollination. The majority of evolutionary transitions to red color proceeded from purple lineages and tend to be genetically simple, almost always involving a few loss-of-function mutations of major phenotypic effect. Here we report on the complex evolution of a novel red floral color in the hummingbird-pollinated Petunia exserta (Solanaceae) from a colorless ancestor. The presence of a red color is remarkable because the genus cannot synthesize red anthocyanins and P. exserta retains a nonfunctional copy of the key MYB transcription factor AN2. We show that moderate upregulation and a shift in tissue specificity of an AN2 paralog, DEEP PURPLE, restores anthocyanin biosynthesis in P. exserta. An essential shift in anthocyanin hydroxylation occurred through rebalancing the expression of three hydroxylating genes. Furthermore, the downregulation of an acyltransferase promotes reddish hues in typically purple pigments by preventing acyl group decoration of anthocyanins. This study presents a rare case of a genetically complex evolutionary transition toward the gain of a novel red color.