Clinal variation in the non-acclimated and cold-acclimated freezing tolerance of Arabidopsis thaliana accessions

Arabidopsis thaliana is a geographically widely spread species consisting of local accessions differing both genetically and phenotypically. These differences may constitute environmental adaptations and a latitudinal cline in freezing tolerance has been shown previously. Many plants, including Arabidopsis, exhibit increased freezing tolerance after cold exposure (cold acclimation). Here we present evidence for geographical clines (both latitudinal and longitudinal) in acclimated (ACC) and non-acclimated (NA) freezing tolerance, estimated from electrolyte leakage measurements on 54 accessions. Leaf Pro contents were not correlated with freezing tolerance, while sugar contents (Glc, Fru, Suc, Raf) were in the ACC, but not the NA state. Expression levels of 14 cold-induced genes were investigated before and after 2 weeks of cold acclimation by quantitative RT-PCR. Expression of the CBF1, 2 and 3 genes was not correlated with freezing tolerance. The expression of some CBF-regulated (COR) genes, however, was correlated specifically with ACC freezing tolerance. A tight correlation between CBF and COR gene expression was only observed under non-acclimating conditions, where CBF and COR expression were also correlated with the expression of PRR5, a component of the circadian clock. Collectively, this study sheds new light on the molecular determinants of plant-freezing tolerance and cold acclimation and their geographical dependence.     

Understanding genetic control of root system architecture in soybean: Insights into the genetic basis of lateral root number

Developing crops with better root systems is a promising strategy to ensure productivity in both optimum and stress environments. Root system architectural traits in 397 soybean accessions were characterized and a high‐density single nucleotide polymorphisms (SNPs)‐based genome‐wide association study was performed to identify the underlying genes associated with root structure. SNPs associated with root architectural traits specific to landraces and elite germplasm pools were detected. Four loci were detected in landraces for lateral root number (LRN) and distribution of root thickness in diameter Class I with a major locus on chromosome 16. This major loci was detected in the coding region of unknown protein, and subsequent analyses demonstrated that root traits are affected with mutated haplotypes of the gene. In elite germplasm pool, 3 significant SNPs in alanine‐glyoxalate aminotransferase, Leucine‐Rich Repeat receptor/No apical meristem, and unknown functional genes were found to govern multiple traits including root surface area and volume. However, no major loci were detected for LRN in elite germplasm. Nucleotide diversity analysis found evidence of selective sweeps around the landraces LRN gene. Soybean accessions with minor and mutated allelic variants of LRN gene were found to perform better in both water‐limited and optimal field conditions.     

Nested core collections maximizing genetic diversity in Arabidopsis thaliana

The successful exploitation of natural genetic diversity requires a basic knowledge of the extent of the variation present in a species. To study natural variation in Arabidopsis thaliana, we defined nested core collections maximizing the diversity present among a worldwide set of 265 accessions. The core collections were generated based on DNA sequence data from a limited number of fragments evenly distributed in the genome and were shown to successfully capture the molecular diversity in other loci as well as the morphological diversity. The core collections are available to the scientific community and thus provide an important resource for the study of genetic variation and its functional consequences in Arabidopsis. Moreover, this strategy can be used in other species to provide a rational framework for undertaking diversity surveys, including single nucleotide polymorphism (SNP) discovery and phenotyping, allowing the utilization of genetic variation for the study of complex traits.     

Towards the genetic architecture of seed lipid biosynthesis and accumulation in Arabidopsis thaliana

We report the quantitative genetic analysis of seed oil quality and quantity in six Arabidopsis thaliana recombinant inbred populations, in which the parent accessions were from diverse geographical origins, and were selected on the basis of variation for seed oil content and lipid composition. Although most of the biochemical steps involved in lipid biosynthesis are known and the key genes have been identified, the regulation of the processes that results in the final oil composition and total amount is not understood. By using physically anchored markers it was possible to compare results across populations. A total of 219 quantitative trait loci (QTLs) were identified, of which 81 were significant at P<0.001. Some of these colocalise with QTLs identified previously, but many novel QTLs were also identified. The results highlight the importance of studying traits in multiple populations, which will lead to a better understanding of the contribution that natural variation makes to the genetic architecture of a phenotype.     

