Molecular Characterization and Population Structure Analysis of Pepper Mild Mottle Virus Infecting Capsicum in Himachal Pradesh, India

Research Article

Ann Agric Crop Sci. 2021; 6(7): 1097.

Molecular Characterization and Population Structure Analysis of Pepper Mild Mottle Virus Infecting Capsicum in Himachal Pradesh, India

Kumari N1,2*, Patel PB1, Chaudhary S1 and Sharma PN1

¹Department of Plant Pathology, CSK HPKV, Himachal Pradesh, India

²Division of Crop Protection, ICAR-Central Institute for Subtropical Horticulture, Uttar Pradesh, India

*Corresponding author: Nidhi Kumari, Scientist, Division of Crop Protection, ICAR-Central Institute for Subtropical Horticulture, Rehmankhera, Lucknow-226101, Uttar Pradesh, India

Received: September 01, 2021; Accepted: October 06, 2021; Published: October 13, 2021

Abstract

Pepper Mild Mottle Virus (PMMoV) is posing a great threat to capsicum cultivation worldwide. Based on the ability of PMMoV in overcoming L mediated resistance in Capsicum spp., it has been categorized in to five pathotypes. The information of pathotype is crucial in framing sustainable management practices. The surveys conducted in seven districts of Himachal Pradesh (H.P) during 2016 and 2017 revealed maximum incidence of PMMoV in Kullu (88.89%) district followed by Mandi (78.57%) after confirmation through DAS-ELISA and RT-PCR. A total 16 isolates were selected for their pathotype characterization by deducing their Amino Acid (aa) sequence and their pathogenicity on differential capsicum cultivars carrying L1, L2, L+, L3 and L4 genes. None of the test isolates produced visible symptoms on differential cultivarsPI-159236 and PI-260429 carrying L3 and L4 genes, respectively. The aa sequences of all test isolates except PMMoV-16.6, 16.10 and PMMoV-17.1 were 100% similar to the P12 pathotypes available at NCBI database. Additionally, the 3-D structure of PMMoV-16.6, 16.10 and 17.1 isolates constructed using SWISSMODEL affirmed the presence of mutations in the short and long hairpin loops of viral CP corresponding to the inner surface and thus do not contribute to the elicitor activity of CP. Thus, all the present isolates were characterized as P12 pathotype. The full genome sequences of three isolates viz., PMMoV-16.7 (MN496153), 16.9 (MN496154) and 17.3 (MN734123) determined in the present study revealed their maximum similarity with PMMoV-HP1 (P12, KJ631123.1) and a Japanese isolate (P12, AB000709.2) and least with a Korean isolate (LC082100.1). The full genome-based recombination analysis did not show any significant recombination event.

Keywords: PMMoV; Capsicum; Pathotype; L gene; Coat protein

Abbreviations

PMMoV: Pepper Mild Mottle Virus; DAS-ELISA: Double Antibody Sandwich Enzyme Linked Immunosorbent Assay; RTPCR: Reverse Transcriptase Polymerase Chain Reaction; CP: Coat Protein; NCBI: National Center for Biotechnology Information

Introduction

