Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach of crop modelling and genome-wide association.

Understanding the genetic basis of phenotypic plasticity is crucial for predicting and managing climate change effects on wild plants and crops. Here, we combined crop modelling and quantitative genetics to study the genetic control of oil yield plasticity for multiple abiotic stresses in sunflower. First, we developed stress indicators to characterize 14 environments for three abiotic stresses (cold, drought and nitrogen) using the SUNFLO crop model and phenotypic variations of three commercial varieties. The computed plant stress indicators better explain yield variation than descriptors at the climatic or crop levels. In those environments, we observed oil yield of 317 sunflower hybrids and regressed it with three selected stress indicators. The slopes of cold stress norm reaction were used as plasticity phenotypes in the following genome-wide association study. Among the 65 534 tested Single Nucleotide Polymorphisms (SNPs), we identified nine quantitative trait loci controlling oil yield plasticity to cold stress. Associated single nucleotide polymorphisms are localized in genes previously shown to be involved in cold stress responses: oligopeptide transporters, lipid transfer protein, cystatin, alternative oxidase or root development. This novel approach opens new perspectives to identify genomic regions involved in genotype-by-environment interaction of a complex traits to multiple stresses in realistic natural or agronomical conditions.

Combined linkage and association mapping of flowering time in Sunflower (Helianthus annuus L.).

Association mapping and linkage mapping were used to identify quantitative trait loci (QTL) and/or causative mutations involved in the control of flowering time in cultivated sunflower Helianthus annuus. A panel of 384 inbred lines was phenotyped through testcrosses with two tester inbred lines across 15 location × year combinations. A recombinant inbred line (RIL) population comprising 273 lines was phenotyped both per se and through testcrosses with one or two testers in 16 location × year combinations. In the association mapping approach, kinship estimation using 5,923 single nucleotide polymorphisms was found to be the best covariate to correct for effects of panel structure. Linkage disequilibrium decay ranged from 0.08 to 0.26 cM for a threshold of 0.20, after correcting for structure effects, depending on the linkage group (LG) and the ancestry of inbred lines. A possible hitchhiking effect is hypothesized for LG10 and LG08. A total of 11 regions across 10 LGs were found to be associated with flowering time, and QTLs were mapped on 11 LGs in the RIL population. Whereas eight regions were demonstrated to be common between the two approaches, the linkage disequilibrium approach did not detect a documented QTL that was confirmed using the linkage mapping approach.

An SSR-SNP Linkage Map of the Parasitic Weed Orobanche cumana Wallr. Including a Gene for Plant Pigmentation.

Sunflower broomrape (Orobanche cumana Wallr.) is a holoparasitic plant that causes major yield losses to sunflower crops in the Old World. Efforts to understand how this parasitic weed recognizes and interacts with sunflowers are important for developing long-term genetic resistance strategies. However, such studies are hampered by the lack of genetic tools for O. cumana. The objectives of this research were to construct a genetic linkage map of this species using SSR and SNP markers, and mapping the Pg locus that is involved in plant pigmentation. The genetic map was developed from the progenies of a cross between the O. cumana inbred lines EK-12 and EK-A1, which originated from populations belonging to two distant and geographically separated gene pools identified in Spain. The inbred lines also differed in plant pigmentation, with EK-A1 lacking anthocyanin pigmentation (pgpg genotype). A genetic map comprising 26 SSR and 701 SNP markers was constructed, which displayed 19 linkage groups (LGs), corresponding to the 19 chromosome pairs of O. cumana. The total length of the map was 1795.7 cM, with an average distance between two adjacent positions of 2.5 cM and a maximum map distance of 41.9 cM. The Pg locus mapped to LG19 between the SNP markers OS02468 and OS01653 at 7.5 and 3.4 cM, respectively. This study constitutes the first linkage map and trait mapping study in Orobanche spp., laying a key foundation for further genome characterization and providing a basis for mapping additional traits such as those having a key role in parasitism.

Genomic regions involved in response to grain yield selection at high and low nitrogen fertilization in maize.

In order to validate the role of genomic regions involved in nitrogen use efficiency and detected in a population of recombinant inbred lines (RIL), we have applied from the same population a recurrent selection for adaptation to low N-input (N0) and to high N-input (N1). Variation of allele frequency at neutral marker during the two cycles of recurrent selection may provide information about markers linked to QTLs. Significant temporal variation of allele frequency was investigated using the test of Waples, which tests the hypothesis of genetic drift versus selection. Most genomic regions (12/19) responding to selection were detected for selection at high N-input and only two were common to selection at high and low N-inputs. This was consistent with the greater grain yield response to selection observed for the population selected under high N-input compared with the population selected under low N-input, when they were evaluated at high N-fertilization. In contrast, when they were evaluated at low N-input both types of selection gave similar yield. As was expected, in the first cycle we observed selection of markers linked to grain yield QTLs. In the course of the second cycle three situations were observed: the confirmation of most regions already selected in C1 including all C1 regions overlapping with grain yield QTLs; the non-confirmation of some C1 regions (2/9); and the identification of new genomic zones (10/17). The detected marker-QTL associations revealed the consistency of the involvement of some traits, such as root architecture and glutamine synthetase activity, which would be of major importance for grain yield setting whatever the nitrogen fertilization.

Modelling postsilking nitrogen fluxes in maize (Zea mays) using 15N-labelling field experiments.

In maize (Zea mays), nitrogen (N) remobilization and postflowering N uptake are two processes that provide amino acids for grain protein synthesis. To study the way in which N is allocated to the grain and to the stover, two different 15N-labelling techniques were developed. 15NO(3-) was provided to the soil either at the beginning of stem elongation or after silking. The distribution of 15N in the stover and in the grain was monitored by calculating relative 15N-specific allocation (RSA). A nearly linear relationship between the RSA of the kernels and the RSA of the stover was found as a result of two simultaneous N fluxes: N remobilization from the stover to the grain, and N allocation to the stover and to the grain originating from N uptake. By modelling the 15N fluxes, it was possible to demonstrate that, as a consequence of protein turnover, a large proportion of the amino acids synthesized from the N taken up after silking were integrated into the proteins of the stover, and these proteins were further hydrolysed to provide N to the grain.

Two cytosolic glutamine synthetase isoforms of maize are specifically involved in the control of grain production.

The roles of two cytosolic maize glutamine synthetase isoenzymes (GS1), products of the Gln1-3 and Gln1-4 genes, were investigated by examining the impact of knockout mutations on kernel yield. In the gln1-3 and gln1-4 single mutants and the gln1-3 gln1-4 double mutant, GS mRNA expression was impaired, resulting in reduced GS1 protein and activity. The gln1-4 phenotype displayed reduced kernel size and gln1-3 reduced kernel number, with both phenotypes displayed in gln1-3 gln1-4. However, at maturity, shoot biomass production was not modified in either the single mutants or double mutants, suggesting a specific impact on grain production in both mutants. Asn increased in the leaves of the mutants during grain filling, indicating that it probably accumulates to circumvent ammonium buildup resulting from lower GS1 activity. Phloem sap analysis revealed that unlike Gln, Asn is not efficiently transported to developing kernels, apparently causing reduced kernel production. When Gln1-3 was overexpressed constitutively in leaves, kernel number increased by 30%, providing further evidence that GS1-3 plays a major role in kernel yield. Cytoimmunochemistry and in situ hybridization revealed that GS1-3 is present in mesophyll cells, whereas GS1-4 is specifically localized in the bundle sheath cells. The two GS1 isoenzymes play nonredundant roles with respect to their tissue-specific localization.