Massive haplotypes underlie ecotypic differentiation in sunflowers

'Chimes of resilience': what makes forest trees genetically resilient?

Forest trees are foundation species of many ecosystems and are challenged by global environmental changes. We assemble genetic facts and arguments supporting or undermining resilient responses of forest trees to those changes. Genetic resilience is understood here as the capacity of a species to restore its adaptive potential following environmental changes and disturbances. Importantly, the data come primarily from European temperate tree species with large distributions and consider only marginally species with small distributions. We first examine historical trajectories of trees during repeated climatic changes. Species that survived the Pliocene-Pleistocene transition and underwent the oscillations of glacial and interglacial periods were equipped with life history traits enhancing persistence and resilience. Evidence of their resilience also comes from the maintenance of large effective population sizes across time and rapid microevolutionary responses to recent climatic events. We then review genetic mechanisms and attributes shaping resilient responses. Usually, invoked constraints to resilience, such as genetic load or generation time and overlap, have limited consequences or are offset by positive impacts. Conversely, genetic plasticity, gene flow, introgression, genetic architecture of fitness-related traits and demographic dynamics strengthen resilience by accelerating adaptive responses. Finally, we address the limitations of this review and highlight critical research gaps.

Beyond the Standard GWAS-A Guide for Plant Biologists

Classic genome-wide association studies (GWAS) look for associations between individual single-nucleotide polymorphisms (SNPs) and phenotypes of interest. With the rapid progress of high-throughput genotyping and phenotyping technologies, GWAS have become increasingly powerful for detecting genetic determinants and their molecular mechanisms underpinning natural phenotypic variation. However, GWAS frequently yield results with neither expected nor promising loci, nor any significant associations. This is often because associations between SNPs and a single phenotype are confounded, for example with the environment, other traits or complex genetic structures. Such confounding can mask true genotype-phenotype associations, or inflate spurious associations. To address these problems, numerous methods have been developed that go beyond the standard model. Such advanced GWAS models are flexible and can offer improved statistical power for understanding the genetics underlying complex traits. Despite this advantage, these models have not been widely adopted and implemented compared to the standard GWAS approach, partly because this literature is diverse and often technical. In this review, our aim is to provide an overview of the application and the benefits of various advanced GWAS models for handling complex traits and genetic structures, targeting plant biologists who wish to carry out GWAS more effectively.

Haplotype-resolved genomes provide insights into the origins and functional significance of genome diversity in bivalves

Bivalves are famed for exhibiting vast genetic diversity of poorly understood origins and functional significance. Through comparative genomics, we demonstrate that high genetic diversity in these invertebrates is not directly linked to genome size. Using oysters as a representative clade, we show that despite genome size reduction during evolution, these bivalves maintain remarkable genetic variability. By constructing a haplotype-resolved genome for Crassostrea sikamea, we identify widespread haplotype divergent sequences (HDSs), representing genomic regions unique to each haplotype. We show that HDSs are driven by transposable elements, playing a key role in creating and maintaining genetic diversity during oyster evolution. Comparisons of haplotype-resolved genomes across four bivalve orders uncover diverse HDS origins, highlighting a role in genetic innovation and expression regulation across broad timescales. Further analyses show that, in oysters, haplotype polymorphisms drive gene expression variation, which is likely to promote phenotypic plasticity and adaptation. These findings advance our understanding of the relationships among genome structure, diversity, and adaptability in a highly successful invertebrate group.

Adaptation repeatedly uses complex structural genomic variation

Structural elements are widespread across genomes, but their complexity and role in repeatedly driving local adaptation remain unclear. In this work, we use phased genome assemblies to show that adaptive divergence in cryptic color pattern in a stick insect is repeatedly underlain by structural variation, but not a simple chromosomal inversion. We found that color pattern in populations of stick insects on two mountains is associated with translocations that have also been inverted. These translocations differ in size and origin on each mountain, but they overlap partially and involve some of the same gene regions. Moreover, this structural variation is subject to divergent selection and arose without introgression between species. Our results show how the origin of structural variation provides a mechanism for repeated bouts of adaptation.

Uncovering the genomic basis of phenological traits in Chouardia litardierei (Asparagaceae) through a genome-wide association study (GWAS)

Chouardia litardierei (Asparagaceae) is a non-model, perennial species characterized by exceptional ecological plasticity. In this research, we studied the genetic architecture underlying several phenological traits in selected ecologically diverged populations of this species. We conducted a genome-wide association study (GWAS) to identify genomic regions linked to the following populations-specific phenological traits: Beginning of Sprouting (BOS), Beginning of Flowering (BOF), Flowering Period Duration (FPD), and Vegetation Period Duration (VPD). Combining phenological data from a common garden experiment with an SNP dataset obtained through the ddRAD-seq approach, we identified numerous loci associated with these traits using single- and multi-locus GWAS models. Narrow-sense heritability estimates were high for all traits, with the VPD trait showing the highest estimate (86.95%), emphasizing its importance for local adaptation. Functional annotation of associated genomic regions revealed key protein families involved in flowering time regulation, vegetative growth timing, and stress adaptation. These findings provide insights into the molecular mechanisms of local adaptation in C. litardierei's populations from different habitats, emphasizing the role of genetic factors in phenological trait variation and ecological divergence across populations.

A Segregating Structural Variant Defines Novel Venom Phenotypes in the Eastern Diamondback Rattlesnake

Of all mutational mechanisms contributing to phenotypic variation, structural variants are both among the most capable of causing major effects as well as the most technically challenging to identify. Intraspecific variation in snake venoms is widely reported, and one of the most dramatic patterns described is the parallel evolution of streamlined neurotoxic rattlesnake venoms from hemorrhagic ancestors by means of deletion of snake venom metalloproteinase (SVMP) toxins and recruitment of neurotoxic dimeric phospholipase A2 (PLA2) toxins. While generating a haplotype-resolved, chromosome-level genome assembly for the eastern diamondback rattlesnake (Crotalus adamanteus), we discovered that our genome animal was heterozygous for a ∼225 Kb deletion containing six SVMP genes, paralleling one of the two steps involved in the origin of neurotoxic rattlesnake venoms. Range-wide population-genomic analysis revealed that, although this deletion is rare overall, it is the dominant homozygous genotype near the northwestern periphery of the species' range, where this species is vulnerable to extirpation. Although major SVMP deletions have been described in at least five other rattlesnake species, C. adamanteus is unique in not additionally gaining neurotoxic PLA2s. Previous work established a superficially complementary north-south gradient in myotoxin (MYO) expression based on copy number variation with high expression in the north and low in the south, yet we found that the SVMP and MYO genotypes vary independently, giving rise to an array of diverse, novel venom phenotypes across the range. Structural variation, therefore, forms the basis for the major axes of geographic venom variation for C. adamanteus.

The long and short of hyperdivergent regions

The increasing prevalence of genome sequencing and assembly has uncovered evidence of hyperdivergent genomic regions - loci with excess genetic diversity - in species across the tree of life. Hyperdivergent regions are often enriched for genes that mediate environmental responses, such as immunity, parasitism, and sensory perception. Especially in self-fertilizing species where the majority of the genome is homozygous, the existence of hyperdivergent regions might imply the historical action of evolutionary forces such as introgression and/or balancing selection. We anticipate that the application of new sequencing technologies, broader taxonomic sampling, and evolutionary modeling of hyperdivergent regions will provide insights into the mechanisms that generate and maintain genetic diversity within and between species.

