Genomic determinants, architecture, and constraints in drought-related traits in Corymbia calophylla.

BACKGROUND: Drought adaptation is critical to many tree species persisting under climate change, however our knowledge of the genetic basis for trees to adapt to drought is limited. This knowledge gap impedes our fundamental understanding of drought response and application to forest production and conservation. To improve our understanding of the genomic determinants, architecture, and trait constraints, we assembled a reference genome and detected ~ 6.5 M variants in 432 phenotyped individuals for the foundational tree Corymbia calophylla. RESULTS: We found 273 genomic variants determining traits with moderate heritability (h2SNP = 0.26-0.64). Significant variants were predominantly in gene regulatory elements distributed among several haplotype blocks across all chromosomes. Furthermore, traits were constrained by frequent epistatic and pleiotropic interactions. CONCLUSIONS: Our results on the genetic basis for drought traits in Corymbia calophylla have several implications for the ability to adapt to climate change: (1) drought related traits are controlled by complex genomic architectures with large haplotypes, epistatic, and pleiotropic interactions; (2) the most significant variants determining drought related traits occurred in regulatory regions; and (3) models incorporating epistatic interactions increase trait predictions. Our findings indicate that despite moderate heritability drought traits are likely constrained by complex genomic architecture potentially limiting trees response to climate change.

Interspecies genome divergence is predominantly due to frequent small scale rearrangements in Eucalyptus.

Synteny, the ordering of sequences within homologous chromosomes, must be maintained within the genomes of sexually reproducing species for the sharing of alleles and production of viable, reproducing offspring. However, when the genomes of closely related species are compared, a loss of synteny is often observed. Unequal homologous recombination is the primary mechanism behind synteny loss, occurring more often in transposon rich regions, and resulting in the formation of chromosomal rearrangements. To examine patterns of synteny among three closely related, interbreeding, and wild Eucalyptus species, we assembled their genomes using long-read DNA sequencing and de novo assembly. We identify syntenic and rearranged regions between these genomes and estimate that ~48% of our genomes remain syntenic while ~36% is rearranged. We observed that rearrangements highly fragment microsynteny. Our results suggest that synteny between these species is primarily lost through small-scale rearrangements, not through sequence loss, gain, or sequence divergence. Further examination of identified rearrangements suggests that rearrangements may be altering the phenotypes of Eucalyptus species. Our study also underscores that the use of single reference genomes in genomic variation studies could lead to reference bias, especially given the scale at which we show potentially adaptive loci have highly diverged, deleted, duplicated and/or rearranged. This study provides an unbiased framework to look at potential speciation and adaptive loci among a rapidly radiating foundation species of woodland trees that are free from selective breeding seen in most crop species.

Genomic consequences of artificial selection during early domestication of a wood fibre crop.

From its origins in Australia, Eucalyptus grandis has spread to every continent, except Antarctica, as a wood crop. It has been cultivated and bred for over 100 yr in places such as South Africa. Unlike most annual crops and fruit trees, domestication of E. grandis is still in its infancy, representing a unique opportunity to interrogate the genomic consequences of artificial selection early in the domestication process. To determine how a century of artificial selection has changed the genome of E. grandis, we generated single nucleotide polymorphism genotypes for 1080 individuals from three advanced South African breeding programmes using the EUChip60K chip, and investigated population structure and genome-wide differentiation patterns relative to wild progenitors. Breeding and wild populations appeared genetically distinct. We found genomic evidence of evolutionary processes known to have occurred in other plant domesticates, including interspecific introgression and intraspecific infusion from wild material. Furthermore, we found genomic regions with increased linkage disequilibrium and genetic differentiation, putatively representing early soft sweeps of selection. This is, to our knowledge, the first study of genomic signatures of domestication in a timber species looking beyond the first few generations of cultivation. Our findings highlight the importance of intra- and interspecific hybridization during early domestication.

Genomic evidence of introgression and adaptation in a model subtropical tree species, Eucalyptus grandis.

The genetic consequences of adaptation to changing environments can be deciphered using population genomics, which may help predict species' responses to global climate change. Towards this, we used genome-wide SNP marker analysis to determine population structure and patterns of genetic differentiation in terms of neutral and adaptive genetic variation in the natural range of Eucalyptus grandis, a widely cultivated subtropical and temperate species, serving as genomic reference for the genus. We analysed introgression patterns at subchromosomal resolution using a modified ancestry mapping approach and identified provenances with extensive interspecific introgression in response to increased aridity. Furthermore, we describe potentially adaptive genetic variation as explained by environment-associated SNP markers, which also led to the discovery of what is likely a large structural variant. Finally, we show that genes linked to these markers are enriched for biotic and abiotic stress responses.

