Introducción

La interacción genotipo-ambiente provoca variaciones en la expresión genética de los cerdos. Esto puede generar diferencias en cada etapa fisiológica, influidas por la fuente alimenticia utilizada, por lo que resulta necesario conocer estas variaciones para potencializar la producción de esta especie animal (^{Albuquerque et al., 2020}; ^{Malgwi et al., 2022}; ^{Óvilo et al., 2014}; ^{Shang-Qiao et al., 2019}; ^{Tao et al., 2017}). Con la secuenciación directa del ARN (RNA-Seq), se han identificado genes relacionados con los procesos biológicos, que han permitido comprender los mecanismos fisiológicos y el efecto que tiene sobre estos mecanismos la fuente nutricional empleada (^{Benítez et al., 2021}; ^{Muñoz et al., 2018}; ^{Núñez et al., 2021}; ^{Shang-Qiao et al., 2019}).

El cerdo pasa por varias etapas fisiológicas durante su engorda, en las que se producen cambios en su crecimiento, ocasionados por la expresión de genes. Es importante conocer los procesos biológicos que se ven afectados durante este periodo, con el fin de identificar los efectos, mejorar la producción y modular la calidad de la carne (^{Wang et al., 2017}; ^{Wang et al., 2020}). Además, un conocimiento profundo de la genética del cerdo podría permitir la selección de ingredientes según su carga genética (^{Malgwi et al., 2022}). Por tal motivo, el objetivo de esta investigación fue identificar los genes con expresión genética diferencial en los procesos biológicos del cerdo, en los tejidos Longissimus dorsi (músculo dorsal ancho) e hígado, por medio del análisis de transcriptoma durante las etapas de crecimiento (55 ± 1,05 kg) y la engorda final (101 ± 7,8 kg).

Materiales y métodos

Resultados

De acuerdo con los resultados obtenidos en la comparación entre grupos y tipo de tejido (Cuadro 2), en el hígado se identificó la mayor cantidad de genes con expresión genética diferencial: 81 genes para el grupo de crecimiento (GC) y 286 para el grupo de engorda final (GF). En el tejido L. dorsi del GC, se observaron 35 genes con expresión y 25 sin expresión. En el hígado sucede lo contrario, la mayor cantidad de genes (286) con expresión se presentó en el GF, y 25 genes no se expresaron en el GC.

La mayor cantidad de genes sin expresión se observó en el tejido hígado: 2175 genes en el GC y 2646 en el GF (Cuadro 3). Estos mismos genes fueron los que no mostraron expresión en el tejido L. dorsi. También se identificaron 1913 genes sin expresión en el tejido hígado del GC, así como 1660 genes que se expresaron en el tejido L. dorsi del GF.

En el cromosoma 6, al comparar la expresión genética diferencial entre los tejidos hígado y L. dorsi, se observó que en ambos grupos (GC y GF) existió el mayor número de genes con expresión (205 y 258, respectivamente). En el cromosoma 2, ocurrió algo similar a lo descrito en el cromosoma 6. En cuanto al cromosoma 1, la mayor cantidad de genes con expresión genética diferencial baja se encontró en el tejido hígado, tanto en el GC como en el GF, con 183 y 157 genes, respectivamente (Cuadro 4).

Los genes con expresión genética diferencial extrema y mayores cambios de pliegue Log2 en L. dorsi e hígado se presentan en el Cuadro 5. Se identificaron con expresión los genes FUT1, SESN2 y FGF21 en el GC. También se identificaron genes con expresión en el mismo grupo (GC) relacionados con detoxificación (GPX2, MT4 y SESN2). En el grupo de crecimiento, se identificaron los genes NR4A3, PDK4, PER1 y PTPRO, los cuales están involucrados en funciones biológicas como procesos del sistema inmune y ritmo circadiano. Sin embargo, en los cerdos del GC no se identificaron genes con expresión que pudieran estar involucrados con las siguientes funciones biológicas: respuesta al estrés, sistema inmune, locomoción, regulación de lipoproteínas, ritmo circadiano, comportamiento y detoxificación. Los genes IHH y MYL7 fueron los más involucrados en funciones biológicas con expresión en el GC frente al GF.

En el hígado (Cuadro 6), se observó una mayor expresión genética diferencial extrema y mayores cambios de pliegue Log2, mientras que en el GC no hubo expresión. En este sentido, cuando no se observó expresión en el hígado del GC, sí hubo expresión en el GF. Los genes MFSD2A, LIPG, THBS1, TGFB2, LTF y APOA4 se identificaron como los más involucrados en procesos biológicos.

Discusión

Se ha señalado en la literatura que el músculo L. dorsi tiene una mayor expresión de genes lipogénicos, lo cual es fundamental para la modulación de la grasa intramuscular. Sin embargo, también se reconoce que el hígado tiene una mayor función biológica y, en general, una mayor expresión de genes, por lo que resulta un órgano esencial en todos los procesos metabólicos durante la vida del cerdo (^{Muñoz et al., 2021}). Otros estudios han demostrado resultados similares, reconociendo que existe una mayor expresión de genes en el hígado y que este órgano es el más sensible al cambio durante las diferentes etapas de desarrollo del cerdo. Además, se ha señalado que el hígado es muy influido por las fuentes alimenticias utilizadas durante el proceso de producción (^{Benítez et al., 2015}; ^{Muñoz et al., 2021}).

