Avaliação do efeito de produtos biológicos e químicos após a germinação de Colletotrichum spp.

Autores

Prof. Dr. Walmes Zeviani

Débora Petermann

1 Sobre o experimento

Introdução

Este estudo tem como objetivo avaliar o efeito de diferentes produtos biológicos e um produto químico após a germinação dos conídios de Colletotrichum nymphaeae e C. chrysophilum, incluindo cepas sensíveis e com sensibilidade reduzida a fungicidas. A variável de interesse é o diâmetro das colônias, medido em milímetros.

Materiais e Métodos

Os experimentos 1 e 2 foram realizados no Laboratório de Epidemiologia para o Manejo Integrado de Doenças de Plantas (LEMID) sob a supervisão da aluna de mestrado Débora Petermann, orientada pela professora doutora Louise Larissa May De Mio.

  • Experimento 2:
    • Data de montagem: 10 de junho de 2022.
    • Data de avaliação: 15 de junho de 2022.
    • Estrutura do experimento: 8 tratamentos, 4 isolados e 3 repetições (placas), com dois valores por placa devido à adição de duas gotas do produto + esporo.

Foram utilizados os seguintes tratamentos:

  1. Bacillus amyloliquefaciens na dosagem de 0,5g/L (Bam 0,5g/L) - produto biológico.
  2. Bacillus amyloliquefaciens na dosagem de 1,0g/L (Bam 1,0g/L) - produto biológico.
  3. Bacillus amyloliquefaciens na dosagem de 2,0g/L (Bam 2,0g/L) - produto biológico.
  4. Bacillus subtilis (Bs) - produto biológico.
  5. Sobrenadante de Bacillus alcalophilus (SBa) - produto biológico.
  6. Biomassa de Bacillus alcalophilus (BBa) - produto biológico.
  7. CabrioTop - produto químico.
  8. Testemunha composta por água esterilizada (Controle).

Os isolados utilizados foram:

  • MdCn-142: C. nymphaeae sensível a fungicida.
  • MdCn-181: C. nymphaeae sensibilidade reduzida a fungicida.
  • MdCc-110: C. chrysophillum sensível a fungicida.
  • MdCc-55PR: C. chrysophillum sensibilidade reduzida a fungicida.

Metodologia

Neste experimento, após a germinação dos conídios e antes de paralisar a germinação com lactofenol, foram retiradas gotas contendo o produto e os esporos dos fungos após 12 horas de contato. Essas gotas foram então colocadas em placas de vidro contendo um meio de cultura rico em nutrientes. Essa metodologia foi aplicada apenas no experimento 2. A decisão de adotar essa abordagem foi tomada devido às altas porcentagens de germinação observadas no experimento 1, o que levantou preocupações sobre a eficácia dos produtos. A literatura indicou que, mesmo com o contato com alguns produtos, os esporos podem germinar, mas a germinação pode não ser viável. Portanto, ao transferir os esporos para o meio de cultura, pôde-se observar que, embora a germinação tenha ocorrido em muitos tratamentos, o crescimento subsequente das colônias foi inibido, indicando a inviabilidade da germinação.

Análise Estatística

A análise estatística será conduzida por meio de uma análise fatorial, em que um dos fatores será composto pelos tratamentos (8 tratamentos) e o segundo fator pelos isolados (4 isolados). A ênfase será na comparação entre os dois isolados de C. nymphaeae (sensível e sensibilidade reduzida) e os dois isolados de C. chrysophilum (sensível e sensibilidade reduzida). É importante ressaltar que os experimentos 1 e 2 serão analisados separadamente.

O crescimento apresentou vários zeros indicando que não houve germinação. Para contornar o muito provável afastamento dos pressupostos decorrente do excesso de valores iguais a zero, optou-se por cortar os valores em classe e análisar como resposta de sobreviência, ainda mantendo a suposição de distribuição normal. Dessa maneira, valores iguais a zero ou muito próximo de, que caírem na primeira classe, serão considerados censuras intervalares. Valores positivos que estiverem além da primeira classe, não serão considerados censurados. Para amplitude de classe considerou-se 2.5 mm.

Ao final, pode-se apresentar as médias amostrais de diâmetro ao lado das estimativas obtidas com o modelo de sobreviência.

2 Análise de dados

Código
#-----------------------------------------------------------------------
# Pacotes.

library(car)
library(emmeans)
library(survival)
library(tidyverse)

#-----------------------------------------------------------------------
# Importação e preparo.

# Importação.
tb <- readxl::read_excel(
    "analises_artigo_mestrado.xlsx",
    sheet = "5 plaqueamento")