Copy number variation contributes to parallel local adaptation in an invasive plant

Adaptation is a critical determinant of the diversification, persistence, and geographic range limits of species. Yet the genetic basis of adaptation is often unknown and potentially underpinned by a wide range of mutational types-from single nucleotide changes to large-scale alterations of chromosome structure. Copy number variation (CNV) is thought to be an important source of adaptive genetic variation, as indicated by decades of candidate gene studies that point to CNVs underlying rapid adaptation to strong selective pressures. Nevertheless, population-genomic studies of CNVs face unique logistical challenges not encountered by other forms of genetic variation. Consequently, few studies have systematically investigated the contributions of CNVs to adaptation at a genome-wide scale. We present a genome-wide analysis of CNV contributing to the adaptation of an invasive weed, Ambrosia artemisiifolia. CNVs show clear signatures of parallel local adaptation between North American (native) and European (invaded) ranges, implying widespread reuse of CNVs during adaptation to shared heterogeneous patterns of selection. We used a local principal component analysis (PCA) to genotype CNV regions in whole-genome sequences of samples collected over the last two centuries. We identified 16 large CNV regions of up to 11.85 megabases in length, eight of which show signals of rapid evolutionary change, with pronounced frequency shifts between historic and modern populations. Our results provide compelling genome-wide evidence that CNV underlies rapid adaptation over contemporary timescales of natural populations.

Repeatable Selection on Large Ancestry Blocks in an Avian Hybrid Zone

Hybrid zones create natural tests of genetic incompatibilities by combining loci from 2 species in the same genetic background in the wild, making them useful for identifying loci involved in both intrinsic and ecological (extrinsic) isolation. Two Swainson's thrush subspecies form a hybrid zone in western North America. These coastal and inland subspecies exhibit dramatic differences in migration routes; their hybrids exhibit poor migratory survival, suggesting that ecological incompatibilities maintain this zone. We used a panel of ancestry informative markers to identify repeated patterns of selection and introgression across 4 hybrid populations that span the entire length of the Swainson's thrush hybrid zone. Two repeatable patterns consistent with selection against incompatibilities-steep genomic clines and few transitions between ancestry states-were found in large genetic blocks on chromosomes 1 and 5. The block on chromosome 1 showed evidence for inland subspecies introgression while the block on chromosome 5 exhibited coastal subspecies introgression. Some regions previously associated with migratory phenotypes, including migratory orientation, or exhibiting misexpression between the subspecies exhibited signatures of selection in the hybrid zone. Both selection and introgression across the genome were shaped by genomic structural features and evolutionary history, with stronger selection and reduced introgression in regions of low recombination, high subspecies differentiation, positive selection within the subspecies, and on macrochromosomes. Cumulatively, these results suggest that linkage among loci interacts with divergent selection and past divergent evolution between species to strengthen barriers to gene flow within hybrid zones.

Rapid evolution of recombination landscapes during the divergence of cichlid ecotypes in Lake Masoko

Variation of recombination rate along the genome is of crucial importance to rapid adaptation and organismal diversification. Many unknowns remain regarding how and why recombination landscapes evolve in nature. Here, we reconstruct recombination maps based on linkage disequilibrium and use subsampling and simulations to derive a new measure of recombination landscape evolution: the Population Recombination Divergence Index (PRDI). Using PRDI, we show that fine-scale recombination landscapes differ substantially between two cichlid fish ecotypes of Astatotilapia calliptera that diverged only ~2,500 generations ago. Perhaps surprisingly, recombination landscape differences are not driven by divergence in terms of allele frequency (FST) and nucleotide diversity (Δ(π)): although there is some association, we observe positive PRDI in regions where FST and Δ(π) are zero. We found a stronger association between the evolution of recombination and 47 large haplotype blocks that are polymorphic in Lake Masoko, cover 21% of the genome, and appear to include multiple inversions. Among haplotype blocks, there is a strong and clear association between the degree of recombination divergence and differences between ecotypes in heterozygosity, consistent with recombination suppression in heterozygotes. Overall, our work provides a holistic view of changes in population recombination landscapes during the early stages of speciation with gene flow.

Plant sperm cell sequencing for genome phasing and determination of meiotic crossover points

Haplotype phasing represents a pivotal procedure in genome analysis, entailing the identification of specific genetic variant combinations on each chromosome. Achieving chromosome-level genome phasing constitutes a considerable challenge, particularly in organisms with large and complex genomes. To address this challenge, we have developed a robust, gamete cell-based phasing pipeline, including wet-laboratory processes for plant sperm cell isolation, short-read sequencing and a bioinformatics workflow to generate chromosome-level phasing. The bioinformatics workflow is applicable for both plant and other sperm cells, for example, those of mammals. Our pipeline ensures high-quality single-nucleotide polymorphism (SNP) calling for each sperm cell and the subsequent construction of a high-density genetic map. The genetic map facilitates accurate chromosome-level genome phasing, enables crossover event detection and could be used to correct potential assembly errors. Our bioinformatics pipeline runs on a Linux system and most of its steps can be executed in parallel, expediting the analysis process. The entire workflow can be performed over the course of 1 d. We provide a practical example from our previous research using this protocol and provide the whole bioinformatics pipeline as a Docker image to ensure its easy adaptability to other studies.

Agricultural landscape genomics to increase crop resilience

Populations are continually adapting to their environment. Knowledge of which populations and individuals harbor unique and agriculturally useful variations has the potential to accelerate crop adaptation to the increasingly challenging environments predicted for the coming century. Landscape genomics, which identifies associations between environmental and genomic variation, provides a means for obtaining this knowledge. However, despite extensive efforts to assemble and characterize ex situ collections of crops and their wild relatives, gaps remain in the genomic and environmental datasets needed to robustly implement this approach. This article outlines the history of landscape genomics, which, to date, has mainly been used in conservation and evolutionary studies, provides an overview of crops and wild relative collections that have the necessary data for implementation and identifies areas where new data generation is needed. We find that 60% of the crops covered by the International Treaty on Plant Genetic Resources for Food and Agriculture lack the data necessary to conduct this kind of analysis, necessitating identification of crops in need of more collections, sequencing, or phenotyping. By highlighting these aspects, we aim to help develop agricultural landscape genomics as a sub-discipline that brings together evolutionary genetics, landscape ecology, and plant breeding, ultimately enhancing the development of resilient and adaptable crops for future environmental challenges.

Incomplete recombination suppression fuels extensive haplotype diversity in a butterfly colour pattern supergene

Supergenes can evolve when recombination-suppressing mechanisms like inversions promote co-inheritance of alleles at two or more polymorphic loci that affect a complex trait. Theory shows that such genetic architectures can be favoured under balancing selection or local adaptation in the face of gene flow, but they can also bring costs associated with reduced opportunities for recombination. These costs may in turn be offset by rare 'gene flux' between inverted and ancestral haplotypes, with a range of possible outcomes. We aimed to shed light on these processes by investigating the 'BC supergene', a large genomic region comprising multiple rearrangements associated with three distinct wing colour morphs in Danaus chrysippus, a butterfly known as the African monarch, African queen and plain tiger. Using whole-genome resequencing data from 174 individuals, we first confirm the effects of BC on wing colour pattern: background melanism is associated with SNPs in the promoter region of yellow, within an inverted subregion of the supergene, while forewing tip pattern is most likely associated with copy-number variation in a separate subregion of the supergene. We then show that haplotype diversity within the supergene is surprisingly extensive: there are at least six divergent haplotype groups that experience suppressed recombination with respect to each other. Despite high divergence between these haplotype groups, we identify an unexpectedly large number of natural recombinant haplotypes. Several of the inferred crossovers occurred between adjacent inversion 'modules', while others occurred within inversions. Furthermore, we show that new haplotype groups have arisen through recombination between two pre-existing ones. Specifically, an allele for dark colouration in the promoter of yellow has recombined into distinct haplotype backgrounds on at least two separate occasions. Overall, our findings paint a picture of dynamic evolution of supergene haplotypes, fuelled by incomplete recombination suppression.

