Skip to main content

Comprehensive analysis of little leaf disease incidence and resistance in eggplant

Abstract

Background

Little leaf disease caused by phytoplasma infection is a significant threat to eggplant (also known as brinjal) cultivation in India. This study focused on the molecular characterisation of the phytoplasma strains and insect vectors responsible for its transmission and screening of brinjal germplasm for resistance to little leaf disease.

Results

Surveys conducted across districts in the Tamil Nadu state of India during 2021–2022 showed a higher incidence of phytoplasma during the Zaid (March to June), followed by Kharif (June to November) and Rabi (November to March) seasons with mean incidence ranging from 22 to 27%. As the name indicates, phytoplasma infection results in little leaf (reduction in leaf size), excessive growth of axillary shoots, virescence, phyllody, stunted growth, leaf chlorosis and witches’ broom symptoms. PCR amplification with phytoplasma-specific primers confirmed the presence of this pathogen in all symptomatic brinjal plants and in Hishimonus phycitis (leafhopper), providing valuable insights into the role of leafhoppers in disease transmission. BLAST search and phylogenetic analysis revealed the phytoplasma strain as “Candidatus Phytoplasma trifolii”. Insect population and disease dynamics are highly influenced by environmental factors such as temperature, relative humidity and rainfall. Further, the evaluation of 22 eggplant accessions revealed immune to highly susceptible responses where over 50% of the entries were highly susceptible. Finally, additive main effect and multiplicative interaction (AMMI) and won-where biplot analyses identified G18 as a best-performing accession for little leaf resistance due to its consistent responses across multiple environments.

Conclusions

This research contributes essential information on little leaf incidence, symptoms, transmission and resistance profiles of different brinjal genotypes, which together ensure effective and sustainable management of this important disease of eggplants.

Peer Review reports

Background

Solanum melongena L., commonly called eggplant or brinjal, holds a pivotal position in solanaceous vegetable crops, cultivated across the globe, primarily in subtropical and tropical regions. China is the top producer, contributing more than 60% of the world’s brinjal production, followed by India with 22.58%, while Egypt, Turkey and Indonesia are the other significant producers [3]. Being rich in vitamins, minerals, soluble and free-reducing sugars, proteins, phenols and anthocyanin and with its remarkable productivity and accessibility have earned it the endearing mark of “poor man’s vegetable” in Asian countries [43].

However, despite its prominence, eggplant cultivation is disturbed significantly due to biotic and abiotic stresses. Among them, little leaf disease, caused by a phytoplasma, emerged as a formidable rival, inflicting substantial economic losses [27, 41]. Infected plants exhibit distinctive characteristics such as the development of undersized leaves, excessive shoot proliferation, phyllody and stunted growth [37]. In India, the disease was first documented by Thomas and Krishnaswami in 1939, followed by a series of comprehensive investigations specifically on the biological aspects of disease development [18, 19, 27].

The etiology of eggplant little leaf disease in India has been firmly established through symptomatology, electron microscopy and PCR assays [4, 9, 19]. Furthermore, weed species Datura stramonium, D. inoxia, Cannabis sativa subsp. sativa, Portulaca oleracea and P. grandiflora were identified as alternate hosts for the phytoplasma [19, 39, 44] while leafhopper Hishimonas phycitis was predicted as the likely vectors for the disease transmission [4, 17, 19].

The current approach to managing little leaf disease relies heavily on controlling its insect vector population through pesticide applications and cultural practices [25]. However, these methods have significant environmental and economic drawbacks when host resistance emerges as a promising alternative [8]. However, the process of identifying resistant genotypes across diverse environments remains challenging and requires appropriate assessment techniques and analytical tools. Utilizing disease incidence percentages [16, 19] and the apparent rate of disease progress between two observations has offered solutions for determining resistance in genotypes. Moreover, the advent of additive main effect and multiplicative interaction (AMMI) analysis has provided an effective approach for studying the performance of genotypes in diverse environmental conditions, thereby enhancing our understanding of genotype-environment interactions (G x E) [28, 35, 36, 40].

Hence, in addition to the molecular characterisation of the phytoplasma strains and insect vectors responsible for disease transmission, in this study, we screened eggplant germplasm for resistance to little leaf disease and assessed the stability of the identified resistant genotypes across multiple environments through AMMI analysis.

Results

Symptomatology and disease incidence survey

A comprehensive survey of eggplant fields across nine districts in Tamil Nadu state of India during 2021–2022 revealed varying incidences of phytoplasma infection. Based on the Percent Disease Incidence (PDI) mean value, the highest incidence (27.05%) occurred during the Zaid (March to June), followed by Kharif (June to November) (25.08%) and Rabi (November to March) (22.07%) seasons (Fig. 1). Phytoplasma-associated symptoms primarily included a drastic reduction in leaf size and excessive growth of axillary shoots, resulting in a bushy appearance of affected plants. Common symptoms observed included little leaf, virescence, phyllody, severe growth stunting, leaf chlorosis, and witches’ broom (Fig. 2A-F).

Fig. 1
figure 1

Incidence of brinjal little leaf in nine districts of Tamil Nadu state in India. Each bar is mean of three replications

Fig. 2
figure 2

Typical phytoplasma symptoms were observed in infected eggplants. (A) little leaf, (B) chlorosis, (C) severe stunting and necrosis, (D) phyllody (E) virescence and (F) giant calyx

