Use of propidium monoazide for live/dead distinction in microbial ecology

One of the prerequisites of making ecological conclusions derived from genetic fingerprints is that bacterial community profiles reflect the live portion of the sample of interest. Propidium monoazide is a membrane-impermeant dye that selectively penetrates cells with compromised membranes, which can be considered dead. Once inside the cells, PMA intercalates into the DNA and can be covalently cross-linked to it, which strongly inhibits PCR amplification. By using PCR after PMA treatment, the analysis of bacterial communities can theoretically be limited to cells with intact cell membranes. Four experiments were performed to study the usefulness of PMA treatment of mixed bacterial communities comprising both intact and compromised cells in combination with end-point PCR by generating community profiles from the following samples: (i) defined mixtures of live and isopropanol-killed cells from pure cultures of random environmental isolates, (ii) wastewater treatment plant influent spiked with defined ratios of live and dead cells, (iii) selected environmental communities, and (iv) a water sediment sample exposed to increasing heat stress. Regions of 16S rRNA genes were PCR amplified from extracted genomic DNA, and PCR products were analyzed by using denaturing gradient gel electrophoresis (DGGE). Results from the first two experiments show that PMA treatment can be of value with end-point PCR by suppressing amplification of DNA from killed cells. The last two experiments suggest that PMA treatment can affect banding patterns in DGGE community profiles and their intensities, although the intrinsic limitations of end-point PCR have to be taken into consideration.     

DNase I and Proteinase K eliminate DNA from injured or dead bacteria but not from living bacteria in microbial reference systems and natural drinking water biofilms for subsequent molecular biology analyses

Molecular techniques, such as polymerase chain reaction (PCR) and quantitative PCR (qPCR), are very sensitive, but may detect total DNA present in a sample, including extracellular DNA (eDNA) and DNA coming from live and dead cells. DNase I is an endonuclease that non-specifically cleaves single- and double-stranded DNA. This enzyme was tested in this study to analyze its capacity of digesting DNA coming from dead cells with damaged cell membranes, leaving DNA from living cells with intact cell membranes available for DNA-based methods. For this purpose, an optimized DNase I/Proteinase K (DNase/PK) protocol was developed. Intact Staphylococcus aureus cells, heat-killed Pseudomonas aeruginosa cells, free genomic DNA of Salmonella enterica, and a mixture of these targets were treated according to the developed DNase/PK protocol. In parallel, these samples were treated with propidium monoazide (PMA) as an already described assay for live-dead discrimination. Quantitative PCR and PCR-DGGE of the eubacterial 16S rDNA fragment were used to test the ability of the DNase/PK and PMA treatments to distinguish DNA coming from cells with intact cell membranes in the presence of DNA from dead cells and free genomic DNA. The methods were applied to three months old autochthonous drinking water biofilms from a pilot facility built at a German waterworks. Shifts in the DNA patterns observed after DGGE analysis demonstrated the applicability of DNase/PK as well as of the PMA treatment for natural biofilm investigation. However, the DNase/PK treatment demonstrated some practical advantages in comparison with the PMA treatment for live/dead discrimination of bacterial targets in drinking water systems.     

Comparison of propidium monoazide with ethidium monoazide for differentiation of live vs. dead bacteria by selective removal of DNA from dead cells

The differentiation between live and dead bacterial cells presents an important challenge in many microbiological applications. Due to the persistence of DNA in the environment after cells have lost viability, DNA-based detection methods cannot differentiate whether positive signals originate from live or dead bacterial targets. We present here a novel chemical, propidium monoazide (PMA), that (like propidium iodide) is highly selective in penetrating only into 'dead' bacterial cells with compromised membrane integrity but not into live cells with intact cell membranes/cell walls. Upon intercalation in the DNA of dead cells, the photo-inducible azide group allows PMA to be covalently cross-linked by exposure to bright light. This process renders the DNA insoluble and results in its loss during subsequent genomic DNA extraction. Subjecting a bacterial population comprised of both live and dead cells to PMA treatment thus results in selective removal of DNA from dead cells. We provide evidence that this chemical can be applied to a wide range of species across the bacterial kingdom presenting a major advantage over ethidium monoazide (EMA). The general application of EMA is hampered by the fact that the chemical can also penetrate live cells of some bacterial species. Transport pumps actively export EMA out of metabolically active cells, but the remaining EMA level can lead to substantial loss of DNA. The higher charge of PMA might be the reason for the higher impermeability through intact cell membranes, thus avoiding DNA loss.     