Pepper Mild Mottle Virus (PMMoV), a Tobamovirus of Virgaviridae family is a major threat to Capsicum spp. cultivation across the world due to its wide prevalence, long term stability, highly contagious and seed borne nature. PMMoV was first described as a latent strain of Tobacco Mosaic Virus (TMV) in literature [1] but later isolated from peppers cultivated in Italy and designated as a distinct virus [2]. The occurrence of PMMoV is ubiquitous as its presence in abundance has been reported in different capsicum based processed food products, soil, drinking water, river water, sewage water, human feces yet in infective form [3,4,5]. In recent years this virus has also acquired attention from different group of workers worldwide as a potential water quality indicator [5]. PMMoV is a rigid ~312nm long rod-shaped positive sense single stranded RNA virus with ~6.3kb genome size. On pepper, this virus causes significant economic losses and the symptoms produced include mosaic and mottling on leaves and fruits, leaf puckering, thickening of vein, stunting, upward cupping of leaves, fruit deformations although sometimes remain unnoticed due to asymptomatic or mild leaf infection. PMMoV transmission takes place vertically through seed, plant to plant via contact, soil (a major source of inoculum for epidemics), and surface water, however, no insect transmission has been reported yet in case of PMMoV [6]. For the management of plant viruses, use of genetic resistance if durable is the most economic, target specific, highly efficient and best alternative approach with negligible deleterious effects to the environment [7]. The genetic resistance has been employed for plant virus management since more than 90 years and more than 200 genes conferring resistance have been identified in monocots, dicots, wild relative of crops and Arabidopsis spp. [8]. Pepper (Capsicum spp.) harbors L gene that conditions plant’s resistance against tobamoviruses. There are 4 classical L gene alleles named L1, L2, L3 and L4 which provide increasing recognition spectra and resistance against tobamoviruses [9,10]. In addition to these classical L alleles, a new allele L1a designated as Hk [11]. However, often virus mutant capable of overcoming existing resistance in host cultivar are encountered [12]. Different PMMoV pathotypes viz., P0, P1, P12, P123 and P1234 exist on the basis of its ability to overcome L1, L2, L3 and L4 mediated resistance [10]. Alleles L1, L2, L3 and L4 characterized with localization of P0 pathotype of TMV in “Bruinsma Wonder” and “Verbeterde glas”, both P0 and P1 pathotype in C. frutescense cv Tabasco, P0, P1 and P12 pathotypes in C. chinense and P0, P1, P12 and P123 strains of TMV in C. chacoense accession PI260429, respectively [13]. However, L1c allele is characterized with localization of P0 pathotype at high temperature only. In addition to these classical L alleles, a new allele L1a designated as Hk gene showing resistance to not only P0 pathotype at any temperature but also Paprika mild mottle virus Japanese strain (PaMMV-J, P1) at 24oC [11]. With the increasing numbers of full genome determination in PMMoV, different pathotypes capable of breaking L2, L3 and even L4 mediated resistance are also reported [2,14-17]. Earlier workers have shown the role of single Amino Acid (aa) mutations in the Coat Protein (CP) gene of the virus in breaking L mediated resistance in pepper, thus CP is designated as elicitor of L mediated resistance in pepper against tobamoviruses [18,19]. Based on the chimera and mutation analysis it has been suggested that there exists a hierarchical interaction between tobamoviruses and L gene due to the interaction between multiple Leucine Rich Repeats (LRR) domain’s sub-regions of L gene with different viral CP themselves or some protein complexes including them [20]. PMMoV has emerged as a great threat to the capsicum cultivation both in protected and open conditions in Himachal Pradesh (H.P.), India [21]. All the cultivars/capsicum hybrids like California wonder, Indira, Yolo wonder etc., famous among the farmers are susceptible to PMMoV. For undertaking breeding programme to introgress resistant genes in high yielding commercial cultivars as well as for the sustainable use of genetic resistance for its management, the knowledge of pathogen populations is of utmost importance. Thus, in this work we have analyzed and characterized the PMMoV populations from major capsicum growing areas of H.P. and identified their pathotypes(s) based on their pathogenic behavior on differential capsicum cultivars as well as CP gene and whole genome sequences.

Materials and Methods

Surveys and collection of virus isolates

In order to collect the virus isolates and assess the prevalence of pepper mild mottle disease, surveys of capsicum/chili crop cultivated under polyhouse and open field conditions were conducted in various districts of H.P. during 2016 and 2017. The districts surveyed included Kangra, Hamirpur, Bilaspur, Mandi, Kullu, Shimla and Solan. The plants were observed for the presence of virus like symptoms and insect vectors. The observation on disease incidence was recorded by counting the number of diseased plants in 2m² area at four locations in each polyhouse/field and per cent incidence was calculated as:

Total 97 diseased samples were collected from plants expressing virus like symptoms. The leaf samples properly wrapped in aluminum foil and packed in plastic bag with proper labeling were brought to the laboratory and stored at -20oC till future use.

Virus identification through DAS-ELISA and RT-PCR

The presence of PMMoV in the diseased samples was determined through DAS-ELISA using commercially available PMMoV specific antibodies (Bioreba, Switzerland) [22]. Further confirmation of PMMoV association with the infected samples was done through RT-PCR using CP specific primers (CPF: 5’ CCAATGGCTGACAGATTACG 3’, CPR: 5’ CAACGACAA CCCTTCGATTT 3’) [6]. Once the presence of virus has been confirmed as PMMoV, the virus culture was maintained on susceptible capsicum variety “California wonder” by inoculating through standard leaf rub method and kept in insect free greenhouse at 23±2oC following proper precautions to avoid any contamination.

Pathogenic Variability

Identification of PMMoV pathotypes

