pISA-tree - a data management framework for life science research projects using a standardised directory tree.

We developed pISA-tree, a straightforward and flexible data management solution for organisation of life science project-associated research data and metadata. pISA-tree was initiated by end-user requirements thus its strong points are practicality and low maintenance cost. It enables on-the-fly creation of enriched directory tree structure (project/Investigation/Study/Assay) based on the ISA model, in a standardised manner via consecutive batch files. Templates-based metadata is generated in parallel at each level enabling guided submission of experiment metadata. pISA-tree is complemented by two R packages, pisar and seekr. pisar facilitates integration of pISA-tree datasets into bioinformatic pipelines and generation of ISA-Tab exports. seekr enables synchronisation with the FAIRDOMHub repository. Applicability of pISA-tree was demonstrated in several national and international multi-partner projects. The system thus supports findable, accessible, interoperable and reusable (FAIR) research and is in accordance with the Open Science initiative. Source code and documentation of pISA-tree are available at https://github.com/NIB-SI/pISA-tree .

Hay meadow vibroscape and interactions within insect vibrational community.

Our experiences shape our knowledge and understanding of the world around us. The natural vibrational environment (vibroscape) is hidden to human senses but is nevertheless perceived and exploited by the majority of animals. Here, we show that the vibroscape recorded on plants in a temperate hay meadow is a dynamic low-frequency world, rich in species-specific vibrational signals. The overall vibroscape composition changed throughout the season and also depended on the plant species, as well as on the spatial position of individual plants within the meadow. Within the studied community, vibrationally signaling species sharing this communication channel avoided interference primarily by partitioning vibrational space on a fine temporal scale. The vibroscape is a reliable source of information in the environment and expands our understanding of ecological and evolutionary processes.

Precision transcriptomics of viral foci reveals the spatial regulation of immune-signaling genes and identifies RBOHD as an important player in the incompatible interaction between potato virus Y and potato.

Whereas the activation of resistance (R) proteins has been intensively studied, the downstream signaling mechanisms leading to the restriction of the pathogen remain mostly unknown. We studied the immunity network response conditioned by the potato Ny-1 gene against potato virus Y. We analyzed the processes in the cell death zone and surrounding tissue on the biochemical and gene expression levels in order to reveal the spatiotemporal regulation of the immune response. We show that the transcriptional response in the cell death zone and surrounding tissue is dependent on salicylic acid (SA). For some genes the spatiotemporal regulation is completely lost in the SA-deficient line, whereas other genes show a different response, indicating multiple connections between hormonal signaling modules. The induction of NADPH oxidase RBOHD expression occurs specifically on the lesion border during the resistance response. In plants with silenced RBOHD, the functionality of the resistance response is perturbed and the spread of the virus is not arrested at the site of infection. RBOHD is required for the spatial accumulation of SA, and conversely RBOHD is under the transcriptional regulation of SA. Using spatially resolved RNA-seq, we also identified spatial regulation of an UDP-glucosyltransferase, another component in feedback activation of SA biosynthesis, thus deciphering a novel aspect of resistance signaling.

Localization patterns of cathepsins K and X and their predictive value in glioblastoma.

Background Glioblastoma is a highly aggressive central nervous system neoplasm characterized by extensive infiltration of malignant cells into brain parenchyma, thus preventing complete tumor eradication. Cysteine cathepsins B, S, L and K are involved in cancer progression and are overexpressed in glioblastoma. We report here for the first time that cathepsin X mRNA and protein are also abundantly present in malignant glioma. Materials and methods Gene expression of cathepsins K and X was analyzed using publically-available tran-scriptomic datasets and correlated with glioma grade and glioblastoma subtype. Kaplan-Maier survival analysis was performed to evaluate the predictive value of cathepsin K and X mRNA expression. Cathepsin protein expression was localized and semi-quantified in tumor tissues by immunohistochemistry. Results Highest gene expression of cathepsins K and X was found in glioblastoma, in particular in the mesenchymal subtype. Overall, high mRNA expression of cathepsin X, but not that of cathepsin K, correlated with poor patients' survival. Cathepsin K and X proteins were abundantly and heterogeneously expressed in glioblastoma tissue. Immuno-labeling of cathepsins K and X was observed in areas of CD133-positive glioblastoma stem cells, localized around arterioles in their niches that also expressed SDF-1α and CD68. mRNA levels of both cathepsins K and X correlated with mRNA levels of markers of glioblastoma stem cells and their niches. Conclusions The presence of both cathepsins in glioblastoma stem cell niche regions indicates their possible role in regulation of glioblastoma stem cell homing in their niches. The clinical relevance of this data needs to be elaborated in further prospective studies.

