Tube-Wise Diagnostic Microarray for the Multiplex Characterization of the Complex Plant Pathogen Ralstonia solanacearum

Ralstonia solanacearum is a well-known agricultural and ecological threat worldwide. The complexity of the R. solanacearum species complex (Rssc) represents a challenge for the accurate characterization of epidemiological strains by official services and research laboratories. The majority of protocols only focus on a narrow range of strains; however, this species complex includes strains that represent major constraints and are under strict regulation. The main drawback associated with the current methods of detecting and characterizing Rssc strains is their reliance on combining different protocols to properly characterize the strains at the ecotype level, which require time and money. Therefore, we used microarray technology (ArrayTube) to develop a standard protocol, which characterizes 17 major groups of interest in the Rssc, in a single multiplex reaction. These 17 majors groups are linked with a phylogenetic assignation (phylotypes, sequevars), but also with an ecotype assignation associated with a range of hosts (e.g., brown rot, Moko). Probes were designed with a 50-mer length constraint and thoroughly evaluated for any flaws or secondary structures. The strains are characterized based on a DNA extraction from pure culture. Validation data showed strong intra-repeatability, inter-repeatability, and reproducibility as well as good specificity. A hierarchical analysis of the probe groups is suitable for an accurate characterization. Compared with single marker detection tests, the method described in this paper addresses efficiently the issue of combining several tests by testing a large number of phylogenetic markers in a single reaction assay. This custom microarray (RsscAT) represents a significant improvement in the epidemiological monitoring of Rssc strains worldwide, and it has the potential to provide insights for phylogenetic incongruence of Rssc strains based on the host of isolation and may be used to indicate potentially emergent strains.


