Phytophthora Infection Reorients the Composition of Rhizospheric Microbial Assembly in Khasi Mandarin ( Citrus reticulata Blanco)

: Phytophthora gummosis, foot rot, and root rot are considered major challenges to the citrus industry worldwide. Little is known about the Phytophthora –microbiome interaction, despite several studies demonstrating changes in the microbial composition of the rhizosphere following challenges by a pathogen. In the present study, we studied the microbial diversity and community structure in healthy rhizospheres and Phytophthora -infected rhizospheres of Khasi mandarin ( Citrus reticulata Blanco), a commercial cultivar extensively grown in the northeast of India. An exploratory study was conducted to identify Phytophthora -infected orchards of Khasi mandarin, and the isolated pathogen was confirmed as P. nicotianae based on its morpho-cultural and molecular characteristics coupled with pathogenicity tests. This study on culturable microbes established the dominance of Trichoderma spp. in the healthy rhizosphere, while the diseased rhizosphere showed the presence of Fusarium spp. A metagenomic study further revealed that the rhizospheres of Phytophthora -infected plants were dominated by species such as Bacteroidia spp., Patescibacteria spp., and Pythium spp., while the healthy Khasi mandarin rhizospheres had a more diverse community predominantly represented by Trichoderma , Penicillium , Linnemannia , Mortierella , Talaromyces , Saitozyma , Bacteroidetes , Pseudomonas , Cytophagia , Cyanobacteria , Bacteroidia , Sphingobacteriia , Burkholderia , Bacillus , and Bradyrhizobium . Terrabacteria and FCB ( Fibrobacterota , Chlorobiota , and Bacteroidota groups) were found to exist in higher relative abundance in disease-free soils than in Phytophthora-infected soils, while phylum Proteobacteria showed identical relative abundance in all soil types. The phyla represented by Pseudomonas , Flavobacteriia , Candidatus , Mycobacterium , Rhizobium , Mesorhizobium , Sphingomonas , and Cytophagia were the most common bacterial phyla in all soil samples, but healthy soil exhibited a greater abundance of Bacteroidetes , Pseudomonas, Cytophagia, Cyanobacteria, Bacteroidia, Sphingobac-teriia, Burkholderia, Bacillus, and Bradyrhizobium . Our study suggests that the presence of Phytophthora spp. in the rhizosphere alters microbial community structure, having potentially strong implications for plant health and productivity. These rhizosphere microbiome-derived citrus responses shed light on exploring effective management strategies for Phytophthora gummosis disease ailing Khasi mandarin.


Introduction
The rhizosphere is "the narrow region of soil that is directly influenced by root secretions and associated soil microorganisms" [1], a critical interface between plant roots and soil microorganisms playing a crucial role in plant growth and health.It is a dynamic and complex environment influenced by various factors such as plant genotype, soil type, and environmental conditions.The concept of the rhizosphere was first introduced by Hiltner in 1904 [2].The rhizosphere microbiome, comprising diverse microbial communities, is reported to have strong implications for plant health, disease resistance, and the ecosystem [3,4].Microorganisms in the rhizosphere interact with plants through various chemical signals and mechanisms, altering the soil's physical and chemical properties, and ultimately influencing plant productivity.Recent studies revealed the role of the rhizosphere in shaping plant-associated microbial communities and their interactions with pathogens [5,6].For example, a study by Yuan et al. [7] showed that the rhizosphere microbiome of maize plants influenced the colonization and virulence of fungal pathogen Fusarium verticillioides by modulating the expression of plant genes involved in jasmonic acid (JA) and salicylic acid (SA) signalling.Similarly, Berg et al. [8] found that the rhizosphere microbiome of Arabidopsis thaliana promoted the growth and biocontrol activities of fungal endophyte Serendipita indica, and enhanced the plant's resistance against the foliar pathogen Pseudomonas syringae.Consequently, any alteration in the microbial makeup of the rhizosphere has the potential to influence its ecological roles, relating to the physiological condition of the plant, and ultimately, the plant's productivity.
Root rot or gummosis is a devastating disease of citrus plants, caused by Phytophthora spp., a group of soilborne pathogens [9] with sensitivity against all major fungicides [10].This disease is characterized by the formation of gum pockets on the bark, trunk, and branches, eventually leading to the dieback and decline of the tree.The pathogen infects the roots, trunk, and branches of the tree, and the infection is often associated with wounds or injuries caused by mechanical damage or adverse environmental conditions [11].In recent years, research has focused on the effects of Phytophthora infection on the composition of rhizospheric microbial communities, playing an important role in plant health and productivity.A study by Yang et al. [12] examined the effects of Phytophthora infection on the rhizospheric microbial community composition of citrus plants.
Khasi mandarin (Citrus reticulata Blanco) is an important citrus fruit crop grown in the northeastern region of India, particularly in the states of Meghalaya, Assam, and Arunachal Pradesh [13], having significant economic and cultural value in the region.However, the production of Khasi mandarin is constrained by various biotic and abiotic factors, including pest and disease pressures.Phytophthora spp., a soil-borne pathogen, is a major risk to Khasi mandarin, causing root rot and gummosis in citrus plants.The disease leads to significant yield losses and poses a major constraint on citrus production in the region.For instance, a report from Singh et al. [14] stated that Khasi mandarin is the most preferred citrus fruit in the eastern Himalayan region, the most notable northeastern region of India.However, it was also notable that the crop is vulnerable to various diseases, including Phytophthora root rot, frequently leading to significant losses in fruit yield and quality.Similarly, Khasi mandarin is reported to be susceptible to a range of pests and diseases, including citrus psyllid, citrus greening, foot rot, etc., which could further adversely affect the productivity and quality of the crop [15,16].
Despite these challenges, India has made significant progress towards the development of management strategies for Phytophthora-induced foot and root rot, comprising the development of resistant varieties and the use of integrated disease management approaches to ensure the sustainability of Khasi mandarin production in the region.For example, Deka et al. (2018) [17], in their studies, assessed three modules, the bio-intensive IPM module (BIPM), IPM module, and farmers' practice (FP), based on previous field studies, which showed IPM displaying significant promise with regard to yielding higher net returns, followed by BIPM.However, the sustainable field management of Phytophthora is still a formidable challenge.Therefore, culturable and metagenomic-based studies of the rhizosphere of citrus are crucial in understanding citrus' microbial diversity and its functional contribution towards the health and productivity of citrus crops.The objectives of such studies were to identify microbial diversity and its functional roles, assess the effect of soil management practices on the microbial community, and evaluate the interaction between the plant and the microbial community.These studies will pave the way towards developing more effective strategies for the microbial antagonist-mediated management of the Phytophthora disease complex, a problem researchers are as yet clueless about, except for the chemical control strategy.

