SciELO - Scientific Electronic Library Online

 
vol.8Composición de especies e ingreso económico por la pesca ribereña de Moluscos en la Costa Grande de Guerrero MéxicoLas Nanopartículas de SiO2 Mejoran la Absorción de Nutrientes en Plantas de Tomate Desarrolladas en Presencia de Arsénico índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista bio ciencias

versión On-line ISSN 2007-3380

Revista bio ciencias vol.8  Tepic  2021  Epub 04-Oct-2021

https://doi.org/10.15741/revbio.08.e1052 

Artículos originales

Frecuencia del alelo menor en predicción genómica para características de crecimiento en bovinos Suizo Europeo

M. Z. Trujano-Chavez1 

J. E. Valerio-Hernández1 

R. López-Ordaz1 

A. Ruíz-Flores1  * 

1 Universidad Autónoma Chapingo, Posgrado en Producción Animal, Departamento de Zootecnia, Estado de México, México.


Resumen

En el estudio se compararon umbrales de frecuencia del alelo menor (MAF) durante el control de calidad de los genotipos, utilizando polimorfismos de nucleótido simple (SNPs) para predecir valores genómicos de características de crecimiento en bovinos Suizo Europeo. El análisis se realizó con 28,973 y 18,994 registros de pesos al nacimiento y destete, 12,835 polimorfismos de nucleótido simple y 300 animales. Las características se analizaron bajo el enfoque de evaluación genómica en un solo paso. Los criterios de comparación fueron: 1) jerarquización de valores genómicos: se estimaron coeficientes de correlación Spearman y Pearson entre valores genómicos predichos con frecuencia estándar del alelo menor (0.05) y niveles menores; 2) habilidad de predicción de los modelos: se realizó validación cruzada con cuatro repeticiones y se obtuvo el promedio de los coeficientes de correlación de los fenotipos reales contra predichos; y 3) coeficiente de regresión simple: los valores genómicos obtenidos con frecuencia del alelo menor estándar e inferiores, fueron las variables independiente y dependiente, respectivamente. El efecto del nivel de frecuencia del alelo menor en la jerarquización de los valores genómicos no fue importante (correlaciones mayores que 0.99, <2.2x10-16). La habilidad predictiva del modelo fue similar: 0.7 promedio. Los coeficientes de regresión fueron cercanos a uno (0.99, <2x10-16). Los resultados sugieren que una frecuencia del alelo menor entre 0.0 y 0.05, no altera las predicciones de valores genómicos en poblaciones pequeñas. Sin embargo, estos resultados son preliminares y susceptibles a cambios si el estudio se repite con un número de marcadores y animales mayor.

Palabras clave: Polimorfismo de nucleótido simple; frecuencia del alelo menor; control de calidad de genotipos; predicción de valores genómicos; Pardo Suizo

Abstract

The present research aimed to compare the effect of minor allele frequency (MAF) threshold during the quality control process of genotypes, using single nucleotide polymorphisms (SNPs) to predict genomic values of growth traits in Braunvieh cattle. The analysis was performed including 28,973 and 18,994 phenotypic records of birth and weaning weights, 12,835 single nucleotide polymorphisms and 300 animals. The traits were analyzed with the single-step genomic evaluation approach. The criteria for comparison of the minor allele frequency thresholds compared were: 1) ranking of genomic values: Spearman and Pearson correlation coefficient estimates were obtained for genomic values predicted with standard (0.05) and lower levels of minor allele frequency; 2) prediction ability of models: cross- validation with four replicates was performed, average of correlation between real and predicted phenotypes was obtained; and 3) regression coefficient estimates: the genomic values obtained with standard and lower levels of minor allele frequency were used as independent and dependent variables, respectively. The effect of minor allele frequency threshold on the ranking of genomic values was not important (correlation estimates greater than 0.99, <2.2x10-16). The prediction ability model was similar: around 0.7. the regression coefficient estimates were close to one (0.99, <2x10-16). The results suggest that a minor allele frequency between 0.0 and 0.05, does not influence the genomic value prediction in small populations. However, these results should be regarded as preliminary and susceptible to changes if the study is repeated with a greater number of markers and animals.

Keywords: Single nucleotide polymorphism; minor allele frequency; genotype quality control; prediction of genetic values; Braunvieh

Introducción

