The necrotrophic effector ToxA from Parastagonospora nodorum interacts with wheat NHL proteins to facilitate Tsn1-mediated necrosis.

The plant pathogen Parastagonospora nodorum secretes necrotrophic effectors to promote disease. These effectors induce cell death on wheat cultivars carrying dominant susceptibility genes in an inverse gene-for-gene manner. However, the molecular mechanisms underpinning these interactions and resulting cell death remain unclear. Here, we used a yeast two-hybrid library approach to identify wheat proteins that interact with the necrotrophic effector ToxA. Using this strategy, we identified an interaction between ToxA and a wheat transmembrane NDR/HIN1-like protein (TaNHL10) and confirmed the interaction using in planta co-immunoprecipitation and confocal microscopy co-localization analysis. We showed that the C-terminus of TaNHL10 is extracellular whilst the N-terminus is localized in the cytoplasm. Further analyses using yeast two-hybrid and confocal microscopy co-localization showed that ToxA interacts with the C-terminal LEA2 extracellular domain of TaNHL10. Random mutagenesis was then used to identify a ToxA mutant, ToxAN109D , which was unable to interact with TaNHL10 in yeast two-hybrid assays. Subsequent heterologous expression and purification of ToxAN109D in Nicotiania benthamiana revealed that the mutated protein was unable to induce necrosis on Tsn1-dominant wheat cultivars, confirming that the interaction of ToxA with TaNHL10 is required to induce cell death. Collectively, these data advance our understanding on how ToxA induces cell death during infection and further highlight the importance of host cell surface interactions in necrotrophic pathosystems.

The Necrotrophic Pathogen Parastagonospora nodorum Is a Master Manipulator of Wheat Defense.

Parastagonospora nodorum is a necrotrophic pathogen of wheat that is particularly destructive in major wheat-growing regions of the United States, northern Europe, Australia, and South America. P. nodorum secretes necrotrophic effectors that target wheat susceptibility genes to induce programmed cell death (PCD), resulting in increased colonization of host tissue and, ultimately, sporulation to complete its pathogenic life cycle. Intensive research over the last two decades has led to the functional characterization of five proteinaceous necrotrophic effectors, SnTox1, SnToxA, SnTox267, SnTox3, and SnTox5, and three wheat susceptibility genes, Tsn1, Snn1, and Snn3D-1. Functional characterization has revealed that these effectors, in addition to inducing PCD, have additional roles in pathogenesis, including chitin binding that results in protection from wheat chitinases, blocking defense response signaling, and facilitating plant colonization. There are still large gaps in our understanding of how this necrotrophic pathogen is successfully manipulating wheat defense to complete its life cycle. This review summarizes our current knowledge, identifies knowledge gaps, and provides a summary of well-developed tools and resources currently available to study the P. nodorum-wheat interaction, which has become a model for necrotrophic specialist interactions. Further functional characterization of the effectors involved in this interaction and work toward a complete understanding of how P. nodorum manipulates wheat defense will provide fundamental knowledge about this and other necrotrophic interactions. Additionally, a broader understanding of this interaction will contribute to the successful management of Septoria nodorum blotch disease on wheat. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A Revised Nomenclature for ToxA Haplotypes Across Multiple Fungal Species.

ToxA is one of the most studied proteinaceous necrotrophic effectors produced by plant pathogens. It has been identified in four pathogens (Pyrenophora tritici-repentis, Parastagonospora nodorum, Parastagonospora pseudonodorum [formerly Parastagonospora avenaria f. sp. tritici], and Bipolaris sorokiniana) causing leaf spot diseases on cereals worldwide. To date, 24 different ToxA haplotypes have been identified. Some P. tritici-repentis and related species also express ToxB, another small protein necrotrophic effector. We present here a revised and standardized nomenclature for these effectors, which could be extended to other poly-haplotypic genes found across multiple species.

Victorin, the host-selective cyclic peptide toxin from the oat pathogen Cochliobolus victoriae, is ribosomally encoded.