Optimising the identification of causal variants across varying genetic architectures in crops

Association studies use statistical links between genetic markers and the phenotype variation across many individuals to identify genes controlling variation in the target phenotype. However, this approach, particularly conducted on a genome‐wide scale (GWAS), has limited power to identify the genes responsible for variation in traits controlled by complex genetic architectures. In this study, we employ real‐world genotype datasets from four crop species with distinct minor allele frequency distributions, population structures and linkage disequilibrium patterns. We demonstrate that different GWAS statistical approaches provide favourable trade‐offs between power and accuracy for traits controlled by different types of genetic architectures. FarmCPU provides the most favourable outcomes for moderately complex traits while a Bayesian approach adopted from genomic prediction provides the most favourable outcomes for extremely complex traits. We assert that by estimating the complexity of genetic architectures for target traits and selecting an appropriate statistical approach for the degree of complexity detected, researchers can substantially improve the ability to dissect the genetic factors controlling complex traits such as flowering time, plant height and yield component.     

Metabolome‐genome‐wide association study dissects genetic architecture for generating natural variation in rice secondary metabolism

Plants produce structurally diverse secondary (specialized) metabolites to increase their fitness for survival under adverse environments. Several bioactive compounds for new drugs have been identified through screening of plant extracts. In this study, genome‐wide association studies (GWAS) were conducted to investigate the genetic architecture behind the natural variation of rice secondary metabolites. GWAS using the metabolome data of 175 rice accessions successfully identified 323 associations among 143 single nucleotide polymorphisms (SNPs) and 89 metabolites. The data analysis highlighted that levels of many metabolites are tightly associated with a small number of strong quantitative trait loci (QTLs). The tight association may be a mechanism generating strains with distinct metabolic composition through the crossing of two different strains. The results indicate that one plant species produces more diverse phytochemicals than previously expected, and plants still contain many useful compounds for human applications.     

Genome‐wide signatures of flowering adaptation to climate temperature: Regional analyses in a highly diverse native range of Arabidopsis thaliana

Current global change is fueling an interest to understand the genetic and molecular mechanisms of plant adaptation to climate. In particular, altered flowering time is a common strategy for escape from unfavourable climate temperature. In order to determine the genomic bases underlying flowering time adaptation to this climatic factor, we have systematically analysed a collection of 174 highly diverse Arabidopsis thaliana accessions from the Iberian Peninsula. Analyses of 1.88 million single nucleotide polymorphisms provide evidence for a spatially heterogeneous contribution of demographic and adaptive processes to geographic patterns of genetic variation. Mountains appear to be allele dispersal barriers, whereas the relationship between flowering time and temperature depended on the precise temperature range. Environmental genome‐wide associations supported an overall genome adaptation to temperature, with 9.4% of the genes showing significant associations. Furthermore, phenotypic genome‐wide associations provided a catalogue of candidate genes underlying flowering time variation. Finally, comparison of environmental and phenotypic genome‐wide associations identified known (Twin Sister of FT, FRIGIDA‐like 1, and Casein Kinase II Beta chain 1) and new (Epithiospecifer Modifier 1 and Voltage‐Dependent Anion Channel 5) genes as candidates for adaptation to climate temperature by altered flowering time. Thus, this regional collection provides an excellent resource to address the spatial complexity of climate adaptation in annual plants.     

Clinal variation in freezing tolerance among natural accessions of Arabidopsis thaliana

Low temperature represents a form of abiotic stress that varies predictably with latitude and altitude and to which organisms have evolved multiple physiological responses. Plants provide an especially useful experimental system for investigating the ecological and evolutionary dynamics of tolerance to low temperature because of their sessile lifestyle and inability to escape ambient atmospheric conditions. Here, intraspecific variation in freezing tolerance was investigated in Arabidopsis thaliana by conducting freezing tolerance assays on 71 accessions collected from across the native range of the species. Assays were performed at multiple minimum temperatures and on both cold-acclimated and non-cold-acclimated individuals. Considerable variation in freezing tolerance was observed among accessions both with and without a prior cold-acclimation treatment, suggesting that differences among accessions in cold-acclimation capacity as well as differences in intrinsic physiology contribute to variation in this phenotype. A highly significant positive relationship was observed between freezing tolerance and latitude of origin of accessions, consistent with a major role for natural selection in shaping variation in this phenotype. Clinal variation in freezing tolerance in A. thaliana coupled with considerable knowledge of the underlying genetics and physiology of this phenotype should allow evolutionary genetic analysis at multiple levels.     

