Heritable differences in abundance of bacterial rhizosphere taxa are correlated with fungal necrotrophic pathogen resistance.

Host-microbe interactions are increasingly recognized as important drivers of organismal health, growth, longevity and community-scale ecological processes. However, less is known about how genetic variation affects hosts' associated microbiomes and downstream phenotypes. We demonstrate that sunflower (Helianthus annuus) harbours substantial, heritable variation in microbial communities under field conditions. We show that microbial communities co-vary with heritable variation in resistance to root infection caused by the necrotrophic pathogen Sclerotinia sclerotiorum and that plants grown in autoclaved soil showed almost complete elimination of pathogen resistance. Association mapping suggests at least 59 genetic locations with effects on both microbial relative abundance and Sclerotinia resistance. Although the genetic architecture appears quantitative, we have elucidated previously unexplained genetic variation for resistance to this pathogen. We identify new targets for plant breeding and demonstrate the potential for heritable microbial associations to play important roles in defence in natural and human-altered environments.

Gene expression and alternative splicing contribute to adaptive divergence of ecotypes.

Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.

The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution.

The domesticated sunflower, Helianthus annuus L., is a global oil crop that has promise for climate change adaptation, because it can maintain stable yields across a wide variety of environmental conditions, including drought. Even greater resilience is achievable through the mining of resistance alleles from compatible wild sunflower relatives, including numerous extremophile species. Here we report a high-quality reference for the sunflower genome (3.6 gigabases), together with extensive transcriptomic data from vegetative and floral organs. The genome mostly consists of highly similar, related sequences and required single-molecule real-time sequencing technologies for successful assembly. Genome analyses enabled the reconstruction of the evolutionary history of the Asterids, further establishing the existence of a whole-genome triplication at the base of the Asterids II clade and a sunflower-specific whole-genome duplication around 29 million years ago. An integrative approach combining quantitative genetics, expression and diversity data permitted development of comprehensive gene networks for two major breeding traits, flowering time and oil metabolism, and revealed new candidate genes in these networks. We found that the genomic architecture of flowering time has been shaped by the most recent whole-genome duplication, which suggests that ancient paralogues can remain in the same regulatory networks for dozens of millions of years. This genome represents a cornerstone for future research programs aiming to exploit genetic diversity to improve biotic and abiotic stress resistance and oil production, while also considering agricultural constraints and human nutritional needs.

Estimating phylogenies from genomes: A beginners review of commonly used genomic data in vertebrate phylogenomics.

Despite the increasing feasibility of sequencing whole genomes from diverse taxa, a persistent problem in phylogenomics is the selection of appropriate genetic markers or loci for a given taxonomic group or research question. In this review, we aim to streamline the decision-making process when selecting specific markers to use in phylogenomic studies by introducing commonly used types of genomic markers, their evolutionary characteristics, and their associated uses in phylogenomics. Specifically, we review the utilities of ultraconserved elements (including flanking regions), anchored hybrid enrichment loci, conserved nonexonic elements, untranslated regions, introns, exons, mitochondrial DNA, single nucleotide polymorphisms, and anonymous regions (nonspecific regions that are evenly or randomly distributed across the genome). These various genomic elements and regions differ in their substitution rates, likelihood of neutrality or of being strongly linked to loci under selection, and mode of inheritance, each of which are important considerations in phylogenomic reconstruction. These features may give each type of marker important advantages and disadvantages depending on the biological question, number of taxa sampled, evolutionary timescale, cost effectiveness, and analytical methods used. We provide a concise outline as a resource to efficiently consider key aspects of each type of genetic marker. There are many factors to consider when designing phylogenomic studies, and this review may serve as a primer when weighing options between multiple potential phylogenomic markers.

Adaptation to distinct habitats is maintained by contrasting selection at different life stages in sunflower ecotypes.

