Complex Polyploids: Origins, Genomic Composition, and Role of Introgressed Alleles

Introgression allows polyploid species to acquire new genomic content from diploid progenitors or from other unrelated diploid or polyploid lineages, contributing to genetic diversity and facilitating adaptive allele discovery. In some cases, high levels of introgression elicit the replacement of large numbers of alleles inherited from the polyploid's ancestral species, profoundly reshaping the polyploid's genomic composition. In such complex polyploids it is often difficult to determine which taxa were the progenitor species and which taxa provided additional introgressive blocks through subsequent hybridization. Here, we use population-level genomic data to reconstruct the phylogenetic history of Betula pubescens (downy birch), a tetraploid species often assumed to be of allopolyploid origin and which is known to hybridize with at least four other birch species. This was achieved by modeling of polyploidization and introgression events under the multispecies coalescent and then using an approximate Bayesian computation (ABC) rejection algorithm to evaluate and compare competing polyploidization models. We provide evidence that B. pubescens is the outcome of an autoploid genome doubling event in the common ancestor of B. pendula and its extant sister species, B. platyphylla, that took place approximately 178,000-188,000 generations ago. Extensive hybridization with B. pendula, B. nana, and B. humilis followed in the aftermath of autopolyploidization, with the relative contribution of each of these species to the B. pubescens genome varying markedly across the species' range. Functional analysis of B. pubescens loci containing alleles introgressed from B. nana identified multiple genes involved in climate adaptation, while loci containing alleles derived from B. humilis revealed several genes involved in the regulation of meiotic stability and pollen viability in plant species.

Early-acting inbreeding depression can evolve as an inbreeding avoidance mechanism

Despite the potential for mechanical, developmental and/or chemical mechanisms to prevent self-fertilization, incidental self-fertilization is inevitable in many predominantly outcrossing species. In such cases, inbreeding can compromise individual fitness. Unquestionably, much of this inbreeding depression is maladaptive. However, we show that when reproductive compensation allows for the replacement of inviable embryos lost early in development, selection can favour deleterious recessive variants that induce 'self-sacrificial' death of inbred embryos. Our theoretical results provide numerous testable predictions which could challenge the assumption that inbreeding depression is always maladaptive. Our work is applicable any species that cannot fully avoid inbreeding, exhibits substantial inbreeding depression, and has the potential to compensate embryos lost early in development. In addition to its general applicability, our theory suggests that self-sacrificial variants might be responsible for the remarkably low realized selfing rates of gymnosperms with high primary selfing rates, as gymnosperms exhibit strong inbreeding depression, have effective reproductive compensation mechanisms, and cannot evolve chemical self-incompatibility.

Between but Not Within-Species Variation in the Distribution of Fitness Effects

New mutations provide the raw material for evolution and adaptation. The distribution of fitness effects (DFE) describes the spectrum of effects of new mutations that can occur along a genome, and is, therefore, of vital interest in evolutionary biology. Recent work has uncovered striking similarities in the DFE between closely related species, prompting us to ask whether there is variation in the DFE among populations of the same species, or among species with different degrees of divergence, that is whether there is variation in the DFE at different levels of evolution. Using exome capture data from six tree species sampled across Europe we characterized the DFE for multiple species, and for each species, multiple populations, and investigated the factors potentially influencing the DFE, such as demography, population divergence, and genetic background. We find statistical support for the presence of variation in the DFE at the species level, even among relatively closely related species. However, we find very little difference at the population level, suggesting that differences in the DFE are primarily driven by deep features of species biology, and those evolutionarily recent events, such as demographic changes and local adaptation, have little impact.

BELMM: Bayesian model selection and random walk smoothing in time-series clustering

MOTIVATION

Due to advances in measuring technology, many new phenotype, gene expression, and other omics time-course datasets are now commonly available. Cluster analysis may provide useful information about the structure of such data.

RESULTS

In this work, we propose BELMM (Bayesian Estimation of Latent Mixture Models): a flexible framework for analysing, clustering, and modelling time-series data in a Bayesian setting. The framework is built on mixture modelling: first, the mean curves of the mixture components are assumed to follow random walk smoothing priors. Second, we choose the most plausible model and the number of mixture components using the Reversible-jump Markov chain Monte Carlo. Last, we assign the individual time series into clusters based on the similarity to the cluster-specific trend curves determined by the latent random walk processes. We demonstrate the use of fast and slow implementations of our approach on both simulated and real time-series data using widely available software R, Stan, and CU-MSDSp.

AVAILABILITY AND IMPLEMENTATION