All sheeps and sizes: a genetic investigation of mature body size across sheep breeds reveals a polygenic nature

Mature body size is genetically correlated with growth rate, an important economic trait in the sheep industry. Mature body size has been studied extensively in humans as well as cattle and other domestic animal populations but not in sheep. Six‐hundred and sixteen ewes, across 22 breeds, were measured for 28 linear measurements representing various skeletal parts. PCA from these measures generated principal components 1 and 2 which represented 66 and 7% of the phenotypic variation respectively. Two‐hundred and twenty sheep were genotyped on the Illumina Ovine HD beadchip for a GWAS investigating mature body size and linear body measurements. Forty‐six (Bonferroni P < 0.05) SNP associations across 14 chromosomes were identified utilizing principal component 1, representing overall body size, revealing mature body size to have fewer loci of large effect than other domestic species such as dogs and horses. Genome‐wide associations for individual linear measures identified major quantitative trait loci for withers height and ear length. Withers height was associated (Bonferroni P < 0.05) with 12 SNPs across six chromosomes whereas ear length was associated with a single locus on chromosome 3, containing MSRB3. This analysis identified several loci known to be associated with mature body size in other species such as NCAPG, LCORL, and HMGA2. Mature body size is more polygenic in sheep than other domesticated species, making the development of genomic selection for the trait the most efficient option for maintaining or reducing mature body size in sheep.     

Heterogeneity of genetic architecture of body size traits in a free‐living population

Knowledge of the underlying genetic architecture of quantitative traits could aid in understanding how they evolve. In wild populations, it is still largely unknown whether complex traits are polygenic or influenced by few loci with major effect, due to often small sample sizes and low resolution of marker panels. Here, we examine the genetic architecture of five adult body size traits in a free‐living population of Soay sheep on St Kilda using 37 037 polymorphic SNPs. Two traits (jaw and weight) show classical signs of a polygenic trait: the proportion of variance explained by a chromosome was proportional to its length, multiple chromosomes and genomic regions explained significant amounts of phenotypic variance, but no SNPs were associated with trait variance when using GWAS. In comparison, genetic variance for leg length traits (foreleg, hindleg and metacarpal) was disproportionately explained by two SNPs on chromosomes 16 (s23172.1) and 19 (s74894.1), which each explained >10% of the additive genetic variance. After controlling for environmental differences, females heterozygous for s74894.1 produced more lambs and recruits during their lifetime than females homozygous for the common allele conferring long legs. We also demonstrate that alleles conferring shorter legs have likely entered the population through a historic admixture event with the Dunface sheep. In summary, we show that different proxies for body size can have very different genetic architecture and that dense SNP helps in understanding both the mode of selection and the evolutionary history at loci underlying quantitative traits in natural populations.     

Evidence for parallel adaptation to climate across the natural range of Arabidopsis thaliana

How organisms adapt to different climate habitats is a key question in evolutionary ecology and biological conservation. Species distributions are often determined by climate suitability. Consequently, the anthropogenic impact on earth's climate is of key concern to conservation efforts because of our relatively poor understanding of the ability of populations to track and evolve to climate change. Here, we investigate the ability of Arabidopsis thaliana to occupy climate space by quantifying the extent to which different climate regimes are accessible to different A. thaliana genotypes using publicly available data from a large‐scale genotyping project and from a worldwide climate database. The genetic distance calculated from 149 single‐nucleotide polymorphisms (SNPs) among 60 lineages of A. thaliana was compared to the corresponding climate distance among collection localities calculated from nine different climatic factors. A. thaliana was found to be highly labile when adapting to novel climate space, suggesting that populations may experience few constraints when adapting to changing climates. Our results also provide evidence of a parallel or convergent evolution on the molecular level supporting recent generalizations regarding the genetics of adaptation.     

