Transcriptome Analysis of Persian Oak ( Quercus brantii L.) Decline Using RNA-seq Technology

Since the late 1980s, the oak decline has affected the Zagros oak forests in western Iran. Persian oak ( Quercus brantii L.) the most important tree species of these forests has been damaged more than any other plant species. In the present study, the RNA sequencing technique was used for the first time to identify key genes and molecular mechanisms involved in Persian oak decline. The RNA was extracted from the leaves of healthy and declined oak trees, and sequenced using the Illumina HiSeq 2500 platform (2 × 150 bp paired-end reads). De novo transcriptome assembly of Persian oak revealed 56,743 unigenes and 6049 differentially expressed genes (DEGs) between declined and control samples. The results of gene ontology analysis showed that most of the DEGs involved in oak decline belong to the group of stress-responsive genes. In general, oak decline samples showed significant reductions in gene expression associated with “photosynthesis and storage of sugar” and “protein synthesis and related processes.” Additionally, DEGs related to the starch degradation pathway were up-regulated, whereas DEGs associated with acetate-mevalonate (MVA), biosynthesis of lignin, and lignases pathways were down-regulated. The present study's findings can be an effective step in identifying the genes involved in oak decline and deciphering the relationship between this phenomenon and biotic and abiotic stresses.