Genome-wide identification of long non-coding RNAs suggests a potential association with effector gene transcription in Phytophthora sojae

Numerous long non-coding RNAs (lncRNAs) identified and characterized in mammals, plants and fungi have been found to play critical regulatory roles in biological processes. However, little is known about the role of lncRNAs in oomycete plant pathogens, which cause devastating damage to the economy and ecosystems. We used strand-specific RNA sequencing (RNA-seq) to generate a computational pipeline to identify lncRNAs in Phytophthora sojae, a model oomycete plant pathogen. In total, 940 lncRNAs with 1010 isoforms were identified from RNA-seq data obtained from four representative stages of P.sojae. The lncRNAs had shorter transcript lengths, longer exon lengths, fewer numbers of exons, lower GC content and higher minimum free energy values compared with protein-coding genes. lncRNAs in P.sojae exhibited low sequence conservation amongst oomycetes and P.sojae isolates. Transcriptional data indicated that P.sojae lncRNAs tended to be transcribed in a stage-specific manner; representative lncRNAs were validated by semi-quantitative reverse transcription-polymerase chain reaction. Phytophthora sojae lncRNAs were concentrated in gene-sparse regions, and lncRNAs were associated with secreted protein and effector coding genes. The neighbouring genes of lncRNAs encoded various effector family members, and RNA-seq data revealed a correlation between the transcription level of lncRNAs and their neighbouring genes. Our results provide the first comprehensive identification of lncRNAs in oomycetes and suggest a potential association between lncRNAs and effector genes.