The French mortality dataset is available at http://www.mortality.org, the Drosophila melanogaster embryogenesis gene expression data at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE121160. Details on our simulated datasets are available in the Supplementary Material, and R scripts and a detailed tutorial on GitHub at https://github.com/ollisa/BELMM. The software CU-MSDSp is available on GitHub at https://github.com/jtchavisIII/CU-MSDSp.

Strong Purifying Selection in Haploid Tissue-Specific Genes of Scots Pine Supports the Masking Theory

The masking theory states that genes expressed in a haploid stage will be under more efficient selection. In contrast, selection will be less efficient in genes expressed in a diploid stage, where the fitness effects of recessive deleterious or beneficial mutations can be hidden from selection in heterozygous form. This difference can influence several evolutionary processes such as the maintenance of genetic variation, adaptation rate, and genetic load. Masking theory expectations have been confirmed in single-cell haploid and diploid organisms. However, in multicellular organisms, such as plants, the effects of haploid selection are not clear-cut. In plants, the great majority of studies indicating haploid selection have been carried out using male haploid tissues in angiosperms. Hence, evidence in these systems is confounded with the effects of sexual selection and intraspecific competition. Evidence from other plant groups is scarce, and results show no support for the masking theory. Here, we have used a gymnosperm Scots pine megagametophyte, a maternally derived seed haploid tissue, and four diploid tissues to test the strength of purifying selection on a set of genes with tissue-specific expression. By using targeted resequencing data of those genes, we obtained estimates of genetic diversity, the site frequency spectrum of 0-fold and 4-fold sites, and inferred the distribution of fitness effects of new mutations in haploid and diploid tissue-specific genes. Our results show that purifying selection is stronger for tissue-specific genes expressed in the haploid megagametophyte tissue and that this signal of strong selection is not an artifact driven by high expression levels.

Reproductive compensation and selection among viable embryos drive the evolution of polyembryony

AbstractSimple polyembryony, where one gametophyte produces multiple embryos with different sires but the same maternal haplotype, is common among vascular plants. We develop an infinite-sites, forward population genetics model showing that together polyembryony's two benefits-"reproductive compensation" achieved by providing a backup for inviable embryos and the opportunity to favor the fitter of surviving embryos-can favor its evolution. Our model tests how these factors can favor the evolution of polyembryony and how these underlying benefits of polyembryony shape the genetic load under a range of biological parameters. While these two benefits are difficult to disentangle in nature, we construct variant models of polyembryony that either only include or only exclude the opportunity for reproductive compensation. We find that reproductive compensation strongly favors the evolution of polyembryony and that polyembryony is favored much more weekly in its absence, suggesting that the benefit of a backup embryo is the major force favoring polyembryony. Remarkably, we find nearly identical results in cases in which mutations impact either embryonic or postembryonic fitness (no pleiotropy) and in cases in which mutations have identical effects on embryonic and postembryonic fitness (extreme pleiotropy). Finally, we find that the consequences of polyembryony depend on its function-polyembryony results in a decrease in mean embryonic fitness when acting as a mechanism of embryo compensation and ultimately increases mean embryonic fitness when we exclude this potential benefit.

Taming the massive genome of Scots pine with PiSy50k, a new genotyping array for conifer research

Pinus sylvestris (Scots pine) is the most widespread coniferous tree in the boreal forests of Eurasia, with major economic and ecological importance. However, its large and repetitive genome presents a challenge for conducting genome-wide analyses such as association studies, genetic mapping and genomic selection. We present a new 50K single-nucleotide polymorphism (SNP) genotyping array for Scots pine research, breeding and other applications. To select the SNP set, we first genotyped 480 Scots pine samples on a 407 540 SNP screening array and identified 47 712 high-quality SNPs for the final array (called 'PiSy50k'). Here, we provide details of the design and testing, as well as allele frequency estimates from the discovery panel, functional annotation, tissue-specific expression patterns and expression level information for the SNPs or corresponding genes, when available. We validated the performance of the PiSy50k array using samples from Finland and Scotland. Overall, 39 678 (83.2%) SNPs showed low error rates (mean = 0.9%). Relatedness estimates based on array genotypes were consistent with the expected pedigrees, and the level of Mendelian error was negligible. In addition, array genotypes successfully discriminate between Scots pine populations of Finnish and Scottish origins. The PiSy50k SNP array will be a valuable tool for a wide variety of future genetic studies and forestry applications.

Genomics and adaptation in forest ecosystems

