Genome-wide association study and transcriptome analysis reveal the genetic basis underlying the environmental adaptation of plant height in a woody plant

Studies on plant height have been conducted in several crops. However, the underlying genetic mechanisms in woody plants remain unclear. To improve the genetic understanding of plant height, the genome-wide association study (GWAS) was conducted on the 298 individuals of paper mulberry (Broussonetia papyrifera), and the individuals with the highest and lowest plant heights were selected for comparative transcriptome analysis. The analysis of phenotypic data showed that plant height decreased from low latitude (N: 24°30′) to high latitude (N: 41°00′), ranging from 372 to 150 cm. Furthermore, the plant height of paper mulberry was significantly correlated with environmental factors, such as latitude, frost-free period, hours of sunshine and so on, indicating adaptive phenotypic divergence across environmental gradients. A total of 228 candidate genes were identified through the GWAS, including three genes (Bp10g0547, Bp10g0551 and Bp10g0817) that contained nonsynonymous SNP variations significantly affecting plant height. A total of 2554 differentially expressed genes (DEGs) were identified through RNA sequencing (RNA-seq) analysis, including 28, 5, 3, 20 and 138 DEGs involved in auxin, gibberellin, cytokinin, ubiquitylation and transcription factors, respectively. Besides, there were 13 common genes identified by integrating GWAS and RNA-seq analysis, including Bp10g0817, which encodes COP1 (CONSTITUTIVELY PHOTOMORPHOGENIC 1) and belongs to the RING type E3 ubiquitin ligase gene family. Collectively, this study provides valuable insights into the genetic mechanisms underlying plant height and adaptation of woody plants to diverse environments.