Spatial, climate and ploidy factors drive genomic diversity and resilience in the widespread grass Themeda triandra.

Global climate change poses a significant threat to natural communities around the world, with many plant species showing signs of climate stress. Grassland ecosystems are not an exception, with climate change compounding contemporary pressures such as habitat loss and fragmentation. In this study, we assess the climate resilience of Themeda triandra, a foundational species and the most widespread plant in Australia, by assessing the relative contributions of spatial, environmental and ploidy factors to contemporary genomic variation. Reduced-representation genome sequencing on 472 samples from 52 locations was used to test how the distribution of genomic variation, including ploidy polymorphism, supports adaptation to hotter and drier climates. We explicitly quantified isolation by distance (IBD) and isolation by environment (IBE) and predicted genomic vulnerability of populations to future climates based on expected deviation from current genomic composition. We found that a majority (54%) of genomic variation could be attributed to IBD, while an additional 22% (27% when including ploidy information) could be explained by two temperature and two precipitation climate variables demonstrating IBE. Ploidy polymorphisms were common within populations (31/52 populations), indicating that ploidy mixing is characteristic of T. triandra populations. Genomic vulnerabilities were found to be heterogeneously distributed throughout the landscape, and our analysis suggested that ploidy polymorphism, along with other factors linked to polyploidy, reduced vulnerability to future climates by 60% (0.25-0.10). Our data suggests that polyploidy may facilitate adaptation to hotter climates and highlight the importance of incorporating ploidy in adaptive management strategies to promote the resilience of this and other foundation species.

Rapid Recovery Gene Downregulation during Excess-Light Stress and Recovery in Arabidopsis.

Stress recovery may prove to be a promising approach to increase plant performance and, theoretically, mRNA instability may facilitate faster recovery. Transcriptome (RNA-seq, qPCR, sRNA-seq, and PARE) and methylome profiling during repeated excess-light stress and recovery was performed at intervals as short as 3 min. We demonstrate that 87% of the stress-upregulated mRNAs analyzed exhibit very rapid recovery. For instance, HSP101 abundance declined 2-fold every 5.1 min. We term this phenomenon rapid recovery gene downregulation (RRGD), whereby mRNA abundance rapidly decreases promoting transcriptome resetting. Decay constants (k) were modeled using two strategies, linear and nonlinear least squares regressions, with the latter accounting for both transcription and degradation. This revealed extremely short half-lives ranging from 2.7 to 60.0 min for 222 genes. Ribosome footprinting using degradome data demonstrated RRGD loci undergo cotranslational decay and identified changes in the ribosome stalling index during stress and recovery. However, small RNAs and 5'-3' RNA decay were not essential for recovery of the transcripts examined, nor were any of the six excess light-associated methylome changes. We observed recovery-specific gene expression networks upon return to favorable conditions and six transcriptional memory types. In summary, rapid transcriptome resetting is reported in the context of active recovery and cellular memory.

HOME: a histogram based machine learning approach for effective identification of differentially methylated regions.

BACKGROUND: The development of whole genome bisulfite sequencing has made it possible to identify methylation differences at single base resolution throughout an entire genome. However, a persistent challenge in DNA methylome analysis is the accurate identification of differentially methylated regions (DMRs) between samples. Sensitive and specific identification of DMRs among different conditions requires accurate and efficient algorithms, and while various tools have been developed to tackle this problem, they frequently suffer from inaccurate DMR boundary identification and high false positive rate. RESULTS: We present a novel Histogram Of MEthylation (HOME) based method that takes into account the inherent difference in the distribution of methylation levels between DMRs and non-DMRs to discriminate between the two using a Support Vector Machine. We show that generated features used by HOME are dataset-independent such that a classifier trained on, for example, a mouse methylome training set of regions of differentially accessible chromatin, can be applied to any other organism's dataset and identify accurate DMRs. We demonstrate that DMRs identified by HOME exhibit higher association with biologically relevant genes, processes, and regulatory events compared to the existing methods. Moreover, HOME provides additional functionalities lacking in most of the current DMR finders such as DMR identification in non-CG context and time series analysis. HOME is freely available at https://github.com/ListerLab/HOME . CONCLUSION: HOME produces more accurate DMRs than the current state-of-the-art methods on both simulated and biological datasets. The broad applicability of HOME to identify accurate DMRs in genomic data from any organism will have a significant impact upon expanding our knowledge of how DNA methylation dynamics affect cell development and differentiation.