# Nomes em caixa baixa.
tb <- tb |>
    rename_with(tolower) |>
    rename("diameter" = "colony diameter (mm)")
str(tb)
tibble [192 × 8] (S3: tbl_df/tbl/data.frame)
 $ experiment: num [1:192] 2 2 2 2 2 2 2 2 2 2 ...
 $ treatment : chr [1:192] "Bam 0,5g/L" "Bam 0,5g/L" "Bam 0,5g/L" "Bam 0,5g/L" ...
 $ isolated  : chr [1:192] "MdCc-110" "MdCc-110" "MdCc-110" "MdCc-110" ...
 $ species   : chr [1:192] "C. chrysophillum" "C. chrysophillum" "C. chrysophillum" "C. chrysophillum" ...
 $ phenotype : chr [1:192] "Sensitive" "Sensitive" "Sensitive" "Sensitive" ...
 $ colony    : num [1:192] 1 2 1 2 1 2 1 2 1 2 ...
 $ repetition: num [1:192] 1 1 2 2 3 3 1 1 2 2 ...
 $ diameter  : num [1:192] 0 0 0 0 0 0 0 0 0 0 ...
Código
# Alguns erros de digitação da `specie` precisam ser corrigidos.
tb_fix <-
    tb |>
    count(isolated, species, phenotype, sort = TRUE) |>
    slice_head(n = 4)

# Corrigindo valores de `specie` e `phenotype` para `isolated`.
tb <-
    tb |>
    select(-species, -phenotype) |>
    left_join(tb_fix[, 1:3], by = "isolated")

# Variáveis `specie` e `phenotype` são características do isolado.
tb |>
    count(isolated, species, phenotype, sort = TRUE)
# A tibble: 4 × 4
  isolated  species          phenotype               n
  <chr>     <chr>            <chr>               <int>
1 MdCc-110  C. chrysophillum Sensitive              48
2 MdCc-55PR C. chrysophillum Reduced sensitivity    48
3 MdCn-142  C. nymphaeae     Sensitive              48
4 MdCn-181  C. nymphaeae     Reduced sensitivity    48
Código
# Variável `experiment` é a repetição do experimento.
tb |>
    xtabs(formula = ~treatment + experiment)
            experiment
treatment     2
  Bam 0,5g/L 24
  Bam 1,0g/L 24
  Bam 2,0g/L 24
  BBa        24
  Bs         24
  CabrioTop® 24
  Control    24
  SBa        24
Código
# A frequência de cada combinação experimental.
tb |>
    xtabs(formula = ~experiment + isolated + treatment) |>
    ftable()
                     treatment Bam 0,5g/L Bam 1,0g/L Bam 2,0g/L BBa Bs CabrioTop® Control SBa
experiment isolated                                                                          
2          MdCc-110                     6          6          6   6  6          6       6   6
           MdCc-55PR                    6          6          6   6  6          6       6   6
           MdCn-142                     6          6          6   6  6          6       6   6
           MdCn-181                     6          6          6   6  6          6       6   6
Código
# Passa para fator.
tb <- tb |>
    mutate_at(c("treatment", "isolated", "experiment"), factor)

# Coloca em um ordem mais lógica.
# levels(tb$treatment) |> dput()
tb$treatment <- fct_relevel(tb$treatment,
                            "Control",
                            "Bam 0,5g/L", "Bam 1,0g/L", "Bam 2,0g/L",
                            "SBa", "BBa", "Bs",
                            "CabrioTop®")

#-----------------------------------------------------------------------
# Visualização.

ggplot(data = tb,
       mapping = aes(x = treatment, y = diameter,
                     color = isolated, group = isolated)) +
    geom_point() +
    stat_summary(fun = "mean", geom = "line")

Código
# ATTENTION: Presença de zeros que corresponde a esporos que não
# germinaram.

# Qual a frequência de zeros?
tb |>
    mutate(zero = ifelse(diameter == 0, "zero", "positive")) |>
    count(treatment, experiment, isolated, zero) |>
    pivot_wider(names_from = zero,
                values_from = n,
                values_fill = 0) |>
    arrange(zero) |>
    print(n = Inf)
# A tibble: 32 × 5
   treatment  experiment isolated  positive  zero
   <fct>      <fct>      <fct>        <int> <int>
 1 Control    2          MdCc-110         6     0
 2 Control    2          MdCn-142         6     0
 3 Control    2          MdCn-181         6     0
 4 SBa        2          MdCc-110         6     0
 5 SBa        2          MdCn-142         6     0
 6 BBa        2          MdCc-110         6     0
 7 BBa        2          MdCc-55PR        6     0
 8 BBa        2          MdCn-142         6     0
 9 BBa        2          MdCn-181         6     0
