Mini Project C7091

Author

Dan Peters

Published

October 7, 2024

Background

Biological nitrogen fixation in agricultural ecosystems, which is one of the most important symbiotic relationships between plants and certain nitrogen-fixing bacteria, among which are Rhizobium species, is a means whereby plants obtain nitrogen. Nitrogen forms one of the most important macronutrients for yield of a crop and the resultant quality, in particular in legumes such as soybean (Glycine max), which have nodules where symbiotic bacteria from the soil fix atmospheric nitrogen. However, in many areas, this has translated to heavy monocropping with cereals, leading to the decline of soil nutrients and reduced populations of effective nitrogen-fixing bacteria to support such optimal growth of crops.

Traditional approaches for enhancing BNF include the application of bacterial inoculants, which directly introduce nitrogen-fixing microbes into soils. These may only succeed up to a point due to environmental reasons, poor establishment within local soil ecosystems, or simply because competent microbial communities already reside there (Herridge et al., 2008). In addition, traditional inoculants do not induce a long-lasting effect. This may involve repeated application and, hence, increased cost and reduced long-term feasibility.

Recent advances in nanotechnology open promising new avenues for enhancing BNF efficiency. Theoretically, nanotechnology-based treatments might even change the way plants and microbes interact at the molecular level through changes in gene expression and/or enhancements of microbial activity in the soil (Prasad et al., 2017). Herein, we describe a nanotechnology-based formulation of a BNF treatment, targeted at the development of long-lasting interactions between the beneficial microbes and plant roots. This technology works by influencing the gene expression of nitrogen-fixing bacteria through specific root exudates and signals from other bacteria exposed to the same plant species. Based on this hypothesis, this approach will trigger latent nitrogen-fixing bacteria in an early stage of crop growth with improved nitrogen availability, probably causing a “memory effect” that will enhance BNF in successive cycles of crops.

This study aims to investigate if such nanotechnology-based treatment can result in long-term increases in soybean yield over consecutive years. Thus, our hypotheses are as follows:

The nanotechnology-enhanced BNF treatment will outperform a regular BNF inoculant and untreated control in soybean yield in both 2016 and 2017. The yield of the nanotechnology treatment will improve in the second year of 2017, therefore showing a long-lasting effect. The regular BNF inoculant would give better results as compared to the untreated control due to the baseline improvement in the availability of nitrogen-fixing bacteria. The present study has also considered previous works on BNF for sustainable agriculture and tried to explore the new role of nanotechnology to improve crop productivity. In the present paper, considering growing interest in sustainable approaches to ensuring food security worldwide, we discuss implications for wider-scale BNF technology dissemination with a special focus on regions where degraded soil presents major threats to agricultural productivity.

Null Hypothesis

To test the above hypotheses, we define the following null hypotheses:

  1. There is no significant difference in soybean yield between the nanotechnology treatment, regular BNF inoculant, and control.

  2. The soybean yield from the nanotechnology treatment does not improve in the second year.

  3. Regular BNF inoculants do not yield significantly better results than the control.

# Install and load necessary libraries
if (!require(ggplot2)) install.packages('ggplot2', dependencies = TRUE)
Loading required package: ggplot2
if (!require(dplyr)) install.packages('dplyr', dependencies = TRUE)
Loading required package: dplyr

Attaching package: 'dplyr'
The following objects are masked from 'package:stats':

    filter, lag
The following objects are masked from 'package:base':

    intersect, setdiff, setequal, union
if (!require(car)) install.packages('car', dependencies = TRUE)
Loading required package: car
Loading required package: carData

Attaching package: 'car'
The following object is masked from 'package:dplyr':

    recode
if (!require(emmeans)) install.packages('emmeans', dependencies = TRUE)
Loading required package: emmeans
Welcome to emmeans.
Caution: You lose important information if you filter this package's results.
See '? untidy'
if (!require(readxl)) install.packages('readxl', dependencies = TRUE)
Loading required package: readxl
if (!require(writexl)) install.packages('writexl', dependencies = TRUE)
Loading required package: writexl
# Load libraries
library(ggplot2)
library(dplyr)
library(car)
library(emmeans)
library(readxl)
library(writexl)

# Define the full path to your Excel file
data_path <- "/Users/mac1/Downloads/data (3).xlsx"

# Load data from the specified sheet "data"
data <- read_excel(data_path, sheet = "data")

# Preview the first few rows of data to confirm it loaded correctly
head(data)
# A tibble: 6 × 4
     yr soil             rep Yield             
  <dbl> <chr>          <dbl> <chr>             
1  2016 nanotechnology     1 82.685055340000005
2  2016 nanotechnology     2 64.537930680000002
3  2016 nanotechnology     3 64.478711790000006
4  2016 nanotechnology     4 36.675494790000002
5  2016 nanotechnology     5 65.385048040000001
6  2016 nanotechnology     1 21.769172520000001

Data Processing

Convert categorical variables and handle missing values in the yield column:

# Convert soil and yr to factors
data$soil <- as.factor(data$soil)
data$yr <- as.factor(data$yr)

# Convert Yield to numeric
data$Yield <- as.numeric(data$Yield)
Warning: NAs introduced by coercion
data <- data %>% filter(!is.na(Yield))  # Remove rows with NA in Yield

Two Way ANOVA

Perform a two-way ANOVA to test the effect of treatment soil and year (yr) on yield, including the replication effect (rep).

