Delineating the Tnt1 Insertion Landscape of the Model Legume Medicago truncatula cv. R108 at the Hi-C Resolution Using a Chromosome-Length Genome Assembly.

Legumes are of great interest for sustainable agricultural production as they fix atmospheric nitrogen to improve the soil. is a well-established model legume, and extensive studies in fundamental molecular, physiological, and developmental biology have been undertaken to translate into trait improvements in economically important legume crops worldwide. However, reference genome was generated in the accession Jemalong A17, which is highly recalcitrant to transformation. R108 is more attractive for genetic studies due to its high transformation efficiency and -insertion population resource for functional genomics. The need to perform accurate synteny analysis and comprehensive genome-scale comparisons necessitates a chromosome-length genome assembly for cv. R108. Here, we performed in situ Hi-C (48×) to anchor, order, orient scaffolds, and correct misjoins of contigs in a previously published genome assembly (R108 v1.0), resulting in an improved genome assembly containing eight chromosome-length scaffolds that span 97.62% of the sequenced bases in the input assembly. The long-range physical information data generated using Hi-C allowed us to obtain a chromosome-length ordering of the genome assembly, better validate previous draft misjoins, and provide further insights accurately predicting synteny between A17 and R108 regions corresponding to the known chromosome 4/8 translocation. Furthermore, mapping the insertion landscape on this reference assembly presents an important resource for functional genomics by supporting efficient mutant gene identification in insertion lines. Our data provide a much-needed foundational resource that supports functional and molecular research into the Leguminosae for sustainable agriculture and feeding the future.

Kaur P ，Lui C ，Dudchenko O ，Nandety RS ，Hurgobin B ，Pham M ，Lieberman Aiden E ，Wen J ，Mysore K ... -  《INTERNATIONAL JOURNAL OF MOLECULAR SCIENCES》

Comparison of structural variants in the whole genome sequences of two Medicago truncatula ecotypes: Jemalong A17 and R108.

Structural variants (SVs) constitute a large proportion of the genomic variation that results in phenotypic variation in plants. However, they are still a largely unexplored feature in most plant genomes. Here, we present the whole-genome landscape of SVs between two model legume Medicago truncatula ecotypes-Jemalong A17 and R108- that have been extensively used in various legume biology studies. To catalogue SVs, we first resolved the previously published R108 genome assembly (R108 v1.0) to chromosome-scale using 124 × Hi-C data, resulting in a high-quality genome assembly. The inter-chromosomal reciprocal translocations between chromosomes 4 and 8 were confirmed by performing syntenic analysis between the two genomes. Combined with the Hi-C data, it appears that these translocation events had a significant effect on chromatin organization. Using both whole-genome and short-read alignments, we identified the genomic landscape of SVs between the two genomes, some of which may account for several phenotypic differences, including their differential responses to aluminum toxicity and iron deficiency, and the development of different anthocyanin leaf markings. We also found extensive SVs within the nodule-specific cysteine-rich gene family which encodes antimicrobial peptides essential for terminal bacteroid differentiation during nitrogen-fixing symbiosis. Our results provide a near-complete R108 genome assembly and the first genomic landscape of SVs obtained by comparing two M. truncatula ecotypes. This may provide valuable genomic resources for the functional and molecular research of legume biology in the future.

Li A ，Liu A ，Wu S ，Qu K ，Hu H ，Yang J ，Shrestha N ，Liu J ，Ren G ... -  《BMC PLANT BIOLOGY》

Strategies for optimizing BioNano and Dovetail explored through a second reference quality assembly for the legume model, Medicago truncatula.

Third generation sequencing technologies, with sequencing reads in the tens- of kilo-bases, facilitate genome assembly by spanning ambiguous regions and improving continuity. This has been critical for plant genomes, which are difficult to assemble due to high repeat content, gene family expansions, segmental and tandem duplications, and polyploidy. Recently, high-throughput mapping and scaffolding strategies have further improved continuity. Together, these long-range technologies enable quality draft assemblies of complex genomes in a cost-effective and timely manner. Here, we present high quality genome assemblies of the model legume plant, Medicago truncatula (R108) using PacBio, Dovetail Chicago (hereafter, Dovetail) and BioNano technologies. To test these technologies for plant genome assembly, we generated five assemblies using all possible combinations and ordering of these three technologies in the R108 assembly. While the BioNano and Dovetail joins overlapped, they also showed complementary gains in continuity and join numbers. Both technologies spanned repetitive regions that PacBio alone was unable to bridge. Combining technologies, particularly Dovetail followed by BioNano, resulted in notable improvements compared to Dovetail or BioNano alone. A combination of PacBio, Dovetail, and BioNano was used to generate a high quality draft assembly of R108, a M. truncatula accession widely used in studies of functional genomics. As a test for the usefulness of the resulting genome sequence, the new R108 assembly was used to pinpoint breakpoints and characterize flanking sequence of a previously identified translocation between chromosomes 4 and 8, identifying more than 22.7 Mb of novel sequence not present in the earlier A17 reference assembly. Adding Dovetail followed by BioNano data yielded complementary improvements in continuity over the original PacBio assembly. This strategy proved efficient and cost-effective for developing a quality draft assembly compared to traditional reference assemblies.

Moll KM ，Zhou P ，Ramaraj T ，Fajardo D ，Devitt NP ，Sadowsky MJ ，Stupar RM ，Tiffin P ，Miller JR ，Young ND ，Silverstein KAT ，Mudge J ... -  《BMC GENOMICS》

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.

Nandety RS ，Oh S ，Lee HK ，Krom N ，Gupta R ，Mysore KS ... -  《-》

Genome variations account for different response to three mineral elements between Medicago truncatula ecotypes Jemalong A17 and R108.

Resequencing can be used to identify genome variations underpinning many morphological and physiological phenotypes. Legume model plant Medicago truncatula ecotypes Jemalong A17 (J. A17) and R108 differ in their responses to mineral toxicity of aluminum and sodium, and mineral deficiency of iron in growth medium. The difference may result from their genome variations, but no experimental evidence supports this hypothesis. A total of 12,750 structure variations, 135,045 short insertions/deletions and 764,154 single nucleotide polymorphisms were identified by resequencing the genome of R108. The suppressed expression of MtAACT that encodes a putative aluminum-induced citrate efflux transporter by deletion of partial sequence of the second intron may account for the less aluminum-induced citrate exudation and greater accumulation of aluminum in roots of R108 than in roots of J. A17, thus rendering R108 more sensitive to aluminum toxicity. The higher expression-level of MtZpt2-1 encoding a TFIIIA-related transcription factor in J. A17 than R108 under conditions of salt stress can be explained by the greater number of stress-responsive elements in its promoter sequence, thus conferring J. A17 more tolerant to salt stress than R108 plants by activating the expression of downstream stress-responsive genes. YSLs (Yellow Stripe-Likes) are involved in long-distance transport of iron in plants. We found that an YSL gene was deleted in the genome of R108 plants, thus rendering R108 less tolerance to iron deficiency than J. A17 plants. The deletion or change in several genes may account for the different responses of M. truncatula ecotypes J. A17 and R108 to mineral toxicity of aluminum and sodium as well as iron deficiency. Uncovering genome variations by resequencing is an effective method to identify different traits between species/ecotypes that are genetically related. These findings demonstrate that analyses of genome variations by resequencing can shed important light on differences in responses of M. truncatula ecotypes to abiotic stress in general and mineral stress in particular.

Wang TZ ，Tian QY ，Wang BL ，Zhao MG ，Zhang WH ... -  《BMC PLANT BIOLOGY》