Conspecific populations living in adjacent but contrasting microenvironments represent excellent systems for studying natural selection. These systems are valuable because gene flow is expected to force genetic homogeneity except at loci experiencing divergent selection. A history of reciprocal transplant and common garden studies in such systems, and a growing number of genomic studies, have contributed to understanding how selection operates in natural populations. While selection can vary across different fitness components and life stages, few studies have investigated how this ultimately affects allele frequencies and the maintenance of divergence between populations. Here, we study two sunflower ecotypes in distinct, adjacent habitats by combining demographic models with genome-wide sequence data to estimate fitness and allele frequency change at multiple life stages. This framework allows us to estimate that only local ecotypes are likely to experience positive population growth (λ > 1) and that the maintenance of divergent adaptation appears to be mediated via habitat- and life stage-specific selection. We identify genetic variation, significantly driven by loci in chromosomal inversions, associated with different life history strategies in neighbouring ecotypes that optimize different fitness components and may contribute to the maintenance of distinct ecotypes.

Rapid evolutionary changes in gene expression in response to climate fluctuations.

There is now abundant evidence of rapid evolution in natural populations, but the genetic mechanisms of these changes remain unclear. One possible route to rapid evolution is through changes in the expression of genes that influence traits under selection. We examined contemporary evolutionary gene expression changes in plant populations responding to environmental fluctuations. We compared genome-wide gene expression, using RNA-seq, in two populations of Brassica rapa collected over four time points between 1997 and 2014, during which precipitation in southern California fluctuated dramatically and phenotypic and genotypic changes occurred. By combining transcriptome profiling with the resurrection approach, we directly examined evolutionary changes in gene expression over time. For both populations, we found a substantial number of differentially expressed genes between generations, indicating rapid evolution in the expression of many genes. Using existing gene annotations, we found that many changes occurred in genes involved in regulating stress responses and flowering time. These appeared related to the fluctuations in precipitation and were potentially adaptive. However, the evolutionary changes in gene expression differed across generations within and between populations, indicating largely independent evolutionary trajectories across populations and over time. Our study provides strong evidence for rapid evolution in gene expression, and indicates that changes in gene expression can be one mechanism of rapid evolutionary responses to selection episodes. This study also illustrates that combining resurrection studies with transcriptomics is a powerful approach for investigating evolutionary changes at the gene regulatory level, and will provide new insights into the genetic basis of contemporary evolution.

Breeding for sustainable oilseed crop yield and quality in a changing climate.

As the effects of climate change continue to alter crop-growing conditions year-to-year on both prime and marginal agricultural landscapes, we must consider the effects not only on yield but also on quality. This is particularly true for oilseed crops. In this review, we explore the importance of oilseeds in general and the specific uses of major oilseed crops including soybean, sunflower, canola, peanut, and cottonseed. We review the physiology of seed oil production, from the perspective of the plant's adaptation to environmental changes. Of particular importance is the role of temperature and water availability on oil synthesis. We then discuss how this influences genetic variation, phenotype variability due to environment, and the interaction of genetics and environment to affect composition and yield of vegetable oils. The ability to predict these effects using genomics and bioinformatics is an important new frontier for breeders to maximize stability of a desired fatty acid composition for their crop over increasingly extreme agricultural environments.

Genetic loci underlying quantitative resistance to necrotrophic pathogens Sclerotinia and Diaporthe (Phomopsis), and correlated resistance to both pathogens.

KEY MESSAGE: We provide results rooted in quantitative genetics, which combined with knowledge of candidate gene function, helps us to better understand the resistance to two major necrotrophic pathogens of sunflower. Necrotrophic pathogens can avoid or even benefit from plant defenses used against biotrophic pathogens, and thus represent a distinct challenge to plant populations in natural and agricultural systems. Sclerotinia and Phomopsis/Diaporthe are detrimental pathogens for many dicotyledonous plants, including many economically important plants. With no well-established methods to prevent infection in susceptible plants, host-plant resistance is currently the most effective strategy. Despite knowledge of a moderate, positive correlation in resistance to the two diseases in sunflower, detailed analysis of the genetics, in the same populations, has not been conducted. We present results of genome-wide analysis of resistance to both pathogens in a diversity panel of 218 domesticated sunflower genotypes of worldwide origin. We identified 14 Sclerotinia head rot and 7 Phomopsis stem canker unique QTLs, plus 1 co-located QTL for both traits, and observed extensive patterns of linkage disequilibrium between sites for both traits. Most QTLs contained one credible candidate gene, and gene families were common for the two disease resistance traits. These results suggest there has been strong, simultaneous selection for resistance to these two diseases and that a generalized mechanism for defense against these necrotrophic pathogens exists.