El descubrimiento de cantidades grandes de marcadores polimorfismos de un solo nucleótido (SNPs); la tecnología de rendimiento alto para genotipar animales en miles de SNPs de manera rentable y los avances en metodologías estadísticas para análisis, ayudaron a generalizar el uso de selección genómica (SG) en el mundo (Meuwissen et al., 2016). Actualmente se estudia el efecto en conjunto de marcadores moleculares en la predicción de valores genéticos (VG) de características de interés económico en los animales (De los Campos et al., 2013).

El éxito del aprovechamiento de marcadores SNP en SG reside en lograr predicciones precisas de los VG para las características estudiadas. Según Meuwissen et al. (2016) se requieren individuos con información fenotípica confiable y con la mayor cantidad posible de SNPs. Para mejorar la precisión de las predicciones de los VG, se han propuesto variantes de la metodología innovadora de (Meuwissen et al., 2001), como los enfoques semiparamétricos (Gianola et al., 2006) o no paramétricos (Gianola et al., 2009). Además, continúan surgiendo estudios sobre modificaciones y aportaciones metodológicas, como el uso del agrupamiento de K-means para validación cruzada contra la metodología común (Saatchi et al., 2011). La evaluación genómica en un solo paso (ssGBLUP) es un análisis estadístico novedoso en SG que permite utilizar información genómica y de pedigrí, con precisión superior en comparación con otros métodos de varios pasos (Wang et al., 2014). El análisis ssGBLUP requiere la construcción de varias matrices: la matriz A (de relaciones genéticas aditivas), la matriz G (de relaciones genómicas) y la matriz H (de relaciones entre animales genotipados y no genotipados) (Pértile et al., 2016).

Independientemente del uso de metodologías estadísticas para el análisis de la información, se requiere que ésta sea de buena calidad. Anderson et al. (2010) detallan la metodología para realizar un buen control de calidad del genotipo (GQC), una de las etapas es la eliminación de los marcadores con frecuencia del alelo menor o Minor Allele Frequency (MAF) muy baja. La MAF se define como todos aquellos alelos, en la secuencia de ADN, cuya frecuencia varía de 0.01 a 0.2 (Chouraki & Seshadri, 2014). La MAF también se ha relacionado con la prevalencia de enfermedades genéticas raras (Kido et al., 2018). Estos autores reportaron que 50 % de MAF se asocia con enfermedades complejas, como Alzheimer, diabetes, cáncer y esquizofrenia.

En la edición de la información genética, al eliminar los alelos de menor frecuencia se pretende aumentar la frecuencia de los alelos que tienen un efecto importante en las características en estudio. Para asegurar la calidad del genotipo en las predicciones de valores genéticos se eliminan los marcadores con MAF menor que 5 % o hasta 1 %. Por ejemplo, Minozzi et al. (2012) trabajaron en el GQC de genotipos bovinos, eliminando marcadores con MAF menor que 1 %. Anderson et al. (2010) indican que el no eliminar los alelos de baja frecuencia provoca el uso de información falsa, por dos razones 1) las asociaciones observadas en estos SNPs son pequeñas porque son impulsadas por los genotipos de pocos individuos; y 2) pueden provenir de errores de genotipado en los marcadores que en realidad son monomórficos en la población.

Para características de crecimiento en bovinos no existe claridad completa del efecto del umbral de MAF en la precisión de las predicciones en poblaciones pequeñas, aunque se entiende que al aumentar la calidad del genotipo se incrementará también la precisión en las predicciones de los VG, aun cuando utilizar valores altos de MAF implique eliminar gran parte de la información de marcadores SNP. Por lo anterior, el objetivo del presente estudio fue comparar el efecto de diferentes umbrales MAF en la edición de GQC utilizando SNPs, para la predicción genómica en una sola etapa de características de crecimiento en bovinos Suizo Europeo.

Material y Métodos

Origen de la información

La información fenotípica provino de la base de datos de la Asociación Mexicana de Criadores de Ganado Suizo de Registro. Se tomaron los registros genealógicos, identificación y los pesos al nacimiento (BW) y al destete (WW) de los animales. El pedigrí lo conformaron 184,788 animales nacidos entre 1901 y 2016.

