Genome-Wide Association Study in Arabidopsis thaliana of Natural Variation in Seed Oil Melting Point: A Widespread Adaptive Trait in Plants

Extensive simulations assess the performance of genome-wide association mapping in various Saccharomyces cerevisiae subpopulations

With the advent of high throughput sequencing technologies, genome-wide association studies (GWAS) have become a powerful paradigm for dissecting the genetic origins of the observed phenotypic variation. We recently completely sequenced the genome of 1011 Saccharomyces cerevisiae isolates, laying a strong foundation for GWAS. To assess the feasibility and the limits of this approach, we performed extensive simulations using five selected subpopulations as well as the total set of 1011 genomes. We measured the ability to detect the causal genetic variants involved in Mendelian and more complex traits using a linear mixed model approach. The results showed that population structure is well accounted for and is not the main problem when the sample size is high enough. While the genetic determinant of a Mendelian trait is easily mapped in all studied subpopulations, discrepancies are seen between datasets when performing GWAS on a complex trait in terms of detection, false positive and false negative rate. Finally, we performed GWAS on the different defined subpopulations using a real quantitative trait (resistance to copper sulfate) and showed the feasibility of this approach. The performance of each dataset depends simultaneously on several factors such as sample size, relatedness and population evolutionary history.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Interactions between genetics and environment shape Camelina seed oil composition

BACKGROUND

Camelina sativa (gold-of-pleasure) is a traditional European oilseed crop and emerging biofuel source with high levels of desirable fatty acids. A twentieth century germplasm bottleneck depleted genetic diversity in the crop, leading to recent interest in using wild relatives for crop improvement. However, little is known about seed oil content and genetic diversity in wild Camelina species.

RESULTS

We used gas chromatography, environmental niche assessment, and genotyping-by-sequencing to assess seed fatty acid composition, environmental distributions, and population structure in C. sativa and four congeners, with a primary focus on the crop's wild progenitor, C. microcarpa. Fatty acid composition differed significantly between Camelina species, which occur in largely non-overlapping environments. The crop progenitor comprises three genetic subpopulations with discrete fatty acid compositions. Environment, subpopulation, and population-by-environment interactions were all important predictors for seed oil in these wild populations. A complementary growth chamber experiment using C. sativa confirmed that growing conditions can dramatically affect both oil quantity and fatty acid composition in Camelina.

CONCLUSIONS

Genetics, environmental conditions, and genotype-by-environment interactions all contribute to fatty acid variation in Camelina species. These insights suggest careful breeding may overcome the unfavorable FA compositions in oilseed crops that are predicted with warming climates.

Detection of significant SNP associated with production and oil quality traits in interspecific oil palm hybrids using RARSeq

A RARSeq based Association mapping study was performed in a population of 104 Elaeis oleifera x E. guineensis hybrids of five origins with the aim of finding functional markers associated to six productive and 19 oil quality traits. For this purpose mRNA of each genotype was isolated and double stranded cDNA was synthesized. Following digestion with two restriction enzymes and adapter ligation, a size selected pool of barcoded amplicons was produced and sequenced using Illumina MiSeq. The obtained sequences were processed with a "snakemake" pipeline, filtered and missing values were imputed. For all traits except two significant effects of the origin was observed. Genetic diversity analyses revealed high variability within origins and an excess of heterozygosity in the population. Two GLM models with Q matrix or PCA matrix as covariates and two MLM models incorporating in addition a Kinship matrix were tested for genotype-phenotype associations using GAPIT software. Using unadjusted p values (< 0.01) 78 potential associations were detected involving 25 SNP and 20 traits. When applying FDR multiple testing with p <  0.05, 25 significant associations remained involving eight SNP and six quality traits. Four SNP were located in genes with a potential relevant biological meaning.

Natural variation in acyl editing is a determinant of seed storage oil composition