10 Control    2          MdCc-55PR        5     1
11 Bam 0,5g/L 2          MdCc-110         0     6
12 Bam 0,5g/L 2          MdCc-55PR        0     6
13 Bam 0,5g/L 2          MdCn-142         0     6
14 Bam 0,5g/L 2          MdCn-181         0     6
15 Bam 1,0g/L 2          MdCc-110         0     6
16 Bam 1,0g/L 2          MdCc-55PR        0     6
17 Bam 1,0g/L 2          MdCn-142         0     6
18 Bam 1,0g/L 2          MdCn-181         0     6
19 Bam 2,0g/L 2          MdCc-110         0     6
20 Bam 2,0g/L 2          MdCc-55PR        0     6
21 Bam 2,0g/L 2          MdCn-142         0     6
22 Bam 2,0g/L 2          MdCn-181         0     6
23 SBa        2          MdCc-55PR        0     6
24 SBa        2          MdCn-181         0     6
25 Bs         2          MdCc-110         0     6
26 Bs         2          MdCc-55PR        0     6
27 Bs         2          MdCn-142         0     6
28 Bs         2          MdCn-181         0     6
29 CabrioTop® 2          MdCc-110         0     6
30 CabrioTop® 2          MdCc-55PR        0     6
31 CabrioTop® 2          MdCn-142         0     6
32 CabrioTop® 2          MdCn-181         0     6
Código
# Algumas combinações experimentais tem todas as repetições com zero de
# crescimento micelial. Isso corresponde a problemas com pressupostos.
# Seria interessante cortar em classes de crescimento micelial e
# analisar como uma resposta de sobreviência.

# Faça um histograma da variável `diameter` por `experiment`.
ggplot(data = tb,
       aes(x = diameter)) +
    geom_histogram(binwidth = 2.5, color = "black") +
    geom_rug()

Código
#-----------------------------------------------------------------------
# Corta a resposta em classe.

# Cortar a resposta em classes.
binwidth <- 2.5
breaks <- seq(from = 0,
              to = max(tb$diameter) + binwidth,
              by = binwidth)
tb$diameter_class <- cut(tb$diameter,
                         breaks = breaks,
                         include.lowest = TRUE)

# Separar valores para criar os limites do intervalo.
tb <-
    tb |>
    mutate(newvalue = as.character(diameter_class),
           newvalue = trimws(
               str_replace_all(newvalue,
                               pattern = "[\\[(\\]]",
                               replacement = " "))) |>
    separate(newvalue,
             into = c("left", "right"),
             sep = ",",
             convert = TRUE)

# IMPORTANT: Aqui foi feito censura intervalar para a primeira classe
# que é o intervalo que contém o 0. Valores positivos que cairam nas
# demais classes são declarados com eventos sem censura.

# `i` indica se é a classe que contém o 0 ou não.
tb <-
    tb |>
    mutate(i = diameter_class == levels(diameter_class)[1],
           status = ifelse(i, yes = 3, no = 1),
           right = ifelse(i, yes = right, no = diameter),
           left = ifelse(i, yes = left, no = diameter))

2.1 Experimento 2

Código
#-----------------------------------------------------------------------
# Experimento 2 --------------------------------------------------------

tbi <- tb |>
    filter(experiment == "2")

# Ajuste do modelo.
m0 <- survreg(
    formula = Surv(
        time = left + 0.01,
        time2 = right,
        event = status,
        type = "interval") ~ treatment * isolated,
    dist = "gaussian",
    # dist = "loggaussian",
    # dist = "loglogistic",
    # dist = "weibull",
    data = tbi)

# logLik(m0)
# AIC(m0)
# BIC(m0)

# Quadro de análise de deviance.
# anova(m0)
Anova(m0)
Analysis of Deviance Table (Type II tests)
                   LR Chisq Df Pr(>Chisq)    
treatment            331.02  7    < 2e-16 ***
isolated              11.07  3    0.01137 *  
treatment:isolated   260.45 21    < 2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Código
# Comparar entre isolados de mesma espécie.
tb_iso <-
    tbi |>
    distinct(species, isolated) |>
    with({
        split(as.character(isolated), species)
    })
tb_iso
$`C. chrysophillum`
[1] "MdCc-110"  "MdCc-55PR"

