Summary The diploid wild cotton species Gossypium australe possesses excellent traits including resistance to disease and delayed gland morphogenesis, and has been successfully used for distant breeding programmes to incorporate disease resistance traits into domesticated cotton. Here, we sequenced the G. australe genome by integrating PacBio, Illumina short read, BioNano (DLS) and Hi‐C technologies, and acquired a high‐quality reference genome with a contig N50 of 1.83 Mb and a scaffold N50 of 143.60 Mb. We found that 73.5% of the G. australe genome is composed of various repeat sequences, differing from those of G. arboreum (85.39%), G. hirsutum (69.86%) and G. barbadense (69.83%). The G. australe genome showed closer collinear relationships with the genome of G. arboreum than G. raimondii and has undergone less extensive genome reorganization than the G. arboreum genome. Selection signature and transcriptomics analyses implicated multiple genes in disease resistance responses, including GauCCD7 and GauCBP1, and experiments revealed induction of both genes by Verticillium dahliae and by the plant hormones strigolactone (GR24), salicylic acid (SA) and methyl jasmonate (MeJA). Experiments using a Verticillium‐resistant domesticated G. barbadense cultivar confirmed that knockdown of the homologues of these genes caused a significant reduction in resistance against Verticillium dahliae. Moreover, knockdown of a newly identified gland‐associated gene GauGRAS1 caused a glandless phenotype in partial tissues using G. australe. The G. australe genome represents a valuable resource for cotton research and distant relative breeding as well as for understanding the evolutionary history of crop genomes.

In plants, long non-coding RNAs regulate disease resistance against fungi and other pathogens. However, the specific mechanism behind this regulation remains unclear. In this study, we identified disease resistance-related lncRNAs as well as their regulating genes and assessed their functions by infection of cotton (Gossypium) chromosome segment substitution lines with Verticillium dahliae. Our results demonstrated that lncRNA7 and its regulating gene Pectin methylesterase inhibitor 13 (GbPMEI13) positively regulated disease resistance via the silencing approach, while ectopic overexpression of GbPMEI13 in Arabidopsis (Arabidopsis thaliana) promoted growth and enhanced resistance to V. dahliae. In contrast, lncRNA2 and its regulating gene Polygalacturonase 12 (GbPG12) negatively regulated resistance to V. dahliae. We further found that fungal disease-related agents, including the pectin-derived oligogalacturonide (OG), could down-regulate the expression of lncRNA2 and GbPG12, leading to pectin accumulation. Conversely, OG up-regulated the expression of lncRNA7, which encodes a plant peptide phytosulfokine (PSK-α), which was confirmed by lncRNA7 overexpression and UPLC-MS experiments. We showed that PSK-α promoted IAA accumulation and activated GbPMEI13 expression through Auxin Response Factor 5 (ARF5). Since it is an inhibitor of pectin methylesterase (PME), GbPMEI13 promotes pectin methylation and therefore increases the resistance to V. dahliae. Consistently, we also demonstrated that GbPMEI13 inhibits the mycelial growth and spore germination of V. dahliae in vitro. In this study, we demonstrated that lncRNA7, lncRNA2, and their regulating genes modulate cell wall defense against V. dahliae via auxin-mediated signaling, providing a strategy for cotton breeding.

