A chickpea genetic variation map based on the sequencing of 3,366 genomes.

Zero hunger and good health could be realized by 2030 through effective conservation, characterization and utilization of germplasm resources1. So far, few chickpea (Cicer arietinum) germplasm accessions have been characterized at the genome sequence level2. Here we present a detailed map of variation in 3,171 cultivated and 195 wild accessions to provide publicly available resources for chickpea genomics research and breeding. We constructed a chickpea pan-genome to describe genomic diversity across cultivated chickpea and its wild progenitor accessions. A divergence tree using genes present in around 80% of individuals in one species allowed us to estimate the divergence of Cicer over the last 21 million years. Our analysis found chromosomal segments and genes that show signatures of selection during domestication, migration and improvement. The chromosomal locations of deleterious mutations responsible for limited genetic diversity and decreased fitness were identified in elite germplasm. We identified superior haplotypes for improvement-related traits in landraces that can be introgressed into elite breeding lines through haplotype-based breeding, and found targets for purging deleterious alleles through genomics-assisted breeding and/or gene editing. Finally, we propose three crop breeding strategies based on genomic prediction to enhance crop productivity for 16 traits while avoiding the erosion of genetic diversity through optimal contribution selection (OCS)-based pre-breeding. The predicted performance for 100-seed weight, an important yield-related trait, increased by up to 23% and 12% with OCS- and haplotype-based genomic approaches, respectively.

Proteomics analysis of Medicago truncatula response to infection by the phytopathogenic bacterium Ralstonia solanacearum points to jasmonate and salicylate defence pathways.

The infection of the model legume Medicago truncatula with Ralstonia solanacearum GMI1000 gives rise to bacterial wilt disease via colonisation of roots. The root and leaf responses to early infection (1 and 3 days post infection) were characterised to investigate the molecular mechanisms of plant resistance or susceptibility. A proteomics approach based on pools of susceptible and resistant recombinant inbred lines was used to specifically target the mechanisms for tolerance. Differential abundances were evidenced for proteins involved in defence (e.g., PR5, PR10, or Kunitz protease inhibitors) and signalling pathways (such as cyclophilin). R. solanacearum inoculation modifies expression levels of those genes, either in both genotypes (AOS1, LOX4, and proteinase inhibitors) or specifically in the resistant line (PR proteins). Exogenous application of salicylic acid (SA) enhanced tolerance to the bacteria, whereas methyl jasmonate (MeJA) enhanced short-term tolerance then promoted disease in the susceptible ecotype, suggesting that they may mediate defence responses. Conversely, proteomics-identified genes were also shown to be SA or MeJA responsive. This is the first description of differential response to R. solanacearum in M. truncatula. Our results suggest that root basal defence is activated at 1 dpi, together with the JA pathway. Specific resistance is then evidenced at three dpi, with the up-regulation of SA-dependent PR proteins.

Medicago truncatula Oleanolic-Derived Saponins Are Correlated with Caterpillar Deterrence.

Plant resistance mechanisms to insect herbivory can potentially be bred into crops as an important strategy for integrated pest management. Medicago truncatula ecotypes inoculated with the rhizobium Ensifer medicae (Sinorhizobium medica) WSM419 were screened for resistance to herbivory by caterpillars of the beet armyworm, Spodoptera exigua, through leaf and whole plant choice studies; TN1.11 and F83005.5 are identified as the least and most deterrent ecotypes, respectively. In response to caterpillar herbivory, both ecotypes mount a robust burst of plant defensive jasmonate phytohormones. Restriction of caterpillars to either of these ecotypes does not adversely affect pest performance. This argues for an antixenosis (deterrence) resistance mechanism associated with the F83005.5 ecotype. Unbiased metabolomic profiling identified strong ecotype-specific differences in metabolite profile, particularly in the content of oleanolic-derived saponins that may act as antifeedants. Compared to the more susceptible ecotype, F83005.5 has higher levels of oleanolic-type zanhic acid- and medicagenic acid-derived compounds. Together, these data support saponin-mediated deterrence as a resistance mechanism of the F83005.5 ecotype and implicates these compounds as potential antifeedants that could be used in agricultural sustainable pest management strategies.

Cell length instead of cell number becomes the predominant factor contributing to hypocotyl length genotypic differences under abiotic stress in Medicago truncatula.