Rapid human-induced environmental changes like climate warming represent a challenge for forest ecosystems. Due to their biological complexity and the long generation time of their keystone tree species, genetic adaptation in these ecosystems might not be fast enough to keep track with conditions changing at such a fast pace. The study of adaptation to environmental change and its genetic mechanisms is therefore key for ensuring a sustainable support and management of forests. The 4-day conference of the European Research Group EvolTree (https://www.evoltree.eu) on the topic of "Genomics and Adaptation in Forest Ecosystems" brought together over 130 scientists to present and discuss the latest developments and findings in forest evolutionary research. Genomic studies in forest trees have long been hampered by the lack of high-quality genomics resources and affordable genotyping methods. This has dramatically changed in the last few years; the conference impressively showed how such tools are now being applied to study past demography, adaptation and interactions with associated organisms. Moreover, genomic studies are now finally also entering the world of conservation and forest management, for example by measuring the value or cost of interspecific hybridization and introgression, assessing the vulnerability of species and populations to future change, or accurately delineating evolutionary significant units. The newly launched conference series of EvolTree will hopefully play a key role in the exchange and synthesis of such important investigations.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11295-022-01542-1.

Adaptive responses to temperature and precipitation variation at the early-life stages of Pinus sylvestris

Early-stage fitness variation has been seldom evaluated at broad scales in forest tree species, despite the long tradition of studying climate-driven intraspecific genetic variation. In this study, we evaluated the role of climate in driving patterns of population differentiation at early-life stages in Pinus sylvestris and explored the fitness and growth consequences of seed transfer within the species range. We monitored seedling emergence, survival and growth over a 2-yr period in a multi-site common garden experiment which included 18 European populations and spanned 25° in latitude and 1700 m in elevation. Climate-fitness functions showed that populations exhibited higher seedling survival and growth at temperatures similar to their home environment, which is consistent with local adaptation. Northern populations experienced lower survival and growth at warmer sites, contrary to previous studies on later life stages. Seed mass was higher in populations from warmer areas and was positively associated with survival and growth at more southern sites. Finally, we did not detect a survival-growth trade-off; on the contrary, bigger seedlings exhibited higher survival probabilities under most climatic conditions. In conclusion, our results reveal that contrasting temperature regimes have played an important role in driving the divergent evolution of P. sylvestris populations at early-life stages.

Atlas of tissue-specific and tissue-preferential gene expression in ecologically and economically significant conifer Pinus sylvestris

Despite their ecological and economical importance, conifers genomic resources are limited, mainly due to the large size and complexity of their genomes. Additionally, the available genomic resources lack complete structural and functional annotation. Transcriptomic resources have been commonly used to compensate for these deficiencies, though for most conifer species they are limited to a small number of tissues, or capture only a fraction of the genes present in the genome. Here we provide an atlas of gene expression patterns for conifer Pinus sylvestris across five tissues: embryo, megagametophyte, needle, phloem and vegetative bud. We used a wide range of tissues and focused our analyses on the expression profiles of genes at tissue level. We provide comprehensive information of the per-tissue normalized expression level, indication of tissue preferential upregulation and tissue-specificity of expression. We identified a total of 48,001 tissue preferentially upregulated and tissue specifically expressed genes, of which 28% have annotation in the Swiss-Prot database. Even though most of the putative genes identified do not have functional information in current biological databases, the tissue-specific patterns discovered provide valuable information about their potential functions for further studies, as for example in the areas of plant physiology, population genetics and genomics in general. As we provide information on tissue specificity at both diploid and haploid life stages, our data will also contribute to the understanding of evolutionary rates of different tissue types and ploidy levels.

Model guided trait-specific co-expression network estimation as a new perspective for identifying molecular interactions and pathways

A wide variety of 1) parametric regression models and 2) co-expression networks have been developed for finding gene-by-gene interactions underlying complex traits from expression data. While both methodological schemes have their own well-known benefits, little is known about their synergistic potential. Our study introduces their methodological fusion that cross-exploits the strengths of individual approaches via a built-in information-sharing mechanism. This fusion is theoretically based on certain trait-conditioned dependency patterns between two genes depending on their role in the underlying parametric model. Resulting trait-specific co-expression network estimation method 1) serves to enhance the interpretation of biological networks in a parametric sense, and 2) exploits the underlying parametric model itself in the estimation process. To also account for the substantial amount of intrinsic noise and collinearities, often entailed by expression data, a tailored co-expression measure is introduced along with this framework to alleviate related computational problems. A remarkable advance over the reference methods in simulated scenarios substantiate the method's high-efficiency. As proof-of-concept, this synergistic approach is successfully applied in survival analysis, with acute myeloid leukemia data, further highlighting the framework's versatility and broad practical relevance.

The GenTree Platform: growth traits and tree-level environmental data in 12 European forest tree species

BACKGROUND

Progress in the field of evolutionary forest ecology has been hampered by the huge challenge of phenotyping trees across their ranges in their natural environments, and the limitation in high-resolution environmental information.

FINDINGS