Inversions contribute disproportionately to parallel genomic divergence in dune sunflowers

The probability of parallel genetic evolution is a function of the strength of selection and constraints imposed by genetic architecture. Inversions capture locally adapted alleles and suppress recombination between them, which limits the range of adaptive responses. In addition, the combined phenotypic effect of alleles within inversions is likely to be greater than that of individual alleles; this should further increase the contributions of inversions to parallel evolution. We tested the hypothesis that inversions contribute disproportionately to parallel genetic evolution in independent dune ecotypes of Helianthus petiolaris. We analysed habitat data and identified variables underlying parallel habitat shifts. Genotype-environment association analyses of these variables indicated parallel responses of inversions to shared selective pressures. We also confirmed larger seed size across the dunes and performed quantitative trait locus mapping with multiple crosses. Quantitative trait loci shared between locations fell into inversions more than expected by chance. We used whole-genome sequencing data to identify selective sweeps in the dune ecotypes and found that the majority of shared swept regions were found within inversions. Phylogenetic analyses of shared regions indicated that within inversions, the same allele typically was found in the dune habitat at both sites. These results confirm predictions that inversions drive parallel divergence in the dune ecotypes.

The recombination landscape of the barn owl, from families to populations

Homologous recombination is a meiotic process that generates diversity along the genome and interacts with all evolutionary forces. Despite its importance, studies of recombination landscapes are lacking due to methodological limitations and limited data. Frequently used approaches include linkage mapping based on familial data that provides sex-specific broad-scale estimates of realized recombination and inferences based on population linkage disequilibrium that reveal a more fine-scale resolution of the recombination landscape, albeit dependent on the effective population size and the selective forces acting on the population. In this study, we use a combination of these 2 methods to elucidate the recombination landscape for the Afro-European barn owl (Tyto alba). We find subtle differences in crossover placement between sexes that lead to differential effective shuffling of alleles. Linkage disequilibrium-based estimates of recombination are concordant with family-based estimates and identify large variation in recombination rates within and among linkage groups. Larger chromosomes show variation in recombination rates, while smaller chromosomes have a universally high rate that shapes the diversity landscape. We find that recombination rates are correlated with gene content, genetic diversity, and GC content. We find no conclusive differences in the recombination landscapes between populations. Overall, this comprehensive analysis enhances our understanding of recombination dynamics, genomic architecture, and sex-specific variation in the barn owl, contributing valuable insights to the broader field of avian genomics.

Shared Selection and Genetic Architecture Drive Strikingly Repeatable Evolution in Long-Term Experimental Hybrid Populations

The degree to which evolution repeats itself has implications regarding the major forces driving evolution and the potential for evolutionary biology to be a predictive (vs. solely historical) science. To understand the factors that control evolutionary repeatability, we experimentally evolved four replicate hybrid populations of sunflowers at natural sites for up to 14 years and tracked ancestry across the genome. We found that there was very strong negative selection against introgressed ancestry in several chromosomes, but positive selection for introgressed ancestry in one chromosome. Further, the strength of selection was influenced by recombination rate. High recombination regions had lower selection against introgressed ancestry due to more frequent recombination away from incompatible backgrounds. Strikingly, evolution was highly parallel across replicates, with shared selection driving 88% of variance in introgressed allele frequency change. Parallel evolution was driven by both high levels of sustained linkage in introgressed alleles and strong selection on large-effect quantitative trait loci. This work highlights the repeatability of evolution through hybridization and confirms the central roles that natural selection, genomic architecture, and recombination play in the process.

A High-Quality Phased Genome Assembly of Stinging Nettle (Urtica dioica ssp. dioica)

Stinging nettles (Urtica dioica) have a long history of association with human civilization, having been used as a source of textile fibers, food and medicine. Here, we present a chromosome-level, phased genome assembly for a diploid female clone of Urtica dioica from Romania. Using a combination of PacBio HiFi, Oxford Nanopore, and Illumina sequencing, as well as Hi-C long-range interaction data (using a novel Hi-C protocol presented here), we assembled two haplotypes of 574.9 Mbp (contig N50 = 10.9 Mbp, scaffold N50 = 44.0 Mbp) and 521.2 Mbp (contig N50 = 13.5 Mbp, scaffold N50 = 48.0 Mbp), with assembly BUSCO scores of 92.6% and 92.2%. We annotated 20,333 and 20,140 genes for each haplotype, covering over 90% of the complete BUSCO genes and including two copies of a gene putatively encoding the neurotoxic peptide urthionin, which could contribute to nettle's characteristic sting. Despite its relatively small size, the nettle genome displays very high levels of repetitiveness, with transposable elements comprising more than 60% of the genome, as well as considerable structural variation. This genome assembly represents an important resource for the nettle community and will enable the investigation of the genetic basis of the many interesting characteristics of this species.

Population Genomics of Adaptive Radiation

Adaptive radiations are rich laboratories for exploring, testing, and understanding key theories in evolution and ecology because they offer spectacular displays of speciation and ecological adaptation. Particular challenges to the study of adaptive radiation include high levels of species richness, rapid speciation, and gene flow between species. Over the last decade, high-throughput sequencing technologies and access to population genomic data have lessened these challenges by enabling the analysis of samples from many individual organisms at whole-genome scales. Here we review how population genomic data have facilitated our knowledge of adaptive radiation in five key areas: (1) phylogenetics, (2) hybridization, (3) timing and rates of diversification, (4) the genomic basis of trait evolution, and (5) the role of genome structure in divergence. We review current knowledge in each area, highlight outstanding questions, and focus on methods that facilitate detection of complex patterns in the divergence and demography of populations through time. It is clear that population genomic data are revolutionising the ability to reconstruct evolutionary history in rapidly diversifying clades. Additionally, studies are increasingly emphasising the central role of gene flow, re-use of standing genetic variation during adaptation, and structural genomic elements as facilitators of the speciation process in adaptive radiations. We highlight hybridization-and the hypothesized processes by which it shapes diversification-and questions seeking to bridge the divide between microevolutionary and macroevolutionary processes as rich areas for future study. Overall, access to population genomic data has facilitated an exciting era in adaptive radiation research, with implications for deeper understanding of fundamental evolutionary processes across the tree of life.

Pangenome graphs and their applications in biodiversity genomics

Complete datasets of genetic variants are key to biodiversity genomic studies. Long-read sequencing technologies allow the routine assembly of highly contiguous, haplotype-resolved reference genomes. However, even when complete, reference genomes from a single individual may bias downstream analyses and fail to adequately represent genetic diversity within a population or species. Pangenome graphs assembled from aligned collections of high-quality genomes can overcome representation bias by integrating sequence information from multiple genomes from the same population, species or genus into a single reference. Here, we review the available tools and data structures to build, visualize and manipulate pangenome graphs while providing practical examples and discussing their applications in biodiversity and conservation genomics across the tree of life.