The necrotrophic fungal pathogen Cochliobolus victoriae produces victorin, a host-selective toxin (HST) essential for pathogenicity to certain oat cultivars with resistance against crown rust. Victorin is a mixture of highly modified heterodetic cyclic hexapeptides, previously assumed to be synthesized by a nonribosomal peptide synthetase. Herein, we demonstrate that victorin is a member of the ribosomally synthesized and posttranslationally modified peptide (RiPP) family of natural products. Analysis of a newly generated long-read assembly of the C. victoriae genome revealed three copies of precursor peptide genes (vicA1-3) with variable numbers of "GLKLAF" core peptide repeats corresponding to the victorin peptide backbone. vicA1-3 are located in repeat-rich gene-sparse regions of the genome and are loosely clustered with putative victorin biosynthetic genes, which are supported by the discovery of compact gene clusters harboring corresponding homologs in two distantly related plant-associated Sordariomycete fungi. Deletion of at least one copy of vicA resulted in strongly diminished victorin production. Deletion of a gene encoding a DUF3328 protein (VicYb) abolished the production altogether, supporting its predicted role in oxidative cyclization of the core peptide. In addition, we uncovered a copper amine oxidase (CAO) encoded by vicK, in which its deletion led to the accumulation of new glycine-containing victorin derivatives. The role of VicK in oxidative deamination of the N-terminal glycyl moiety of the hexapeptides to the active glyoxylate forms was confirmed in vitro. This study finally unraveled the genetic and molecular bases for biosynthesis of one of the first discovered HSTs and expanded our understanding of underexplored fungal RiPPs.

Optimized Production of Disulfide-Bonded Fungal Effectors in Escherichia coli Using CyDisCo and FunCyDisCo Coexpression Approaches.

Effectors are a key part of the arsenal of plant-pathogenic fungi and promote pathogen virulence and disease. Effectors typically lack sequence similarity to proteins with known functional domains and motifs, limiting our ability to predict their functions and understand how they are recognized by plant hosts. As a result, cross-disciplinary approaches involving structural biology and protein biochemistry are often required to decipher and better characterize effector function. These approaches are reliant on high yields of relatively pure protein, which often requires protein production using a heterologous expression system. For some effectors, establishing an efficient production system can be difficult, particularly those that require multiple disulfide bonds to achieve their naturally folded structure. Here, we describe the use of a coexpression system within the heterologous host Escherichia coli, termed CyDisCo (cytoplasmic disulfide bond formation in E. coli) to produce disulfide bonded fungal effectors. We demonstrate that CyDisCo and a naturalized coexpression approach termed FunCyDisCo (Fungi CyDisCo) can significantly improve the production yields of numerous disulfide-bonded effectors from diverse fungal pathogens. The ability to produce large quantities of functional recombinant protein has facilitated functional studies and crystallization of several of these reported fungal effectors. We suggest this approach could be broadly useful in the investigation of the function and recognition of a broad range of disulfide bond-containing effectors.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

PR1-mediated defence via C-terminal peptide release is targeted by a fungal pathogen effector.

The effector SnTox3 from Parastagonospora nodorum elicits a strong necrotic response in susceptible wheat and also interacts with wheat pathogenesis-related protein 1 (TaPR-1), although the function of this interaction in disease is unclear. Here, we dissect TaPR1 function by studying SnTox3-TaPR1 interaction and demonstrate the dual functionality of SnTox3. We utilized site-directed mutagenesis to identify an SnTox3 variant, SnTox3P173S , that was unable to interact with TaPR1 in yeast-two-hybrid assays. Additionally, using recombinant proteins we established a novel protein-mediated phenotyping assay allowing functional studies to be undertaken in wheat. Wheat leaves infiltrated with TaPR1 proteins showed significantly less disease compared to control leaves, correlating with a strong increase in defence gene expression. This activity was dependent on release of the TaCAPE1 peptide embedded within TaPR1 by an unidentified serine protease. The priming activity of TaPR1 was compromised by SnTox3 but not the noninteracting variant SnTox3P173S , and we demonstrate that SnTox3 prevents TaCAPE1 release from TaPR1 in vitro. SnTox3 independently functions to induce necrosis through recognition by Snn3 and also suppresses host defence through a direct interaction with TaPR1 proteins. Importantly, this study also advances our understanding of the role of PR1 proteins in host-microbe interactions as inducers of host defence signalling.

The crystal structure of SnTox3 from the necrotrophic fungus Parastagonospora nodorum reveals a unique effector fold and provides insight into Snn3 recognition and pro-domain protease processing of fungal effectors.