Genetics of alternative splicing evolution during sunflower domestication.

Alternative splicing enables organisms to produce the diversity of proteins necessary for multicellular life by using relatively few protein-coding genes. Although differences in splicing have been identified among divergent taxa, the shorter-term evolution of splicing is understudied. The origins of novel splice forms, and the contributions of alternative splicing to major evolutionary transitions, are largely unknown. This study used transcriptomes of wild and domesticated sunflowers to examine splice differentiation and regulation during domestication. We identified substantial splicing divergence between wild and domesticated sunflowers, mainly in the form of intron retention. Transcripts with divergent splicing were enriched for seed-development functions, suggesting that artificial selection impacted splicing patterns. Mapping of quantitative trait loci (QTLs) associated with 144 differential splicing cases revealed primarily trans-acting variation affecting splicing patterns. A large proportion of identified QTLs contain known spliceosome proteins and are associated with splicing variation in multiple genes. Examining a broader set of wild and domesticated sunflower genotypes revealed that most differential splicing patterns in domesticated sunflowers likely arose from standing variation in wild Helianthus annuus and gained frequency during the domestication process. However, several domesticate-associated splicing patterns appear to be introgressed from other Helianthus species. These results suggest that sunflower domestication involved selection on pleiotropic regulatory alleles. More generally, our findings indicate that substantial differences in isoform abundances arose rapidly during a recent evolutionary transition and appear to contribute to adaptation and population divergence.

Stochastic processes drive rapid genomic divergence during experimental range expansions.

Range expansions are crucibles for rapid evolution, acting via both selective and neutral mechanisms. While selection on traits such as dispersal and fecundity may increase expansion speed, neutral mechanisms arising from repeated bottlenecks and genetic drift in edge populations (i.e. gene surfing) could slow spread or make it less predictable. Thus, it is necessary to disentangle the effects of selection from neutral mechanisms to robustly predict expansion dynamics. This is difficult to do with expansions in nature, as replicated expansions are required to distinguish selective and neutral processes in the genome. Using replicated microcosms of the red flour beetle ( Tribolium castaneum), we identify a robust signature of stochastic, neutral mechanisms in genomic changes arising over only eight generations of expansion and assess the role of standing variation and de novo mutations in driving these changes. Average genetic diversity was reduced within edge populations, but with substantial among-replicate variability in the changes at specific genomic windows. Such variability in genomic changes is consistent with a large role for stochastic, neutral processes. This increased genomic divergence among populations was mirrored by heightened variation in population size and expansion speed, suggesting that stochastic variation in the genome could increase unpredictability of range expansions.

A taxonomically broad metagenomic survey of 339 species spanning 57 families suggests cystobasidiomycete yeasts are not ubiquitous across all lichens.

