Medicago truncatula: Genetic and Genomic Resources

Genome-wide methylation landscape during somatic embryogenesis in Medicago truncatula reveals correlation between Tnt1 retrotransposition and hyperactive methylation regions

Medicago truncatula is a model legume for fundamental research on legume biology and symbiotic nitrogen fixation. Tnt1, a retrotransposon from tobacco, was used to generate insertion mutants in M. truncatula R108. Approximately 21 000 insertion lines have been generated and publicly available. Tnt1 retro-transposition event occurs during somatic embryogenesis (SE), a pivotal process that triggers massive methylation changes. We studied the SE of M. truncatula R108 using leaf explants and explored the dynamic shifts in the methylation landscape from leaf explants to callus formation and finally embryogenesis. Higher cytosine methylation in all three contexts of CG, CHG, and CHH patterns was observed during SE compared to the controls. Higher methylation patterns were observed in assumed promoter regions (~2-kb upstream regions of transcription start site) of the genes, while lowest was recorded in the untranslated regions. Differentially methylated promoter region analysis showed a higher CHH methylation in embryogenesis tissue samples when compared to CG and CHG methylation. Strong correlation (89.71%) was identified between the differentially methylated regions (DMRs) and the site of Tnt1 insertions in M. truncatula R108 and stronger hypermethylation of genes correlated with higher number of Tnt1 insertions in all contexts of CG, CHG, and CHH methylation. Gene ontology enrichment and KEGG pathway enrichment analysis identified genes and pathways enriched in the signal peptide processing, ATP hydrolysis, RNA polymerase activity, transport, secondary metabolites, and nitrogen metabolism pathways. Combined gene expression analysis and methylation profiling showed an inverse relationship between methylation in the DMRs (regions spanning genes) and the expression of genes. Our results show that a dynamic shift in methylation happens during the SE process in the context of CG, CHH and CHG methylation, and the Tnt1 retrotransposition correlates with the hyperactive methylation regions.

Mutation of BAM2 rescues the sunn hypernodulation phenotype in Medicago truncatula, suggesting that a signaling pathway like CLV1/BAM in Arabidopsis affects nodule number

The unique evolutionary adaptation of legumes for nitrogen-fixing symbiosis leading to nodulation is tightly regulated by the host plant. The autoregulation of nodulation (AON) pathway negatively regulates the number of nodules formed in response to the carbon/nitrogen metabolic status of the shoot and root by long-distance signaling to and from the shoot and root. Central to AON signaling in the shoots of Medicago truncatula is SUNN, a leucine-rich repeat receptor-like kinase with high sequence similarity with CLAVATA1 (CLV1), part of a class of receptors in Arabidopsis involved in regulating stem cell populations in the root and shoot. This class of receptors in Arabidopsis includes the BARELY ANY MERISTEM family, which, like CLV1, binds to CLE peptides and interacts with CLV1 to regulate meristem development. M. truncatula contains five members of the BAM family, but only MtBAM1 and MtBAM2 are highly expressed in the nodules 48 hours after inoculation. Plants carry mutations in individual MtBAMs, and several double BAM mutant combinations all displayed wild-type nodule number phenotypes. However, Mtbam2 suppressed the sunn-5 hypernodulation phenotype and partially rescued the short root length phenotype of sunn-5 when present in a sunn-5 background. Grafting determined that bam2 suppresses supernodulation from the roots, regardless of the SUNN status of the root. Overexpression of MtBAM2 in wild-type plants increases nodule numbers, while overexpression of MtBAM2 in some sunn mutants rescues the hypernodulation phenotype, but not the hypernodulation phenotypes of AON mutant rdn1-2 or crn. Relative expression measurements of the nodule transcription factor MtWOX5 downstream of the putative bam2 sunn-5 complex revealed disruption of meristem signaling; while both bam2 and bam2 sunn-5 influence MtWOX5 expression, the expression changes are in different directions. We propose a genetic model wherein the specific root interactions of BAM2/SUNN are critical for signaling in nodule meristem cell homeostasis in M. truncatula.

