Cysteine-rich receptor-like-kinases (CRKs), a transmembrane subfamily of receptor-like kinase, play crucial roles in plant adaptation. As such cotton is the major source of fiber for the textile industry, but environmental stresses are limiting its growth and production. Here, we have performed a deep computational analysis of CRKs in five Gossypium species, including G. arboreum (60 genes), G. raimondii (74 genes), G. herbaceum (65 genes), G. hirsutum (118 genes), and G. barbadense (120 genes). All identified CRKs were classified into 11 major classes and 43 subclasses with the finding of several novel CRK-associated domains including ALMT, FUSC_2, Cript, FYVE, and Pkinase. Of these, DUF26_DUF26_Pkinase_Tyr was common and had elevated expression under different biotic and abiotic stresses. Moreover, the 35 land plants comparison identified several new CRKs domain-architectures. Likewise, several SNPs and InDels were observed in CLCuD resistant G. hirsutum. The miRNA target side prediction and their expression profiling in different tissues predicted miR172 as a major CRK regulating miR. The expression profiling of CRKs identified multiple clusters with co-expression under certain stress conditions. The expression analysis under CLCuD highlighted the role of GhCRK057, GhCRK059, GhCRK058, and GhCRK081 in resistant accession. Overall, these results provided primary data for future potential functional analysis as well as a reference study for other agronomically important crops.

Cotton (Gossypium hirsutum) is an economically potent crop in many countries including Pakistan, India, and China. For the last three decades, cotton production is under the constant stress of cotton leaf curl disease (CLCuD) caused by begomoviruses/satellites complex that is transmitted through the insect pest, whitefly (Bemisia tabaci). In 2018, we identified a highly recombinant strain; Cotton leaf curl Multan virus-Rajasthan (CLCuMuV-Raj), associated with the Cotton leaf curl Multan betasatellite-Vehari (CLCuMuB^Veh). This strain is dominant in cotton-growing hub areas of central Punjab, Pakistan, causing the third epidemic of CLCuD. In the present study, we have explored the CLCuD diversity from central to southern districts of Punjab (Faisalabad, Lodhran, Bahawalpur, Rahimyar Khan) and the major cotton-growing region of Sindh (Tandojam), Pakistan for 2 years (2020–2021). Interestingly, we found same virus (CLCuMuV-Raj) and associated betasatellite (CLCuMuB^Veh) strain that was previously reported with the third epidemic in the central Punjab region. Furthermore, we found minor mutations in two genes of CLCuMuV-Raj C4 and C1 in 2020 and 2021 respectively as compared to its isolates in 2018, which exhibited virus evolution. Surprisingly, we did not find these mutations in CLCuMuV-Raj isolates identified from Sindh province. The findings of the current study represent the stability of CLCuMuV-Raj and its spread toward the Sindh province where previously Cotton leaf curl Kokhran virus (CLCuKoV) and Cotton leaf curl Shahdadpur virus (CLCuShV) have been reported. The findings of the current study demand future research on CLCuD complex to explore the possible reasons for prevalence in the field and how the virus-host-vector compatible interaction can be broken to develop resistant cultivars.