A comparative analysis of small RNAs between two Upland cotton backcross inbred lines with different fiber length: Expression and distribution

Abstract

The cotton fiber is the most important raw material for the textile industry and an ideal model system for studying cell elongation. However, the genetic variation of fiber elongation in relation to miRNA is poorly understood. A high-throughput comparative RNA-seq of two lines differing in fiber length (FL) from a backcross inbred line (BIL) population of G. hirsutum × G. barbadense revealed differentially expressed (DE) miRNAs and their targets in rapidly elongating fibers. A real-time quantitative PCR analysis was further performed to validate the results. A total of 463 (including 47 DE) miRNAs were identified, and seven DE miRNAs were co-localized with seven FL quantitative trait loci (QTL) identified in the G. hirsutum × G. barbadense population. Of 82 (including 21 DE) targets identified, nine (including one DE) were also co-localized with the seven FL QTL. The relationship between the allopolyploid and its diploid ancestral species with respect to miRNAs and their targets was also characterized. These results will facilitate the understanding of the molecular genetic mechanism of fiber elongation with regards to miRNAs in cotton. Keywords: Allopolyploid Gossypium, Fiber elongation, miRNA, QTL hotspots