Combining Individual-Based Radio-Tracking With Whole-Genome Sequencing Data Reveals Candidate for Genetic Basis of Partial Migration in a Songbird

Partial migration is a phenomenon where migratory and resident individuals of the same species co-exist within a population, and has been linked to both intrinsic (e.g., genetic) as well as environmental factors. Here we investigated the genomic architecture of partial migration in the common blackbird, a songbird that comprises resident populations in the southern distribution range, partial migratory populations in central Europe, and exclusively migratory populations in northern and eastern Europe. We generated whole-genome sequencing data for 60 individuals, each of which was phenotyped for migratory behavior using radio-telemetry tracking. These individuals were sampled across the species' distribution range, including resident populations (Spain and France), obligate migrants (Russia), and a partial migratory population with equal numbers of migratory and resident individuals in Germany. We estimated genetic differentiation (FST) of single-nucleotide variants (SNVs) in 2.5 kb windows between all possible population and migratory phenotype combinations, and focused our characterization on birds from the partial migratory population in Germany. Despite overall low differentiation within the partial migratory German population, we identified several outlier regions with elevated differentiation on four distinct chromosomes. The region with the highest relative and absolute differentiation was located on chromosome 9, overlapping PER2, which has previously been shown to be involved in the control of the circadian rhythm across vertebrates. While this region showed high levels of differentiation, no fixed variant could be identified, supporting the notion that a complex phenotype such as migratory behavior is likely controlled by a large number of genetic loci.

Centromeres are hotspots for chromosomal inversions and breeding traits in mango

Chromosomal inversions can preserve combinations of favorable alleles by suppressing recombination. Simultaneously, they reduce the effectiveness of purifying selection enabling deleterious alleles to accumulate. This study explores how areas of low recombination, including centromeric regions and chromosomal inversions, contribute to the accumulation of deleterious and favorable loci in 225 Mangifera indica genomes from the Australian Mango Breeding Program. Here, we identify 17 chromosomal inversions that cover 7.7% (29.7 Mb) of the M. indica genome: eight pericentric (inversion includes the centromere) and nine paracentric (inversion is on one arm of the chromosome). Our results show that these large pericentric inversions are accumulating deleterious loci, while the paracentric inversions show deleterious levels above and below the genome wide average. We find that despite their deleterious load, chromosomal inversions contain small effect loci linked to variation in crucial breeding traits. These results indicate that chromosomal inversions have likely facilitated the evolution of key mango breeding traits. Our study has important implications for selective breeding of favorable combinations of alleles in regions of low recombination.

The evolution of genomic organization through chromosomal rearrangements in Arctic charr (Salvelinus alpinus)

Chromosomal rearrangements (CRs) can play an important role in evolutionary diversification by preserving linkage among favourable alleles through reduced recombination and/or by reducing hybrid fitness due to genomic incompatibilities. Our goal was to determine to what extent CRs contribute to known patterns of genetic variation in Arctic charr (Salvelinus alpinus). To address this goal, we compared genetic linkage maps to identify whole arm CRs and smaller scale structural variants (SVs) such as translocations/transpositions and inversions found in groups of populations that reflect the temporal sequence of geographic isolation events. If CRs contribute to genetic differentiation, we expected that CRs would be specific to glacial lineages, geographic clusters of populations within lineages, and sympatric morphs. We detected fusions and fissions of whole chromosome arms and SV involving translocations/transpositions of the sex-determining gene (sdY) and inversions. Several CRs were shared across populations from the Arctic and Atlantic glacial lineages, Canadian and Icelandic populations within the Atlantic lineage, between two Icelandic populations and sympatric morphs within Icelandic populations, suggesting that their origin predates geographic isolation in glacial refugia. Other CRs were specific to single populations, which suggests a more recent origin of these variants in refugia, during post-glacial recolonization and/or in contemporary populations. Thus, CRs contribute relatively little to known patterns of genetic differentiation at different geographic scales but represent a pool of standing genetic variation for evolution.

Complex Genomic Landscape of Inversion Polymorphism in Europe's Most Destructive Forest Pest

In many species, polymorphic genomic inversions underlie complex phenotypic polymorphisms and facilitate local adaptation in the face of gene flow. Multiple polymorphic inversions can co-occur in a genome, but the prevalence, evolutionary significance, and limits to complexity of genomic inversion landscapes remain poorly understood. Here, we examine genome-wide genetic variation in one of Europe's most destructive forest pests, the spruce bark beetle Ips typographus, scan for polymorphic inversions, and test whether inversions are associated with key traits in this species. We analyzed 240 individuals from 18 populations across the species' European range and, using a whole-genome resequencing approach, identified 27 polymorphic inversions covering ∼28% of the genome. The inversions vary in size and in levels of intra-inversion recombination, are highly polymorphic across the species range, and often overlap, forming a complex genomic architecture. We found no support for mechanisms such as directional selection, overdominance, and associative overdominance that are often invoked to explain the presence of large inversion polymorphisms in the genome. This suggests that inversions are either neutral or maintained by the combined action of multiple evolutionary forces. We also found that inversions are enriched in odorant receptor genes encoding elements of recognition pathways for host plants, mates, and symbiotic fungi. Our results indicate that the genome of this major forest pest of growing social, political, and economic importance harbors one of the most complex inversion landscapes described to date and raise questions about the limits of intraspecific genomic architecture complexity.

The Role of Hybridization in Species Formation and Persistence

Hybridization, or interbreeding between different taxa, was traditionally considered to be rare and to have a largely detrimental impact on biodiversity, sometimes leading to the breakdown of reproductive isolation and even to the reversal of speciation. However, modern genomic and analytical methods have shown that hybridization is common in some of the most diverse clades across the tree of life, sometimes leading to rapid increase of phenotypic variability, to introgression of adaptive alleles, to the formation of hybrid species, and even to entire species radiations. In this review, we identify consensus among diverse research programs to show how the field has progressed. Hybridization is a multifaceted evolutionary process that can strongly influence species formation and facilitate adaptation and persistence of species in a rapidly changing world. Progress on testing this hypothesis will require cooperation among different subdisciplines.

Distinct patterns of genetic variation at low-recombining genomic regions represent haplotype structure

Genomic regions sometimes show patterns of genetic variation distinct from the genome-wide population structure. Such deviations have often been interpreted to represent effects of selection. However, systematic investigation of whether and how non-selective factors, such as recombination rates, can affect distinct patterns has been limited. Here, we associate distinct patterns of genetic variation with reduced recombination rates in a songbird, the Eurasian blackcap (Sylvia atricapilla), using a new reference genome assembly, whole-genome resequencing data and recombination maps. We find that distinct patterns of genetic variation reflect haplotype structure at genomic regions with different prevalence of reduced recombination rate across populations. At low-recombining regions shared in most populations, distinct patterns reflect conspicuous haplotypes segregating in multiple populations. At low-recombining regions found only in a few populations, distinct patterns represent variance among cryptic haplotypes within the low-recombining populations. With simulations, we confirm that these distinct patterns evolve neutrally by reduced recombination rate, on which the effects of selection can be overlaid. Our results highlight that distinct patterns of genetic variation can emerge through evolutionary reduction of local recombination rate. The recombination landscape as an evolvable trait therefore plays an important role determining the heterogeneous distribution of genetic variation along the genome.

Prezygotic barriers effectively limit hybridization in a rapid evolutionary radiation