Rhizospheric Soil Sampling and Isolation of Pathogen
We conducted sampling activity in a 10-year-old Khasi mandarin orchard (showing signs of branch dieback due to root rot) located in Khaman Pathar (Sivasagar), Assam, India, in December 2021.We further confirmed the presence of root rot by visually examining the youngest feeder roots.Briefly, 100 g of soil samples was collected from each of the four sides of the tree canopy using a 25 mm soil auger from inside of the canopy skirt, ensuring the soil contained citrus major feeder roots.The samples were collected with a shovel that was washed with 70% ethanol between each sampling point.Each sample, including rhizospheric soil and feeder roots, individually packed in zip-locked polyethene bags within an ice pack, properly labelled with sample descriptions, were finally brought to the laboratory, in Department of Plant Pathology, Assam Agricultural University, Jorhat (India).Subsequent dilutions were prepared using 1 g of bulked rhizosphere soil and 99 mL of sterile distilled water, and homogenized via vigorous shaking.Dilutions of 1:10 and 1:100 were then streaked onto Petri dishes with Corn Meal Agar-PARPH medium (Cat.60786; Sigma, St. Louis, MO, USA) in triplicates.Plates were incubated at 24 • C for three days.The purified isolates were examined for their cultural traits, such as colony colour, topography, margin, opacity, and mycelial growth pattern on the PDA medium, five days after inoculation.Morphological features such as sporangial morphology, including shape and size, were studied through microscopic observation (ZEISS Axiolab 5, Oberkochen, Germany).

Molecular Characterization
DNA extraction was carried out from the isolated Phytophthora spp.following the CTAB protocol originally described by Prabha et al. in 2013 [18] and further amplified using ITS-1 and ITS-4 primers.The DNA was checked for quality and quantity using agarose gel electrophoresis.Amplification was performed with an initial denaturation step at 94 • C for 4 min followed by 30 cycles of denaturation, annealing, and extension.The amplified DNA was visualized using gel electrophoresis.The sequence data were analysed using the Bio-Edit sequence alignment editor and searched against the NCBI database using nucleotide BLAST.A phylogenetic tree was constructed using MEGA version 11 with the maximum likelihood method and 500 bootstrap replicates.

Pathogenicity Test
The pathogenic potential of identified fungal isolates was assessed as previously described by Das et al. (2017) [19] by taking on detached leaves and excised stems using two methods, the mycelial bit inoculation technique (MBIT) and the pin prick method [20].Under the MBIT, healthy Khasi mandarin leaves were surface-sterilized and then injured with a sterilized pin.Mycelial discs of 7-day-old Phytophthora spp.were cut and placed on the injured leaves, and kept in a controlled environment for up to 120 h, while under the excised inoculation method, small stems were surface-sterilized and a block of Phytophthora spp. was placed on the exposed portion of the stems kept in a controlled environment for up to 10 days.For comparisons, control leaves and stems were inoculated with sterile water agar discs.Necrotic progression on the leaves and stems was visually monitored every day.Five repetitions were carried out for better visualization.

