Plant Pathol J > Volume 40(1); 2024 > Article
Abbas, Nahid, Nawaz ul Rehman, Shaheen, and Liaquat: Assessment of Resistance Induction in Mungbean against Alternaria alternata through RNA Interference


A comprehensive survey of mungbean-growing areas was conducted to observe leaf spot disease caused by Alternaria alternata. Alternaria leaf spot symptoms were observed on the leaves. Diversity of 50 genotypes of mungbean was assessed against A. alternata and data on pathological traits was subjected to cluster analysis. The results showed that genotypes of mungbean were grouped into four clusters based on resistance parameters under the influence of disease. The principal component biplot demonstrated that all the disease-related parameters (% disease incidence, % disease intensity, lesion area, and % of infection) were strongly correlated with each other. Alt a 1 gene that is precisely found in Alternaria species and is responsible for virulence and pathogenicity. Alt a 1 gene was amplified using gene specific primers. The isolated pathogen produced similar symptoms when inoculated on mungbean and tobacco. The sequence analysis of the internal transcribed spacer (ITS) region, a 600 bp fragment amplified using specific primers, ITS1 and ITS2 showed 100% identity with A. alternata. Potato virus X (PVX)-based silencing vector expressing Alt a 1 gene was constructed to control this pathogen through RNA interference in tobacco. Out of 50 inoculated plants, 9 showed delayed onset of disease. Furthermore, to confirm our findings at molecular level semi-quantitative reverse transcriptase polymerase chain reaction was used. Both phenotypic and molecular investigation indicated that RNAi induced through the VIGS vector was efficacious in resisting the pathogen in the model host, Tobacco (Nicotiana tabacum). To the best of our knowledge, this study has been reported for the first time.

The population of the world is increasing day by day and therefore, food security is being challenged in many parts of the world. Most of the suffering is due to imbalances and malnutrition (Rasul et al., 2012). The substrates for the body’s regular physiological processes are provided by the consumption of a balanced diet of nutrients (Kiani et al., 2022). Legumes serve as a nutritional complement to cereals; they are a good source of protein, essential amino acids, and minerals, and therefore, can be a replacement for meat. Mungbean (Vigna radiata) is commonly known as a green gram in most countries, especially in Pakistan. It is a well-known pulse contributing to human food consumption. Mungbean is certainly a primary food and rich in source of protein and amino acids (Habibullah et al., 2007).
The growth of pulses such as mungbean is recognized to fulfill the demands for food. Several biotic challenges for instance pathogenic bacterial and fungal organisms pose a negative impact on mungbean production, therefore, limiting the maximum output. Some of the fungal diseases of mungbean include Alternaria leaf spot, Cercospora leaf spot, powdery mildew, and anthracnose which contribute to a range of losses in production. A wide variety of research based on epidemiological studies has revealed that Alternaria alternata spores are the most commonly recognized fungal spore in the environment (Woudenberg et al., 2015). Alternaria black rot, Alternaria leaf spot, and Alternaria brown spot are the most significant phenotypes of plant diseases caused by A. alternata spores. Alternaria leaf spot damage is more prominent on aerial parts and approximately 80% of losses were recorded (Singh, 1987). Due to the pathogen’s ability to spread via seed market globalization, long-distance airborne spore dissemination, and the effects of climate change, A. alternata represent a major risk of infection to horticulture crops all over the world (Matić et al., 2020). The first case of Alternaria leaf spot in Korea was observed when researchers observed many mungbean plants with dark-brown leaf spots. This study was effective in studying the presence and control of this disease in mungbean (Jung et al., 2019).
A variety of control measures including integrated pest management are being in use to control the mungbean disease (Batzer et al., 2022). Utilization of host resistance is one of the operative and economically approachable techniques of disease management in crops, especially mungbean. Screening of resistant sources is the first step towards the exploitation of host resistance (Yadav et al., 2014). Different valuable and reliable methods are in use for the screening of mungbean reactions towards diseases like Cercospora leaf spot, powdery mildew, and anthracnose. First, these diseases (the foliar diseases) are screened where the pressure of the disease is high and then it is tested by the artificial inoculation of pathogen spores (Kasettranan et al., 2010). Resistance sources are screened by the artificial inoculation of spores of pathogens. For Africa, Asia, and Australia, several disease-resistant mungbean lines have been introduced (Batzer et al., 2022). One of the most useful strategies for coping with fungal infection is the use of fungicides. Fungicides are chemicals that are employed to stop or get rid of fungal infections in plants or seeds (Price et al., 2015). To control plant infections brought on by phytopathogenic fungus, chemical-based fungicides have been utilized (Deresa and Diriba, 2023).
The majority of main diseases in growing crops meant for the market have been successfully controlled by various fungicide types (Lonergan et al., 2015). The damaging effects of overusing these chemical disease control methods on agricultural fields result in a decline in soil fertility (Carrascal-Hernández et al., 2022).To avoid chemical use, the natural defensive mechanism known as gene silencing is commonly employed which both protects against viral infection and controls endogenous gene expression (Koch et al., 2016; Šečić and Kogel, 2021).
RNAi is known as one of the important and promising approaches for genomic analysis. Research has shown that it has a significant impact on the improvement of crops. Without altering the expression of other genes, it permits the downregulation of gene activity. This technology can be used to exploit resistance thus preventing the economic loss due to crop pathogens every year. This technology is also useful in plant stress management (Mouyna et al., 2004). Target gene suppression through RNAi-generated signals in the plant, RNA silencing plays a crucial part in the research of host-induced gene silencing and the control of phytopathogenic fungi (Carrascal-Hernández et al., 2022). The expression of specific dsRNA sequences in the host plant and the silencing of target genes are necessary for the method of host-induced gene silencing through RNAi-generated signals in the host plant (Sang and Kim, 2020). The expression of dsRNA suggests that the target pathogen will transform the host plant by creating a hairpin dsRNA that contains its gene. Small RNAi homologous to the target RNA is produced during the target pathogen’s mRNA breakdown, protecting plants from illnesses spread by phytopathogenic fungi (Carrascal-Hernández et al., 2022; Koch et al., 2013).
Therefore keeping in view the importance of Alternaria leaf spot, a comprehensive survey was done to access the infection of A. alternata on mungbean, both morphological and molecular characterization was done and pathogenicity was confirmed as well as the current study aimed to develop an RNAi-based resistance strategy against A. alternata in model host (tobacco, Nicotiana tabacum).