PREMISE: Lichens are fungi that enter into obligate symbioses with photosynthesizing organisms (algae, cyanobacteria). Traditional narratives of lichens as binary symbiont pairs have given way to their recognition as dynamic metacommunities. Basidiomycete yeasts, particularly of the genus Cyphobasidium, have been inferred to be widespread and important components of lichen metacommunities. Yet, the presence of basidiomycete yeasts across a wide diversity of lichen lineages has not previously been tested. METHODS: We searched for lichen-associated cystobasidiomycete yeasts in newly generated metagenomic data from 413 samples of 339 lichen species spanning 57 families and 25 orders. The data set was generated as part of a large-scale project to study lichen biodiversity gradients in the southern Appalachian Mountains Biodiversity Hotspot of southeastern North America. RESULTS: Our efforts detected cystobasidiomycete yeasts in nine taxa (Bryoria nadvornikiana, Heterodermia leucomelos, Lecidea roseotincta, Opegrapha vulgata, Parmotrema hypotropum, P. subsumptum, Usnea cornuta, U. strigosa, and U. subgracilis), representing 2.7% of all species sampled. Seven of these taxa (78%) are foliose (leaf-like) or fruticose (shrubby) lichens that belong to families where basidiomycete yeasts have been previously detected. In several of the nine cases, cystobasidiomycete rDNA coverage was comparable to, or greater than, that of the primary lichen fungus single-copy nuclear genomic rDNA, suggesting sampling artifacts are unlikely to account for our results. CONCLUSIONS: Studies from diverse areas of the natural sciences have led to the need to reconceptualize lichens as dynamic metacommunities. However, our failure to detect cystobasidiomycetes in 97.3% (330 species) of the sampled species suggests that basidiomycete yeasts are not ubiquitous in lichens.

Diversity and evolution of the repetitive genomic content in Cannabis sativa.

BACKGROUND: The repetitive content of the genome, once considered to be "junk DNA", is in fact an essential component of genomic architecture and evolution. In this study, we used the genomes of three varieties of Cannabis sativa, three varieties of Humulus lupulus and one genotype of Morus notabilis to explore their repetitive content using a graph-based clustering method, designed to explore and compare repeat content in genomes that have not been fully assembled. RESULTS: The repetitive content in the C. sativa genome is mainly composed of the retrotransposons LTR/Copia and LTR/Gypsy (14% and 14.8%, respectively), ribosomal DNA (2%), and low-complexity sequences (29%). We observed a recent copy number expansion in some transposable element families. Simple repeats and low complexity regions of the genome show higher intra and inter species variation. CONCLUSIONS: As with other sequenced genomes, the repetitive content of C. sativa's genome exhibits a wide range of evolutionary patterns. Some repeat types have patterns of diversity consistent with expansions followed by losses in copy number, while others may have expanded more slowly and reached a steady state. Still, other repetitive sequences, particularly ribosomal DNA (rDNA), show signs of concerted evolution playing a major role in homogenizing sequence variation.

Reductions in complexity of mitochondrial genomes in lichen-forming fungi shed light on genome architecture of obligate symbioses.

Symbioses among co-evolving taxa are often marked by genome reductions such as a loss of protein-coding genes in at least one of the partners as a means of reducing redundancy or intergenomic conflict. To explore this phenomenon in an iconic yet under-studied group of obligate symbiotic organisms, mitochondrial genomes of 22 newly sequenced and annotated species of lichenized fungi were compared to 167 mitochondrial genomes of nonlichenized fungi. Our results demonstrate the first broad-scale loss of atp9 from mitochondria of lichenized fungi. Despite key functions in mitochondrial energy production, we show that atp9 has been independently lost in three different lineages spanning 10 of the 22 studied species. A search for predicted, functional copies of atp9 among genomes of other symbionts involved in each lichen revealed the full-length, presumably functional copies of atp9 in either the photosynthetic algal partner or in other symbiotic fungi in all 10 instances. Together, these data yield evidence of an obligate symbiotic relationship in which core genomic processes have been streamlined, likely due to co-evolution.

Demographic inference in barn swallows using whole-genome data shows signal for bottleneck and subspecies differentiation during the Holocene.

