Search

Sunflower is a hybrid crop that is considered moderately drought-tolerant and adapted to new cropping systems required for the agro-ecological transition. Here, we studied the impact of hybridity status (hybrids vs. inbred lines) on the responses to drought at the molecular and eco-physiological level exploiting publicly available datasets. Eco-physiological traits and leaf proteomes were measured in eight inbred lines and their sixteen hybrids grown in the high-throughput phenotyping platform Phenotoul-Heliaphen. Hybrids and parental lines showed different growth strategies: hybrids grew faster in the absence of water constraint and arrested their growth more abruptly than inbred lines when subjected to water deficit. We identified 471 differentially accumulated proteins, of which 256 were regulated by drought. The amplitude of up- and downregulations was greater in hybrids than in inbred lines. Our results show that hybrids respond more strongly to water deficit at the molecular and eco-physiological levels. Because of presence/absence polymorphism, hybrids potentially contain more genes than their parental inbred lines. We propose that detrimental homozygous mutations and the lower number of genes in inbred lines lead to a constitutive defense mechanism that may explain the lower growth of inbred lines under well-watered conditions and their lower reactivity to water deficit.

Genomic prediction is a useful tool for plant and animal breeding programs and is starting to be used to predict human diseases as well. A shortcoming that slows down the genomic selection deployment is that the accuracy of the prediction is not known a priori. We propose EthAcc (Estimated THeoretical ACCuracy) as a method for estimating the accuracy given a training set that is genotyped and phenotyped. EthAcc is based on a causal quantitative trait loci model estimated by a genome-wide association study. This estimated causal model is crucial; therefore, we compared different methods to find the one yielding the best EthAcc. The multilocus mixed model was found to perform the best. We compared EthAcc to accuracy estimators that can be derived via a mixed marker model. We showed that EthAcc is the only approach to correctly estimate the accuracy. Moreover, in case of a structured population, in accordance with the achieved accuracy, EthAcc showed that the biggest training set is not always better than a smaller and closer training set. We then performed training set optimization with EthAcc and compared it to CDmean. EthAcc outperformed CDmean on real datasets from sugar beet, maize, and wheat. Nonetheless, its performance was mainly due to the use of an optimal but inaccessible set as a start of the optimization algorithm. EthAcc's precision and algorithm issues prevent it from reaching a good training set with a random start. Despite this drawback, we demonstrated that a substantial gain in accuracy can be obtained by performing training set optimization.

Heliaphen is an outdoor platform designed for high-throughput phenotyping. It allows the automated management of drought scenarios and monitoring of plants throughout their lifecycles. A robot moving between plants growing in 15-L pots monitors the plant water status and phenotypes the leaf or whole-plant morphology. From these measurements, we can compute more complex traits, such as leaf expansion (LE) or transpiration rate (TR) in response to water deficit. Here, we illustrate the capabilities of the platform with two practical cases in sunflower (Helianthus annuus): a genetic and genomic study of the response of yield-related traits to drought, and a modeling study using measured parameters as inputs for a crop simulation. For the genetic study, classical measurements of thousand-kernel weight (TKW) were performed on a biparental population under automatically managed drought stress and control conditions. These data were used for an association study, which identified five genetic markers of the TKW drought response. A complementary transcriptomic analysis identified candidate genes associated with these markers that were differentially expressed in the parental backgrounds in drought conditions. For the simulation study, we used a crop simulation model to predict the impact on crop yield of two traits measured on the platform (LE and TR) for a large number of environments. We conducted simulations in 42 contrasting locations across Europe using 21 years of climate data. We defined the pattern of abiotic stresses occurring at the continental scale and identified ideotypes (i.e., genotypes with specific trait values) that are more adapted to specific environment types. This study exemplifies how phenotyping platforms can assist the identification of the genetic architecture controlling complex response traits and facilitate the estimation of ecophysiological model parameters to define ideotypes adapted to different environmental conditions.

Prediction of hybrid performance using incomplete factorial mating designs is widely used in breeding programs including different heterotic groups. Based on the general combining ability (GCA) of the parents, predictions are accurate only if the genetic variance resulting from the specific combining ability is small and both parents have phenotyped descendants. Genomic selection (GS) can predict performance using a model trained on both phenotyped and genotyped hybrids that do not necessarily include all hybrid parents. Therefore, GS could overcome the issue of unknown parent GCA. Here, we compared the accuracy of classical GCA-based and genomic predictions for oil content of sunflower seeds using several GS models. Our study involved 452 sunflower hybrids from an incomplete factorial design of 36 female and 36 male lines. Re-sequencing of parental lines allowed to identify 468,194 non-redundant SNPs and to infer the hybrid genotypes. Oil content was observed in a multi-environment trial (MET) over 3 years, leading to nine different environments. We compared GCA-based model to different GS models including female and male genomic kinships with the addition of the female-by-male interaction genomic kinship, the use of functional knowledge as SNPs in genes of oil metabolic pathways, and with epistasis modeling. When both parents have descendants in the training set, the predictive ability was high even for GCA-based prediction, with an average MET value of 0.782. GS performed slightly better (+0.2%). Neither the inclusion of the female-by-male interaction, nor functional knowledge of oil metabolism, nor epistasis modeling improved the GS accuracy. GS greatly improved predictive ability when one or both parents were untested in the training set, increasing GCA-based predictive ability by 10.4% from 0.575 to 0.635 in the MET. In this scenario, performing GS only considering SNPs in oil metabolic pathways did not improve whole genome GS prediction but increased GCA-based prediction ability by 6.4%. Our results show that GS is a major improvement to breeding efficiency compared to the classical GCA modeling when either one or both parents are not well-characterized. This finding could therefore accelerate breeding through reducing phenotyping efforts and more effectively targeting for the most promising crosses.