Isolation and Identification of Culturable Fungi from Healthy and Diseased Plant Rhizosphere
In order to retrieve the potentially pathogenic and non-pathogenic fungi associated with the rhizosphere, fungal isolations were carried out from the soil samples collected previously from both healthy and diseased Khasi mandarin rhizospheres.For the soil samples, upon serial dilution, 0.5 mL of the soil water sample was added to a Petri plate containing nutrient agar medium, and a dilution factor of 1:100,000 was used to isolate bacteria from the soil samples.Similarly, 1 ml of the 1:1000 dilution soil water sample was added to Petri dishes containing sterilized PDA medium for fungi and incubated upside-down at 25-30 • C in the BOD incubator.Bacterial colony-forming units (CFU) were estimated by counting the number of colonies using the method described by Chowdhry and Varshney (2000) [21].Fungal morphological observations of mycelia, conidiophores, and conidia were taken using Zeiss Axiolab 5. Recovered rhizospheric microbes from healthy and diseased plant rhizospheres were characterized based on their morphological and reproductive structures, consulting various monographs [21][22][23][24][25]. Purified bacteria were maintained in NA slants for further experiments.
Relative abundance (RA): This method calculates the proportion of individuals in a sample belonging to a particular species or group compared with the total number of individuals in the sample.RA is expressed as a percentage, and represents the contribution of a species to the overall abundance of the community; it is calculated as follows: RA (%) = Ci/Ct × 100, where Ci is the Total number of colonies of individual species and Ct = Total number of colonies of all species.

Bioinformatics Analysis 2.5.1. Metagenomic Analysis
Metagenomic DNA was extracted from healthy and diseased plants' rhizospheric soil samples collected from Khaman Pathar using a DNA concentration kit, and purity was determined using NanoDrop Spectrophotometer (Thermo Fischer Scientific, Waltham, MA, USA).Extracted DNA was sent for Metagenome Shotgun Sequencing at Centyle Biotech Private Limited, Rudrapur, Uttarakhand, India.After DNA fragmentation, Illumina sequencing adapters were added, and products were amplified using PCR, during which unique indices were added to each sample.After library amplification, libraries were sequenced using the NovaSeq 6000 system (Illumina, San Diego, CA, USA) following standard methods provided by the manufacturer for pair-end sequencing (2 × 150 bp).Shotgun metagenomic sequencing was performed using the NOVASEQ6000 sequencer (Illumina) in accordance with the standard protocol (https://www.illumina.com/,accessed on 26 February 2022).
Reads' adapter sequences were filtered and trimmed using Fastp v0.23.2.Reads were treated to have an average Phred score of 20; we removed bases with low complexity, and trimmed poly-tails if there were any.After quality control, the healthy sequences had 39,951,198 metagenomic reads whereas the infected sequences had 35,301,808 metagenomic reads.These were assembled into contigs (metagenomes were assembled individually) using the software MEGAHIT, version 1.2.8 [26], with kmer values from 21 to 127.
We used a protein space-based annotator called Kaiju [27], which essentially performs a 6-frame translation of the datasets and uses a non-redundant taxonomy database comprising bacteria, archaea, viruses, and micro-eukaryotes like fungi and dikarya.To this end, we used a modified Greedy classification approach to classify and assign taxonomies to as many contigs/reads as possible, by using a minimal match of 7 amino acids, and a match score of 50, with an e value of 0.01 [27].

Quantitative Analysis Shannon-Weiner Diversity Index
This method measures the diversity of a community by taking into account both the number of species present and their relative abundances.The Shannon-Weiner index [28] considers the information content of each species and provides a measure of the uncertainty associated with predicting the identity of a randomly selected individual from the community.The higher the Shannon index value, the greater the diversity of the community.
where H and ρi represent the Shannon index and relative abundance, respectively.

Shannon-Weiner Evenness Index
This method measures how evenly individuals are distributed among the species in a community [29].It considers the relative abundances of species and the total number of species present.The Shannon-Weiner evenness index ranges from 0 to 1, with 1 indicating a perfectly even distribution of individuals among species.
where E, i, and m represent species evenness, relative abundance, and total number of species, respectively.

Identification and Characterization of Pathogen
We conducted a thorough survey of the productive Khasi mandarin orchards (Figure 1).where H and ρi represent the Shannon index and relative abundance, respectively.

Shannon-Weiner Evenness Index
This method measures how evenly individuals are distributed among the species in a community [29].It considers the relative abundances of species and the total number of species present.The Shannon-Weiner evenness index ranges from 0 to 1, with 1 indicating a perfectly even distribution of individuals among species.
where E, i, and m represent species evenness, relative abundance, and total number of species, respectively.

Identification and Characterization of Pathogen
We conducted a thorough survey of the productive Khasi mandarin orchards (Figure 1).Phytophthora-induced gummosis disease was identified and confirmed based on field observations of external symptoms including sticky brown fluid and charred bark on the plant's trunk, brown rot, or sores on the skin of the diseased fruits (Figure 2).Phytophthora-induced gummosis disease was identified and confirmed based on field observations of external symptoms including sticky brown fluid and charred bark on the plant's trunk, brown rot, or sores on the skin of the diseased fruits (Figure 2).