Para la información genómica se colectaron muestras de pelo de 300 animales, 236 hembras y 64 machos, seleccionados aleatoriamente de ranchos ubicados en los estados de Colima, Jalisco y Veracruz, México. Posteriormente las muestras se enviaron a genotipar a GeneSeek (Lincoln, NE, USA). El chip usado fue Genomic Profile Bovine LD, de 30,000 y 50,000 marcadores SNP para 150 animales cada uno.

Fenotipos

Para BW, se eliminaron los registros de individuos con BW fuera del intervalo promedio ± tres desviaciones estándar, sin información de hato o edad de la madre. Los grupos contemporáneos se definieron combinando los efectos de sexo, hato, año y época de nacimiento. Las épocas de nacimiento se definieron considerando el día juliano: del día 80-171, primavera; 172-264, verano; 265-354, otoño; 355-366, invierno; y 1-79, invierno. La información de animales en grupos contemporáneos con menos de dos individuos y animales en grupos de tres o más individuos con varianza cero también se eliminó. Finalmente se obtuvieron 28,973 registros fenotípicos para análisis posteriores.

Para WW se eliminaron los registros con edades al destete inconsistentes, que estuvieron tres desviaciones estándar por arriba o por debajo de la media, sin información de manejo, hato o edad de la madre. Los grupos contemporáneos se definieron combinando los efectos de manejo, sexo, hato, año y época de nacimiento. El manejo para WW se definió de tres maneras 1) alimentación de la cría con leche de la madre; 2) leche de la madre más alimento balanceado; y 3) leche de la madre y de una nodriza más alimento balanceado. Las épocas de nacimiento se definieron igual que para BW, así como el criterio de eliminación dentro de grupos contemporáneos. Finalmente se obtuvieron 18,994 registros fenotípicos para análisis subsecuentes.

Los fenotipos seleccionados para el análisis con el método ssGBLUP fueron de los animales genotipados y de sus ancestros cinco generaciones atrás. La cantidad de datos con que se realizaron estos análisis se muestran en el Tabla 1.

Tabla 1 Número de animales genotipados y no genotipados utilizados en el estudio para los pesos al nacimiento (BW) y al destete (WW). 

Trait Animals Genotyped Non-genotyped
BW 325 226 99
WW 255 206 49

Genotipos

La información de los genotipos se editó de la siguiente manera: información genómica: se tomaron los SNPs en común del chip de 30K y del de 50K, obteniéndose 12,835 SNPs de los 300 animales; recodificación: los efectos aditivos se recodificaron como AA=0, AB=1 y BB=2; imputación: los genotipos faltantes se imputaron dentro de tamaño de chip utilizando las muestras de distribución marginal de los genotipos de los marcadores; y calidad del genotipo: se eliminaron los SNPs monomórficos y con diferentes umbrales de MAF. En la Tabla 2 se observan los números de marcadores SNP disponibles una vez que se aplicaron los umbrales MAF estudiados.

Tabla 2 Frecuencia de alelo menor (MAF) y número de marcadores SNP disponibles para el análisis de los pesos al nacimiento y al destete. 

MAF SNP Markers
0 12,557
0.02 11,496
0.04 8,629
0.05 6,648

Cálculo de la matriz H

La matriz de relaciones entre animales genotipados y no genotipados (H) se obtuvo de manera similar a como lo descrito por Christensen & Lund (2010) y Legarra et al. (2009). La matriz H, que incluye información de pedigrí y genómica, se define como:

H=Ann+A'gnAgg-1Ga-AggAgg-1AgnA'gnAgg-1GaGaAgg-1AgnGa

donde Ann, Agn y Agg son submatrices de A que contienen las relaciones entre animales no genotipados, entre genotipados con no genotipados y entre genotipados, respectivamente.

Ga es una matriz de relación ajustada obtenida de la matriz G descrita anteriormente:

Ga=βG+α

Ga=βG+ y se obtienen al resolver el sistema de ecuaciones:

En la primera ecuación el promedio de los elementos de la

AvgGβ+α=AvgAggAvgdiag(G)β+α=AvgdiagAgg

diagonal de G es igual al promedio de los elementos de la diagonal de Agg. En la segunda ecuación el promedio de todos los elementos de G es igual al promedio de todos los elementos de Agg. Para más detalles sobre el cálculo de H consultar (Christensen et al., 2012).

Análisis estadístico

Las características de crecimiento BW y WW se analizaron bajo el enfoque de análisis genómico en un solo paso, modificación del método de mejor predicción genómica insesgada genómica (GBLUP) propuesto por (VanRaden et al., 2009), siendo el cambio más importante la permutación de la matriz G, por la matriz H. En notación matricial el modelo puede escribirse como:

y=xβ+z1c+z2u+e

donde y es el vector de fenotipos; Z1 y Z2 son matrices de incidencia que relacionan las mediciones fenotípicas con los efectos fijos y los efectos aleatorios y ; y es el vector de efectos aleatorios de error residual. Las varianzas aditivas genéticas, de grupo contemporáneo y residual se supusieron:

Varu=Hσa2Varc=Iσg2 yVare=Iσe2

respectivamente.

El modelo para BW (1) incluyó los efectos fijos de la covariable edad de la madre lineal (EM) y cuadrático (EM2); el efecto de grupo contemporáneo (GC) fue considerado aleatorio. El modelo para WW (2) incluyó los mismos efectos fijos que (1) más el de la covariable edad al destete (ED).

BW=EM+EM2+GC+e (1)

WW=ED+EM+EM2+GC+e (2)

Validación cruzada

La habilidad predictiva de los modelos utilizados se probó mediante validación cruzada. La población se dividió en cinco grupos, uno de los cuales sirvió como de entrenamiento y el resto como de prueba. La validación cruzada se repitió cuatro ocasiones, con cada uno de los grupos de datos de prueba y finalmente se obtuvo la media aritmética de los resultados de cada iteración obteniendo un único resultado para cada análisis.

Criterios de comparación de los análisis Efecto de MAF en la jerarquización de animales

Para determinar el efecto de MAF sobre el grado de similitud entre la jerarquización de animales de acuerdo con los valores genómicos predichos (GEBVs) para BW y WW, se obtuvieron los estimadores de correlación Spearman y Pearson entre los GEBVs obtenidos con MAF = 0.05 y con los correspondientes GEBVs obtenidos con MAF = 0, 0.02 y 0.04.

Habilidad de predicción de los modelos

Para evaluar el efecto de MAF en la habilidad predictiva de los modelos se realizó una validación cruzada. Se dividió la información en conjuntos de entrenamiento y validación, utilizando el conjunto de validación para evaluar la capacidad predictiva del modelo entrenado. Se eliminó aleatoriamente 20 % de los valores fenotípicos de las características de crecimiento BW y WW. Posteriormente para cada modelo con MAF de 0, 0.02, 0.04 y 0.05 se obtuvieron los GEBVs para dichas características y se obtuvo el promedio de cuatro repeticiones de los coeficientes de correlación de los fenotipos eliminados con los predichos.

Coeficiente de regresión lineal simple

Los estimadores de los coeficientes de regresión lineal simple se utilizaron para determinar el efecto de MAF en los cambios en magnitud de los GEBVs, utilizando como variables independientes los GEBVs con MAF 0.05 y como variables dependientes los valores genómicos predichos con MAF 0, 0.02 y 0.04. La hipótesis nula probada fue que el estimador del coeficiente de regresión fue igual que cero.

Para todos los análisis estadísticos se utilizó el programa R (RStudio Team, 2019).

Resultados y discusión

Jerarquización de los animales

Los coeficientes de correlación de Spearman y Pearson entre los GEBVs con los análisis considerando diferentes umbrales de MAF se muestran en la Tabla 3. Los estimadores en ambos tipos de correlaciones fueron mayores que 0.999 (p ≤ 2x10-16), hubo pocos cambios en la jerarquización de los animales de acuerdo con los GEBVs obtenidos con los diferentes umbrales probados. Esto indica que el efecto de MAF no fue importante en la jerarquización de los GEBVs.

Tabla 3 Estimadores de los coeficientes de correlación Pearson y Spearman, en paréntesis nivel de significancia para la prueba de hipótesis del estimador igual que cero, entre valores genómicos predichos con frecuencia del alelo menor (MAF) = 0.05 contra MAF = 0, 0.02 y 0.04 para pesos al nacimiento (BW) y al destete (WW) en bovinos Suizo Europeo. 

MAF BW WW
Pearson Spearman Pearson Spearman
0 0.99993
(<2×10-16)
0.99980
(<2×10-16)
0.99976
(<2.2×10-16)
0.99960
(<2.2×10-16)
0.02 0.99990
(<2×10-16)
0.99976
(<2×10-16)
0.99974
(<2.2×10-16)
0.99957
(<2.2×10-16)
0.04 0.99991
(<2×10-16)
0.99978
(<2×10-16)
0.99983
(<2.2×10-16)
0.99968
(<2.2×10-16)