Confirmation of phytoplasma infection through PCR amplification

Amplification of the nine individual symptomatic samples collected from nine districts of Tamil Nadu, India using phytoplasma-specific primer pairs (P1/P6) [10] amplified a 1.5 kb product (Supplementary figure S1). Subsequent evaluation using R16F2n/R16R2 [45] in the nested PCR assays yielded specific amplicons of 1.2 kb from each of the nine symptomatic eggplant plant samples (Supplementary figure S2). These results indicated the presence of phytoplasma in all evaluated symptomatic eggplant samples and sequenced for further confirmation. The sequence data was also submitted to NCBI Bank (Table 1). No amplification was observed in DNAs extracted from healthy plant samples used as negative controls in both the first round and nested PCR experiments with the same set of primers.

Table 1 Sample ID, GenBank accession numbers, Percentage ID with reference strain from data set as obtained by iPhyClassifier and ‘Candidatus Phytoplasma’ species assignment

Phytoplasma detection in insects by PCR assays

In the PCR analysis of the various insect species (10 insects per species) associated with eggplant, including Amrasca biguttula biguttula, Hishimonus phycitis, Leucinodes arbonalis, Bemisia tabaci and Henosepilachna vigintioctopunctata, 1.25 kb PCR amplicon was observed only from H. phycitis, confirming its prominent role as the vector for the disease transmission (Supplementary figure S3). The 16 S rRNA of H. phycitis was sequenced for further confirmation, and the sequence data was also submitted to GenBank (HP: Hp CBE; Accession No: ON870850). Additionally, greenhouse transmission assays were carried out utilizing positive H. phycitis populations. Following 2 months of feeding under the cage, 50% of brinjal plants in all replications exhibited characteristic little leaf disease symptoms, while control cages without leafhopper inoculation showed no signs of disease.

BLAST and phylogenetic analyses of H. phycitis 16 S rDNA sequences

BLAST analysis of the 16S rRNA partial sequences obtained from all H. phycitis samples indicated similarity with the 16Sr VI group of Candidatus Phytoplasma trifolii. Subsequently, in the phylogenetic tree analysis, the 16S rRNA sequence from the eggplant little leaf isolates BLLP-1 to BLLP-9 and the H. phycitis isolate Hp - CBE clustered together with the reference strain Candidatus Phytoplasma trifolii, derived from clover (Fig. 3).

Fig. 3
figure 3

Phylogenetic tree based on 16SrDNA constructed by neighbour-joining method showing the relationships among eggplant, leafhopper and the reference phytoplasma strains

Characteristics of H. phycitis.

Hishimonus phycitis was characterised at the species level in the Department of Agricultural Entomology, Tamil Nadu Agricultural University (TNAU), Coimbatore, India. Species-specific features consist of compound eyes black with a whitish apical margin and ocelli at the anterior margin of the vertex contiguous with compound eyes. The vertex is devoid of spots and the cypleal suture is distinct. The scutellum is triangular with black traverse suture in the middle with yellowish brown. Wings dark brown mottling spots are present all over the wing, and a black semi-circular spot is present at the commissural margin of each wing. when wings are folded at rest, this spot looks like a prominent circular spot. Pictographic representations of H. phycitis are given in (Fig. 4).

Fig. 4
figure 4

External Morphology of Hishimonus phycitis adult male: (A). Dorsal view, (B). Lateral view, (C). Ventral view, (D). Head with pronotum, dorsal view, (E). Head, Ventral view

Population dynamics of H. phycitis

Population dynamics of H. phycitis were investigated during the Kharif, Rabi, and Zaid seasons (Table 2). The population of H. phycitis reached a maximum during Zaid, followed by the Kharif and Rabi seasons. During the Kharif season, the population of H. phycitis showed positive correlations with maximum and minimum temperatures, morning and evening relative humidity, and wind speed, and negative correlations with rainfall. During the Rabi season, positive correlations with maximum temperature (28-30.5 °C) and wind speed (3.2-8 km/hr), while negative correlations with minimum temperature (19.5–27.5 °C), morning and evening relative humidity (84–91 & 45–79%) and rainfall (0.5–63 mm) were obvious. However, in the Zaid season, the population of H. phycitis displayed positive correlations with minimum temperature, evening relative humidity, rainfall, and wind speed, and negative correlations with maximum temperature and morning relative humidity (Table 3). The population dynamics of H. phycitis were also found to be influenced by weather parameters, accounting for 99, 95 and 96% of the variation (Table 4).

Table 2 Population dynamics of Hishimonus phycitis at Coimbatore during 2021–2022
Table 3 Correlation coefficient (r) and regression equation Hishimonus phycitis vs. weather parameters for Kharif, Rabi and Zaid seasons
Table 4 Multiple regression equation Hishimonus phycitis vs. total weather parameters for Kharif, Rabi and Zaid seasons

Resistance responses of eggplant germplasm against little leaf disease