Morpho-Cultural Characteristics and Molecular Identification of Pathogen
One isolate associated with Phytophthora gummosis was recovered and purified in Corn Meal Agar-PARPH media.The pure cultures showed distinctive fluffy to cottony growth patterns after being subjected to 28 °C incubation on Potato Dextrose Agar media (Figures 3 and 4).Careful observations made under a ZESIS Axiolab 5 microscope revealed distinctive ovoid-to-obpyriform-shaped sporangia (L × B: 44.8-45.3× 26.6-27.2microns) that were notably papillate, were occasionally spherical, and often had a distorted shape (Table 1).The morpho-cultural characteristics indicated that the pathogen was Phytophthora spp.

Morpho-Cultural Characteristics and Molecular Identification of Pathogen
One isolate associated with Phytophthora gummosis was recovered and purified in Corn Meal Agar-PARPH media.The pure cultures showed distinctive fluffy to cottony growth patterns after being subjected to 28 • C incubation on Potato Dextrose Agar media (Figures 3 and 4).Careful observations made under a ZESIS Axiolab 5 microscope revealed distinctive ovoid-to-obpyriform-shaped sporangia (L × B: 44.8-45.3× 26.6-27.2microns) that were notably papillate, were occasionally spherical, and often had a distorted shape (Table 1).The morpho-cultural characteristics indicated that the pathogen was Phytophthora spp.

Morpho-Cultural Characteristics and Molecular Identification of Pathogen
One isolate associated with Phytophthora gummosis was recovered and purified in Corn Meal Agar-PARPH media.The pure cultures showed distinctive fluffy to cottony growth patterns after being subjected to 28 °C incubation on Potato Dextrose Agar media (Figures 3 and 4).Careful observations made under a ZESIS Axiolab 5 microscope revealed distinctive ovoid-to-obpyriform-shaped sporangia (L × B: 44.8-45.3× 26.6-27.2microns) that were notably papillate, were occasionally spherical, and often had a distorted shape (Table 1).The morpho-cultural characteristics indicated that the pathogen was Phytophthora spp.Three PCR products (Takara PCR Master mix) obtained after PCR amplification, labelled as P1, P2, P3, showed a band size of approx.943 bp aligned with a 100 bp ladder (Takara).The PCR products were sequenced by Eurofins Ltd. (Bothell, WA, USA), and the sequences were used to construct a phylogenetic tree, while 33 closely related species showed the nearest similarity with the Phytophthora nicotianae strains of West Bengal (P2) and Italy (P1 and P3) (Figure 5).Three PCR products (Takara PCR Master mix) obtained after PCR amplification, labelled as P1, P2, P3, showed a band size of approx.943 bp aligned with a 100 bp ladder (Takara).The PCR products were sequenced by Eurofins Ltd. (Bothell, WA, USA), and the sequences were used to construct a phylogenetic tree, while 33 closely related species showed the nearest similarity with the Phytophthora nicotianae strains of West Bengal (P2) and Italy (P1 and P3) (Figure 5).

In Vitro Studies on Pathogenicity on Different Plant Parts
The identified isolate, P. nicotianae, was tested for its pathogenicity using different plant parts: the leaf, excised shoot, and root (Figure 6).Though the symptoms varied, all three plant parts were observed to be affected with the pathogen, displaying its aggressive behavioural polymorphism with regard to the host-pathogen interaction.

Identification of Microbes from Diseased and Healthy Rhizospheres
Colonies with varied morphology and a dominant growth pattern were isolated from the microbial-rich rhizosphere of Khasi mandarin.In total, 18 fungal isolates and 10 bacterial isolates were isolated from the diseased plant's rhizosphere, showing typical symptoms of Phytophthora infection.Similarly, from the healthy rhizosphere, 16 isolates of fungi and 28 isolates of bacteria with distinct morphology were isolated for further characterization.

Morphological Characterization and Relative Abundance of Isolated Fungal Species
The results of the determination of the morphological and reproductive characteristics of the fungal species plus a microscopic study, besides the observations of [22][23][24][25], collectively showed that 16 isolates from the healthy rhizosphere soil belonged to eight different fungal genera, including Fusarium, Penicillium, Aspergillus, Acremonium, Mucor, Geotrichum, Culvularia, and Trichoderma (Figure 7).In contrast, the as many 18 fungal isolates from the diseased rhizosphere were represented by seven different fungal genera: Penicillium, Aspergillus, Geotrichum, Culvularia, Fusarium, and Mucor (Figure 8).
Comparing the relative abundance of different culture-dependent fungal isolates between the health and diseased rhizosphere (Table 2) showed a complete absence of Trichoderma spp. on one hand and the increasing dominance of Aspergillus spp.and Fusarium spp. in the diseased rhizosphere on the other; however, a distinct reduction in the abundance of Penicillium in the healthy rhizosphere was observed compared with that in the diseased rhizosphere.These observations lend strong support in favour of structural and compositional (as well as possibly functional) changes in the rhizospheres, depending upon their status, i.e., healthy or disease-infected.