Hypocotyl elongation in the dark is a crucial process to ensure seedling emergence. It relies both on the cell number and cell length. The contribution of these two factors to the maximal hypocotyl length and the impact of environmental conditions on this contribution are not known. This is surprising considering the agronomic and economical importance of seedling emergence in crop establishment. Using 14 genotypes from a nested core collection representing Medicago truncatula (barrel medic) natural variation, we investigated how epidermal cell number and cell length contribute to hypocotyl length under optimal, low temperature (8°C) and water deficit (-0.50 MPa) conditions. Both cell number and length vary according to genotypes and contribute to maximal hypocotyl length differences between genotypes. This contribution, however, depends on growth conditions. Cell number is the major contributor under optimal conditions (60%) whereas cell length becomes the major determinant under stress. Maximal hypocotyl length is correlated with hypocotyl elongation rate under both stresses but not under optimal condition, revealing contrasted genotypes for cell elongation capacity under stress. To identify the genetic regulators determining cell number and cell length, quantitative trait loci (QTLs) were detected using a recombinant inbred lines population exhibiting segregation in maximal hypocotyl length. Two QTLs controlling cell number and three QTLs controlling cell length at low temperature were detected. One QTL for cell number and two for cell length were found to be associated with hypocotyl length under low temperature. This study provides new information to improve seedling emergence under abiotic stress.

Whole-genome nucleotide diversity, recombination, and linkage disequilibrium in the model legume Medicago truncatula.

Medicago truncatula is a model for investigating legume genetics, including the genetics and evolution of legume-rhizobia symbiosis. We used whole-genome sequence data to identify and characterize sequence polymorphisms and linkage disequilibrium (LD) in a diverse collection of 26 M. truncatula accessions. Our analyses reveal that M. truncatula harbors both higher diversity and less LD than soybean (Glycine max) and exhibits patterns of LD and recombination similar to Arabidopsis thaliana. The population-scaled recombination rate is approximately one-third of the mutation rate, consistent with expectations for a species with a high selfing rate. Linkage disequilibrium, however, is not extensive, and therefore, the low recombination rate is likely not a major constraint to adaptation. Nucleotide diversity in 100-kb windows was negatively correlated with gene density, which is expected if diversity is shaped by selection acting against slightly deleterious mutations. Among putative coding regions, members of four gene families harbor significantly higher diversity than the genome-wide average. Three of these families are involved in resistance against pathogens; one of these families, the nodule-specific, cysteine-rich gene family, is specific to the galegoid legumes and is involved in control of rhizobial differentiation. The more than 3 million SNPs that we detected, approximately one-half of which are present in more than one accession, are a valuable resource for genome-wide association mapping of genes responsible for phenotypic diversity in legumes, especially traits associated with symbiosis and nodulation.

Molecular Markers Help with Breeding for Agronomic Traits of Spring Wheat in Kazakhstan and Siberia.

The Kazakhstan-Siberia Network for Spring Wheat Improvement (KASIB) was established in 2000, forming a collaboration between breeding and research programs through biannual yield trials. A core set of 142 genotypes from 15 breeding programs was selected, genotyped for 81 DNA functional markers and phenotyped for 10 agronomic traits at three sites in Kazakhstan (Karabalyk, Shortandy and Shagalaly) and one site in Russia (Omsk) in 2020-2022. The study aim was to identify markers demonstrating significant effects on agronomic traits. The average grain yield of individual trials varied from 118 to 569 g/m2. Grain yield was positively associated with the number of days to heading, plant height, number of grains per spike and 1000-kernel weight. Eight DNA markers demonstrated significant effects. The spring-type allele of the Vrn-A1 gene accelerated heading by two days (5.6%) and was present in 80% of the germplasm. The winter allele of the Vrn-A1 gene significantly increased grain yield by 2.7%. The late allele of the earliness marker per se, TaMOT1-D1, delayed development by 1.9% and increased yield by 4.5%. Translocation of 1B.1R was present in 21.8% of the material, which resulted in a 6.2% yield advantage compared to 1B.1B germplasm and a reduction in stem rust severity from 27.6 to 6.6%. The favorable allele of TaGS-D1 increased both kernel weight and yield by 2-3%. Four markers identified in ICARDA germplasm, ISBW2-GY (Kukri_c3243_1065, 3B), ISBW3-BM (TA004946-0577, 1B), ISBW10-SM2 (BS00076246_51, 5A), ISBW11-GY (wsnp_Ex_c12812_20324622, 4A), showed an improved yield in this study of 3-4%. The study recommends simultaneous validation and use of selected markers in KASIB's network.