Plant pathogens cause disease through secreted effector proteins, which act to promote infection. Typically, the sequences of effectors provide little functional information and further targeted experimentation is required. Here, we utilized a structure/function approach to study SnTox3, an effector from the necrotrophic fungal pathogen Parastagonospora nodorum, which causes cell death in wheat-lines carrying the sensitivity gene Snn3. We developed a workflow for the production of SnTox3 in a heterologous host that enabled crystal structure determination and functional studies. We show this approach can be successfully applied to study effectors from other pathogenic fungi. The ß-barrel fold of SnTox3 is a novel fold among fungal effectors. Structure-guided mutagenesis enabled the identification of residues required for Snn3 recognition. SnTox3 is a pre-pro-protein, and the pro-domain of SnTox3 can be cleaved in vitro by the protease Kex2. Complementing this, an in silico study uncovered the prevalence of a conserved motif (LxxR) in an expanded set of putative pro-domain-containing fungal effectors, some of which can be cleaved by Kex2 in vitro. Our in vitro and in silico study suggests that Kex2-processed pro-domain (designated here as K2PP) effectors are common in fungi and this may have broad implications for the approaches used to study their functions.

The global regulator of pathogenesis PnCon7 positively regulates Tox3 effector gene expression through direct interaction in the wheat pathogen Parastagonospora nodorum.

To investigate effector gene regulation in the wheat pathogenic fungus Parastagonospora nodorum, the promoter and expression of Tox3 was characterised through a series of complementary approaches. Promoter deletion and DNase I footprinting experiments identified a 25 bp region in the Tox3 promoter as being required for transcription. Subsequent yeast one-hybrid analysis using the DNA sequence as bait identified that interacting partner as the C2H2 zinc finger transcription factor PnCon7, a putative master regulator of pathogenesis. Silencing of PnCon7 resulted in the down-regulation of Tox3 demonstrating that the transcription factor has a positive regulatory role on gene expression. Analysis of Tox3 expression in the PnCon7 silenced strains revealed a strong correlation with PnCon7 transcript levels, supportive of a direct regulatory role. Subsequent pathogenicity assays using PnCon7-silenced isolates revealed that the transcription factor was required for Tox3-mediated disease. The expression of two other necrotrophic effectors (ToxA and Tox1) was also affected but in a non-dose dependent manner suggesting that the regulatory role of PnCon7 on these genes was indirect. Collectively, these data have advanced our fundamental understanding of the Con7 master regulator of pathogenesis by demonstrating its positive regulatory role on the Tox3 effector in P. nodorum through direct interaction.

Rapid Parallel Evolution of Azole Fungicide Resistance in Australian Populations of the Wheat Pathogen Zymoseptoria tritici.

Zymoseptoria tritici is a globally distributed fungal pathogen which causes Septoria tritici blotch on wheat. Management of the disease is attempted through the deployment of resistant wheat cultivars and the application of fungicides. However, fungicide resistance is commonly observed in Z. tritici populations, and continuous monitoring is required to detect breakdowns in fungicide efficacy. We recently reported azole-resistant isolates in Australia; however, it remained unknown whether resistance was brought into the continent through gene flow or whether resistance emerged independently. To address this question, we screened 43 isolates across five Australian locations for azole sensitivity and performed whole-genome sequencing on 58 isolates from seven locations to determine the genetic basis of resistance. Population genomic analyses showed extremely strong differentiation between the Australian population recovered after azoles began to be used and both Australian populations recovered before azoles began to be used and populations on different continents. The apparent absence of recent gene flow between Australia and other continents suggests that azole fungicide resistance has evolved de novo and subsequently spread within Tasmania. Despite the isolates being distinct at the whole-genome level, we observed combinations of nonsynonymous substitutions at the CYP51 locus identical to those observed elsewhere in the world. We observed nine previously reported nonsynonymous mutations as well as isolates that carried a combination of the previously reported L50S, S188N, A379G, I381V, Y459DEL, G460DEL, and N513K substitutions. Assays for the 50% effective concentration against a subset of isolates exposed to the tebuconazole and epoxiconazole fungicides showed high levels of azole resistance. The rapid, parallel evolution of a complex CYP51 haplotype that matches a dominant European haplotype demonstrates the enormous potential for de novo resistance emergence in pathogenic fungi.IMPORTANCE Fungicides are essential to control diseases in agriculture because many crops are highly susceptible to pathogens. However, many pathogens rapidly evolve resistance to fungicides. A large body of studies have described specific mutations conferring resistance and have often made inferences about the origins of resistance based on sequencing data from the target gene alone. Here, we show the de novo acquisition of resistance to the ubiquitously used azole fungicides in genetically isolated populations of the wheat pathogen Zymoseptoria tritici in Tasmania, Australia. We confirm evidence for parallel evolution through genome-scale analyses of representative worldwide populations. The emergence of complex resistance haplotypes following a well-documented recent introduction of azoles into Australian farming practices demonstrates how rapidly chemical resistance evolves in agricultural ecosystems.