In Vitro Studies on Pathogenicity on Different Plant Parts
The identified isolate, P. nicotianae, was tested for its pathogenicity using different plant parts: the leaf, excised shoot, and root (Figure 6).Though the symptoms varied, all three plant parts were observed to be affected with the pathogen, displaying its aggressive behavioural polymorphism with regard to the host-pathogen interaction.The results of the determination of the morphological and reproductive characteristics of the fungal species plus a microscopic study, besides the observations of [22][23][24][25], collectively showed that 16 isolates from the healthy rhizosphere soil belonged to eight different fungal genera, including Fusarium, Penicillium, Aspergillus, Acremonium, Mucor, Geotrichum, Culvularia, and Trichoderma (Figure 7).In contrast, the as many 18 fungal isolates from the diseased rhizosphere were represented by seven different fungal genera: Penicillium, Aspergillus, Geotrichum, Culvularia, Fusarium, and Mucor (Figure 8).Comparing the relative abundance of different culture-dependent fungal isolates between the health and diseased rhizosphere (Table 2) showed a complete absence of Trichoderma spp. on one hand and the increasing dominance of Aspergillus spp.and Fusarium spp. in the diseased rhizosphere on the other; however, a distinct reduction in the abundance of Penicillium in the healthy rhizosphere was observed compared with that in the diseased rhizosphere.These observations lend strong support in favour of structural and compositional (as well as possibly functional) changes in the rhizospheres, depending upon their status, i.e., healthy or disease-infected.

Diversity of Rhizospheric Bacterial and Fungal Communities
Soil samples were further evaluated for CFUs (colony forming units) per gram and showed that the CFU counts of microbes varied significantly between healthy and diseased rhizosphere soil.The average CFU values for the bacterial population ranged from (2.3-4.6 × 10 6 ) to (4.8-6.1 × 10 6 ) in healthy and diseased rhizospheric soil, respectively, while the fungal counts ranged between 2.9-4.1 × 10 4 and 3.6-4.8× 10 4 in healthy and diseased rhizosphere soil, respectively (Table 3).The microbial diversity in the rhizosphere soils, determined using the Shannon-Weiner index, showed comparatively higher diversity values in the disease-infected rhizosphere, ranging from 1.06 to 1.49 against values of 0.97 to 1.18 in the healthy rhizosphere.Additionally, for the mean species healthy soil samples, these values varied between 0.53 and 0.71.These showed that the diversity index in the diseased soil was far higher compared with that in the healthy soil (Table 4).

Metagenomics Analysis
From the metagenomic sequencing, raw reads of 6.15 Gbp (DAM_H) and 5.82 Gbp (DAM_D) were obtained from healthy and infected soil after quality control, and totals of 6.01 Gbp (DAM_H) and 5.31 Gbp (DAM_D) were retained with an average of 16.193 Gbp for each sample.The clean reads were assembled into 14,645 contigs (>1000 bp) and 2137 contigs (>1000 bp) for the healthy and infected rizosphere, respectively.The detailed sequencing data and metagenome assembly statistics in each sample are listed (Table 5).Metagenomics-related abundance analysis showed the occurrence of nine prominent phyla, including Proteobacteria, Terrabacteria, Actinobacteria, Sphingobacteria, Flavobacteria, Patescibacteria, Cyanobacteria, Planctobacteria, and Acidobacteria, regardless of the nature of the rhizosphere, whether derived from healthy or disease-infected trees, although the proportionate distribution of these phyla differed between the two types of rhizosphere samples (Figure 9).The diseased rhizosphere soil (DAM_D) had a higher number of bacterial communities represented by seven phyla, while the healthy rhizosphere soil samples (DAM_H) showed a far lower number of bacterial phyla comparatively, marked by the absence of Flavobacteria, Patescibacteria, PVC, and Acidobacteria, providing some strong clues with which to exploit upon such biochemical changes towards microbially sustainable management options.Interestingly, the comparative study of DAM_D (with an abundance range of 0.409-33.118%and a mean value of 8.334% with a standard deviation of 11.196%) versus DAM_H (with an abundance range of 0.000-33.06%and a mean value of 14.076 with a standard deviation of 14.076%) further revealed a marginal difference when parameters such as abundance range and its mean value were coupled with the standard deviation.These observations are indicative of DAM_H displaying a comparatively higher mean abundance variability of bacterial phyla.
At the genus level, the soils of Khaman Pathar had higher occurrences of Bacteroidetes and Pseudomonas.Other notable genera present in the diseased soils included Streptomyces, Corynebacteria, Candidatus, Patescibacteria, Mycobacterium, Rhizobium, Mesorhizobium, Sphingomonas, Bacteroidia, Propionibacteria, Acinetobacter, Cytophagia, Cyanobacteria, Bacteroidia, and Sphingobacter (Figure 10).The diseased rhizosphere soil (DAM_D) had a higher number of bacterial communities represented by seven phyla, while the healthy rhizosphere soil samples (DAM_H) showed a far lower number of bacterial phyla comparatively, marked by the absence of Flavobacteria, Patescibacteria, PVC, and Acidobacteria, providing some strong clues with which to exploit upon such biochemical changes towards microbially sustainable management options.Interestingly, the comparative study of DAM_D (with an abundance range of 0.409-33.118%and a mean value of 8.334% with a standard deviation of 11.196%) versus DAM_H (with an abundance range of 0.000-33.06%and a mean value of 14.076 with a standard deviation of 14.076%) further revealed a marginal difference when parameters such as abundance range and its mean value were coupled with the standard deviation.These observations are indicative of DAM_H displaying a comparatively higher mean abundance variability of bacterial phyla.
At the genus level, the soils of Khaman Pathar had higher occurrences of Bacteroidetes and Pseudomonas.Other notable genera present in the diseased soils included Streptomyces, Corynebacteria, Candidatus, Patescibacteria, Mycobacterium, Rhizobium, Mesorhizobium, Sphingomonas, Bacteroidia, Propionibacteria, Acinetobacter, Cytophagia, Cyanobacteria, Bacteroidia, and Sphingobacter (Figure 10).The minimum abundance value for DAM_D is 0.000%, and the maximum abundance value is 14.467%, with a mean of 2.145% and a standard deviation of 3.185%.For DAM_H the minimum value is 0.000%, and the maximum value is 22.276%, with a mean of 2.595% and a standard deviation of 4.937%.The relative abundance in diseased soil (DAM_D) is lower than the relative abundance of other microorganisms in healthy soil (DAM_H) (Figure 11).The minimum abundance value for DAM_D is 0.000%, and the maximum abundance value is 14.467%, with a mean of 2.145% and a standard deviation of 3.185%.For DAM_H, the minimum value is 0.000%, and the maximum value is 22.276%, with a mean of 2.595% and a standard deviation of 4.937%.The relative abundance in diseased soil (DAM_D) is lower than the relative abundance of other microorganisms in healthy soil (DAM_H).(Figure 11).