The GenTree Platform contains phenotypic and environmental data from 4,959 trees from 12 ecologically and economically important European forest tree species: Abies alba Mill. (silver fir), Betula pendula Roth. (silver birch), Fagus sylvatica L. (European beech), Picea abies (L.) H. Karst (Norway spruce), Pinus cembra L. (Swiss stone pine), Pinus halepensis Mill. (Aleppo pine), Pinus nigra Arnold (European black pine), Pinus pinaster Aiton (maritime pine), Pinus sylvestris L. (Scots pine), Populus nigra L. (European black poplar), Taxus baccata L. (English yew), and Quercus petraea (Matt.) Liebl. (sessile oak). Phenotypic (height, diameter at breast height, crown size, bark thickness, biomass, straightness, forking, branch angle, fructification), regeneration, environmental in situ measurements (soil depth, vegetation cover, competition indices), and environmental modeling data extracted by using bilinear interpolation accounting for surrounding conditions of each tree (precipitation, temperature, insolation, drought indices) were obtained from trees in 194 sites covering the species' geographic ranges and reflecting local environmental gradients.

CONCLUSION

The GenTree Platform is a new resource for investigating ecological and evolutionary processes in forest trees. The coherent phenotyping and environmental characterization across 12 species in their European ranges allow for a wide range of analyses from forest ecologists, conservationists, and macro-ecologists. Also, the data here presented can be linked to the GenTree Dendroecological collection, the GenTree Leaf Trait collection, and the GenTree Genomic collection presented elsewhere, which together build the largest evolutionary forest ecology data collection available.

Selection patterns on early-life phenotypic traits in Pinus sylvestris are associated with precipitation and temperature along a climatic gradient in Europe

Understanding the dynamics of selection is key to predicting the response of tree species to new environmental conditions in the current context of climate change. However, selection patterns acting on early recruitment stages and their climatic drivers remain largely unknown in most tree species, despite being a critical period of their life cycle. We measured phenotypic selection on Pinus sylvestris seed mass, emergence time and early growth rate over 2 yr in four common garden experiments established along the latitudinal gradient of the species in Europe. Significant phenotypic plasticity and among-population genetic variation were found for all measured phenotypic traits. Heat and drought negatively affected fitness in the southern sites, but heavy rainfalls also decreased early survival in middle latitudes. Climate-driven directional selection was found for higher seed mass and earlier emergence time, while the form of selection on seedling growth rates differed among sites and populations. Evidence of adaptive and maladaptive phenotypic plasticity was found for emergence time and early growth rate, respectively. Seed mass, emergence time and early growth rate have an adaptive role in the early stages of P. sylvestris and climate strongly influences the patterns of selection on these fitness-related traits.

New model species for arctic-alpine plant molecular ecology

Arctic and alpine, high latitude and high elevation environments are one of the most stressful environments for species to inhabit. This harshness manifests itself in lower species richness in comparison to more southern vegetation zones (Francis & Currie, 2003). Furthermore, the climatic oscillations-past and predicted-have the most dramatic effect on these ecosystems. For example, in regions of continental ice sheets-the northernmost part of Western Europe and North America-the Arctic species assemblages are no older than a few thousands of years, which is a relatively short period from an evolutionary perspective. Although similar environments may have existed further south during the Ice Age, allowing some preadaptation for the Arctic species, the current habitat is a unique combination of environmental factors such as the climate, soil, bedrock, and photoperiod. Hence, understanding the evolutionary forces shaping Arctic-alpine species will be important for predicting these vulnerable environments' population viability and adaptive potential in the future. In this issue of Molecular Ecology Resources, Nowak et al. (Molecular Ecology Resources) present extensive genome-wide resources for an Arctic-alpine plant Draba nivalis. This adds a valuable new member into the cabbage family models for evolutionary genetics and adaptation studies, to accompany e.g., Arabidopsis (Nature Genetics, 43, 476; Nature, 408, 796), Arabis (Nature Plants, 1, 14023) and Capsella (Nature Genetics, 45, 831). A whole new avenue will open up for molecular ecological studies not only for D. nivalis, but the whole large Draba genus with its diverse ecological and evolutionary characteristics.

The GenTree Dendroecological Collection, tree-ring and wood density data from seven tree species across Europe

The dataset presented here was collected by the GenTree project (EU-Horizon 2020), which aims to improve the use of forest genetic resources across Europe by better understanding how trees adapt to their local environment. This dataset of individual tree-core characteristics including ring-width series and whole-core wood density was collected for seven ecologically and economically important European tree species: silver birch (Betula pendula), European beech (Fagus sylvatica), Norway spruce (Picea abies), European black poplar (Populus nigra), maritime pine (Pinus pinaster), Scots pine (Pinus sylvestris), and sessile oak (Quercus petraea). Tree-ring width measurements were obtained from 3600 trees in 142 populations and whole-core wood density was measured for 3098 trees in 125 populations. This dataset covers most of the geographical and climatic range occupied by the selected species. The potential use of it will be highly valuable for assessing ecological and evolutionary responses to environmental conditions as well as for model development and parameterization, to predict adaptability under climate change scenarios.