Analysis of the structure and function of the LYK cluster of Medicago truncatula A17 and R108

The establishment of the Legume-Rhizobia symbiosis is generally dependent on the production of rhizobial lipochitooligosaccharidic Nod factors (NFs) and their perception by plant Lysin Motif Receptor-Like Kinases (LysM-RLKs). In this study, we characterized a cluster of LysM-RLK genes implicated in strain-specific recognition in two highly divergent and widely-studied Medicago truncatula genotypes, A17 and R108. We then used reverse genetic approaches and biochemical analyses to study the function of selected genes in the clusters and the ability of their encoded proteins to bind NFs. Our study has revealed that the LYK cluster exhibits a high degree of variability among M. truncatula genotypes, which in A17 and R108 includes recent recombination events within the cluster and a transposon insertion in A17. The essential role of LYK3 in nodulation in A17 is not conserved in R108 despite similar sequences and good nodulation expression profiles. Although, LYK2, LYK5 and LYK5bis are not essential for nodulation of the two genotypes, some evidence points to accessory roles in nodulation, but not through high-affinity NF binding. This work shows that recent evolution in the LYK cluster provides a source of variation for nodulation, and potential robustness of signaling through genetic redundancy.

Transcript profiling of hairy vetch (Vicia villosa Roth) identified interesting genes for seed dormancy

Hairy vetch, a diploid annual legume species, has a robust growth habit, high biomass yield, and winter hardy characteristics. Seed hardness is a major constraint for growing hairy vetch commercially. Hard seeded cultivars are valuable as forages, whereas soft seeded and shatter resistant cultivars have advantages for their use as a cover crop. Transcript analysis of hairy vetch was performed to understand the genetic mechanisms associated with important hairy vetch traits. RNA was extracted from leaves, flowers, immature pods, seed coats, and cotyledons of contrasting soft and hard seeded "AU Merit" plants. A range of 31.22-79.18 Gb RNA sequence data per tissue sample were generated with estimated coverage of 1040-2639×. RNA sequence assembly and mapping of the contigs against the Medicago truncatula (V4.0) genome identified 76,422 gene transcripts. A total of 24,254 transcripts were constitutively expressed in hairy vetch tissues. Key genes, such as KNOX4 (a class II KNOTTED-like homeobox KNOXII gene), qHs1 (endo-1,4-β-glucanase), GmHs1-1 (calcineurin-like metallophosphoesterase), chitinase, shatterproof 1 and 2 (SHP1, SHP2), shatter resistant 1-5 (SHAT1-5)(NAC transcription factor), PDH1 (prephenate dehydrogenase 1), and pectin methylesterases with a potential role in seed hardness and pod shattering, were further explored based on genes involved in seed hardness from other species to query the hairy vetch transcriptome data. Identification of interesting candidate genes in hairy vetch can facilitate the development of improved cultivars with desirable seed characteristics for use as a forage and as a cover crop.

Medicago truncatula resources to study legume biology and symbiotic nitrogen fixation

Medicago truncatula is a chosen model for legumes towards deciphering fundamental legume biology, especially symbiotic nitrogen fixation. Current genomic resources for M. truncatula include a completed whole genome sequence information for R108 and Jemalong A17 accessions along with the sparse draft genome sequences for other 226 M. truncatula accessions. These genomic resources are complemented by the availability of mutant resources such as retrotransposon (Tnt1) insertion mutants in R108 and fast neutron bombardment (FNB) mutants in A17. In addition, several M. truncatula databases such as small secreted peptides (SSPs) database, transporter protein database, gene expression atlas, proteomic atlas, and metabolite atlas are available to the research community. This review describes these resources and provide information regarding how to access these resources.

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.

Plastid phylogenomics uncovers multiple species in Medicago truncatula (Fabaceae) germplasm accessions