Diversity and Abundance of Fungi in Soil Samples
Metagenomic analysis of the data based on the relative abundance at the phylum level revealed that Ascomycota, Zygomycota, and Basidiomycota were the dominant phyla (Figure 12).

Diversity and Abundance of Fungi in Soil Samples
Metagenomic analysis of the data based on the relative abundance at the phylum level revealed that Ascomycota, Zygomycota, and Basidiomycota were the dominant phyla (Figure 12).

Diversity and Abundance of Fungi in Soil Samples
Metagenomic analysis of the data based on the relative abundance at the phylum level revealed that Ascomycota, Zygomycota, and Basidiomycota were the dominant phyla (Figure 12).Over 60% of fungal communities in general belonged to the Ascomycota phylum in all soil samples.The mean value for DAM_D is 8.334, and the standard deviation is 11.196.It is quite noticeable that the mean value for DAM_H is higher than that of DAM_D, with a value of 9.512, and that the former had a higher standard deviation of 14.076, showing a dominant pattern.(Figure 10).Similar with that of the phylum, genusspecific diversity further showed the higher abundance of the Phytophthora genus in the diseased rhizosphere (DAM_D) compared with that in the healthy rhizosphere (DAM_H).The DAM_D samples showed a 0.86% abundance of Phytophthora, while DAM_H samples showed an abundance of 0.30%, thereby showing a net reduction in the presence of the pathogen of 0.56% (Figure 10).
It is noteworthy that the metagenome of the diseased soil samples had a much higher content of Fusarium (13.38%) compared with that of the healthy soil samples.The healthy soil sample from Khaman Pathar also had a particularly high abundance of Penicillium (21.02%).The mean for DAM_D is 0.027, and the standard deviation is 0.025.The mean for DAM_H is higher than that of DAM_D, with a value of 0.053, and the former also has a higher standard deviation of 0.034.

