Confirmation of Fusarium root rot resistance QTL Fsp-Ps 2.1 of pea under controlled conditions.

BACKGROUND: Dry pea production has increased substantially in North America over the last few decades. With this expansion, significant yield losses have been attributed to an escalation in Fusarium root rots in pea fields. Among the most significant rot rotting pathogenic fungal species, Fusarium solani fsp. pisi (Fsp) is one of the main causal agents of root rot of pea. High levels of partial resistance to Fsp has been identified in plant genetic resources. Genetic resistance offers one of the best solutions to control this root rotting fungus. A recombinant inbred population segregating for high levels of partial resistance, previously single nucleotide polymorphism (SNP) genotyped using genotyping-by-sequencing, was phenotyped for disease reaction in replicated and repeated greenhouse trials. Composite interval mapping was deployed to identify resistance-associated quantitative trait loci (QTL). RESULTS: Three QTL were identified using three disease reaction criteria: root disease severity, ratios of diseased vs. healthy shoot heights and dry plant weights under controlled conditions using pure cultures of Fusarium solani fsp. pisi. One QTL Fsp-Ps 2.1 explains 44.4-53.4% of the variance with a narrow confidence interval of 1.2 cM. The second and third QTL Fsp-Ps3.2 and Fsp-Ps3.3 are closely linked and explain only 3.6-4.6% of the variance. All of the alleles are contributed by the resistant parent PI 180693. CONCLUSION: With the confirmation of Fsp-Ps 2.1 now in two RIL populations, SNPs associated with this region make a good target for marker-assisted selection in pea breeding programs to obtain high levels of partial resistance to Fusarium root rot caused by Fusarium solani fsp. pisi.

A local score approach improves GWAS resolution and detects minor QTL: application to Medicago truncatula quantitative disease resistance to multiple Aphanomyces euteiches isolates.

Quantitative trait loci (QTL) with small effects, which are pervasive in quantitative phenotypic variation, are difficult to detect in genome-wide association studies (GWAS). To improve their detection, we propose to use a local score approach that accounts for the surrounding signal due to linkage disequilibrium, by accumulating association signals from contiguous single markers. Simulations revealed that, in a GWAS context with high marker density, the local score approach outperforms single SNP p-value-based tests for detecting minor QTL (heritability of 5-10%) and is competitive with regard to alternative methods, which also aggregate p-values. Using more than five million SNPs, this approach was applied to identify loci involved in Quantitative Disease Resistance (QDR) to different isolates of the plant root rot pathogen Aphanomyces euteiches, from a GWAS performed on a collection of 174 accessions of the model legume Medicago truncatula. We refined the position of a previously reported major locus, underlying MYB/NB-ARC/tyrosine kinase candidate genes conferring resistance to two closely related A. euteiches isolates belonging to pea pathotype I. We also discovered a diversity of minor resistance QTL, not detected using p-value-based tests, some of which being putatively shared in response to pea (pathotype I and III) and/or alfalfa (race 1 and 2) isolates. Candidate genes underlying these QTL suggest pathogen effector recognition and plant proteasome as key functions associated with M. truncatula resistance to A. euteiches. GWAS on any organism can benefit from the local score approach to uncover many weak-effect QTL.

Development of two major resources for pea genomics: the GenoPea 13.2K SNP Array and a high-density, high-resolution consensus genetic map.

Single nucleotide polymorphism (SNP) arrays represent important genotyping tools for innovative strategies in both basic research and applied breeding. Pea is an important food, feed and sustainable crop with a large (about 4.45 Gbp) but not yet available genome sequence. In the present study, 12 pea recombinant inbred line populations were genotyped using the newly developed GenoPea 13.2K SNP Array. Individual and consensus genetic maps were built providing insights into the structure and organization of the pea genome. Largely collinear genetic maps of 3918-8503 SNPs were obtained from all mapping populations, and only two of these exhibited putative chromosomal rearrangement signatures. Similar distortion patterns in different populations were noted. A total of 12 802 transcript-derived SNP markers placed on a 15 079-marker high-density, high-resolution consensus map allowed the identification of ohnologue-rich regions within the pea genome and the localization of local duplicates. Dense syntenic networks with sequenced legume genomes were further established, paving the way for the identification of the molecular bases of important agronomic traits segregating in the mapping populations. The information gained on the structure and organization of the genome from this research will undoubtedly contribute to the understanding of the evolution of the pea genome and to its assembly. The GenoPea 13.2K SNP Array and individual and consensus genetic maps are valuable genomic tools for plant scientists to strengthen pea as a model for genetics and physiology and enhance breeding.

