Conjoint Analysis of Genome-Wide lncRNA and mRNA Expression of Heteromorphic Leavesin Response to Environmental Heterogeneityin Populus euphratica

Heterophylly is the phenomenon of leaf forms varying along the longitudinal axis within a single plant. Populus euphratica, a heterophyllous woody plant, develops lanceolate leaves and dentate broad-ovate leaves on the bottom and top of the canopy, respectively, which are faced with different intensities of ambient solar radiation. However, the mechanism of the heteromorphic leaf response to the microenvironment in P. euphratica remains elusive. Here, we show that the dentate broad-ovate leaves have advantages in tolerating high light intensity, while lanceolate leaves are excellent at capturing light. Compared with lanceolate leaves, more trichomes, higher stomatal density, thicker lamina, and higher specific leaf weight were observed in dentate broad-ovate leaves. Furthermore, high-throughput RNA sequencing analysis revealed that the expression patterns of genes and long noncoding RNAs (lncRNAs) are different between the two heteromorphic leaves. A total of 36,492 genes and 1725 lncRNAs were detected, among which 586 genes and 54 lncRNAs were differentially expressed. Based on targets prediction, lncRNAs and target genes involved in light adaption, protein repair, stress response, and growth and development pathways were differentially expressed in heteromorphic leaves, 10 pairs of which were confirmed by quantitative real-time PCR. Additionally, the analysis of interactions indicated that lncRNA–mRNA interactions were involved in the response to the microenvironment of heteromorphic leaves. Taken together, these results suggest that the morphological features and joint regulation of lncRNA–mRNA in heteromorphic leaves may serve as survival strategies for P. euphratica, which could lead to optimal utilization of environmental factors.