- Data note
- Open Access
Whole genome bisulfite sequencing of Medicago truncatula A17 wild type and lss mutants
BMC Research Notes volume 13, Article number: 192 (2020)
Earlier work in our lab identified a spontaneous mutant (likesunnsupernodulator-lss) in Medicago truncatula, resulting in increased nodulation. Molecular genetic evidence indicated the phenotype was due to an unknown lesion resulting in cis-silencing of the SUNN gene. Altered methylation of the promoter was suspected, but analysis of the SUNN promoter by bisulfite sequencing at the time of publication revealed no significant methylation differences between the SUNN promoter in wild type and lss plants. Using advances in methylome generation we compared the methylome of wild type and the lss mutant in the larger 810 kB area of the genome where lss maps.
The data show the distribution of types of methylation across the entire genome between A17 wild type and lss mutants, the number of differentially methylated cytosines between genotypes, and the overall pattern of gene methylation between genotypes. We expect the wild type data will be especially useful as a reference for other investigations of methylation using M. truncatula.
Legume plants regulate the number of nodules formed through a long distance signal transduction pathway that involves many genes . Earlier work in our lab identified a spontaneous mutant in Medicago truncatula in the Jemalong cultivar, resulting in increased nodulation . Molecular genetic evidence indicated the phenotype of the like-sunn-supernodulator (lss) mutant was due to an unknown lesion resulting in cis-silencing of the SUNN gene, which has a wild type sequence in lss mutants. Expression of SUNN in the shoots is critical to regulation of nodulation in the roots . The lesion is mapped to an 810 kilobase area of the genome on chromosome 4, including the SUNN gene, but the nature of the lesion has not been determined . Altered methylation of the promoter was suspected, but analysis of the SUNN promoter by bisulfite sequencing at the time of publication revealed no significant methylation differences between the SUNN promoter in A17 wild type and lss plants . Using genome capture combined with bi-sulfite sequencing, Satgé et al.  identified 474 regions that were differentially methylated during nodule development in M. truncatula and over 400 genes downregulated in plants with a mutant copy of the DEMETER demethylase gene. Because the lss lesion behaves like a paramutation, including reversion events , we decided to expand our bisulfite sequencing beyond the SUNN promoter and compare the methylomes of the A17 wild type to the lss mutant.
Whole genome bisulfite sequencing
Data consist of sequencing results of two bisulfite libraries made from young leaves of individual 6-week-old greenhouse grown Medicago truncatula A17 wild type and lss mutant plants (one of each). DNA was extracted using the DNeasy Plant Mini Kit (Qiagen). Bisulfite treatment of A17 and lss DNA was performed using the EZ DNA Methylation-Gold Kit (Zymo Research) according to manufacturer instructions. Whole genome bisulfite sequencing (WGBS) libraries were prepared using the Illumina TruSeq DNA Methylation library preparation kit. Library quality was assessed using an Agilent Bioanalyzer and quantitated using a Q-PCR kit from KAPA Biosystems following the manufacturer’s instructions. Sequencing was performed on an Illumina HiSeq 2500 with the HiSeq Cluster Kit v4 (1 × 125 bp single-read; see Data set 1 and Data set 2 for unaligned sequencing reads). Raw sequencing FASTQ files were trimmed to remove adapters, low-quality bases, and short reads using Trim Galore! v0.4.2 (Quality Phred score cutoff 5, min length 20). Alignment was carried out using Bismark v0.16.3  and Bowtie2 v2.2.9 with N-1 for increased sensitivity against the M. truncatula 4.0 genome. Duplicated sequences were removed with Picard MarkDuplicates v2.8.0 and alignments with MAPQ scores of 0 removed with Samtools v1.3.1. Sequencing and alignment results are summarized in data file 1 referenced in Table 1.
Methylation calling was performed using MethylDackel v0.2.1. Statistically significant differentially methylated cytosines (DMCs) were identified using the Bioconductor R package DSS v2.14.0  (p.threshold 1e−5) (Data file 2 referenced in Table 1). Differentially methylated regions (DMRs) were also identified through DSS (pct.sig 0.5, minimum number of DMCs 10, 50 bp minimum length, and DMRs merged if within 100 bp) (Data file 3 referenced in Table 1). There were 307 DMRs in the CG context and 772 DMRs in the CHG context. Displays of the number of DMCs in CG and CHG contexts, the number of DMRs identified by comparison of A17 and lss in CHG and CG contexts and the distributions of methylation across the exons and 1-kb flanking sequences for A17 and lss are graphically displayed in Data file 4 referenced in Table 1. Scripts used to analyze data and produce figures can be found at https://bitbucket.org/nfreese/medicagobseq. Bedtools v2.26.0 was used to associate DMRs with overlapping or closest gene (against M. truncatula MedtrA17_4.0 annotation) (Data file 3 referenced in Table 1) . Methylation data was visualized using the Integrated Genome Browser in the area around the receptor protein kinase SUNN (Medtr4g070970) gene (Data file 5 referenced in Table 1) . There was a single significant CG DMR within the first Medtr4g070970 exon. No DMRs were identified upstream or downstream within the CG or CHG context. The CG DMR identified within Medtr4g070970 displayed decreased methylation in the A17 sample.
The data sets were generated without biological replicates and thus any comparison of the A17 and lss data is limited by a small sample size.
Differentially methylated regions
Differentially methylated cytosines
- lsslike-sunn :
Whole genome bisulfite sequencing
Satge C, Moreau S, Sallet E, Lefort G, Auriac MC, Rembliere C, Cottret L, Gallardo K, Noirot C, Jardinaud MF, Gamas P. Reprogramming of DNA methylation is critical for nodule development in Medicago truncatula. Nat Plants. 2016;2:16166. https://doi.org/10.1038/nplants.2016.166.
Roy S, Liu W, Nandety RS, Crook A, Mysore KS, Pislariu C, Frugoli JA, Dickstein R, Udvardi M. Celebrating 20 years of genetic discoveries in legume nodulation and symbiotic nitrogen fixation. Plant Cell. 2020;32:15–41. https://doi.org/10.1105/tpc.19.00279.
Schnabel EL, Mukherjee A, Smith LS, Kassaw T, Long SR, Frugoli JA. The lss supernodulation mutant of medicago truncatula reduces expression of the SUNN gene. Plant Physiol. 2010;154:139–1402. https://doi.org/10.1104/pp.110.164889.
Penmetsa RV, Frugoli JA, Smith LS, Long SR, Cook DR. Dual genetic pathways controlling nodule number in Medicago truncatula. Plant Physiol. 2003;131:998–1008. https://doi.org/10.1104/pp.015677.
Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011;27:1571–2. https://doi.org/10.1093/bioinformatics/btr167.
Park Y, Wu H. Differential methylation analysis for BS-seq data under general experimental design. Bioinformatics. 2016;32:1446–53. https://doi.org/10.1093/bioinformatics/btw026.
Quinlan AR, Hall IM. IM, BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2. https://doi.org/10.1093/bioinformatics/btq033.
Freese NH, Norris DC, Loraine AE. Integrated genome browser: visual analytics platform for genomics. Bioinformatics. 2016;32:2089–95. https://doi.org/10.1093/bioinformatics/btw069.
Freese N, Schnabel E, Frugoli J. Whole genome bisulfite sequencing of Medicago truncatula A17 wild type and lss mutants. Open Science Framework. 2019. https://doi.org/10.17605/OSF.IO/B9CKM.
Freese N, Schnabel E, Frugoli J. Whole genome bisulfite sequencing of Medicago truncatula A17 wild type and lss mutants. INSDC. SRA SRP155259. 2019. https://identifiers.org/insdc.sra:SRP155259.
Bisulfite-Seq of Medicago truncatula: A17 leaves. NCBI Sequence Read Archive; 2019. https://identifiers.org/insdc.sra:SRX4457708.
Bisulfite-Seq of Medicago truncatula: LSS leaves. NCBI Sequence Read Archive; 2019. https://identifiers.org/insdc.sra:SRX4457709.
We acknowledge the David H. Murdock Research Institute Genomics Laboratory for generation of the sequence data.
The design, tissue collection and sequencing was supported by NSF IOS 1444461 to Frugoli and the creation of the libraries and informatics analysis by NSF NPGI Postdoctoral fellowship 1523814 to Freese.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Freese, N.H., Schnabel, E.L. & Frugoli, J.A. Whole genome bisulfite sequencing of Medicago truncatula A17 wild type and lss mutants. BMC Res Notes 13, 192 (2020). https://doi.org/10.1186/s13104-020-05012-6
- Medicago truncatula
- Autoregulation of nodulation