Seeds exhibit wide variation in the fatty acid composition of their storage oil. However, the genetic basis of this variation is only partially understood. Here we have used a multi-parent advanced generation inter-cross (MAGIC) population to study the genetic control of fatty acid chain length in Arabidopsis thaliana seed oil. We mapped four quantitative trait loci (QTL) for the quantity of the major very long chain fatty acid species 11-eicosenoic acid (20:1), using multiple QTL modelling. Surprisingly, the main-effect QTL does not coincide with FATTY ACID ELONGASE 1 and a parallel genome wide association study suggested that LYSOPHOSPHATIDYLCHOLINE ACYLTRANSFERASE 2 (LPCAT2) is a candidate for this QTL. Regression analysis also suggested that LPCAT2 expression and 20:1 content in seeds of the 19 MAGIC founder accessions are related. LPCAT is a key component of the Lands cycle; an acyl editing pathway that enables acyl-exchange between the acyl-Coenzyme A and phosphatidylcholine precursor pools used for microsomal fatty acid elongation and desaturation, respectively. We Mendelianised the main-effect QTL using biparental chromosome segment substitution lines and carried out complementation tests to show that a single cis-acting polymorphism in the LPCAT2 promoter causes the variation in seed 20:1 content, by altering the LPCAT2 expression level and total LPCAT activity in developing siliques. Our work establishes that oilseed species exhibit natural variation in the enzymic capacity for acyl editing and this contributes to the genetic control of storage oil composition.

The genetic architecture of amylose biosynthesis in maize kernel

Starch is the most abundant storage carbohydrate in maize kernel. The content of amylose and amylopectin confers unique properties in food processing and industrial application. Thus, the resurgent interest has been switched to the study of individual amylose or amylopectin rather than total starch, whereas the enzymatic machinery for amylose synthesis remains elusive. We took advantage of the phenotype of amylose content and the genotype of 9,007,194 single nucleotide polymorphisms from 464 inbred maize lines. The genome-wide association study identified 27 associated loci involving 39 candidate genes that were linked to amylose content including transcription factors, glycosyltransferases, glycosidases, as well as hydrolases. Except the waxy gene that encodes the granule-bound starch synthase, the remaining candidate genes were located in the upstream pathway of amylose synthesis, while the downstream members were already known from prior studies. The linked candidate genes could be transferred to manipulate amylose content and thus add value to maize kernel in the breeding programme.

Genome-wide association study reveals candidate genes influencing lipids and diterpenes contents in Coffea arabica L

Lipids, including the diterpenes cafestol and kahweol, are key compounds that contribute to the quality of coffee beverages. We determined total lipid content and cafestol and kahweol concentrations in green beans and genotyped 107 Coffea arabica accessions, including wild genotypes from the historical FAO collection from Ethiopia. A genome-wide association study was performed to identify genomic regions associated with lipid, cafestol and kahweol contents and cafestol/kahweol ratio. Using the diploid Coffea canephora genome as a reference, we identified 6,696 SNPs. Population structure analyses suggested the presence of two to three groups (K = 2 and K = 3) corresponding to the east and west sides of the Great Rift Valley and an additional group formed by wild accessions collected in western forests. We identified 5 SNPs associated with lipid content, 4 with cafestol, 3 with kahweol and 9 with cafestol/kahweol ratio. Most of these SNPs are located inside or near candidate genes related to metabolic pathways of these chemical compounds in coffee beans. In addition, three trait-associated SNPs showed evidence of directional selection among cultivated and wild coffee accessions. Our results also confirm a great allelic richness in wild accessions from Ethiopia, especially in accessions originating from forests in the west side of the Great Rift Valley.

Genome-wide Association Study Reveals that the Aquaporin NIP1;1 Contributes to Variation in Hydrogen Peroxide Sensitivity in Arabidopsis thaliana