The experiment was conducted during the year 2021–2022 and aimed to identify resistance sources for little leaf disease in eggplant. None of the tested accessions were found to be free from little leaf disease incidence. Observations were statistically analysed and mean values were tabulated (Table 5). Among the 24 genotypes tested, CBE-SM-105 and CBE-SM-106 exhibited the lowest little leaf disease incidence (8.06%), followed by PLR 1 (9.10%), and CBE-SM-119 (18.24%), which had the highest incidence under natural conditions. In artificially inoculated conditions, a few entries, such as CBE-SM-105, CBE-SM-086, CBE-SM-104, CBE-SM-106, CBE-SM-108, and PLR-1, displayed moderately resistant reactions, while the remaining entries were susceptible and highly susceptible (Table 6).

Table 5 Reaction of eggplant genotypes to little leaf disease under artificial and field conditions
Table 6 Resistance categories of different eggplant genotypes against little leaf disease

AMMI analysis

Genotypes such as G18 (CBE-SM-106), G17 (CBE-SM-105), G16 (CBE-SM-104), G14 (CBE-SM-086), G12 (CBE-SM-084), G1(C0-2), G22 (CBE -SM-110), G3 (CBE-SM-113) and G9 (CBE-SM-119), were plotted away from the x axis in AMMI I biplot which shows low interaction or consistent performance across the environment. However, the lines G18 (CBE-SM-106) and G17 (CBE-SM-105) are highly resistant, since they had low mean values for disease incidence. G9 (CBE-SM-119) and G5 (CBE-SM-115) lines with high mean values were considered as susceptible to little leaf. On the other hand, the genotypes on a horizontal line have similar interaction patterns, such as G12 (CBE-SM-084) and G22 (CBE-SM-110) (Fig. 5).

Fig. 5
figure 5

AMMI1 Biplot display for mean little leaf disease incidence (%) and IPCA 1 scores of the 24 genotypes (G) tested across three environments (E). E1: Kharif, E2: Rabi, E3: Zaid

Genotypes or environments on the right side of the y-axis have higher mean disease incidence, while those placed on the left side have a lower mean. As a result, it is quite evident from the AMMI biplot that, the genotypes G9 (CBE-SM-119), G5 (CBE-SM-115), G4 (CBE-SM-114), G3 (CBE-SM-113), G11 (CBE-SM-083), G22 (CBE-SM-110), G1 (C0-2) and G12 (CBE-SM-084) had a low level of resistance. In contrast, G18 (CBE-SM-106), G17 (CBE-SM-105), G16 (CBE-SM-104), G14 (CBE-SM-086) and G20 (CBE-SM-108) were highly resistant to the little leaf of eggplant. However, in the presence of environmental influence, E3 on the right-hand side are favourable for disease incidence, while E1 and E2 on the opposite side were slightly less favourable.

The AMMI2 biplot derived using PC1 and PC2 shows interaction effect, and the genotypes plotted near the origin are considered to be stable. Thus G14 (CBE-SM-086), G12 (CBE-SM-084), G22 (CBE-SM-110), G3 (CBE-SM-113), G13 (CBE-SM-085), G5 (CBE-SM-115), G19 (CBE-SM-107) and G16 (CBE-SM-104) are stable for resistant to little leaf disease. The lines that are far away from the origin, such as G8 (CBE-SM-118), G6 (CBE-SM-116) and G21 (CBE-SM-109), are less stable for disease incidence. The environments E1, E2 and E3 all had higher level projections from the origin, indicating that all environments are highly discriminative for disease incidence (Fig. 6). Also, E2 observed with low PCA1 scores had low interactions and great stability. However, the high PCA1 scores in E3 and E1 showed high interaction, stability and adaptation to specific environments.

Fig. 6
figure 6

AMMI2 biplot for 24 genotypes (G) and three environments (E) displayed on the first and second principal component axis (IPAC 1 vs. IPAC 2). E1: Kharif, E2: Rabi, E3: Zaid

A significantly larger proportion of genotypes recorded low IPCA1 scores and showed a small interaction effect, which led to the clustering of the genotypes on the biplot. Genotypes G8 (CBE-SM-118) recorded the highest IPCA1 score of 0.6 followed by genotypes viz. G6 (CBE-SM-116), G7 (CBE-SM-117), G20 (CBE-SM-108) and G2 (CBE-SM-112). This indicated that the high interaction effect of genotypes had specific stability and adoption to environments (Tables 7 and 8).

Table 7 Environments
Table 8 AMMI analysis of variance (ANOVA) for G x E interactions

GGE biplot

GGE biplot analysis was used to study the performance of genotypes in each environment. The centre of the concentric circles defines the ideal genotype, while those closer to this point are regarded as near the ideal genotype. In this study, the clusters of genotype present in the centre of the concentric ring show the ideal genotype, whereas G17 (CBE-SM-105), G16 (CBE-SM-104) and G24 (PLR-1) are near the ideal genotype (Fig. 7). In the average environment coordination (AEC) analysis, E3 (Zaid) appeared to be the most discriminating environment when compared with E1 (Kharif) and E2 (Rabi) (Fig. 8). The analysis also showed that all environments were ideal for little leaf disease screening. Genotypes clustered closely to the biplot origin of G17 (CBE-SM-105) were similar in little leaf disease-resistant status and stability. The genotypes scattered from the origin were considered to vary in little leaf disease-resistant status and were found to be less stable. The which-won-where biplot was employed to identify the best genotypes suited to various test environments, where G18 (CBE-SM-106) recorded the best-performing genotype in all test environments (Fig. 9).