DiNAR: revealing hidden patterns of plant signalling dynamics using Differential Network Analysis in R.

BACKGROUND: Progress in high-throughput molecular methods accompanied by more complex experimental designs demands novel data visualisation solutions. To specifically answer the question which parts of the specifical biological system are responding in particular perturbation, integrative approach in which experimental data are superimposed on a prior knowledge network is shown to be advantageous. RESULTS: We have developed DiNAR, Differential Network Analysis in R, a user-friendly application with dynamic visualisation that integrates multiple condition high-throughput data and extensive biological prior knowledge. Implemented differential network approach and embedded network analysis allow users to analyse condition-specific responses in the context of topology of interest (e.g. immune signalling network) and extract knowledge concerning patterns of signalling dynamics (i.e. rewiring in network structure between two or more biological conditions). We validated the usability of software on the Arabidopsis thaliana and Solanum tuberosum datasets, but it is set to handle any biological instances. CONCLUSIONS: DiNAR facilitates detection of network-rewiring events, gene prioritisation for future experimental design and allows capturing dynamics of complex biological system. The fully cross-platform Shiny App is hosted and freely available at https://nib-si.shinyapps.io/DiNAR. The most recent version of the source code is available at https://github.com/NIB-SI/DiNAR/ with a DOI 10.5281/zenodo.1230523 of the archived version in Zenodo.

Novel algorithm including CA-125, HE4 and body mass index in the diagnosis of endometrial cancer.

OBJECTIVES: To evaluate the diagnostic and prognostic potential of preoperative serum CA-125 and HE4 levels in patients with endometrial cancer. METHODS: Prospective case-control study of 133 women who underwent surgical treatment at the University Medical Centre Ljubljana (64 patients with endometrial cancer, 69 control patients with prolapsed uterus or myoma). Serum CA-125 and HE4 levels were determined using electrochemiluminescent assays. RESULTS: Serum CA-125 and HE4 levels were significantly higher in patients with endometrial cancer, compared to the controls (p=2.67×10-4, 1.36×10-7, respectively). A diagnostic model that combines serum CA-125 and HE4 levels and body mass index separated patients with endometrial cancer from controls, with AUC of 0.804, sensitivity of 66.7%, and specificity of 84.6%. Serum HE4 levels showed good prognostic potential and stratified the patients according to presence/absence of deep myometrial invasion (p=0.001) or lymphovascular invasion (p=0.003), with AUCs of 0.78 and 0.81, respectively. In low-risk patients with grade 1 and 2 endometrioid cancer for whom lymphadenectomy can be avoided, HE4 allowed stratification according to deep myometrial invasion (p=3.39×10-4), with AUC of 0.84. Although median HE4 levels were higher in patients with lymphovascular invasion, this difference did not reach significance (p=0.06). CONCLUSIONS: A model based on preoperative serum CA-125 and HE4 levels and body mass index has good diagnostic accuracy for separation of patients with endometrial cancer and control patients. In patients with endometrial cancer, serum HE4 levels allow prediction of deep myometrial and lymphovascular invasion.

Inter-laboratory assessment of different digital PCR platforms for quantification of human cytomegalovirus DNA.

Quantitative PCR (qPCR) is an important tool in pathogen detection. However, the use of different qPCR components, calibration materials and DNA extraction methods reduces comparability between laboratories, which can result in false diagnosis and discrepancies in patient care. The wider establishment of a metrological framework for nucleic acid tests could improve the degree of standardisation of pathogen detection and the quantification methods applied in the clinical context. To achieve this, accurate methods need to be developed and implemented as reference measurement procedures, and to facilitate characterisation of suitable certified reference materials. Digital PCR (dPCR) has already been used for pathogen quantification by analysing nucleic acids. Although dPCR has the potential to provide robust and accurate quantification of nucleic acids, further assessment of its actual performance characteristics is needed before it can be implemented in a metrological framework, and to allow adequate estimation of measurement uncertainties. Here, four laboratories demonstrated reproducibility (expanded measurement uncertainties below 15%) of dPCR for quantification of DNA from human cytomegalovirus, with no calibration to a common reference material. Using whole-virus material and extracted DNA, an intermediate precision (coefficients of variation below 25%) between three consecutive experiments was noted. Furthermore, discrepancies in estimated mean DNA copy number concentrations between laboratories were less than twofold, with DNA extraction as the main source of variability. These data demonstrate that dPCR offers a repeatable and reproducible method for quantification of viral DNA, and due to its satisfactory performance should be considered as candidate for reference methods for implementation in a metrological framework.