Medicago truncatula is a model legume that has been extensively investigated in diverse subdisciplines of plant science. Medicago littoralis can interbreed with M. truncatula and M. italica; these three closely related species form a clade, i.e. TLI clade. Genetic studies have indicated that M. truncatula accessions are heterogeneous but their taxonomic identities have not been verified. To elucidate the phylogenetic position of diverse M. truncatula accessions within the genus, we assembled 54 plastid genomes (plastomes) using publicly available next-generation sequencing data and conducted phylogenetic analyses using maximum likelihood. Five accessions showed high levels of plastid DNA polymorphism. Three of these highly polymorphic accessions contained sequences from both M. truncatula and M. littoralis. Phylogenetic analyses of sequences placed some accessions closer to distantly related species suggesting misidentification of source material. Most accessions were placed within the TLI clade and maximally supported the interrelationships of three subclades. Two Medicago accessions were placed within a M. italica subclade of the TLI clade. Plastomes with a 45-kb (rpl20-ycf1) inversion were placed within the M. littoralis subclade. Our results suggest that the M. truncatula accession genome pool represents more than one species due to possible mistaken identities and gene flow among closely related species.

A newly evolved chimeric lysin motif receptor-like kinase in Medicago truncatula spp. tricycla R108 extends its Rhizobia symbiotic partnership

Rhizobial lipochitooligosaccharidic Nod factors (NFs), specified by nod genes, are the primary determinants of host specificity in the legume-Rhizobia symbiosis. We examined the nodulation ability of Medicago truncatula cv Jemalong A17 and M. truncatula ssp. tricycla R108 with the Sinorhizobium meliloti nodF/nodL mutant, which produces modified NFs. We then applied genetic and functional approaches to study the genetic basis and mechanism of nodulation of R108 by this mutant. We show that the nodF/nodL mutant can nodulate R108 but not A17. Using genomics and reverse genetics, we identified a newly evolved, chimeric LysM receptor-like kinase gene in R108, LYK2bis, which is responsible for the phenotype and can allow A17 to gain nodulation with the nodF/nodL mutant. We found that LYK2bis is involved in nodulation by mutants producing nonO-acetylated NFs and interacts with the key receptor protein NFP. Many, but not all, natural S. meliloti and S. medicae strains tested require LYK2bis for efficient nodulation of R108. Our findings reveal that a newly evolved gene in R108, LYK2bis, extends nodulation specificity to mutants producing nonO-acetylated NFs and is important for nodulation by many natural Sinorhizobia. Evolution of this gene may present an adaptive advantage to allow nodulation by a greater variety of strains.

Correction to: Comparison of efficiency and time to regeneration of Agrobacterium-mediated transformation methods in Medicago truncatula

[This corrects the article DOI: 10.1186/s13007-019-0404-1.].

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.

Somatic Embryogenesis in the Medicago truncatula Model: Cellular and Molecular Mechanisms

Medicago truncatula is now widely regarded as a legume model where there is an increasing range of genomic resources. Highly regenerable lines have been developed from the wild-type Jemalong cultivar, most likely due to epigenetic changes. These lines with high rates of somatic embryogenesis (SE) can be compared with wild-type where SE is rare. Much of the research has been with the high SE genotype Jemalong 2HA (2HA). SE can be induced from leaf tissue explants or isolated mesophyll protoplasts. In 2HA, the exogenous phytohormones 1-naphthaleneacetic acid (NAA) and 6-benzylaminopurine (BAP) are central to SE. However, there are interactions with ethylene, abscisic acid (ABA), and gibberellic acid (GA) which produce maximum SE. In the main, somatic embryos are derived from dedifferentiated cells, undergo organellar changes, and produce stem-like cells. There is evidence that the SE is induced as a result of a stress and hormone interaction and this is discussed. In M. truncatula, there are connections between stress and specific up-regulated genes and specific hormones and up-regulated genes during the SE induction phase. Some of the transcription factors have been knocked down using RNAi to show they are critical for SE induction (MtWUSCHEL, MtSERF1). SE research in M. truncatula has utilized high throughput transcriptomic and proteomic studies and the more detailed investigation of some individual genes. In this review, these studies are integrated to suggest a framework and timeline for some of the key events of SE induction in M. truncatula.