275 years of forestry meets genomics in Pinus sylvestris

Pinus sylvestris has a long history of basic and applied research that is relevant for both forestry and evolutionary studies. Its patterns of adaptive variation and role in forest economic and ecological systems have been studied extensively for nearly 275 years, detailed demography for a 100 years and mating system more than 50 years. However, its reference genome sequence is not yet available and genomic studies have been lagging compared to, for example, Pinus taeda and Picea abies, two other economically important conifers. Despite the lack of reference genome, many modern genomic methods are applicable for a more detailed look at its biological characteristics. For example, RNA-seq has revealed a complex transcriptional landscape and targeted DNA sequencing displays an excess of rare variants and geographically homogenously distributed molecular genetic diversity. Current DNA and RNA resources can be used as a reference for gene expression studies, SNP discovery, and further targeted sequencing. In the future, specific consequences of the large genome size, such as functional effects of regulatory open chromatin regions and transposable elements, should be investigated more carefully. For forest breeding and long-term management purposes, genomic data can help in assessing the genetic basis of inbreeding depression and the application of genomic tools for genomic prediction and relatedness estimates. Given the challenges of breeding (long generation time, no easy vegetative propagation) and the economic importance, application of genomic tools has a potential to have a considerable impact. Here, we explore how genomic characteristics of P. sylvestris, such as rare alleles and the low extent of linkage disequilibrium, impact the applicability and power of the tools.

Genome-Wide Analysis of Colonization History and Concomitant Selection in Arabidopsis lyrata

The high climatic variability in the past hundred thousand years has affected the demographic and adaptive processes in many species, especially in boreal and temperate regions undergoing glacial cycles. This has also influenced the patterns of genome-wide nucleotide variation, but the details of these effects are largely unknown. Here we study the patterns of genome-wide variation to infer colonization history and patterns of selection of the perennial herb species Arabidopsis lyrata, in locally adapted populations from different parts of its distribution range (Germany, UK, Norway, Sweden, and USA) representing different environmental conditions. Using site frequency spectra based demographic modeling, we found strong reduction in the effective population size of the species in general within the past 100,000 years, with more pronounced effects in the colonizing populations. We further found that the northwestern European A. lyrata populations (UK and Scandinavian) are more closely related to each other than with the Central European populations, and coalescent based population split modeling suggests that western European and Scandinavian populations became isolated relatively recently after the glacial retreat. We also highlighted loci showing evidence for local selection associated with the Scandinavian colonization. The results presented here give new insights into postglacial Scandinavian colonization history and its genome-wide effects.

No evidence of prenatal diversifying selection at locus or supertype levels in the dog MHC class II loci

BACKGROUND

Despite decades of studying, the mechanisms maintaining high diversity in the genes of the Major Histocompatibility Complex (MHC) are still puzzling scientists. In addition to pathogen recognition and other functions, MHC molecules may act prenatally in mate choice and in maternal-foetal interactions. These interactions are potential selective mechanisms that increase genetic diversity in the MHC. During pregnancy, immune response has a dual role: the foetus represents foreign tissue compared to mother, but histo-incompatibility is required for successful pregnancy. We have studied the prenatal selection in MHC class II loci (DLA-DQA1, DLA-DQB1 and DLA-DRB1) in domestic dogs by comparing the observed and expected offspring genotype proportions in 110 dog families. Several potential selection targets were addressed, including the peptide-binding site, the MHC locus, three-locus haplotype and supertype levels. For the supertype analysis, the first canine supertype classification was created based on in silico analysis of peptide-binding amino-acid polymorphism.

RESULTS

In most loci and levels, no deviation from the expected genotype frequencies was observed. However, one peptide-binding site in DLA-DRB1 had an excess of heterozygotes among the offspring. In addition, if the father shared a DLA-DRB1 allele with the mother, that allele was inherited by the offspring more frequently than expected, suggesting the selective advantage of a histo-compatible foetus, in contrast to our expectations.

CONCLUSIONS

We conclude that there is some evidence of post-copulatory selection at nucleotide site level in the MHC loci of pet dogs. But due to no indication of selection at locus, three-locus, or supertype levels, we estimated that the prenatal selection coefficient is less than 0.3 in domestic dogs and very likely other factors are more important in maintaining the genetic diversity in MHC loci.

Relationship between wild greylag and European domestic geese based on mitochondrial DNA

