Analyzing the genomic and transcriptomic architecture of milk traits in Murciano-Granadina goats

Dailu Guan, Vincenzo Landi, María Gracia Luigi-Sierra, Juan Vicente Delgado, Xavier Such, Anna Castelló, Betlem Cabrera, Emilio Mármol-Sánchez, Javier Fernández-Alvarez, José Luis Ruiz De La Torre Casañas, Amparo Martínez, Jordi Jordana, Marcel Amills*

*Corresponding author for this work

Research output: Contribution to journalArticleResearchpeer-review

10 Citations (Scopus)


Background: In this study, we aimed to investigate the molecular basis of lactation as well as to identify the genetic factors that influence milk yield and composition in goats. To achieve these two goals, we have analyzed how the mRNA profile of the mammary gland changes in seven Murciano-Granadina goats at each of three different time points, i.e. 78 d (T1, early lactation), 216 d (T2, late lactation) and 285 d (T3, dry period) after parturition. Moreover, we have performed a genome-wide association study (GWAS) for seven dairy traits recorded in the 1st lactation of 822 Murciano-Granadina goats. Results: The expression profiles of the mammary gland in the early (T1) and late (T2) lactation were quite similar (42 differentially expressed genes), while strong transcriptomic differences (more than one thousand differentially expressed genes) were observed between the lactating (T1/T2) and non-lactating (T3) mammary glands. A large number of differentially expressed genes were involved in pathways related with the biosynthesis of amino acids, cholesterol, triglycerides and steroids as well as with glycerophospholipid metabolism, adipocytokine signaling, lipid binding, regulation of ion transmembrane transport, calcium ion binding, metalloendopeptidase activity and complement and coagulation cascades. With regard to the second goal of the study, the performance of the GWAS allowed us to detect 24 quantitative trait loci (QTLs), including three genome-wide significant associations: QTL1 (chromosome 2, 130.72-131.01 Mb) for lactose percentage, QTL6 (chromosome 6, 78.90-93.48 Mb) for protein percentage and QTL17 (chromosome 17, 11.20 Mb) for both protein and dry matter percentages. Interestingly, QTL6 shows positional coincidence with the casein genes, which encode 80% of milk proteins. Conclusions: The abrogation of lactation involves dramatic changes in the expression of genes participating in a broad array of physiological processes such as protein, lipid and carbohydrate metabolism, calcium homeostasis, cell death and tissue remodeling, as well as immunity. We also conclude that genetic variation at the casein genes has a major impact on the milk protein content of Murciano-Granadina goats.

Original languageAmerican English
Article number35
JournalJournal of Animal Science and Biotechnology
Issue number1
Publication statusPublished - 11 Mar 2020


  • Casein genes
  • Dairy traits
  • GWAS
  • Lactation
  • QTLs
  • RNA-Seq


Dive into the research topics of 'Analyzing the genomic and transcriptomic architecture of milk traits in Murciano-Granadina goats'. Together they form a unique fingerprint.

Cite this