The differential varieties comprising of Capsicum spp. viz., C. annuum (Yolo wonder), C. frutescense (Tabasco), C. annuum (Doux des Landes), C. chinense (PI-159236) and C. chacoense (PI- 260429) having tobamovirus resistance genes viz., L¹, L², L+, L³ and L4, respectively were procured from The Centro de Investigacion and Tecnologia Agroalimenatria de Aragon (CITA), Spain through National Bureau of Plant Genetic Resources (NBPGR), New Delhi. Five plants of each differential variety were inoculated through standard leaf rub method at 4-5 leaf stage and 2 plants inoculated with plain water were kept as control in greenhouse at temperature 23±2oC and observed for the appearance of symptoms. The plants showing symptoms as well as those which were asymptomatic were subjected to DAS-ELISA to confirm the presence or absence of the virus.

Coat Protein gene amplification, cloning and sequencing

To identify the presence of various PMMoV strains in H.P., CP gene of each test isolates was amplified using CP specific primers [6]. The amplified regions of each test isolates were gel purified using Gel extraction kit (GeNei) following the manufacturer’s instructions. The eluted PCR products representing the CP gene of all present isolates were ligated in pGEMT- Easy vector (Promega) and E. coli strain DH5a cells were transformed using ligated product. Colony PCR was performed to confirm the presence of insert in randomly picked white colonies in case of each test isolate sample. Plasmid was isolated from recombinant transformed colonies through Alkali lysis method. Isolated plasmid was further screened by restriction digestion using restriction endonuclease EcoR1 and custom sequenced using the services of AgriGenome Labs. The sequences obtained were subjected to vecscreen to remove the vector sequences followed by NCBI nBLAST analysis. The sequences were translated through Expasy translate tool and the amino acid sequences of all the present isolates were observed for any single amino acid mutations and compared with the sequences of P12, P123 and P1234 pathotypes.

Phylogenetic analysis

The nucleotide sequences obtained after custom sequencing were first screened for vector sequences using vecscreen software and the sequences obtained were blasted using Blastn program (http:// www.ncbi.nih.gov/blast) to the confirm the identity as PMMoV at molecular level. The CP gene sequences of all the test isolates were compared with each other as well as other CP gene sequences available at NCBI databases. The multiple sequences were aligned by ClustalW program (https://www.genome.jp/tools-bin/clustalw) [23]. The CP gene nucleotide sequences of the present were translated using Expasy Translation Tool (www.expasy.org/tools/dna.html). The sequences were observed for single amino acid mutations to determine the pathotype of PMMoV isolates. The amino acid sequences of present isolates were compared with that of known PMMoV pathotypes obtained from NCBI database and aligned using ClustalW program (https://www.genome.jp/tools-bin/clustalw) [23] and observed for the single amino acid mutations. The amino acid sequences of the present isolates were used to construct the threedimensional structure of coat protein using SWISS Model and viewed using SWISS-Pdbviewer 4.1.0 to locate the position of single amino acid substitutions in the present isolates.

Determination of complete genome sequence

Three isolates PMMoV-16.7, 16.9 and 17.3 from Sundernagar (Mandi), Nagwain (Kullu) and Phagli (Shimla), representing three districts were selected at random and their complete nucleotide sequence was determined [6]. The whole genome was assembled using DNAstar package (DNASTAR® Inc., USA).

Full genome based phylogenetic and recombination analysis

Four full PMMoV genome sequences of HP isolates (3 isolates sequenced in present study and one PMMoV-HP1(KJ631123) sequenced [6] along with other 25 full genomes present in NCBI database were aligned using MUSCLE alignment and a phylogenetic tree was constructed through neighbor joining method in GENEIOUS Prime software. To predict the recombination events and breakpoints based on PMMoV full genomes, RDP4 program with RDP, BOOTSCAN, GENECONV, CHIMAERA, MAXIMUM CHI SQUARE, SISCAN, and 3SEQ recombination detection methods was used. The analysis was carried out with default options for general recombination detection options. The events which were detected with at least four methods were considered plausible evidences for recombination.

Results and Discussion

Surveys and collection of viral isolates

During 2016-17, surveys were conducted in 29 locations of Kangra, Hamirpur, Bilaspur, Mandi, Kullu, Shimla and Solan districts of HP to assess the prevalence of disease and collect the viral isolates. The diseased plants in various polyhouses exhibited symptoms like mosaic, mottling, chlorosis, upward or downward curling and deformation of leaves along with stunting of plants (Figure 1A). The fruits beared by diseased plants were deformed and had lumpy appearance with mosaic pattern on their surface making them unfit for market. Based on the symptomatology, maximum incidence of 70 per cent was observed in Kullu district, followed by Mandi with 60 per cent incidence (Figure 1B). The prevalence of mites and white fly was observed in all the polyhouses. The attack of white fly was much more in those cases where the capsicum and tomato crops were being cultivated together. The prevalence of symptoms like upward cupping of young leaves along with puckering was high where the mite attack was severe.