R-loops are functional structures in chromatin comprising one single-stranded DNA and a DNA:RNA hybrid. Here, we report ssDRIP-seq, a single-strand DNA ligation-based library preparation technique for genome-wide identification of R-loops. When applied in Arabidopsis, ssDRIP-seq exhibits high efficiency, low bias and strand specificity. We found that Arabidopsis R-loops are enriched by both AT and GC skews, and are formed in the sense and antisense orientations. R-loops are strongly enriched in gene promoters and gene bodies, and are highly associated with noncoding RNA and repetitive genomic regions. Furthermore, R-loops are negatively correlated with CG DNA hypermethylation, and are prevalent in regions with multiple chromatin modifications, showing strong correlations with the activated and repressed gene loci. Our analyses indicate that R-loops are common features in the Arabidopsis genome and suggest that the R-loops play diverse roles in genome organization and gene regulation, thereby providing insights into plant nuclear genome formation and function.