Selective quantification of viable Escherichia coli bacteria in biosolids by quantitative PCR with propidium monoazide modification

Quantitative differentiation of live cells in biosolids samples, without the use of culturing-based approaches, is highly critical from a public health risk perspective, as recent studies have shown significant regrowth and reactivation of indicator organisms. Persistence of DNA in the environment after cell death in the range of days to weeks limits the application of DNA-based approaches as a measure of live cell density. Using selective nucleic acid intercalating dyes like ethidium monoazide (EMA) and propidium monoazide (PMA) is one of the alternative approaches to detecting and quantifying viable cells by quantitative PCR. These compounds have the ability to penetrate only into dead cells with compromised membrane integrity and intercalate with DNA via their photoinducible azide groups and in turn inhibit DNA amplification during PCRs. PMA has been successfully used in different studies and microorganisms, but it has not been evaluated sufficiently for complex environmental samples such as biosolids. In this study, experiments were performed with Escherichia coli ATCC 25922 as the model organism and the uidA gene as the target sequence using real-time PCR via the absolute quantification method. Experiments with the known quantities of live and dead cell mixtures showed that PMA treatment inhibits PCR amplification from dead cells with over 99% efficiency. The results also indicated that PMA-modified quantitative PCR could be successfully applied to biosolids when the total suspended solids (TSS) concentration is at or below 2,000 mg·liter(-1).     

Effect of PCR amplicon length on suppressing signals from membrane-compromised cells by propidium monoazide treatment

Treatment of microbiological samples with viability dyes prior to extraction of DNA and PCR amplification for downstream analysis has evolved into a commonly applied method. The addition of this easy-to-perform step to the sample analysis procedure inhibits the amplification of DNA from dead cells with compromised cell membranes. The method is currently used both in combination with quantitative PCR (qPCR), end-point PCR, and isothermal amplification. We present here a detailed study of the effect of amplicon size on amplification signals from unstressed and heat-exposed cells after treatment with propidium monoazide (PMA). PMA treatment was shown to be more efficient in excluding dead cells from the analysis both in combination with qPCR (PMA-qPCR) and denaturing gradient gel electrophoresis (PMA-DGGE), when longer amplicons were used. When applied to pure cultures of the fish pathogens Vibrio anguillarum and Flavobacterium psychrophilum exposed to a heat gradient ranging from mild to lethal, qPCR product lengths did not influence PMA-qPCR results at low temperatures, whereas an increasingly strong impact was seen at higher temperatures. Membrane permeability as a result of heat exposure might however have to be considered a conservative parameter for cell death for these pathogens as culturability and redox activity were lost at lower stress intensities than membrane integrity. When applying PMA-DGGE to an environmental water sample which was either left untreated or was exposed to heat, differences to non-PMA treated samples tended to slightly increase when amplified fragments in the first round of the nested PCR were longer, whereas the impact of 1st-round cycle numbers remains unclear.     

Method to detect only viable cells in microbial ecology

Propidium monoazide can limit the analysis of microbial communities derived from genetic fingerprints to viable cells with intact cell membranes. However, PMA treatment cannot completely suppress polymerase chain reaction (PCR) amplification when the targeted gene is too short. PMA treatment in combination with two-step nested PCR was designed to overcome this problem. Four experiments were performed to determine the limitation of PMA treatment and to evaluate the suitability of the method by applying the following samples: (1) pure cultures of Escherichia coli O157:H7, Enterobacter aerogenes, and Alcaligenes faecalis; (2) pond water samples spiked with heat-killed E. coli O157:H7 and E. aerogenes; (3) anaerobic sludge samples exposed to increasing heat stress; and (4) selected natural samples of estuarine sediment and lake mud. Results from the first two experiments show that PMA treatment cannot efficiently suppress dead cells from PCR amplification when the targeted gene is as short as 190 bp, however, the two-step nested PCR can overcome this problem. The last two experiments indicate the method that PMA treatment in combination with two-step nested PCR is useful for viable cells detection in microbial ecology.     

Molecular monitoring of disinfection efficacy using propidium monoazide in combination with quantitative PCR