Applicability of human osteoarthritic chondrocytes for in vitro efficacy testing of anti-TNFα drugs.

In vitro cell-based models are important tools for assessing efficacies of new leads in early phases of drug development. Human osteoarthritic chondrocytes (OACs), obtained from biomedical waste material, represent a valuable, relatively accessible cellular source that could be used for this purpose. By employing reverse transcription-polymerase chain reaction (qRT-PCR) we compared gene expression profiles of key anabolic, catabolic and inflammatory genes of freshly isolated vs. monolayer cultured OACs (passages P0-P2) and non-stimulated vs. tumor necrosis factor alpha (TNF-α) stimulated P2 OACs. After expansion of OACs in monolayer cultures, the expression of almost all analyzed genes significantly decreased. The subsequent addition of TNF-α to OACs at P2 significantly increased expressions of all catabolic and inflammatory genes, leaving the anabolic profile almost unchanged. TNF-α-treated OACs were later utilized for efficacy testing of anti-TNF-α drugs infliximab and etanercept and both significantly reduced the expressions of all catabolic and inflammatory genes tested.

Multiplex quantification of four DNA targets in one reaction with Bio-Rad droplet digital PCR system for GMO detection.

The advantages of the digital PCR technology are already well documented until now. One way to achieve better cost efficiency of the technique is to use it in a multiplexing strategy. Droplet digital PCR platforms, which include two fluorescence filters, support at least duplex reactions and with some developments and optimization higher multiplexing is possible. The present study not only shows a development of multiplex assays in droplet digital PCR, but also presents a first thorough evaluation of several parameters in such multiplex digital PCR. Two 4-plex assays were developed for quantification of 8 different DNA targets (7 genetically modified maize events and maize endogene). Per assay, two of the targets were labelled with one fluorophore and two with another. As current analysis software does not support analysis of more than duplex, a new R- and Shiny-based web application analysis tool (http://bit.ly/ddPCRmulti) was developed that automates the analysis of 4-plex results. In conclusion, the two developed multiplex assays are suitable for quantification of GMO maize events and the same approach can be used in any other field with a need for accurate and reliable quantification of multiple DNA targets.

An international comparability study on quantification of mRNA gene expression ratios: CCQM-P103.1.

Measurement of RNA can be used to study and monitor a range of infectious and non-communicable diseases, with profiling of multiple gene expression mRNA transcripts being increasingly applied to cancer stratification and prognosis. An international comparison study (Consultative Committee for Amount of Substance (CCQM)-P103.1) was performed in order to evaluate the comparability of measurements of RNA copy number ratio for multiple gene targets between two samples. Six exogenous synthetic targets comprising of External RNA Control Consortium (ERCC) standards were measured alongside transcripts for three endogenous gene targets present in the background of human cell line RNA. The study was carried out under the auspices of the Nucleic Acids (formerly Bioanalysis) Working Group of the CCQM. It was coordinated by LGC (United Kingdom) with the support of National Institute of Standards and Technology (USA) and results were submitted from thirteen National Metrology Institutes and Designated Institutes. The majority of laboratories performed RNA measurements using RT-qPCR, with datasets also being submitted by two laboratories based on reverse transcription digital polymerase chain reaction and one laboratory using a next-generation sequencing method. In RT-qPCR analysis, the RNA copy number ratios between the two samples were quantified using either a standard curve or a relative quantification approach. In general, good agreement was observed between the reported results of ERCC RNA copy number ratio measurements. Measurements of the RNA copy number ratios for endogenous genes between the two samples were also consistent between the majority of laboratories. Some differences in the reported values and confidence intervals ('measurement uncertainties') were noted which may be attributable to choice of measurement method or quantification approach. This highlights the need for standardised practices for the calculation of fold change ratios and uncertainties in the area of gene expression profiling.

Antibacterial Activity of Wild Mushroom Extracts on Bacterial Wilt Pathogen Ralstonia solanacearum.

In total, 150 protein extracts from 94 different basidiomycete and ascomycete wild mushroom species were tested for antibacterial activity against the quarantine plant-pathogen bacterium Ralstonia solanacearum. In in vitro microtiter plate assays, 15 extracts with moderate to high antibacterial activities were identified: 11 completely inhibited bacterial growth and 4 showed partial inhibition. Of these 15 extracts, 5 were further tested and 3 extracts slowed disease progression and reduced disease severity in artificially inoculated tomato and potato plants. However, the in vitro activities of the extracts did not always correlate with their in vivo activities, which emphasizes the importance of performing early screening tests also in vivo. Testing of selected extracts against 12 R. solanacearum strains identified 6 with potential for broader applicability. Further analysis of extracts from Amanita phalloides and Clitocybe geotropa showed that the active substances are proteins with an approximate size of 180 kDa. To our knowledge, this is the first in vitro and in vivo study that demonstrates that mushroom protein extracts can be promising for treatment of bacterial wilt caused by R. solanacearum.

Transcriptome study and identification of potential marker genes related to the stable expression of recombinant proteins in CHO clones.

BACKGROUND: Chinese hamster ovary (CHO) cells have become the host of choice for the production of recombinant proteins, due to their capacity for correct protein folding, assembly, and posttranslational modifications. The most widely used system for recombinant proteins is the gene amplification procedure that uses the CHO-Dhfr expression system. However, CHO cells are known to have a very unstable karyotype. This is due to chromosome rearrangements that can arise from translocations and homologous recombination, especially when cells with the CHO-Dhfr expression system are treated with methotrexate hydrate. The present method used in the industry for testing clones for their long-term stability of recombinant protein production is empirical, and it involves their cultivation over extended periods of time prior to the selection of the most suitable clone for further bioprocess development. The aim of the present study was the identification of marker genes that can predict stable expression of recombinant genes in particular clones early in the development stage. RESULTS: The transcriptome profiles of CHO clones with stable and unstable recombinant protein production were investigated over 10-weeks of cultivation, using a DNA microarray. We identified 14 genes that were differentially expressed between the stable and unstable clones already at 2 weeks from the beginning of the cultivation. Their expression was validated by reverse-transcription quantitative real-time PCR (RT-qPCR). Furthermore, the k-nearest neighbour algorithm approach shows that the combination of the gene expression patterns of only five of these 14 genes is sufficient to predict stable recombinant protein production in clones in the early phases of cell-line development. CONCLUSIONS: The exact molecular mechanisms that cause unstable recombinant protein production are not fully understood. However, the expression profiles of some genes in clones with stable and unstable recombinant protein production allow prediction of such instability early in the cell-line development stage. We have thus developed a proof-of-concept for a novel approach to eliminate unstable clones in the CHO-Dhfr expression system, which saves time and labour-intensive work in cell-line development.

The Effect of Timing of Female Vibrational Reply on Male Signalling and Searching Behaviour in the Leafhopper Aphrodes makarovi.

Sexual communication in animals often involves duetting characterized by a coordinated reciprocal exchange of acoustic signals. We used playback experiments to study the role of timing of a female reply in the species-specific duet structure in the leafhopper Aphrodes makarovi (Hemiptera: Cicadellidae). In leafhoppers, mate recognition and location is mediated exclusively by species- and sex-specific substrate-borne vibrational signals and a female signal emitted in reply to male advertisement calls is essential for recognition and successful location of the female. In A. makarovi, males have to initiate each exchange of vibrational signals between partners, and in a duet the beginning of a female reply overlaps the end of the male advertisement call. Results of playback treatments in which female replies were delayed and did not overlap with the male call revealed that in order to trigger an appropriate behavioural response of the male, female reply has to appear in a period less than 400 ms after the end of the initiating male call. Results also suggest that males are not able to detect a female reply while calling, since female reply that did not continue after the end of male call triggered male behaviour similar to behaviour observed in the absence of female reply. Together, our results show that vibrational duets are tightly coordinated and that the species-specific duet structure plays an important role in mate recognition in location processes.

Bimodal dynamics of primary metabolism-related responses in tolerant potato-Potato virus Y interaction.

BACKGROUND: Potato virus Y (PVY) is a major pathogen that causes substantial economic losses in worldwide potato production. Different potato cultivars differ in resistance to PVY, from severe susceptibility, through tolerance, to complete resistance. The aim of this study was to better define the mechanisms underlying tolerant responses of potato to infection by the particularly aggressive PVY(NTN) strain. We focused on the dynamics of the primary metabolism-related processes during PVY(NTN) infection. RESULTS: A comprehensive analysis of the dynamic changes in primary metabolism was performed, which included whole transcriptome analysis, nontargeted proteomics, and photosynthetic activity measurements in potato cv. Désirée and its transgenic counterpart depleted for accumulation of salicylic acid (NahG-Désirée). Faster multiplication of virus occurred in the NahG-Désirée, with these plants developing strong disease symptoms. We show that while the dynamics of responses at the transcriptional level are extensive and bimodal, this is only partially translated to the protein level, and to the final functional outcome. Photosynthesis-related genes are transiently induced before viral multiplication is detected and it is down-regulated later on. This is reflected as a deficiency of the photosynthetic apparatus at the onset of viral multiplication only. Interestingly, specific and constant up-regulation of some RuBisCO transcripts was detected in Désirée plants, which might be important, as these proteins have been shown to interact with viral proteins. In SA-deficient and more sensitive NahG-Désirée plants, consistent down-regulation of photosynthesis-related genes was detected. A constant reduction in the photochemical efficiency from the onset of viral multiplication was identified; in nontransgenic plants this decrease was only transient. The transient reduction in net photosynthetic rate occurred in both genotypes with the same timing, and coincided with changes in stomatal conductivity. CONCLUSIONS: Down-regulation of photosynthesis-related gene expression and decreased photosynthetic activity is in line with other studies that have reported the effects of biotic stress on photosynthesis. Here, we additionally detected induction of light-reaction components in the early stages of PVY(NTN) infection of tolerant interaction. As some of these components have already been shown to interact with viral proteins, their overproduction might contribute to the absence of symptoms in cv. Désirée.

Transcriptomic variation between different Chinese hamster ovary cell lines.

OBJECTIVES: To identify transcription markers that uniquely determine specific Chinese hamster ovary (CHO) cell lines and can be used for the identification of cell lines in the process of biopharmaceutical cell-line development. RESULTS: Five CHO cell lines with different origins were extensively characterised at the transcriptomic level and the results were compared to their karyotype characterisation. The analysed cell lines differ in their karyotype but, due to the genome instability observed during parental and recombinant cell-line establishment, karyotyping is not the preferred method for accurate identification of the various CHO cell lines. Marker genes unique to a specific cell line were identified by microarrays, and their expression was validated by reverse-transcription quantitative real-time PCR. The analysed cell lines can be differentiated by the presence/absence of detectable marker gene expression. Additionally, the similarity of the transcriptional profiles is dependent on cell-line history but independent of the manipulation steps involved in the recombinant cell-line development process. CONCLUSIONS: Certain transcripts can be used as markers for the identification of a CHO cell line undergoing recombinant development and thus represent a powerful tool for ensuring the maintenance of high quality standards.

Analysis of Glioblastoma Patients' Plasma Revealed the Presence of MicroRNAs with a Prognostic Impact on Survival and Those of Viral Origin.

BACKGROUND: Glioblastoma multiforme (GBM) is among the most aggressive cancers with a poor prognosis in spite of a plethora of established diagnostic and prognostic biomarkers and treatment modalities. Therefore, the current goal is the detection of novel biomarkers, possibly detectable in the blood of GBM patients that may enable an early diagnosis and are potential therapeutic targets, leading to more efficient interventions. EXPERIMENTAL PROCEDURES: MicroRNA profiling of 734 human and human-associated viral miRNAs was performed on blood plasma samples from 16 healthy individuals and 16 patients with GBM, using the nCounter miRNA Expression Assay Kits. RESULTS: We identified 19 miRNAs with significantly different plasma levels in GBM patients, compared to the healthy individuals group with the difference limited by a factor of 2. Additionally, 11 viral miRNAs were found differentially expressed in plasma of GBM patients and 24 miRNA levels significantly correlated with the patients' survival. Moreover, the overlap between the group of candidate miRNAs for diagnostic biomarkers and the group of miRNAs associated with survival, consisted of ten miRNAs, showing both diagnostic and prognostic potential. Among them, hsa miR 592 and hsa miR 514a 3p have not been previously described in GBM and represent novel candidates for selective biomarkers. The possible signalling, induced by the revealed miRNAs is discussed, including those of viral origin, and in particular those related to the impaired immune response in the progression of GBM. CONCLUSION: The GBM burden is reflected in the alteration of the plasma miRNAs pattern, including viral miRNAs, representing the potential for future clinical application. Therefore proposed biomarker candidate miRNAs should be validated in a larger study of an independent cohort of patients.

Gene expression of cultured human chondrocytes as a model for assessing neutralization efficacy of soluble TNFα by TNFα antagonists.

Tumor necrosis factor-alpha (TNFα) antagonists are efficacious in the treatment of various immune-mediated inflammatory diseases. Because of rapidly growing demand for developing new or biosimilar versions of these biologicals, the need to create in vitro testing models that best represent physiological conditions is increasing. Primary human chondrocytes were used for potency evaluation and comparison between the molecular effects of anti-TNFα biologicals. Infliximab and etanercept were chosen to assess the suitability of chondrocyte cell culture for determination of anti-TNFα neutralization efficacy employing quantitative reverse transcription-polymerase chain reaction (qRT-PCR) technology. Use of both anti-TNFα biologics resulted in decrease of TNFα-stimulated expression of various matrix metalloproteinases, interleukins and other inflammation-related genes in our cell model. Significant differences in inhibition efficacy of etanercept and infliximab were observed, which were confirmed also on protein level. To evaluate the potency of anti-TNFα biologicals, a selection of TNFα-responsive target genes was made from the gene array data. The selected genes were employed in development of statistical model, which enables comparability of anti-TNFα biologicals. The presented analytical approach is suitable for assessment of the neutralization efficacy of various anti-TNFα biologicals. As such, it can be used for additional comprehensive characterization and comparability of TNF antagonists in preclinical drug testing.

Identification of plasma biomarker candidates in glioblastoma using an antibody-array-based proteomic approach.

BACKGROUND: Glioblastoma multiforme (GBM) is a brain tumour with a very high patient mortality rate, with a median survival of 47 weeks. This might be improved by the identification of novel diagnostic, prognostic and predictive therapy-response biomarkers, preferentially through the monitoring of the patient blood. The aim of this study was to define the impact of GBM in terms of alterations of the plasma protein levels in these patients. MATERIALS AND METHODS: We used a commercially available antibody array that includes 656 antibodies to analyse blood plasma samples from 17 healthy volunteers in comparison with 17 blood plasma samples from patients with GBM. RESULTS: We identified 11 plasma proteins that are statistically most strongly associated with the presence of GBM. These proteins belong to three functional signalling pathways: T-cell signalling and immune responses; cell adhesion and migration; and cell-cycle control and apoptosis. Thus, we can consider this identified set of proteins as potential diagnostic biomarker candidates for GBM. In addition, a set of 16 plasma proteins were significantly associated with the overall survival of these patients with GBM. Guanine nucleotide binding protein alpha (GNAO1) was associated with both GBM presence and survival of patients with GBM. CONCLUSIONS: Antibody array analysis represents a useful tool for the screening of plasma samples for potential cancer biomarker candidates in small-scale exploratory experiments; however, clinical validation of these candidates requires their further evaluation in a larger study on an independent cohort of patients.

Reverse transcriptase droplet digital PCR shows high resilience to PCR inhibitors from plant, soil and water samples.

BACKGROUND: Detection and quantification of plant pathogens in the presence of inhibitory substances can be a challenge especially with plant and environmental samples. Real-time quantitative PCR has enabled high-throughput detection and quantification of pathogens; however, its quantitative use is linked to standardized reference materials, and its sensitivity to inhibitors can lead to lower quantification accuracy. Droplet digital PCR has been proposed as a method to overcome these drawbacks. Its absolute quantification does not rely on standards and its tolerance to inhibitors has been demonstrated mostly in clinical samples. Such features would be of great use in agricultural and environmental fields, therefore our study compared the performance of droplet digital PCR method when challenged with inhibitors common to plant and environmental samples and compared it with quantitative PCR. RESULTS: Transfer of an existing Pepper mild mottle virus assay from reverse transcription real-time quantitative PCR to reverse transcription droplet digital PCR was straight forward. When challenged with complex matrices (seeds, plants, soil, wastewater) and selected purified inhibitors droplet digital PCR showed higher resilience to inhibition for the quantification of an RNA virus (Pepper mild mottle virus), compared to reverse transcription real-time quantitative PCR. CONCLUSIONS: This study confirms the improved detection and quantification of the PMMoV RT-ddPCR in the presence of inhibitors that are commonly found in samples of seeds, plant material, soil, and wastewater. Together with absolute quantification, independent of standard reference materials, this makes droplet digital PCR a valuable tool for detection and quantification of pathogens in inhibition prone samples.