DOI: 10.5281/zenodo.6458498
title: Data from: High-resolution methylome analysis uncovers stress-responsive genomic hotspots and drought-sensitive TE superfamilies in the clonal Lombardy poplar

description: The following dataset contains the processed data presented in the article<strong> "High-resolution methylome analysis uncovers stress-responsive genomic hotspots and drought-sensitive TE superfamilies in the clonal Lombardy poplar"</strong> <strong>Supplementary_methods.docx: </strong>contain detailed information for the experimental stress treatments, sequencing library preparation, sequencing and DMR calling. <strong>BedGraph files (CpG.bed, CHG.bed, CHH.bed):</strong> contain methylation levels (%) for each cytosine in the Lombardy poplar genome, in the respective sequence context. The first three columns represent the genomic coordinates of the cytosine, the 56 following columns indicate the methylation levels for each of the samples. Missing values are represented with NA (when particular cytosines were not captured by the sequencing method). <strong>DMR_annotation_Populus_nigra_Italica_after_biotic_and_abiotic_treatments.txt: </strong>contains all the identified regions that showed significant stress-induced differential methylation (DMR). The file include all annotation for each single DMR: genomic location, genomic feature, gene, TE, sequence context and stress treatment, besides other specific relevant information. <strong>sample_IDs_basic_metadata.txt</strong>: contains the sample ID and the associated metadata (stress treatment and ortet location and ID) for all samples used in the analysis. <strong>Supplementary_file_1_metadata_samples.xlsx: </strong>contains the metadata associated to each sample including: sequencing statistics before and after quality and adapter trimming, read mapping and coverage statistics, and number of interrogated cytosines on each sequence context (CpG, CHG, CHH). <strong>Supplementary_file_2_GO_enrichments.xlsx: </strong>contains the complete results for the GO enrichment analysis for different gene datasets associated to: drought-CHH-DMRs, SINEs, MITEs, SINEs + MITEs. <strong>italica_denovo_TE_280920.gff</strong>: contains the predicted TEs using the following methodology. First, TEs were de-novo annotated using the Extensive de-novo TE Annotator (EDTA) (version 1.8.3) (https://github.com/oushujun/EDTA) with default parameters, except for option --sensitive: 1, which uses RepeatModeler (version 2.0.1) to identify remaining TEs. All the steps in EDTA pipeline were selected (filter, final and anno) in order to perform whole-genome annotation/analysis after the TE library was constructed. Then, in the annotated library from EDTA, we merged overlapping fragments and fragments located at a close distance (&lt;10bp) in a strand wise manner. The merged fragment was annotated as the family of longer merged fragment. Structural variants derived from nanopore data were used to redefine the boundaries of overlapping TE fragments to be more precise with actual predictions. LINE elements were identified independently by RepeatModeler in order to construct a more comprehensive de-novo TE library. <strong>SaliS.fasta:</strong> contains the consensus sequences of <strong>Sali</strong>caceae <strong>S</strong>INE families (SaliS), the file was built by extracting information from the supplementary table 2 of the publication: "Divergence of 3′ ends as a driver of short interspersed nuclear element (SINE) evolution in the Salicaceae" (https://doi.org/10.1111/tpj.14721) <strong>Pnigra_Italica_SaliS.bed: </strong>the file contains the annotated SaliS found by blastn over the P. nigra Italica reference genome (-qcov_hsp_perc 90 -perc_identity 70 -word_size 7). Column headers: chr, start, end, length, strand, perc_identity, SaliS family. <strong>Pnigra_Italica_all_TEs_for_anno.bed: </strong>contains the merged information from <strong>italica_denovo_TE_280920.gff </strong>and<strong> Pnigra_Italica_SaliS.bed.</strong> Column headers: chr, start, end, length, strand, perc_identity (only for SaliS), TE superfamily. <strong>CXX_ortet_DMRs_merged.bed</strong>: contains DMRs merged from all pairwise DMR callings between two ortets. One file per context. Column headers: chr, start, end, number of comparisons where the DMR occur, avg number of cytosines (when called in multiple DMR callings), avg differential methylation vs. control (when called in multiple DMR callings), avg adjusted p value (when called in multiple DMR callings), avg DMR length (when called in multiple DMR callings). <strong>SCRIPTS</strong> <strong>cov_filtering.sh:</strong> to filter individual positions according to a custom threshold. <strong>unionbedg_with_NAs.sh:</strong> to merge information from different samples in a single file taking into account the percentage of missing values per position across the given samples. <strong>anovas_and_contrasts_boxplots_barplots_cld.r</strong>: to perform statistical tests for the effect of treatments and ortets on the average global methylation. Each sequence context was analyzed separately. <strong>CHH_noise_filter.sh:</strong> to remove cytosines with invariable methylation values across 90% of the samples. <strong>GlobalMethAvg_calculation.r</strong>: to calculate global average methylation given a methylation file (CpG.bed, CHG/bed or CHH.bed) and sample file. <strong>Hclustering_and_PCAs_analysis.r:</strong> to perform hierarchical clustering, principal component analysis and plot the respective figures. <strong>ICC_matrices_analysis.r:</strong> to calculate intraclass correlation coefficients among all pairwise combinations and plot colored grids Annotations are based on the de novo reference genome of the Populus nigra cv. Italica clone uploaded in the ENA project: PRJEB44889 (www.ebi.ac.uk/ena/browser/view/GCA_950102115). Bisulfite sequencing data can be found under the ENA project: PRJEB51831