Hybridization is increasingly recognized as an important evolutionary process across the tree of life. In many clades, phylogenomic approaches have permitted unparalleled insight into the extent and frequency of hybridization. However, we continue to lack a deep understanding of the factors that limit and shape patterns of hybridization, especially in evolutionary radiations. In this study, we characterized patterns of introgression across Costus (Costaceae), a young evolutionary radiation of tropical understory plants that maintain widespread interfertility despite exhibiting strong prezygotic reproductive isolation. We analyzed a phylogenomic dataset of 756 genes from 54 Costus species using multiple complementary approaches - D-statistics, gene-tree-based tests, and phylogenetic network analyses - to detect and characterize introgression events throughout the evolutionary history of the radiation. Our results identified a moderate number of introgression events, including a particularly ancient, well-supported event spanning one of the deepest divergences in the clade. Most introgression events occurred between taxa or ancestral lineages that shared the same pollination syndrome (bee-pollinated or hummingbird-pollinated). These findings suggest that prezygotic barriers, including pollinator specialization, have been key to the balance between introgression and reproductive isolation in Costus.

Fitness consequences of structural variation inferred from a House Finch pangenome

Genomic structural variants (SVs) play a crucial role in adaptive evolution, yet their average fitness effects and characterization with pangenome tools are understudied in wild animal populations. We constructed a pangenome for House Finches (Haemorhous mexicanus), a model for studies of host-pathogen coevolution, using long-read sequence data on 16 individuals (32 de novo-assembled haplotypes) and one outgroup. We identified 887,118 SVs larger than 50 base pairs, mostly (60%) involving repetitive elements, with reduced SV diversity in the eastern US as a result of its introduction by humans. The distribution of fitness effects of genome-wide SVs was estimated using maximum likelihood approaches and revealed that SVs in both coding and noncoding regions were on average more deleterious than smaller indels or single nucleotide polymorphisms. The reference-free pangenome facilitated identification of a > 10-My-old, 11-megabase-long pericentric inversion on chromosome 1. We found that the genotype frequencies of the inversion, estimated from 135 birds widely sampled temporally and geographically, increased steadily over the 25 y since House Finches were first exposed to the bacterial pathogen Mycoplasma gallisepticum and showed signatures of balancing selection, capturing genes related to immunity and telomerase activity. We also observed shorter telomeres in populations with a greater number of years exposure to Mycoplasma. Our study illustrates the utility of long-read sequencing and pangenome methods for understanding wild animal populations, estimating fitness effects of genome-wide SVs, and advancing our understanding of adaptive evolution through structural variation.

Partial unidirectional translocation from 5AL to 7BS leads to dense spike in an EMS-induced wheat mutant

BACKGROUND

As the inflorescence of wheat, spike architecture largely determines grain productivity. Dissecting the genetic basis for the spike morphology of wheat can contribute to the designation of ideal spike morphology to improve grain production.

RESULTS

The present study characterizes a dense spike1 (ds1) mutant, derived from Nongda3753, induced by EMS treatment, which exhibits a dense spike and reduced plant height. Through bulked segregant analysis sequencing (BSA-Seq) of two segregating populations, ds1 was mapped to the short arm of chromosome 7B. Further genotypic and phenotypic analyses of the residual heterozygous lines from F3 to F6 of Yong3002×ds1 revealed that there was a 0-135 Mb deletion in chromosome 7B associated with the dense spike phenotype. The reads count analysis of the two bulks in BSA-Seq, along with the cytological analysis of ds1, ND3753, NIL-ds1 and NIL-Y3002, confirmed that the partial unidirectional translocation of 5AL (543-713 Mb) to 7BS (0-135 Mb) exists in ds1. This translocation led to an increase in both copy number and expression of the Q gene, which is one of the reasons for the dense spike phenotype observed in ds1.

CONCLUSION

Partial unidirectional translocation from 5AL to 7BS was identified in the EMS-induced mutant ds1, which exhibits dense spike phenotype. This research illustrates the effect of one chromosome structure variation on wheat spike morphology, and provides new materials with several chromosome structure variations for future wheat breeding.

Admixture With Cultivated Sunflower Likely Facilitated Establishment and Spread of Wild Sunflower (Helianthus annuus) in Argentina

A better understanding of the genetic and ecological factors underlying successful invasions is critical to mitigate the negative impacts of invasive species. Here, we study the invasion history of Helianthus annuus populations from Argentina, with particular emphasis on the role of post-introduction admixture with cultivated sunflower (also H. annuus) and climate adaptation driven by large haploblocks. We conducted genotyping-by-sequencing of samples of wild populations as well as Argentinian cultivars and compared them with wild (including related annual Helianthus species) and cultivated samples from the native range. We also characterised samples for 11 known haploblocks associated with environmental variation in native populations to test whether haploblocks contributed to invasion success. Population genomics analyses supported two independent geographic sources for Argentinian populations, the central United States and Texas, but no significant contribution of related annual Helianthus species. We found pervasive admixture with cultivated sunflower, likely as result of post-introduction hybridization. Genomic scans between invasive populations and their native sources identified multiple genomic regions of divergence, possibly indicative of selection, in the invaded range. These regions significantly overlapped between the two native-invasive comparisons and showed disproportionally high crop ancestry, suggesting that crop alleles contributed to invasion success. We did not find evidence of climate adaptation mediated by haploblocks, yet outliers of genome scans were enriched in haploblock regions and, for at least two haploblocks, the cultivar haplotype was favoured in Argentina. Our results show that admixture with cultivated sunflower played a major role in the establishment and spread of H. annuus populations in Argentina.

Historic breeding practices contribute to germplasm divergence in leaf specialized metabolism and ecophysiology in cultivated sunflower (Helianthus annuus)

PREMISE

The use of hybrid breeding systems to increase crop yields has been the cornerstone of modern agriculture and is exemplified in the breeding and improvement of cultivated sunflower (Helianthus annuus). However, it is poorly understood what effect supporting separate breeding pools in such systems, combined with continued selection for yield, may have on leaf ecophysiology and specialized metabolite variation.

METHODS

We analyzed 288 lines of cultivated H. annuus to examine the genomic basis of several specialized metabolites and agronomically important traits across major heterotic groups.

RESULTS

Heterotic group identity supports phenotypic divergences between fertility restoring and cytoplasmic male-sterility maintainer lines in leaf ecophysiology and specialized metabolism. However, the divergence is not associated with physical linkage to nuclear genes that support current hybrid breeding practices in cultivated H. annuus. Additionally, we identified four genomic regions associated with leaf ecophysiology and specialized metabolism that colocalize with previously identified QTLs for quantitative self-compatibility traits and with S-protein homolog (SPH) proteins, a recently discovered family of proteins associated with self-incompatibility and self/nonself recognition in Papaver rhoeas (common poppy) with suggested conserved downstream mechanisms among eudicots.

CONCLUSIONS

Further work is necessary to confirm the self-incompatibility mechanisms in cultivated H. annuus and their relationship to the integrative and polygenic architecture of leaf ecophysiology and specialized metabolism in cultivated sunflower. However, because self-compatibility is a derived quantitative trait in cultivated H. annuus, trait linkage to divergent phenotypic traits may have partially arisen as a potential unintended consequence of historical breeding practices and selection for yield.

Inversions encounter relaxed genetic constraints and balance birth and death of TPS genes in Curcuma

