In Silico Recombination Analysis of DNA-A Sequence From
Total Page:16
File Type:pdf, Size:1020Kb
Int.J.Curr.Microbiol.App.Sci (2014) 3(2): 489-497 ISSN: 2319-7706 Volume 3 Number 2 (2014) pp. 489-497 http://www.ijcmas.com Original Research Article In Silico Recombination Analysis of DNA-A sequence from Begomovirus reported in India: This identified recombinant is the evolution from other viruses prevailing at different geographical region of Pakistan and China. Avinash Marwal, Rajneesh Prajapat and R K Gaur* Department of Science, Faculty of Arts, Science and Commerce, Mody Institute of Technology and Science, Lakshmangarh, Sikar - 332311, Rajasthan, India *Corresponding author A B S T R A C T K e y w o r d s The Begomovirus DNA-A genome was isolated from an ornamental plant Marigold and is identified as a new recombinant species, sharing nucleotide Begomovirus; identity with other isolates reported from China and Pakistan. The present study Marigold; remarkably suggests that the exchange of DNA-A with other begomoviruses would Recombinant; create a new disease complex posing a serious threat to agriculture crops and Ornamental. horticulture ornamental plants production. Introduction Recombination has played, and continues currently unknown as to whether the to play, a pivotal role in geminiviral sequences in particular parts of the evolution and may be contributing to the begomovirus genomes are exchangeable emergence of new forms of geminiviruses between different species and/or members because the high frequency of mixed of the same genus from different infections of begomoviruses provides an geographical locations. In addition, opportunity for the emergence of new Northern India seems to be unusually rich viruses arising from recombination among in virus biodiversity; we investigated the strains and / or species (Harrison and extent of recombination events and Robinson, 1999). In some cases, the examined their role in the evolution of recombinants exhibited a new pathogenic virus in India and its neighboring phenotype which is often more virulent countries. than the parents (Zhou et al., 1997). The earlier reports shown the molecular The specific recombination events characterization of complete genome of a including the recombination breaks and begomovirus isolated from an ornamental hot spots have not been reported so far in plant Marigold (Marwal et al., 2013) i.e. Marigold infecting begomovirus. It is also Ageratum enation virus (AEV: 489 Int.J.Curr.Microbiol.App.Sci (2014) 3(2): 489-497 KC589699). Here we are presenting and recombination scan and the manual highlighting the In Silico Recombination checking of automated analysis results. analysis approach for in depth study about Analysis was allowed by employing the nature of the virus. Thus, in order to Bonferroni correction with confidence take a step forward to find a cure against greater than 95% (P value 0.05). In RDP such viruses that causes major crop loss analysis, the length of the window was set worldwide. to 10 variable sites, and the step size was set to one nucleotide. P values were Materials and Methods estimated by randomizing the alignment 1,000 times. Recombination between divergent genomes is believed to be a major Results and Discussion mechanism by which diversity amongst Recombination positions were observed in viruses is generated (Robertson et al., the Ageratum enation virus (AEV: 1995). To detect the possibility of KC589699) sequence (Figure 1). The recombination in begomovirus, schematic sequence display is where the Recombination Detection Program (RDP) results of automated recombination scans was utilized, which is based on a pair wise are presented and it is the part of the scanning approach. It usually runs under program that is used to drive the manual Windows 95/98/NT/XP/VISTA/7 and checking of automated analysis results. couples a high degree of analysis The colored rectangles correspond to automation with an interactive and detailed sequence fragments, thus representing the graphical user interface (Posada and recombinant, major and minor parents in a Crandall, 2001). Using various graphical representation of a sequence recombination detection methods the fragments that have been potentially conclusion of recombination studies were derived through recombination from a evaluated (Posada, 2002). The sequence resembling the one named to the recombination breakpoint could be right of the rectangle. identified by using Recombination detection program [RDP], GENECONV, The RDP plot of Ageratum enation virus Maximum-Chi, BOOTSCAN, (AEV: KC589699) sequence conservation CHIMAERA, and 3SEQ methods. All displayed a graphical overview of the these methods were implemented in RDP sequence alignment that also indicates the v.3.44 (Martin et al., 2005). portion of the alignment. Within the Recombination positions were recognized sequence part of the display, individual by only RDP method and other methods nucleotides are color coded according to such as GENECONV, Bootscan, Maxchi their degree of conservation. When a and Chimera were not found suitable for recombination event selected, the toggle recombination analysis because of lowest sequence display button can be used to recombination breakpoint detection highlight nucleotide polymorphisms that accuracy. contribute to the recombination signals depicted in the plot display. As per the The Ageratum enation virus (AEV: schematic sequence display six KC589699) DNA-A sequence were recombination break point positions were subjected to recombination analyses using identified in the begomovirus infecting RDP method used to drive automated Marigold by the RDP method. 490 Int.J.Curr.Microbiol.App.Sci (2014) 3(2): 489-497 The first evidence is given in Figure 2 E-26 and region probability (MC where breakpoint begin from 2665th corrected) was 2.512 E-23. Hence [position 2751 in alignments] position and recommended recombination was ending breakpoint ends at 287th [position observed at the AV2 ORF, spanning 309 in alignments] position. Approximate completely AV1 ORF and the starting of p-value for this region was 3.281 x 10-11. AC3 ORF fragment of the sequence. The region probability (MC Uncorrected) The third recombination event was was 1.125 E-14 and region probability discovered at the downstream of the (MC corrected) was 7.652 E-12. It second recombination event in the suggested recombination in the common begomovirus DNA-A sequence, where region (IR) and AV2 ORF fragment of the breakpoint begin from 938th [position 965 sequence. The major parent was identified in alignments] position and ending as Pedilanthus leaf curl virus (JQ012790) breakpoint ends at 1042th [position 1094 in identified in India and were found alignments] position. Approximate p-value infecting Cestrum nocturnum. Whereas the for this region was 1.335 x 10-12. In this minor parent was Croton yellow vein case the major parent was identified as mosaic virus (FN678906) found infecting Ageratum enation virus (JQ911765) an ornamental plant Alcea rosea in causing disease in Papaver somniferum Pakistan. This clearly indicates that this and was reported from India (Figure 4). portion of recombinant fragment of The minor parent in the RDP plot was common region (IR) and AV2 ORF is found to be Tomato leaf curl Ranchi virus contributed from the two viruses (GQ994095) reported from India infecting prevailing at different geographical region, Tomato sp. The region probability (MC undoubtedly pointing towards the Uncorrected) was 1.963 E-15 and region begomovirus evolution. probability (MC corrected) was 1.335 E- 12. The recombination was detected in the RDP looks for regions within a sequence AC3 ORF fragment of the DNA-A alignment in which sequence pairs are sequence. sufficiently similar to suspect that they may have arisen through recombination. The fourth recombination event was The second recombination was detected discovered at the downstream of the third downstream of the first recombinant recombination event in the begomovirus sequence where breakpoint begin from DNA-A sequence, where breakpoint begin 388th [position 310 in alignments] position from 1042th [position 1095 in alignments] and ending breakpoint ends at 935th position and ending breakpoint ends at [position 962 in alignments] position 1200th [position 1258 in alignments] (Figure 3). Approximate p-value for this position. Approximate p-value for this region was 5.787 x 10-24. Here the region was 1.660 x 10-08. In this case the contribution of major parent in the RDP major parent was identified as Ageratum plot was by Tomato leaf curl Karnataka enation virus (JQ911765) causing disease virus (FJ514798) infecting Mentha viridis in Papaver somniferum and was reported in India. The minor parent was identified from India (Figure 5). The minor parent in as Croton yellow vein mosaic virus the RDP plot was found to be Tomato leaf (AJ507777) infecting Croton curl Ranchi virus (GQ994095) reported bonplandianum in India. The region from India infecting Tomato sp. The probability (MC Uncorrected) was 3.694 region probability (MC Uncorrected) was 491 Int.J.Curr.Microbiol.App.Sci (2014) 3(2): 489-497 Figure.1 Diagram of the schematic sequence display representing the RDP recombination map of the recombinant fragments for the Ageratum enation virus (AEV: KC589699) sequence. Each color/pattern represents a sequence specific of a virus. The virus genome organization is represented under the diagram, positioning the different viral genes named according to the begomovirus convention. Figure.2 An RDP pairwise identity plot for the piece of sequence from the major parent (JQ012790_Pedil)