MtQRRS1, an R-locus required for Medicago truncatula quantitative resistance to Ralstonia solanacearum.

Ralstonia solanacearum is a major soilborne pathogen that attacks > 200 plant species, including major crops. To characterize MtQRRS1, a major quantitative trait locus (QTL) for resistance towards this bacterium in the model legume Medicago truncatula, genetic and functional approaches were combined. QTL analyses together with disease scoring of heterogeneous inbred families were used to define the locus. The candidate region was studied by physical mapping using a bacterial artificial chromosome (BAC) library of the resistant line, and sequencing. In planta bacterial growth measurements, grafting experiments and gene expression analysis were performed to investigate the mechanisms by which this locus confers resistance to R. solanacearum. The MtQRRS1 locus was localized to the same position in two recombinant inbred line populations and was narrowed down to a 64 kb region. Comparison of parental line sequences revealed 15 candidate genes with sequence polymorphisms, but no evidence of differential gene expression upon infection. A role for the hypocotyl in resistance establishment was shown. These data indicate that the quantitative resistance to bacterial wilt conferred by MtQRRS1, which contains a cluster of seven R genes, is shared by different accessions and may act through intralocus interactions to promote resistance.

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.

Temperature increase modifies susceptibility to Verticillium wilt in Medicago spp and may contribute to the emergence of more aggressive pathogenic strains.

Global warming is expected to have a direct impact on plant disease patterns in agro-eco-systems. However, few analyses report the effect of moderate temperature increase on disease severity due to soil-borne pathogens. For legumes, modifications of root plant-microbe interactions either mutualistic or pathogenic due to climate change may have dramatic effects. We investigated the effect of increasing temperature on the quantitative disease resistance to Verticillium spp., a major soil-borne fungal pathogen, in the model legume Medicago truncatula and the crop M. sativa. First, twelve pathogenic strains isolated from various geographical origin were characterized with regard to their in vitro growth and pathogenicity at 20°C, 25°C and 28°C. Most of them exhibited 25°C as the optimum temperature for in vitro parameters, and between 20°C and 25°C for pathogenicity. Second, a V. alfalfae strain was adapted to the higher temperature by experimental evolution, i.e. three rounds of UV mutagenesis and selection for pathogenicity at 28°C on a susceptible M. truncatula genotype. Inoculation of monospore isolates of these mutants on resistant and susceptible M. truncatula accessions revealed that at 28°C they were all more aggressive than the wild type strain, and that some had acquired the ability to cause disease on resistant genotype. Third, one mutant strain was selected for further studies of the effect of temperature increase on the response of M. truncatula and M. sativa (cultivated alfalfa). The response of seven contrasted M. truncatula genotypes and three alfalfa varieties to root inoculation was followed using disease severity and plant colonization, at 20°C, 25°C and 28°C. With increasing temperature, some lines switched from resistant (no symptoms, no fungus in the tissues) to tolerant (no symptoms but fungal growth into the tissues) phenotypes, or from partially resistant to susceptible. Further studies in greenhouse evidence the reduction in plant fitness due to disease in susceptible lines. We thus report that root pathogenic interactions are affected by anticipated global warming, with trends towards increased plant susceptibility and larger virulence for hot-adapted strains. New threats due to hot-adapted strains of soil-borne pathogens, with possibly wider host range and increased aggressiveness, might occur.

Medicago truncatula quantitative resistance to a new strain of Verticillium alfalfae from Iran revealed by a genome-wide association study.

Verticillium wilt is a major threat to many crops, among them alfalfa (Medicago sativa). The model plant Medicago truncatula, a close relative of alfalfa was used to study the genetic control of resistance towards a new Verticillium alfalfae isolate. The accidental introduction of pathogen strains through global trade is a threat to crop production and such new strains might also be better adapted to global warming. Isolates of V. alfalfae were obtained from alfalfa fields in Iran and characterized. The Iranian isolate AF1 was used in a genome-wide association study (GWAS) involving 242 accessions from the Mediterranean region. Root inoculations were performed with conidia at 25°C and symptoms were scored regularly. Maximum Symptom Score and Area under Disease Progess Curve were computed as phenotypic traits to be used in GWAS and for comparison to a previous study with French isolate V31.2 at 20°C. This comparison showed high correlation with a shift to higher susceptibility, and similar geographical distribution of resistant and susceptible accessions to AF1 at 25°C, with resistant accessions mainly in the western part. GWAS revealed 30 significant SNPs linked to resistance towards isolate AF1. None of them were common to the previous study with isolate V31.2 at 20°C. To confirm these loci, the expression of nine underlying genes was studied. All genes were induced in roots following inoculation, in susceptible and resistant plants. However, in resistant plants induction was higher and lasted longer. Taken together, the use of a new pathogen strain and a shift in temperature revealed a completely different genetic control compared to a previous study that demonstrated the existence of two major QTLs. These results can be useful for Medicago breeding programs to obtain varieties better adapted to future conditions.