Fig. 7
figure 7

Genotype view of GGE biplot analysis of genotypes relative to ideal genotype (centre of the concentric circles)

Fig. 8
figure 8

Environment view of GGE biplot analysis: ranking of environments relative to an ideal test location (represented by the centre of the concentric circles); E1: Kharif, E2: Rabi, E3: Zaid

Fig. 9
figure 9

The which-won-where view of the GGE biplot showing which genotypes performed best in which environments: E1: Kharif, E2: Rabi, E3: Zaid

Discussion

Eggplant or brinjal is a globally important vegetable crop, cultivated year-round in tropical regions. Among the numerous challenges faced by eggplant cultivation, the emergence of eggplant little leaf disease, caused by phytoplasma, poses a significant threat, potentially leading to 100% crop loss [19]. Phytoplasmas are non-helical bacteria devoid of cell walls, known to affect thousands of plant species [15, 21, 23, 26, 38, 42]. They are considered one of the most molecularly enigmatic genera of plant pathogens due to their inability to be cultured in vitro, irregular host distribution, and the challenges associated with serology and electron microscopy-based identification. However, biotechnological interventions, such as nested PCR targeting the 16S rRNA gene, have overcome these limitations, aiding in the identification and classification of phytoplasmas [7, 29]. Phytoplasmas are increasingly recognized as emerging plant pathogens with significant economic impacts, comparable to viruses.

Eggplant little leaf disease was characterized by a conspicuous reduction in leaf size, leading to the gradual shrinkage of newly produced leaves. This was accompanied by a notable reduction in petiole size, resulting in slender, hairless, yellow, and closely appressed leaves with shorter internodes. Axillary buds exhibited excessive growth, developing into stunted branches with tiny leaves, giving the plants a rosette-like appearance. Phyllody, transforming floral structures into vegetative-like structures, was another typical symptom. In severe cases, plants became sterile as flowers remained green, upright, and generally unproductive [24].

These symptoms closely resembled those reported in previous studies of phytoplasma-related diseases. For instance, reduced leaves, virescence, and witches’ broom symptoms in phytoplasma-infected Catharanthus plants in Egypt [32]. Symptoms include stunting, yellow to purple leaves, severe lateral bud proliferation, hypertrophic calyx development, virescence, suppression of ovary and anther production, and petal fusion [33]. Increased branching, reduction in leaf size, and phyllody were also noticed in infected plants. Fruiting was also affected, with malformed fruits and shrivelled seeds in late-stage infections [27].

Molecular detection of phytoplasma in this study utilized nested PCR assays targeting the 16S rRNA gene. The 16S rRNA gene region was amplified with universal primers P1/P6, followed by nested primer pair R16F2n and R16R2, resulting in DNA amplification of approximately 1.5 and 1.2 kb, respectively. These findings align with previous research, including the work of [19], who reported similar results in both plants and vectors, using the same primer pairs. The highly conserved 16S rRNA gene region was employed as the primary molecular tool for phytoplasma identification, genotyping, taxonomic assignment, and group/subgroup classification [22]. All nine phytoplasma samples from plants and one from vectors in this study were identified as ‘Candidatus Phytoplasma trifolii’ related strains, focusing on the conserved region of the 16S rDNA. This study validates the findings of [19], who associated ‘Candidatus Phytoplasma trifolii’ with eggplant little leaf. Also, the phytoplasma associated with eggplant little leaf disease belongs to the 16SrVI clover proliferation phytoplasma group [4].

Phytoplasmas primarily reside in the sieve cells of phloem tissue, inducing various disease symptoms and hormonal imbalances. Insect vectors acquire phytoplasmas from infected plants while feeding on phloem sap. These vectors transport phytoplasmas through their gastrointestinal tracts, hemocoel, and various internal organs, ultimately transmitting them to the plant phloem via their salivary glands [11]. Two leafhopper species, H. phycitis and A. biguttula biguttula, have been implicated in the transmission of little leaf in eggplant. However, their efficiency as vectors may vary. In the eggplant ecosystem of West Bengal [25], reported the presence of both insect vectors, but only H. phycitis was confirmed as a carrier and natural vector of the 16SrVI-D subgroup of phytoplasmas using nested PCR assays, sequence comparison, phylogeny, virtual RFLP analysis, and transmission assays. This study reaffirms the role of H. phycitis as a vector for eggplant little leaf disease [25]. The phloem-sap feeding behaviour of H. phycitis likely contributes to the transmission of eggplant little leaf [19].

A total of 24 genotypes of germplasm were evaluated for the little leaf resistance in which CBE-SM-105 and CBE-SM-106 were moderately resistant also reported earlier [6, 31]. The resistant genotypes could be used as a parent for developing a little leaf-resistant variety to promote cultivation in disease-prone and endemic areas of Tamil Nadu, India. The AMMI analysis model first accounts for the main effect, while the subsequent PCA predicted the interactions [13]. Previous findings also confirmed that AMMI analysis is an efficient tool for screening the resistance level of various genotypes [28, 36]. The additional years of testing the germplasm could confirm a specific test location and environment that will be most representative, discriminating, and ideal. Further research at the molecular level might be useful in identifying the genes and proteins responsible for resistance and susceptibility genotypes. A genotype with higher mean performance and stability across the test locations is said to be an ideal genotype. The which-won-where biplot recorded G18 (CBE-SM-106) as ideal genotypes based on the presence in the centre of the concentric ring. Similarly, a study on genotype G4 (Karuna) was closer to the centre of the concentric ring, followed by G1 (Harasoya) and G5 (Himso) in which Karuna may be chosen as an ideal genotype [30]. Cumulatively, PCA1 and PCA2 contributed to 60% of the total interaction. Thus, the interaction between the 24 genotypes across three test environments, in which the best genotype was predicted by the first two principal compounds. Sufficient research reports have also stated that the genotype interactions were precisely presented by two PCAs [28, 35].