Biosynthesis of a Tricyclo[6.2.2.02,7 ]dodecane System by a Berberine Bridge Enzyme-Like Aldolase.

The aldol reaction is one of the most fundamental stereocontrolled carbon-carbon bond-forming reactions and is mainly catalyzed by aldolases in nature. Despite the fact that the aldol reaction has been widely proposed to be involved in fungal secondary metabolite biosynthesis, a dedicated aldolase that catalyzes stereoselective aldol reactions has only rarely been reported in fungi. Herein, we activated a cryptic polyketide biosynthetic gene cluster that was upregulated in the fungal wheat pathogen Parastagonospora nodorum during plant infection; this resulted in the production of the phytotoxic stemphyloxin II (1). Through heterologous reconstruction of the biosynthetic pathway and in vitro assay by using cell-free lysate from Aspergillus nidulans, we demonstrated that a berberine bridge enzyme (BBE)-like protein SthB catalyzes an intramolecular aldol reaction to establish the bridged tricyclo[6.2.2.02,7 ]dodecane skeleton in the post-assembly tailoring step. The characterization of SthB as an aldolase enriches the catalytic toolbox of classic reactions and the functional diversities of the BBE superfamily of enzymes.

Bipolenins K-N: New sesquiterpenoids from the fungal plant pathogen Bipolaris sorokiniana.

Chemical investigation of the barley and wheat fungal pathogen Bipolaris sorokiniana BRIP10943 yielded four new sativene-type sesquiterpenoid natural products, bipolenins K-N (1-4), together with seven related known analogues (5-11), and a sesterterpenoid (12). Their structures were determined by detailed analysis of spectroscopic data, supported by TDDFT calculations and comparison with previously reported analogues. These compounds were evaluated for their phytotoxic activity against wheat seedlings and wheat seed germination. The putative biosynthetic relationships between the isolated sesquiterpenoids were also explored.

Genetic mapping of Stb19, a new resistance gene to Zymoseptoria tritici in wheat.

KEY MESSAGE: A new and dominant R gene Stb19 is identified from a soft wheat cultivar 'Lorikeet' and was mapped on the distal region of chromosome 1DS. Two tightly linked KASP markers were also discovered and validated for molecular-assisted breeding programs. A new R gene, designated as Stb19, provides resistance to Zymoseptoria tritici in wheat. This new dominant gene resides on the short arm of chromosome 1D, exhibiting complete resistance to three Z. tritici isolates, WAI332, WAI251, and WAI161, at the seedling stage. A genetic linkage map, based on an F2:3 population of 'Lorikeet' and 'Summit,' found the Stb19 gene at a 9.3 cM region on 1DS, closely linked with two Kompetitive Allele-Specific PCR markers, snp_4909967 and snp_1218021. Further, the two markers were tested and validated in another F2:3 population and 266 different wheat accessions, which gave over 95% accuracy of resistance/susceptibility prediction. Combined with the physical location of the identified SNPs and the previous evidence of gene order on chromosome 1DS (centromere-Sr45-Sr33-Lr21-telomere), Stb19 is proposed to be located between Sr33 and Lr21. Thus, the newly discovered Stb19 along with the KASP markers represents an increase in genetic resources available for wheat breeding resistance to Z. tritici.

Wheat PR-1 proteins are targeted by necrotrophic pathogen effector proteins.