$`C. nymphaeae`
[1] "MdCn-142" "MdCn-181"
Código
# Para armazenar as comparações de médias.
tb_means <- list(isolated = NULL, treatment = NULL)

# Aplica o teste de comparação de médias separado por espécie.
tb_means$isolated[["C. nymphaeae"]] <-
    emmeans(m0,
            specs = ~isolated | treatment,
            at = list(isolated = tb_iso[["C. nymphaeae"]])) |>
    multcomp::cld(Letters = letters, reversed = TRUE) |>
    as.data.frame() |>
    mutate(.group = trimws(.group))
tb_means$isolated[["C. chrysophillum"]] <-
    emmeans(m0,
            specs = ~isolated | treatment,
            at = list(isolated = tb_iso[["C. chrysophillum"]])) |>
    multcomp::cld(Letters = letters, reversed = TRUE) |>
    as.data.frame() |>
    mutate(.group = trimws(.group))
tb_means$isolated <- bind_rows(tb_means$isolated, .id = "Specie")

if (m0$dist == "weibull") {
    tb_means$isolated <-
        tb_means$isolated |>
        mutate(weibull_shape = 1/m0[["scale"]]) |>
        mutate_at(c("emmean", "lower.CL", "upper.CL"),
                  ~exp(.) * gamma(1 + 1/weibull_shape)) |>
        mutate(weibull_shape = NULL)
}

ggplot(data = tb_means$isolated,
       mapping = aes(y = isolated, x = emmean)) +
    facet_grid(facets = treatment ~ .) +
    geom_errorbarh(mapping = aes(xmin = lower.CL, xmax = upper.CL),
                   height = 0) +
    geom_point(mapping = aes(color = Specie)) +
    geom_text(mapping = aes(label = sprintf("%0.1f %s", emmean, .group)),
              vjust = 0, nudge_y = 0.10, size = 2.5) +
    labs(y = "Isolated",
         x = "Colony diameter (mm)")

Código
# Comparação entre os tratamentos.
tb_means$treatment <-
    emmeans(m0,
            specs = ~treatment | isolated) |>
    multcomp::cld(Letters = letters, reversed = TRUE) |>
    as.data.frame() |>
    mutate(.group = trimws(.group)) |>
    identity()
tb_means$treatment
    treatment  isolated   emmean        SE  df   lower.CL  upper.CL .group
2     Control  MdCc-110 21.61750 0.8179044 159 20.0021420 23.232858      a
3         BBa  MdCc-110 19.93083 0.8179044 159 18.3154753 21.546191      a
4         SBa  MdCc-110 19.19583 0.8179044 159 17.5804753 20.811191      a
7  CabrioTop®  MdCc-110  1.25500 0.8729915 159 -0.4691549  2.979155      b
8          Bs  MdCc-110  1.25500 0.8729915 159 -0.4691549  2.979155      b
5  Bam 2,0g/L  MdCc-110  1.25500 0.8729915 159 -0.4691549  2.979155      b
6  Bam 1,0g/L  MdCc-110  1.25500 0.8729915 159 -0.4691549  2.979155      b
1  Bam 0,5g/L  MdCc-110  1.25500 0.8729915 159 -0.4691549  2.979155      b
13        BBa MdCc-55PR 28.33167 0.8179044 159 26.7163087 29.947025      a
10    Control MdCc-55PR 23.16571 0.8185183 159 21.5491408 24.782282      b
11 CabrioTop® MdCc-55PR  1.25500 0.8729915 159 -0.4691549  2.979155      c
12         Bs MdCc-55PR  1.25500 0.8729915 159 -0.4691549  2.979155      c
15 Bam 2,0g/L MdCc-55PR  1.25500 0.8729915 159 -0.4691549  2.979155      c
16 Bam 1,0g/L MdCc-55PR  1.25500 0.8729915 159 -0.4691549  2.979155      c
9  Bam 0,5g/L MdCc-55PR  1.25500 0.8729915 159 -0.4691549  2.979155      c
14        SBa MdCc-55PR  1.25500 0.8729915 159 -0.4691549  2.979155      c
18        BBa  MdCn-142 19.70250 0.8179044 159 18.0871420 21.317858      a
19        SBa  MdCn-142 19.28750 0.8179044 159 17.6721420 20.902858      a
20    Control  MdCn-142 18.78833 0.8179044 159 17.1729753 20.403691      a
23 CabrioTop®  MdCn-142  1.25500 0.8729915 159 -0.4691549  2.979155      b
24         Bs  MdCn-142  1.25500 0.8729915 159 -0.4691549  2.979155      b
17 Bam 2,0g/L  MdCn-142  1.25500 0.8729915 159 -0.4691549  2.979155      b
21 Bam 1,0g/L  MdCn-142  1.25500 0.8729915 159 -0.4691549  2.979155      b
22 Bam 0,5g/L  MdCn-142  1.25500 0.8729915 159 -0.4691549  2.979155      b
29    Control  MdCn-181 22.46750 0.8179044 159 20.8521420 24.082858      a
27        BBa  MdCn-181 16.75000 0.8179044 159 15.1346420 18.365358      b
31 CabrioTop®  MdCn-181  1.25500 0.8729915 159 -0.4691549  2.979155      c
26 Bam 2,0g/L  MdCn-181  1.25500 0.8729915 159 -0.4691549  2.979155      c
28 Bam 0,5g/L  MdCn-181  1.25500 0.8729915 159 -0.4691549  2.979155      c
32         Bs  MdCn-181  1.25500 0.8729915 159 -0.4691549  2.979155      c
30 Bam 1,0g/L  MdCn-181  1.25500 0.8729915 159 -0.4691549  2.979155      c
25        SBa  MdCn-181  1.25500 0.8729915 159 -0.4691549  2.979155      c
Código
if (m0$dist == "weibull") {
    tb_means$treatment <-
        tb_means$treatment |>
        mutate(weibull_shape = 1/m0[["scale"]]) |>
        mutate_at(c("emmean", "lower.CL", "upper.CL"),
                  ~exp(.) * gamma(1 + 1/weibull_shape)) |>
        mutate(weibull_shape = NULL)
}