Eggplant cultivation in India often coincides with the cultivation of various other crops, raising concerns about the natural spread of phytoplasmas across different plant species. Efficient vector species facilitate the transmission of phytoplasmas between eggplants and other plants. Thus, it is crucial to comprehensively assess the epidemiological factors contributing to the natural dissemination of phytoplasmas associated with eggplant little leaf disease in India [19]. This study substantiates the presence of a new natural vector and identifies resistant and susceptible genotypes harbouring phytoplasma. Eggplant little leaf disease poses a significant challenge in regions where eggplant is cultivated in India. To effectively combat this disease, a multifaceted approach is needed, including the cultivation of resistant eggplant cultivars, vector management, and regulation of natural hosts [37]. Identifying this vector and the associated conditions for phytoplasma occurrence, along with recognising potential genotypes, as reported in this study, opens new avenues for research. This information holds promise for developing more effective strategies to manage eggplant little leaf disease in India.

Conclusion

In conclusion, this comprehensive study of eggplant little leaf disease, caused by phytoplasma, in Tamil Nadu, India, has provided valuable insights into the symptomatology, molecular detection, phytoplasma characterization, vector species, and evaluation of germplasm resistance. The disease’s symptoms, including leaf reduction, petiole thinning, axillary bud proliferation, and phyllody, closely resemble those observed in phytoplasma-related diseases in other plant species, aiding in its early identification. Molecular detection using nested PCR targeting the 16S rRNA gene confirmed the presence of ‘Candidatus Phytoplasma trifolii’ related strains, facilitating precise diagnosis and taxonomic assignment. The study reaffirmed the role of H. phycitis as a vector for eggplant little leaf disease and identified potential resistant genotypes, providing valuable resources for breeding programs. These findings have significant implications for managing eggplant little leaf disease in India and underscore the importance of a multifaceted approach, including developing resistant cultivars and vector control, to ensure the sustainability of eggplant cultivation in the region.

Methods

Disease survey and collection of samples

Between May to September 2022, leaf samples in triplicate were collected from symptomatic eggplant plants located in nine districts of Tamil Nadu state in India with the farmers’ permission. Samples were collected from 45 to 60 days old eggplant crops where a minimum of three to four samples were from diseased plants when two samples from healthy plants served as the negative control.

DNA extraction and PCR assays

Genomic DNA from leaf samples, as well as for leafhopper species, whiteflies, hadda beetle, and shoot and fruit borer larvae, were extracted using the CTAB method outlined by [1]. Amplification of phytoplasma ribosomal DNA was achieved using the phytoplasma-universal primer (P1/P6) was designed on basis of 16S rDNA sequence which amplified a 1.5 kb product at initial denaturation 94ºC for 2 min, followed by 35 cycles of denaturation at 94ºC for 1 min, annealing at 55ºC for 1 min, extension at 72ºC for 1 min and final extension 72ºC for 10 min [10], Subsequent evaluation using R16F2n/R16R2 in the nested PCR assays by PCR products amplified by the universal primers P1/P6 were diluted 0 to 50 times in sterile deionized water and used as templates which yielded a specific amplicons of 1.2 kb at initial denaturation 94ºC for 2 min, followed by 35 cycles of denaturation at 94ºC for 1 min, annealing at 50ºC for 1 min, extension at 72ºC for 1 min and final extension 72ºC for 10 min [45]. For the nested PCR, the product from the first round of PCR assay was diluted to 1:4 ratio with sterile water and 2 µl of the diluted product was used as template. Subsequently, 5 µl of the nested PCR product were subjected to electrophoresis in a 1.0% (w/v) agarose gel, stained with ethidium bromide, and visualized under a UV transilluminator.

Sequencing, pairwise sequence comparison, and phylogenetic analysis

