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.Non peer reviewe

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

Author summary Here we built up a mathematically justified bridge between 1) parametric approaches and 2) co-expression networks in light of identifying molecular interactions underlying complex traits. We first shared our concern that methodological improvements around these schemes, adjusting only their power and scalability, are bounded by more fundamental scheme-specific limitations. Subsequently, our theoretical results were exploited to overcome these limitations to find gene-by-gene interactions neither of which can capture alone. We also aimed to illustrate how this framework enables the interpretation of co-expression networks in a more parametric sense to achieve systematic insights into complex biological processes more reliably. The main procedure was fit for various types of biological applications and high-dimensional data to cover the area of systems biology as broadly as possible. In particular, we chose to illustrate the method's applicability for gene-profile based risk-stratification in cancer research using public acute myeloid leukemia datasets. 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.Peer reviewe

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 theHopscotch that predates domestication led us to investigate its prevalence and potential role in teosinte. We assessed the prevalence of the Hopscotchelement 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 Hopscotchelement 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

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.Peer reviewe

Utilization of Tissue Ploidy Level Variation in de Novo Transcriptome Assembly of Pinus sylvestris

Compared to angiosperms, gymnosperms lag behind in the availability of assembled and annotated genomes. Most genomic analyses in gymnosperms, especially conifer tree species, rely on the use of de novo assembled transcriptomes. However, the level of allelic redundancy and transcript fragmentation in these assembled transcriptomes, and their effect on downstream applications have not been fully investigated. Here, we assessed three assembly strategies for short-reads data, including the utility of haploid megagametophyte tissue during de novo assembly as single-allele guides, for six individuals and five different tissues in Pinus sylvestris. We then contrasted haploid and diploid tissue genotype calls obtained from the assembled transcriptomes to evaluate the extent of paralog mapping. The use of the haploid tissue during assembly increased its completeness without reducing the number of assembled transcripts. Our results suggest that current strategies that rely on available genomic resources as guidance to minimize allelic redundancy are less effective than the application of strategies that cluster redundant assembled transcripts. The strategy yielding the lowest levels of allelic redundancy among the assembled transcriptomes assessed here was the generation of SuperTranscripts with Lace followed by CD-HIT clustering. However, we still observed some levels of heterozygosity (multiple gene fragments per transcript reflecting allelic redundancy) in this assembled transcriptome on the haploid tissue, indicating that further filtering is required before using these assemblies for downstream applications. We discuss the influence of allelic redundancy when these reference transcriptomes are used to select regions for probe design of exome capture baits and for estimation of population genetic diversity.Peer reviewe

Genomics of clinal local adaptation in Pinus sylvestris under continuous environmental and spatial genetic setting

Understanding the consequences of local adaptation at the genomic diversity is a central goal in evolutionary genetics of natural populations. In species with large continuous geographical distributions the phenotypic signal of local adaptation is frequently clear, but the genetic basis often remains elusive. We examined the patterns of genetic diversity inPinus sylvestris, a keystone species in many Eurasian ecosystems with a huge distribution range and decades of forestry research showing that it is locally adapted to the vast range of environmental conditions. MakingP. sylvestrisan even more attractive subject of local adaptation study, population structure has been shown to be weak previously and in this study. However, little is known about the molecular genetic basis of adaptation, as the massive size of gymnosperm genomes has prevented large scale genomic surveys. We generated a both geographically and genomically extensive dataset using a targeted sequencing approach. By applying divergence-based and landscape genomics methods we identified several loci contributing to local adaptation, but only few with large allele frequency changes across latitude. We also discovered a very large (ca. 300 Mbp) putative inversion potentially under selection, which to our knowledge is the first such discovery in conifers. Our results call for more detailed analysis of structural variation in relation to genomic basis of local adaptation, emphasize the lack of large effect loci contributing to local adaptation in the coding regions and thus point out the need for more attention toward multi-locus analysis of polygenic adaptation

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

Peer reviewe

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.Peer reviewe

Complex patterns of local adaptation in teosinte

Populations of widely distributed species often encounter and adapt to specific environmental conditions. However, comprehensive characterization of the genetic basis of adaptation is demanding, requiring genome-wide genotype data, multiple sampled populations, and a good understanding of population structure. We have used environmental and high-density genotype data to describe the genetic basis of local adaptation in 21 populations of teosinte, the wild ancestor of maize. We found that altitude, dispersal events and admixture among subspecies formed a complex hierarchical genetic structure within teosinte. Patterns of linkage disequilibrium revealed four mega-base scale inversions that segregated among populations and had altitudinal clines. Based on patterns of differentiation and correlation with environmental variation, inversions and nongenic regions play an important role in local adaptation of teosinte. Further, we note that strongly differentiated individual populations can bias the identification of adaptive loci. The role of inversions in local adaptation has been predicted by theory and requires attention as genome-wide data become available for additional plant species. These results also suggest a potentially important role for noncoding variation, especially in large plant genomes in which the gene space represents a fraction of the entire genome

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 were particularly 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 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 informing 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