Occurrence of Cotton leaf curl Multan virus and associated betasatellites with leaf curl disease of Bhut-Jolokia chillies (Capsicum chinense Jacq.) in India.

Geminiviridae comprises the largest family of plant viruses which causes severe crop losses in India. The highest pungency chilli Bhut-Jolokia or ghost pepper (Capsicum chinense Jaqc.) hails from North-East region of India and is used in many dishes to add flavors and also for its medicinal value. However, this chilli variety is also affected by viruses leading to crop and economic losses. The present study reports the identification of begomoviruses in the infected chilli Bhut-Jolokia leaf samples collected from eight different places of North-East region (Manipur) of India. The infected leaf samples were screened for the presence of viral genome by rolling circle amplification (RCA) followed by PCR using degenerate primer pairs. The subsequent analyses using restriction fragment length polymorphism and sequencing revealed the presence of Cotton leaf curl Multan virus (CLCuMuV), and Tomato leaf curl Patna betasatellite (ToLCPaB). The findings focus on the phylogenetic relatedness, probable recombinational hot-spots and evolutionary divergence of the viral DNA sequences with the current reported begomoviral genome. To the best of our knowledge, this is the first report showing the presence of CLCuMuV, and associated non-cognate ToLCPaB with leaf curl disease of Bhut-Jolokia chillies. The study reveals potential recombination sites on both viral genome and betsatellite which, during the course of evolution, may have aided the virus to progress and successfully establish infection in chilli plants. Taken together, our results suggest a possible spread of CLCuMuV to the hitherto non-host crop in the North-East region of India.

Diversity and recombination analysis of Cotton leaf curl Multan virus: a highly emerging begomovirus in northern India.

BACKGROUND: Cotton leaf curl disease (CLCuD), caused by begomoviruses in association with satellite molecules, is a major threat to cotton production causing enormous losses to cotton crop in most of the cotton growing countries including Indian subcontinent. In this study, isolates of begomovirus and satellite molecules associated with CLCuD were collected from North India (Haryana, New Delhi). They were amplified employing rolling circle replication mechanism, cloned, sequenced and, their phylogenetic and recombination analysis was performed. RESULTS: The five Cotton leaf curl Multan virus (CLCuMuV) isolates investigated in this study showed monopartite organization of the genome typical of Old World begomoviruses. Nucleotide sequence analyses assigned them as the strains of CLCuMuV and were designated as CLCuMuV-SR13, CLCuMuV-SR14, CLCuMuV-ND14, CLCuMuV-ND15 and CLCuMuV-SR15. The genome of CLCuMuV-SR13 shared a highest level of nucleotide sequence identity (98%) with CLCuMuV (JN678804), CLCuMuV-SR14 and CLCuMuV-SR15 exhibited 96% with CLCuMuV (KM096471), while isolates CLCuMuV-ND15 and CLCuMuV-SR15 revealed 96% sequence identity with CLCuMuV (AY765253). The four betasatellite molecules investigated in this study shared 95-99% nucleotide sequence identity with Cotton leaf curl Multan betasatellite (CLCuMB) from India. The betasatellite molecules were designated as CLCuMB-SR13, CLCuMB-SR14, CLCuMB-ND14 and CLCuMB-ND15. Alphasatellite molecules in this study, designated as GLCuA-SR14, GLCuA-ND14 and GLCuA-SR15, revealed 98% identity with Guar leaf curl alphasatellite (GLCuA) reported from Pakistan. CONCLUSION: The phylogenetic and recombination studies concluded that the isolates of CLCuMuV genomes undertaken in this study have a potential recombinant origin. Remarkably, significant recombination was detected in almost all the genes with contribution of Cotton leaf curl Kokhran Virus (CLCuKoV) in IR, V1, V2, C1, C4 and C5 regions and of CLCuMuV in C2 region of CLCuMuV-SR14. CLCuKoV also donated in C2, C3 regions of CLCuMuV-ND14; V1, V2, C2 and C3 regions of CLCuMuV-ND15 and C1 of CLCuMuV-SR15. Altogether, these observations signify the uniqueness in Indian CLCuMuV isolates showing contribution of CLCuKoV in all the genes. An interesting observation was frequent identification of GLCuA in CLCuD leaf samples.

Dominance of Cotton leaf curl Multan virus-Rajasthan strain associated with third epidemic of cotton leaf curl disease in Pakistan.

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 (CLCuMuBVeh). 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 (CLCuMuBVeh) 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.

Gossypium hirsutum calmodulin-like protein (CML 11) interaction with geminivirus encoded protein using bioinformatics and molecular techniques.