PCR products of R16F2n/R16R2 were outsourced and sequenced bidirectionally at Eurofins Genomics India Pvt Ltd, Bangalore, India. The BioEdit sequence alignment editor tool was employed to assemble and generate consensus sequences by aligning sequences in both forward and reverse directions [14]. The 16S rRNA gene sequences of the samples were analyzed using the iPhy Classifier program (http://www.plantpathology.ba.ars.usda.gov/cgibin/resource/iphyclassifier.cgi), and the consensus sequences were submitted to GenBank. Additionally, 16S rRNA gene sequences of various phytoplasma groups were obtained from GenBank. Sequence homology percentages with other phytoplasma isolates were noted and subsequently utilized to construct a phylogenetic tree. Sequence alignment was performed using the CLUSTALW2 algorithm [20], and the phylogenetic tree was constructed using the neighbour-joining (NJ) method with MEGA version 11 [46], supported by 1000 bootstrap replications. The phylogenetic tree was rooted using the Acholeplasma laidlawii phytoplasma 16S rDNA sequence (Ac. No. NR074448).

Insect sampling, identification, and population dynamics

The observation of the H. phycitis population involved counting the number of leafhoppers per 10 sweeps. The experiments spanned the Kharif, Rabi, and Zaid seasons of 2021–2022. Weather parameter data, including temperature, rainfall, relative humidity, daylight hours, and wind velocity, for the three seasons were obtained from the Agro Climate Research Centre, TNAU, Coimbatore, Tamil Nadu, India (Tables 9, 10 and 11). The collected insects were preserved in plastic vials containing 70% ethanol and stored at 4ºC for subsequent analysis. Taxonomic identification of the collected insects was carried out by the Department of Entomology, TNAU, Coimbatore, India. In accordance with Kumar et al.‘s protocol [18], transmission assays were conducted using the positive leafhopper samples. Five plants were maintained in both treated and control cages for the experiment.

Table 9 Meteorological data during the crop growth period (Kharif 2021–2022)
Table 10 Meteorological data during the crop growth period (Rabi 2021–2022)
Table 11 Meteorological data during the crop growth period (Zaid 2021–2022)

Germplasm evaluation for resistance to little leaf disease

Germplasm for resistance evaluation comprised 22 accessions and two check varieties, Co2 and PLR1. Of these, three were interspecific hybrid derivatives collected and maintained at the Department of Vegetable Crops, Horticultural College and Research Institute, TNAU, Coimbatore. The evaluation was conducted from February to June 2022 in a randomized block design and with three replications. Observations were recorded for 12 different traits from five randomly selected plants of each replication. The mean percentage of disease incidence was calculated based on the ratio of diseased to total plants across all replications. Five categories were established based on percent disease incidence: (i) immune (0%), (ii) resistant (0.1–10%), (iii) moderately resistant (10.1–20%), (iv) susceptible (20.1–50%), and (v) highly susceptible (> 50%) [8]. Germplasm was also screened for resistance to little leaf disease under artificial conditions using an insect-proof net house. Fifteen plants per line were raised in 25 cm pots, and each plant was inoculated with phytoplasma through side grafting using young shoots from diseased plants [8].

AMMI and genotype by environment (G x E) analyses

AMMI statistical model and computational methods were used in this study as described in [5, 12, 28, 48]. The Analysis of variance (ANOVA) was generated using the computer software program MATMODEL version 3.0 [13]. The ‘genotype and genotype by environment’ (GGE) biplot analysis was utilized to check the stability of genotypes over the Kharif, Rabi and Zaid seasons [34, 47]. The AMMI biplot and PC analysis were carried out using Plant Breeding Tools Version 1.4 to find the susceptibility of the genotypes. As the data contains zero due to non-germination (NG) of seeds, 0.05 has been included for all the individual observations so that the results can be directly used with the biplot graph generated.

Data availability

The 16S rRNA of H. phycitis generated during the current study is available in NCBI GenBank (Accession No: ON870850).

References

  1. Ahrens U, Seemuller E. Detection of DNA of plant pathogenic mycoplasmalike organisms by a polymerase chain reaction that amplifies a sequence of the 16 S rRNA gene. Phytopathology. 1992;82:828–32.

    Article  CAS  Google Scholar 

  2. Mello APOA, Eckstein B, Flores D, Kreyci PF, Bedendo IP. Identification by computer-simulated RFLP of phytoplasmas associated with eggplant giant calyx representative of two subgroups, a lineage of 16SrIII-J and the new subgroup 16SrIII-U. Int J Systemic Evolutionary Microbiol. 2011;61:1454–61.

    Article  Google Scholar 

  3. Anonymous. 2022a. Agriculture data, Agricultural Statistics databases. Organisation of United Nations, Rome, Italyhttp://faostat.fao.org.

  4. Azadvar M, Baranwal VK. Multilocus sequence analysis of phytoplasma associated with brinjal little leaf disease and its detection in Hishimonus phycitis in India. Phytopathogenic Mollicutes. 2012;2:15–21.

    Article  Google Scholar 

  5. Bose LK, Jambhulkar NN, Pande K, Singh ON. Use of AMMI and other stability statistics in the simultaneous selection of rice genotypes for yield and stability under direct-seeded conditions. Chil J Agricultural Res. 2014;74:3–9.

    Article  Google Scholar 

  6. Chakrabarti AK, Choudhury B. Breeding brinjal resistant to little leaf disease. Proceedings of the Indian National Science Academy, 1975, 41: 379–385.

  7. Clarridge IIIJE. Impact of 16S rRNA gene sequence analysis for identification of bacteria on clinical microbiology and infectious diseases. Clin Microbiol Rev. 2004;17:840–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Das AK, Mitrn DK. Screening of brinjal varieties, lines and wild Solanum spp for resistance to little leaf disease. Indian J Agricultural Sci. 1999;69:726–8.

    Google Scholar 

  9. Das AK, Mitra DK, Ahlawat YS. Detection of brinjal little leaf phytoplasma in phloem elements isolated by enzyme treatment. Indian Phytopathol. 2000;53:146–8.

    CAS  Google Scholar 

  10. Deng S, Hiruki C. Amplification of 16S rRNA genes from culturable and nonculturable mollicutes. J Microbiol Methods. 1991;14:53–61.

    Article  CAS  Google Scholar 

  11. Firrao G, Garcia-Chapa M, Marzachì C. Phytoplasmas: genetics, diagnosis and relationships with the plant and insect host. Front Biosci. 2007;12:1353–75.

    Article  CAS  PubMed  Google Scholar 

  12. Gauch HG Jr. A simple protocol for AMMI analysis of yield trials. Crop Sci. 2013;53:1860–9.

    Article  Google Scholar 

  13. Gauch HG. MATMODEL version 3.0: open source software for AMMI and related analyses. Crop and Soil sciences. Ithaca, NY: Cornell University; 2007.

    Google Scholar 

  14. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT, In: Nucleic acids symposium series. Oxford University Press, USA, 1999, pp. 95–98.

  15. Hogenhout S, Oshima K, Ammar ED, Kakizawa S, Kingdom HN, Namba S. Phytoplasmas: bacteria that manipulate plants and insects. Mol Plant Pathol. 2008;9:403–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Kumar J, Gunapati S, Singh SP, Lalit A, Sharma NC, Tuli R. First report of a ‘Candidatus Phytoplasma asteris’ (16SrI group) associated with little leaf disease of Solanum melongena (brinjal) in India. New Disease Rep. 2012;26:p21.

    Article  Google Scholar 

  17. Kumar M. Genetic diversity and natural spread sources of brinjal little leaf phytoplasma. M.Sc. Thesis. Division of Plant Pathology, Indian Agricultural Research Institute, New Delhi, India, 2015 p.103.

  18. Kumar M, Katiyar A, Madhupriya, Rao GP. First report of association of potato virus X and potato virus Y and ‘Candidatus Phytoplasma trifolii’ in brinjal in India. Virus Disease. 2016;27:207–8.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Kumar M, Madhupriya, Rao GP. Molecular characterization, vector identification and sources of phytoplasmas associated with brinjal little leaf disease in India. 3 Biotech. 2017;7:1–11.

    Article  Google Scholar 

  20. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.

    Article  CAS  PubMed  Google Scholar 

  21. Lee IM, Davis RE, Gundersen-Rindal DE. Phytoplasma: phytopathogenic mollicutes. Annu Rev Microbiol. 2000;54:221–55.

    Article  CAS  PubMed  Google Scholar 

  22. Lee IM, Gundersen-Rindal DE, Bertaccini A. Phytoplasma: ecology and genomic diversity. Phytopathology. 1998;88:1359–66.

    Article  CAS  PubMed  Google Scholar 

  23. Li ZN, Zhang L, Bai YB, Liu P, Wu YF. Detection and identification of the elm yellows group phytoplasma associated with Puna chicory flat stem in China. Can J Plant Pathol. 2012;34:34–41.

    Article  Google Scholar 

  24. Maheshwari M, Kumar M, Rao GP. Identification and characterization of phytoplasma associated with brinjal little leaf, phyllody and witches’ broom disease in India. Indian Phytopathol. 2017;70:258–61.

    Article  Google Scholar 

  25. Majumdar S, Das BK. Studies on little leaf of brinjal and morpho taxonomy of the leafhopper species associated from Bengal. J Entomol Zool Stud. 2020;8:514–21.

    Google Scholar 

  26. Marcone C, Neimark H, Ragozzino A, Lauer U, Seemuller E. Chromosome sizes of phytoplasmas composing major phylogenetic groups and subgroups. Phytopathology. 1999;89:805–10.

    Article  CAS  PubMed  Google Scholar 

  27. Mitra DK. In: Raychaudhuri SP, Teakle, editors. Little leaf, a serious disease of eggplant (Solanum melongena), management of Plant diseases caused by fastidious prokaryotes. New Delhi, India: DS, Associated Publishing Co; 1993. pp. 73–8.

    Google Scholar 

  28. Mukherjee AK, Mohapatra NK, Bose LK, Jambhulkar NN, Nayak P. Additive main effects and multiplicative interaction (AMMI) analysis of GxE interactions in rice-blast pathosystem to identify stable resistant genotypes. Afr J Agric Res. 2013;8:5492–507.

    Google Scholar 

  29. Namba S. Molecular and biological properties of phytoplasmas. Proceedings of the Japan Academy, 2019, Series B, 95: 401–418.

  30. Nataraj V, Pandey N, Ramteke R, Verghese P, Reddy R, Onkarappa T, Mehtre SP, Gupta V, Satpute GK, Mohan Y, Shivakumar M, Chandraand S, Rajesh V. GGE biplot analysis of vegetable type soybean genotypes under multi-environmental conditions in India. J Environ Biol. 2021;42:247–53.

    Article  Google Scholar 

  31. Nirmala Devi S, Mathew K, Abdul Vahab M. Evaluation of brinjal for resistance to little leaf disease. South Indian Hortic. 1995;43:120–1.

    Google Scholar 

  32. Omar AF, Emeran AA, Abass JM. Detection of phytoplasma associated with periwinkle virescence in Egypt. Plant Pathol J (Faisalabad). 2008;7:92–7.

    Article  CAS  Google Scholar 

  33. Omar AF, Foissac X. Occurrence and incidence of phytoplasmas of the 16SrII-D subgroup on solanaceous and cucurbit crops in Egypt. Eur J Plant Pathol. 2012;133:353–60.

    Article  CAS  Google Scholar 

  34. Parihar AK, Basandrai AK, Sirari A, Dinakaran D, Singh D, Kannan K, Kushawaha KP, Adinarayan M, Akram M, Latha TKS, Paranidharan V. Assessment of mungbean genotypes for durable resistance to Yellow Mosaic Disease: Genotype× environment interactions. Plant Breeding. 2017;136:94–100.

    Article  CAS  Google Scholar 

  35. Persaud R, Saravanakumar D, Persaud M. Identification of resistant cultivars for sheath blight and use of AMMI models to understand genotype and environment interactions. Plant Dis. 2019;103:2204–11.

    Article  CAS  PubMed  Google Scholar 

  36. Persaud R, Saravanakumar D. Screening for blast resistance in rice using AMMI models to understand G x E interaction in Guyana. Phytoparasitica. 2018;46:551–68.

    Article  CAS  Google Scholar 

  37. Rao G, Mall S, Raj S, Snehi S. Phytoplasma diseases affecting various plant species in India. Acta Phytopathologica et Entomologica Hungarica. 2011;46:59–99.

    Article  CAS  Google Scholar 

  38. Rao GP. Our understanding about phytoplasma research scenario in India. Indian Phytopathol. 2021;74:371–401.

    Article  Google Scholar 

  39. Rao GP, Kumar M. World status of phytoplasma diseases associated with eggplant. Crop Prot. 2017;96:22–9.

    Article  Google Scholar 

  40. Rodrigues PC, Malosetti M, Gauch HG Jr, Eeuwijk FA. A weighted AMMI algorithm to study genotype by environment interaction and QTL by environment interaction. Crop Sci. 2014;54:1555–70.

    Article  Google Scholar 

  41. Schneider B, Seemueller E, Smart CD, Kirkpatrick BC. Phylogenetic classification of plant pathogenic mycoplasma-like organisms or phytoplasmas. Mol Diagn Procedures Mycoplasmology. 1995;1:369–80.

    Article  CAS  Google Scholar 

  42. Seemüller E, Garnier M, Schneider B. Mycoplasmas of plants and insects. Molecular Biology and pathogenicity of Mycoplasmas. Boston, MA: Springer; 2002. pp. 91–115.

    Book  Google Scholar 

  43. Shukla V, Naik LB. Agro-techniques of solanaceous vegetables. Adv Hortic. 1993;5:365.

    Google Scholar 

  44. Singh N, Rao GP, Upadhayaya PP. Candidatus Phytoplasma Trifolii’ associated with little leaf and witches’ broom disease of Datura stramonium L. in India. Phytopathogenic Mollicutes. 2012;2:69–71.

    Article  CAS  Google Scholar 

  45. Smart CD, Schneider B, Blomquist CL, Guerra LJ, Harrison NA, Ahrens U, Lorenz KH, Seemüller E, Kirkpatrick BC. Phytoplasma-specific PCR primers based on sequences of the 16S-23S rRNA spacer region. Appl Environ Microbiol. 1996;62:2988–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Tamura K, Stecher G, Kumar S. MEGA11: molecular evolutionary genetics analysis version 11. Mol Biol Evol. 2021;38:3022–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Yan W. GGE biplot - A Windows application for graphical analysis of multi environment trial data and other types of two-way data. Agron J. 2001;93:1111–8.

    Article  Google Scholar 

  48. Gauch HG, Zobel RW. AMMI analysis of yield trails. Genotype by Environment Interaction. Boca Raton, Florida: CRC Publication; 1996. pp. 88–122.

    Book  Google Scholar 

Download references

Acknowledgements

We express our gratitude to the Department of Plant Pathology, Centre for Plant Protection Studies, Tamil Nadu Agricultural University, Tamil Nadu and the Centre for Crop Health, University of Southern Queensland, Toowoomba, Australia for providing the facilities and data analytical support during the research process.

Funding

We acknowledge FMC India Private Limited, Mumbai and National Agricultural Higher Education Project - Institutional Development Plan, Indian Council of Agricultural Research (ICAR), India and University of Southern Queensland’s Strategic Research Scheme for funding support.

Author information

Authors and Affiliations

Authors

Contributions

MK, SP, GPY and TE conceived the study and designed the experiments. MK, IJ, GPY, SBJ, SPD, NMB and SM conducted the survey, collection of samples and PCR analysis. TE conducted insect identification and population dynamics. BAJJ for germplasm evolution studies. NMB for molecular and sequence analysis. SM for meteorological data analysis. HSN, DLA MK, and SP for statistical analysis. MK, GPY, SBJ and SPD wrote the initial draft. SP revised the draft and all authors provided comments.

Corresponding authors

Correspondence to Muthusamy Karthikeyan or Sambasivam Periyannan.

Ethics declarations

Ethics approval and consent to participate

Experiments were conducted in accordance with the relevant institutional guidelines. Doesn’t involve any experiments with animals.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Karthikeyan, M., Yogiraj, G.P., Elaiyabharathi, T. et al. Comprehensive analysis of little leaf disease incidence and resistance in eggplant. BMC Plant Biol 24, 576 (2024). https://doi.org/10.1186/s12870-024-05257-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05257-4

Keywords