Exploring the role of polymorphic interspecies structural variants in reproductive isolation and adaptive divergence in Eucalyptus.

Structural variations (SVs) play a significant role in speciation and adaptation in many species, yet few studies have explored the prevalence and impact of different categories of SVs. We conducted a comparative analysis of long-read assembled reference genomes of closely related Eucalyptus species to identify candidate SVs potentially influencing speciation and adaptation. Interspecies SVs can be either fixed differences or polymorphic in one or both species. To describe SV patterns, we employed short-read whole-genome sequencing on over 600 individuals of Eucalyptus melliodora and Eucalyptus sideroxylon, along with recent high-quality genome assemblies. We aligned reads and genotyped interspecies SVs predicted between species reference genomes. Our results revealed that 49,756 of 58,025 and 39,536 of 47,064 interspecies SVs could be typed with short reads in E. melliodora and E. sideroxylon, respectively. Focusing on inversions and translocations, symmetric SVs that are readily genotyped within both populations, 24 were found to be structural divergences, 2,623 structural polymorphisms, and 928 shared structural polymorphisms. We assessed the functional significance of fixed interspecies SVs by examining differences in estimated recombination rates and genetic differentiation between species, revealing a complex history of natural selection. Shared structural polymorphisms displayed enrichment of potentially adaptive genes. Understanding how different classes of genetic mutations contribute to genetic diversity and reproductive barriers is essential for understanding how organisms enhance fitness, adapt to changing environments, and diversify. Our findings reveal the prevalence of interspecies SVs and elucidate their role in genetic differentiation, adaptive evolution, and species divergence within and between populations.

Phylogenomics and the rise of the angiosperms.

Angiosperms are the cornerstone of most terrestrial ecosystems and human livelihoods1,2. A robust understanding of angiosperm evolution is required to explain their rise to ecological dominance. So far, the angiosperm tree of life has been determined primarily by means of analyses of the plastid genome3,4. Many studies have drawn on this foundational work, such as classification and first insights into angiosperm diversification since their Mesozoic origins5-7. However, the limited and biased sampling of both taxa and genomes undermines confidence in the tree and its implications. Here, we build the tree of life for almost 8,000 (about 60%) angiosperm genera using a standardized set of 353 nuclear genes8. This 15-fold increase in genus-level sampling relative to comparable nuclear studies9 provides a critical test of earlier results and brings notable change to key groups, especially in rosids, while substantiating many previously predicted relationships. Scaling this tree to time using 200 fossils, we discovered that early angiosperm evolution was characterized by high gene tree conflict and explosive diversification, giving rise to more than 80% of extant angiosperm orders. Steady diversification ensued through the remaining Mesozoic Era until rates resurged in the Cenozoic Era, concurrent with decreasing global temperatures and tightly linked with gene tree conflict. Taken together, our extensive sampling combined with advanced phylogenomic methods shows the deep history and full complexity in the evolution of a megadiverse clade.

Species-specific basecallers improve actual accuracy of nanopore sequencing in plants.

BACKGROUND: Long-read sequencing platforms offered by Oxford Nanopore Technologies (ONT) allow native DNA containing epigenetic modifications to be directly sequenced, but can be limited by lower per-base accuracies. A key step post-sequencing is basecalling, the process of converting raw electrical signals produced by the sequencing device into nucleotide sequences. This is challenging as current basecallers are primarily based on mixtures of model species for training. Here we utilise both ONT PromethION and higher accuracy PacBio Sequel II HiFi sequencing on two plants, Phebalium stellatum and Xanthorrhoea johnsonii, to train species-specific basecaller models with the aim of improving per-base accuracy. We investigate sequencing accuracies achieved by ONT basecallers and assess accuracy gains by training single-species and species-specific basecaller models. We also evaluate accuracy gains from ONT's improved flowcells (R10.4, FLO-PRO112) and sequencing kits (SQK-LSK112). For the truth dataset for both model training and accuracy assessment, we developed highly accurate, contiguous diploid reference genomes with PacBio Sequel II HiFi reads. RESULTS: Basecalling with ONT Guppy 5 and 6 super-accurate gave almost identical results, attaining read accuracies of 91.96% and 94.15%. Guppy's plant-specific model gave highly mixed results, attaining read accuracies of 91.47% and 96.18%. Species-specific basecalling models improved read accuracy, attaining 93.24% and 95.16% read accuracies. R10.4 sequencing kits also improve sequencing accuracy, attaining read accuracies of 95.46% (super-accurate) and 96.87% (species-specific). CONCLUSIONS: The use of a single mixed-species basecaller model, such as ONT Guppy super-accurate, may be reducing the accuracy of nanopore sequencing, due to conflicting genome biology within the training dataset and study species. Training of single-species and genome-specific basecaller models improves read accuracy. Studies that aim to do large-scale long-read genotyping would primarily benefit from training their own basecalling models. Such studies could use sequencing accuracy gains and improving bioinformatics tools to improve study outcomes.

Annual and perennial Medicago show signatures of parallel adaptation to climate and soil in highly conserved genes.

Human induced environmental change may require rapid adaptation of plant populations and crops, but the genomic basis of environmental adaptation remain poorly understood. We analysed polymorphic loci from the perennial crop Medicago sativa (alfalfa or lucerne) and the annual legume model species M. truncatula to search for a common set of candidate genes that might contribute to adaptation to abiotic stress in both annual and perennial Medicago species. We identified a set of candidate genes of adaptation associated with environmental gradients along the distribution of the two Medicago species. Candidate genes for each species were detected in homologous genomic linkage blocks using genome-environment (GEA) and genome-phenotype association analyses. Hundreds of GEA candidate genes were species-specific, of these, 13.4% (M. sativa) and 24% (M. truncatula) were also significantly associated with phenotypic traits. A set of 168 GEA candidates were shared by both species, which was 25.4% more than expected by chance. When combined, they explained a high proportion of variance for certain phenotypic traits associated with adaptation. Genes with highly conserved functions dominated among the shared candidates and were enriched in gene ontology terms that have shown to play a central role in drought avoidance and tolerance mechanisms by means of cellular shape modifications and other functions associated with cell homeostasis. Our results point to the existence of a molecular basis of adaptation to abiotic stress in Medicago determined by highly conserved genes and gene functions. We discuss these results in light of the recently proposed omnigenic model of complex traits.

Multiple chromosomal inversions contribute to adaptive divergence of a dune sunflower ecotype.

Both models and case studies suggest that chromosomal inversions can facilitate adaptation and speciation in the presence of gene flow by suppressing recombination between locally adapted alleles. Until recently, however, it has been laborious and time-consuming to identify and genotype inversions in natural populations. Here we apply RAD sequencing data and newly developed population genomic approaches to identify putative inversions that differentiate a sand dune ecotype of the prairie sunflower (Helianthus petiolaris) from populations found on the adjacent sand sheet. We detected seven large genomic regions that exhibit a different population structure than the rest of the genome and that vary in frequency between dune and nondune populations. These regions also show high linkage disequilibrium and high heterozygosity between, but not within, arrangements, consistent with the behaviour of large inversions, an inference subsequently validated in part by comparative genetic mapping. Genome-environment association analyses show that key environmental variables, including vegetation cover and soil nitrogen, are significantly associated with inversions. The inversions colocate with previously described "islands of differentiation," and appear to play an important role in adaptive divergence and incipient speciation within H. petiolaris.

Landscape drivers of genomic diversity and divergence in woodland Eucalyptus.

Spatial genetic patterns are influenced by numerous factors, and they can vary even among coexisting, closely related species due to differences in dispersal and selection. Eucalyptus (L'Héritier 1789; the "eucalypts") are foundation tree species that provide essential habitat and modulate ecosystem services throughout Australia. Here we present a study of landscape genomic variation in two woodland eucalypt species, using whole-genome sequencing of 388 individuals of Eucalyptus albens and Eucalyptus sideroxylon. We found exceptionally high genetic diversity (π ≈ 0.05) and low genome-wide, interspecific differentiation (FST  = 0.15) and intraspecific differentiation between localities (FST  ≈ 0.01-0.02). We found no support for strong, discrete population structure, but found substantial support for isolation by geographic distance (IBD) in both species. Using generalized dissimilarity modelling, we identified additional isolation by environment (IBE). Eucalyptus albens showed moderate IBD, and environmental variables have a small but significant amount of additional predictive power (i.e. IBE). Eucalyptus sideroxylon showed much stronger IBD and moderate IBE. These results highlight the vast adaptive potential of these species and set the stage for testing evolutionary hypotheses of interspecific adaptive differentiation across environments.

Chloroplast genome analysis of box-ironbark Eucalyptus.

Eucalyptus L'Hérit. (Myrtaceae) is a taxonomically complex and highly speciose genus that dominates much of Australia's woody vegetation. However, very little information is available about the molecular biology and chloroplast diversity of certain groups, such as Eucalyptus section Adnataria, which is found in many woodland habitats of eastern Australia. We report four new complete chloroplast genomes of Eucalyptus, including three genomes from species previously lacking any chloroplast reference sequences. Plastomes of E. albens, E. conica, E. crebra and E. melliodora assembled using a de novo approach were shown to be largely identical to each other, and similar in size and structure to previously published chloroplast genomes from Eucalyptus. A total of 132 genes (114 single-copy genes and 18 duplicated genes in the IR regions) were identified, and shown to be highly conserved in terms of gene order, content and organization. Slightly higher divergence in the intergenic spacers was identified through comparative genomic analyses. Chloroplast sequences of 35 additional individuals representing 12 species were assembled using a reference guided approach. Rates of nucleotide substitution varied among the protein coding genes, with 17 genes under possible positive selection, and 29 invariant genes. Phylogenetic analysis of either the whole reconstructed plastome sequences or the individual genes revealed extreme discordance with expected species boundaries or higher-level relationships. Plastome relationships were better predicted by geography than by nuclear DNA or taxonomic relationships, suggesting a substantial influence of gene flow over and above the effects of incomplete lineage sorting. These results provide resources for future research and valuable insights into the prevalence of interspecific gene flow among Eucalyptus species.

An evaluation of alternative explanations for widespread cytonuclear discordance in annual sunflowers (Helianthus).

Cytonuclear discordance is commonly observed in phylogenetic studies, yet few studies have tested whether these patterns reflect incomplete lineage sorting or organellar introgression. Here, we used whole-chloroplast sequence data in combination with over 1000 nuclear single-nucleotide polymorphisms to clarify the extent of cytonuclear discordance in wild annual sunflowers (Helianthus), and to test alternative explanations for such discordance. Our phylogenetic analyses indicate that cytonuclear discordance is widespread within this group, both in terms of the relationships among species and among individuals within species. Simulations of chloroplast evolution show that incomplete lineage sorting cannot explain these patterns in most cases. Instead, most of the observed discordance is better explained by cytoplasmic introgression. Molecular tests of evolution further indicate that selection may have played a role in driving patterns of plastid variation - although additional experimental work is needed to fully evaluate the importance of selection on organellar variants in different parts of the geographic range. Overall, this study represents one of the most comprehensive tests of the drivers of cytonuclear discordance and highlights the potential for gene flow to lead to extensive organellar introgression in hybridizing taxa.

Landscape genomic prediction for restoration of a Eucalyptus foundation species under climate change.

As species face rapid environmental change, we can build resilient populations through restoration projects that incorporate predicted future climates into seed sourcing decisions. Eucalyptus melliodora is a foundation species of a critically endangered community in Australia that is a target for restoration. We examined genomic and phenotypic variation to make empirical based recommendations for seed sourcing. We examined isolation by distance and isolation by environment, determining high levels of gene flow extending for 500 km and correlations with climate and soil variables. Growth experiments revealed extensive phenotypic variation both within and among sampling sites, but no site-specific differentiation in phenotypic plasticity. Model predictions suggest that seed can be sourced broadly across the landscape, providing ample diversity for adaptation to environmental change. Application of our landscape genomic model to E. melliodora restoration projects can identify genomic variation suitable for predicted future climates, thereby increasing the long term probability of successful restoration.

The K = 2 conundrum.

Assessments of population genetic structure have become an increasing focus as they can provide valuable insight into patterns of migration and gene flow. structure, the most highly cited of several clustering-based methods, was developed to provide robust estimates without the need for populations to be determined a priori. structure introduces the problem of selecting the optimal number of clusters, and as a result, the ΔK method was proposed to assist in the identification of the "true" number of clusters. In our review of 1,264 studies using structure to explore population subdivision, studies that used ΔK were more likely to identify K = 2 (54%, 443/822) than studies that did not use ΔK (21%, 82/386). A troubling finding was that very few studies performed the hierarchical analysis recommended by the authors of both ΔK and structure to fully explore population subdivision. Furthermore, extensions of earlier simulations indicate that, with a representative number of markers, ΔK frequently identifies K = 2 as the top level of hierarchical structure, even when more subpopulations are present. This review suggests that many studies may have been over- or underestimating population genetic structure; both scenarios have serious consequences, particularly with respect to conservation and management. We recommend publication standards for population structure results so that readers can assess the implications of the results given their own understanding of the species biology.

Voyage of discovery? A comment on Koch et al. "A voyage to Terra Australis: human-mediated dispersal of cats".

The origins of feral cats in Australia may be understood with the help of molecular studies, but it is important that hypotheses be tested with appropriate sampling and methodology. We point out several shortcomings in the analysis by Koch et al. (BMC Evol Biol 15:262, 2015; A voyage to Terra Australis: human-mediated dispersal of cats. Dryad Digital Repository, 2015), present a reanalysis of part of the study and discuss the challenges of elucidating the early history of feral cats.

Multiple reproductive barriers separate recently diverged sunflower ecotypes.

Measuring reproductive barriers between groups of organisms is an effective way to determine the traits and mechanisms that impede gene flow. However, to understand the ecological and evolutionary factors that drive speciation, it is important to distinguish between the barriers that arise early in the speciation process and those that arise after speciation is largely complete. In this article, we comprehensively test for reproductive isolation between recently diverged (<10,000 years bp) dune and nondune ecotypes of the prairie sunflower, Helianthus petiolaris. We find reproductive barriers acting at multiple stages of hybridization, including premating, postmating-prezygotic, and postzygotic barriers, despite the recent divergence. Barriers include extrinsic selection against immigrants and hybrids, a shift in pollinator assemblage, and postpollination assortative mating. Together, these data suggest that multiple barriers can be important for reducing gene flow in the earliest stages of speciation.

Assessing the reproducibility of discriminant function analyses.

Data are the foundation of empirical research, yet all too often the datasets underlying published papers are unavailable, incorrect, or poorly curated. This is a serious issue, because future researchers are then unable to validate published results or reuse data to explore new ideas and hypotheses. Even if data files are securely stored and accessible, they must also be accompanied by accurate labels and identifiers. To assess how often problems with metadata or data curation affect the reproducibility of published results, we attempted to reproduce Discriminant Function Analyses (DFAs) from the field of organismal biology. DFA is a commonly used statistical analysis that has changed little since its inception almost eight decades ago, and therefore provides an opportunity to test reproducibility among datasets of varying ages. Out of 100 papers we initially surveyed, fourteen were excluded because they did not present the common types of quantitative result from their DFA or gave insufficient details of their DFA. Of the remaining 86 datasets, there were 15 cases for which we were unable to confidently relate the dataset we received to the one used in the published analysis. The reasons ranged from incomprehensible or absent variable labels, the DFA being performed on an unspecified subset of the data, or the dataset we received being incomplete. We focused on reproducing three common summary statistics from DFAs: the percent variance explained, the percentage correctly assigned and the largest discriminant function coefficient. The reproducibility of the first two was fairly high (20 of 26, and 44 of 60 datasets, respectively), whereas our success rate with the discriminant function coefficients was lower (15 of 26 datasets). When considering all three summary statistics, we were able to completely reproduce 46 (65%) of 71 datasets. While our results show that a majority of studies are reproducible, they highlight the fact that many studies still are not the carefully curated research that the scientific community and public expects.

Patterns of domestication in the Ethiopian oil-seed crop noug (Guizotia abyssinica).

Noug (Guizotia abyssinica) is a semidomesticated oil-seed crop, which is primarily cultivated in Ethiopia. Unlike its closest crop relative, sunflower, noug has small seeds, small flowering heads, many branches, many flowering heads, and indeterminate flowering, and it shatters in the field. Here, we conducted common garden studies and microsatellite analyses of genetic variation to test whether high levels of crop-wild gene flow and/or unfavorable phenotypic correlations have hindered noug domestication. With the exception of one population, analyses of microsatellite variation failed to detect substantial recent admixture between noug and its wild progenitor. Likewise, only very weak correlations were found between seed mass and the number or size of flowering heads. Thus, noug's 'atypical' domestication syndrome does not seem to be a consequence of recent introgression or unfavorable phenotypic correlations. Nonetheless, our data do reveal evidence of local adaptation of noug cultivars to different precipitation regimes, as well as high levels of phenotypic plasticity, which may permit reasonable yields under diverse environmental conditions. Why noug has not been fully domesticated remains a mystery, but perhaps early farmers selected for resilience to episodic drought or untended environments rather than larger seeds. Domestication may also have been slowed by noug's outcrossing mating system.

Genomic variation across landscapes: insights and applications.

The distribution of genomic variation across landscapes can provide insights into the complex interactions between the environment and the genome that influence the distribution of species, and mediate phenotypic adaptation to local conditions. High throughput sequencing technologies now offer unprecedented power to explore these interactions, allowing powerful inferences about historical processes of colonization, gene flow and divergence, as well as the identification of loci that mediate local adaptation. These 'landscape genomic' approaches have been validated in model species and are now being applied to nonmodel organisms, including foundation species that have substantial effects on ecosystem processes. Here we review the growing field of landscape genomics from a very broad perspective. In particular, we describe the inferential power that is gained by taking a genome-wide view of genetic variation, strategies for study design to best capture adaptive variation, and how to apply this information to practical challenges, such as restoration.

Genomic analysis of a migratory divide reveals candidate genes for migration and implicates selective sweeps in generating islands of differentiation.

Differential gene flow, reductions in diversity following linked selection and/or features of the genome can structure patterns of genomic differentiation during the process of speciation. Possible sources of reproductive isolation are well studied between coastal and inland subspecies groups of Swainson's thrushes, with differences in seasonal migratory behaviour likely playing a key role in reducing hybrid fitness. We assembled and annotated a draft reference genome for this species and generated whole-genome shotgun sequence data for populations adjacent to the hybrid zone between these groups. We documented substantial genomewide heterogeneity in relative estimates of genetic differentiation between the groups. Within population diversity was lower in areas of high relative differentiation, supporting a role for selective sweeps in generating this pattern. Absolute genetic differentiation was reduced in these areas, further suggesting that recurrent selective sweeps in the ancestral population and/or between divergent populations following secondary contact likely occurred. Relative genetic differentiation was also higher near centromeres and on the Z chromosome, suggesting that features of the genome also contribute to genomewide heterogeneity. Genes linked to migratory traits were concentrated in islands of differentiation, supporting previous suggestions that seasonal migration is under divergent selection between Swainson's thrushes. Differences in migratory behaviour likely play a central role in the speciation of many taxa; we developed the infrastructure here to permit future investigations into the role several candidate genes play in reducing gene flow between not only Swainson's thrushes but other species as well.

On the adaptive value of cytoplasmic genomes in plants.

Is DNA variation maintained in organelle genomes selectively neutral? The answer to this question has important implications for many aspects of ecology and evolution. While traditionally the answer has been 'yes', recent studies in animals have shown that, on the contrary, mitochondrial DNA polymorphism is frequently adaptive. In plants, however, the neutrality assumption has not been strongly challenged. Here, we begin with a critical evaluation of arguments in favour of this long-held view. We then discuss the latest empirical evidence for the opposing prediction that sequence variation in plant cytoplasmic genomes is frequently adaptive. While outstanding research progress is being made towards understanding this fundamental topic, we highlight the need for studies that combine information ranging from field experiments to physiology to molecular evolutionary biology. Such an interdisciplinary approach provides a means for determining the frequency, drivers and evolutionary significance of adaptive organelle DNA variation.

Explaining intraspecific diversity in plant secondary metabolites in an ecological context.

Plant secondary metabolites (PSMs) are ubiquitous in plants and play many ecological roles. Each compound can vary in presence and/or quantity, and the composition of the mixture of chemicals can vary, such that chemodiversity can be partitioned within and among individuals. Plant ontogeny and environmental and genetic variation are recognized as sources of chemical variation, but recent advances in understanding the molecular basis of variation may allow the future deployment of isogenic mutants to test the specific adaptive function of variation in PSMs. An important consequence of high intraspecific variation is the capacity to evolve rapidly. It is becoming increasingly clear that trait variance linked to both macro- and micro-environmental variation can also evolve and may respond more strongly to selection than mean trait values. This research, which is in its infancy in plants, highlights what could be a missing piece of the picture of PSM evolution. PSM polymorphisms are probably maintained by multiple selective forces acting across many spatial and temporal scales, but convincing examples that recognize the diversity of plant population structures are rare. We describe how diversity can be inherently beneficial for plants and suggest fruitful avenues for future research to untangle the causes and consequences of intraspecific variation.

The availability of research data declines rapidly with article age.

Policies ensuring that research data are available on public archives are increasingly being implemented at the government [1], funding agency [2-4], and journal [5, 6] level. These policies are predicated on the idea that authors are poor stewards of their data, particularly over the long term [7], and indeed many studies have found that authors are often unable or unwilling to share their data [8-11]. However, there are no systematic estimates of how the availability of research data changes with time since publication. We therefore requested data sets from a relatively homogenous set of 516 articles published between 2 and 22 years ago, and found that availability of the data was strongly affected by article age. For papers where the authors gave the status of their data, the odds of a data set being extant fell by 17% per year. In addition, the odds that we could find a working e-mail address for the first, last, or corresponding author fell by 7% per year. Our results reinforce the notion that, in the long term, research data cannot be reliably preserved by individual researchers, and further demonstrate the urgent need for policies mandating data sharing via public archives.

Divergence is focused on few genomic regions early in speciation: incipient speciation of sunflower ecotypes.

Early in speciation, as populations undergo the transition from local adaptation to incipient species, is when a number of transient, but potentially important, processes appear to be most easily detected. These include signatures of selective sweeps that can point to asymmetry in selection between habitats, divergence hitchhiking, and associations of adaptive genes with environments. In a genomic comparison of ecotypes of the prairie sunflower, Helianthus petiolaris, occurring at Great Sand Dunes National Park and Preserve (Colorado), we found that selective sweeps were mainly restricted to the dune ecotype and that there was variation across the genome in whether proximity to the nondune population constrained or promoted divergence. The major regions of divergence were few and large between ecotypes, in contrast with an interspecific comparison between H. petiolaris and a sympatric congener, Helianthus annuus. In general, the large regions of divergence observed in the ecotypic comparison swamped locus-specific associations with environmental variables. In both comparisons, regions of high divergence occurred in portions of the genetic map with high marker density, probably reflecting regions of low recombination. The difference in genomic distributions of highly divergent regions between ecotypic and interspecific comparisons highlights the value of studies spanning the spectrum of speciation in related taxa.