Plant viruses are the most abundant destructive agents that exist in every ecosystem, causing severe diseases in multiple crops worldwide. Currently, a major gap is present in computational biology determining plant viruses interaction with its host. We lay out a strategy to extract virus-host protein interactions using various protein binding and interface methods for Geminiviridae, a second largest virus family. Using this approach, transcriptional activator protein (TrAP/C2) encoded by Cotton leaf curl Kokhran virus (CLCuKoV) and Cotton leaf curl Multan virus (CLCuMV) showed strong binding affinity with calmodulin-like (CML) protein of Gossypium hirsutum (Gh-CML11). Higher negative value for the change in Gibbs free energy between TrAP and Gh-CML11 indicated strong binding affinity. Consensus from gene ontology database and in-silico nuclear localization signal (NLS) tools identified subcellular localization of TrAP in the nucleus associated with Gh-CML11 for virus infection. Data based on interaction prediction and docking methods present evidences that full length and truncated C2 strongly binds with Gh-CML11. This computational data was further validated with molecular results collected from yeast two-hybrid, bimolecular fluorescence complementation system and pull down assay. In this work, we also show the outcomes of full length and truncated TrAP on plant machinery. This is a first extensive report to delineate a role of CML protein from cotton with begomoviruses encoded transcription activator protein.

Temporal Dynamic of the Ratio between Monopartite Begomoviruses and Their Associated Betasatellites in Plants, and Its Modulation by the Viral Gene ßC1.

The begomovirus-betasatellite complex constantly threatens crops in Asia. However, the quantitative relationship between begomoviruses and betasatellites remains largely unknown. The quantities of tobacco curly shoot virus (TbCSV) and its betasatellite (TbCSB) and their ratio varied significantly in initial infection, and thereafter, the ratio tended to become constant. The TbCSB/TbCSV ratio in agrobacteria inoculum significantly affected that in plants in the initial infection but not thereafter. Null-mutation of ßC1 that encodes a multifunctional protein important for pathogenesis in TbCSB significantly reduced the TbCSB/TbCSV ratio in plants. Viral inoculum plants with higher TbCSB/TbCSV ratios promoted whitefly transmission of the virus. The expression of AV1 encoded by TbCSV, ßC1 encoded by TbCSB and the ßC1/AV1 ratio varied significantly in the initial infection and thereafter the ratio tended to become constant. Additionally, the temporal dynamics of the ratio between another begomovirus and its betasatellite was similar to that of TbCSV and was positively regulated by ßC1. These results indicate that the ratio between monopartite begomoviruses and betasatellites tend to become constant as infection progresses, and is modulated by ßC1, but a higher betasatellite/begomovirus ratio in virally inoculated plants promotes virus transmission by whiteflies. Our findings provide novel insights into the association between begomoviruses and betasatellites.

Cotton leaf curl Multan virus differentially regulates innate antiviral immunity of whitefly (Bemisia tabaci) vector to promote cryptic species-dependent virus acquisition.

Begomoviruses represent the largest group of economically important, highly pathogenic, DNA plant viruses that contribute a substantial amount of global crop disease burden. The exclusive transmission of begomoviruses by whiteflies (Bemisia tabaci) requires them to interact and efficiently manipulate host responses at physiological, biological and molecular scales. However, the molecular mechanisms underlying complex begomovirus-whitefly interactions that consequently substantiate efficient virus transmission largely remain unknown. Previously, we found that whitefly Asia II 7 cryptic species can efficiently transmit cotton leaf curl Multan virus (CLCuMuV) while MEAM1 cryptic species is a poor carrier and incompetent vector of CLCuMuV. To investigate the potential mechanism/s that facilitate the higher acquisition of CLCuMuV by its whitefly vector (Asia II 7) and to identify novel whitefly proteins that putatively interact with CLCuMuV-AV1 (coat protein), we employed yeast two-hybrid system, bioinformatics, bimolecular fluorescence complementation, RNA interference, RT-qPCR and bioassays. We identified a total of 21 Asia II 7 proteins putatively interacting with CLCuMuV-AV1. Further analyses by molecular docking, Y2H and BiFC experiments validated the interaction between a whitefly innate immunity-related protein (BTB/POZ) and viral AV1 (coat protein). Gene transcription analysis showed that the viral infection significantly suppressed the transcription of BTB/POZ and enhanced the accumulation of CLCuMuV in Asia II 7, but not in MEAM1 cryptic species. In contrast to MEAM1, the targeted knock-down of BTB/POZ substantially reduced the ability of Asia II 7 to acquire and accumulate CLCuMuV. Additionally, antiviral immune signaling pathways (Toll, Imd, Jnk and Jak/STAT) were significantly suppressed following viral infection of Asia II 7 whiteflies. Taken together, the begomovirus CLCuMuV potentiates efficient virus accumulation in its vector B. tabaci Asia II 7 by targeting and suppressing the transcription of an innate immunity-related BTB/POZ gene and other antiviral immune responses in a cryptic species-specific manner.