Evolutionary dynamics of inversion and its impact on biochemical traits are a puzzling question. Here, we show abundance of inversions in three Curcuma species (turmeric, hidden ginger and Siam tulip). Genes within inversions display higher long terminal repeat content and lower expression level compared with genomic background, suggesting inversions in Curcuma experience relaxed genetic constraints. It is corroborated by depletion of selected SNPs and enrichment of deleterious mutations in inversions detected among 56 Siam tulip cultivars. Functional verification of tandem duplicated terpene synthase (TPS) genes reveals that genes within inversions become pseudogenes, while genes outside retain catalytic function. Our findings suggest that inversions act as a counteracting force against tandem duplication in balancing birth and death of TPS genes and modulating terpenoid contents in Curcuma. This study provides an empirical example that inversions are likely not adaptive but affect biochemical traits.

Population-aware permutation-based significance thresholds for genome-wide association studies

Motivation

Permutation-based significance thresholds have been shown to be a robust alternative to classical Bonferroni significance thresholds in genome-wide association studies (GWAS) for skewed phenotype distributions. The recently published method permGWAS introduced a batch-wise approach to efficiently compute permutation-based GWAS. However, running multiple univariate tests in parallel leads to many repetitive computations and increased computational resources. More importantly, traditional permutation methods that permute only the phenotype break the underlying population structure.

Results

We propose permGWAS2, an improved method that does not break the population structure during permutations and uses an elegant block matrix decomposition to optimize computations, thereby reducing redundancies. We show on synthetic data that this improved approach yields a lower false discovery rate for skewed phenotype distributions compared to the previous version and the commonly used Bonferroni correction. In addition, we re-analyze a dataset covering phenotypic variation in 86 traits in a population of 615 wild sunflowers (Helianthus annuus L.). This led to the identification of dozens of novel associations with putatively adaptive traits, and removed several likely false-positive associations with limited biological support.

Availability and implementation

permGWAS2 is open-source and publicly available on GitHub for download: https://github.com/grimmlab/permGWAS.

Haplotype-resolved T2T genome assemblies and pangenome graph of pear reveal diverse patterns of allele-specific expression and the genomic basis of fruit quality traits

Hybrid crops often exhibit increased yield and greater resilience, yet the genomic mechanism(s) underlying hybrid vigor or heterosis remain unclear, hindering our ability to predict the expression of phenotypic traits in hybrid breeding. Here, we generated haplotype-resolved T2T genome assemblies of two pear hybrid varieties, 'Yuluxiang' (YLX) and 'Hongxiangsu' (HXS), which share the same maternal parent but differ in their paternal parents. We then used these assemblies to explore the genome-scale landscape of allele-specific expression (ASE) and create a pangenome graph for pear. ASE was observed for close to 6000 genes in both hybrid cultivars. A subset of ASE genes related to aspects of fruit quality such as sugars, organic acids, and cuticular wax were identified, suggesting their important contributions to heterosis. Specifically, Ma1, a gene regulating fruit acidity, is absent in the paternal haplotypes of HXS and YLX. A pangenome graph was built based on our assemblies and seven published pear genomes. Resequencing data for 139 cultivated pear genotypes (including 97 genotypes sequenced here) were subsequently aligned to the pangenome graph, revealing numerous structural variant hotspots and selective sweeps during pear diversification. As predicted, the Ma1 allele was found to be absent in varieties with low organic acid content, and this association was functionally validated by Ma1 overexpression in pear fruit and calli. Overall, these results reveal the contributions of ASE to fruit-quality heterosis and provide a robust pangenome reference for high-resolution allele discovery and association mapping.

Mutation potentiates migration swamping in polygenic local adaptation

Locally adapted traits can exhibit a wide range of genetic architectures, from pronounced divergence at a few loci to small frequency divergence at many loci. The type of architecture that evolves depends strongly on the migration rate, as weakly selected loci experience swamping and do not make lasting contributions to divergence. Simulations from previous studies showed that even when mutations are strongly selected and should resist migration swamping, the architecture of adaptation can collapse and become transient at high mutation rates. Here, we use an analytical two-population model to study how this transition in genetic architecture depends upon population size, strength of selection, and parameters describing the mutation process. To do this, we develop a mathematical theory based on the diffusion approximation to predict the threshold mutation rate above which the transition occurs. We find that this performs well across a wide range of parameter space, based on comparisons with individual-based simulations. The threshold mutation rate depends most strongly on the average effect size of mutations, weakly on the strength of selection, and marginally on the population size. Across a wide range of the parameter space, we observe that the transition to a transient architecture occurs when the trait-wide mutation rate is 10-3-10-2, suggesting that this phenomenon is potentially relevant to complex traits with a large mutational target. On the other hand, based on the apparent stability of genetic architecture in many classic examples of local adaptation, our theory suggests that per-trait mutation rates are often relatively low.

36-year study reveals stability of a wild wheat population across microhabitats

Long-term genetic studies of wild populations are very scarce, but are essential for connecting ecological and population genetics models, and for understanding the dynamics of biodiversity. We present a study of a wild wheat population sampled over a 36-year period at high spatial resolution. We genotyped 832 individuals from regular sampling along transects during the course of the experiment. Genotypes were clustered into ecological microhabitats over scales of tens of metres, and this clustering was remarkably stable over the 36 generations of the study. Simulations show that it is difficult to determine whether this spatial and temporal stability reflects extremely limited dispersal or fine-scale local adaptation to ecological parameters. Using a common-garden experiment, we showed that the genotypes found in distinct microhabitats differ phenotypically. Our results provide a rare insight into the population genetics of a natural population over a long monitoring period.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these 2 fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we lay out a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., genome-wide association studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur analytically and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study, we re-examine an analysis testing for coevolution of expression levels between genes across a fungal phylogeny and show that including eigenvectors of the covariance matrix as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

Dissecting the genetic architecture of sunflower disc diameter using genome-wide association study

Sunflower (Helianthus annuus L.) plays an essential role in meeting the demand for edible oil worldwide. The yield of sunflower seeds encompasses several component traits, including the disc diameter. Over three consecutive years, 2019, 2020, and 2022, we assessed phenotypic variation in disc diameter across a diverse set of sunflower accessions (N = 342) in replicated field trials. Upon aggregating the phenotypic data from multiple years, we estimated the broad sense heritability (H 2) of the disc diameter trait to be 0.88. A subset of N = 274 accessions was genotyped by using the tunable genotyping-by-sequencing (tGBS) method, resulting in 226,779 high-quality SNPs. Using these SNPs and the disc diameter phenotype, we conducted a genome-wide association study (GWAS) employing two statistical approaches: the mixed linear model (MLM) and the fixed and random model circulating probability unification (farmCPU). The MLM and farmCPU GWAS approaches identified 106 and 8 significant SNPs located close to 53 and 21 genes, respectively. The MLM analysis identified two significant peaks: a prominent signal on chromosome 10 and a relatively weaker signal on chromosome 16, both of which were also detected by farmCPU. The genetic loci associated with disc diameter, as well as the related candidate genes, present promising avenues for further functional validation and serve as a basis for sunflower oil yield improvement.

Structural variation discovery in wheat using PacBio high-fidelity sequencing