QTL mapping of oleic acid content in modern VNIIMK sunflower (Helianthus annuus L.) lines by using GBS-based SNP map.

Oleic acid is a monounsaturated fatty acid increasing oil oxidative stability. High content of oleic acid is thus a valuable trait in oilseed crops. Sunflower (Helianthus annuus L.) normally accumulates linoleic acid as a major fatty acid, but a mutant expressing a high oleic phenotype form was previously obtained by chemical mutagenesis and mapped on the sunflower genome. Several studies suggest the presence of additional genes involved in the control of the high content of oleic acid, with their expression possibly depending on the genetic background. To test this hypothesis, we performed a QTL mapping of the high oleic acid trait within two independent F2 crosses involving lines with contrasting oleic acid content from the Pustovoit All-Russia Research Institute of Oil Crops (VNIIMK) collection. We applied genotyping-by-sequencing (GBS) to construct single nucleotide polymorphism-based genetic maps and performed QTL mapping using quantitative and qualitative encoding for oleic acid content. Our results support the finding that the oleic acid content in the assessed crosses is controlled by one major effect locus. However, different dominant/recessive effects of the major locus were reported for both crosses. Additionally, a possible translocation between chromosome 7 and 14 was reported in one assessed cross. We defined a set of single nucleotide polymorphism markers for each cross which could be used for marker-assisted selection.

Comparative proteomic analysis of the thermotolerant plant Portulaca oleracea acclimation to combined high temperature and humidity stress.

Elevated temperature and humidity are major environmental factors limiting crop yield and distribution. An understanding of the mechanisms underlying plant tolerance to high temperature and humidity may facilitate the development of cultivars adaptable to warm or humid regions. Under conditions of 90% humidity and 35 °C, the thermotolerant plant Portulaca oleracea exhibits excellent photosynthetic capability and relatively little oxidative damage. To determine the proteomic response that occurs in leaves of P. oleracea following exposure to high temperature and high humidity, a proteomic approach was performed to identify protein changes. A total of 51 differentially expressed proteins were detected and characterized functionally and structurally; these identified proteins were involved in various functional categories, mainly including material and energy metabolism, the antioxidant defense responses, protein destination and storage, and transcriptional regulation. The subset of antioxidant defense-related proteins demonstrated marked increases in activity with exposure to heat and humidity, which led to lower accumulations of H(2)O(2) and O(2)(-) in P. oleracea compared with the thermosensitive plant Arabidopsis thaliana. The quickly accumulations of proline content and heat-shock proteins, and depleting abscisic acid (ABA) via increasing ABA-8'-hydroxylase were also found in P. oleracea under stress conditions, that resulted into greater stomata conductance and respiration rates. On the basis of these findings, we propose that P. oleracea employs multiple strategies to enhance its adaptation to high-temperature and high-humidity conditions.

Using a physiological framework for improving the detection of quantitative trait loci related to nitrogen nutrition in Medicago truncatula.

Medicago truncatula is used as a model plant for exploring the genetic and molecular determinants of nitrogen (N) nutrition in legumes. In this study, our aim was to detect quantitative trait loci (QTL) controlling plant N nutrition using a simple framework of carbon/N plant functioning stemming from crop physiology. This framework was based on efficiency variables which delineated the plant's efficiency to take up and process carbon and N resources. A recombinant inbred line population (LR4) was grown in a glasshouse experiment under two contrasting nitrate concentrations. At low nitrate, symbiotic N(2) fixation was the main N source for plant growth and a QTL with a large effect located on linkage group (LG) 8 affected all the traits. Significantly, efficiency variables were necessary both to precisely localize a second QTL on LG5 and to detect a third QTL involved in epistatic interactions on LG2. At high nitrate, nitrate assimilation was the main N source and a larger number of QTL with weaker effects were identified compared to low nitrate. Only two QTL were common to both nitrate treatments: a QTL of belowground biomass located at the bottom of LG3 and another one on LG6 related to three different variables (leaf area, specific N uptake and aboveground:belowground biomass ratio). Possible functions of several candidate genes underlying QTL of efficiency variables could be proposed. Altogether, our results provided new insights into the genetic control of N nutrition in M. truncatula. For instance, a novel result for M. truncatula was identification of two epistatic interactions in controlling plant N(2) fixation. As such this study showed the value of a simple conceptual framework based on efficiency variables for studying genetic determinants of complex traits and particularly epistatic interactions.

