Genetic variations and miRNA-target interactions contribute to natural phenotypic variations in Populus

Variation in regulatory factors, including microRNAs (miRNAs), contributes to variation in quantitative and complex traits. However, in plants, variants in miRNAs and their target genes that contribute to natural phenotypic variation, and the underlying regulatory networks, remain poorly characterized. We investigated the associations and interactions of single-nucleotide polymorphisms (SNPs) in miRNAs and their target genes with phenotypes in 435 individuals from a natural population of Populus. We used RNA-seq to identify 217 miRNAs differentially expressed in a tension wood system, and identified 1196 candidate target genes; degradome sequencing confirmed 60 of the target sites. In addition, 72 miRNA-target pairs showed significant co-expression. Gene ontology (GO) term analysis showed that most of the genes in the co-regulated pairs participate in biological regulation. Genome resequencing found 5383 common SNPs (frequency ≥ 0.05) in 139 miRNAs and 31 037 SNPs in 819 target genes. Single-SNP association analyses identified 232 significant associations between wood traits (P ≤ 0.05) and SNPs in 102 miRNAs and 1387 associations with 478 target genes. Among these, 102 miRNA-target pairs associated with the same traits. Multi-SNP associations found 102 epistatic pairs associated with traits. Furthermore, a reconstructed regulatory network contained 12 significantly co-expressed pairs, including eight miRNAs and nine targets associated with traits. Lastly, both expression and genetic association showed that miR156i, miR156j, miR396a and miR6445b were involved in the formation of tension wood. This study shows that variants in miRNAs and target genes contribute to natural phenotypic variation and annotated roles and interactions of miRNAs and their target genes by genetic association analysis.