Structural variations (SVs) pervade plant genomes and contribute substantially to the phenotypic diversity. However, most SVs were ineffectively assayed due to their complex nature and the limitations of early genomic technologies. By applying the PacBio high-fidelity (HiFi) sequencing for wheat genomes, we performed a comprehensive evaluation of mainstream long-read aligners and SV callers in SV detection. The results indicated that the accuracy of deletion discovery is markedly influenced by callers, accounting for 87.73% of the variance, whereas both aligners (38.25%) and callers (49.32%) contributed substantially to the accuracy variance for insertions. Among the aligners, Winnowmap2 and NGMLR excelled in detecting deletions and insertions, respectively. For SV callers, SVIM achieved the best performance. We demonstrated that combining the aligners and callers mentioned above is optimal for SV detection. Furthermore, we evaluated the effect of sequencing depth on the accuracy of SV detection, revealing that low-coverage HiFi sequencing is sufficiently robust for high-quality SV discovery. This study thoroughly evaluated SV discovery approaches and established optimal workflows for investigating structural variations using low-coverage HiFi sequencing in the wheat genome, which will advance SV discovery and decipher the biological functions of SVs in wheat and many other plants.

Phenotypic plasticity rather than ecotypic differentiation explains the broad realized niche of a Neotropical orchid species

Local adaptation is common in plant species, and knowing whether a population is locally adapted has fundamental and applied relevance. However, local adaptation in tropical plants remains largely less studied, and covering this gap is not simple since reciprocal transplantation - the gold standard for detecting local adaptation - is not feasible for most species. Here, we combined genetic, climatic and phenotypic data to investigate ecotypic differentiation, an important aspect of local adaptation, in coastal and inland populations of the orchid Epidendrum fulgens Brongn., a long-lived tropical plant for which reciprocal transplantation would not be feasible. We used nine microsatellite markers to estimate genetic divergence between inland and coastal populations. Moreover, occurrence data and climate data were used to test for differences in the realized niche of those populations. Finally, we assessed saturated water content, leaf specific area, height, and stomatal density in common garden and in situ to investigate the effects of ecotypic differentiation and plasticity on the phenotype. Coastal and inland groups' niches do not overlap, the former occupying a wetter and warmer area. However, this differentiation does not seem to be driven by ecotypic differentiation since there was no positive correlation between genetic structure and climate dissimilarity. Moreover, specific leaf area and leaf saturated water content, which are important phenotypic traits related to soil fertility and drought stress, were rather plastic. We conclude that ecotypic differentiation is absent, since phenotypic plasticity is an important mechanism explaining the niche broadness of this species.

The genetic architecture of repeated local adaptation to climate in distantly related plants

Closely related species often use the same genes to adapt to similar environments. However, we know little about why such genes possess increased adaptive potential and whether this is conserved across deeper evolutionary lineages. Adaptation to climate presents a natural laboratory to test these ideas, as even distantly related species must contend with similar stresses. Here, we re-analyse genomic data from thousands of individuals from 25 plant species as diverged as lodgepole pine and Arabidopsis (~300 Myr). We test for genetic repeatability based on within-species associations between allele frequencies in genes and variation in 21 climate variables. Our results demonstrate significant statistical evidence for genetic repeatability across deep time that is not expected under randomness, identifying a suite of 108 gene families (orthogroups) and gene functions that repeatedly drive local adaptation to climate. This set includes many orthogroups with well-known functions in abiotic stress response. Using gene co-expression networks to quantify pleiotropy, we find that orthogroups with stronger evidence for repeatability exhibit greater network centrality and broader expression across tissues (higher pleiotropy), contrary to the 'cost of complexity' theory. These gene families may be important in helping wild and crop species cope with future climate change, representing important candidates for future study.

Local adaptation to climate facilitates a global invasion

Local adaptation may facilitate range expansion during invasions, but the mechanisms promoting destructive invasions remain unclear. Cheatgrass (Bromus tectorum), native to Eurasia and Africa, has invaded globally, with particularly severe impacts in western North America. We sequenced 307 genotypes and conducted controlled experiments. We found that diverse lineages invaded North America, where long-distance gene flow is common. Ancestry and phenotypic clines in the native range predicted those in the invaded range, indicating pre-adapted genotypes colonized different regions. Common gardens showed directional selection on flowering time that reversed between warm and cold sites, potentially maintaining clines. In the Great Basin, genomic predictions of strong local adaptation identified sites where cheatgrass is most dominant. Preventing new introductions that may fuel adaptation is critical for managing ongoing invasions.

Genotype of pioneer plant Miscanthus is not a key factor in the structure of rhizosphere bacterial community in heavy metal polluted sites

Miscanthus is a common pioneer plant with abundant genetic variation in abandoned mines in southern China. However, the extent to which genetic differentiation among species modulates rhizosphere bacterial communities remains unclear. Miscanthus samples were collected from 26 typical abandoned heavy-metal mines with different soil types in southern China, tested using 14 pairs of simple sequence repeats (SSR) primers, and classified into two genotypes based on Nei's genetic distance. The structure and diversity of rhizosphere bacterial communities were examined using 16 S rRNA sequencing. The results showed that among the factors affecting the rhizosphere bacterial community structure of Miscanthus samples, the role of genotype was not significant, and geographical conditions were the most important factors, followed by pH and total organic carbon (TOC). The process of rhizospheric community assembly varied among different genotypes; however, the recruited species and their abundances were similar. Collectively, we provided an approach based on genetic differentiation to quantify the relative contribution of genotypes to the rhizosphere bacterial community, demonstrating that genotypes contribute less than soil conditions. Our findings provide new insights into the role of host genetics in the ecological processes of plant rhizosphere bacterial communities in abandoned mines and provide theoretical support for microbe-assisted phytoremediation.

Genomic Insights into Disease Resistance in Sunflower (Helianthus annuus): Identifying Key Regions and Candidate Genes for Verticillium dahliae Resistance

Sunflower (Helianthus annuus) is a globally significant field crop, and disease resistance is crucial for ensuring yield stability and crop quality. Verticillium dahliae is a notorious soilborne pathogen that causes Verticillium Wilt (VW) and threatens sunflower production worldwide. In this study, we conducted a comprehensive assessment of sunflower resistance to V. dahliae across 231 sunflower cultivar lines, from the Sunflower Association Mapping (SAM) population. We employed EMMAX and ridge regression best linear unbiased prediction (rrBLUP) and identified 148 quantitative trait loci (QTLs) and 23 putative genes associated with V. dahliae resistance, including receptor like kinases, cell wall modification, transcriptional regulation, plant stress signalling and defense regulation genes. Our enrichment and quantitative real-time PCR validation results highlight the importance of membrane vesicle trafficking in the sunflower immune system for efficient signaling and defense upon activation by V. dahliae. This study also reveals the polygenic architecture of V. dahliae resistance in sunflowers and provides insights for breeding sunflower cultivars resistant to VW. This research contributes to ongoing efforts to enhance crop resilience and reduce yield losses due to VW, ultimately benefiting sunflower growers and the agricultural sector.

PNNGS, a multi-convolutional parallel neural network for genomic selection