Discussion
Khasi mandarin is an important fruit crop grown in tropical and subtropical parts of Arunachal Pradesh, Meghalaya, and Assam, India, but it is often challenged by fungal pathogens such as Phytophthora spp., the major cause of gummosis and root rot in growing plantations [19].Due to the polycyclic nature of the pathogen and the perennial nature of the crop, chemical management strategies could not determine the root cause of the disease [7,8].However, it was found that some soils of the rhizosphere showed suppressive effect on Phytophthora spp., suggesting that differences in the microbial community structure of healthy vs. pathogen-rich soils may play a pivotal role in finding the microbial source of disease development as well as disease suppressiveness [29]; it may unlock potential strategies for controlling Phytophthora diseases in the crop using beneficial microbes with antagonistic traits.
The major symptoms of Phytophthora spp.-induced disease usually manifest as the oozing of gums followed by the gradual blackening of the trunk [30].Among the nine surveyed orchards, Khaman Pathar and Motapung (Tinsukia, Assam, India), and Kamalabari (Sivasagar, Assam, India) showed typical symptoms like those mentioned before in these three locations, where we were able to isolate three possible Phytophthora spp., giving them the codes P1, P2, and P3.These kinds of characteristic symptoms of Phytophthora infection were recorded independently in different time periods by Mekonen et al. [29], Das et al. [30], and El-Guilli et al. [31] for sweet orange, Kinnow mandarin, and sour orange, respectively [29].
Morpho-cultural characteristics of these three isolates, P1, P2 and P3, showed a topography of a flat to fluffy colony with a smooth margin, without any pigmentation.These cultural features were also reported earlier by the authors of [32,33], who identified that these characteristics of Phytophthora spp.Were associated with Khasi mandarin as well.Microscopic features revealed that sporangia that were papillate, ovoid, elongated, occasionally spherical, and often distorted in shape, were detected in all the isolates, as reported through earlier studies as well [30,31].
These morphological characteristics were earlier associated with Phytophthora nicotianae [34,35].Phytophthora causing gummosis, and root and foot rot diseases were recorded in citrus orchards of Maharashtra and Punjab due to P. nicotianae.[19].The genomic DNA from three Phytophthora isolates, P1, P2, and P3, were amplified with ITS1 and ITS4 primers.Agarose gel amplification showed amplification products measuring ~943 bp in size using a universal primer (Figure 2).After amplification, PCR products were outsourced for sequencing, and the results were checked against the GenBank reference sequence database at the National Centre for Biotechnology Information (NCBI, www.ncbi.nlm.nih.gov,accessed on 15 December 2023).Phylogenetic analysis revealed that all the isolates displayed the closest similarity with Phytophthora nicotianae.(Figures 3-5).
Pathogenicity tests showed isolated Phytophthora spp.causing symptoms of brown water-soaked discoloration and white-to-whitish mycelial growth on detached leaves, excised stems, and feeder roots (Figure 6), depending upon the plant parts used.The appearance of these symptoms over a period of 2-7 days following inoculation indicated the ability of the identified pathogen to produce disease symptoms in Khasi mandarin plants.The control leaves showed no discoloration or mycelial growth after inoculation, confirming the specificity of the symptoms against the Phytophthor nicotianae spp.Isolates, with the gummosis, and root and foot rot symptoms observed in infected Khasi mandarin orchards.These results agree with the observations made by Das et al. (2017) during their study [10].Similarly, excised stems and shoots with inoculation showed brown lesions over the inoculated areas compared with the uninoculated shoots showing whitish-coloured vascular bundles.These symptoms are in accordance with findings from the extensive work of Das et al. [19].Das et al. [19] reported that within 7 days after treatment, Kinnow mandarin leaves developed brownish-black necrotic lesions, indicating that Phytophthora spp. is the root cause of the citrus gummosis disease.Similarly, the pathogenicity assay employing excised stems conducted by Das et al. [19].(2017) in citrus revealed that the isolate could cause a progressive necrotic lesion in the Kinnow mandarin stem from the site of inoculation [15].
The population dynamics of microbes in diseased and healthy rhizospheres varied significantly, with the latter maintaining a marginally higher population, comprising both bacterial as well as fungal counts.However, it was not possible to establish which group of microbes modulated the pathosystem or host resistance (Table 4).
We could identify as many as ten (10) bacterial isolates and eighteen (18) fungal isolates in the rhizosphere of diseased-infected plants, while the healthy rhizosphere displayed thirteen (13) fungal isolates and twenty-eight (28) bacterial isolates with distinct morphology.However, the diseased rhizosphere showed a much higher Shannon-Wiener diversity index, indicating its higher diversity index, besides that for mean species evenness due to the presence of opportunistic microorganisms in the diseased rhizosphere in response to plant secretomes carrying root exudates.Root exudates can have a positive or negative effect on fungal community structure, depending upon the nature of substrates or signalling molecules, to support an increase in soil fungal taxa for a higher diversity index [34], while in the case of disease suppressiveness, beneficial microbes are more likely to remain abundant in the root zone of plants, suppressing the growth of harmful and other opportunistic microorganisms, thereby lowering the diversity index of healthy rhizospheric soil [35].Interestingly, culturable microbes, when compared on the basis of species distribution between healthy and diseased rhizospheres, revealed no differentiating clues, but the pattern of distribution of fungal species led us to identify the source of suppressiveness against Phytophthora nicotianae.
The suppressive ability of different microbial strains is facilitated through a variety of mechanisms of action, including direct antagonistic behaviour, competition for resources, and indirect antagonistic behaviour through the induction of plant resistance [36,37].However, to fully understand the citrus plant-microbiome interaction and its resultant susceptibility/suppressiveness, a detailed study on the community composition and nature of root exudates is crucial.Some important clues for exploiting the citrus microbiome for disease resistance against soil-borne pathogens can be found in reviews by Srivastava et al. [6] and Zhang et al. [38].
The metagenomic analysis results of the healthy and diseased plants are in sync with the cultural microbiome data establishing the lower intensity of Phytophthora spp. in the rhizosphere of disease-free plants compared with that in the infected Khasi mandarin rhizosphere.This established the presence of Phytophthora spp., leading to the disease.This study also identified some of the microbes and fungal populations varying in accordance with the health status of the plant, further supporting the hypothesis that community composition will affect or becomes affected by the disease in the same field and cultivar.Further analysis such as that on differences in the enzymes activated and the metabolic network will give a better picture of how bacterial community plays a role in disease suppression.With the current data size of each sample (5 GB), despite them revealing the dominating microbes, arriving to any conclusions will be difficult.Taking clues from here, further in-depth metagenomic sequencing will be conducive to gaining insights into either the microbe-microbe interaction or the host-microbe interaction at the molecular and metabolic levels for better elucidations.
Our analysis found that Ascomycota and Basidiomycota were the most prevalent fungal phyla in all the soil samples, but the quantity of Ascomycota was dramatically reduced in diseased soils.The genera Fusarium, Aspergillus, Saitozyma, Rhizopus, Exophiala, Rhizophagus, Collietotrichum, and Purpureocillium were found to have the highest relative abundance in Phytophthora-infected soil, while Trichoderma, Penicillium, Linnemannia, Mortierella, Talaromyces, Globisporangium, Saitozyma, and Serendipita were most abundant in the diseasefree rhizosphere, suggesting their putative role in the natural suppression of Phytophthora spp.These studies will strengthen our understanding of the already available community structure associated with Phytophthora [39,40].
Although increasing bacterial diversity is a better predictor of disease suppression than is the presence of a specific bacterial taxon, our study found that the most abundant phyla, such as Proteobacteria, Terrabacteria, Actinobacteria, FCB, Flavobacteria, Patescibacteria, Cyanobacteria, PVC, and Acidobacteria, were present in both types of soil samples.Beneficial effects of the plant microbiota are frequently described as being provided by synergistic interactions between microbes, and it is important to identify and understand the microbiome that has a synergistic effect to explore the maximum expression of crop resilience to disease [41].However, it is difficult to determine which functional groups contribute to suppressiveness.