Materials and Methods

Survey and sample collection in mungbean-growing areas

Different surveys of almost three growing season of different mungbean-growing areas of Faisalabad, Layyah, and Bahawalpur districts, Pakistan were done comprehensively (map attached for sampling locations) (Fig. 1A) from July to September as the disease was most prevalent in that time duration. Most of the areas have seeds sowing through the drill method. The mungbean-growing areas are plain lands as this crop needed warm climatic conditions. Thirty samples were collected from each location. After collection the infected areas of mungbean leaves were cut and leaf discs were carefully shifted to potato dextrose agar (PDA) media plates after sterilization (Devappa and Thejakumar, 2016) and incubated at 28 ± 2°C for 6 to 7 days (Nagrale et al., 2013).

Morphological and molecular identification

The characteristics of the colony, colony growth rate, the occurrence of sporangia, and the dimensions of conidia were observed after seven days of incubation. The extraction of genomic DNA involves a variety of techniques and methods (Zhang et al., 2010). Fungal DNA was isolated using the extraction buffer method (Masoodi et al., 2021).

Amplification of ITS and Alt a 1 gene

The ITS regions of all the isolates were amplified by the polymerase chain reaction (PCR) using primers ITS1 and ITS4 (White et al., 1990). Briefly, the primers ITS1 (5′-TCCGTAGGTGAACCTGCGG-3′) and ITS4 (5′-TCCTCCGCTTATTGATATGC-3′) were used. The PCR protocol followed Chen et al. (2000) with some modifications, the final PCR reaction mixture of 50 μl contained 15 μM each of the two primers, 50 ng of ITS DNA template, 100 μM dNTPs, 2 mM MgSO4, and 1 U Thermo Scientific Taq DNA polymerase (Waltham, MA, USA).
To further confirm that the fungal isolates belong to the Alternaria species, the major allergen gene Alt a 1 (Deards and Montague, 1991) was amplified. The primers Alt-for (5′-ATGCAGTTCACCACCATCGC-3′) and Alt-rev (5′-ACGAGGGTGAYGTAGGCGTC-3′) were designed according to Hong et al. (2005).

Phylogenetic analysis

