Genome-Wide Expression Profiles of Hemp (Cannabis sativa L.) in Response to Drought Stress

Drought is the main environmental factor impairing hemp growth and yield. In order to decipher the molecular responses of hemp to drought stress, transcriptome changes of droughtstressed hemp (DS1 and DS2), compared to wellwatered control hemp (CK1 and CK2), were studied with RNASeq technology. RNASeq generated 9.83, 11.30, 11.66, and 11.31M clean reads in the CK1, CK2, DS1, and DS2 libraries, respectively. A total of 1292 differentially expressed genes (DEGs), including 409 (31.66%) upregulated and 883 (68.34%) downregulated genes, were identified. The expression patterns of 12 selected genes were validated by qRT-PCR, and the results were accordant with Illumina analysis. Gene Ontology (GO) and KEGG analysis illuminated particular important biological processes and pathways, which enriched many candidate genes such as NAC, B3, peroxidase, expansin, and inositol oxygenase that may play important roles in hemp tolerance to drought. Eleven KEGG pathways were significantly influenced, the most influenced being the plant hormone signal transduction pathway with 15 differentially expressed genes. A similar expression pattern of genes involved in the abscisic acid (ABA) pathway under drought, and ABA induction, suggested that ABA is important in the drought stress response of hemp. These findings provide useful insights into the drought stress regulatory mechanism in hemp.