Conclusions
Our study provided important insights into the changes in the rhizosphere microbiome of Khasi mandarin infected with Phytophthora nicotianae, the causative agent of foot and root rot disease, as established via pathogenicity tests.With the comparatively low depth of our metagenomic data, most of the results are in accordance with some of the metagenomic revelations reported earlier.Regardless of the health status of a rhizosphere (though higher culturable microbial counts prevailed in the healthy rhizosphere compared with the diseaseinfected rhizosphere), unknown pathogenic species currently referred to as Candidatus predominantly dictated the rhizosphere functions, but a higher abundance of Bacteriodetes and Pseudomonas in the healthy rhizosphere provided strong clues with which to identify the microbial genes, the source of disease resistance.Likewise, identifying functional fungal microbial communities as a source of disease suppression could potentially pave the way for developing newer methods for the consortium-mediated management of Phytophthora-induced diseases in citrus crops, supported by robust screening efforts with focus on the endophytic behaviour of identified microbes.An in-depth metagenomic study of healthy versus disease-infected rhizospheres is further suggested to derive better results with regard to developing microbial antagonists as aids in citrus health management.

Figure 3 .
Figure 3. Pure cultures of the isolates of Phytophthora spp.(A-C) with different growth patterns.(A) Light cottony mycelial growth; (B) dense, fluffy, and cottony mycelial growth; and (C) light cottony mycelial growth.

Figure 3 .
Figure 3. Pure cultures of the isolates of Phytophthora spp.(A-C) with different growth patterns.(A) Light cottony mycelial growth; (B) dense, fluffy, and cottony mycelial growth; and (C) light cottony mycelial growth.

Figure 3 .
Figure 3. Pure cultures of the isolates of Phytophthora spp.(A-C) with different growth patterns.(A) Light cottony mycelial growth; (B) dense, fluffy, and cottony mycelial growth; and (C) light cottony mycelial growth.

Figure 6 .
Figure 6.Pathogenicity of Phytophthora nicotianae on leaf inoculation ((A) as the control; (B) as the treated sample), excised shoot inoculation ((C) as the control; (D) as the treated sample), and root inoculation ((E) as the control; (F) as the treated sample).The effectiveness of P. nicotianae was marked by fungal growth on the inoculated leaf after 7 days.(B) Overlapping mycelial growth on excised shoot after 15 days.(D) Loss of feeder roots after 72 h of pathogen inoculation, (F) marked by the appearance of blackish-brown necrotic lesions over the roots, while the untreated control (E) remained symptomless.

Figure 6 .
Figure 6.Pathogenicity of Phytophthora nicotianae on leaf inoculation ((A) as the control; (B) as the treated sample), excised shoot inoculation ((C) as the control; (D) as the treated sample), and root inoculation ((E) as the control; (F) as the treated sample).The effectiveness of P. nicotianae was marked by fungal growth on the inoculated leaf after 7 days.(B) Overlapping mycelial growth on excised shoot after 15 days.(D) Loss of feeder roots after 72 h of pathogen inoculation, (F) marked by the appearance of blackish-brown necrotic lesions over the roots, while the untreated control (E) remained symptomless.

Table 3 .
Population enumeration of culturable rhizospheric microbes in Phytophthora-infected and healthy Khasi mandarin rhizospheres.

Table 4 .
Percentage relative abundance of isolated rhizospheric fungi in Phytophthora-infected versus healthy Khasi mandarin rhizospheres.

Table 5 .
Contig assembly statistics at a minimal length of 1000 bp/1 kbp.