One of the major drawbacks of DNA-based microbial diagnostics is its inability to discriminate between live and dead bacteria. Due to the persistence of DNA in the environment after cells have lost their viability, DNA-based assays cannot assess pathogenic risk since signals can originate from both live and dead cells. Presented here is a potential application of the novel chemical propidium monoazide (PMA), which results in the selective suppression of DNA detection from dead cells. PMA can only penetrate dead cells with permeabilized cell membranes. Upon intercalation into the DNA, covalent crosslinkage of PMA to DNA is achieved through light exposure. This modification prevents the DNA from being amplified by PCR. The method, in combination with quantitative PCR as a diagnostic tool, successfully monitored the disinfection efficacy of hypochlorite, benzalkonium and heat on several model pathogens. Threshold cycle numbers increased with increasing disinfection strength after PMA treatment of samples compared to non-PMA treated samples. With some disinfectant-specific differences, monitoring viability loss with membrane integrity as an indicator seemed to be more conservative than monitoring viability loss with plate counts. Loss of viability after short UV-exposure could not be monitored with PMA as UV light affects viability by inducing DNA damage without directly affecting membrane permeability.     

Molecular monitoring of inactivation efficiencies of bacteria during pulsed electric field treatment of clinical wastewater

Aims: The applicability of an alternative wastewater disinfection concept based on the pulsed electric field (PEF) treatment is tested with molecular biology techniques using clinical wastewaters. Methods and Results: Hospital wastewater was treated with the PEF technology. The inactivation efficiencies of bacteria were successfully monitored with real‐time polymerase chain reaction (PCR). As the differentiation between living and dead bacterial cells is important for the determination of the disinfection efficiency, propidium monoazide (PMA) was applied. PMA selectively penetrates cells with compromised membranes and intercalates into the DNA inhibiting a subsequent PCR amplification. The rates of reduction were examined for specific pathogens and wastewater populations using PCR‐denaturing gradient gel electrophoresis. The results showed that the main part of the bacterial population could be inactivated efficiently with the PEF treatment. Moreover, it was demonstrated that naturally occurring nuclease activities were not affected by the PEF treatment in contrast to a thermal treatment. Conclusions: The results indicated that the PEF treatment is an appropriate alternative disinfection concept for the treatment of clinical wastewaters and surpass the disadvantages of other disinfection methods. Significance and Impact of the Study: With the use of propidium monoazide for live–dead distinction, a new concept could be developed for the evaluation of disinfection methods.     

Distribution of bacterioplankton in meromictic Lake Saelenvannet, as determined by denaturing gradient gel electrophoresis of PCR-amplified gene fragments coding for 16S rRNA.

The community structure of bacterioplankton in meromictic Lake Saelenvannet was examined by PCR amplification of the V3 region of 16S rRNA from microbial communities recovered from various depths in the water column. Two different primer sets were used, one for amplification of DNA from the domain Bacteria and another specific for DNA from the domain Archaea. Amplified DNA fragments were resolved by denaturing gradient gel electrophoresis (DGGE), and the resulting profiles were reproducible and specific for the communities from different depths. Bacterial diversity estimated from the number and intensity of specific fragments in DGGE profiles decreased with depth. The reverse was true for the Archaea, with the diversity increasing with depth. Hybridization of DGGE profiles with oligonucleotide probes specific for phylogenetic groups of microorganisms showed the presence of both sulfate-reducing bacteria and methanogens throughout the water column, but they appeared to be most abundant below the chemocline. Several dominant fragments in the DGGE profiles were excised and sequenced. Among the dominant populations were representatives related to Chlorobium phaeovibrioides, chloroplasts from eukaryotic algae, and unidentified Archaea.     

Use of propidium monoazide and increased amplicon length reduce false-positive signals in quantitative PCR for bioburden analysis

Rapid microbiological methods (RMMs) as an alternative to conventional cultivation-based bioburden analysis are receiving increasing attention although no single technology is currently able to satisfy the needs of the health care industry. Among the RMMs, quantitative PCR (qPCR) seems particularly suited. Its implementation is, however, hampered by false-positive signals originating from free DNA in PCR reagents or from dead cells in the samples to be analysed. In this study, we assessed the capability of propidium monoazide (PMA) to inactivate exogenous DNA in PCR reagents and thus to minimise its impact in bioburden analysis. PMA is a membrane-impermeant dye that intercalates into DNA and covalently binds to it upon photoactivation leading to strong inhibition of PCR amplification. PMA is currently used mainly for treatment of microbiological samples to exclude signals from membrane-compromised cells, but is also very useful for suppression of exogenous DNA signals. In addition to testing the effect of different PMA concentrations on non-template controls and target DNA, we demonstrate the effect of amplicon length on the exclusion of background amplification. Targeting a 1,108-bp 16S rRNA gene fragment using universal bacterial primers and PCR reagents treated with 5 μM PMA resulted in complete suppression of signals from exogenous DNA within 50 cycles of amplification, while a limit of detection of 10 copies of Escherichia coli genomic DNA per PCR reaction was achieved. A combined PMA treatment of sample and PCR reagents furthermore improved the selective detection of live cells making this method appear a highly attractive RMM.     

Polymerase chain reaction amplification length-dependent ethidium monoazide suppression power for heat-killed cells of Enterobacteriaceae

The polymerase chain reaction (PCR) can confirm the presence of bacteria, but it is unable to differentiate between live and dead bacteria. Although ethidium monoazide (EMA)- and propidium monoazide (PMA)-based PCR have been evaluated, a quantity of ≥ 10(3)cells/ml dead cells produces a false-positive reading at 40 to 50 cycles (K. Rudi et al., Appl. Environ. Microbiol. 71 (2005) 1018-1024). After confirming the precision of real-time PCR of a long DNA target (16S or 23S ribosomal RNA [rRNA] gene, 1490 or 2840 bp), we evaluated the degree of suppression of an EMA treatment on the 16S/23S PCR using various amplification lengths (110-2840 bp) with heat-killed cells of Enterobacteriaceae (e.g., Salmonella enteritidis). We found that the inhibition rate was proportional to the PCR amplification length; short DNA (110 bp) amplification slightly delayed the threshold cycle (C(T)) of heat-killed cells of Enterobacteriaceae when compared with no EMA treatment. Regardless of the amplification length, the C(T) delay using live cells of Enterobacteriaceae with EMA was negligible. Thus, our real-time PCR of a long DNA (16S or 23S) template following EMA treatment is a rapid viable bacterial assay, which can potentially target all genera, for testing pasteurized milk that may have originally been contaminated with high levels of dead bacteria.     

应用DGGE法对青海相邻两盐湖中细菌多样性的快速检测

分别对青海相邻两盐湖柯柯盐湖、茶卡盐湖土样、泥样进行富集培养后,从中提取的DNA用两套不同的细菌通用引物进行扩增,分别得到包含V8和V9高变区的16SrDNA片断。经变性梯度凝胶电泳(DGGE)分析,结果显示这两个盐湖富集样品中细菌多样性具有较大的差异,而且同一盐湖不同性质样品中的细菌多样性差异也较大,两湖的泥样富集样品中均表现出了稍丰富的多样性。     

应用DGGE研究微生物群落时的常见问题分析

变性梯度凝胶电泳(DGGE)是通过核酸片段对微生物群落进行研究,可以监测未培养细菌及其功能基因,被广泛地应用于微生物群落多样性和动态分析,并成为微生物分子生态学研究中的重要手段之一。文中论述了DGGE操作过程中遇到的常见问题,并提出了相应的解决方法。全面分析了样品预处理过程和PCR扩增效果对DGGE分析的影响,探讨了DGGE图谱的优化过程和图谱分析方法,并对DGGE的应用前景进行了综述。     

Selective removal of DNA from dead cells of mixed bacterial communities by use of ethidium monoazide

The distinction between viable and dead bacterial cells poses a major challenge in microbial diagnostics. Due to the persistence of DNA in the environment after cells have lost viability, DNA-based quantification methods overestimate the number of viable cells in mixed populations or even lead to false-positive results in the absence of viable cells. On the other hand, RNA-based diagnostic methods, which circumvent this problem, are technically demanding and suffer from some drawbacks. A promising and easy-to-use alternative utilizing the DNA-intercalating dye ethidium monoazide bromide (EMA) was published recently. This chemical is known to penetrate only into "dead" cells with compromised cell membrane integrity. Subsequent photoinduced cross-linking was reported to inhibit PCR amplification of DNA from dead cells. We provide evidence here that in addition to inhibition of amplification, most of the DNA from dead cells is actually lost during the DNA extraction procedure, probably together with cell debris which goes into the pellet fraction. Exposure of bacteria to increasing stress and higher proportions of dead cells in defined populations led to increasing loss of genomic DNA. Experiments were performed using Escherichia coli O157:H7 and Salmonella enterica serovar Typhimurium as model pathogens and using real-time PCR for their quantification. Results showed that EMA treatment of mixed populations of these two species provides a valuable tool for selective removal of DNA of nonviable cells by using conventional extraction protocols. Furthermore, we provide evidence that prior to denaturing gradient gel electrophoresis, EMA treatment of a mature mixed-population drinking-water biofilm containing a substantial proportion of dead cells can result in community fingerprints dramatically different from those for an untreated biofilm. The interpretation of such fingerprints can have important implications in the field of microbial ecology.     

Selective detection of viable Helicobacter pylori using ethidium monoazide or propidium monoazide in combination with real-time polymerase chain reaction

Because Helicobacter pylori has a role in the pathogenesis of gastric cancer, chronic gastritis and peptic ulcer disease, detection of its viable form is very important. The objective of this study was to optimize a PCR method using ethidium monoazide (EMA) or propidium monoazide (PMA) for selective detection of viable H. pylori cells in mixed samples of viable and dead bacteria. Before conducting the real-time PCR using SodB primers of H. pylori, EMA or PMA was added to suspensions of viable and/or dead H. pylori cells at concentrations between 1 and 100 μM. PMA at a concentration of 50 μM induced the highest DNA loss in dead cells with little loss of genomic DNA in viable cells. In addition, selective detection of viable cells in the mixtures of viable and dead cells at various ratios was possible with the combined use of PMA and real-time PCR. In contrast, EMA penetrated the membranes of both viable and dead cells and induced degradation of their genomic DNA. The findings of this study suggest that PMA, but not EMA, can be used effectively to differentiate viable H. pylori from its dead form.     

Numerical Analysis of Grassland Bacterial Community Structure under Different Land Management Regimens by Using 16S Ribosomal DNA Sequence Data and Denaturing Gradient Gel Electrophoresis Banding Patterns

Bacterial diversity in unimproved and improved grassland soils was assessed by PCR amplification of bacterial 16S ribosomal DNA (rDNA) from directly extracted soil DNA, followed by sequencing of ~45 16S rDNA clones from each of three unimproved and three improved grassland samples (A. E. McCaig, L. A. Glover, and J. I. Prosser, Appl. Environ. Microbiol. 65:1721–1730, 1999) or by denaturing gradient gel electrophoresis (DGGE) of total amplification products. Semi-improved grassland soils were analyzed only by DGGE. No differences between communities were detected by calculation of diversity indices and similarity coefficients for clone data (possibly due to poor coverage). Differences were not observed between the diversities of individual unimproved and improved grassland DGGE profiles, although considerable spatial variation was observed among triplicate samples. Semi-improved grassland samples, however, were less diverse than the other grassland samples and had much lower within-group variation. DGGE banding profiles obtained from triplicate samples pooled prior to analysis indicated that there was less evenness in improved soils, suggesting that selection for specific bacterial groups occurred. Analysis of DGGE profiles by canonical variate analysis but not by principal-coordinate analysis, using unweighted data (considering only the presence and absence of bands) and weighted data (considering the relative intensity of each band), demonstrated that there were clear differences between grasslands, and the results were not affected by weighting of data. This study demonstrated that quantitative analysis of data obtained by community profiling methods, such as DGGE, can reveal differences between complex microbial communities.     

Numerical analysis of grassland bacterial community structure under different land management regimens by using 16S ribosomal DNA sequence data and denaturing gradient gel electrophoresis banding patterns

Bacterial diversity in unimproved and improved grassland soils was assessed by PCR amplification of bacterial 16S ribosomal DNA (rDNA) from directly extracted soil DNA, followed by sequencing of ~45 16S rDNA clones from each of three unimproved and three improved grassland samples (A. E. McCaig, L. A. Glover, and J. I. Prosser, Appl. Environ. Microbiol. 65:1721-1730, 1999) or by denaturing gradient gel electrophoresis (DGGE) of total amplification products. Semi-improved grassland soils were analyzed only by DGGE. No differences between communities were detected by calculation of diversity indices and similarity coefficients for clone data (possibly due to poor coverage). Differences were not observed between the diversities of individual unimproved and improved grassland DGGE profiles, although considerable spatial variation was observed among triplicate samples. Semi-improved grassland samples, however, were less diverse than the other grassland samples and had much lower within-group variation. DGGE banding profiles obtained from triplicate samples pooled prior to analysis indicated that there was less evenness in improved soils, suggesting that selection for specific bacterial groups occurred. Analysis of DGGE profiles by canonical variate analysis but not by principal-coordinate analysis, using unweighted data (considering only the presence and absence of bands) and weighted data (considering the relative intensity of each band), demonstrated that there were clear differences between grasslands, and the results were not affected by weighting of data. This study demonstrated that quantitative analysis of data obtained by community profiling methods, such as DGGE, can reveal differences between complex microbial communities.     

Quantifying fungal viability in air and water samples using quantitative PCR after treatment with propidium monoazide (PMA)

A method is described to discriminate between live and dead cells of the infectious fungi Aspergillus fumigatus, Aspergillus flavus, Aspergillus terreus, Mucor racemosus, Rhizopus stolonifer and Paecilomyces variotii. To test the method, conidial suspensions were heat inactivated at 85 degrees C or held at 5 degrees C (controls) for 1 h. Polycarbonate filters (25 mm diameter, 0.8 microm pore size) were placed on "welled" slides (14 mm diameter) and the filters treated with either PBS or PMA. Propidium monoazide (PMA), which enters dead cells but not live cells, was incubated with cell suspensions, exposed to blue wavelength light-emitting diodes (LED) to inactivate remaining PMA and secure intercalation of PMA with DNA of dead cells. Treated cells were extracted and the live and dead cells evaluated with quantitative PCR (QPCR). After heat treatment and DNA modification with PMA, all fungal species tested showed an approximate 100- to 1000-fold difference in cell viability estimated by QPCR analysis which was consistent with estimates of viability based on culturing.     

Differentiation of genes extracted from non-viable versus viable micro-organisms in environmental samples using ethidium monoazide bromide

Differentiation of DNA derived from viable or non-viable microorganisms within mixed microbial communities continues to be one of the greatest challenges in molecular studies of environmental samples. A novel method developed for microbial food pathogens is tested here on environmental samples. This technique involves the use of ethidium monoazide bromide (EMA) for the distinction of live/dead cells. In non-viable cells EMA intercalates into the DNA which prevents amplification by PCR. We adapted and evaluated the EMA technique for soil, elemental sulfur and river biofilm samples. Quantitative PCR determined that EMA suppressed 99.99% of E. coli LKI gfp+ signal in non-viable cultures and 100.00% when the cultures were added to soil samples. The same technique was also successful at suppressing DNA amplification from spiked non-viable cells in elemental sulfur samples by 100.00%, but not in three Saskatchewan River biofilms. In sub Antarctic soil, EMA-Q-PCR was used to detect the prevalence of a functional gene, amoA, and this was closely correlated to nitrification activity measurements. The ability of EMA to differentiate between viable and non-viable populations in soil was confirmed by the similarity of the 16S rRNA denaturing-gradient-gel electrophoresis DNA fingerprint of EMA treated soil and the 16S rRNA cDNA fingerprint of non-EMA treated soil. The EMA technique effectively suppressed amplification of non-viable spiked controls, closely mirrored activity assays and yielded community composition profiles similar to rRNA techniques. The use of EMA in soil effectively suppressed amplification of non-viable organism DNA, however it was not effective in biofilm samples and EMA partially inhibited amplification of viable organism DNA in elemental sulfur samples.     

Variable selection in high-dimensional multivariate binary data with application to the analysis of microbial community DNA fingerprints

In order to understand the relevance of microbial communities on crop productivity, the identification and characterization of the rhizosphere soil microbial community is necessary. Characteristic profiles of the microbial communities are obtained by denaturing gradient gel electrophoresis (DGGE) of polymerase chain reaction (PCR) amplified 16S rDNA from soil extracted DNA. These characteristic profiles, commonly called community DNA fingerprints, can be represented in the form of high-dimensional binary vectors. We address the problem of modeling and variable selection in high-dimensional multivariate binary data and present an application of our methodology in the context of a controlled agricultural experiment.