BackgroundNucleotide binding site (NBS) genes encode a large family of disease resistance (R) proteins in plants. The availability of genomic data of the two diploid cotton species, Gossypium arboreum and Gossypium raimondii, and the two allotetraploid cotton species, Gossypium hirsutum (TM-1) and Gossypium barbadense allow for a more comprehensive and systematic comparative study of NBS-encoding genes to elucidate the mechanisms of cotton disease resistance.ResultsBased on the genome assembly data, 246, 365, 588 and 682 NBS-encoding genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. The distribution of NBS-encoding genes among the chromosomes was nonrandom and uneven, and was tended to form clusters. Gene structure analysis showed that G. arboreum and G. hirsutum possessed a greater proportion of CN, CNL, and N genes and a lower proportion of NL, TN and TNL genes compared to that of G. raimondii and G. barbadense, while the percentages of RN and RNL genes remained relatively unchanged. The percentage changes among them were largest for TNL genes, about 7 times. Exon statistics showed that the average exon numbers per NBS gene in G. raimondii and G. barbadense were all greater than that in G. arboretum and G. hirsutum. Phylogenetic analysis revealed that the TIR-NBS genes of G. barbadense were closely related with that of G. raimondii. Sequence similarity analysis showed that diploid cotton G. arboreum possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. hirsutum, while diploid G. raimondii possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. barbadense. The synteny analysis showed that more NBS genes in G. raimondii and G. arboreum were syntenic with that in G. barbadense and G. hirsutum, respectively.ConclusionsThe structural architectures, amino acid sequence similarities and synteny of NBS-encoding genes between G. arboreum and G. hirsutum, and between G. raimondii and G. barbadense were the highest among comparisons between the diploid and allotetraploid genomes, indicating that G. hirsutum inherited more NBS-encoding genes from G. arboreum, while G. barbadense inherited more NBS-encoding genes from G. raimondii. This asymmetric evolution of NBS-encoding genes may help to explain why G. raimondii and G. barbadense are more resistant to Verticillium wilt, whereas G. arboreum and G. hirsutum are more susceptible to Verticillium wilt. The disease resistances of the allotetraploid cotton were related to their NBS-encoding genes especially in regard from which diploid progenitor they were derived, and the TNL genes may have a significant role in disease resistance to Verticillium wilt in G. raimondii and G. barbadense.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-017-3682-x) contains supplementary material, which is available to authorized users.

Circular RNAs (circRNAs), a class of recently discovered non-coding RNAs, play a role in biological and developmental processes. A recent study showed that circRNAs exist in plants and play a role in their environmental stress responses. However, cotton circRNAs and their role in Verticillium wilt response have not been identified up to now. In this study, two CSSLs (chromosome segment substitution lines) of G.barbadense introgressed into G. hirsutum, CSSL-1 and CSSL-4 (a resistant line and a susceptible line to Verticillium wilt, respectively), were inoculated with V. dahliae for RNA-seq library construction and circRNA analysis. A total of 686 novel circRNAs were identified. CSSL-1 and CSSL-4 had similar numbers of circRNAs and shared many circRNAs in common. However, CSSL-4 differentially expressed approximately twice as many circRNAs as CSSL-1, and the differential expression levels of the common circRNAs were generally higher in CSSL-1 than in CSSL-4. Moreover, two C-RRI comparisons, C-RRI-vs-C-RRM and C-RRI-vs-C-RSI, possessed a large proportion (approximately 50%) of the commonly and differentially expressed circRNAs. These results indicate that the differentially expressed circRNAs may play roles in the Verticillium wilt response in cotton. A total of 280 differentially expressed circRNAs were identified. A Gene Ontology analysis showed that most of the ‘stimulus response’ term source genes were NBS family genes, of which most were the source genes from the differentially expressed circRNAs, indicating that NBS genes may play a role in Verticillium wilt resistance and might be regulated by circRNAs in the disease-resistance process in cotton.

Elucidating the mechanism of resistance to biotic and abiotic stress is of great importance in cotton. In this study, a gene containing the NAC domain, designated GbNAC1, was identified from Gossypium barbadense L. Homologous sequence alignment indicated that GbNAC1 belongs to the TERN subgroup. GbNAC1 protein localized to the cell nucleus. GbNAC1 was expressed in roots, stems, and leaves, and was especially highly expressed in vascular bundles. Functional analysis showed that cotton resistance to Verticillium wilt was reduced when the GbNAC1 gene was silenced using the virus-induced gene silencing (VIGS) method. GbNAC1-overexpressing Arabidopsis showed enhanced resistance to Verticillium dahliae compared to wild-type. Thus, GbNAC1 is involved in the positive regulation of resistance to Verticillium wilt. In addition, analysis of GbNAC1-overexpressing Arabidopsis under different stress treatments indicated that it is involved in plant growth, development, and response to various abiotic stresses (ABA, mannitol, and NaCl). This suggests that GbNAC1 plays an important role in resistance to biotic and abiotic stresses in cotton. This study provides a foundation for further study of the function of NAC genes in cotton and other plants.