DNA methylation profiles of diverse Brachypodium distachyon align with underlying genetic diversity.

DNA methylation, a common modification of genomic DNA, is known to influence the expression of transposable elements as well as some genes. Although commonly viewed as an epigenetic mark, evidence has shown that underlying genetic variation, such as transposable element polymorphisms, often associate with differential DNA methylation states. To investigate the role of DNA methylation variation, transposable element polymorphism, and genomic diversity, whole-genome bisulfite sequencing was performed on genetically diverse lines of the model cereal Brachypodium distachyon Although DNA methylation profiles are broadly similar, thousands of differentially methylated regions are observed between lines. An analysis of novel transposable element indel variation highlighted hundreds of new polymorphisms not seen in the reference sequence. DNA methylation and transposable element variation is correlated with the genome-wide amount of genetic variation present between samples. However, there was minimal evidence that novel transposon insertions or deletions are associated with nearby differential methylation. This study highlights unique relationships between genetic variation and DNA methylation variation within Brachypodium and provides a valuable map of DNA methylation across diverse resequenced accessions of this model cereal species.

Axe: rapid, competitive sequence read demultiplexing using a trie.

Summary: We describe a rapid algorithm for demultiplexing DNA sequence reads with in-read indices. Axe selects the optimal index present in a sequence read, even in the presence of sequencing errors. The algorithm is able to handle combinatorial indexing, indices of differing length and several mismatches per index sequence. Availability and implementation: Axe is implemented in C, and is used as a command-line program on Unix-like systems. Axe is available online at https://github.com/kdmurray91/axe, and is available in Debian/Ubuntu distributions of GNU/Linux as the package axe-demultiplexer. Supplementary information: Supplementary data are available at Bioinformatics online.

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.

Adaptation genomics: The angel is in the details.

This article comments on: Genome-wide signatures of flowering adaptation to climate temperature: Regional analyses in a highly diverse native range of Arabidopsis thaliana.

Genomic diversity guides conservation strategies among rare terrestrial orchid species when taxonomy remains uncertain.

Background and Aims: Species are often used as the unit for conservation, but may not be suitable for species complexes where taxa are difficult to distinguish. Under such circumstances, it may be more appropriate to consider species groups or populations as evolutionarily significant units (ESUs). A population genomic approach was employed to investigate the diversity within and among closely related species to create a more robust, lineage-specific conservation strategy for a nationally endangered terrestrial orchid and its relatives from south-eastern Australia. Methods: Four putative species were sampled from a total of 16 populations in the Victorian Volcanic Plain (VVP) bioregion and one population of a sub-alpine outgroup in south-eastern Australia. Morphological measurements were taken in situ along with leaf material for genotyping by sequencing (GBS) and microsatellite analyses. Key Results: Species could not be differentiated using morphological measurements. Microsatellite and GBS markers confirmed the outgroup as distinct, but only GBS markers provided resolution of population genetic structure. The nationally endangered Diuris basaltica was indistinguishable from two related species ( D. chryseopsis and D. behrii ), while the state-protected D. gregaria showed genomic differentiation. Conclusions: Genomic diversity identified among the four Diuris species suggests that conservation of this taxonomically complex group will be best served by considering them as one ESU rather than separately aligned with species as currently recognized. This approach will maximize evolutionary potential among all species during increased isolation and environmental change. The methods used here can be applied generally to conserve evolutionary processes for groups where taxonomic uncertainty hinders the use of species as conservation units.

Population and phylogenomic decomposition via genotyping-by-sequencing in Australian Pelargonium.

Species delimitation has seen a paradigm shift as increasing accessibility of genomic-scale data enables separation of lineages with convergent morphological traits and the merging of recently diverged ecotypes that have distinguishing characteristics. We inferred the process of lineage formation among Australian species in the widespread and highly variable genus Pelargonium by combining phylogenomic and population genomic analyses along with breeding system studies and character analysis. Phylogenomic analysis and population genetic clustering supported seven of the eight currently described species but provided little evidence for differences in genetic structure within the most widely distributed group that containing P. australe. In contrast, morphometric analysis detected three deep lineages within Australian Pelargonium; with P. australe consisting of five previously unrecognized entities occupying separate geographic ranges. The genomic approach enabled elucidation of parallel evolution in some traits formerly used to delineate species, as well as identification of ecotypic morphological differentiation within recognized species. Highly variable morphology and trait convergence each contribute to the discordance between phylogenomic relationships and morphological taxonomy. Data suggest that genetic divergence among species within the Australian Pelargonium may result from allopatric speciation while morphological differentiation within and among species may be more strongly driven by environmental differences.

Natural allelic variation underlying a major fitness trade-off in Arabidopsis thaliana.

Plants can defend themselves against a wide array of enemies, from microbes to large animals, yet there is great variability in the effectiveness of such defences, both within and between species. Some of this variation can be explained by conflicting pressures from pathogens with different modes of attack. A second explanation comes from an evolutionary 'tug of war', in which pathogens adapt to evade detection, until the plant has evolved new recognition capabilities for pathogen invasion. If selection is, however, sufficiently strong, susceptible hosts should remain rare. That this is not the case is best explained by costs incurred from constitutive defences in a pest-free environment. Using a combination of forward genetics and genome-wide association analyses, we demonstrate that allelic diversity at a single locus, ACCELERATED CELL DEATH 6 (ACD6), underpins marked pleiotropic differences in both vegetative growth and resistance to microbial infection and herbivory among natural Arabidopsis thaliana strains. A hyperactive ACD6 allele, compared to the reference allele, strongly enhances resistance to a broad range of pathogens from different phyla, but at the same time slows the production of new leaves and greatly reduces the biomass of mature leaves. This allele segregates at intermediate frequency both throughout the worldwide range of A. thaliana and within local populations, consistent with this allele providing substantial fitness benefits despite its marked impact on growth.

A methyl transferase links the circadian clock to the regulation of alternative splicing.

Circadian rhythms allow organisms to time biological processes to the most appropriate phases of the day-night cycle. Post-transcriptional regulation is emerging as an important component of circadian networks, but the molecular mechanisms linking the circadian clock to the control of RNA processing are largely unknown. Here we show that PROTEIN ARGININE METHYL TRANSFERASE 5 (PRMT5), which transfers methyl groups to arginine residues present in histones and Sm spliceosomal proteins, links the circadian clock to the control of alternative splicing in plants. Mutations in PRMT5 impair several circadian rhythms in Arabidopsis thaliana and this phenotype is caused, at least in part, by a strong alteration in alternative splicing of the core-clock gene PSEUDO RESPONSE REGULATOR 9 (PRR9). Furthermore, genome-wide studies show that PRMT5 contributes to the regulation of many pre-messenger-RNA splicing events, probably by modulating 5'-splice-site recognition. PRMT5 expression shows daily and circadian oscillations, and this contributes to the mediation of the circadian regulation of expression and alternative splicing of a subset of genes. Circadian rhythms in locomotor activity are also disrupted in dart5-1, a mutant affected in the Drosophila melanogaster PRMT5 homologue, and this is associated with alterations in splicing of the core-clock gene period and several clock-associated genes. Our results demonstrate a key role for PRMT5 in the regulation of alternative splicing and indicate that the interplay between the circadian clock and the regulation of alternative splicing by PRMT5 constitutes a common mechanism that helps organisms to synchronize physiological processes with daily changes in environmental conditions.

Genome-wide association study of 107 phenotypes in Arabidopsis thaliana inbred lines.

Although pioneered by human geneticists as a potential solution to the challenging problem of finding the genetic basis of common human diseases, genome-wide association (GWA) studies have, owing to advances in genotyping and sequencing technology, become an obvious general approach for studying the genetics of natural variation and traits of agricultural importance. They are particularly useful when inbred lines are available, because once these lines have been genotyped they can be phenotyped multiple times, making it possible (as well as extremely cost effective) to study many different traits in many different environments, while replicating the phenotypic measurements to reduce environmental noise. Here we demonstrate the power of this approach by carrying out a GWA study of 107 phenotypes in Arabidopsis thaliana, a widely distributed, predominantly self-fertilizing model plant known to harbour considerable genetic variation for many adaptively important traits. Our results are dramatically different from those of human GWA studies, in that we identify many common alleles of major effect, but they are also, in many cases, harder to interpret because confounding by complex genetics and population structure make it difficult to distinguish true associations from false. However, a-priori candidates are significantly over-represented among these associations as well, making many of them excellent candidates for follow-up experiments. Our study demonstrates the feasibility of GWA studies in A. thaliana and suggests that the approach will be appropriate for many other organisms.

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.

New Arabidopsis advanced intercross recombinant inbred lines reveal female control of nonrandom mating.

Female control of nonrandom mating has never been genetically established, despite being linked to inbreeding depression and sexual selection. In order to map the loci that control female-mediated nonrandom mating, we constructed a new advanced intercross recombinant inbred line (RIL) population derived from a cross between Arabidopsis (Arabidopsis thaliana) accessions Vancouver (Van-0) and Columbia (Col-0) and mapped quantitative trait loci (QTLs) responsible for nonrandom mating and seed yield traits. We genotyped a population of 490 RILs. A subset of these lines was used to construct an expanded map of 1,061.4 centimorgans with an average interval of 6.7±5.3 centimorgans between markers. QTLs were then mapped for female- and male-mediated nonrandom mating and seed yield traits. To map the genetic loci responsible for female-mediated nonrandom mating and seed yield, we performed mixed pollinations with genetically marked Col-0 pollen and Van-0 pollen on RIL pistils. To map the loci responsible for male-mediated nonrandom mating and seed yield, we performed mixed pollinations with genetically marked Col-0 and RIL pollen on Van-0 pistils. Composite interval mapping of these data identified four QTLs that control female-mediated nonrandom mating and five QTLs that control female-mediated seed yield. We also identified four QTLs that control male-mediated nonrandom mating and three QTLs that control male-mediated seed yield. Epistasis analysis indicates that several of these loci interact. To our knowledge, the results of these experiments represent the first time female-mediated nonrandom mating has been genetically defined.

Genetic architecture of regulatory variation in Arabidopsis thaliana.

Studying the genetic regulation of expression variation is a key method to dissect complex phenotypic traits. To examine the genetic architecture of regulatory variation in Arabidopsis thaliana, we performed genome-wide association (GWA) mapping of gene expression in an F(1) hybrid diversity panel. At a genome-wide false discovery rate (FDR) of 0.2, an associated single nucleotide polymorphism (SNP) explains >38% of trait variation. In comparison with SNPs that are distant from the genes to which they were associated, locally associated SNPs are preferentially found in regions with extended linkage disequilibrium (LD) and have distinct population frequencies of the derived alleles (where Arabidopsis lyrata has the ancestral allele), suggesting that different selective forces are acting. Locally associated SNPs tend to have additive inheritance, whereas distantly associated SNPs are primarily dominant. In contrast to results from mapping of expression quantitative trait loci (eQTL) in linkage studies, we observe extensive allelic heterogeneity for local regulatory loci in our diversity panel. By association mapping of allele-specific expression (ASE), we detect a significant enrichment for cis-acting variation in local regulatory variation. In addition to gene expression variation, association mapping of splicing variation reveals both local and distant genetic regulation for intron and exon level traits. Finally, we identify candidate genes for 59 diverse phenotypic traits that were mapped to eQTL.

Population genomic variation reveals roles of history, adaptation and ploidy in switchgrass.

Geographic patterns of genetic variation are shaped by multiple evolutionary processes, including genetic drift, migration and natural selection. Switchgrass (Panicum virgatum L.) has strong genetic and adaptive differentiation despite life history characteristics that promote high levels of gene flow and can homogenize intraspecific differences, such as wind-pollination and self-incompatibility. To better understand how historical and contemporary factors shape variation in switchgrass, we use genotyping-by-sequencing to characterize switchgrass from across its range at 98 042 SNPs. Population structuring reflects biogeographic and ploidy differences within and between switchgrass ecotypes and indicates that biogeographic history, ploidy incompatibilities and differential adaptation each have important roles in shaping ecotypic differentiation in switchgrass. At one extreme, we determine that two Panicum taxa are not separate species but are actually conspecific, ecologically divergent types of switchgrass adapted to the extreme conditions of coastal sand dune habitats. Conversely, we identify natural hybrids among lowland and upland ecotypes and visualize their genome-wide patterns of admixture. Furthermore, we determine that genetic differentiation between primarily tetraploid and octoploid lineages is not caused solely by ploidy differences. Rather, genetic diversity in primarily octoploid lineages is consistent with a history of admixture. This suggests that polyploidy in switchgrass is promoted by admixture of diverged lineages, which may be important for maintaining genetic differentiation between switchgrass ecotypes where they are sympatric. These results provide new insights into the mechanisms shaping variation in widespread species and provide a foundation for dissecting the genetic basis of adaptation in switchgrass.