Entre los genes con expresión en el GC en L. dorsi, se identificó al gen FUT1, que es importante en la glucosilación intestinal y desempeña un rol fundamental en la salud digestiva del cerdo, según lo reportado por ^{Hesselager et al. (2016)}. Estos autores sugieren que la activación de este gen podría mejorar la salud intestinal de los cerdos, lo cual tendría un impacto positivo en la producción animal, pues podría darles cierta resistencia frente a patógenos como Escherichia coli, que tiene una presencia significativa en esta especie. En este contexto, mejorar la salud intestinal, y por ende la producción en general, sería de gran relevancia. Además, ^{Hesselager et al. (2016)} mencionan que una alternativa es la cría selectiva de animales que presenten una alta expresión del gen FUT1, lo que potenciaría aún más la mejora.

El gen SESN2 también se relaciona con la calidad de la carne y a la regulación de la grasa intramuscular (^{Puig-Oliveras et al., 2016}; ^{Wang et al., 2015}). Por su parte, el gen FGF21 se asocia al crecimiento, pero también se ha reportado que es un supresor de la adipogénesis en células grasas intramusculares (^{Wang et al., 2015}). Otros autores como ^{Wang et al. (2017)}, señalan que después de 120 días de edad hay una mayor acumulación de grasa, por lo que en el músculo L. dorsi del GC fue mayor la expresión de genes involucrados en el crecimiento, lo que concuerda con los hallazgos previos. La edad de los cerdos que se utilizaron en este estudio estaba dentro de ese rango y sí se observó una mayor expresión de genes.

En el músculo L. dorsi, los genes sin expresión en el GC y con expresión en el GF fueron NR4A3, PDK4, PER1 y PTPRO. El gen NR4A3 está involucrado en la respuesta inflamatoria y presenta expresión en cerdos con mucha grasa intramuscular (^{Wang et al., 2021}). Este gen tiene mucha similitud con el PDK4, que se asocia a la grasa intramuscular y a la calidad de la carne. Por su parte, el gen PER1 participa en la regulación del ciclo de alimentación, lo que permite la ingesta de comida en los periodos más convenientes para su metabolismo; además, regula el ciclo del sueño e influye en la actividad locomotora y en el comportamiento (^{Figueiredo et al., 2017}). En cuanto al gen PTPRO, este es un receptor de la tirosina fosfatasa que se involucra en procesos del sistema inmune y en el desarrollo de la estructura anatómica y la morfogénesis (Cuadro 5).

Con base en la expresión de los genes identificados en el músculo L. dorsi, se puede diferenciar que en el GC suceden procesos biológicos más asociados al crecimiento, mientras que en el GF los cambios están más relacionados con la inmunidad, el crecimiento y la calidad de la carne.

El gen IHH se identificó con expresión en el hígado (Cuadro 6). Este gen codifica proteínas de la familia Hedgehog, moléculas esenciales de señalización que regulan una variedad de procesos asociados al desarrollo, como el crecimiento y la morfogénesis (^{The Human Gene Database (THGD), n. d.}). En cuanto al gen MYL7, conocido por codificar proteínas, se ha reportado que entre sus vías relacionadas se encuentran la respuesta inmune, la remodelación del citoesqueleto y la regulación del citoesqueleto de actina por Rho GTPasas (^{THGD, n. d.).} Este funciona de manera específica para el músculo, en su etapa de diferenciación. Además, su baja expresión se asocia con distrofia muscular, miocardiopatías, entre otros problemas (^{Tamiyakul et al., 2020}); por lo tanto, la expresión de estos genes puede representar una mejor respuesta productiva e inmune.

En el hígado, sin expresión en el GC, se identificó al gen MFSD2A, el cual está involucrado en el metabolismo de nutrientes y se localiza en el cromosoma 6. Este gen se asocia al crecimiento con alta frecuencia genética en la raza Duroc (^{Lemus-Flores et al., 2020}). También se ha reportado que está relacionado con el transporte de moléculas, el metabolismo de los lípidos y la homeostasis (^{Wong et al., 2022}; ^{Zhang et al., 2019}). Por su parte, el gen LIPG actúa en el metabolismo de lípidos en el músculo L. dorsi (^{Tao et al., 2017}). En cuanto al gen THBS1, este se relaciona con glicoproteínas y formación de heparina, la coagulación de la sangre y la formación de vasos sanguíneos nuevos; si existen alteraciones negativas, puede verse asociado a una alteración en la proliferación celular. Al igual que TGFB2, está presente en la morfogénesis y formación de la estructura anatómica (^{THGD, n. d.}; ^{The Human Protein Atlas, n. d.}). Por último, los genes TGFB2 y LTF se relacionan con la respuesta inflamatoria, mientras que el gen APOA4 se vincula con la detoxificación, la modulación de glucosa y el metabolismo de lípidos y lipoproteínas (^{Wilkinson et al., 2010}; ^{Wong et al., 2022}).

Conclusiones

Al comparar los cerdos en etapa de crecimiento con los de etapa de finalización en cuanto a la expresión genética, el primer grupo presentó genes relacionados con el crecimiento, mientras que el segundo grupo mostró expresión de genes asociados a inmunidad, crecimiento, metabolismo de nutrientes y lípidos, lipoproteínas y detoxificación. Estas diferencias fueron consistentes independientemente del tipo de tejido analizado.

Conflicto de interés

El manuscrito fue preparado y revisado con la participación de todos los autores, quienes declaran que no existe ningún conflicto de intereses que ponga en riesgo la validez de los resultados presentados.