A. Preparación de los sistemas de fermentación:
Para la realización de este proyecto, se tomaron 4 botellas de plástico de 600 mililitros (mL) a las cuales se les agregó 300 mL de agua destilada, 2 gramos (g) de Saccharomyces cerevisiae, y 45 g de azúcar comercial. A cada botella se le aplicó un tratamiento diferente basado en pH y temperatura, a 2 de ellas se les agregó ácido Clorhírico concentrado (HCL conc) hasta obtener un pH entre 2 y 3 y una de ellas se sometió a una temperatura ambiente y la otra a frío, a las botellas restantes se les añadió escamas de Hidróxido de Sodio (NaOH) llegando a tener un pH entre 9 y 11 y se repitió el procedimiento antes mencionado realizado con las botellas ácidas. Se replicó 6 veces cada tratamiento y se obtuvo en total 28 sistemas de fermentación, catorce en condiciones ácidas y catorce en condiciones básicas, por lo que al final se contaron con 4 tratamientos, los cuales son AA=ácido Ambiente, AF= ácido Frío, BA= Básico Ambiente, BF= Básico Frío. Por otro lado, las botellas se sellaron y se les hizo un agujero en la tapa por donde se introdujo el extremo de una manguera, el cual no tocó la disolución, el otro extremo de la manguera se dirigió hacia un recipiente con agua para fomentar la salida de dióxido de carbono (CO2). Con los sistemas de fermentación listos se midió el pH inicial de cada uno, luego 14 de ellos se colocaron en una temperatura fría en un rango de 3 a 4.4ºC en una refrigeradora en el Laboratorio de Biotecnología con Microalgas de la Escuela de Biología y los restantes se colocaron a temperatura ambiente en el laboratorio LAFIT en la Escuela de Química. En estas condiciones los 28 sistemas se dejaron fermentando por una semana.
B. Cuantificación de las variables:
Porcentaje de Etanol (v/v): Se destilaron las muestras con un sistema de destilación simple y se utilizó un refractómetro para medir el índice de refracción del etanol destilado y utilizando gráficas, se determinó el grado de alcohol de la mezcla obtenida. pH: Se utilizó una interfase y pH-metro para determinar el pH de las muestras iniciales y finales para luego ser comparadas. Glucosa (µg/mL): Se utilizó el método Fenol-Sulfúrico y un espectrofotómetro para determinar la cantidad de glucosa en las muestras.
C. Análisis Estadístico:
Se utilizó el análisis de varianza ANDEVA para probar la hipótesis de la diferencia de las medias de los 4 tratamientos, para las variables medidas de glucosa (µg/mL), porcentaje de alcohol (v/v) y pH. Se revisó la normalidad de los datos con la prueba Shapiro y la homocedasticidad con la prueba Barlett para escoger si usar un ANDEVA paramétrico o no paramétrico con la prueba de Fisher o Kruskal-Wallis respectivamente, se hicieron pruebas posteriori dependiendo de los resultados de normalidad y homocedasticidad, siendo la prueba TukeyHSD para el ANDEVA de Fisher y Wilcoxon con método de Bonferroni para el ANDEVA de Kruskal-Wallis. Además, se utilizó el método de regresión lineal y múltiple, para así poder obtener el mejor modelo de trabajo, primeramente se realizó shapiro.test para comprobar la normalidad de las variables, y de esta manera se pudo determinar el método de correlación a usar, si “Spearman” o “Pearson”, seguidamente se trabajó con cor.test para obtener el valor de correlación y se identificó cuales se relacionan mejor entre sí con respecto al tratamiento. Posterior a eso, se usó el modelo de regresión lineal entre dos variables (pH-Azúcar) para los tratamientos iniciales y con respecto a las regresiones finales se utilizó un modelo de regresión múltiple entre tres variables( pH-Azúcar-Etanol) con tratamientos finales.
Resultados
Se desarrollaron pruebas de ANDEVA al inicio para corroborar que todas las medias de los datos iniciales fueran homogéneas y posteriormente se volvió a realizar ANDEVA pero a las mediciones finales de cada tratamiento para observar los cambios. Cabe destacar que se desarrolló la prueba ANOVA en los tratamientos iniciales de etanol para corroborar que este último no estuviera presente y se obtuvo que todos los valores de etanol resultaron ser cero como se esperaba. (Figura 5.) Por otro lado, tanto en las ANDEVAS paramétricas como no paramétricas, todas las muestras resultaron balanceadas con un “n” igual a 7 en todos los tratamientos y se obtuvo resultados normales al aplicar la prueba de Shapiro a los residuos de las pruebas y la homocedasticidad se realizó con la prueba Barlett, la mayoría de estas variables son homocedásticos, la única que resulta ser heterocedástica es la variable del pH, tanto inicial como final, por lo que a este ensayo se le aplicó la ANDEVA no paramétrica, la cual consiste en aplicar la prueba Kruskal-Wallis para seguidamente aplicar la prueba posteriori que consta de Wilcoxon con método de Bonferroni. Con las variables que resultaron ser homocedásticas se usó la prueba de Fisher con su prueba posteriori que sería TukeyHSD. Al aplicar las pruebas posteriori se muestran cuáles datos son significativos o no significativos entre ellos como se observa en el siguiente cuadro.
Cuadro 1. Comparación entre tratamientos según la cantidad inicial de glucosa. (AA=Ácido Ambiente, AF= Ácido Frío, BA= Básico Ambiente, BF= Básico Frío).
Por otro lado, para obtener una comparación más visual, se realizaron gráficas tipo diagramas de cajas y bigotes, tanto en las variables iniciales como finales para asegurarse de que todas cumplían con lo establecido, para luego poder compararlas entre sí y observar el cambio que tuvo el proceso de fermentación después de una semana. Los diagramas mencionados se observan a continuación.
library(readxl)
Glucosa <- read_excel("C:/Users/Lenovo/Desktop/Glucosa.xls",
col_types = c("text", "numeric", "numeric"))
View(Glucosa)
Anova inicial
boxplot(Azucar$Glucosa_inicial~Azucar$Categorica, main = " Glucosa Inicial", xlab = "Tratamiento", ylab = "Cantidad Inicial µg/ml", col = "pink")
annovainicial<- aov(Azucar$Glucosa_inicial~Azucar$Categorica)
summary(annovainicial)
shapiro.test(annovainicial$residuals)
bartlett.test(Azucar$Glucosa_inicial~Azucar$Categorica)
anova<- aov(Azucar$Glucosa_inicial~Azucar$Categorica)
anova
summary(anova)
TukeyHSD(anova)
Anova final
boxplot(Azucar$Glucosa_final~Azucar$Categorica, main= "Glucosa Finales", xlab = "Tratamiento", ylab = "Cantidad Final µg/ml", col =" purple")
annovafinal<- aov(Azucar$Glucosa_final~Azucar$Categorica)
summary(annovafinal)
shapiro.test(annovafinal$residuals)
bartlett.test(Azucar$Glucosa_final~Azucar$Categorica)
anova<- aov(Azucar$Glucosa_final~Azucar$Categorica)
anova
summary(anova)
TukeyHSD(anova)
Figura 1. Comparación entre tratamientos según la cantidad inicial de glucosa. (AA=Ácido Ambiente, AF= Ácido Frío, BA= Básico Ambiente, BF= Básico Frío)
Figura 2. Comparación entre tratamientos según la cantidad final de glucosa. (AA=Ácido Ambiente, AF= Ácido Frío, BA= Básico Ambiente, BF= Básico Frío)
library(readxl)
Tabla_pH <- read_excel("C:/Users/Lenovo/Desktop/Tabla pH.xls",
col_types = c("text", "numeric", "numeric"))
View(Tabla_pH)
pH inicial
boxplot(Tabla_pH$pH_inicial~Tabla_pH$Categorica, main = "pH Inicial", xlab = "Tratamiento", ylab = "Cantidad Inicial PH", col = "red"))
annovainicialpH<- aov(Tabla_pH$pH_inicial~Tabla_pH$Categorica)
summary(annovainicialpH)
shapiro.test(annovainicialpH$residuals)
bartlett.test(Tabla_pH$pH_inicial~Tabla_pH$Categorica)
anova<-kruskal.test(Tabla_pH$pH_inicial~Tabla_pH$Categorica)
anova
pairwise.wilcox.test(Tabla_pH$pH_inicial,Tabla_pH$Categorica, p.adjust.method = "b", exact= F)
pH final
boxplot(Tabla_pH$pH_final~Tabla_pH$Categorica, main = "pH Final", xlab = " Tratamiento", ylab = "Cantidad Final pH", col = "blue")
annovafinalpH<- aov(Tabla_pH$pH_final~Tabla_pH$Categorica)
summary(annovafinalpH)
shapiro.test(annovafinalpH$residuals)
bartlett.test(Tabla_pH$pH_final~Tabla_pH$Categorica)
anova<-kruskal.test(Tabla_pH$pH_final~Tabla_pH$Categorica)
anova
pairwise.wilcox.test(Tabla_pH$pH_final,Tabla_pH$Categorica, p.adjust.method = "b", exact= )
Figura 3. Comparación entre tratamientos según la cantidad inicial de pH. (AA=Ácido Ambiente, AF= Ácido Frío, BA= Básico Ambiente, BF= Básico Frío)
Figura 4. Comparación entre tratamientos según la cantidad final de pH. (AA=Ácido Ambiente, AF= Ácido Frío, BA= Básico Ambiente, BF= Básico Frío)
Etanol
library(readxl)
Etanol <- read_excel("C:/Users/Lenovo/Desktop/Etanol.xls",
col_types = c("text", "numeric"))
View(Etanol)
Boxplot inicial
boxplot(etanol$etanol_final~etanol$Categorica, main ="Porcentaje inicial de etanol (v/v)",xlab = "Tratamiento", ylab = "Cantidad Final % (v/v)", col = "yellow")
Boxplot final
boxplot(Etanol$Etanol~Etanol$Categorica, main ="Porcentaje inicial de etanol (v/v)", xlab = "Tratamiento", ylab = "Cantidad Inicial % (v/v)", col = "green"))
annovaetanol<- aov(Etanol$Etanol~Etanol$Categorica)
summary(annovaetanol)
shapiro.test(annovaetanol$residuals) #p-value = 0.5499
bartlett.test(Etanol$Etanol~Etanol$Categorica) #p-value = 0.2585
anova<- aov(Etanol$Etanol~Etanol$Categorica)
anova
summary(anova)
TukeyHSD(anova)
Figura 5. Comparación entre tratamientos según la cantidad inicial de etanol. (AA=Ácido Ambiente, AF= Ácido Frío, BA= Básico Ambiente, BF= Básico Frío)
Figura 6. Comparación entre tratamientos según la cantidad final de etanol. (AA=Ácido Ambiente, AF= Ácido Frío, BA= Básico Ambiente, BF= Básico Frío)
Regresiones Lineales Se trabajó con el método de regresiones lineales después de haber aplicado las correlaciones, en este se pudo obtener los grados de libertad, valores de p-value del pH y de la glucosa, así mismo el p-value total, todo esto dependiendo del tratamiento aplicado inicialmente, ya que con los finales se añade la variable de etanol y se realizan regresiones múltiples. Se pudo rescatar que el mejor modelo fue el tratamiento inicial ácido frío en el que las dos variables fueron significativas y con un p-value total de 0.004615.
En este primer tratamiento inicial si se pudo realizar una regresión con las variables de pH y Glucosa. Se hizo realizó una regresión lineal en la cual se notó que fue positiva, pero no significativa, como queda evidenciada en la figura 7, con una cantidad de grados de libertad de 5.
IacidoF
library(readxl)
IacidoF <- read_excel("C:/Users/Lenovo/Desktop/IacidoF.xls",
col_types = c("text", "numeric", "numeric"))
View(IacidoF)
shapiro.test(IacidoF$pH)
shapiro.test(IacidoF$Glucosa)
cor.test(IacidoF$pH,IacidoF$Glucosa, method = "spearman")
reg<-lm(IacidoF$pH~IacidoF$Glucosa)
summary(reg)
library(ggplot2)
ggplot(IacidoF, aes(x=pH, y=Glucosa)) +
ggtitle ("Tratamiento Inicial Acido Ambiente") +
theme (plot.title = element_text(family="Comic Sans MS",
size=rel(1), #Tamaño relativo de la letra del título
vjust=2, #Justificación vertical, para separarlo del gráfico
face="bold", #Letra negrilla. Otras posibilidades "plain", "italic", "bold" y "bold.italic"
color="black", #Color del texto
lineheight=1.5)) + #Separación entre líneas
geom_point(shape=1) +
geom_smooth(method=lm)
Figura 7. Regresión lineal del tratamiento inicial ácido frío, del % pH con respecto a la cantidad de glucosa.
+ Inicial básico frío
Lo sucedido en este caso fue que los valores de todas las variables fue mayor a 0.05 por ende no se pudo aplicar regresiones del tratamiento final básico frío.Fue hecha una regresión lineal la cual dio negativa y no significativa y con una cantidad de grados de libertad de 5 como se evidencia en la figura 8.
IbasicoF
library(readxl)
IbasicoF <- read_excel("C:/Users/Lenovo/Desktop/IbasicoF.xls",
col_types = c("text", "numeric", "numeric"))
View(IbasicoF)
shapiro.test(IbasicoF$pH)
shapiro.test(IbasicoF$Glucosa)
cor.test(IbasicoF$pH,IbasicoF$Glucosa, method = "pearson")
reg<-lm(IbasicoF$pH~IbasicoF$Glucosa)
summary(reg)
library(ggplot2)
ggplot(IbasicoF, aes(x=pH, y=Glucosa)) +
ggtitle ("Tratamiento Inicial Basico Frio") +
theme (plot.title = element_text(family="Comic Sans MS",
size=rel(1), #Tamaño relativo de la letra del título
vjust=2, #Justificación vertical, para separarlo del gráfico
face="bold", #Letra negrilla. Otras posibilidades "plain", "italic", "bold" y "bold.italic"
color="black", #Color del texto
lineheight=1.5)) + #Separación entre líneas
geom_point(shape=1) +
geom_smooth(method=lm)
Figura 8. Regresión lineal del tratamiento inicial básico frío, del % pH con respecto a la cantidad de azúcar glucosa.
+ Inicial ácido ambiente
En esta regresión se pudo notar que las variables de glucosa y de pH no es un valor menor a 0.05.Se realizo una regrecion lineal en la que los resultados dieron no significativos con una pendiente positiva, de igual forma con 5 grados de libertad como se ve en la figura 9.
IacidoA
library(readxl)
IacidoA <- read_excel("C:/Users/Lenovo/Desktop/IacidoA.xls",
col_types = c("text", "numeric", "numeric"))
View(IacidoA)
shapiro.test(IacidoA$pH)
shapiro.test(IacidoA$Glucosa)
cor.test(IacidoA$pH,IacidoA$Glucosa, method = "spearman")
reg<-lm(IacidoA$pH~IacidoA$Glucosa)
summary(reg)
library(ggplot2)
ggplot(IacidoA, aes(x=pH, y=Glucosa)) +
ggtitle ("Tratamiento Inicial Acido Ambiente") +
theme (plot.title = element_text(family="Comic Sans MS",
size=rel(1), #Tamaño relativo de la letra del título
vjust=2, #Justificación vertical, para separarlo del gráfico
face="bold", #Letra negrilla. Otras posibilidades "plain", "italic", "bold" y "bold.italic"
color="black", #Color del texto
lineheight=1.5)) + #Separación entre líneas
geom_point(shape=1) +
geom_smooth(method=lm)
Figura 9. Regresión lineal del tratamiento inicial ácido ambiente, de % pH con respecto a la cantidad de azúcar Glucosa.
En este último tratamiento se pudo resaltar que las dos variables obtuvieron resultados mayores a 0.05. En este caso también realizamos una regresión lineal para evidenciar este resultado como se ve en la figura 10 donde la pendiente fue positiva y tuvo datos no significativos, con 5 grados de libertad.
lbasicoA
library(readxl)
IbasicoA <- read_excel("C:/Users/Lenovo/Desktop/IbasicoA.xls",
col_types = c("text", "numeric", "numeric"))
View(IbasicoA)
shapiro.test(IbasicoA$pH)
cor.test(IbasicoA$pH,IbasicoA$Glucosa, method = "spearman")
reg<-lm(IbasicoA$pH~IbasicoA$Glucosa)
summary(reg)
library(ggplot2)
ggplot(IbasicoA, aes(x=pH, y=Glucosa)) +
ggtitle ("Tratamiento Inicial Basico Ambiente") +
theme (plot.title = element_text(family="Comic Sans MS",
size=rel(1), #Tamaño relativo de la letra del título
vjust=2, #Justificación vertical, para separarlo del gráfico
face="bold", #Letra negrilla. Otras posibilidades "plain", "italic", "bold" y "bold.italic"
color="black", #Color del texto
lineheight=1.5)) + #Separación entre líneas
geom_point(shape=1) +
geom_smooth(method=lm)
Figura 10. Regresión lineal del tratamiento inicial básico ambiente, de % pH con respecto a la cantidad de Glucosa.
Regresiones Múltiples
+ Final ácido frío
En este primer tratamiento final si se pudo realizar una regresión en la que pudimos demostrar que hubo un regresión positiva y significativa , con una cantidad de grados de libertada de 5, con las variables de pH y Etanol, no obstante fue la única regresión que se pudo hacer ya que el valor de la variable de Glucosa dio mayor a 0.05, esto claro hablando del tratamiento final ácido frío.
FacidoF
library(readxl)
IacidoF <- read_excel("C:/Users/Lenovo/Desktop/IacidoF.xls",
col_types = c("text", "numeric", "numeric"))
View(IacidoF)
attach(FacidoF)
reg2 <- lm(pH~Azucar+Etanol)
summary(reg2)
library(ggplot2)
ggplot(FacidoF, aes(x=pH, y=Azucar)) +
ggtitle ("Tratamiento Inicial Basico Ambiente") +
theme (plot.title = element_text(family="Comic Sans MS",
size=rel(1), #Tamaño relativo de la letra del titulo
vjust=2,
face="bold", #Letra negrilla. Otras posibilidades "plain", "italic", "bold" y "bold.italic"
color="black", #Color del texto
lineheight=1.5)) +
geom_point(shape=1) +
geom_smooth(method=lm)
Figura 11. Regresión múltiple del tratamiento final ácido frío, del % pH con respecto al % de etanol. + Final básico frío
Lo acontecido en este caso fue que los valores de todas las variables fue mayor a 0.05 por ende no se pudo aplicar regresiones del tratamiento final básico frío. No obstante si se realizó la regresión múltiple pero como no hubo resultados significativos no la pusimos.
FbasicoF
library(readxl)
IbasicoF <- read_excel("C:/Users/Lenovo/Desktop/IbasicoF.xls",
col_types = c("text", "numeric", "numeric"))
attach(FbsicoF)
reg3 <- lm(pH~Azucar+Etanol)
summary(reg3)
+ Final ácido ambiente
En esta situación podemos notar que las variables de glucosa y etanol dieron un valor mayor a 0.05 y aunque el pH si diera un valor menor no se pudieron realizar regresiones con el tratamiento final acido ambiente. Se realizó la regresión múltiple pero como no hubo resultados significativos no la agregamos el trabajo.
FacidoA
library(readxl)
IacidoA <- read_excel("C:/Users/Lenovo/Desktop/IacidoA.xls",
col_types = c("text", "numeric", "numeric"))
attach(FacidoA)
reg1 <- lm(pH~Azucar+Etanol)
summary(reg1)
+ Final basico ambiente
En este último tratamiento se pudo resaltar que todos los valores de las variables fueron mayor a 0.05 por lo que no se pudo realizar ninguna regresión en el tratamiento final basico ambiente. Al realizar la regresión múltiple no obtuvimos resultados significativos por ende no la agregamos el trabajo.
FbasicoA
library(readxl)
IbasicoA <- read_excel("C:/Users/Lenovo/Desktop/IbasicoA.xls",
col_types = c("text", "numeric", "numeric"))
attach(FbsicoA)
reg4 <- lm(pH~Azucar+Etanol)
summary(reg4)
### Discusión y conclusiones
Arango, R., & Peña, C. (2008). Evaluación de la producción de etanol utilizando cepas recombinantes de Saccharomyces cerevisiae a partir de melaza de caña de azúcar. Medellín, Colombia: Corporación para Investigaciones Biológicas-CIB. Campbell, N. & Reece, J. (2007). “Biología” Séptima Edición. Madrid, España: Editorial Médica Panamericana.
Consejo Superior de Investigaciones Científicas (CSIC). (s.f).La levadura de la cerveza y… del laboratorio. Recuperado de http://seresmodelicos.csic.es/llevat.html [ Consultado el 14 de octubre del 2019].
Fernández, G. (s. f). Acidez y basicidad de alcoholes. Recuperado de https://www.quimicaorganica.org/alcoholes/409-acidez-y-basicidad-de-alcoholes.html [Consultado el 14 de octubre del 2019]. Ospina, A & Palacios, M. (1994). Efectos del cultivo de levaduras sobre la carga orgánica de los efluentes de SUCROMILES S.A. Tesis pregrado Microbiología. Facultad de Ingeniería. Universidad del Valle. Cali, Colombia. 23-29p.
Rubín, A. (2015). Saccharomyces cerevisiae: características, morfología, ciclo vital. Recuperado de https://www.lifeder.com/saccharomyces-cerevisiae/ [Consultado el 14 de octubre del 2019].
Suárez, C., Garrido, N., Guevara, C. (2016). Levadura Saccharomyces cerevisiae y la producción de alcohol. La Habana, Cuba: Instituto Cubano de Investigaciones de los Derivados de la Caña de Azúcar, 50, 20-28.(1)
Tomasso, M. (2004). Tolerancia de las levaduras al etanol. Montevideo, Uruguay: Universidad de la República Uruguay. Facultad de Química.