Abstract
Long intergenic non-coding RNAs (lincRNAs) are emerging as regulators of protein-coding genes (PCGs) in many plant and animal developmental processes and stress responses. In this study, we characterize the genome-wide lincRNAs in potatoes responsive to a vascular bacterial disease presumably caused by Candidatus Liberibacter solanacearum (CLso). Approximately 4397 lincRNAs were detected in healthy and infected potato plants at various stages of zebra chip (ZC) disease progression. Of them, ~65% (2844) were novel lincRNAs, and less than 1% (9) were orthologs of Arabidopsis and rice based on reciprocal BLAST analysis, suggesting species-specific expansion. Among the proximal lincRNAs within 50 kbp from a PCG, ~49% were transcribed from the same strand, while ~39% and ~15% followed convergent (head-to-head) and divergent (tail-to-tail) orientations, respectively. Approximately 30% (1308) were differentially expressed following CLso infection, with substantial changes occurring 21 days after infection (DAI). Weighted Gene Co-expression Network Analysis (WGCNA) of lincRNAs and PCGs identified 46 highly correlated lincRNA-PCG pairs exhibiting co-up or co-downregulation. Furthermore, overexpression of selected lincRNAs in transgenic potato hairy roots resulted in perturbation of neighboring PCG expression and conferred tolerance to CLso infection. Our results provide novel insights into potato lincRNAs' identity, expression dynamics, and functional relevance to CLso infection.