The origins of the European domestic goose are uncertain. The available information comes from archaeological findings and historical literature, but genetic evidence has hitherto been scarce. The domestic goose in Europe is derived from the greylag goose (Anser anser), but it is not known where the initial domestication took place and which of the two subspecies of greylag goose was ancestral. We aimed to determine the amount and geographical distribution of genetic diversity in modern populations of greylag geese as well as in different breeds of the domestic goose to make inferences about goose domestication. We studied DNA sequence variation in the mitochondrial control region of greylag geese from multiple populations across Europe and western Asia as well as specimens of domestic geese representing 18 modern breeds and individuals not belonging to any recognised breed. Our results show notable differences in genetic diversity between different greylag goose populations and the presence of six mitochondrial haplogroups which show a degree of geographical partitioning. The genetic diversity of the domestic goose is low, with 84% of sampled individuals having one of two major closely related haplotypes, suggesting that modern European domestic geese may derive from a narrow genetic base. The site of domestication remains unresolved, but domestic geese in Turkey were unusually diverse, indicating the importance of further sampling in the vicinity of the eastern Mediterranean and the Near East. There appears to be past or ongoing hybridisation between greylags and domestic geese in particular areas, consistent with field observations.

Natural variation in teosinte at the domestication locus teosinte branched1 (tb1)

The teosinte branched1(tb1) gene is a major QTL controlling branching differences between maize and its wild progenitor, teosinte. The insertion of a transposable element (Hopscotch) upstream of tb1 is known to enhance the gene’s expression, causing reduced tillering in maize. Observations of the maize tb1 allele in teosinte and estimates of an insertion age of the Hopscotch that predates domestication led us to investigate its prevalence and potential role in teosinte. We assessed the prevalence of the Hopscotch element across an Americas-wide sample of 837 maize and teosinte individuals using a co-dominant PCR assay. Additionally, we calculated population genetic summaries using sequence data from a subset of individuals from four teosinte populations and collected phenotypic data using seed from a single teosinte population where Hopscotch was found segregating at high frequency. Genotyping results indicate the Hopscotch element is found in a number of teosinte populations and linkage disequilibrium near tb1 does not support recent introgression from maize. Population genetic signatures are consistent with selection on the tb1 locus, revealing a potential ecological role, but a greenhouse experiment does not detect a strong association between the Hopscotch and tillering in teosinte. Our findings suggest the role of Hopscotch differs between maize and teosinte. Future work should assess tb1 expression levels in teosinte with and without the Hopscotch and more comprehensively phenotype teosinte to assess the ecological significance of the Hopscotch insertion and, more broadly, the tb1 locus in teosinte.

Complex patterns of local adaptation in teosinte

Populations of widely distributed species encounter and must adapt to local environmental conditions. However, comprehensive characterization of the genetic basis of adaptation is demanding, requiring genome-wide genotype data, multiple sampled populations, and an understanding of population structure and potential selection pressures. Here, we used single-nucleotide polymorphism genotyping and data on numerous environmental variables to describe the genetic basis of local adaptation in 21 populations of teosinte, the wild ancestor of maize. We found complex hierarchical genetic structure created by altitude, dispersal events, and admixture among subspecies, which complicated identification of locally beneficial alleles. Patterns of linkage disequilibrium revealed four large putative inversion polymorphisms showing clinal patterns of frequency. Population differentiation and environmental correlations suggest that both inversions and intergenic polymorphisms are involved in local adaptation.

A recent local sweep at the PHYA locus in the Northern European Spiterstulen population of Arabidopsis lyrata

Northern and central European Arabidopsis lyrata ssp. petraea populations are locally adapted to prevailing climatic conditions through differences in timing of life history events. The timing of flowering and, in perennials, the timing of growth cessation influence fitness. Phytochrome A may have an important role in regulating these life history traits as it perceives changes in daylength. We asked whether PHYA has contributed to local adaptation to the northern conditions in A. l. petraea. To search for signals of directional selection at the PHYA locus, we resequenced PHYA and 9 short fragments around PHYA from a 57-kb region from a German (Plech) and a Norwegian (Spiterstulen) population and compared patterns of differentiation and diversity to a set of 19 reference loci around the genome. First, we found that the populations were highly differentiated: there were three nonsynonymous fixed differences at the PHYA locus, which was in stark contrast with the total four fixed differences in the 19 reference loci. Compatible with a sweep hypothesis, variation was almost completely removed from the 9.4-kb region around PHYA in the northern Spiterstulen population. The overall level of linkage disequilibrium (LD) was higher in Spiterstulen, but there was no LD across the PHYA locus in the population, which is also a known consequence of a selective sweep. The sweep has likely occurred after the last glacial maximum, which suggests that it has contributed to adaptation to the northern conditions.

