Exploring the molecular basis of dairy production and viability traits in Murciano-Granadina goats with genomic, metagenomic and transcriptomic approaches

Tesis presentada el día 19 de diciembre de 2022 para la obtención del título de doctor en Producción Animal por la Universidad Autónoma de Barcelona.

[EN] In this thesis, we aimed to investigate the genetic basis of milk and morphological traits in Murciano-Granadina goats as well as to explore the microbial composition of several types of caprine cheeses. Moreover, we wanted to assess the consequences of inbreeding on dairy performance and to study the molecular basis of traits related with the biological viability of goats. In study 1, we have identified several genetic determinants of milk yield and composition traits in Murciano-Granadina goats with records for three lactations by performing three independent GWAS (one for each lactation) and a longitudinal GWAS. Both approaches consistently revealed a genome-wide significant QTL for protein percentage on a chromosome 6 (74. 8-94. 6 Mb) region which harbours the four casein genes. Another QTL on chromosome 2 (129. 77-131. 01 Mb) for lactose percentage was also consistently detected. More QTL were identified in the longitudinal GWAS than in the three independent GWAS, possibly because the former has more statistical power than the latter. In study 2, we also used a GWAS approach to identify genomic regions associated with the variation of morphological traits in Murciano-Granadina goats. In this analysis, few significant associations were found and positional coincidence with morphology QTL detected in other breeds was low. These results are compatible with a highly polygenic determinism of morphological traits in goats. Since dairy production in goats is focused to the cheesemaking and the technological and organoleptic properties of this food are largely determined by its microbiota, in study 3, we characterised the microbiota of six different goat cheeses by sequencing hyper-variable regions of the bacterial 16S rRNA gene. We have identified lactic acid bacteria as the predominant microbial community in all six cheeses. Besides, we detected an abundance of psychrophilic bacteria, which are common post-pasteurisation contaminants of milk, in fresh cheeses. Regarding viability traits, in study 4, we detected low levels of inbreeding in the Murciano-Granadina population under study (FROH = 0. 053±0. 04). Despite this, a significant inbreeding depression on the logarithm of the somatic cell score (lnSCC), an important trait for udder health and milk quality, was detected and its determinants were finely mapped to a few chromosome 8 and 25 regions containing immunity genes. To ascertain genomic regions harbouring potentially harmful mutations (that are not expected to show a Mendelian segregation in the offspring of carrier individuals), in study 5 we carried out a scan to detected transmission ratio distortion (TRD) in Murciano-Granadina goats at a genome-wide scale.

[ES] En esta tesis, hemos querido investigar las bases genéticas de los caracteres lecheros y morfológicos de cabras Murciano-Granadina, así como explorar la composición microbiana de diversos tipos de quesos caprinos. Adicionalmente, hemos querido evaluar las consecuencias de la consanguinidad sobre la producción y composición de la leche de cabra y estudiar la base molecular de distintos caracteres relacionados con la viabilidad biológica. En el estudio 1, hemos identificado varios marcadores genéticos asociados a la producción y composición de la leche en la raza Murciano-Granadina utilizando registros de tres lactaciones y aplicando tres GWAS independientes (por lactación) y un GWAS longitudinal. Ambas aproximaciones revelaron un QTL significativo para el porcentaje de proteínas en una región del cromosoma 6 (74,8-94,6 Mb) que alberga los cuatro genes de las caseínas. Se detectó un QTL en el cromosoma 2 (129,77-131,01 Mb) para el porcentaje de la lactosa. Se identificaron más QTL en el GWAS longitudinal que en los tres GWAS independientes, posiblemente porque el GWAS longitudinal tiene mayor poder estadístico. En el estudio 2 también aplicamos GWAS para identificar regiones genómicas asociadas a la variación de caracteres morfológicos en cabras de la raza Murciano-Granadina. Se hallaron pocas asociaciones significativas y la coincidencia posicional con QTL morfológicos detectados en otras razas fue baja. Estos resultados sugieren un determinismo poligénico de los caracteres morfológicos en caprino. Dado que la producción de leche en cabra está enfocada a la elaboración de quesos y las propiedades tecnológicas y organolépticas de este alimento están en gran parte determinadas por su microbiota, en el estudio 3, hemos caracterizado la microbiota de seis quesos de cabra secuenciando regiones hipervariables del gen bacteriano 16S ARNr. Las bacterias acido-lácticas son las más predominantes en los seis quesos. En los quesos frescos hemos detectado cierta abundancia de bacterias psicrófilas, contaminantes habituales de la leche post-pasteurización. Respecto a los rasgos de viabilidad, en el estudio 4 se detectaron niveles bajos de consanguinidad en la población Murciano-Granadina estudiada (FROH = 0. 053±0. 04). Sin embargo, se detectó depresión consanguínea significativa para el logaritmo del recuento de células somáticas (lnSCC), importante para la salud de la glándula mamaria y la calidad de la leche, y se mapeó finamente los determinantes de esta depresión en regiones de los cromosomas 8 y 25 que contienen genes inmunitarios. Para determinar las regiones genómicas que contienen mutaciones potencialmente adversas sobre la viabilidad en el estudio 5 investigamos la existencia de distorsiones del ratio de transmisión alélica (TRD) en cabras Murciano-Granadinas.

