El incremento en las emisiones de CO2 durante el último siglo, ha producido cambios en los océanos, produciendo aumento de la temperatura y disminución del pH en ellos (Intergovernmental Panel on Climate Change, 2007). En el caso particular del mar Caribe, Chollett, Mumby, Müller-Karger, y Hu (2012), reportan que la temperatura ha aumentado 1.41 °C durante el último siglo (de 0.19 °C a 0.27 °C por decenio), alcanzando una temperatura media de 28 °C. Por otra parte, Magaña Rueda, Graizbord, Buenfil Friedman y Gómez Mendoza (2009) sugieren que la temperatura en esta región podría ascender hasta los 31 °C hacia el año 2100.
Modelaciones de acidificación señalan que el pH de los océanos podría descender de 8.1 (media actual) a 7.6 hacia el año 2100 (Caldeira & Wickett, 2003; Orr et al., 2005). De acuerdo con Gledhill, Wanninkhof, Millero y Eakin (2008), se calcula que la presión parcial de CO2 (pCO2) en el Caribe ha aumentado + 2.2 µatm.año-1 durante el último siglo, produciendo una disminución de 0.05 unidades de pH.
Cambios en estos parámetros afectan directamente la dinámica de las poblaciones, especialmente a organismos calcificadores, los cuales son más vulnerables durante su fase larvaria (Sale et al., 2010).
Roller y Strickle (1989), observaron que al aumentar la temperatura (de 20 a 30 °C) se presenta un incremento en las tasas de desarrollo, crecimiento y mortalidad para Thais haemastoma canaliculata, por otra parte, Watson, Southgate, Tyler, y Peck (2009), mencionan que larvas de Saccostrea glomerata expuestas a reducción de pH (de -0.1 a -0.2 unidades) presentan disminución en sus tasas de crecimiento y desarrollo, así como una esqueletogénesis anormal.
Strombus gigas (Linnaeus, 1758) es un organismo con distribución metapoblacional en el Caribe, presenta un estadio larval, cuya cinética de desarrollo ha sido estudiada en condiciones normales de temperatura (28 - 29 °C), por lo que no se conoce el efecto que el incremento de temperatura y la disminución del pH ejercen durante su etapa larval. En este sentido, la finalidad del presente estudio fue evaluar el efecto del incremento de la temperatura y la acidificación en la eclosión, desarrollo, crecimiento, mortalidad y calcificación de la concha durante su fase larvaria.
MATERIALES Y MÉTODOS
Cultivo larval: Se recolectó un segmento de masa ovígera de S. gigas recién desovada en Quintana Roo (20°19’00”N - 87°21’25” W, agosto 2015), la cual fue transportada para su cultivo a la estación Marina del Cinvestav-IPN (Telchac, Yucatán, 21°20’28” N - 89°18’25” W). La incubación de embriones y el cultivo larval se realizaron de agosto a octubre 2015, empleando la condición 28 °C - pH 8.1 como control y tres tratamientos de temperatura y pH (T1 = 28 °C - pH 7.6, T2 = 31 °C - pH 8.1, T3 = 31 °C - pH 7.6.), por triplicado. Los niveles de temperatura y pH utilizados se establecieron de acuerdo a las condiciones actuales y las proyecciones hacia 2100 para la región Norte del Caribe Mexicano (Caldeira & Wickett, 2003; Chollett et al., 2012; Magaña Rueda et al., 2009). La temperatura se controló utilizando calentadores para acuarios. El pH de los tratamientos se reguló mediante burbujeo, suministrando CO2 cada 4 hr. Temperatura y pH se monitorearon con un potenciómetro Oakton Mod. pH 5+.
En cada unidad de cultivo se colocó un segmento de masa ovígera de 15 cm de longitud (± 106 huevos por cm lineal). Se registró el tiempo a la eclosión (hr) para cada tratamiento. A partir de la eclosión, las larvas fueron distribuidas manualmente a unidades de cultivo de 5 L a razón de 160 larvas.L-1.
La alimentación se realizó por las mañanas (frecuencia: 24 hr), suministrando una mezcla de microalgas (Nannochloropsis oculata, Isochrysis galbana, Chaetoceros mulleri, Pavlovav sp. y Navicula sp.), a una concentración total de 1 000 cel.mL-1 por unidad de cultivo (200 cel.mL-1 de cada microalga por ración diaria).
Se realizaron recambios totales cada 5 días empleando agua de mar filtrada mediante filtros de polipropileno (5 µm) y esterilizada con UV. Se utilizaron mallas de 180 µm para la retención de las larvas.
Eclosión, desarrollo, crecimiento y mortalidad: Se registró el tiempo inicial y final de eclosión para cada tratamiento, posteriormente, se contabilizó el número de huevos no eclosionados por segmento de masa ovígera. El porcentaje de eclosión se calculó de la siguiente manera:
A fin de evaluar el desarrollo, cada 24 hr usando una pipeta Pasteur se extrajeron 30 larvas por tratamiento, las cuales se observaron (in vivo) a través un microscopio óptico LEICA ICC50. El desarrollo se observó de acuerdo a Brito-Manzano, Aldana-Aranda, de la Cruz-Lázaro, y Estrada-Botello (2006), los indicadores fueron: número de lóbulos del velum, número de espiras de la concha, corazón funcional (larval y adulto), opérculo, branquias (Fig. 1 A) y conducta nado-reptante. El crecimiento (N = 30 larvas por edad/tratamiento) y mortalidad fueron evaluados en larvas de 0, 10, 20 y 30 días. Se midió la longitud sifonal (µm) de la concha (Fig. 1 B) utilizando el software Leica LAS Interactive Measurements (resolución 10x). La mortalidad fue calculada a partir de las conchas vacías obtenidas durante los recambios de agua. El porcentaje de mortalidad se obtuvo con la siguiente fórmula:
Calcificación: Se determinó la proporción en peso de Calcio (% wt) para larvas de 30 días (N = 3 larvas por tratamiento), se eligió esta edad dado que la teleoconcha se forma entre los 30 y 35 días (Enriquez-Díaz, Volland, Chavez-Villegas, Aldana-Aranda, & Gros, 2015). Para la cuantificación se utilizó Microscopia Electrónica de Barrido empleando el Análisis de Espectroscopía de Energía Dispersiva “EDX” (Microscopio PHILIPS ESEM-XL30). La medición se hizo en cuatro zonas de la concha siguiendo la longitud sifonal (Fig. 1B). Además, se utilizó Microscopía RAMAN (Microscopio WITEC ALPHA 300, Láser 488nm), empleando las distancias raman (cm-1) propuestas por Nehrke, Poigner, Wilhelms-Dick, Brey, & Abele (2012) para la identificación y cuantificación de aragonita (206, 705 y 1 085 cm-1) y calcita (282, 711 y 1 085 cm-1) (N = 3 - 5 larvas por tratamiento). El efecto de fluorescencia registrado en las muestras se eliminó mediante el software CrystalSleuth–RRUFF Project. La intensidad RAMAN se calculó con el software OriginPro 8, y se reportó en unidades arbitrarias (a.u.).
Se utilizó el software Infostat versión estudiantil, para calcular media y desviación estándar, así como realizar prueba de normalidad (Shapiro-Wilks), de homocedasticidad (Kolmogorov) y análisis de Kruskal-Wallis (P (α0.05)) para eclosión, crecimiento, mortalidad y calcificación (EDX) entre tratamientos de temperatura y pH.
RESULTADOS
Eclosión: Se presentó a partir del día cinco de incubación, iniciando en larvas de 31 °C a partir de las 110 hr posteriores al desove, y las 118 hr para los tratamientos de 28 °C. La duración de la eclosión fue de 18 hr para 31 °C y de 30 hr para 28 °C (Fig. 2A). En general, se registró ~90 % de eclosión, con excepción del T2 (31 °C - pH 8.1) donde se obtuvo un 96.67 ± 1.39 % de larvas eclosionadas (Fig. 2B). No hubo variación significativa en la proporción de larvas eclosionadas entre tratamientos (H (g.l. 3) = 6.38, P (α0.05) = 0.0943), pero si en el tiempo de eclosión (H (g.l. 3) = 8.31, P (α0.05) < 0.0001).
Desarrollo larvario: La cinética de desarrollo del velum fue mayor para larvas de 31 °C - 8.1pH (T2), donde la reducción de lóbulos velares (indicador del inicio de la metamorfosis) comenzó a parir del día 38, mientras que en el control (28 °C - pH 8.1) se observó a partir del día 41. La concha mostró menor desarrollo en larvas del T1 (28 °C - pH 7.6), las cuales no desarrollaron 4.5 espiras, característica observada en el resto de los tratamientos. La elongación del opérculo, la aparición del corazón adulto funcional, el desarrollo de las branquias e inicio de la conducta nado-reptante, se presentó en menor tiempo y mayor proporción en los tratamientos de 31 °C (Cuadro 1).
Crecimiento: Las larvas presentaron tallas de 230 ± 4.12 µm (28 °C - pH 7.6) a 313.27 ± 11.34 µm (31 °C - pH 8.1) durante la eclosión. Se registró un crecimiento similar entre el control y los tratamientos hacia el día 10 (L.S. = ~650 µm), hacia el día 20, larvas de los tratamientos 31 °C fueron más grandes (28 °C = ~710 µm y 31 °C = ~820 µm), mismo patrón se registró durante todo el desarrollo, alcanzando las mayores tallas en los tratamientos de 31 °C (pH 7.6 = 1 031.50 ± 11.03 µm; pH 8.1 = 1 054.50 ± 11.13 µm) (Fig. 3). El análisis de Kruskal-Wallis mostró variación significativa en el crecimiento entre tratamientos (H (g.l. 3) = 26.19, P (α0.05) < 0.0001).
Mortalidad: Se registró de 0.59 ± 0.19 % a 3.01 ± 0.25 % de mortalidad durante el día 10 (Fig. 4A), incrementando de 1.38 ± 0.54 % a 5.67 ± 1.16 % hacia el día 20 (Fig. 4B). Se presentó un pico de mortalidad el día 30 para todos los tratamientos (de 19.13 % a 67.38), mostrando mayores registros en el T3 con 58.42 ± 12.03 % (Fig. 4C). El mayor porcentaje de mortalidad acumulada fue en los tratamientos de 31 °C (pH 8.1 = 65.64 ± 8.49 % y pH 7.6 = 65.10 ± 15.49 %). La mortalidad no mostró variación significativa entre tratamientos (H (g.l. 3) = 5.41, P (α0.05) = 0.1441), mientras que en el análisis por edad, se observó variación solo para las observaciones realizadas durante el día 20 (H (g.l. 3) = 8.74, P (α0.05) = 0.0324).
Calcificación:
EDX: La proporción de Calcio (% wt) fue similar entre tratamientos ~35 % wt, siendo mayor en el control con 37.29 ± 16.81 % wt (Fig. 5). No se presentó variación significativa para la proporción de Calcio entre tratamientos (H (g.l. 3) = 0.57, P (α0.05) = 0.9038).
Raman: Se detectaron tres picos para aragonita (206, 705 y 1 085 cm-1) y Calcita (282, 711 y 1 085 cm-1) (Fig. 6). El pico de 1 085 cm-1 mostró mayor intensidad en 31 °C (pH 7.6 = 3 338.77 ± 1 046.86 a.u. y pH 8.1 = 1 259.15 ± 608.35 a.u.). El segundo pico (de 705 a 711 cm-1) registró más intensidad en pH 7.6 (28 °C = 259.15 ± 121.02 a.u. y 31 °C = 606.50 ± 111.80 a.u.). Se detectó aragonita en todos los tratamientos, con mayor ocurrencia en 28 °C (100 %, independientemente del pH). La intensidad de aragonita (206 cm-1) fue mayor en pH 8.1 (Control: 1 039.54 ± 780.26 a.u. y 31 °C: 1 228.87 ± 0.00 a.u.) (Fig. 7). Se registró calcita (282 cm-1) solo en larvas cultivadas a 31 °C (pH 8.1 = 80 %; pH 7.6 = 17 %), la intensidad de este carbonato fue mayor en pH 8.1 (206.04 ± 122.33 a.u.) (Fig. 7).
DISCUSIÓN
El calentamiento y acidificación del océano podrían tener efectos severos en organismos calcificadores (Gazeau et al., 2007). Las larvas de organismos marinos, al ser más vulnerables a cambios ambientales, son empleadas como indicador del efecto del incremento de la temperatura y reducción del pH a través del estudio de las tasas de eclosión, desarrollo, crecimiento, mortalidad y calcificación en larvas de invertebrados (Sale et al., 2010; Talmage & Gobler, 2009; Ushakova, 2003).
En general, los estudios de cambio climático se han desarrollado con especies de ambientes fríos a templados dentro de las cuales destacan Haliotis sp. y Mytilus sp. (Gaylord et al., 2011; Gazeau et al., 2010, 2011; Miller, Reynolds, Sobrino, & Riedel, 2009; Sunday, Crim, Harley, & Hart, 2011), siendo escasas las investigaciones con moluscos de ambientes tropicales.
Al evaluar la tasa de eclosión Parker, Ross, y O’Connor (2010), en tratamientos combinados de temperatura (18, 22, 26 y 30 °C) y pH (8.2, 8.0, 7.9 y 7.8), para Saccostrea glomerata y C. gigas, obtuvieron 100 % de larva D en condiciones normales (26 °C - pH 8.2); sin embargo, observaron menor proporción de larvas al disminuir la temperatura y el pH (18 °C - pH 7.8: C. gigas = 35 % y S. glomerata = 58 %). Por otra parte, Gazeau et al. (2011) reportan menor porcentaje de eclosión para embriones de C. gigas al disminuir el pH de 8.03 a 7.41 (de 95 a 20 %). En el caso de S. gigas, la tasa de eclosión no se vio influenciada por ningún parámetro, sin embargo, se puede determinar que el tiempo de incubación se reduce, al incrementar la temperatura, independientemente del pH.
Aldana, Lucas, Brulé, Salguero, y Rendon (1989) reportan el crecimiento para Strombus costatus en función de la temperatura, y registraron tallas de 812 µm en 28 °C y 685 µm en 24 °C (para larvas de 31 días). En cuanto a la acidificación, se reportaron menores tasas de crecimiento para larvas de Mercenaria mercenaria, Argopecten irradians y C. virginica cultivadas en pH 7.5 (Talmage & Gobler, 2009), en el presente estudio observamos que si la temperatura se mantiene constante (28 °C) y el pH disminuye, el crecimiento de S. gigas es menor, en contraparte, al aumentar la temperatura, el efecto del pH es contrarrestado, los que beneficia el crecimiento del organismo.
Se conoce que el incremento de la temperatura y la disminución del pH tienen un efecto negativo en las tasas de mortalidad larval de moluscos. Aldana Aranda et al. (1989) registraron mortalidades de 83 a 90 % para S. costatus ante incrementos de temperatura (de 24 a 28 °C, respectivamente). En cuanto a la acidificación, Talmage y Glober (2009, 2010) observaron que a menor pH (7.5) la mortalidad diaria incrementa para M. mercenaria (de 1.8 a 6.7 %) y A. irradians, (de 3.8 a 7.8 %). Resultados similares fueron presentados por Kurihara, Kato, e Ishimatsu (2007) y Gazeau et al. (2011) para C. gigas, donde la mortalidad aumentó de 20 a 90 % al reducirse el pH (de 8.1 a 7.6, respectivamente). Los resultados de este estudio mostraron un mayor porcentaje de mortalidad hacia el día 20 de desarrollo, observando que las larvas cultivadas a 31 °C fueron más afectadas a nivel fisiológico por la temperatura.
La acidificación del medio marino produce baja saturación de aragonita y calcita, lo cual afecta directamente a organismos calcificadores al reducirse la tasa de incorporación de calcio en sus concha (Talmage & Gobler, 2009). Kurihara et al. (2007) y Gazeau et al. (2011), señalaron que la tasa de calcificación en larvas de C. gigas disminuye de 72 a 30 % ante menor pH (de 8.1 a 7.4, respectivamente); contrariamente, en este estudio se observó que la proporción de calcio en el estadio planctónico de S. gigas, no mostró variación en ninguna condición experimental.
Weiss, Tuross, Addadi, y Weiner (2002) informaron que en estadios iniciales, la concha de M. mercenaria y C. gigas se compone de carbonato de calcio amorfo, sustituyéndose por aragonita al final de la vida larvaria. Auzoux-Bordenave et al. (2010) presentaron una revisión para siete especies de moluscos, entre ellas Ostrea edulis y Pinctada fucata, en la que señalaron que la protoconcha se compone por aragonita y trazas de carbonato amorfo, mientras que la teleoconcha por calcita. En este estudio, las conchas de las larvas de 30 días (premetamórficas) de todos los tratamientos, presentaron aragonita, mientras que calcita solo se observó en organismos cultivados a 31 °C. Ambos polimorfos de carbonato de calcio mostraron menores registros en función del pH.
Se sabe que el incremento de temperatura puede aumentar los efectos negativos de la acidificación del océano, cuando pasan el umbral de tolerancia térmica, ya que puede alterar la membrana celular y modificar la homeóstasis ácido-base; también puede influir en la actividad de las proteínas implicadas en la regulación ácido-base (O’Donnell, Hammond, & Hofmann, 2008; Padilla-Gamiño, Kelly, Evans, & Hofmann, 2013; Pörtner, 2010). De acuerdo con Le Moullac et al. (2016), las larvas pueden adaptarse a condiciones de acidificación, siendo la temperatura el principal estresor en el desarrollo larvario. En este estudio se observó que S. gigas podría adaptarse al calentamiento y acidificación de los océanos, sin embargo, los procesos de biomineralización de la concha fueron afectados, razón por la cual las conchas de S. gigas mostraron mayor proporción de aragonita (carbonato presente en las fases tempranas del desarrollo larvario), retardándose así el proceso de formación de calcita. Aunque la aragonita es el polimorfo de carbonato de calcio menos estable, su presencia en todos los tratamientos podría atribuirse a la dureza de este carbonato, por lo cual, se detectó aún en condiciones de pH menor, no obstante con menores intensidades, lo cual, se traduce en la formación de conchas con menor densidad o estructuras menos compactas, lo que podría generar repercusiones a nivel fisiológico.
Si bien, se determinaron las fases (presencia/ausencia) e intensidades de aragonita y calcita de la concha de S. gigas, se debe considerar realizar un análisis químico a fin de conocer la variación o estratificación de estos dos carbonatos de calcio, a fin de conocer a fondo los efectos de la acidificación, empleando diferentes ambientes de pH, a través de la vida larvaria de S. gigas.
AGRADECIMIENTOS
Agradecemos al Consejo Nacional de Ciencia y Tecnología por financiar el presente estudio (Proyecto: “El caracol rosa como indicador de Cambio climático en el Caribe: Calentamiento y Acidificación Oceánica: CB-2011-01/181329). Al Centro de Investigación y de Estudios Avanzados del IPN, Unidad Mérida, en especial a Ana Ruth Cristóbal Ramos y José Bante Guerra por su apoyo técnico en la obtención de datos mediante EDX y Microscopía RAMAN. A Geonel Rodríguez Gattorno por su contribución en el análisis RAMAN.
RESUMEN
El aumento de las emisiones de CO2 produce calentamiento y reducción del pH en los océanos, lo cual puede generar efectos negativos en muchos organismos marinos, particularmente en aquellos con estructuras calcáreas (i.e. moluscos), afectando principalmente a sus estadios larvarios (Le Moullac et al., 2016), en este sentido, se estudió a S. gigas, gasterópodo de importancia comercial en el mar Caribe, con el fin de conocer el efecto de la temperatura y la acidificación en el desarrollo, el crecimiento, la mortalidad y la calcificación durante su fase larvaria. Se realizó un cultivo larvario empleando cuatro tratamientos de temperatura y pH (control = 28 °C - pH 8.1, T1 = 28 °C pH 7.6, T2 = 31 °C pH 8.1 y T3 = 31 °C - pH 7.6) por triplicado. La eclosión se registró al inicio del experimento (No. de huevos – No. de larvas nacidas), por otra parte, el desarrollo de órganos, el crecimiento de la concha y la mortalidad se evaluaron a través del tiempo. La calcificación fue estudiada mediante análisis EDX y Raman para larvas de 30 días de edad. Se observó que el desarrollo y el crecimiento de órganos fue mayor a 31 °C (talla inicial = 230 ± 4.12 a 313.27 ± 11.34 µm, talla final = 829.50 ± 11.33 a 1 054.50 ± 11.13 µm; para T1 y T2, respectivamente), mismo patrón se presentó para el tiempo de eclosión (18 hr) y la tasa de mortalidad (~ 57 %). La proporción de calcio (% wt) fue similar entre tratamientos (de 34.37 ± 10.05 a 37.29 ± 16.81 % wt). El análisis Raman mostró aragonita para las conchas de todas las condiciones experimentales, con valores más altos en el control (1 039.54 ± 780.26 a.u.). Calcita solo se detectó en los tratamientos de 31 °C (174.56 ± 127.19 a.u.), mientras que, a menor pH (7.6), la intensidad de aragonita y calcita fue menor. En conclusión, S. gigas podría adaptarse a escenarios futuros de temperatura y acidificación, sin embargo, puede verse afectado durante los procesos de biomineralización de la concha.
Palabras clave: desarrollo, temperatura, acidificación, calcificación, Strombus gigas.
REFERENCIAS
Aldana-Aranda, D., Lucas, A., Brule, T., Salguero, E., & Rendon, F. (1989). Effects of temperature, algal food, feeding rate and density on the larval growth of the milk conch (Strombus costatus) in Mexico. Aquaculture, 76(3), 361-371. http://doi.org/10.1016/0044-8486(89)90087-2
Auzoux-Bordenave, S., Badou, A., Gaume, B., Berland, S., Helléouet, M. N., Milet, C., & Huchette, S. (2010). Ultrastructure, chemistry and mineralogy of the growing shell of the European abalone Haliotis tuberculata. Journal of Structural Biology, 171(3), 277-290. http://doi.org/10.1016/j.jsb.2010.05.012
Brito-Manzano, N., Aldana-Aranda, D., de la Cruz-Lázaro, E., & Estrada-Botello, M. A. (2006). Organogénesis larvaria de Strombus gigas (Mesogastropoda: Strombidae) en el Arrecife Alacranes durante el período máximo de su época reproductiva. Revista Universidad y Ciencia, 22(1), 75-82.
Caldeira, K., & Wickett, M. E. (2003). Oceanography: Anthropogenic carbon and ocean pH. Nature, 425(6956), 365-365. http://doi.org/10.1038/425365ª
Chollett, I., Mumby, P. J., Müller-Karger, F. E., & Hu, C. (2012). Physical environments of the Caribbean Sea. Limnology and Oceanography, 57(4), 1233-1244. http://doi.org/10.4319/lo.2012.57.4.1233
Enriquez-Diaz, M. R., Volland, J. M., Chavez-Villegas, J. F., Aldana-Aranda, D., & Gros, O. (2015). Development of the planktotrophic veligers and plantigrades of Strombus pugilis (Gastropoda). Journal of Molluscan Studies, eyv011. http://doi.org/10.1093/mollus/eyv011
Gaylord, B., Hill, T. M., Sanford, E., Lenz, E. A., Jacobs, L. A., Sato, K. N., … Hettinger, A. (2011). Functional impacts of ocean acidification in an ecologically critical foundation species. Journal of Experimental Biology, 214(15), 2586-2594. http://doi.org/10.1242/jeb.055939
Gazeau, F., Gattuso, J. P., Dawber, C., Pronker, A. E., Peene, F., Peene, J., … Middelburg, J. J. (2010). Effect of ocean acidification on the early life stages of the blue mussel Mytilus edulis. Biogeosciences, 7(7), 2051-2060. http://doi.org/10.5194/bg-7-2051-2010
Gazeau, F., Gattuso, J. P., Greaves, M., Elderfield, H., Peene, J., Heip, C. H. R., & Middelburg, J. J. (2011). Effect of Carbonate Chemistry Alteration on the Early Embryonic Development of the Pacific Oyster (Crassostrea gigas). PLOS One, 6(8), e23010. http://doi.org/10.1371/journal.pone.0023010
Gazeau, F., Quiblier, C., Jansen, J. M., Gattuso, J. P., Middelburg, J. J., & Heip, C. H. R. (2007). Impact of elevated CO2 on shellfish calcification. Geophysical Research Letters, 34(7), L07603. http://doi.org/10.1029/2006GL028554
Gledhill, D. K., Wanninkhof, R., Millero, F. J., & Eakin, M. (2008). Ocean acidification of the Greater Caribbean Region 1996–2006. Journal of Geophysical Research: Oceans, 113(C10), C10031. http://doi.org/10.1029/2007JC004629
Intergovernmental Panel on Climate Change. (2007). Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press. Retrieved from https://www.ipcc.ch/pdf/assessment-report/ar4/wg2/ar4_wg2_full_report.pdf
Kurihara, H., Kato, S., & Ishimatsu, A. (2007). Effects of increased seawater pCO2 on early development of the oyster Crassostrea gigas. Aquatic Biology, 1(1), 91-98. http://doi.org/10.3354/ab00009
Le Moullac, G., Soyez, C., Latchere, O., Vidal-Dupiol, J., Fremery, J., Saulnier, D., … Gueguen, Y. (2016). Pinctada margaritifera responses to temperature and pH: Acclimation capabilities and physiological limits. Estuarine, Coastal and Shelf Science. http://doi.org/10.1016/j.ecss.2016.04.011
Magaña Rueda, V., Graizbord, B., Buenfil Friedman, J., & Gómez Mendoza, L. (2009). Escenarios de cambio climático y tendencias en la zona del Golfo de México. In J. Buenfil (Ed.). Adaptación a los impactos del cambio climático en los humedales costeros del Golfo de México (Vol. 2, pp. 569-673). México D.F.: Secretaría de Medio Ambiente y Recursos Naturales / Instituto Nacional de Ecología (INE-SEMARNAT). Retrieved from http://www2.inecc.gob.mx/publicaciones/libros/610/cap6.pdf
Miller, A. W., Reynolds, A. C., Sobrino, C., & Riedel, G. F. (2009). Shellfish Face Uncertain Future in High CO2 World: Influence of Acidification on Oyster Larvae Calcification and Growth in Estuaries. PLOS One, 4(5), e5661. http://doi.org/10.1371/journal.pone.0005661
Nehrke, G., Poigner, H., Wilhelms-Dick, D., Brey, T., & Abele, D. (2012). Coexistence of three calcium carbonate polymorphs in the shell of the Antarctic clam Laternula elliptica. Geochemistry, Geophysics, Geosystems - Wiley Online Library, 13(5), Q05014. http://doi.org/10.1029/2011GC003996
O’Donnell, M. J., Hammond, L. M., & Hofmann, G. E. (2008). Predicted impact of ocean acidification on a marine invertebrate: elevated CO2 alters response to thermal stress in sea urchin larvae. Marine Biology, 156(3), 439-446. http://doi.org/10.1007/s00227-008-1097-6
Orr, J. C., Fabry, V. J., Aumont, O., Bopp, L., Doney, S. C., Feely, R. A., … Yool, A. (2005). Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms. Nature, 437(7059), 681-686. http://doi.org/10.1038/nature04095
Padilla-Gamiño, J. L., Kelly, M. W., Evans, T. G., & Hofmann, G. E. (2013). Temperature and CO2 additively regulate physiology, morphology and genomic responses of larval sea urchins, Strongylocentrotus purpuratus. Proceedings of the Royal Society of London B: Biological Sciences, 280(1759), 20130155. http://doi.org/10.1098/rspb.2013.0155
Parker, L. M., Ross, P. M., & O’Connor, W. A. (2010). Comparing the effect of elevated pCO2 and temperature on the fertilization and early development of two species of oysters. Marine Biology, 157(11), 2435-2452. http://doi.org/10.1007/s00227-010-1508-3
Pörtner, H. O. (2010). Oxygen- and capacity-limitation of thermal tolerance: a matrix for integrating climate-related stressor effects in marine ecosystems. Journal of Experimental Biology, 213(6), 881-893. http://doi.org/10.1242/jeb.037523
Roller, R. A., & Stickle, W. B. (1989). Temperature and salinity effects on the intracapsular development, metabolic rates, and survival to hatching of Thais haemastoma canaliculata (Gray) (Prosobranchia:Muricidae) under laboratory conditions. Journal of Experimental Marine Biology and Ecology, 125(3), 235-251. http://doi.org/10.1016/0022-0981(89)90099-3
Sale, P. F., Van Lavieren, H., Ablan, M. C., Atema, J., Butler, M., Fauvelot, C., ... & Stewart, H. L. (2010). Conservando la Conectividad de los Arrecifes: Guía Para los Administradores de las Áreas Marinas Protegidas. Grupo de Trabajo de Conectividad, Programa de Investigación Dirigido a los Arrecifes de Coral y a la Creación de Capacidades para la Gestión, UNU-INWEH. Retrieved from http://ccres.net/images/uploads/publications/2/crtr_connectivityhandbook_web_spanish.pdf
Sunday, J. M., Crim, R. N., Harley, C. D. G., & Hart, M. W. (2011). Quantifying rates of evolutionary adaptation in response to ocean acidification. PloS One, 6(8), e22881. http://doi.org/10.1371/journal.pone.0022881
Talmage, S. C., & Gobler, C. J. (2009). The effects of elevated carbon dioxide concentrations on the metamorphosis, size, and survival of larval hard clams (Mercenaria mercenaria), bay scallops (Argopecten irradians), and Eastern oysters (Crassostrea virginica). Limnology and Oceanography, 54(6), 2072-2080. http://doi.org/10.4319/lo.2009.54.6.2072
Talmage, S. C., & Gobler, C. J. (2010). Effects of past, present, and future ocean carbon dioxide concentrations on the growth and survival of larval shellfish. Proceedings of the National Academy of Sciences of the United States of America, 107(40), 17246-17251. http://doi.org/10.1073/pnas.0913804107
Ushakova, O. O. (2003). Combined effect of salinity and temperature on Spirorbis spirorbis L. and Circeus spirillum L. larvae from the White Sea. Journal of Experimental Marine Biology and Ecology, 296(1), 23-33.
Watson, S. A., Southgate, P. C., Tyler, P. A., & Peck, L. S. (2009). Early Larval Development of the Sydney Rock Oyster Saccostrea glomerata Under Near-Future Predictions of CO2-Driven Ocean Acidification. Journal of Shellfish Research, 28(3), 431-437. http://doi.org/10.2983/035.028.0302
Weiss, I. M., Tuross, N., Addadi, L., & Weiner, S. (2002). Mollusc larval shell formation: amorphous calcium carbonate is a precursor phase for aragonite. Journal of Experimental Zoology, 293(5), 478-491. http://doi.org/10.1002/jez.90004