The genomic signature of crop-wild introgression in maize

The evolutionary significance of hybridization and subsequent introgression has long been appreciated, but evaluation of the genome-wide effects of these phenomena has only recently become possible. Crop-wild study systems represent ideal opportunities to examine evolution through hybridization. For example, maize and the conspecific wild teosinte Zea mays ssp. mexicana (hereafter, mexicana) are known to hybridize in the fields of highland Mexico. Despite widespread evidence of gene flow, maize and mexicana maintain distinct morphologies and have done so in sympatry for thousands of years. Neither the genomic extent nor the evolutionary importance of introgression between these taxa is understood. In this study we assessed patterns of genome-wide introgression based on 39,029 single nucleotide polymorphisms genotyped in 189 individuals from nine sympatric maize-mexicana populations and reference allopatric populations. While portions of the maize and mexicana genomes appeared resistant to introgression (notably near known cross-incompatibility and domestication loci), we detected widespread evidence for introgression in both directions of gene flow. Through further characterization of these genomic regions and preliminary growth chamber experiments, we found evidence suggestive of the incorporation of adaptive mexicana alleles into maize during its expansion to the highlands of central Mexico. In contrast, very little evidence was found for adaptive introgression from maize to mexicana. The methods we have applied here can be replicated widely, and such analyses have the potential to greatly inform our understanding of evolution through introgressive hybridization. Crop species, due to their exceptional genomic resources and frequent histories of spread into sympatry with relatives, should be particularly influential in these studies.

Diversity and abundance of the abnormal chromosome 10 meiotic drive complex in Zea mays

Maize Abnormal chromosome 10 (Ab10) contains a classic meiotic drive system that exploits the asymmetry of meiosis to preferentially transmit itself and other chromosomes containing specialized heterochromatic regions called knobs. The structure and diversity of the Ab10 meiotic drive haplotype is poorly understood. We developed a bacterial artificial chromosome (BAC) library from an Ab10 line and used the data to develop sequence-based markers, focusing on the proximal portion of the haplotype that shows partial homology to normal chromosome 10. These molecular and additional cytological data demonstrate that two previously identified Ab10 variants (Ab10-I and Ab10-II) share a common origin. Dominant PCR markers were used with fluorescence in situ hybridization to assay 160 diverse teosinte and maize landrace populations from across the Americas, resulting in the identification of a previously unknown but prevalent form of Ab10 (Ab10-III). We find that Ab10 occurs in at least 75% of teosinte populations at a mean frequency of 15%. Ab10 was also found in 13% of the maize landraces, but does not appear to be fixed in any wild or cultivated population. Quantitative analyses suggest that the abundance and distribution of Ab10 is governed by a complex combination of intrinsic fitness effects as well as extrinsic environmental variability.

Complex patterns of local adaptation in teosinte

Populations of widely distributed species encounter and must adapt to local environmental conditions. However, comprehensive characterization of the genetic basis of adaptation is demanding, requiring genome-wide genotype data, multiple sampled populations, and an understanding of population structure and potential selection pressures. Here, we used single-nucleotide polymorphism genotyping and data on numerous environmental variables to describe the genetic basis of local adaptation in 21 populations of teosinte, the wild ancestor of maize. We found complex hierarchical genetic structure created by altitude, dispersal events, and admixture among subspecies, which complicated identification of locally beneficial alleles. Patterns of linkage disequilibrium revealed four large putative inversion polymorphisms showing clinal patterns of frequency. Population differentiation and environmental correlations suggest that both inversions and intergenic polymorphisms are involved in local adaptation.

Time scales of divergence and speciation among natural populations and subspecies of Arabidopsis lyrata (Brassicaceae)

PREMISE OF THE STUDY

Plant populations that face new environments adapt and diverge simultaneously, and both processes leave footprints in their genetic diversity. Arabidopsis lyrata is an excellent species for studying these processes. Pairs of populations and subspecies of A. lyrata represent different stages of divergence. These populations are also known to be locally adapted and display various stages of emerging reproductive isolation.

METHODS

We used nucleotide diversity data from 19 loci to estimate divergence times and levels of diversity among nine A. lyrata populations. Traditional distance-based methods and model-based clustering analysis were used to supplement pairwise coalescence-based analysis of divergence.

KEY RESULTS

Estimated divergence times varied from 130,000 generations between North American and European subspecies to 39,000 generations between central European and Scandinavian populations. In concordance with previous studies, the highest level of diversity was found in Central Europe and the lowest in North America and a diverged Russian Karhumäki population. Local adaptation among Northern and central European populations has emerged during the last 39,000 generations. Populations that are reproductively isolated by prezygotic mechanisms have been separated for a longer time period of ∼70,000 generations but still have shared nucleotide polymorphism.

CONCLUSIONS

In A. lyrata, reproductively isolated populations started to diverge ∼70,000 generations ago and more closely related, locally adapted populations have been separate lineages for ∼39,000 generations. However, based on the posterior distribution of divergence times, the processes leading to reproductive isolation and local adaptation are likely to temporally coincide.

Comparative population genomics of maize domestication and improvement

Domestication and plant breeding are ongoing 10,000-year-old evolutionary experiments that have radically altered wild species to meet human needs. Maize has undergone a particularly striking transformation. Researchers have sought for decades to identify the genes underlying maize evolution, but these efforts have been limited in scope. Here, we report a comprehensive assessment of the evolution of modern maize based on the genome-wide resequencing of 75 wild, landrace and improved maize lines. We find evidence of recovery of diversity after domestication, likely introgression from wild relatives, and evidence for stronger selection during domestication than improvement. We identify a number of genes with stronger signals of selection than those previously shown to underlie major morphological changes. Finally, through transcriptome-wide analysis of gene expression, we find evidence both consistent with removal of cis-acting variation during maize domestication and improvement and suggestive of modern breeding having increased dominance in expression while targeting highly expressed genes.

Maize HapMap2 identifies extant variation from a genome in flux

Whereas breeders have exploited diversity in maize for yield improvements, there has been limited progress in using beneficial alleles in undomesticated varieties. Characterizing standing variation in this complex genome has been challenging, with only a small fraction of it described to date. Using a population genetics scoring model, we identified 55 million SNPs in 103 lines across pre-domestication and domesticated Zea mays varieties, including a representative from the sister genus Tripsacum. We find that structural variations are pervasive in the Z. mays genome and are enriched at loci associated with important traits. By investigating the drivers of genome size variation, we find that the larger Tripsacum genome can be explained by transposable element abundance rather than an allopolyploid origin. In contrast, intraspecies genome size variation seems to be controlled by chromosomal knob content. There is tremendous overlap in key gene content in maize and Tripsacum, suggesting that adaptations from Tripsacum (for example, perennialism and frost and drought tolerance) can likely be integrated into maize.

Selection on nuclear genes in a Pinus phylogeny

In this study, we investigate natural selection in a pine phylogeny. DNA sequences from 18 nuclear genes were used to construct a very well-supported species tree including 10 pine species. This tree is in complete agreement with a previously reported supertree constructed from morphological and molecular data, but there are discrepancies with previous chloroplast phylogenies within the section Pinus. A significant difference in evolutionary rate between Picea and Pinus was found, which could potentially indicate a lower mutation rate in Picea, but other scenarios are also possible. Several approaches were used to study selection patterns in a set of 21 nuclear genes in pines and in some cases in Picea and Pseudotsuga. The overall pattern suggests efficient purifying selection resulting in low branch-specific d(n)/d(s) ratios with an average of 0.22, which is similar to other higher plants. Evidence for purifying selection was common and found on at least 55% of the branches. Evidence of positive selection at several sites was found in a phytocyanin homolog and significant differences in d(n)/d(s) among the branches in the gene tree in dehydrin 1. Several genes suitable for further phylogenetic analysis at various levels of divergence were identified.

Demographic history has influenced nucleotide diversity in European Pinus sylvestris populations

To infer the role of natural selection in shaping standing genetic diversity, it is necessary to assess the genomewide impact of demographic history on nucleotide diversity. In this study we analyzed sequence diversity of 16 nuclear loci in eight Pinus sylvestris populations. Populations were divided into four geographical groups on the basis of their current location and the geographical history of the region: northern Europe, central Europe, Spain, and Turkey. There were no among-group differences in the level of silent nucleotide diversity, which was approximately 0.005/bp in all groups. There was some evidence that linkage disequilibrium extended further in northern Europe than in central Europe: the estimates of the population recombination rate parameter, rho, were 0.0064 and 0.0294, respectively. The summary statistics of nucleotide diversity in central and northern European populations were compatible with an ancient bottleneck rather than the standard neutral model.

Genomic diversity in forest trees

Forest trees in general are out-crossing, long-lived, and at early stages of domestication. Molecular evolution at neutral sites is very slow because of the long generation times. Transferring information between closely related conifer species is facilitated by high sequence similarity. At the nucleotide level, trees have at most intermediate levels of variation relative to other plants. Importantly, in many species linkage disequilibrium within genes declines within less than 1000 bp. In contrast to the slow rate of neutral evolution, large tree populations respond rapidly to natural selection. Detecting traces of selection may be easier in tree populations than in many other species. Association studies between genotypes and phenotypes are proving to be useful tools for functional genomics.