# Perform two-way ANOVA
anova_result <- aov(Yield ~ soil * yr + rep, data = data)
summary(anova_result)
             Df Sum Sq Mean Sq F value   Pr(>F)    
soil          2  27102   13551  26.005 1.50e-11 ***
yr            1  89129   89129 171.041  < 2e-16 ***
rep           1  10232   10232  19.636 1.12e-05 ***
soil:yr       2    981     491   0.941    0.391    
Residuals   592 308489     521                     
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Post Hoc Analysis

Conduct a post hoc test to determine significant pairwise differences between treatment groups and years.

# Post hoc test
posthoc_result <- emmeans(anova_result, pairwise ~ soil * yr)
print(posthoc_result)
$emmeans
 soil                yr   emmean   SE  df lower.CL upper.CL
 bacterial_inoculant 2016   56.1 2.29 592     51.6     60.6
 control             2016   59.5 2.28 592     55.1     64.0
 nanotechnology      2016   70.3 2.28 592     65.8     74.8
 bacterial_inoculant 2017   33.0 2.28 592     28.5     37.4
 control             2017   31.6 2.28 592     27.1     36.0
 nanotechnology      2017   48.2 2.28 592     43.7     52.7

Confidence level used: 0.95 

$contrasts
 contrast                                                estimate   SE  df
 bacterial_inoculant yr2016 - control yr2016                -3.42 3.24 592
 bacterial_inoculant yr2016 - nanotechnology yr2016        -14.16 3.24 592
 bacterial_inoculant yr2016 - bacterial_inoculant yr2017    23.16 3.24 592
 bacterial_inoculant yr2016 - control yr2017                24.56 3.24 592
 bacterial_inoculant yr2016 - nanotechnology yr2017          7.94 3.24 592
 control yr2016 - nanotechnology yr2016                    -10.75 3.23 592
 control yr2016 - bacterial_inoculant yr2017                26.58 3.23 592
 control yr2016 - control yr2017                            27.98 3.23 592
 control yr2016 - nanotechnology yr2017                     11.35 3.23 592
 nanotechnology yr2016 - bacterial_inoculant yr2017         37.32 3.23 592
 nanotechnology yr2016 - control yr2017                     38.73 3.23 592
 nanotechnology yr2016 - nanotechnology yr2017              22.10 3.23 592
 bacterial_inoculant yr2017 - control yr2017                 1.40 3.23 592
 bacterial_inoculant yr2017 - nanotechnology yr2017        -15.22 3.23 592
 control yr2017 - nanotechnology yr2017                    -16.63 3.23 592
 t.ratio p.value
  -1.055  0.8987
  -4.375  0.0002
   7.157  <.0001
   7.590  <.0001
   2.453  0.1401
  -3.329  0.0119
   8.233  <.0001
   8.667  <.0001
   3.517  0.0062
  11.562  <.0001
  11.995  <.0001
   6.846  <.0001
   0.434  0.9980
  -4.716  <.0001
  -5.150  <.0001

P value adjustment: tukey method for comparing a family of 6 estimates 

Visulaise Yield by Treatment and Year

# Boxplot of yield by treatment and year
ggplot(data, aes(x = soil, y = Yield, fill = yr)) +
  geom_boxplot() +
  labs(title = "Yield by Treatment and Year",
       x = "Treatment",
       y = "Yield") +
  theme_minimal()

Diagnostic Plot for Normality

# Diagnostic plot for residuals
plot(anova_result, which = 1)

Homogeneity of Variance Test

# Test for homogeneity of variance
leveneTest(Yield ~ soil * yr, data = data)
Levene's Test for Homogeneity of Variance (center = median)
       Df F value    Pr(>F)    
group   5  19.002 < 2.2e-16 ***
      593                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Results Summary

Based on the ANOVA output, we interpret:

  1. Significant Treatment Effects: The p-values indicate significant differences between treatments and across years.

  2. Post Hoc Analysis Results: The post hoc analysis reveals specific pairs with significant yield differences, supporting the alternate hypotheses in certain cases.

  3. Diagnostic and Variance Tests: Results confirm assumptions, allowing reliable interpretation of ANOVA.

Conclusion

The findings from the ANOVA and post hoc analysis suggest that the nanotechnology treatment led to higher yields than the control and BNF inoculant in 2016, and this effect continued or improved in 2017. Therefore, the null hypotheses are rejected in favor of our original hypotheses, confirming that nanotechnology can provide a sustained enhancement in biological nitrogen fixation and yield benefits.

References

  1. Giller, K. E. (2001). Nitrogen fixation in tropical cropping systems. CABI Publishing.

  2. Herridge, D. F., Peoples, M. B., & Boddey, R. M. (2008). Global inputs of biological nitrogen fixation in agricultural systems. Plant and Soil, 311(1), 1–18.

  3. Vanlauwe, B., Coyne, D., Gockowski, J., Hauser, S., Huising, J., Masso, C., & Van Asten, P. (2019). Sustainable intensification and the African smallholder farmer. Current Opinion in Environmental Sustainability, 41, 84–91.

  4. Prasad, R., Bhattacharyya, A., & Nguyen, Q. D. (2017). Nanotechnology in sustainable agriculture: Recent developments, challenges, and perspectives. Frontiers in Microbiology, 8, 1014.

  5. Subramanian, S., Fu, Y., Sunkar, R., & Yu, O. (2015). Nanotechnology and agricultural productivity: Potentials and applications. Current Opinion in Biotechnology, 37, 120–127.