Genomic selection is focused on prediction of breeding values of selection candidates by means of high density of markers. It relies on the assumption that all quantitative trait loci (QTLs) tend to be in strong linkage disequilibrium (LD) with at least one marker. In this context, we present theoretical results regarding the accuracy of genomic selection, i.e., the correlation between predicted and true breeding values. Typically, for individuals (so-called test individuals), breeding values are predicted by means of markers, using marker effects estimated by fitting a ridge regression model to a set of training individuals. We present a theoretical expression for the accuracy; this expression is suitable for any configurations of LD between QTLs and markers. We also introduce a new accuracy proxy that is free of the QTL parameters and easily computable; it outperforms the proxies suggested in the literature, in particular, those based on an estimated effective number of independent loci (Me). The theoretical formula, the new proxy, and existing proxies were compared for simulated data, and the results point to the validity of our approach. The calculations were also illustrated on a new perennial ryegrass set (367 individuals) genotyped for 24,957 single nucleotide polymorphisms (SNPs). In this case, most of the proxies studied yielded similar results because of the lack of markers for coverage of the entire genome (2.7 Gb).

Adaptation to pollinators is a key factor of diversification in angiosperms. The Caribbean sister genera Rhytidophyllum and Gesneria present an important diversification of floral characters. Most of their species can be divided in two major pollination syndromes. Large-open flowers with pale colours and great amount of nectar represent the generalist syndrome, while the hummingbird-specialist syndrome corresponds to red tubular flowers with a less important nectar volume. Repeated convergent evolution toward the generalist syndrome in this group suggests that such transitions rely on few genes of moderate to large effect. To test this hypothesis, we built a linkage map and performed a QTL detection for divergent pollination syndrome traits by crossing one specimen of the generalist species Rhytidophyllum auriculatum with one specimen of the hummingbird pollinated R. rupincola. Using geometric morphometrics and univariate traits measurements, we found that floral shape among the second-generation hybrids is correlated with morphological variation observed between generalist and hummingbird-specialist species at the genus level. The QTL analysis showed that colour and nectar volume variation between syndromes involve each one major QTL while floral shape has a more complex genetic basis and rely on few genes of moderate effect. Finally, we did not detect any genetic linkage between the QTLs underlying those traits. This genetic independence of traits could have facilitated evolution toward optimal syndromes.

Bacterial wilt caused by the soilborne bacterium Ralstonia solanacearum attacks hundreds of plant species, including many agriculturally important crops. Natural resistance to this disease has been found in some species and is usually inherited as a polygenic trait. In tomato, a model crop plant, genetic analysis previously revealed the involvement of several QTL (quantitative trait loci) controlling resistance and, in all of these studies with different strains of the pathogen, loci on chromosome 6 played the predominant role in controlling this trait. Using quantitative data collected from a greenhouse test F3 population, we identified a new locus on chromosome 12 that appears to be active specifically against a race 1 biovar 3 Pss4 bacterial strain endemic to Taiwan. Chromosome 6 still contributes significantly to the control of the resistance, and weaker associations of the trait to other regions of the genome are observed. These results are discussed in the context of current molecular knowledge about the strain specificity of disease resistance genes.

Modern technologies and especially next generation sequencing facilities are giving a cheaper access to genotype and genomic data measured on the same sample at once. This creates an ideal situation for multifactorial experiments designed to infer gene regulatory networks. The fifth "Dialogue for Reverse Engineering Assessments and Methods" (DREAM5) challenges are aimed at assessing methods and associated algorithms devoted to the inference of biological networks. Challenge 3 on "Systems Genetics" proposed to infer causal gene regulatory networks from different genetical genomics data sets. We investigated a wide panel of methods ranging from Bayesian networks to penalised linear regressions to analyse such data, and proposed a simple yet very powerful meta-analysis, which combines these inference methods. We present results of the Challenge as well as more in-depth analysis of predicted networks in terms of structure and reliability. The developed meta-analysis was ranked first among the 16 teams participating in Challenge 3A. It paves the way for future extensions of our inference method and more accurate gene network estimates in the context of genetical genomics.