Recent studies have identified that proteinaceous effectors secreted by Parastagonospora nodorum are required to cause disease on wheat. These effectors interact in a gene-for-gene manner with host-dominant susceptibilty loci, resulting in disease. However, whilst the requirement of these effectors for infection is clear, their mechanisms of action remain poorly understood. A yeast-two-hybrid library approach was used to search for wheat proteins that interacted with the necrotrophic effector SnTox3. Using this strategy we indentified an interaction between SnTox3 and the wheat pathogenicity-related protein TaPR-1-1, and confirmed it by in-planta co-immunprecipitation. PR-1 proteins represent a large family (23 in wheat) of proteins that are upregulated early in the defence response; however, their function remains ellusive. Interestingly, the P. nodorum effector SnToxA has recently been shown to interact specifically with TaPR-1-5. Our analysis of the SnTox3-TaPR-1 interaction demonstrated that SnTox3 can interact with a broader range of TaPR-1 proteins. Based on these data we utilised homology modeling to predict, and validate, regions on TaPR-1 proteins that are likely to be involved in the SnTox3 interaction. Precipitating from this work, we identified that a PR-1-derived defence signalling peptide from the C-terminus of TaPR-1-1, known as CAPE1, enhanced the infection of wheat by P. nodorum in an SnTox3-dependent manner, but played no role in ToxA-mediated disease. Collectively, our data suggest that P. nodorum has evolved unique effectors that target a common host-protein involved in host defence, albeit with different mechanisms and potentially outcomes.

Differential effector gene expression underpins epistasis in a plant fungal disease.

Fungal effector-host sensitivity gene interactions play a key role in determining the outcome of septoria nodorum blotch disease (SNB) caused by Parastagonospora nodorum on wheat. The pathosystem is complex and mediated by interaction of multiple fungal necrotrophic effector-host sensitivity gene systems. Three effector sensitivity gene systems are well characterized in this pathosystem; SnToxA-Tsn1, SnTox1-Snn1 and SnTox3-Snn3. We tested a wheat mapping population that segregated for Snn1 and Snn3 with SN15, an aggressive P. nodorum isolate that produces SnToxA, SnTox1 and SnTox3, to study the inheritance of sensitivity to SnTox1 and SnTox3 and disease susceptibility. Interval quantitative trait locus (QTL) mapping showed that the SnTox1-Snn1 interaction was paramount in SNB development on both seedlings and adult plants. No effect of the SnTox3-Snn3 interaction was observed under SN15 infection. The SnTox3-Snn3 interaction was however, detected in a strain of SN15 in which SnTox1 had been deleted (tox1-6). Gene expression analysis indicates increased SnTox3 expression in tox1-6 compared with SN15. This indicates that the failure to detect the SnTox3-Snn3 interaction in SN15 is due - at least in part - to suppressed expression of SnTox3 mediated by SnTox1. Furthermore, infection of the mapping population with a strain deleted in SnToxA, SnTox1 and SnTox3 (toxa13) unmasked a significant SNB QTL on 2DS where the SnTox2 effector sensitivity gene, Snn2, is located. This QTL was not observed in SN15 and tox1-6 infections and thus suggesting that SnToxA and/or SnTox3 were epistatic. Additional QTLs responding to SNB and effectors sensitivity were detected on 2AS1 and 3AL.

Functional genomics-guided discovery of a light-activated phytotoxin in the wheat pathogen Parastagonospora nodorum via pathway activation.

Parastagonospora nodorum is an important pathogen of wheat. The contribution of secondary metabolites to this pathosystem is poorly understood. A biosynthetic gene cluster (SNOG_08608-08616) has been shown to be upregulated during the late stage of P. nodorum wheat leaf infection. The gene cluster shares several homologues with the Cercospora nicotianae CTB gene cluster encoding the biosynthesis of cercosporin. Activation of the gene cluster by overexpression (OE) of the transcription factor gene (SNOG_08609) in P. nodorum resulted in the production of elsinochrome C, a perelyenequinone phytotoxin structurally similar to cercosporin. Heterologous expression of the polyketide synthase gene elcA from the gene cluster in Aspergillus nidulans resulted in the production of the polyketide precursor nortoralactone common to the cercosporin pathway. Elsinochrome C could be detected on wheat leaves infected with P. nodorum, but not in the elcA disruption mutant. The compound was shown to exhibit necrotic activity on wheat leaves in a light-dependent manner. Wheat seedling infection assays showed that ΔelcA exhibited reduced virulence compared with wild type, while infection by an OE strain overproducing elsinochrome C resulted in larger lesions on leaves. These data provided evidence that elsinochrome C contributes to the virulence of P. nodorum against wheat.

Transition from heterothallism to homothallism is hypothesised to have facilitated speciation among emerging Botryosphaeriaceae wheat-pathogens.