[CA] En aquesta tesi, hem volgut investigar les bases genètiques dels caràcters lleters i morfològics de les cabres Murciano-Granadina, així com explorar la composició microbiana de diversos tipus de formatges caprins. Addicionalment, hem volgut avaluar les conseqüències de la consanguinitat sobre la producció i composició de la llet de cabra i estudiar la base molecular de diversos caràcters relacionats amb la viabilitat biològica. En l'estudi 1, hem identificat diversos marcadors genètics associats a la producció i composició de la llet en la raça Murciano-Granadina emprant registres de tres lactacions i aplicant tres GWAS independents (per lactació) i un GWAS longitudinal. Ambdues aproximacions van revelar un QTL significatiu per al percentatge de proteïnes en una regió del cromosoma 6 (74,8-94,6 Mb) que alberga els quatre gens de les caseïnes. Es va detectar un QTL al cromosoma 2 (129,77-131,01 Mb) per al percentatge de la lactosa. Es van identificar més QTL al GWAS longitudinal que als tres GWAS independents, possiblement perquè el GWAS longitudinal té més poder estadístic. A l'estudi 2, també hem aplicat GWAS per identificar regions genòmiques associades a la variació de caràcters morfològics en cabres de la raça Murciano-Granadina. Es van trobar poques associacions significatives i la coincidència posicional amb QTL morfològics detectats en altres races va ser baixa. Aquests resultats suggereixen un determinisme poligènic dels caràcters morfològics en cabrum. Atès que la producció de llet en cabra està enfocada a l'elaboració de formatges i les propietats tecnològiques i organolèptiques d'aquest aliment estan en gran part determinades per la seva microbiota, a l'estudi 3, hem caracteritzat la microbiota de sis formatges de cabra mitjançant la seqüenciació de regions hipervariables del gen bacterià 16S ARNr. Els bacteris acido-làctics són els més predominants en els sis formatges. En els formatges frescos hem detectat certa abundància de bacteris psicròfils, contaminants habituals de la llet després de la pasteurització. Respecte als trets de viabilitat, a l'estudi 4 s'han detectat nivells baixos d'endogàmia a la població Murciano-Granadina estudiada (FROH = 0. 053±0. 04). Malgrat això, es va detectar una depressió endogàmica significativa pel logaritme del recompte de cèl·lules somàtiques (lnSCC), important per a la salut de la glàndula mamària i la qualitat de la llet, i es va mapejar finament els determinants d'aquesta depressió a regions dels cromosomes 8 i 25 que contenen gens immunitaris. Per determinar les regions genòmiques que contenen mutacions potencialment adverses sobre la viabilitat a l'estudi 5 vam investigar l'existència de distorsions de la ratio de transmissió aŀlèlica (TRD) en cabres Murciano-Granadines.