MAF = umbral de frecuencia de alelo menor; BW = peso al nacimiento; WW = peso al destete.

Estos resultados sugieren que con fines de selección de animales como progenitores, cuando las evaluaciones se hacen utilizando información de SNPs, se puede utilizar cualquier umbral MAF entre 0.0 y 0.05 en el control de calidad durante la edición de genotipos, aunque 0.05 es el estándar en las evaluaciones genómicas realizadas con animales de acuerdo con (VanRaden et al., 2009).

Las precisiones obtenidas mediante los estimadores de correlación fueron ligeramente mayores para BW que para WW, esto puede deberse al mayor tiempo de influencia del ambiente sobre los animales para WW y que puede afectar las predicciones.

En el control de calidad del genotipo en bovinos HolsteinFriesian australianos, usando el chip Illumina Bovine SNP50TM y 798 animales, se eliminó MAF <2.5 % (Hayes et al., 2009). Estos autores obtuvieron para peso vivo exactitudes de los GEBVs de 50 a 67 %, similares a las obtenidas en el presente estudio. Por su parte, Wiggans et al. (2009) eliminaron registros utilizando MAF menores que 2 % en genotipos de 941 bovinos Jersey y 344 toros Suizo Europeo, las dos poblaciones de animales utilizadas en estos estudios no pueden considerarse grandes de acuerdo con VanRaden et al. (2009). Estos autores sugieren usar MAF no menor que 5 %, aunque en poblaciones grandes se puede disminuir sin afectar la precisión.

Wiggans et al. (2009) argumentan que el número de animales disponibles por raza en su estudio hizo que sus resultados se consideraran preliminares, ya que el tamaño de muestra no permitió obtener resultados definitivos. Lo anterior puede explicar los resultados obtenidos en la presente investigación, pues la cantidad de animales utilizados no se considera una población grande, en comparación con los tamaños poblacionales utilizados en otros estudios (VanRaden et al., 2009).

Habilidad de predicción del modelo con diferente umbral de MAF

La Tabla 4 muestra que establecer un MAF mayor en el control de calidad en la edición de los genotipos, mejora la precisión de predicción de los GEBVs, aunque la diferencia en precisión, para BW, entre utilizar MAF = 0.0 contra MAF = 0.05 es muy baja, 0.0027, y es más baja para WW donde la diferencia fue 0.0007. Esto significa que, en el GQC de una muestra de 300 animales con 12,835 SNPs, el umbral de MAF a utilizar de hasta 0.05 no es importante en la predicción de GEBVs para BW y WW.

Tabla 4 Promedios de los estimadores de correlación Pearson entre los valores genómicos predichos mediante análisis con diferentes umbrales de frecuencia del alelo menor (MAF) y los valores fenotípicos reales para los pesos al nacimiento (BW) y al destete de cuatro repeticiones de validación cruzada (WW). 

MAF BW WW
0 0.7130 0.6444
0.02 0.7127 0.6447
0.04 0.7136 0.6451
0.05 0.7136 0.6451

MAF = umbral de frecuencia de alelo menor; BW = peso al nacimiento; WW = peso al destete.

Linck & Battey (2019) realizaron un estudio con información de SNPs simulada, utilizando umbrales de MAF entre 0.017 y 0.25; estos autores reportaron que utilizar un MAF mayor que 0.03 implica un cambio en la estructura de la distribución de los datos por la caída en el tamaño total de la matriz después del filtrado y genera un aumento en la precisión de los predictores. Lo anterior puede explicar que en el presente estudio utilizar MAF = 0.05 para BW, implique obtener ligeramente mejores predicciones.

La habilidad de predicción de los modelos para BW y WW para cada uno de los umbrales estudiados no representa cambios significativos. Lo anterior puede explicarse por la poca cantidad de SNPs eliminados dada la información utilizada comparada con la información que se usó en otros estudios como el realizado por (Zhu et al., 2017) con 82,594 marcadores por (De la Cruz & Raska, 2014), 40,000 en el estudio de (Linck & Battey, 2019), y 600,000 en el estudio con aves de (Abdollahi-Arpanahi et al., 2014). Estos autores reportaron un efecto importante del MAF en la predicción genómica de las características estudiadas. Un aspecto adicional es el bajo número de animales que se utilizó en el presente estudio, lo que no permite obtener resultados concluyentes. Se recomienda realizar el estudio con un número mayor de animales.