Testing for adaptation to climate in Arabidopsis thaliana: a calibrated common garden approach

BACKGROUND AND AIMS: A recent method used to test for local adaptation is a common garden experiment where analyses are calibrated to the environmental conditions of the garden. In this study the calibrated common garden approach is used to test for patterns of adaptation to climate in accessions of Arabidopsis thaliana. METHODS: Seedlings from 21 accessions of A. thaliana were planted outdoors in College Park, MD, USA, and development was monitored during the course of a growing season. ANOVA and multiple regression analysis were used to determine if development traits were significant predictors of plant success. Previously published data relating to accessional differences in genetic and physiological characters were also examined. Historical records of climate were used to evaluate whether properties of the site of origin of an accession affected the fitness of plants in a novel environment. KEY RESULTS: By calibrating the analysis to the climatic conditions of the common garden site, performance differences were detected among the accessions consistent with a pattern of adaptation to latitude and climatic conditions. Relatively higher accession fitness was predicted by a latitude and climatic history similar to that of College Park in April and May during the main growth period of this experiment. The climatic histories of the accessions were better predictors of performance than many of the life-history and growth measures taken during the experiment. CONCLUSIONS: It is concluded that the calibrated common garden experiment can detect local adaptation and guide subsequent reciprocal transplant experiments.     

Paths to selection on life history loci in different natural environments across the native range of Arabidopsis thaliana

Selection on quantitative trait loci (QTL) may vary among natural environments due to differences in the genetic architecture of traits, environment‐specific allelic effects or changes in the direction and magnitude of selection on specific traits. To dissect the environmental differences in selection on life history QTL across climatic regions, we grew a panel of interconnected recombinant inbred lines (RILs) of Arabidopsis thaliana in four field sites across its native European range. For each environment, we mapped QTL for growth, reproductive timing and development. Several QTL were pleiotropic across environments, three colocalizing with known functional polymorphisms in flowering time genes (CRY2, FRI and MAF2‐5), but major QTL differed across field sites, showing conditional neutrality. We used structural equation models to trace selection paths from QTL to lifetime fitness in each environment. Only three QTL directly affected fruit number, measuring fitness. Most QTL had an indirect effect on fitness through their effect on bolting time or leaf length. Influence of life history traits on fitness differed dramatically across sites, resulting in different patterns of selection on reproductive timing and underlying QTL. In two oceanic field sites with high prereproductive mortality, QTL alleles contributing to early reproduction resulted in greater fruit production, conferring selective advantage, whereas alleles contributing to later reproduction resulted in larger size and higher fitness in a continental site. This demonstrates how environmental variation leads to change in both QTL effect sizes and direction of selection on traits, justifying the persistence of allelic polymorphism at life history QTL across the species range.     

Analysis of the genetic architecture of maize kernel size traits by combined linkage and association mapping

Kernel size‐related traits are the most direct traits correlating with grain yield. The genetic basis of three kernel traits of maize, kernel length (KL), kernel width (KW) and kernel thickness (KT), was investigated in an association panel and a biparental population. A total of 21 single nucleotide polymorphisms (SNPs) were detected to be most significantly (P < 2.25 × 10^?6) associated with these three traits in the association panel under four environments. Furthermore, 50 quantitative trait loci (QTL) controlling these traits were detected in seven environments in the intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, of which eight were repetitively identified in at least three environments. Combining the two mapping populations revealed that 56 SNPs (P < 1 × 10^?3) fell within 18 of the QTL confidence intervals. According to the top significant SNPs, stable‐effect SNPs and the co‐localized SNPs by association analysis and linkage mapping, a total of 73 candidate genes were identified, regulating seed development. Additionally, seven miRNAs were found to situate within the linkage disequilibrium (LD) regions of the co‐localized SNPs, of which zma‐miR164e was demonstrated to cleave the mRNAs of Arabidopsis CUC1, CUC2 and NAC6 in vitro. Overexpression of zma‐miR164e resulted in the down‐regulation of these genes above and the failure of seed formation in Arabidopsis pods, with the increased branch number. These findings provide insights into the mechanism of seed development and the improvement of molecular marker‐assisted selection (MAS) for high‐yield breeding in maize.     