ggplot(data = tb_means$treatment,
       mapping = aes(y = treatment, x = emmean)) +
    facet_grid(facets = isolated ~ .) +
    geom_errorbarh(mapping = aes(xmin = lower.CL, xmax = upper.CL),
                   height = 0) +
    geom_point() +
    geom_text(mapping = aes(label = sprintf("%0.1f %s", emmean, .group)),
              vjust = 0, nudge_y = 0.10, size = 2.5) +
    labs(y = "Treatment",
         x = "Colony diameter (mm)")

Código
#-----------------------------------------------------------------------
# Comparando com estimativas amostrais.

tbi |>
    group_by(treatment, isolated) |>
    summarise_at("diameter", "mean") |>
    left_join(tb_means$treatment,
              by = c("treatment", "isolated")) |>
    select(treatment, isolated, diameter, emmean) |>
    ungroup() |>
    # arrange(diameter) |>
    print(n = Inf)
# A tibble: 32 × 4
   treatment  isolated  diameter emmean
   <fct>      <fct>        <dbl>  <dbl>
 1 Control    MdCc-110      21.6  21.6 
 2 Control    MdCc-55PR     22.8  23.2 
 3 Control    MdCn-142      18.8  18.8 
 4 Control    MdCn-181      22.5  22.5 
 5 Bam 0,5g/L MdCc-110       0     1.25
 6 Bam 0,5g/L MdCc-55PR      0     1.25
 7 Bam 0,5g/L MdCn-142       0     1.25
 8 Bam 0,5g/L MdCn-181       0     1.25
 9 Bam 1,0g/L MdCc-110       0     1.25
10 Bam 1,0g/L MdCc-55PR      0     1.25
11 Bam 1,0g/L MdCn-142       0     1.25
12 Bam 1,0g/L MdCn-181       0     1.25
13 Bam 2,0g/L MdCc-110       0     1.25
14 Bam 2,0g/L MdCc-55PR      0     1.25
15 Bam 2,0g/L MdCn-142       0     1.25
16 Bam 2,0g/L MdCn-181       0     1.25
17 SBa        MdCc-110      19.2  19.2 
18 SBa        MdCc-55PR      0     1.25
19 SBa        MdCn-142      19.3  19.3 
20 SBa        MdCn-181       0     1.25
21 BBa        MdCc-110      19.9  19.9 
22 BBa        MdCc-55PR     28.3  28.3 
23 BBa        MdCn-142      19.7  19.7 
24 BBa        MdCn-181      16.7  16.8 
25 Bs         MdCc-110       0     1.25
26 Bs         MdCc-55PR      0     1.25
27 Bs         MdCn-142       0     1.25
28 Bs         MdCn-181       0     1.25
29 CabrioTop® MdCc-110       0     1.25
30 CabrioTop® MdCc-55PR      0     1.25
31 CabrioTop® MdCn-142       0     1.25
32 CabrioTop® MdCn-181       0     1.25