INTRODUCTION
Bacterial strains of the Ralstonia solanacearum species complex (Rssc) are considered priority plant pathogens in many countries in the world and are classified as quarantine organisms as well as bioterrorism and double usage agents in the USA and in Europe (2000; 2008). Plant protection across borders and territories requires the availability of rapid and reliable protocols to prevent the accidental or intentional introduction of exotic strains. The genetic diversity of Rssc strains is distributed into four major phylotypes and 53 sequevars based on partial sequences of the endoglucanase (egl) gene (Fegan and Prior, 2005). The number of sequevars or lineages can be extended with the discovery of new strains, as demonstrated in sequevars that were identified and characterized in Martinique (French West Indies; Wicker et al., 2009) and Brazil (Albuquerque et al., 2014). A recent polyphasic taxonomic approach based primarily on DNA-DNA hybridizations, Internal transcribed spacer (ITS) and egl partial sequence phylogenetic analyses, and phenotypic data merged the Rssc into three species (Safni et al., 2014): R. solanacearum clustering strains from phylotype II; R. pseudosolanacearum clustering strains from phylotypes I and III; and R. syzygii subsp. syzygii comb. nov. clustering strains from phylotype IV, which included two subspecies, R. syzygii subsp. celebesensis subsp. nov. and R. syzygii subsp. indonesiensis subsp. nov. These data were supported by comparative analyses of whole bacterial genomes (Prior et al., 2016).
Such complexity represents a challenge for any diagnostician attempting to develop protocols and tools to specifically discriminate Rssc strains. Moreover, the ability to identify and distinguish Rssc ecotypes, a group of bacteria sharing the same ecological niche (Cohan, 2002), is important because they represent the most epidemiological active strains in the Rssc causing "famous" diseases, such as potato brown rot or the banana Moko, in valuable agronomic plants. Certain ecotypes can be split into different phylogenetic groups, such as the paraphyletic Moko strains pathogenic to Musaceae and Solanaceae (Cellier et al., 2012;Albuquerque et al., 2014), which have been phylogenetically categorized into sequevars 3, 4, 6, 24, 25, 41, and 53. Moreover, the difficulty of discriminating among ecotypes increases when high proximity is observed in the genome content of the ecotypes, which has been observed for the pathological variant IIB-4NPB strains (not pathogenic on banana) and Moko sequevar 4 strains (Cellier et al., 2012;Ailloud et al., 2015).
Diagnostic protocols mainly rely on techniques that offer compatible routine application procedures and provide for userfriendly environments, good industrial production capacities, and a high level of reliability. Techniques that primarily employ immunoassays and immune-strip tests (ImmunoStrip R ; Agdia, Elkhart, IN) are mainly used for the detection of Rssc strains at the species level (Danks and Barker, 2000), and although these assays are simple and affordable, they are known to produce false positives (Narayanasamy, 2010). DNA-based approaches using conventional PCR (Table 1) can identify the Rssc at the species level (Huang and Schell, 1990;Seal et al., 1993;Opina et al., 1997;Lee and Wang, 2000;Glick et al., 2002;Schonfeld et al., 2003), the phylotype level (Fegan and Prior, 2005); and when these approaches are coupled with sequencing capabilities, they can be used to produce phylogenetic trees (Lane, 1991;Taghavi et al., 1996;Fegan and Prior, 2005;Prior and Fegan, 2005b). General protocols that rely on the evolution of PCR techniques are also available, such as qPCR (Weller et al., 2000;Ozakman and Schaad, 2003;Smith and De Boer, 2009;Inoue and Nakaho, 2014) or LAMP PCR . Specific protocols have been designed to target a particular group of strains or ecotypes of the Rssc, such as brown rot strains (Fegan et al., 1998;Weller et al., 2000;Ozakman and Schaad, 2003;Smith and De Boer, 2009;Kubota et al., 2011;Ha et al., 2012;Li et al., 2014;Kubota and Jenkins, 2015;; Moko strains (Prior and Fegan, 2005a;Cellier et al., 2015); and Blood Disease Bacterium (BDB) strains (Kubota et al., 2011). These diagnostic methods are primarily limited to the detection of strains associated with the brown rot ecotype (Li et al., 2014); however, the Rssc presents a wide diversity of strains.
The main drawback associated with the current methods of detecting and characterizing Rssc strains relies in combining different protocols to properly characterize the strains at their ecotype level. The development of a single method to fully characterize the main groups/ecotypes of the Rssc would greatly benefit the diagnostic field and streamline the workflow for strain identification in an efficient and affordable manner, and it would also facilitate the detection of emerging or unknown strains.
In this study, we employed a multiplex approach to manage the complexity of the Rssc within the framework of microarray technology, which has the potential to test for multiple biomolecular targets in a single reaction. In recent years, applications for the Rssc have been developed for both fundamental (Guidot et al., 2007(Guidot et al., , 2009Cellier et al., 2012;Lefeuvre et al., 2013) and applied research (Aittamaa et al., 2008;Pelludat et al., 2009;Dobnik et al., 2014). To provide a costeffective strategy that meets the need for user-friendly processing via conventional lab equipment and high-volume manufacturing capacities that comply with in vitro diagnostic (IVD) regulations, the ArrayTube (AT) platform (Alere Technologies GmbH, Jena, Germany) was selected as the final protocol. This array consists of a custom microarray integrated into a micro reaction vial, which simplifies handling and is used for routing testing in the medical field (Braun et al., 2012;Schneeberg et al., 2015).
Recently, new genomes of the Rssc were fully sequenced and annotated , thus providing new insights for the development of more reliable markers for diagnostic purposes. This paper provides the first description of the production of a microarray (RsscAT) dedicated to characterize strains within 17 major groups of interest in the Rssc with respect to their phylogenetic assignation and ecotype. This portable custom RsscAT provides significant improvements that allow for the rapid assignment of the phylogenetic position of a strain, which can be used to predict the pathogenicity of a strain. RsscAT has the potential to be used for the epidemiological monitoring of heterogeneous Rssc strains worldwide.

Probe Design and Selection
The genomes used to design the probes are freely available on the MicroScope platform (available online at https://www. genoscope.cns.fr/agc/microscope/home/index.php; Vallenet et al., 2009;Genoscope, Evry, France) and the EMBL nucleotide sequence database . A first bioinformatics screening of candidate CDSs was performed using both the comparative genomics tool of the MicroScope platform and a BLAST search against NCBI databases to select CDSs that were conserved in the different phylogenetic subgroups and presented limited or absent identities with DNA sequences from non-target genomes present in the databases. A second bioinformatics selection eliminated the CDSs that were too close to mobile elements (e.g., transposases, integrases, etc.) using a parsing R (v3.0.0) script (R Development Core Team, 2011). The probes (1-9 per CDS) were designed using two software programs, the Array Oligo Selector (AOS; Bozdech et al., 2003) and OligoArray (OA; Rouillard et al., 2003), according to the following specifications: 50-mer length, 40-70% GC content, 2 probes per CDS, binding energy ranging between −20 and −35 kcal/mol (AOS), Tm 68-88 • C (OA), exclusion of probes including 20 AT or more (AOS), and exclusion of probes including five consecutive repetitions of one of the four bases (OA). The probes were selected from a file that included all of the previously selected CDSs vs. a file that compiled all of the R. solanacearum genomes. A total of 256 candidate probes were generated, and their specificity was verified using a BLAST search against NCBI databases and a specific database that includes Rssc genomes and nontarget genomes (strains genetically close to the Rssc or strains occurring in the same ecological niche (Supplementary Table  1). Primer pairs flanking the candidate probes were designed using Primer 3 (Untergasser et al., 2012), and PCR assays were performed with the selected strains representative of the Rssc diversity ( Table 3) to validate the occurrence of the targeted DNA in the different subgroups. Following the PCR screening that removed the candidate DNA regions that were improperly amplified by the target strains, a batch of 100 probes was selected to be implemented and cribbed in the AT technology. This last selection removed probes that yielded unexpected results and generated the final set of 32 validated specific biological probes (Table 4), including a negative control probe that verifies the hybridization process. This set of 32 probes/primers was used to multiplex the 17 groups of interest (RsscAT).
The RsscAT were manufactured by Alere Technologies and consisted of oligonucleotide probes with a 3 ′ amino modification and C6 spacer. The probes were spotted in duplicate (first batch-AT1) along with other non-specific probes that were not analyzed; and the probes were also spotted in triplicate (second batch-AT2).

Multiplex Linear DNA Amplification and Labeling
Strains were grown overnight in Luria-Bertani broth (LB) at 28 • C with agitation at 250 rpm, and then mild centrifugation was performed at 5,100 rpm for 10 min at 10 • C. To prepare the cells for a DNA extraction, the pellets were washed twice with 500 µL of 1 M NaCl solution, centrifuged at 5,100 rpm for 10 min at 10 • C, and then processed using the Wizard Genomic DNA Purification Kit (Promega, Madison WI, USA) following the manufacturer's instructions.
Following this DNA extraction, a linear PCR amplification was performed in which only one reverse primer was used to linearly amplify each targeted sequence ( Table 4). The labeling of the genomic DNA was accomplished during the linear amplification step by using dUTP linked biotin as a marker, which allowed site-specific internal labeling of the corresponding target region, thus leading to the production of single-stranded biotin-labeled products. A total of 25 individual reverse primers were used in a multiplex linear PCR amplification for each Rssc strain. Using the HybridisationPlus Kit (Alere Technologies), 5 µL of DNA normalized with 200 ng/µL of genomic DNA were labeled according to the manufacturer's recommendations (3.9 µL Labeling Buffer B1, 0.1 µL Polymerase B2) and 1.35 µM of each of the 25 primers. The linear PCR amplification was performed using a Veriti Thermal Cycler (Applied Biosystems,

Microarray Hybridization and Data Acquisition
The reagents required for the hybridization steps were provided by the HybridisationPlus Kit (Alere Technologies), and the protocols were performed according to the manufacturer's recommendations except for certain conditions that were optimized according to the probe design requirements. The prewashing steps of the RsscAT consisted of adding 500 µL of 60 • C preheated ultrapure water for a 5 min incubation at 60 • C and stirring at 550 rpm. The flow was discarded, and then 200 µL of hybridization buffer C1 preheated at 60 • C was added and incubated for 2 min at 60 • C with stirring at 550 rpm. The flow was discarded, and then 100 µL of hybridization mix consisting of 10 µL of labeled DNA and 90 µL of the hybridization buffer (C1) preheated at 60 • C was transferred into a prewashed RsscAT and then incubated for 1 h incubation at 60 • C with stirring at 550 rpm.
The washing steps consisted of heating the hybridized RsscAT at 42 • C and discarding the hybridization mix, and then 500 µL of preheated washing buffer 1 (C2) at 42 • C was added and incubated for 5 min at 42 • C with stirring at 550 rpm. The flow was then discarded, and these steps were repeated two more times, with the last repetition performed at 30 • C. The conjugate mix, which consisted of 1 µL of preheated HRP conjugate 100x (C3) and 99 µL of preheated conjugate buffer (C4) at 30 • C, was added to the washed RsscAT and then incubated for 10 min at 30 • C with stirring at 550 rpm. Another washing step consisted of discarding the flow, adding 500 µL of washing buffer 2 (C5), and then incubating for 5 min at 30 • C with stirring at 550 rpm. The final coloration step consisted of discarding the flow, adding 100 µL of HRP substrate buffer (D1) that was already acclimated to room temperature, and incubating for 5 min at 30 • C without stirring.
Data acquisition was performed with the software provided with the ATR03 scanner (Alere Technologies) following the manufacturer's instructions.

Data Analysis
Hybridization signals were processed using the IconoClust software version 3.3, and all of the spots were automatically normalized by the software according to the equation NI = 1 − (M/BG), where NI is the normalized intensity, M is the average intensity of the automatically recognized spot, and BG is the intensity of the local background. The output range of the signals was between 0 and 1, with 0 being negative and 1 being the maximal possible signal value. These raw data were processed using a home-made R script (provided upon request) that could set a threshold to discriminate positive probe signals from negative probe signals. The threshold was manually assigned to 0.6, which allowed for clear discrimination of positive and negative signals. The threshold for the control probe GC_07 (Rssc strains) was set to 0.4. Signals with intensities higher than the assigned threshold were considered positive and set as "1." Signals lower than the assigned threshold were regarded as negative and set as "0." Statistics Different replicates were performed to assess the repeatability of the protocol. Each probe was spotted in duplicate (batch AT1) or triplicate (batch AT2) to evaluate the intra-assay repeatability, which was calculated based on the consistency rate within each batch or for both batches. Several strains (2 for batch AT1 and 6 for batch AT2) were tested twice or three times in the same batch to evaluate the inter-assay repeatability, which was calculated based on the proportion of consistent results among the different hybridization repeats.
Twenty-seven strains were tested in both the AT1 and AT2 batches. For these strains, the relative specificity, SP = TN TN + FP ; the relative sensitivity, SE = TP TP + FN ; and the Diagnostic Odds Ratio, DOR = TP/FP FN/TN and the corresponding 95% confidence interval (95% CI), were calculated for each batch, where TP, FP, FN, and TN refer to true positive, false positive, false negative and true negative, respectively. The DOR value ranged from 0 to infinity, with higher values indicating better discriminatory test performance. The DOR values were compared between the two batches using the Breslow-Day test of the Homogeneity of Odds Ratios (p = 0.05). Finally, the SP and SE were calculated for each group of interest using the data acquired for the 75 target strains and the 12 outgroups in the two batches.

RESULTS AND DISCUSSION
The RsscAT was developed to evaluate diverse bacterial strains and represents an affordable, accurate, fast, and user-friendly diagnostic tool with a high potential for standardization and routine diagnostic testing use. Microarrays developed using the conventional glass slide format with fluorescence detection (e.g., Cy3/Cy5) have a high operational burden because of the high laboratory skill required to handle the material and perform the bioinformatics computations, as well as the high cost of the equipment. Although high-density microarrays are more likely to meet research laboratory expectations, routine laboratory diagnostics require more user-friendly and efficient protocols. AT format microarrays provide a valuable option for diagnosticians who need to efficiently test a large number of markers at the same time. The ability of the AT microarray to identify and discriminate strains of the Rssc was evaluated using a set of 75 target strains representing the 17 main groups of the Rssc and 12 outgroup strains that are found in the same environment and show a high likelihood of interfering with the identification of Rssc strains (pathogenic or saprophytic strains).

Probe Design
The first step in developing the RsscAT was to design accurate probes specific to the 17 phylogenetic groups of the Rssc. This design step was performed in several stages during which our strong focus on achieving the specificity of the selected templates (CDS or probes) was maintained via in silico analyses and experimental data.
Probes were designed with a 50-mer length constraint, which provided sufficient length for a probe that balances specificity and flexibility and is capable of identifying the best sequences markers (Palmer et al., 2006).
The probability of secondary structures (e.g., hairpins) was checked during the probe design and used to select candidate probes and then enhance the probe signals. Additionally, the use of single-stranded DNA as a template prevented competition between the probe and antisense strand and increased the probability of single-stranded DNA hybridization to the probe.
To our knowledge, we have designed the first probe set that can characterize Rssc ecotypes and strains regardless of the array technology, including macroarrays, microarrays, glass slides, nylon, or AT. Although various strain typing designs used in the medical field provide detection at levels below the species level (Schmoock et al., 2011;Braun et al., 2012), most of the designs available for Rssc strains have focused on detection at the species level (Aittamaa et al., 2008;Pelludat et al., 2009;Dobnik et al., 2014); moreover, these designs are constrained by the sensitivity of the tests, which is a drawback of array technology. Here, RsccAT was presented as a characterization protocol, and it can be integrated at the second step of strain identification after the steps for isolating and detecting the target Rssc strain (by ELISA, PCR, etc.). Therefore, the sensitivity constraint is not applicable because a significant amount of DNA can be extracted from the pure culture strains.

Control Responses
The internal positive controls from Alere Technologies generated the expected positive results, and the negative control probe NEG_Cupri07 generated the expected negative results for each hybridization. The control probe GC_07 for the Rssc strains generated the expected results for most of the strains; however, the following 4 phylotype III strains showed a lower intensity signal: 447, 594, MG464 (Madagascar), and CFBP4964 (Reunion). This result suggested the presence of mismatches between the probe and the amplified targets. Because a signal was observed between the expected positive and expected negative strains, we lowered the threshold for this specific probe to 0.4. Using this threshold, all of the strains of the Rssc tested positive (100% inclusivity). This last result was expected because the probe GC_07 is included in the UDP-3-O-acyl-GlcNAc deacetylase (lpxC) gene, which has been successfully used in the design and validation of specific primers for the Rssc (Opina et al., 1997;Villa et al., 2003;Fegan and Prior, 2005). Additionally, the 12 outgroup strains did not hybridize on any of the 32 developed probes except the R. pickettii strain LMG5942, which was detected by the probes III_04 and BDB_02 (83.3% exclusivity). Nevertheless, this false positive data was not considered because the R. solanacearum positive control probe GC_07 did not test positive for this strain. The results obtained with the control probes confirmed the reliability of the hybridization step and supported the high specificity of the RsscAT.
No significant results were observed between the odds ratios for the batches (Breslow-Day test, p = 0.4375), which indicated that the different strains presented globally similar responses when the two batches were tested; however, the common probes of the two arrays were spotted under different conditions for these two batches (different locations and repetitions), which indicates the high reproducibility of this AT microarray protocol.
The ability of the AT microarray to identify and discriminate strains of the Rssc was evaluated using a set of 75 target strains representing the 17 main groups of the Rssc and 12 outgroup strains that are found in the same environment and show a high likelihood of interfering with the identification of Rssc strains (pathogenic or saprophytic strains).

Phylotype Identification
High relative specificity of 98.07% (95% CI 96.65-99) and sensitivity values of 100% (95% CI 90.26-100) were obtained for both probes that typed phylotype I. Indeed, the probes designed to detect phylotype I yielded the expected specific positive results for all of the tested strains of phylotype I, and they only yielded a false positive signal for the strains CFBP1418 (IIB-4NPB) and JT663 (SZY).
However, high relative specificity and sensitivity values were obtained for the probes that detected phylotypes IIA and IIB, with 99.24% (95% CI 97.27-99.91) and 100% (95% CI 99.12-100) relative specificity for phylotypes IIA and IIB,respectively, relative sensitivity for phylotypes IIA and IIB, respectively. Thus, the target strains for phylotypes IIA and IIB were successfully detected except for strain CIP239 belonging to IIA-40.
All the strains from phylotype III were successfully detected by the two dedicated probes, and a 100% relative sensitivity value (95% CI 95.49-100) was observed. A high relative specificity value of 99.49% (95% CI 98.51-99.89) was also obtained, and false positive signals were only registered for probe III_04 when testing LMG 5942 (see Section Controls).
All of the phylotype IV target strains except strain JT663 were detected, and false positive values were not registered, thereby resulting in a 100% relative specificity value (95% CI 99.38-100) and a high relative sensitivity value of 93.67% (95% CI 85.84-97.91).
To conclude, the specific probes yielded the expected results and showed good specificity. The only group that showed lower reliability was the phylotype II probes, which did not yield true positive signals for the Moko IIB strains and epidemiological IIB-4NPB strains. The phylogenetic and genomic characteristics of the Moko IIB-4 strains and IIB-4NPB strains are closely correlated , whereas these characteristics of the Moko IIB-3 strains are less correlated. The lack of sensitivity of the phylotype II probes for this particular group of strains was not predicted by the computer analysis or PCR amplification of the target. Nevertheless, the results obtained with the probes targeting the two subgroups IIA and IIB inside phylotype II largely overcame the inclusivity problems shown with the phylotype II probes.

Ecotype and Lineage Identification
The ecotype of the Rssc strains (e.g., brown rot, Moko, NPB, BDB, and RSY) is the lowest level of genetic diversity that can be detected with the RsscAT, and an accurate ecotype analysis must be performed with additional groups of probes. Hence, a hierarchical analysis of the probe groups must be performed to develop an accurate characterization of the target strain. Once the status of the Rssc strain is validated, the phylotype can be determined and then, if possible, the ecotype is also determined accordingly to a particular ecotype. This process can circumvent the interpretation of false negative/positive signals.
The results ( Table 5) showed that the Grandville wilt ecotype from phylotype IIA-7, the epidemiological variant IIB-4NPB, and the BDB ecotype were fully and specifically detected.
Brown rot strains were fully characterized by the dedicated specific probes by considering the complementary responses displayed by the four probes. Indeed, false negative responses obtained for the probes IIB1_05 and IIB1_06 were overcome by the IIB1_10 and IIB1_12 responses and vice versa. This result highlighted the importance to have several probes for a given group of interest, located if possible on different CDS, which was the case for the Brown-rot strains-targeting probes. Interestingly, the IIB1_10 and IIB1_12 probes matched with a specific genome region selected with a different in silico approach and used to develop a Brown rot strains specific real-time quantitative PCR assay . The paraphyletic Moko strains are sorted into different clusters within phylotype II (IA-6, IIA-24, IIB-3, and IIB-4), which highlights the genomic complexity and heterogeneity of this ecotype. Although these strains are epidemiologically adapted to wilt Musaceae, these strains also show pathogenicity to Solanaceae in controlled conditions (Cellier and Prior, 2010). Moko sequevars 4 and 6 were specifically detected, whereas sequevars 3 and 24 showed false negative signals despite the amplification of the respective target sequences by PCR. The hybridization process between the probe and the amplified DNA from these three strains might have been compromised by substitution of the DNA sequence, thereby leading to false negative signals. The hybridization parameters have been fully discussed in Tomlinson et al. (2014). Several site substitutions at specific regions of the amplified DNA can lead to a significant decrease in the positive signal, and the same false negative characterization at the ecotype level was observed for strain JT663 in the SZY ecotype and MAFF301558 in the RsIV ecotype. Complementary specific PCRs could be used when information at the ecotype level is lacking ( Table 1).
Several false negative signals were obtained for some ecotypes; however, as previously noted, a hierarchical analysis of the probe groups could help to resolve such strain assignments. For example, strain CMR15 has been typed on AT as an IIB-1 strain by the IIB1_06 probe; however, this strain has been identified as a phylotype III strain and could not be an IIB-1 strain.
Additionally, strain CFBP6783 has been identified as a probable IIB-1 strain; however, it belongs to the IIB-4NPB group. Nevertheless, this result is not reliable (non-homogeneous repetitions of two IIB-1 probes and no identification by the two other IIB-1 probes) because all IIB-4NPB probes indicate that this strain is actually an IIB-4NPB strain.

Intra-Repeatability, Inter-Repeatability, and Reproducibility
A consistency rate of 99.65% was obtained for the duplicate spotted probes for the AT1 batch and 99.97% was obtained for the triplicate spotted probes in the AT2 batch; thus, the RsscAT method showed a high degree of intra-array repeatability. In addition, high inter-array repeatability values were also observed, with consistency rates of 97.06% (batch AT1) and 99.51% (batch AT2). For the set of 27 strains analyzed from both batches (AT1 & AT2), the major indicators of diagnostic test performance presented high values. Thus, the relative specificity, relative sensitivity, odds ratios and the corresponding 95% confidence intervals were 98.81% (95% CI 98. 16-99.27 CI 869.92-2474.9), respectively, when the strains were hybridized on batch AT2.
No significant results were observed between the odds ratios for the batches (Breslow-Day test, p = 0.4375), which indicated that the different strains presented globally similar responses when the two batches were tested; however, the common probes of the two arrays were spotted under different conditions for these two batches (different locations and repetitions), which indicates the high reproducibility of this AT microarray protocol.

CONCLUSION
The design presented here is highly reliable, and the repetitions of probes as well as the different levels of characterization were able to differentiate among the target strains within the 17 different groups of high importance within the Rssc. This design is also able to detect phylogenetic incongruence among strains according to the host isolation, and then can also detect emergent strains that will require further confirmation tests. Different outcomes and scenarios could emerge from this RsccAT analysis: for example, strains characterized as Moko on tomato (Solanum lycopersicum) represent serious cases and require immediate action; however, tests that rely on species-level analyses would give the positive result for the detection of "R. solanacearum" strains and a negative result for the specific detection of brown rot strains. Additionally, the ecotype level represents strain groups that are epidemiologically active worldwide; hence, the strains identified at the phylotype level but not at the ecotype level represent anomalies that require further investigation.
Currently, only limited eradication and containment strategies are available to control Rssc strains once they have been introduced into a territory. Hence, preventing such introductions and developing tools to identify and characterize the quarantined organisms remain priorities for the agricultural industry. The developed RsscAT dedicated to identifying Rssc strains was designed to minimize time and costs and maximize marker multiplexing and reliability, and it is only applied on pure culture strains. Multiplexing the markers into a compact, affordable, and user-friendly technology for infra-species-level characterizations of Rssc strains is the goal of the RsccAT design. A total of 17 groups can be interpreted in a phylogenetic hierarchical fashion that covers the main phylogenetic groups and ecotypes of the Rssc. Although ELISA, conventional or real time PCR, LAMP-PCR, and other molecular tools are designed to test for a multitude of strains in one reaction based on one or very few markers, the RsscAT is design to test for a multitude of markers in one reaction for one strain. This profile makes the RsscAT as a strong complementary diagnostic protocol within a diagnostic scheme. The RsscAT diagnostic array presents a number of benefits, such as an industry-standard format and a high level of reliability, and it also provides diagnostic laboratories with a user-friendly, cost-efficient and time-efficient option that can be used for epidemiological monitoring of Rssc strains.