Isolated ITS PCR products were sequenced from Macrogen Korea, and their sequences were specifically analyzed by BLAST to identify them. The sequences from the GenBank database that indicated a substantial similarity with both isolates during BLAST were downloaded. Phylogenetic analysis was carried out MEGA X ( was used to conduct the phylogenetic analysis. The Maximum Composite Likelihood approach was used to build a neighbor-joining tree and bootstrap confidence values were obtained after 1,000 replications.

Diversity of mungbean germplasm for resistance against A. alternata

Fifty genotypes of mungbean (14104, 14108, 14113, 14114, 14115, 14117, 14123, 14128, 14129, 14133, 14134, 14141, 14150, 14165, 14180, 14198, 14203, 14205, 14208, 14215, 14243, 14250, 14256, 14261, 14265, 14266, 14279, 14289, 14291, 14292, 14294, 14295, 14299, 14300, 14304, 14305, 14306, 14310, 14311, 14314, 14349, 14368, 14369, 14378, 14438, 14443, 14473, 14478, 14480, and 14540) were evaluated for resistance to A. alternata under field conditions. The data on various pathological traits; disease incidence, disease intensity, leaf area, lesion area and infection percentage was observed.

Construction of PVX-based VIGS vector

Specific PCR amplification and cloning in PVX

Oligonucleotide primers were designed with the addition of restriction sites. SalI restriction site G/TCGAC was added at 5′ end of Alt forward primer Alt-SalI-For (5′-GTCGACATGCAGTTCACCACCATCGC-3′) and ClaI AT/CGAT was added at 5′ end of Alt Reverse primer Alt-ClaI-Rev (5′-ATCGATACGAGGGTGAYGTAGGCGTC-3′) as our desired cloning vector PVX (VIGS pgR107) contains ClaI-SmalI-SalI.
The disease-causing genes of the virus have been deleted and the modified cDNA of the virus genome is ligated into a binary vector (pGreen 0000). A PCR fragment amplified by Alt-SalI-For and Alt-ClaI-Rev primers was cloned at respective sites in the vector pgR107for the construction of the PXV-based RNAi construct used for VIGS in tobacco plants. Double restriction digestion initially with ClaI 1× Tango buffer followed by incubation at 37°C for 1 h. A final 2× concentration of Tango buffer along with a 2-fold excess of SalI was added after completion of initial digestion.

Construction of Agrobacterium tumefaciens electro-competent cells

Luria-Berani (LB) broth medium was cultured followed by agitation for 48 h at 28°C with the appropriate antibiotics, i.e., rifampicin and tetracycline (Wise et al., 2006). The culture was further transferred to fresh LB media and incubated at 28°C until the cell attained an OD600 value of 0.5 to 1.0 (1010 cells/ml). Cells were pelleted and re-suspended thrice and were finally suspended in 1.5 ml of cold filter-sterilized glycerol, aliquots were made and kept at −80°C.

Electroporation-based clone transformation in Agrobacterium tumefaciens

Agrobacterium-competent cells were electroporated using the method reported by (Sambrook et al., 1989). The electro-cell manipulator 600 (BTX, San Diego, CA, USA) was used to carry out the electroporation.

Agrobacterium-mediated infiltration of host plants

Following the procedure outlined by Yang et al. (2000), Agrobacterium cultures were introduced into the Nicotiana tabacum and mungbean plants. A. tumefaciens strain GV3101 was electroporated with the infectivity constructs. The bacteria were suspended in 10 mm MgCl2 after being pelleted down at 4,000 rpm for 15 min. One hundred and fifty g/ml of acetosyringone was administered to the bacterial suspension. After 3 h of incubation, the cells were then infiltrated into the fresh, properly matured leaves of N. tabacum that were 3 to 4 weeks old.

Statistical analysis

The experiments were repeated and the average of data collected was subjected to cluster analysis and principal component analysis (PCA) by using statistical software SPSS version 12 and version 5.0 of Statisca (Sneath and Sokal, 1973). K-means clustering was used to perform cluster analysis while the elucidation distances was done by using Ward’s method.

Semi-quantitative reverse transcriptase polymerase chain reaction

Based on the manufacturer’s recommendations, total RNA was extracted from 10 samples (S1 to S2) and 50 more samples (PH1 to PH 50) using TRIzol (Invitrogen, Carlsbad, CA, USA). Superscript III reverse transcriptase (Invitrogen) and random primers were used for cDNA synthesis in accordance with the manufacturer’s instructions (Ali et al., 2014). β-Actin gene of Alterneria alternata was used as an internal control (Li et al., 2021).


Collection and fungus culturing through diseased leaves of mungbean. A thorough study of the various mungbean-growing regions was conducted in Punjab, Pakistan (Fig. 1A). Small, circular, and atypically shaped necrotic spots were seen (Fig. 1C). Dark to gray-brown color spots were present (Fig. 1B). The spots typically had a yellow halo that surrounded them and had brown rings surrounding their pale centers.

Morphological identification

In order to verify and identify the exact fungal pathogen the infected leaf areas sterilized with 15% household sodium hypochlorite were cultured on PDA plates (Fig. 2A). Dull white woolly mycelia growth was observed on leaf discs (Fig. 2B). Mycelia quickly matured, changed into gray and eventually turned greenish-brown to black (Fig. 2C). The initial stages of the mycelium were translucent, later it changed to a greenish-brown color. It was randomly branched, multicelled, and septate (Fig. 2D). Conidiophores were both short and long, arranged in single and small groups. On conidiophores, conidia were arranged in chains of five or more. Smooth, elongated conidia ranged in color from light olivaceous to dark brown with overall length of 19 μm and typically 5 μm wide (Fig. 2E). Based on its morphological traits, the fungus was suspected to be a species from the genus Alternaria.

Diversity of mungbean germplasm for resistance against A. alternata

Fifty genotypes of mungbean were evaluated for resistance to A. alternata under field conditions (Table 1). The data on various pathological traits were subjected to cluster analysis. The results demonstrated that 50 genotypes of mungbean were grouped into four clusters based on resistance parameters under the influence of disease. Cluster analysis showed that the genotypes in cluster 4 have maximum disease intensity, incidence, and infection percentage (Table 2). However, on the other hand, minimum disease parameters were observed in members grouped into cluster 3. The most diverse members were of cluster 4 and cluster 3 derivatives followed by members of cluster 2 and cluster 4. The most similar clusters were cluster 2 and cluster 3. Similar results were shown by the distance metric between the clusters where clusters 3 and 4 had maximum value (61.5108). The members of different clusters in the dendrogram are displayed in (Fig. 3). The PCA biplot demonstrated that all the disease-related parameters were strongly correlated with each other while they showed a weak correlation with leaf area (Fig. 4).

Molecular identification

The genomic DNA of the fungus was isolated using the extraction buffer technique and further examined using agarose gel (Fig. 5A). To further confirm and characterize the fungal isolate the ITS region was amplified. A primer pair of ITS1 and ITS4 was employed for the PCR amplification. Approximately a 600 bp DNA fragment was amplified from the DNA of fungus samples and resolved on 1% agarose gel (Fig. 5B). Alt a 1 gene amplification was carried out by using Alt-for and Alt-rev primers. The PCR products were approximately 534 bp (Fig. 5C) typical of the size of Alt a 1 gene. Selected PCR products were purified and sequenced. The nucleotide sequence from ITS region of the isolate was submitted to NCBI with GenBank accession no. OR349656.

Sequencing and phylogenetics

The sequences were analyzed using BLASTn and the results revealed 100% similarity to A. alternata isolates. The nucleotide sequences of these isolates were downloaded from the GenBank and analyzed with HA-34-ITS, the nucleotide sequence of a PCR product selected from this study. A phylogenetic relationship established through a neighbour-joining statistical approach supported the results that the isolate HA-34-ITS is an isolate of A. alternata (Fig. 6). The morphological information with both visual observations of fungus growth and microscopic analysis revealed that the fungus pathogen infecting mungbean in the field is Alternaria sp. Amplification of Alt a 1 gene specific to the Alternaria spp. further supported the findings. Furthermore, the sequence analysis of the ITS region showed that leaf spot disease of mungbean is caused by the A. alternata.

Construction of PVX-based VIGS vector

To use RNAi technique for gene silencing, PVX-based VIGS (pgR107) vector was selected. ClaI and SalI were used to ligate the PCR product in pgR107 vector. An amplified fragment of 534 bp was observed. DNA fragment specific to Alt a 1 gene was amplified using Alt-SalI-For and Alt-ClaI-Rev primers and ligated into PVX-based VIGS vector pgR 107 at SalI and ClaI site. The ligated plasmids were transformed into competent E. coli cells (DH5α) and clones were selected on kanamycin. Plasmids were isolated from selected colonies (Fig. 7A) and restriction digestion was performed using SalI and ClaI enzymes. Restriction digestion generated two fragments each corresponding to vector pgR107 (about 10,000 bp) and Alt a 1 (534 bp). A clone designated as HA-107-3 (lane 3, Fig. 7B) was selected for transformation in Agrobacterium and VIGS analysis.

Transformation of Agrobacterium

For inoculation on the plant host and the VIGS analysis, the clone HA-107-3 was transformed into A. tumefaciens strain LBA 4404 and was grown on LB plates containing kanamycin and rifampicin as a selection marker. Bacterial colonies containing HA-107-3 were confirmed through colony PCR using Alt a 1 gene primers. Colony PCR produced a 534 bp fragment representative of the size of the Alt a 1 gene. The individual colonies containing HA-107-3 plasmid were grown in LB broth containing antibiotics.

Infiltration of VIGS construct on tobacco (Nicotianatabacum)

Agrobacterium culture containing PVX VIGS construct HA-107-3 was infiltrated on 20-day-old 60 tobacco plants. After 3 days, the inoculum of Alternaria fungus was sprayed on 50 plants (labeled from PVH1 to PVH50) keeping the remaining 10 plants as control (designated as treated control i.e., TC1, TC2, TC3, TC4, TC5, TC6, TC7, TC8, TC9, and TC10). Ten more plants were used as untreated control (neither infiltrated with PVX VIGS nor the fungus inoculum i.e., UC1, UC2, UC3, UC4, UC5, UC6, UC7, UC8, UC9, and UC10). Ten plants were used as inoculated control (only fungus inoculum was applied i.e., S1, S2, S3, S4, S5, S6, S7, S8, S9, and S10). The plants were covered with polythene bags and kept under controlled conditions. The plants were monitored on alternate days. The symptoms of fungus infection started appearing on the inoculated control plants after 5 days. All the 10 plants displayed visible symptoms of fungal infection. Moving toward the 50 plants (PH1 to PH50) infiltrated with PVX VIGS and fungus inoculum, the symptoms were delayed with mild infection in some plants.
Multiple observations were derived from the infectivity analysis (Table 3). Control plants (TC1 to TC10) displayed no symptoms, TC1 is shown in (Fig. 8A). All fungal inoculated plants (S1 to S10) showed clear symptoms of fungal infection, S1 is shown in (Fig. 8C). Plants infiltrated only with PVX VIGS (UC1 to UC10) displayed no symptoms. In 50 plants infiltrated with PVX VIGS and fungus inoculum (PH1 to PH50), the symptoms were delayed up to 15 days and very mild infection was observed compare to the control plants. Out of 50 inoculated plants, 12 plants (PH1, PH6, PH7, PH9, PH11, PH14, PH15, PH30, PH33, PH41, PH45, and PH50) showed no infection and 9 plants (PH2, PH4, PH5, PH8, PH19, PH21, PH31, PH35, and PH42) showed delayed on set of disease, PH1and PH2 is shown in Fig. 8B and D, respectively. Remaining 29 plants (PH3, PH10, PH12, PH13, PH16, PH17, PH18, PH20, PH22, PH23, PH24, PH25, PH26, PH27, PH28, PH29, PH32, PH34, PH36, PH37, PH38, PH39, PH40, PH43, PH44, PH46, PH47, PH48, and PH49) displayed the symptoms of fungal infection. These results showed that the symptom development was associated with the expression of Alt a 1 gene that is precisely responsible for virulence and pathogenicity of Alterneria alternata, hence giving evidence that the target gene (Alt a 1) was silenced. Support for this evidence comes from the observation that fungal infection was reduced to a considerable extent only in plants that were infiltrated with PVX VIGS.

Expression analysis using semi-quantitative reverse transcriptase polymerase chain reaction

Semi-quantitative reverse transcriptase polymerase chain reaction was performed to further confirm the efficacy of PVX VIGS. Samples from the inoculated control on which only fungus inoculum was applied (S1, S2, S3, S4, S5, S6, S7, S8, S9, and S10) displayed high expression of Alt a 1 gene (Fig. 9A). Results from the analysis of 50 plants infiltrated with PVX VIGS were according to the phenotypic impression. Twelve plant samples (PH1, PH6, PH7, PH9, PH11, PH14, PH15, PH30, PH33, PH41, PH45, and PH50) displayed almost negligible Alt a 1 expression (Fig. 9B). While 9 plant samples (PH2, PH4, PH5, PH8, PH19, PH21, PH31, PH35, and PH42) visibly shown relatively low Alt a 1 expression (Fig. 9C). Furthermore, there was no alteration in the expression of actin throughout the analysis.


Mungbean is an economically important crop with numerous practical applications (Brar et al., 2022). Every year, economic loss occurs as a result of fungal and viral infections that lead to reduced mungbean production as a whole (Schreinemachers et al., 2015; Tak and Kumar, 2020). Fungal diseases pose the biggest threat to crops grown for human consumption (Bankole and Adebanjo, 2003; Fisher et al., 2012). The fungus A. alternata, which causes Alternaria leaf spot, can lower plant grains yield and quality (Basım et al., 2017; Nauanova, 2020). Each year, this disease reduces economic output. It is a foliar mungbean disease that can affect the crop at various growth stages.
A comprehensive survey of different mungbean-growing areas of Punjab, Pakistan was done. Leaf spot disease was noticed in various mungbean varieties. The symptoms on numerous different plants were marked with dark-brown leaf spots. There were abnormally shaped tiny, round, and irregular necrotic patches. There were concentric bands of dark to gray-brown color dots. The spots often had brown rings around their light-colored cores and a yellow halo around them. A study carried out in Korea noticed the same brown spots on the leaves of mungbean that were brown to black (Jung et al., 2019). Similar brown spots were also noticed in kidney bean plants in the fields in China (Jin et al., 2022). Previously from Pakistan, there is a lack of information available about the pathogen causing leaf spot diseases in mungbean. However, in the present study, the causal organism responsible for leaf spot disease is identified as A. alternata. The leaf spots recorded in this study were similar to the previous reports of A. alternata infection on mungbean from Korea (Jung et al., 2019) and kidney beans in China (Jin et al., 2022).
To verify and identify the accurate fungal species the infected leaf areas were cultured on PDA plates. The initial stages of the mycelium were translucent, later it changed to a greenish-brown color. It was randomly branched, multicelled, and septate. Conidiophores were both short and long, arranged in single and small groups, and typically 2 to 6 nm in length. On conidiophores, conidia were arranged in chains of five or more. Smooth, elongated conidia ranged in color from light olivaceous to dark brown. The morphological characters observed were identical to the description of A. alternata explained by Woudenberg et al. (2015) and Jin et al. (2022). Microscopic observations carried out by Aung et al. (2020) and Nagrale et al. (2013) further confirmed that the characteristics of the isolates mentioned above were similar to Alternaria sp.
Mungbean genotypes were assessed for their resistance against A. alternata. The results demonstrated that transgressive segregates could be developed between the members of cluster 4 and the members of either cluster 2 or cluster 3 to find some genotypes with enhanced resistance against A. alternata. Cluster analysis and distances between them determined diversity among germplasm (Adhikari et al., 1994, 1999; Mukherjee et al., 1999). The PC biplot showed positively closed correlation. The biplot display provides differential interactions among various traits (Gabriel, 1971; Onasanya et al., 2004). To ensure correct identification of the isolates, molecular analysis of the genomic DNA was performed (Cenis, 1992; Saleem and El-Shahir, 2022; Watanabe et al., 2010). For amplification and molecular characterization of fungal pathogens, the ITS sequences were used according to Toledo et al. (2013). A 600 bp amplified fragment was seen on the gel. The same PCR amplification analysis with the ITS region was used in the study conducted by Aloi et al. (2021) and Tarini et al. (2010). To further confirm that the fungal isolates belong to the Alternaria species, the major allergen gene Alt a 1 was amplified. PCR products were approximately 534 bp long. Results were similar to the studies carried out by (Aung et al., 2020; Gabriel et al., 2017; Hong et al., 2005; Nayyar et al., 2017; Skóra et al., 2015). RNAi technique was used for gene silencing. Agrobacterium-mediated VIGS system could potentially be effective in silencing genes (Zhou et al., 2012). Alt a 1 gene was silenced by following VIGS using PVX in tobacco plants. PVX-based VIGS (pgR107) vector was selected. This PVX-based RNAi construct when agro-infiltrated in the tobacco plant was able to induce resistance against fungus pathogens. Results showed that 12 out of the 50 plants remained symptomless, whereas 9 plants exhibited a delayed onset of fungal infection. These findings revealed confirmation that the target gene (Alt a 1) was silenced as the symptom development is associated with the expression of the Alt a 1 gene, which is specifically responsible for virulence and pathogenicity of A. alternata.
Furthermore, to confirm the phenotypic expression, semi-quantitative reverse transcriptase PCR was performed to investigate the expression of Alt a 1 gene at molecular level. S1, S2, S3, S4, S5, S6, S7, S8, S9, S10 (from inoculated control plants that were sprayed with fungus inoculum only) displayed the presence and adequate expression of Alt a 1 gene. While out of 50 plants infiltrated with (PVX(VIGS) + Inoculum), 12 plant samples (PH1, PH6, PH7, PH9, PH11, PH14, PH15, PH30, PH33, PH41, PH45, and PH50) displayed almost negligible Alt a 1 expression that is also confirmed by phenotypic expression of plant samples that was showing no infection symptoms. While 9 plant samples (PH2, PH4, PH5, PH8, PH19, PH21, PH31, PH35, and PH42) visibly showed relatively low Alt a 1 expression, these 9 plants showed delayed and minor infection symptoms apparently. β-Actin was used as control and showed no difference in the expression. Both the phenotypic results and molecular evidence hence confirm that PVX(VIGS) played a promising role in RNAi-mediated gene silencing. Conclusively, these findings suggest evidence that VIGS was effective in controlling the Alterneria alternata infection in tobacco.
The stable expression of this RNAi construct in host can provide more effective and durable resistance against A. alternata. Furthermore, as the construct is designed to silence the Alt a 1 gene, a major allergen-producing gene the construct will be helpful in minimizing or restricting the expression of allergen in plants. Therefore, could act as a model system in plant protection.


Conflicts of Interest

No potential conflict of interest relevant to this article was reported.


We the authors acknowledge the University of Agriculture, Faisalabad.

Fig. 1
Sampling and surveys of mungbean leaf spot disease. (A) The mungbean fields of districts Layyah, Faisalabad, and Bahawalpur were surveyed to check the disease infestation level. (B) Infected mungbean leaf with brown spots. (C) Clusters of infected brown spots.
Fig. 2
(A) Culturing of leaf discs on potato dextrose agar (PDA). (B) Start of fungus growth on PDA plates. White mycelium can be observed on the surface of leaf discs after three to four days. (C) Mature fungus developing on leaf disc surfaces. Growth that ranges from gray to black in color can be seen (after 7 to 8 days post-culturing of leaf discs). (D) Microscopic view of greenish-brown branched, multicelled, and septate conidiophore and conidia. (E) Smooth, elongated conidia ranged in color from light olivaceous to dark brown. Scale bars = 19 μm (D), 5 μm (E).
Fig. 3
Elucidation of members of different clusters in the dendrogram.
Fig. 4
Graphical representation of principal component analysis among the germplasm entries and disease assessment parameters (lesion area and % disease infection overlapped) in mungbean.
Fig. 5
Molecular identification. (A) Genomic DNA extraction of fungus more than 10,000 bp in size. (B) Polymerase chain reaction (PCR) amplification of internal transcribed spacer region from the genomic DNA of fungi (L, 1 kb DNA ladder; N, negative control). (C) PCR amplification of Alt a 1 gene approximately 534 bp in size.
Fig. 6
HA-34-ITS isolate’s internal transcribed spacer (ITS) region of the rRNA gene and Alternaria alternata isolates’ nucleotide sequences were used to generate a phylogenetic tree. Scale length = 0.050.
Fig. 7
Construction of PVX-based VIGS vector. (A) The recombinant plasmid isolated from bacteria approximately 10,000 bp in size. (B) Restriction digestion of plasmid pgR107 with SalI and ClaI for clone confirmation. Restriction digestion generated two fragments each corresponding to vector pgR107 (about 10,000 bp) and Alt a 1 (534 bp).
Fig. 8
Symbolic reflection of results of Infiltration of VIGS construct on tobacco plants. (A) TC1 plant sample infiltrated only with PVX VIGS showing no symptoms of infection. (B) PH1 plant sample treated with PVX(VIGS) + Inoculum showing no symptoms of fungal infection. (C) S1 plant sample sprayed with fungus inoculum only showing visible fungal infection with irregularly shaped, dark gray-brown necrotic spots. (D) PH2 treated with PVX(VIGS) + Inoculum showing delayed infection with mild symptoms.
Fig. 9
Expression Analysis using semi-quantitative reverse transcriptase polymerase chain reaction. (A) Samples from the inoculated control on which only fungus inoculum was applied S1, S2, S3, S4, S5, S6, S7, S8, S9, and S10 displayed high expression of Alt a 1 gene. No alteration was observed in the expression of β-actin gene. (B) Plant samples PH1, PH6, PH7, PH9, PH11, PH14, PH15, PH30, PH33, PH41, PH45, and PH50 displayed almost imperceptible Alt a 1 expression, the expression of the β-actin gene remained unchanged. (C) Plant samples PH2, PH4, PH5, PH8, PH19, PH21, PH31, PH35, and PH42 visibly showed relatively low Alt a 1 expression while β-actin gene’s expression remained unaffected.
Table 1
Mungbean germplasm data on various pathological parameters on which cluster analysis performed
Sr. no. Germplasm Percent disease incidence Percent disease intensity Lesion area (cm) Leaf area (cm) % of infection
1 14104 47.31 21.90 0.15 1.47 10.20
2 14108 35.53 15.82 0.13 1.54 8.44
3 14113 14.15 6.22 0.17 1.50 11.33
4 14114 65.93 25.65 0.21 1.53 13.72
5 14115 27.05 11.42 0.12 1.52 7.89
6 14117 59.23 24.84 0.17 1.55 10.96
7 14123 34.25 14.38 0.10 1.48 6.75
8 14128 54.18 24.13 0.11 1.51 7.28
9 14129 17.25 8.15 0.10 1.49 6.71
10 14133 25.46 10.28 0.11 1.53 7.18
11 14134 5.55 1.95 0.05 1.51 3.31
12 14141 34.66 7.56 0.10 1.54 6.49
13 14150 14.39 6.43 0.09 1.51 5.96
14 14165 47.15 10.31 0.11 1.48 7.43
15 14180 21.53 9.17 0.10 1.52 6.57
16 14198 74.23 29.13 0.19 1.56 12.17
17 14203 54.89 34.06 0.14 1.53 9.15
18 14205 55.41 21.84 0.13 1.47 8.84
19 14208 19.52 8.97 0.10 1.51 6.62
20 14215 7.51 2.33 0.07 1.53 4.57
21 14243 47.48 19.47 0.12 1.55 7.74
22 14250 67.81 30.65 0.21 1.54 13.63
23 14256 31.15 12.47 0.11 1.49 7.38
24 14261 44.63 21.11 0.13 1.53 8.49
25 14265 64.85 31.23 0.14 1.51 9.27
26 14266 54.23 23.14 0.12 1.49 8.05
27 14279 7.56 3.37 0.10 1.50 6.66
28 14289 14.53 4.45 0.09 1.52 5.92
29 14291 25.87 9.77 0.11 1.47 7.48
30 14292 18.25 8.96 0.10 1.53 6.53
31 14294 45.92 27.41 0.15 1.49 10.06
32 14295 69.82 29.97 0.19 1.51 12.58
33 14299 14.57 3.43 0.10 1.48 6.75
34 14300 22.13 9.44 0.15 1.52 9.86
35 14304 34.57 14.53 0.11 1.55 7.09
36 14305 58.21 26.37 0.17 1.51 11.25
37 14306 78.32 36.11 0.23 1.53 15.03
38 14310 31.62 14.72 0.13 1.48 8.78
39 14311 22.15 9.76 0.10 1.52 6.57
40 14314 7.86 2.69 0.09 1.51 5.96
41 14349 39.56 17.43 0.14 1.47 9.52
42 14368 65.24 29.86 0.19 1.49 12.75
43 14369 35.23 21.94 0.16 1.53 10.45
44 14378 15.29 6.14 0.10 1.51 6.62
45 14438 67.52 31.23 0.17 1.48 12.16
46 14443 47.23 22.32 0.15 1.51 9.93
47 14473 28.54 11.21 0.12 1.47 8.16
48 14478 15.13 5.89 0.10 1.49 6.71
49 14480 48.21 22.17 0.16 1.55 10.32
50 14540 25.78 13.33 0.12 1.51 7.94
Table 2
Cluster analysis of different pathological traits of Alternaria alternata infecting mungbean genotypes
Cluster 1 Cluster 2 Cluster 3 Cluster 4
Disease incidence 48.3653 27.7396 12.6162 68.3700
Disease intensity 23.4960 12.4307 5.8946 31.1250
Lesion area (cm) 0.1410 0.1189 0.0969 0.1888
Leaf area (cm) 1.5110 1.5100 1.5073 1.5175
% of infection 9.7477 8.0646 6.8558 13.0612
Table 3
Infectivity analysis for thes assessment of RNAi based resistance using PVX VIGS vector
Inoculum No. of plants

Treated/used Uninfected Infected
Control 10 10 0
Control + Inoculum 10 0 10
PVX(VIGS) 10 10 0
PVX(VIGS) + Inoculum 50 12 + 9a 29

a Nine plants showed delayed infection.


Adhikari, T. B., Mew, T. W. and Leach, J. E. 1999. Genotypic and pathotypic diversity in Xanthomonas oryzae pv.oryzae in Nepal. Phytopathology 89:687-694.
crossref pmid
Adhikari, T. B., Mew, T. W. and Teng, P. S. 1994. Phenotypic diversity of Xanthomonas oryzae pv. oryzae in Nepal. Plant Dis. 78:68-72.
Ali, M. A., Wieczorek, K., Kreil, D. P. and Bohlmann, H. 2014. The beet cyst nematode Heterodera schachtii modulates the expression of WRKY transcription factors in syncytia to favour its development in Arabidopsis roots. PLoS ONE 9:e102360.
crossref pmid pmc
Aloi, F., Riolo, M., Sanzani, S. M., Mincuzzi, A., Ippolito, A., Siciliano, I., Pane, A., Gullino, M. L. and Cacciola, S. O. 2021. Characterization of Alternaria species associated with heart rot of pomegranate fruit. J. Fungi 7:172.
crossref pmid pmc
Aung, S. L. L., Liu, H. F., Pei, D. F., Lu, B. B., Oo, M. M. and Deng, J. X. 2020. Morphology and molecular characterization of a fungus from the Alternaria alternata species complex causing black spots on Pyrus sinkiangensis (Koerle pear). Mycobiology 48:233-239.
crossref pmid pmc
Bankole, S. A. and Adebanjo, A. 2003. Mycotoxins in food in West Africa: current situation and possibilities of controlling it. Afr. J. Biotechnol. 2:254-263.
Basım, E., Basım, H., Abdulai, M., Baki, D. and Öztürk, N. 2017. Identification and characterization of Alternaria alternata causing leaf spot of olive tree (Olea europaea) in Turkey. Crop Prot. 92:79-88.
Batzer, J. C., Singh, A., Rairdin, A., Chiteri, K. and Mueller, D. S. 2022. Mungbean: a preview of disease management challenges for an alternative U.S. cash crop. J. Integr. Pest Manag. 13:4.
crossref pdf
Brar, A. S., Kaur, K., Sindhu, V. K., Tsolakis, N. and Srai, J. S. 2022. Sustainable water use through multiple cropping systems and precision irrigation. J. Cleaner Prod 333:130117.
Carrascal-Hernández, D. C., Flórez-López, E., Peralta-Ruiz, Y., Chaves-López, C. and Grande-Tovar, C. D. 2022. Eco-friendly biocontrol strategies of Alternaria phytopathogen fungus: a focus on gene-editing techniques. Agriculture 12:1722.
Cenis, J. L. 1992. Rapid extraction of fungal DNA for PCR amplification. Nucleic Acids Res. 20:2380.
crossref pmid pmc
Chen, Y. C., Eisner, J. D., Kattar, M. M., Rassoulian-Barrett, S. L., LaFe, K., Yarfitz, S. L., Limaye, A. P. and Cookson, B. T. 2000. Identification of medically important yeasts using PCR-based detection of DNA sequence polymorphisms in the internal transcribed spacer 2 region of the rRNA genes. J. Clin. Microbiol. 38:2302-2310.
crossref pmid pmc pdf
Deards, M. J. and Montague, A. E. 1991. Purification and characterisation of a major allergen of Alternaria alternata. Mol. Immunol 28:409-415.
crossref pmid
Deresa, E. M. and Diriba, T. F. 2023. Phytochemicals as alternative fungicides for controlling plant diseases: a comprehensive review of their efficacy, commercial representatives, advantages, challenges for adoption, and possible solutions. Heliyon 9:e13810.
crossref pmid pmc
Devappa, V. and Thejakumar, M. B. 2016. Morphological and physiological studies of Alternaria alternata causing leaf spot disease of Chilli (Capsicum annuum L.). Int. J. Appl. Pure Sci. Agric. 2:762-773.
Fisher, M. C., Henk, D. A., Briggs, C. J., Brownstein, J. S., Madoff, L. C., McCraw, S. L. and Gurr, S. J. 2012. Emerging fungal threats to animal, plant and ecosystem health. Nature 484:186-194.
crossref pmid pmc pdf
Gabriel, K. R. 1971. The biplot graphic display of matrices with application to principal component analysis. Biometrika 58:453-467.
Gabriel, M. F., Uriel, N., Teifoori, F., Postigo, I., Suñén, E. and Martínez, J. 2017. The major Alternaria alternata allergen, Alt a 1: a reliable and specific marker of fungal contamination in citrus fruits. Int. J. Food Microbiol. 257:26-30.
crossref pmid
Habibullah Abbas, M. and Shah, H. U. 2007. Proximate and mineral composition of mung bean. Sarhad J. Agric. 23:463-466.
Hong, S. G., Cramer, R. A., Lawrence, C. B. and Pryor, B. M. 2005. Alt a 1 allergen homologs from Alternaria and related taxa: analysis of phylogenetic content and secondary structure. Fungal Genet. Biol. 42:119-129.
crossref pmid
Jin, Y., Xiong, Y., Xu, C., Ren, J., Guo, Y., Zuo, Y., Zhang, Y. and Geng, X. 2022. First report of Alternaria alternata causing leaf spot on kidney bean in China. Plant Dis. 106:1531.
Jung, J. H., Ha, J., Kim, M. Y., Kim, J.-E., Son, H. and Lee, S.-H. 2019. First report of Alternaria alternata causing leaf spot on mungbean (Vigna radiata) in the Republic of Korea. Plant Dis. 103:1028.
Kasettranan, W., Somta, P. and Srinives, P. 2010. Mapping of quantitative trait loci controlling powdery mildew resistance in mungbean (Vigna radiata (L.) Wilczek). J. Crop Sci. Biotechnol. 13:155-161.
crossref pdf
Kiani, A. K., Dhuli, K., Donato, K., Aquilanti, B., Velluti, V., Matera, G., Iaconelli, A., Connelly, S. T., Bellinato, F., Gisondi, P. and Bertelli, M. 2022. Main nutritional deficiencies. J. Prev. Med. Hyg. 63(2 Suppl 3):E93-E101.
pmid pmc
Koch, A., Biedenkopf, D., Furch, A., Weber, L., Rossbach, O., Abdellatef, E., Linicus, L., Johannsmeier, J., Jelonek, L., Goesmann, A., Cardoza, V., McMillan, J., Mentzel, T. and Kogel, K.-H. 2016. An RNAi-based control of Fusarium graminearum infections through spraying of long dsRNAs involves a plant passage and is controlled by the fungal silencing machinery. PLoS Pathog. 12:e1005901.
crossref pmid pmc
Koch, A., Kumar, N., Weber, L., Keller, H., Imani, J. and Kogel, K.-H. 2013. Host-induced gene silencing of cytochrome P450 lanosterol C14α-demethylase-encoding genes confers strong resistance to Fusarium species. Proc. Natl. Acad. Sci. U. S. A 110:19324-19329.
crossref pmid pmc
Li, L., Ma, H., Zheng, F., Chen, Y., Wang, M., Jiao, C., Li, H. and Gai, Y. 2021. The transcription regulator ACTR controls ACT-toxin biosynthesis and pathogenicity in the tangerine pathotype of Alternaria alternata. Microbiol. Res. 248:126747.
crossref pmid
Lonergan, E., Pasche, J., Skoglund, L. and Burrows, M. 2015. Sensitivity of Ascochyta species infecting pea, lentil, and chickpea to boscalid, fluxapyroxad, and prothioconazole. Plant Dis. 99:1254-1260.
crossref pmid
Masoodi, K. Z., Lone, S. M. and Rasool, R. S. 2021. Genomic DNA extraction from Fungal Mycelium using the modified Cenis. 1992 method. In: Advanced methods in molecular biology and biotechnology, eds. by K. Z. Masoodi, S. M. Lone and R. S. Rasool, pp. 51-56. Academic Press, Cambridge, MA, USA.
Matić, S., Tabone, G., Garibaldi, A. and Gullino, M. L. 2020. Alternaria leaf spot caused by Alternaria species: an emerging problem on ornamental plants in Italy. Plant Dis. 104:2275-2287.
crossref pmid
Mouyna, I., Henry, C., Doering, T. L. and Latgé, J.-P. 2004. Gene silencing with RNA interference in the human pathogenic fungus Aspergillus fumigatus. FEMS Microbiol. Lett. 237:317-324.
crossref pmid
Mukherjee, A. K., Rao, A. V. S., De, R. N. and Nayak, P. 1999. Genetic diversity among slow-blasting rice genotypes. Oryza 36:70-73.
Nagrale, D. T., Gaikwad, A. P. and Sharma, L. 2013. Morphological and cultural characterization of Alternaria alternata (Fr.) Keissler blight of gerbera (Gerbera jamesonii H. Bolus ex J.D. Hook). J. Appl. Nat. Sci. 5:171-178.
crossref pdf
Nauanova, A. 2020. Resistance of barley varieties to the pathogens of leaf spotting on the artificial infectious background. Syst. Rev. Pharm. 11:1589-1592.
Nayyar, B. G., Woodward, S., Mur, L. A. J., Akram, A., Arshad, M., Naqvi, S. M. S. and Akhund, S. 2017. The incidence of Alternaria species associated with infected Sesamum indicum L. seeds from fields of the Punjab, Pakistan. Plant Pathol. J. 33:543-553.
crossref pmid pmc pdf
Onasanya, A., Sere, Y., Nwilene, F., Abo, M. E. and Akator, K. 2004. Reactions and resistance status of differential rice genotypes to rice yellow mottle virus, genus Sobemovirus in Cote d’ Ivoire. Asian J. Plant Sci. 3:718-723.
Price, P. P. 3rd, Purvis, M. A., Cai, G., Padgett, G. B., Robertson, C. L., Schneider, R. W. and Albu, S. 2015. Fungicide resistance in Cercospora kikuchii, a soybean pathogen. Plant Dis. 99:1596-1603.
crossref pmid
Rasul, F., Cheema, M. A., Sattar, A., Saleem, M. F. and Wahid, M. A. 2012. Evaluating the performance of three mungbean varieties grown under varying inter-row spacing. J. Anim. Plant Sci. 22:1030-1035.
Saleem, A. and El-Shahir, A. A. 2022. Morphological and molecular characterization of some Alternaria species isolated from tomato fruits concerning mycotoxin production and polyketide synthase genes. Plants 11:1168.
crossref pmid pmc
Sambrook, J., Fritsch, E. F. and Maniatis, T. 1989. Preparation and transformation of competent E. coli. In: Molecular cloning: a laboratory manual, eds. by J. Sambrook, E. F. Fritsch and T. Maniatis, 1:pp. 1-82. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, New York, USA.
Sang, H. and Kim, J.-I. 2020. Advanced strategies to control plant pathogenic fungi by host-induced gene silencing (HIGS) and spray-induced gene silencing (SIGS). Plant Biotechnol. Rep. 14:1-8.
crossref pdf
Schreinemachers, P., Balasubramaniam, S., Boopathi, N. M., Ha, C. V., Kenyon, L., Praneetvatakul, S., Sirijinda, A., Le, N. T., Srinivasan, R. and Wu, M.-H. 2015. Farmers’ perceptions and management of plant viruses in vegetables and legumes in tropical and subtropical Asia. Crop Prot. 75:115-123.
Šečić, E. and Kogel, K.-H. 2021. Requirements for fungal uptake of dsRNA and gene silencing in RNAi-based crop protection strategies. Curr. Opin. Biotechnol. 70:136-142.
crossref pmid
Singh, S. N. 1987. Response of chilli cultivars to Alternaria alternata and losses under field conditions. Farm Sci. J. 2:96-97.
Skóra, J., Otlewska, A., Gutarowska, B., Leszczyńska, J., Majak, I. and Stępień, Ł 2015. Production of the allergenic protein Alt a 1 by Alternaria isolates from working environments. Int. J. Environ. Res. Public Health 12:2164-2183.
crossref pmid pmc
Sneath, P. H. A. and Sokal, R. R. 1973. Numerical taxonomy: the principles and practice of numerical classification. W.F. Freeman & Co, San Francisco, CA, USA. pp. 573.
Tak, Y. and Kumar, M. 2020. Phenolics: a key defence secondary metabolite to counter biotic stress. In: Plant phenolics in sustainable agriculture, eds. by R. Lone, R. Shuab and A. N. Kamili, pp. 309-329. Springer, Singapore.
Tarini, N. M. A., Wahid, M. H., Ibrahim, F., Yasmon, A. and Djauzi, S. 2010. Development of multiplex-PCR assay for rapid detection of Candida spp. Med. J. Indones. 19:83-87.
Toledo, A. V., Simurro, M. E. and Balatti, P. A. 2013. Morphological and molecular characterization of a fungus, Hirsutella sp., isolated from planthoppers and psocids in Argentina. J. Insect Sci. 13:18.
crossref pmid pmc
Watanabe, M., Lee, K., Goto, K., Kumagai, S., Sugita-Konishi, Y. and Hara-Kudo, Y. 2010. Rapid and effective DNA extraction method with bead grinding for a large amount of fungal DNA. J. Food Prot. 73:1077-1084.
crossref pmid pdf
White, T. J., Bruns, T., Lee, S. and Taylor, J. 1990. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: PCR protocols: a guide to methods and applications, eds. by M. A. Innis, D. H. Gelfand, J. J. Sninsky and T. J. White, pp. 315-322. Academic Press, San Diego, CA, USA.
Wise, A. A., Liu, Z. and Binns, A. N. 2006. Three methods for the introduction of foreign DNA into Agrobacterium. Methods Mol. Biol. 343:43-53.
Woudenberg, J. H. C., Seidl, M. F., Groenewald, J. Z., de Vries, M., Stielow, J. B., Thomma, B. P. H. J. and Crous, P. W. 2015. Alternaria section Alternaria: species, formae speciales or pathotypes? Stud. Mycol. 82:1-21.
crossref pmid pmc
Yadav, D. L., Jaisani, P. and Pandey, R. N. 2014. Identification of sources of resistance in mungbean genotypes and influence of fungicidal application to powdery mildew epidemics. Int. J. Curr. Microbiol. Appl. Sci. 3:513-519.
Yang, Y., Li, R. and Qi, M. 2000. In vivo analysis of plant promoters and transcription factors by agroinfiltration of tobacco leaves. Plant J. 22:543-551.
crossref pmid
Zhang, Y. J., Zhang, S., Liu, X. Z., Wen, H. A. and Wang, M. 2010. A simple method of genomic DNA extraction suitable for analysis of bulk fungal strains. Lett. Appl. Microbiol. 51:114-118.
crossref pmid
Zhou, X-F, Sun, J-D, Zhao, Z., Lv, J., Wei, X.-W. and Cai, R. 2012. The feasibility analysis of PVX and TRV vectors as the VIGS tool for studying the gene function. Phys. Proc. 33:46-54.
METRICS Graph View
  • 0 Crossref
  •  0 Scopus
  • 273 View
  • 41 Download
Related articles

Editorial Office
Rm,904 (New Bldg.) The Korean Science & Technology Center 22,
Teheran-ro 7-Gil, Gangnamgu, Seoul 06130, Korea
Tel: +82-2-557-9360    Fax: +82-2-557-9361    E-mail:                

Copyright © 2024 by Korean Society of Plant Pathology.

Developed in M2PI

Close layer
prev next