Genomic selection (GS) can accomplish breeding faster than phenotypic selection. Improving prediction accuracy is the key to promoting GS. To improve the GS prediction accuracy and stability, we introduce parallel convolution to deep learning for GS and call it a parallel neural network for genomic selection (PNNGS). In PNNGS, information passes through convolutions of different kernel sizes in parallel. The convolutions in each branch are connected with residuals. Four different Lp loss functions train PNNGS. Through experiments, the optimal number of parallel paths for rice, sunflower, wheat, and maize is found to be 4, 6, 4, and 3, respectively. Phenotype prediction is performed on 24 cases through ridge-regression best linear unbiased prediction (RRBLUP), random forests (RF), support vector regression (SVR), deep neural network genomic prediction (DNNGP), and PNNGS. Serial DNNGP and parallel PNNGS outperform the other three algorithms. On average, PNNGS prediction accuracy is 0.031 larger than DNNGP prediction accuracy, indicating that parallelism can improve the GS model. Plants are divided into clusters through principal component analysis (PCA) and K-means clustering algorithms. The sample sizes of different clusters vary greatly, indicating that this is unbalanced data. Through stratified sampling, the prediction stability and accuracy of PNNGS are improved. When the training samples are reduced in small clusters, the prediction accuracy of PNNGS decreases significantly. Increasing the sample size of small clusters is critical to improving the prediction accuracy of GS.

Regional differences in leaf evolution facilitate photosynthesis following severe drought

Characterizing physiological and anatomical changes that underlie rapid evolution following climatic perturbation can broaden our understanding of how climate change is affecting biodiversity. It can also provide evidence of cryptic adaptation despite stasis at higher levels of biological organization. Here, we compared evolutionary changes in populations of Mimulus cardinalis from historically different climates in the north and south of the species' range following an exceptional drought. We grew seeds produced from predrought ancestral plants alongside peak-drought descendants in a common glasshouse and exposed them to wet and dry conditions. Before the drought, northern ancestral populations expressed traits contributing to drought escape, while southern ancestral populations expressed drought avoidance. Following the drought, both regions evolved to reduce water loss and maintain photosynthesis in dry treatments (drought avoidance), but via different anatomical alterations in stomata, trichomes, and palisade mesophyll. Additionally, southern populations lost the ability to take advantage of wet conditions. These results reveal rapid evolution towards drought avoidance at an anatomical level following an exceptional drought, but suggest that differences in the mechanisms between regions incur different trade-offs. This sheds light on the importance of characterizing underlying mechanisms for downstream life-history and macromorphological traits.

Karyotypic and phenotypic condensation in allotetraploid wheats accompanied with reproductive strategy transformation: from natural evolution to domestication

MAIN CONCLUSION

Allotetraploid wheat reflects evolutionary divergence and domestication convergence in the karyotypic and phenotypic evolution, accompanied with the transformation from r- strategy to K- strategy in reproductive fitness. Allotetraploid wheat, the progenitor of hexaploidy bread wheat, has undergone 300,000 years of natural evolution and 10,000 years of domestication. The variations in karyotype and phenotype as well as fertility fitness have not been systematically linked. Here, by combining fluorescent in situ hybridization with the quantification of phenotypic and reproductive traits, we compared the karyotype, vegetative growth phenotype and reproductive fitness among synthesized, wild and domesticated accessions of allotetraploid wheat. We detected that the wild accessions showed dramatically high frequencies of homologous recombination and copy number variations of simple sequence repeats (SSR) comparing with synthetic and domesticated accessions. The phenotypic traits reflected significant differences among the populations shaped by distinct evolutionary processes. The diversity observed in wild accessions was significantly greater than that in domesticated ones, particularly in traits associated with vegetative growth and spike morphology. We found that the active pollen of domesticated accessions exhibited greater potential of germination, despite a lower rate of active pollen compared with the wild accessions, indicating a transformation in reproductive fitness strategy for pollen development in domesticated accessions compared to the wild accessions, from r-strategy to K-strategy. Our results demonstrate the condensation of karyotype and phenotype from natural wild accessions to domesticated accessions in allotetraploid wheats. Ecological strategy transformation should be seriously considered from evolution to domestication in polyploid plants, especially crops, which may provide a perspective on the adaptive evolution of polyploid plants.

Chromosomal inversions harbour excess mutational load in the coral, Acropora kenti, on the Great Barrier Reef

The future survival of coral reefs in the Anthropocene depends on the capacity of corals to adapt as oceans warm and extreme weather events become more frequent. Targeted interventions designed to assist evolutionary processes in corals require a comprehensive understanding of the distribution and structure of standing variation, however, efforts to map genomic variation in corals have so far focussed almost exclusively on SNPs, overlooking structural variants that have been shown to drive adaptive processes in other taxa. Here, we show that the reef-building coral, Acropora kenti, harbours at least five large, highly polymorphic structural variants, all of which exhibit signatures of strongly suppressed recombination in heterokaryotypes, a feature commonly associated with chromosomal inversions. Based on their high minor allele frequency, uniform distribution across habitats and elevated genetic load, we propose that these inversions in A. kenti are likely to be under balancing selection. An excess of SNPs with high impact on protein-coding genes within these loci elevates their importance both as potential targets for adaptive selection and as contributors to genetic decline if coral populations become fragmented or inbred in future.

Association studies of salinity tolerance in sunflower provide robust breeding and selection strategies under climate change

Phytotoxic soil salinity is a global problem, and in the northern Great Plains and western Canada, salt accumulates on the surface of marine sediment soils with high water tables under annual crop cover, particularly near wetlands. Crop production can overcome saline-affected soils using crop species and cultivars with salinity tolerance along with changes in management practices. This research seeks to improve our understanding of sunflower (Helianthus annuus) genetic tolerance to high salinity soils. Genome-wide association was conducted using the Sunflower Association Mapping panel grown for two years in naturally occurring saline soils (2016 and 2017, near Indian Head, Saskatchewan, Canada), and six phenotypes were measured: days to bloom, height, leaf area, leaf mass, oil percentage, and yield. Plot level soil salinity was determined by grid sampling of soil followed by kriging. Three estimates of sunflower performance were calculated: (1) under low soil salinity (< 4 dS/m), (2) under high soil salinity (> 4 dS/m), and (3) plasticity (regression coefficient between phenotype and soil salinity). Fourteen loci were significant, with one instance of co-localization between a leaf area and a leaf mass locus. Some genomic regions identified as significant in this study were also significant in a recent greenhouse salinity experiment using the same panel. Also, some candidate genes underlying significant QTL have been identified in other plant species as having a role in salinity response. This research identifies alleles for cultivar improvement and for genetic studies to further elucidate salinity tolerance pathways.

loco-pipe: an automated pipeline for population genomics with low-coverage whole-genome sequencing

Summary

We developed loco-pipe, a Snakemake pipeline that seamlessly streamlines a set of essential population genomic analyses for low-coverage whole genome sequencing (lcWGS) data. loco-pipe is highly automated, easily customizable, massively parallelized, and thus is a valuable tool for both new and experienced users of lcWGS.

Availability and implementation

loco-pipe is published under the GPLv3. It is freely available on GitHub (github.com/sudmantlab/loco-pipe) and archived on Zenodo (doi.org/10.5281/zenodo.10425920).

Local genetic adaptation to habitat in wild chimpanzees

How populations adapt to their environment is a fundamental question in biology. Yet we know surprisingly little about this process, especially for endangered species such as non-human great apes. Chimpanzees, our closest living relatives, are particularly interesting because they inhabit diverse habitats, from rainforest to woodland-savannah. Whether genetic adaptation facilitates such habitat diversity remains unknown, despite having wide implications for evolutionary biology and conservation. Using 828 newly generated exomes from wild chimpanzees, we find evidence of fine-scale genetic adaptation to habitat. Notably, adaptation to malaria in forest chimpanzees is mediated by the same genes underlying adaptation to malaria in humans. This work demonstrates the power of non-invasive samples to reveal genetic adaptations in endangered populations and highlights the importance of adaptive genetic diversity for chimpanzees.