White grain disorder (WGD) is a recently emerged wheat disease in Australia caused by three Botryosphaeriaceae fungi, from the genus Eutiarosporella. These species are E. tritici-australis, E. darliae, and E. pseudodarliae. Characterisation of the mating type genes for the WGD-species show that the genome sequence of a single E. darliae and E. pseudodarliae isolate both harbour MAT1-2-1 and MAT1-1-1, which suggests that these species are homothallic. However, unlike most other characterised mating-type loci from other homothallic Dothideomycetes, these species' MAT1-1-1 are located at a separate locus, inserted within the coding region of another gene. The sequenced strain of E. tritici-australis analysed did not harbour MAT1-1-1. Including the sequenced strain, we screened the mating type genes present in 16 E. tritici-australis individuals isolated from infected grain from fields in South Australia. Of these 16, 11 harbour MAT1-1-1 and the other five harbour MAT1-2-1. The genome of a MAT1-1-1 harbouring isolate was re-sequenced, which demonstrated that MAT1-1-1 was present at the MAT locus. We examined non-coding DNA surrounding the MAT1-1-1 gene in E. pseudodarliae and observed fragments of the MAT locus both up and downstream. These fragments and their orientation around MAT1-1-1 is similar to characterised heterothallic Botryosphaeriaceae. Based on these gene arrangements, we conclude that the new MAT1-1-1 containing locus likely originated from a cryptic DNA integration event between two heterothallic individuals. We hypothesise that this integration event led to the formation of a homothallic lineage, which is the common ancestor of E. darliae and E. pseudodarliae.

Heterologous expression of cytotoxic sesquiterpenoids from the medicinal mushroom Lignosus rhinocerotis in yeast.

BACKGROUND: Genome mining facilitated by heterologous systems is an emerging approach to access the chemical diversity encoded in basidiomycete genomes. In this study, three sesquiterpene synthase genes, GME3634, GME3638, and GME9210, which were highly expressed in the sclerotium of the medicinal mushroom Lignosus rhinocerotis, were cloned and heterologously expressed in a yeast system. RESULTS: Metabolite profile analysis of the yeast culture extracts by GC-MS showed the production of several sesquiterpene alcohols (C15H26O), including cadinols and germacrene D-4-ol as major products. Other detected sesquiterpenes include selina-6-en-4-ol, ß-elemene, ß-cubebene, and cedrene. Two purified major compounds namely (+)-torreyol and α-cadinol synthesised by GME3638 and GME3634 respectively, are stereoisomers and their chemical structures were confirmed by 1H and 13C NMR. Phylogenetic analysis revealed that GME3638 and GME3634 are a pair of orthologues, and are grouped together with terpene synthases that synthesise cadinenes and related sesquiterpenes. (+)-Torreyol and α-cadinol were tested against a panel of human cancer cell lines and the latter was found to exhibit selective potent cytotoxicity in breast adenocarcinoma cells (MCF7) with IC50 value of 3.5 ± 0.58 µg/ml while α-cadinol is less active (IC50 = 18.0 ± 3.27 µg/ml). CONCLUSIONS: This demonstrates that yeast-based genome mining, guided by transcriptomics, is a promising approach for uncovering bioactive compounds from medicinal mushrooms.

Combining protein ratio p-values as a pragmatic approach to the analysis of multirun iTRAQ experiments.

iTRAQ labeling of peptides is widely used for quantitative comparison of biological samples using mass spectrometry. However, iTRAQ determined protein ratios have varying credibility depending on the number and quality of the peptide ratios used to generate them, and accounting for this becomes problematic particularly in the multirun scenario needed for larger scale biological studies. One approach to this problem relies on the use of sophisticated statistical global models using peptide ratios rather than working directly with the protein ratios, but these yield complex models whose solution relies on computational approaches such as stage-wise regression, which are nontrivial to run and verify. Here we evaluate an alternative pragmatic approach to finding differentially expressed proteins based on combining protein ratio p-values across experiments in a fashion similar to running a meta-analysis across different iTRAQ runs. Our approach uses the well-established Stouffer's Z-transform for combining p-values, alongside a ratio trend consistency measure, which we introduce. We evaluate this method with data from two iTRAQ experiments using plant and animal models. We show that in the specific context of iTRAQ data analysis this method has advantages of simplicity, high tolerance of run variability, low false discovery rate, and emphasis on proteins identified with high confidence.