Hydrogen peroxide (H₂O₂) is a reactive oxygen species that affects cell signaling in various plant defense responses and induces programmed cell death. To identify the new components associated with H₂O₂ signaling and tolerance, we conducted a genome-wide association study (GWAS) on the root growth of 133 Arabidopsis thaliana accessions grown in the presence of toxic H₂O₂ levels. The most significant SNPs were associated with a cluster of chromosome 4 genes encoding an aquaporin NODULIN 26-LIKE INTRINSIC PROTEIN 1; 1 (NIP1;1), an NB-ARC domain-containing disease resistance protein (AT4G19050), and a putative membrane lipoprotein (AT4G19070). The expression level of NIP1;1 was relatively high in A. thaliana accessions sensitive to H₂O₂. Additionally, overexpression of NIP1;1 in a tolerant accession (e.g., Col-0) increased the sensitivity of transgenic plants to H₂O₂. An in planta β-glucuronidase reporter assay revealed that variations in the NIP1;1 promoter were responsible for the differences of its expression level in H₂O₂-tolerant and -sensitive accessions. Cell death was extensive and H₂O₂ levels were high in the roots of H₂O₂-sensitive and NIP1;1-overexpressing accessions. Together, our results indicate that the aquaporin NIP1;1 is a key determinant of the sensitivity of A. thaliana to H₂O₂, and contributes to the phenotypic variations detected by our GWAS.

Genome Wide Analysis of Fatty Acid Desaturation and Its Response to Temperature

Plants modify the polyunsaturated fatty acid content of their membrane and storage lipids in order to adapt to changes in temperature. In developing seeds, this response is largely controlled by the activities of the microsomal ω-6 and ω-3 fatty acid desaturases, FAD2 and FAD3. Although temperature regulation of desaturation has been studied at the molecular and biochemical levels, the genetic control of this trait is poorly understood. Here, we have characterized the response of Arabidopsis (Arabidopsis thaliana) seed lipids to variation in ambient temperature and found that heat inhibits both ω-6 and ω-3 desaturation in phosphatidylcholine, leading to a proportional change in triacylglycerol composition. Analysis of the 19 parental accessions of the multiparent advanced generation intercross (MAGIC) population showed that significant natural variation exists in the temperature responsiveness of ω-6 desaturation. A combination of quantitative trait locus (QTL) analysis and genome-wide association studies (GWAS) using the MAGIC population suggests that ω-6 desaturation is largely controlled by cis-acting sequence variants in the FAD2 5' untranslated region intron that determine the expression level of the gene. However, the temperature responsiveness of ω-6 desaturation is controlled by a separate QTL on chromosome 2. The identity of this locus is unknown, but genome-wide association studies identified potentially causal sequence variants within ∼40 genes in an ∼450-kb region of the QTL.

Correlations of Genotype with Climate Parameters Suggest Caenorhabditis elegans Niche Adaptations

Species inhabit a variety of environmental niches, and the adaptation to a particular niche is often controlled by genetic factors, including gene-by-environment interactions. The genes that vary in order to regulate the ability to colonize a niche are often difficult to identify, especially in the context of complex ecological systems and in experimentally uncontrolled natural environments. Quantitative genetic approaches provide an opportunity to investigate correlations between genetic factors and environmental parameters that might define a niche. Previously, we have shown how a collection of 208 whole-genome sequenced wild Caenorhabditis elegans can facilitate association mapping approaches. To correlate climate parameters with the variation found in this collection of wild strains, we used geographic data to exhaustively curate daily weather measurements in short-term (3 month), middle-term (one year), and long-term (three year) durations surrounding the date of strain isolation. These climate parameters were used as quantitative traits in association mapping approaches, where we identified 11 quantitative trait loci (QTL) for three climatic variables: elevation, relative humidity, and average temperature. We then narrowed the genomic interval of interest to identify gene candidates with variants potentially underlying phenotypic differences. Additionally, we performed two-strain competition assays at high and low temperatures to validate a QTL that could underlie adaptation to temperature and found suggestive evidence supporting that hypothesis.