Estimador de los coeficientes de regresión

En la Tabla 5 se observa que los coeficientes de regresión son similares entre sí; sin embargo, el coeficiente para WW se reduce ligeramente conforme incrementa el MAF para el control de calidad en la edición del genotipo. Conforme los estimadores de los coeficientes de regresión sean diferentes de la unidad, la diferencia en la predicción de los GEBVs con las alternativas de análisis en comparación es mayor (Mäntysaari et al., 2010), todos los estimadores de los coeficientes de regresión se encuentran cerca de la unidad en consecuencia las diferencias entre los efectos de los umbrales de MAF probados son mínimas.

Tabla 5 Estimadores de los coeficientes de regresión (± SE), las variables independientes fueron los valores genómicos predichos con frecuencia del alelo menor (MAF) = 0.05 y como las dependientes los valores genómicos predichos con MAF de 0, 0.02 y 0.04; entre paréntesis nivel de significancia para la prueba de hipótesis del estimador igual que cero, para los pesos al nacimiento (BW) y al destete (WW) en bovinos Suizo Europeo. 

MAF BW WW
0 0.99756±0.00072
(<2×10-16)
0.99323 ± 0.00136
(<2×10-16)
0.02 0.99605±0.00076
(<2×10-16)
0.99255 ± 0.00141
(<2×10-16)
0.04 0.99578 ± 0.00072
(<2×10-16)
0.99400±0.00114
(<2×10-16)

MAF = umbral de frecuencia de alelo menor; BW = peso al nacimiento; WW = peso al destete.

Los resultados del presente estudio son similares a los reportados por (Edriss et al., 2012), quienes utilizando información de 5,500 animales Holstein y Jersey con 54,000 SNPs, estimaron para las características de fertilidad, proteína en leche y mastitis, las correlaciones entre valores reales y predichos usando los umbrales de MAF: 0.075, 0.275, 0.3, 0.325 y 0.35, sus resultados mostraron mínimas diferencias, por ejemplo, en Holstein las correlaciones para fertilidad tuvieron una diferencia máxima de 0.014; concluyendo que los resultados de su estudio sugieren que para una población de referencia pequeña, los marcadores SNP con MAF bajo no dañan la predicción genómica, afirmación que también aplica a lo encontrado en el presente estudio, pero para características de crecimiento.

De acuerdo con lo reportado por VanRaden et al. (2009) y Wiggans et al. (2009), en relación con el número de animales que conforman la población y el número de marcadores, los resultados utilizados en este estudio deberán tomarse con precaución, y como mencionaron Hayes et al. (2009), deberán considerarse preliminares y susceptibles a cambios que en un futuro pudieran presentarse al repetir el estudio con un número mayor de animales y marcadores, ya que el efecto de MAF se debe a la caída en el tamaño total de la matriz después del filtrado en la GQC, caída que puede no ser significativa si la cantidad de información es limitada (Linck & Battey, 2019).

Conclusión

El umbral de frecuencia de alelo menor a utilizar en la evaluación genómica puede encontrarse en el rango de 0 a 0.05, sin alterar de forma importante la predicción de valores genómicos, pues no tiene un predictive efecto significativo sobre la jerarquización de los valores genómicos predichos y la habilidad predictiva del modelo. Se reconoce que 0.05 puede ser un umbral recomendable para utilizar en evaluación de poblaciones pequeñas y con un número de marcadores limitado.

Bibliografía

Abdollahi-Arpanahi, R., Nejati-Javaremi, A., Pakdel, A., Moradi-Shahrbabak, M., Morota, G., Valente, B. D., Kranis, A., Rosa, G. J. M. and Gianola, D. (2014). Effect of allele frequencies, effect sizes and number of markers on prediction of quantitative traits in chickens. Journal of Animal Breeding and Genetics, 131(2): 123-133. https://doi.org/10.1111/jbg.12075 [ Links ]

Anderson, C. A., Pettersson, F. H., Clarke, G. M., Cardon, L. R., Morris, A. P. and Zondervan, K. T. (2010). Data quality control in genetic case-control association studies. Nature Protocols, 5, 1564-1573. https://doi.org/10.1038/nprot.2010.116 [ Links ]

