Proxi-RIMS-seq2 applied to native microbiomes uncovers hundreds of known and novel m5C methyltransferase specificities

Abstract

Methylation patterns in bacteria can be used to study Restriction-Modification (RM) or other defense systems with novel properties. While ^m4C and ^m6A methylation is well characterized mainly through PacBio sequencing, the landscape of ^m5C methylation is under-characterized. To bridge this gap, we performed RIMS-seq2 on microbiomes composed of resolved assemblies of distinct genomes through proximity ligation. This high-throughput approach enables the identification of ^m5C methylated motifs and links them to cognate methyltransferases directly on native microbiomes without the need to isolate bacterial strains. Methylation patterns can also be identified on viral DNA and compared to host DNA, strengthening evidence for virus-host interaction. Applied to three different microbiomes, the method unveils over 1900 motifs that were deposited in REBASE. The motifs include a novel 8-base recognition site (CAT^m5CGATG) that was experimentally validated by characterizing its cognate methyltransferase. Our findings suggest that microbiomes harbor arrays of untapped ^m5C methyltransferase specificities, providing insights to bacterial biology and biotechnological applications.

Figure 1 :

A. Overview of the Proxi-RIMS-seq2 analytic steps : [1] Imbalances in MAGs indicate methylation [2] imbalance at known motifs [3] imbalance is used to de-novo identify methylated motifs [4] linking genotype (MAGs) with phenotype (methylation) associates methyltransferases to their predicted recognition motifs [5] methyltransferase activities can be cloned and expressed in-vivo in strains lacking ^m5C methylation for validation. B. Relationship between MAG length (in base pairs, x-axis) and read coverage (log10 scale, y-axis) within the genome-resolved vermicompost microbiome. Each data point represents a MAG, distinguished by the presence (green) or absence (purple) of predicted methylated motifs. Notably, MAG 17, belonging to the genus Patulibacter, exhibits the highest read coverage mapped to its consensus genome. C. Barplot representing the absolute C-to-T transition rate in read 1 (blue) and 2 (red) at all 16 NCN contexts(with N=A,T,C or G) in MAG 17. D. differential C-to-T transition rate between read 1 and 2 (imbalance value) at all 16 NCN contexts in MAG 17. E. PWM found most significantly associated with C-to-T imbalance in MAG 17 (* represent the methylated cytosine).

Figure 2 :

Imbalance profiles for selected clusters (full dataset in Supplementary Figure 2) in A. Vermicompost. B. Dental and C. TruMatrix microbiome. At an imbalance of 1% (dark red), all cytosine at the underlined position are predicted to be methylated. The x-axis corresponds to the annotated MAGs, y-axis corresponds to known recognition sites of ^m5C methyltransferases (as cataloged in REBASE) and z-axis (gray to red) indicates imbalance values, ranging from 0 to 1%. MAGs are color-coded according to their taxonomic order. Both MAGs and motifs are clustered based on their imbalance profiles. Rows are clustered using the Manhattan distance, and columns are clustered using the Minkowski distance, as implemented in Pheatmap.

Figure 3 :. Validation of a new methyltransferase specificity using Tet-assisted PacBio sequencing.

A. IPD ratios at a specific E.coli locus containing both CATCGATG and GATC motifs. B. average IPD ratio for all possible 8 mers for A modifications (^m6A, upper panel) and C modifications (^m4C or ^m5C, lower panel). Y-axis corresponds to the 3 first nucleotides (oriented 3’ 5’), x-axis corresponds to the 4 last nucleotides (orientated 5’-3’) and z-axis (color) corresponds to the average IPD ratio. C. De-novo identification of motifs (output of DiNAMO and PWM). #MotifPos/#TotalPos :