SNP discovery and genetic mapping using genotyping by sequencing of whole genome genomic DNA from a pea RIL population.

BACKGROUND: Progress in genetics and breeding in pea still suffers from the limited availability of molecular resources. SNP markers that can be identified through affordable sequencing processes, without the need for prior genome reduction or a reference genome to assemble sequencing data would allow the discovery and genetic mapping of thousands of molecular markers. Such an approach could significantly speed up genetic studies and marker assisted breeding for non-model species. RESULTS: A total of 419,024 SNPs were discovered using HiSeq whole genome sequencing of four pea lines, followed by direct identification of SNP markers without assembly using the discoSnp tool. Subsequent filtering led to the identification of 131,850 highly designable SNPs, polymorphic between at least two of the four pea lines. A subset of 64,754 SNPs was called and genotyped by short read sequencing on a subpopulation of 48 RILs from the cross 'Baccara' x 'PI180693'. This data was used to construct a WGGBS-derived pea genetic map comprising 64,263 markers. This map is collinear with previous pea consensus maps and therefore with the Medicago truncatula genome. Sequencing of four additional pea lines showed that 33 % to 64 % of the mapped SNPs, depending on the pairs of lines considered, are polymorphic and can therefore be useful in other crosses. The subsequent genotyping of a subset of 1000 SNPs, chosen for their mapping positions using a KASP™ assay, showed that almost all generated SNPs are highly designable and that most (95 %) deliver highly qualitative genotyping results. Using rather low sequencing coverages in SNP discovery and in SNP inferring did not hinder the identification of hundreds of thousands of high quality SNPs. CONCLUSIONS: The development and optimization of appropriate tools in SNP discovery and genetic mapping have allowed us to make available a massive new genomic resource in pea. It will be useful for both fine mapping within chosen QTL confidence intervals and marker assisted breeding for important traits in pea improvement.

Genome-wide association mapping of partial resistance to Aphanomyces euteiches in pea.

BACKGROUND: Genome-wide association (GWA) mapping has recently emerged as a valuable approach for refining the genetic basis of polygenic resistance to plant diseases, which are increasingly used in integrated strategies for durable crop protection. Aphanomyces euteiches is a soil-borne pathogen of pea and other legumes worldwide, which causes yield-damaging root rot. Linkage mapping studies reported quantitative trait loci (QTL) controlling resistance to A. euteiches in pea. However the confidence intervals (CIs) of these QTL remained large and were often linked to undesirable alleles, which limited their application in breeding. The aim of this study was to use a GWA approach to validate and refine CIs of the previously reported Aphanomyces resistance QTL, as well as identify new resistance loci. METHODS: A pea-Aphanomyces collection of 175 pea lines, enriched in germplasm derived from previously studied resistant sources, was evaluated for resistance to A. euteiches in field infested nurseries in nine environments and with two strains in climatic chambers. The collection was genotyped using 13,204 SNPs from the recently developed GenoPea Infinium® BeadChip. RESULTS: GWA analysis detected a total of 52 QTL of small size-intervals associated with resistance to A. euteiches, using the recently developed Multi-Locus Mixed Model. The analysis validated six of the seven previously reported main Aphanomyces resistance QTL and detected novel resistance loci. It also provided marker haplotypes at 14 consistent QTL regions associated with increased resistance and highlighted accumulation of favourable haplotypes in the most resistant lines. Previous linkages between resistance alleles and undesired late-flowering alleles for dry pea breeding were mostly confirmed, but the linkage between loci controlling resistance and coloured flowers was broken due to the high resolution of the analysis. A high proportion of the putative candidate genes underlying resistance loci encoded stress-related proteins and others suggested that the QTL are involved in diverse functions. CONCLUSION: This study provides valuable markers, marker haplotypes and germplasm lines to increase levels of partial resistance to A. euteiches in pea breeding.

Transcriptome sequencing for high throughput SNP development and genetic mapping in Pea.

BACKGROUND: Pea has a complex genome of 4.3 Gb for which only limited genomic resources are available to date. Although SNP markers are now highly valuable for research and modern breeding, only a few are described and used in pea for genetic diversity and linkage analysis. RESULTS: We developed a large resource by cDNA sequencing of 8 genotypes representative of modern breeding material using the Roche 454 technology, combining both long reads (400 bp) and high coverage (3.8 million reads, reaching a total of 1,369 megabases). Sequencing data were assembled and generated a 68 K unigene set, from which 41 K were annotated from their best blast hit against the model species Medicago truncatula. Annotated contigs showed an even distribution along M. truncatula pseudochromosomes, suggesting a good representation of the pea genome. 10 K pea contigs were found to be polymorphic among the genetic material surveyed, corresponding to 35 K SNPs.We validated a subset of 1538 SNPs through the GoldenGate assay, proving their ability to structure a diversity panel of breeding germplasm. Among them, 1340 were genetically mapped and used to build a new consensus map comprising a total of 2070 markers. Based on blast analysis, we could establish 1252 bridges between our pea consensus map and the pseudochromosomes of M. truncatula, which provides new insight on synteny between the two species. CONCLUSIONS: Our approach created significant new resources in pea, i.e. the most comprehensive genetic map to date tightly linked to the model species M. truncatula and a large SNP resource for both academic research and breeding.

High-density genome-wide association mapping implicates an F-box encoding gene in Medicago truncatula resistance to Aphanomyces euteiches.

• The use of quantitative disease resistance (QDR) is a promising strategy for promoting durable resistance to plant pathogens, but genes involved in QDR are largely unknown. To identify genetic components and accelerate improvement of QDR in legumes to the root pathogen Aphanomyces euteiches, we took advantage of both the recently generated massive genomic data for Medicago truncatula and natural variation of this model legume. • A high-density (≈5.1 million single nucleotide polymorphisms (SNPs)) genome-wide association study (GWAS) was performed with both in vitro and glasshouse phenotyping data collected for 179 lines. • GWAS identified several candidate genes and pinpointed two independent major loci on the top of chromosome 3 that were detected in both phenotyping methods. Candidate SNPs in the most significant locus (σ(A)²= 23%) were in the promoter and coding regions of an F-box protein coding gene. Subsequent qRT-PCR and bioinformatic analyses performed on 20 lines demonstrated that resistance is associated with mutations directly affecting the interaction domain of the F-box protein rather than gene expression. • These results refine the position of previously identified QTL to specific candidate genes, suggest potential molecular mechanisms, and identify new loci explaining QDR against A. euteiches.

QTL meta-analysis provides a comprehensive view of loci controlling partial resistance to Aphanomyces euteiches in four sources of resistance in pea.

BACKGROUND: Development of durable plant genetic resistance to pathogens through strategies of QTL pyramiding and diversification requires in depth knowledge of polygenic resistance within the available germplasm. Polygenic partial resistance to Aphanomyces root rot, caused by Aphanomyces euteiches, one of the most damaging pathogens of pea worldwide, was previously dissected in individual mapping populations. However, there are no data available regarding the diversity of the resistance QTL across a broader collection of pea germplasm. In this study, we performed a meta-analysis of Aphanomyces root rot resistance QTL in the four main sources of resistance in pea and compared their genomic localization with genes/QTL controlling morphological or phenological traits and with putative candidate genes. RESULTS: Meta-analysis, conducted using 244 individual QTL reported previously in three mapping populations (Puget x 90-2079, Baccara x PI180693 and Baccara x 552) and in a fourth mapping population in this study (DSP x 90-2131), resulted in the identification of 27 meta-QTL for resistance to A. euteiches. Confidence intervals of meta-QTL were, on average, reduced four-fold compared to mean confidence intervals of individual QTL. Eleven consistent meta-QTL, which highlight seven highly consistent genomic regions, were identified. Few meta-QTL specificities were observed among mapping populations, suggesting that sources of resistance are not independent. Seven resistance meta-QTL, including six of the highly consistent genomic regions, co-localized with six of the meta-QTL identified in this study for earliness and plant height and with three morphological genes (Af, A, R). Alleles contributing to the resistance were often associated with undesirable alleles for dry pea breeding. Candidate genes underlying six main meta-QTL regions were identified using colinearity between the pea and Medicago truncatula genomes. CONCLUSIONS: QTL meta-analysis provided an overview of the moderately low diversity of loci controlling partial resistance to A. euteiches in four main sources of resistance in pea. Seven highly consistent genomic regions with potential use in marker-assisted-selection were identified. Confidence intervals at several main QTL regions were reduced and co-segregation among resistance and morphological/phenological alleles was identified. Further work will be required to identify the best combinations of QTL for durably increasing partial resistance to A. euteiches.

Natural diversity in the model legume Medicago truncatula allows identifying distinct genetic mechanisms conferring partial resistance to Verticillium wilt.

Verticillium wilt is a major threat to alfalfa (Medicago sativa) and many other crops. The model legume Medicago truncatula was used as a host for studying resistance and susceptibility to Verticillium albo-atrum. In addition to presenting well-established genetic resources, this wild plant species enables to investigate biodiversity of the response to the pathogen and putative crosstalk between disease and symbiosis. Symptom scoring after root inoculation and modelling of disease curves allowed assessing susceptibility levels in recombinant lines of three crosses between susceptible and resistant lines, in a core collection of 32 lines, and in mutants affected in symbiosis with rhizobia. A GFP-expressing V. albo-atrum strain was used to study colonization of susceptible plants. Symptoms and colonization pattern in infected M. truncatula plants were typical of Verticillium wilt. Three distinct major quantitative trait loci were identified using a multicross, multisite design, suggesting that simple genetic mechanisms appear to control Verticillium wilt resistance in M. truncatula lines A17 and DZA45.5. The disease functional parameters varied largely in lines of the core collection. This biodiversity with regard to disease response encourages the development of association genetics and ecological approaches. Several mutants of the resistant line, impaired in different steps of rhizobial symbiosis, were affected in their response to V. albo-atrum, which suggests that mechanisms involved in the establishment of symbiosis or disease might have some common regulatory control points.

Genetic variability and QTL mapping of freezing tolerance and related traits in Medicago truncatula.

Freezing is a major environmental limitation to crop productivity for a number of species including legumes. We investigated the genetic determinism of freezing tolerance in the model legume Medicago truncatula Gaertn (M. truncatula). After having observed a large variation for freezing tolerance among 15 M. truncatula accessions, the progeny of a F6 recombinant inbred line population, derived from a cross between two accessions, was acclimated to low above-freezing temperatures and assessed for: (a) number of leaves (NOL), leaf area (LA), chlorophyll content index (CCI), shoot and root dry weights (SDW and RDW) at the end of the acclimation period and (b) visual freezing damage (FD) during the freezing treatment and 2 weeks after regrowth and foliar electrolyte leakage (EL) 2 weeks after regrowth. Consistent QTL positions with additive effects for FD were found on LG1, LG4 and LG6, the latter being the most explanatory (R (2) ≈ 40 %). QTL for NOL, QTL for EL, NOL and RDW, and QTL for EL and CCI colocalized with FD QTL on LG1, LG4 and LG6, respectively. Favorable alleles for these additive effects were brought by the same parent suggesting that this accession contributes to superior freezing tolerance by affecting plants' capacity to maintain growth at low above-freezing temperatures. No epistatic effects were found between FD QTL, but for each of the studied traits, 3-6 epistatic effects were detected between loci not detected directly as QTL. These results open the way to the assessment of syntenic relationships between QTL for frost tolerance in M. truncatula and cultivated legume species.

Five Regions of the Pea Genome Co-Control Partial Resistance to D. pinodes, Tolerance to Frost, and Some Architectural or Phenological Traits.

Evidence for reciprocal links between plant responses to biotic or abiotic stresses and architectural and developmental traits has been raised using approaches based on epidemiology, physiology, or genetics. Winter pea has been selected for years for many agronomic traits contributing to yield, taking into account architectural or phenological traits such as height or flowering date. It remains nevertheless particularly susceptible to biotic and abiotic stresses, among which Didymella pinodes and frost are leading examples. The purpose of this study was to identify and resize QTL localizations that control partial resistance to D. pinodes, tolerance to frost, and architectural or phenological traits on pea dense genetic maps, considering how QTL colocalizations may impact future winter pea breeding. QTL analysis revealed five metaQTLs distributed over three linkage groups contributing to both D. pinodes disease severity and frost tolerance. At these loci, the haplotypes of alleles increasing both partial resistance to D. pinodes and frost tolerance also delayed the flowering date, increased the number of branches, and/or decreased the stipule length. These results question both the underlying mechanisms of the joint control of biotic stress resistance, abiotic stress tolerance, and plant architecture and phenology and the methods of marker-assisted selection optimizing stress control and productivity in winter pea breeding.

New consistent QTL in pea associated with partial resistance to Aphanomyces euteiches in multiple French and American environments.

Partial resistances, often controlled by quantitative trait loci (QTL), are considered to be more durable than monogenic resistances. Therefore, a precursor to developing efficient breeding programs for polygenic resistance to pathogens should be a greater understanding of genetic diversity and stability of resistance QTL in plants. In this study, we deciphered the diversity and stability of resistance QTL to Aphanomyces euteiches in pea towards pathogen variability, environments and scoring criteria, from two new sources of partial resistance (PI 180693 and 552), effective in French and USA infested fields. Two mapping populations of 178 recombinant inbred lines each, derived from crosses between 552 or PI 180693 (partially resistant) and Baccara (susceptible), were used to identify QTL for Aphanomyces root rot resistance in controlled and in multiple French and USA field conditions using several resistance criteria. We identified a total of 135 additive-effect QTL corresponding to 23 genomic regions and 13 significant epistatic interactions associated with partial resistance to A. euteiches in pea. Among the 23 additive-effect genomic regions identified, five were consistently detected, and showed highly stable effects towards A. euteiches strains, environments, resistance criteria, condition tests and RIL populations studied. These results confirm the complexity of inheritance of partial resistance to A. euteiches in pea and provide good bases for the choice of consistent QTL to use in marker-assisted selection schemes to increase current levels of resistance to A. euteiches in pea breeding programs.

A complex genetic network involving a broad-spectrum locus and strain-specific loci controls resistance to different pathotypes of Aphanomyces euteiches in Medicago truncatula.

A higher understanding of genetic and genomic bases of partial resistance in plants and their diversity regarding pathogen variability is required for a more durable management of resistance genetic factors in sustainable cropping systems. In this study, we investigated the diversity of genetic factors involved in partial resistance to Aphanomyces euteiches, a very damaging pathogen on pea and alfalfa, in Medicago truncatula. A mapping population of 178 recombinant inbred lines, from the cross F83005.5 (susceptible) and DZA045.5 (resistant), was used to identify quantitative trait loci for resistance to four A. euteiches reference strains belonging to the four main pathotypes currently known on pea and alfalfa. A major broad-spectrum genomic region, previously named AER1, was localized to a reduced 440 kb interval on chromosome 3 and was involved in complete or partial resistance, depending on the A. euteiches strain. We also identified 21 additive and/or epistatic genomic regions specific to one or two strains, several of them being anchored to the M. truncatula physical map. These results show that, in M. truncatula, a complex network of genetic loci controls partial resistance to different pea and alfalfa pathotypes of A. euteiches, suggesting a diversity of molecular mechanisms underlying partial resistance.

Partial resistance of Medicago truncatula to Aphanomyces euteiches is associated with protection of the root stele and is controlled by a major QTL rich in proteasome-related genes.

A pathosystem between Aphanomyces euteiches, the causal agent of pea root rot disease, and the model legume Medicago truncatula was developed to gain insights into mechanisms involved in resistance to this oomycete. The F83005.5 French accession and the A17-Jemalong reference line, susceptible and partially resistant, respectively, to A. euteiches, were selected for further cytological and genetic analyses. Microscopy analyses of thin root sections revealed that a major difference between the two inoculated lines occurred in the root stele, which remained pathogen free in A17. Striking features were observed in A17 roots only, including i) frequent pericycle cell divisions, ii) lignin deposition around the pericycle, and iii) accumulation of soluble phenolic compounds. Genetic analysis of resistance was performed on an F7 population of 139 recombinant inbred lines and identified a major quantitative trait locus (QTL) near the top of chromosome 3. A second study, with near-isogenic line responses to A. euteiches confirmed the role of this QTL in expression of resistance. Fine-mapping allowed the identification of a 135-kb sequenced genomic DNA region rich in proteasome-related genes. Most of these genes were shown to be induced only in inoculated A17. Novel mechanisms possibly involved in the observed partial resistance are proposed.

Genetic and Pathogenicity Diversity of Aphanomyces euteiches Populations From Pea-Growing Regions in France.

Aphanomyces euteiches is an oomycete pathogen with a broad host-range on legumes that causes devastating root rot disease in many pea-growing countries and especially in France. Genetic resistance is a promising way to manage the disease since consistent QTL controlling partial resistance have been identified in near isogenic lines of pea. However, there are still no resistant pea varieties cultivated in France. This study aimed to evaluate the phenotypic and genetic diversity of A. euteiches populations from the major pea-growing regions in France. A collection of 205 isolates, from soil samples collected in infested pea fields located in five French regions, was established and genotyped using 20 SSR markers. Thirteen multilocus genotypes were found among the 205 isolates which displayed a low genotypic richness (ranged from 0 to 0.333). Two main clusters of isolates were identified using PCoA and STRUCTURE, including a predominant group comprising 88% of isolates and another group representing 12% of isolates mainly from the Bourgogne region. A subset of 34 isolates, representative of the fields sampled, was phenotyped for aggressiveness on a set of resistant and susceptible varieties of four legume hosts (pea, faba bean, vetch, alfalfa). Significant differences in disease severity were found among isolates and three groups of aggressiveness comprising 16, 17, and 2 isolates, respectively, were identified using HCA analysis. A higher diversity in pathogen aggressiveness was observed among isolates from Bourgogne, which included different legumes in its crop history. Little relationship was observed between genetic clusters and pathogenicity in the subset of 34 isolates, as expected using neutral markers. This study provides useful knowledge on the current state of low to moderate diversity among A. euteiches populations before resistant pea varieties are grown in France. New insights and hypotheses about the major factors shaping the diversity and evolution of A. euteiches are also discussed.

Comparative Genome-Wide-Association Mapping Identifies Common Loci Controlling Root System Architecture and Resistance to Aphanomyces euteiches in Pea.