Christensen, O. F. & Lund, M. S. (2010). Genomic relationship matrix when some animals are not genotyped. Genetics Selection Evolution, 42(2): 1-8. https://doi.org/10.1186/1297-9686-42-2 [ Links ]

Christensen, O. F., Madsen, P., Nielsen, B., Ostersen, T. and Su, G. (2012). Single-step methods for genomic evaluation in pigs. Animal 6(10): 1565-1571. https://doi.org/10.1017/S1751731112000742 [ Links ]

Chouraki, V. & Seshadri, S. (2014). Chapter five - Genetics of Alzheimer’s disease. Advances in Genetics, 87, 246-278. https://doi.org/10.1016/B978-0-12-800149-3.00005-6 [ Links ]

De la Cruz, O. & Raska, P. (2014). Population structure at different minor allele frequency levels. BMC Proc, 8, S55. https://doi.org/10.1186/1753-6561-8-s1-s55 [ Links ]

De los Campos, G., Hickey, J. M., Pong-Wong, R., Daetwyler, H. D. and Calus, M. P. L. (2013). Whole-genome regression and prediction methods applied to plant and animal breeding. Genetics, 193(2): 327-345. https://doi.org/10.1534/genetics.112.143313 [ Links ]

Edriss, V., Guldbrandtsen, B., Lund, M. S. and Su, G. (2012). Effect of marker-data editing on the accuracy of genomic prediction. Journal of Animal Breeding and Genetics, 130(2): 128-135. https://doi.org/10.1111/j.14390388.2012.01015.x [ Links ]

Gianola, D., de los Campos, G., Hill, W. G., Manfredi, E. and Fernando, R. (2009). Additive Genetic Variability and the Bayesian Alphabet. Genetics, 183(1): 347-363. https://doi.org/10.1534/genetics.109.103952 [ Links ]

Gianola, D., Fernando, R. L. and Stella, A. (2006). Genomic-Assisted Prediction of Genetic Value with Semiparametric Procedures. Genetics, 173(3): 1761-1776. https://doi.org/10.1534/genetics.105.049510 [ Links ]

Hayes, B. J., Bowman, P. J., Chamberlain, A. J. and Goddard, M. E. (2009). Invited review: Genomic selection in dairy cattle: Progress and challenges. Journal of Dairy Science, 92(2): 433-443. https://doi.org/10.3168/jds.2008-1646 [ Links ]

Kido, T., Sikora-Wohlfeld, W., Kawashima, M., Kikuchi, S., Kamatani, N., Patwardhan, A., Chen, R., Sirota, M., Kodama, K., Hadley, D. and Butte, A. J. (2018). Are minor alleles more likely to be risk alleles? BMC Medical Genomics, 11, 3. https://doi.org/10.1186/s12920-018-0322-5 [ Links ]

Legarra, A., Aguilar, I. and Misztal, I. (2009). A relationship matrix including full pedigree and genomic information. Journal of Dairy Science, 92(9): 4656-4663. https://doi.org/10.3168/jds.2009-2061 [ Links ]

Linck, E. & Battey, C. J. (2019). Minor allele frequency thresholds strongly affect population structure inference with genomic data sets. Molecular Ecology Resources, 19(3): 639-647. https://doi.org/10.1111/1755-0998.12995 [ Links ]

Lopez-Cruz, M., Crossa, J., Bonnett, D., Dreisigacker, S., Poland, J., Jannink, J.-L., Singh, R. P., Autrique, E. and De los Campos, G. (2015). Increased prediction accuracy in wheat breeding trials using a marker × environment interaction genomic selection model. Genes|Genomes|Genetics, 5(4): 569-582. https://doi.org/10.1534/g3.114.016097 [ Links ]

Mäntysaari, E., Liu, Z. and VanRaden, P. (2010). Interbull validation test for genomic evaluations. INTERBULL BULLETIN No 41. https://journal.interbull.org/index.php/ib/article/view/1134Links ]

Meuwissen, T. H. E., Hayes, B. J. and Goddard, M. E. (2001). Prediction of total genetic value using genome-wide dense marker maps. Genetics, 157(4), 1819-1829. https://www.genetics.org/content/157/4/1819#:~:text=A%20genome%20of%201000%20cM,marker%20spacing%20of%201%20cM.&text=It%20was%20concluded%20that%20selection,to%20shorten%20the%20generation%20intervalLinks ]