Genetic mapping of loci involved in oil tocopherol composition control in Russian sunflower (Helianthus annuus L.) lines.

Tocopherols are antioxidants that preserve oil lipids against oxidation and serve as a natural source of vitamin E in the human diet. Compared with other major oilseeds like rapeseed and soybean, sunflower (Helianthus annuus L.) exhibits low phenotypic diversity of tocopherol composition, both in wild and cultivated accessions from germplasm collections. Two major mutations that alter tocopherol composition were identified in genetic collections, and several studies suggested additional loci controlling tocopherol composition, with their expression possibly depending on the genetic background. In the present study, we performed QTL mapping of tocopherol composition in two independent F2 crosses between lines with contrasting tocopherol composition from the Pustovoit All-Russia Research Institute of Oil Crops (VNIIMK) collection. We used genotyping-bysequencing (GBS) to construct single nucleotide polymorphism-based genetic maps, and performed QTL mapping using quantitative and qualitative encoding for phenotypic traits. Our results support the notion that the tocopherol composition in the assessed crosses is controlled by two loci. We additionally selected and validated two single nucleotide polymorphism markers for each cross which could be used for marker-assisted selection.

Identification of Key Metabolic Pathways and Biomarkers Underlying Flowering Time of Guar (Cyamopsis tetragonoloba (L.) Taub.) via Integrated Transcriptome-Metabolome Analysis.

Guar (Cyamopsis tetragonoloba (L.) Taub.) is an annual legume crop native to India and Pakistan. Seeds of the plant serve as a source of galactomannan polysaccharide (guar gum) used in the food industry as a stabilizer (E412) and as a gelling agent in oil and gas fracturing fluids. There were several attempts to introduce this crop to countries of more northern latitudes. However, guar is a plant of a short photoperiod, therefore, its introduction, for example, to Russia is complicated by a long day length during the growing season. Breeding of new guar varieties insensitive to photoperiod slowed down due to the lack of information on functional molecular markers, which, in turn, requires information on guar genome. Modern breeding strategies, e.g., genomic predictions, benefit from integration of multi-omics approaches such as transcriptome, proteome and metabolome assays. Here we present an attempt to use transcriptome-metabolome integration to understand the genetic determination of flowering time variation among guar plants that differ in their photoperiod sensitivity. This study was performed on nine early- and six delayed-flowering guar varieties with the goal to find a connection between 63 metabolites and 1,067 differentially expressed transcripts using Shiny GAM approach. For the key biomarker of flowering in guar myo-inositol we also evaluated the KEGG biochemical pathway maps available for Arabidopsis thaliana. We found that the phosphatidylinositol signaling pathway is initiated in guar plants that are ready for flowering through the activation of the phospholipase C (PLC) gene, resulting in an exponential increase in the amount of myo-inositol in its free form observed on GC-MS chromatograms. The signaling pathway is performed by suppression of myo-inositol phosphate kinases (phosphorylation) and alternative overexpression of phosphatases (dephosphorylation). Our study suggests that metabolome and transcriptome information taken together, provide valuable information about biomarkers that can be used as a tool for marker-assisted breeding, metabolomics and functional genomics of this important legume crop.

Development of SNP Set for the Marker-Assisted Selection of Guar (Cyamopsis tetragonoloba (L.) Taub.) Based on a Custom Reference Genome Assembly.

Guar gum, a polysaccharide derived from guar seeds, is widely used in a variety of industrial applications, including oil and gas production. Although guar is mostly propagated in India, interest in guar as a new industrial legume crop is increasing worldwide, demanding the development of effective tools for marker-assisted selection. In this paper, we report a wide-ranging set of 4907 common SNPs and 327 InDels generated from RADseq genotyping data of 166 guar plants of different geographical origin. A custom guar reference genome was assembled and used for variant calling. A consensus set of variants was built using three bioinformatic pipelines for short variant discovery. The developed molecular markers were used for genome-wide association study, resulting in the discovery of six markers linked to the variation of an important agronomic trait-percentage of pods matured to the harvest date under long light day conditions. One of the associated variants was found inside the putative transcript sequence homologous to an ABC transporter in Arabidopsis, which has been shown to play an important role in D-myo-inositol phosphates metabolism. Earlier, we suggested that genes involved in myo-inositol phosphate metabolism have significant impact on the early flowering of guar plants. Hence, we believe that the developed SNP set allows for the identification of confident molecular markers of important agrobiological traits.

The RAP1 gene confers effective, race-specific resistance to the pea aphid in Medicago truncatula independent of the hypersensitive reaction.

Plant resistance to pathogens is commonly associated with a hypersensitive response (HR), but the degree to which the HR is responsible for incompatibility is subject to debate. Resistance to aphids is likely to share features with resistance to pathogens but is less well understood. Here, we report effective resistance to the pea aphid Acyrthosiphon pisum in Medicago truncatula. Aphids lost weight and died rapidly (within two days) on the resistant genotype Jemalong, which developed necrotic lesions following infestation. Lesions were induced by nonvascular intracellular stylet punctures by aphids, remained localized to the site of stylet entry, stained for the presence of reactive oxygen species, and were similar to the HR induced by the bacterial pathogen Pseudomonas syringae pv. phaseolicola. The implication that aphid-induced lesions confer resistance was tested by quantitative trait loci analysis using recombinant inbred lines derived from a cross between Jemalong and the susceptible genotype DZA315.16. One major locus, RAP1, was identified that was sufficient to confer race-specific resistance against the pea aphid and was mapped to the middle of chromosome 3. Surprisingly, a separate locus, mapping to the top of chromosome 3, governed aphid-induced HR, indicating that the HR-like lesions are not required for RAP1-mediated aphid resistance.

WhoGEM: an admixture-based prediction machine accurately predicts quantitative functional traits in plants.

The explosive growth of genomic data provides an opportunity to make increased use of sequence variations for phenotype prediction. We have developed a prediction machine for quantitative phenotypes (WhoGEM) that overcomes some of the bottlenecks limiting the current methods. We demonstrated its performance by predicting quantitative disease resistance and quantitative functional traits in the wild model plant species, Medicago truncatula, using geographical locations as covariates for admixture analysis. The method's prediction reliability equals or outperforms all existing algorithms for quantitative phenotype prediction. WhoGEM analysis produces evidence that variation in genome admixture proportions explains most of the phenotypic variation for quantitative phenotypes.

Distinct Lotus japonicus Transcriptomic Responses to a Spectrum of Bacteria Ranging From Symbiotic to Pathogenic.

Lotus japonicus is a well-studied nodulating legume and a model organism for the investigation of plant-microbe interactions. The majority of legume transcriptome studies have focused on interactions with compatible symbionts, whereas responses to non-adapted rhizobia and pathogenic bacteria have not been well-characterized. In this study, we first characterized the transcriptomic response of L. japonicus to its compatible symbiont, Mesorhizobium loti R7A, through RNA-seq analysis of various plant tissues. Early symbiotic signaling was largely Nod factor-dependent and enhanced within root hairs, and we observed large-scale transcriptional reprogramming in nodule primordia and mature nitrogen-fixing nodules. We then characterized root transcriptional responses to a spectrum of L. japonicus interacting bacteria ranging from semi-compatible symbionts to pathogens. M. loti R7A and the semi-compatible strain Sinorhizobium fredii HH103 showed remarkably similar responses, allowing us to identify a small number of genes potentially involved in differentiating between fully and semi-compatible symbionts. The incompatible symbiont Bradyrhizobium elkanii USDA61 induced a more attenuated response, but the weakest response was observed for the foliar pathogen Pseudomonas syringae pv. tomato DC3000, where the affected genes also responded to other tested bacteria, pointing to a small set of common bacterial response genes. In contrast, the root pathogen Ralstonia solanacearum JS763 induced a pronounced and distinct transcriptomic pathogen response, which we compared to the results of the other treatments. This comparative analysis did not support the concept that an early defense-like response is generally evoked by compatible rhizobia during establishment of symbiosis.