Combining plant genetic resistance with architectural traits that are unfavorable to disease development is a promising strategy for reducing epidemics. However, few studies have identified root system architecture (RSA) traits with the potential to limit root disease development. Pea is a major cultivated legume worldwide and has a wide level of natural genetic variability for plant architecture. The root pathogen Aphanomyces euteiches is a major limiting factor of pea crop yield. This study aimed to increase the knowledge on the diversity of loci and candidate genes controlling RSA traits in pea and identify RSA genetic loci associated with resistance to A. euteiches which could be combined with resistance QTL in breeding. A comparative genome wide association (GWA) study of plant architecture and resistance to A. euteiches was conducted at the young plant stage in a collection of 266 pea lines contrasted for both traits. The collection was genotyped using 14,157 SNP markers from recent pea genomic resources. It was phenotyped for ten root, shoot and overall plant architecture traits, as well as three disease resistance traits in controlled conditions, using image analysis. We identified a total of 75 short-size genomic intervals significantly associated with plant architecture and overlapping with 46 previously detected QTL. The major consistent intervals included plant shoot architecture or flowering genes (PsLE, PsTFL1) with putative pleiotropic effects on root architecture. A total of 11 genomic intervals were significantly associated with resistance to A. euteiches confirming several consistent previously identified major QTL. One significant SNP, mapped to the major QTL Ae-Ps7.6, was associated with both resistance and RSA traits. At this marker, the resistance-enhancing allele was associated with an increased total root projected area, in accordance with the correlation observed between resistance and larger root systems in the collection. Seven additional intervals associated with plant architecture overlapped with GWA intervals previously identified for resistance to A. euteiches. This study provides innovative results about genetic interdependency of root disease resistance and RSA inheritance. It identifies pea lines, QTL, closely-linked markers and candidate genes for marker-assisted-selection of RSA loci to reduce Aphanomyces root rot severity in future pea varieties.

Combining gene expression and genetic analyses to identify candidate genes involved in cold responses in pea.

Cold stress affects plant growth and development. In order to better understand the responses to cold (chilling or freezing tolerance), we used two contrasted pea lines. Following a chilling period, the Champagne line becomes tolerant to frost whereas the Terese line remains sensitive. Four suppression subtractive hybridisation libraries were obtained using mRNAs isolated from pea genotypes Champagne and Terese. Using quantitative polymerase chain reaction (qPCR) performed on 159 genes, 43 and 54 genes were identified as differentially expressed at the initial time point and during the time course study, respectively. Molecular markers were developed from the differentially expressed genes and were genotyped on a population of 164 RILs derived from a cross between Champagne and Terese. We identified 5 candidate genes colocalizing with 3 different frost damage quantitative trait loci (QTL) intervals and a protein quantity locus (PQL) rich region previously reported. This investigation revealed the role of constitutive differences between both genotypes in the cold responses, in particular with genes related to glycine degradation pathway that could confer to Champagne a better frost tolerance. We showed that freezing tolerance involves a decrease of expression of genes related to photosynthesis and the expression of a gene involved in the production of cysteine and methionine that could act as cryoprotectant molecules. Although it remains to be confirmed, this study could also reveal the involvement of the jasmonate pathway in the cold responses, since we observed that two genes related to this pathway were mapped in a frost damage QTL interval and in a PQL rich region interval, respectively.

Mapping the genetic basis of symbiotic variation in legume-rhizobium interactions in Medicago truncatula.

Mutualisms are known to be genetically variable, where the genotypes differ in the fitness benefits they gain from the interaction. To date, little is known about the loci that underlie such genetic variation in fitness or whether the loci influencing fitness are partner specific, and depend on the genotype of the interaction partner. In the legume-rhizobium mutualism, one set of potential candidate genes that may influence the fitness benefits of the symbiosis are the plant genes involved in the initiation of the signaling pathway between the two partners. Here we performed quantitative trait loci (QTL) mapping in Medicago truncatula in two different rhizobium strain treatments to locate regions of the genome influencing plant traits, assess whether such regions are dependent on the genotype of the rhizobial mutualist (QTL × rhizobium strain), and evaluate the contribution of sequence variation at known symbiosis signaling genes. Two of the symbiotic signaling genes, NFP and DMI3, colocalized with two QTL affecting average fruit weight and leaf number, suggesting that natural variation in nodulation genes may potentially influence plant fitness. In both rhizobium strain treatments, there were QTL that influenced multiple traits, indicative of either tight linkage between loci or pleiotropy, including one QTL with opposing effects on growth and reproduction. There was no evidence for QTL × rhizobium strain or genotype × genotype interactions, suggesting either that such interactions are due to small-effect loci or that more genotype-genotype combinations need to be tested in future mapping studies.