Genome-Wide Investigation of the MiR166 Family Provides New Insights into Its Involvement in the Drought Stress Responses of Tea Plants (Camellia sinensis (L.) O. Kuntze)

Abstract

MicroRNA166 (miR166) is a highly conserved plant miRNA that plays a crucial role in plant growth and the resistance to various abiotic stresses. However, the miR166s in tea (Camellia sinensis (L.) O. Kuntze) have not been comprehensively identified and analyzed. This study identified 30 mature miR166s and twelve pre-miR166s in tea plants. An evolutionary analysis revealed that csn-miR166s originating from the 3′ arm of their precursors were more conserved than the csn-miR166s derived from the 5′ arm of their precursors. The twelve pre-miR166s in tea were divided into two groups, with csn-MIR166 Scaffold364-2 separated from the other precursors. The Mfold-based predictions indicated that the twelve csn-MIR166s formed typical and stable structures comprising a stem-loop hairpin, with minimum free energy ranging from −110.90 to −71.80 kcal/mol. An analysis of the CsMIR166 promoters detected diverse cis-acting elements, including those related to light responses, biosynthesis and metabolism, abiotic stress defenses, and hormone responses. There was no one-to-one relationship between the csn-miR166s and their targets, but most csn-miR166s targeted HD-Zip III genes. Physiological characterization of tea plants under drought stress showed that leaf water content proportionally decreased with the aggravation of drought stress. In contrast, tea leaves’ malondialdehyde (MDA) content proportionally increased. Moreover, the cleavage site of the ATHB-15-like transcript was identified according to a modified 5′ RNA ligase-mediated rapid amplification of cDNA ends. The RT-qPCR data indicated that the transcription of nine csn-miR166s was negatively correlated with their target gene.