Cotton Verticillium wilt, caused by the notorious fungal phytopathogen Verticillium dahliae (V. dahliae), is a destructive soil-borne vascular disease and severely decreases cotton yield and quality worldwide. Transcriptional and post-transcriptional regulation of genes responsive to V. dahliae are crucial for V. dahliae tolerance in plants. However, the specific microRNAs (miRNAs) and the miRNA/target gene crosstalk involved in cotton resistance to Verticillium wilt remain largely limited. To investigate the roles of regulatory RNAs under V. dahliae induction in upland cotton, mRNA and small RNA libraries were constructed from mocked and infected roots of two upland cotton cultivars with the V. dahliae-sensitive cultivar Jimian 11 (J11) and the V. dahliae-tolerant cultivar Zhongzhimian 2 (Z2). A comparative transcriptome analysis revealed 8330 transcripts were differentially expressed under V. dahliae stress and associated with several specific biological processes. Moreover, small RNA sequencing identified a total of 383 miRNAs, including 330 unique conserved miRNAs and 53 novel miRNAs. Analysis of the regulatory network involved in the response to V. dahliae stress revealed 31 differentially expressed miRNA–mRNA pairs, and the up-regulation of GhmiR395 and down-regulation of GhmiR165 were possibly involved in the response to V. dahliae by regulating sulfur assimilation through the GhmiR395-APS1/3 module and the establishment of the vascular pattern and secondary cell wall formation through GhmiR165-REV module, respectively. The integrative analysis of mRNA and miRNA expression profiles from upland cotton lays the foundation for further investigation of regulatory mechanisms of resistance to Verticillium wilt in cotton and other crops.

The NAC transcription factors (NACs) are among the largest plant-specific gene regulators and play essential roles in the transcriptional regulation of both biotic and abiotic stress responses. Verticillium wilt of cotton caused by Verticillium dahliae (V. dahliae) is a destructive soil-borne disease that severely decreases cotton yield and quality. Although NACs constitute a large family in upland cotton (G. hirsutum L.), there is little systematic investigation of the NACs’ responsive to V. dahliae that has been reported. To further explore the key NACs in response to V. dahliae resistance and obtain a better comprehension of the molecular basis of the V. dahliae stress response in cotton, a genome-wide survey was performed in this study. To investigate the roles of GhNACs under V. dahliae induction in upland cotton, mRNA libraries were constructed from mocked and infected roots of upland cotton cultivars with the V. dahliae-sensitive cultivar “Jimian 11” (J11) and V. dahliae-tolerant cultivar “Zhongzhimian 2” (Z2). A total of 271 GhNACs were identified. Genome analysis showed GhNACs phylogenetically classified into 12 subfamilies and distributed across 26 chromosomes and 20 scaffolds. A comparative transcriptome analysis revealed 54 GhNACs were differentially expressed under V. dahliae stress, suggesting a potential role of these GhNACs in disease response. Additionally, one NAC090 homolog, GhNAC204, could be a positive regulator of cotton resistance to V. dahliae infection. These results give insight into the GhNAC gene family, identify GhNACs’ responsiveness to V. dahliae infection, and provide potential molecular targets for future studies for improving V. dahliae resistance in cotton.

• The N-terminal tail of histone H3 is specifically cleaved during EV71 infection. • Viral protease 3C is identified as a protease responsible for proteolytically processing the N-terminal H3 tail. • Our finding reveals a new epigenetic regulatory mechanism for Enterovirus 71 in virus-host interactions.