Accounting for historical demographic features is vital for many types of evolutionary inferences, including the estimation of divergence times between closely related populations. In barn swallow, Hirundo rustica, inferring historical population sizes and subspecies divergence times can shed light on the recent co-evolution of this species with humans. Pairwise sequentially Markovian coalescent uncovered population growth beginning on the order of one million years ago-which may reflect the radiation of the broader Hirundo genus-and a more recent population decline. Additionally, we used approximate Bayesian computation to evaluate hypotheses about recent timescale barn swallow demography, including population growth due to human commensalism, and a potential founder event associated with the onset of nesting on human structures. We found signal for a bottleneck event approximately 7,700 years ago, near the time that humans began building substantial structures, although there was considerable uncertainty associated with this estimate. Subspecies differentiation and subsequent growth occurred after the bottleneck in the best-supported model, an order of magnitude more recently than previous estimates in this system. We also compared results obtained from whole-genome sequencing versus reduced representation sequencing, finding many similar results despite substantial allelic dropout in the reduced representation data, which may have affected estimates of some parameters. This study presents the first genetic evidence of a potential barn swallow founder effect and subspecies divergence coinciding with the Holocene, which is an important step in analysing the biogeographical history of a well-known human commensal species.

Plant domestication and the assembly of bacterial and fungal communities associated with strains of the common sunflower, Helianthus annuus.

Root and rhizosphere microbial communities can affect plant health, but it remains undetermined how plant domestication may influence these bacterial and fungal communities. We grew 33 sunflower (Helianthus annuus) strains (n = 5) that varied in their extent of domestication and assessed rhizosphere and root endosphere bacterial and fungal communities. We also assessed fungal communities in the sunflower seeds to investigate the degree to which root and rhizosphere communities were influenced by vertical transmission of the microbiome through seeds. Neither root nor rhizosphere bacterial communities were affected by the extent of sunflower domestication, but domestication did affect the composition of rhizosphere fungal communities. In particular, more modern sunflower strains had lower relative abundances of putative fungal pathogens. Seed-associated fungal communities strongly differed across strains, but several lines of evidence suggest that there is minimal vertical transmission of fungi from seeds to the adult plants. Our results indicate that plant-associated fungal communities are more strongly influenced by host genetic factors and plant breeding than bacterial communities, a finding that could influence strategies for optimizing microbial communities to improve crop yields.

Genomic variation across two barn swallow hybrid zones reveals traits associated with divergence in sympatry and allopatry.

Hybrid zones are geographic regions where isolating barriers between divergent populations are challenged by admixture. Identifying factors that facilitate or inhibit hybridization in sympatry can illuminate the processes that maintain those reproductive barriers. We analysed patterns of hybridization and phenotypic variation across two newly discovered hybrid zones between three subspecies of barn swallow (Hirundo rustica). These subspecies differ in ventral coloration and wing length, traits that are targets of sexual and natural selection, respectively, and are associated with genome-wide differentiation in allopatry. We tested the hypothesis that the degree of divergence in these traits is associated with the extent of hybridization in secondary contact. We applied measures of population structure based on >23,000 SNPs to confirm that named subspecies correspond to distinct genomic clusters, and assessed coincidence between geographic clines for ancestry and phenotype. Although gene flow was ongoing across both hybrid zones and pairwise FST between subspecies was extremely low, we found striking differences in the extent of hybridization. In the more phenotypically differentiated subspecies pair, clines for ancestry, wing length and ventral coloration were steep and coincident, suggestive of strong isolation and, potentially, selection associated with phenotype. In the less phenotypically differentiated pair, gene flow and phenotypic variation occurred over a wide geographic span, indicative of weaker isolation. Traits associated with genome-wide differentiation in allopatry may thus also contribute to isolation in sympatry. We discuss potentially important additional roles for evolutionary history and ecology in shaping variation in the extent hybridization between closely related pairs of subspecies.

Effects of assortative mate choice on the genomic and morphological structure of a hybrid zone between two bird subspecies.

Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between the alba and personata subspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437 SNP loci on a 3,000 km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome-wide divergence. Variation in only one trait-head plumage patterning-was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100 km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome-wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.

Rapid genome-wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools.

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome-wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late-season drought in California. These ancestor-descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (FST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome-wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.