Color phenotypes are under similar genetic control in two distantly related species of Timema stick insect

Ecology and genetics are both of general interest to evolutionary biologists as they can influence the phenotypic and genetic response to selection. The stick insects Timema podura and Timema cristinae exhibit a green/melanistic body color polymorphism that is subject to different ecologically based selective regimes in the two species. Here, we describe aspects of the genetics of this color polymorphism in T. podura, and compare this to previous results in T. cristinae. We first show that similar color phenotypes of the two species cluster in phenotypic space. We then use genome‐wide association mapping to show that in both species, color is controlled by few loci, dominance relationships between color alleles are the same, and SNPs associated with color phenotypes colocalize to the same linkage group. Regions within this linkage group that harbor genetic variants associated with color exhibit elevated linkage disequilibrium relative to genome wide expectations, but more strongly so in T. cristinae. We use these results to discuss predictions regarding how the genetics of color could influence levels of phenotypic and genetic variation that segregate within and between populations of T. podura and T. cristinae, drawing parallels with other organisms.     

Genome‐wide association mapping and examination of possible maternal effect for the pace trait of horses

Previously, a single nucleotide polymorphism (SNP) related to gait type was identified at position 22 999 655 of chromosome 23 in the coding region of DMRT3 (DMRT3:Ser301Ter) by showing that a cytosine (C) to adenine (A) mutation of this SNP induced pace in the Icelandic horse. We investigated the effect of DMRT3:Ser301Ter on the gait of Hokkaido Native Horses, a Japanese native breed, and examined genetic factors other than DMRT3 by exploring genome‐wide SNPs related to gait determination. All animals exhibiting pace were AA for DMRT3:Ser301Ter, confirming the association of DMRT3:Ser301Ter with gait determination; however, 14.3% of the animals exhibiting trot also had AA for DMRT3:Ser301Ter, suggesting the presence of another factor(s) cooperatively working with DMRT3:Ser301Ter for gait determination. SNPs on chromosomes 13 and 23 were detected by genome‐wide association analysis (false discovery rate <0.05), although SNPs on chromosome 23 were all located in the vicinity of DMRT3:Ser301Ter, confirming the association with DMRT3. A genome‐wide association study targeting only animals with AA for DMRT3:Ser301Ter to examine genetic factors cooperatively working with DMRT3:Ser301Ter for gait determination suggested associations of 23 SNPs on six chromosomes. In a series of analyses of the effect of a maternal factor (dam's gait) on gait determination, the effect was suggested in comparison of the frequencies of exhibiting pace in gait checks in only two animal groups having dams with different DMRT3:Ser301Ter genotypes (P < 0.05), suggesting that the gait of the dam does not have a major effect on whether progeny homozygous for the DMRT3:Ser301Ter mutation will preferentially pace or trot.     

Meta‐analysis of genome‐wide association from genomic prediction models

Genome‐wide association (GWA) studies based on GBLUP models are a common practice in animal breeding. However, effect sizes of GWA tests are small, requiring larger sample sizes to enhance power of detection of rare variants. Because of difficulties in increasing sample size in animal populations, one alternative is to implement a meta‐analysis (MA), combining information and results from independent GWA studies. Although this methodology has been used widely in human genetics, implementation in animal breeding has been limited. Thus, we present methods to implement a MA of GWA, describing the proper approach to compute weights derived from multiple genomic evaluations based on animal‐centric GBLUP models. Application to real datasets shows that MA increases power of detection of associations in comparison with population‐level GWA, allowing for population structure and heterogeneity of variance components across populations to be accounted for. Another advantage of MA is that it does not require access to genotype data that is required for a joint analysis. Scripts related to the implementation of this approach, which consider the strength of association as well as the sign, are distributed and thus account for heterogeneity in association phase between QTL and SNPs. Thus, MA of GWA is an attractive alternative to summarizing results from multiple genomic studies, avoiding restrictions with genotype data sharing, definition of fixed effects and different scales of measurement of evaluated traits.