Comparative transcriptome profiling reveals a network of differentially expressed genes in Asia II 7 and MEAM1 whitefly cryptic species in response to early infection of Cotton leaf curl Multan virus.

Cotton leaf curl Multan virus (CLCuMuV) is a whitefly-vectored begomovirus that poses ramping threat to several economically important crops worldwide. The differential transmission of CLCuMuV by its vector Bemisia tabaci mainly relies on the type of whitefly cryptic species. However, the molecular responses among different whitefly cryptic species in response to early CLCuMuV infection remain elusive. Here, we compared early-stage transcriptomic profiles of Asia II 7 and MEAM1 cryptic species infected by CLCuMuV. Results of Illumina sequencing revealed that after 6 and 12 h of CLCuMuV acquisition, 153 and 141 genes among viruliferous (VF) Asia II 7, while 445 and 347 genes among VF MEAM 1 whiteflies were differentially expressed compared with aviruliferous (AVF) whiteflies. The most abundant groups of differentially expressed genes (DEGs) among Asia II 7 and MEAM1 were associated with HTH-1 and zf-C2H2 classes of transcription factors (TFs), respectively. Notably, in contrast to Asia II 7, MEAM1 cryptic species displayed higher transcriptional variations with elevated immune-related responses following CLCuMuV infection. Among both cryptic species, we identified several highly responsive candidate DEGs associated with antiviral innate immunity (alpha glucosidase, LSM14-like protein B and phosphoenolpyruvate carboxykinase), lysosome (GPI-anchored protein 58) and autophagy/phagosome pathways (sequestosome-1, cathepsin F-like protease), spliceosome (heat shock protein 70), detoxification (cytochrome P450 4C1), cGMP-PKG signaling pathway (myosin heavy chain), carbohydrate metabolism (alpha-glucosidase), biological transport (mitochondrial phosphate carrier) and protein absorption and digestion (cuticle protein 8). Further validation of RNA-seq results showed that 23 of 28 selected genes exhibited concordant expression both in RT-qPCR and RNA-seq. Our findings provide vital mechanistic insights into begomovirus-whitefly interactions to understand the dynamics of differential begomovirus transmission by different whitefly cryptic species and reveal novel molecular targets for sustainable management of insect-transmitted plant viruses.

Molecular phylogenetics and evolutionary analysis of a highly recombinant begomovirus, Cotton leaf curl Multan virus, and associated satellites.

Cotton leaf curl Multan virus (CLCuMuV) and its associated satellites are a major part of the cotton leaf curl disease (CLCuD) caused by the begomovirus species complex. Despite the implementation of potential disease management strategies, the incessant resurgence of resistance-breaking variants of CLCuMuV imposes a continuous threat to cotton production. Here, we present a focused effort to map the geographical prevalence, genomic diversity, and molecular evolutionary endpoints that enhance disease complexity by facilitating the successful adaptation of CLCuMuV populations to the diversified ecosystems. Our results demonstrate that CLCuMuV populations are predominantly distributed in China, while the majority of alphasatellites and betasatellites exist in Pakistan. We demonstrate that together with frequent recombination, an uneven genetic variation mainly drives CLCuMuV and its satellite's virulence and evolvability. However, the pattern and distribution of recombination breakpoints greatly vary among viral and satellite sequences. The CLCuMuV, Cotton leaf curl Multan alphasatellite, and Cotton leaf curl Multan betasatellite populations arising from distinct regions exhibit high mutation rates. Although evolutionarily linked, these populations are independently evolving under strong purifying selection. These findings will facilitate to comprehensively understand the standing genetic variability and evolutionary patterns existing among CLCuMuV populations across major cotton-producing regions of the world.

Recombinant variants of cotton leaf curl Multan virus is associated with the breakdown of leaf curl resistance in cotton in northwestern India.

Cotton leaf curl disease (CLCuD), caused by a begomovirus species complex, is a major constraint to cotton (Gossypium hirsutum) production in northwestern India. During 2006 to 2010, a surveillance was conducted to monitor the spread of CLCuD in Haryana and Rajasthan. Six different field symptoms, upward curling, downward curling, enation, vein thickening, severe curling and mild curling were documented. Six isolates associated with these symptom types were tested positive in PCR to cotton leaf curl Rajasthan virus. The isolates were successfully transmitted through whitefly (Bemisia tabaci) at the rate up to 73.3% to the resistant cotton cultivar, RS2013. All these six isolates were further characterised based on the complete nucleotide sequences of the viral genome and the associated betasatellites. These virus isolates shared highest sequence identity (86-99%) with the cotton leaf curl Multan virus (CLCuMuV) and the associated betasatellites also shared highest sequence identity (78-92%) with cotton leaf curl Multan betasatellite (CLCuMuB). Based on the sequence identity and phylogenetic analysis of the viral genome and betasatellite, these isolates were identified as variants of CLCuMuV. Recombination analysis revealed significant recombination events in these isolates with the other cotton infecting begomoviruses. The isolate, Mo-Raj-2 has been identified as a resistant breaking strain having a major recombination in the coding regions of both viral genome and betasatellite. The natural occurrence of disease symptoms, transmission of the virus isolates through whitefly and complete genome analysis of the virus revealed the association of recombinant variant of CLCuMuV with the breakdown of resistance in cotton in Rajasthan and Haryana, the major cotton belt of India.

Transmission efficiency of Cotton leaf curl Multan virus by three cryptic species of Bemisia tabaci complex in cotton cultivars.

Cotton leaf curl Multan virus (CLCuMuV) is a serious and economically important viral disease agent in cotton and ornamental plants like Hibiscus in many regions of the world, especially in South Asia. CLCuMuV is transmitted exclusively by Bemisia tabaci cryptic species complex. This virus was recently recorded in southern China, presumably an invasion from South Asia. This study was performed to estimate the efficiency of three species of the B. tabaci whitefly complex (tentatively named as MEAM1, MED and Asia II 7, respectively) to transmit CLCuMuV and Cotton leaf curl multan virus betasatelite (CLCuMuB). Transmission assays and real-time quantitative PCR were conducted using three cultivars of cotton, Gossypium hirsutum, including 112-2, Xinhai-21 and Zhongmian-40. The results indicated that Asia II 7 was able to transmit the virus to two of the cotton cultivars, i.e. 112-2 and Xinhai-21, with the highest transmission efficiencies of 40% and 30%, respectively, but was unable to transmit the virus to the cotton cultivar Zhongmian-40. MEAM1 and MED failed to transmit CLCuMuV and CLCuMuB to any of the three cotton cultivars. After the three cryptic species of whiteflies had fed on virus-infected cotton plants for 48 h, the relative quantity of CLCuMuV in Asia II 7 was detected to be significantly higher than that in both MEAM1 and MED (P < 0.05). These results indicate that among the three species of whiteflies Asia II 7 is likely the most efficient vector for CLCuMuV and CLCuMuB in Malvaceae crops in China. Our findings provide valuable information to the control of viral diseases caused by CLCuMuV in the field.

Suppression of cotton leaf curl disease symptoms in Gossypium hirsutum through over expression of host-encoded miRNAs.

Cotton leaf curl disease (CLCuD), a major factor resulting in the enormous yield losses in cotton crop, is caused by a distinct monopartite begomovirus in association with Cotton leaf curl Multan betasatellite (CLCuMB). Micro(mi)RNAs are known to regulate gene expression in eukaryotes, including antiviral defense in plants. In a previous study, we had computationally identified a set of cotton miRNAs, which were shown to have potential targets in the genomes of Cotton leaf curl Multan virus (CLCuMuV) and CLCuMB at multiple loci. In the current study, effect of Gossypium arboreum-encoded miRNAs on the genome of CLCuMuV and CLCuMB was investigated in planta. Two computationally predicted cotton-encoded miRNAs (miR398 and miR2950) that showed potential to bind multiple Open Reading Frames (ORFs; C1, C4, V1, and non- coding intergenic region) of CLCuMuV, and (ßC1) of CLCuMB were selected. Functional validation of miR398 and miR2950 was done by overexpression approach in G. hirsutum var. HS6. A total of ten in vitro cotton plants were generated from independent events and subjected to biological and molecular analyses. Presence of the respective Precursor (pre)-miRNA was confirmed through PCR and Southern blotting, and their expression level was assessed by semi quantitative RT-PCR, Real Time quantitative PCR and northern hybridization in the PCR-positive lines. Southern hybridization revealed 2-4 copy integration of T-DNA in the genome of the transformed lines. Remarkably, expression of pre-miRNAs was shown up to 5.8-fold higher in the transgenic (T0) lines as revealed by Real Time PCR. The virus resistance was monitored following inoculation of the transgenic cotton lines with viruliferous whitefly (Bemisia tabaci) insect vector. After inoculation, four of the transgenic lines remained apparently symptom free. While a very low titre of viral DNA could be detected by Rolling circle amplification, betasatellite responsible for symptom induction could not be detected in any of the healthy looking transgenic lines. In this study for the first time, efficacy of the host (G. arboreum)-encoded miRNAs against CLCuD symptoms was experimentally demonstrated through overexpression of miR398 and miR2950 in G. hirsutum var. HS6 plants. Computational prediction of miRNAs targeting virus genome and their subsequent implication in translational inhibition or cleavage based suppression of viral mRNA via overexpression could help in generating virus resistant plants.