Fig. 4

GhMYB2 generated phasiRNAs in different cotton tissues. Small RNA mapping results showing the tasiRNA distribution on the GhMYB2 locus in early ovule and anther from accession TM-1 (a), hypocotyl and embryogenic callus (EC) tissues from accession YS (b). Small RNA mapping comparisons were normalized by the counts per million reads (CPM) from each tissue sample. Both A and D subgenome allele for GhMYB2 and GhMYB109 were presented. Due to the very few polymorphism, small RNA mapping densities did not exhibit distinctive difference between either allele. GhMYB109 was the negative control for the miR828 targeting. Gh_D07G1901 (miR390 targeting) were predicted to be equivalent to GhphasiRNA1 in the cotton genome and therefore served as the positive controls for phasiRNA deriving genes. UBL1-Gh was the negative control for non-miRNA targeting loci