Meuwissen, T., Hayes, B. and Goddard, M. (2016). Genomic selection: A paradigm shift in animal breeding. Animal Frontiers, 6(1): 6-14. https://doi.org/10.2527/af.2016-0002 [ Links ]

Minozzi, G., Williams, J. L., Stella, A., Strozzi, F., Luini, M., Settles, M. L., Taylor, J. F., Whitlock, R. H., Zanella, R. and Neibergs, H. L. (2012). Meta-analysis of two genome-wide association studies of bovine paratuberculosis. PLOS ONE, 7(3): e32578. https://doi.org/10.1371/journal.pone.0032578 [ Links ]

Pértile, S. F. N., Silva, F. F., Salvian, M. and Mourão, G. B. (2016). Seleção e associação genômica ampla para o melhoramento genético animal com uso do método ssGBLUP. Pesquisa Agropecuária Brasileira, 51(10): 1729-1736. https://doi.org/10.1590/s0100-204x2016001000004 [ Links ]

RStudio Team. (2019). RStudio: Integrated Development for R. RStudio, Inc., Boston, MA. https://www.rstudio.com/Links ]

Saatchi, M., McClure, M. C., McKay, S. D., Rolf, M. M., Kim, J., Decker, J. E., Taxis, T. M., Chapple, R. H., Ramey, H. R., Northcutt, S. L., Bauck, S., Woodward, B., Dekkers, J. C. M., Fernando, R. L., Schnabel, R. D., Garrick, D. J. and Taylor, J. F. (2011). Accuracies of genomic breeding values in American Angus beef cattle using K-means clustering for cross-validation. Genetics Selection Evolution, 43, 40. https://doi.org/10.1186/1297-9686-43-40 [ Links ]

VanRaden, P. M., Van Tassell, C. P., Wiggans, G. R., Sonstegard, T. S., Schnabel, R. D., Taylor, J. F. and Schenkel, F. S. (2009). Invited Review: Reliability of genomic predictions for North American Holstein bulls. Journal of Dairy Science, 92(1): 16-24. https://doi.org/10.3168/jds.2008-1514 [ Links ]

Wang, H., Misztal, I., Aguilar, I., Legarra, A., Fernando, R. L., Vitezica, Z., Okimoto, R., Wing, T., Hawken, R. and Muir, W. M. (2014). Genome-wide association mapping including phenotypes from relatives without genotypes in a singlestep (ssGWAS) for 6-week body weight in broiler chickens. Frontiers in Genetics, 5, 134. https://doi.org/10.3389/fgene.2014.00134 [ Links ]

Wiggans, G. R., Sonstegard, T. S., VanRaden, P. M., Matukumalli, L. K., Schnabel, R. D., Taylor, J. F., Schenkel, F. S. and Van Tassell, C. P. (2009). Selection of single-nucleotide polymorphisms and quality of genotypes used in genomic evaluation of dairy cattle in the United States and Canada. Journal of Dairy Science, 92, 3431-3436. https://doi.org/10.3168/jds.2008-1758 [ Links ]

Zhu, B., Zhang, J., Niu, H., Guan, L., Guo, P., Xu, L., Chen, Y., Zhang, L., Gao, H., Gao, X. and Li, J. (2017). Effects of marker density and minor allele frequency on genomic prediction for growth traits in Chinese Simmental beef cattle. Journal of Integrative Agriculture, 16(4): 911-920. https://doi.org/10.1016/s2095-3119(16)61474-0 [ Links ]

Como citar este artículo: Trujano-Chavez, M. Z., Valerio-Hernández, J. E., López-Ordaz, R., Ruíz-Flores, A. (2021). Minor allele frequency in genomic prediction for growth traits in Braunvieh cattle. Revista Bio Ciencias 8, e1052. doi: DOI: https://doi.org/10.15741/revbio.08.e1052

Recibido: 29 de Agosto de 2020; Aprobado: 22 de Marzo de 2021

* Corresponding Author: Agustín Ruíz-Flores. Universidad Autónoma Chapingo, Carretera Federal México-Texcoco Km 38.5, C.P. 56227, Texcoco, Estado de México, México. Phone: (595) 952 1683. E-mail: aruizf@chapingo.mx, https://produccionanimal.chapingo.mx/

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons