Background: Assessing genetic biodiversity and population structure of minor breeds through the information provided by neutral molecular markers, allows determination of their extinction risk and to design strategies for their management and conservation. Analysis of microsatellite loci is known to be highly informative in the reconstruction of the historical processes underlying the evolution and differentiation of animal populations. Guadarrama goat is a threatened Spanish breed which actual census (2008) consists of 3057 females and 203 males distributed in 22 populations more or less isolated. The aim of this work is to study the genetic status of this breed through the analysis of molecular data from 10 microsatellites typed in historic and actual live animals. Results: The mean expected heterozygosity across loci within populations ranged from 0.62 to 0.77. Genetic differentiation measures were moderate, with a mean FST of 0.074, GST of 0.081 and RST of 0.085. Percentages of variation among and within populations were 7.5 and 92.5, respectively. Bayesian clustering analyses pointed out a population subdivision in 16 clusters, however, no correlation between geographical distances and genetic differences was found. Management factors such as the limited exchange of animals between farmers (estimated gene flow Nm = 3.08) mostly due to sanitary and social constraints could be the major causes affecting Guadarrama goat population subdivision. Conclusion: Genetic diversity measures revealed a good status of biodiversity in the Guadarrama goat breed. Since diseases are the first cause affecting the census in this breed, population subdivision would be an advantage for its conservation. However, to maintain private alleles present at low frequencies in such small populations minimizing the inbreeding rate, it would necessitate some mating designs of animals carrying such alleles among populations. The systematic use of molecular markers will facilitate the comprehensive management of these populations, which in combination with the actual breeding program to increase milk yield, will constitute a good strategy to preserve the breed.

Mitigating the effects of global warming has become the main challenge for humanity in recent decades. Livestock farming contributes to greenhouse gas emissions, with an important output of methane from enteric fermentation processes, mostly in ruminants. Because ruminal microbiota is directly involved in digestive fermentation processes and methane biosynthesis, understanding the ecological relationships between rumen microorganisms and their active metabolic pathways is essential for reducing emissions. This study analysed whole rumen metagenome using long reads and considering its compositional nature in order to disentangle the role of rumen microbes in methane emissions. The β-diversity analyses suggested a subtle association between methane production and overall microbiota composition (0.01 < R² < 0.02). Differential abundance analysis identified 36 genera and 279 KEGGs as significantly associated with methane production (P_adj < 0.05). Those genera associated with high methane production were Eukaryota from Alveolata and Fungi clades, while Bacteria were associated with low methane emissions. The genus-level association network showed 2 clusters grouping Eukaryota and Bacteria, respectively. Regarding microbial gene functions, 41 KEGGs were found to be differentially abundant between low- and high-emission animals and were mainly involved in metabolic pathways. No KEGGs included in the methane metabolism pathway (ko00680) were detected as associated with high methane emissions. The KEGG network showed 3 clusters grouping KEGGs associated with high emissions, low emissions, and not differentially abundant in either. A deeper analysis of the differentially abundant KEGGs revealed that genes related with anaerobic respiration through nitrate degradation were more abundant in low-emission animals. Methane emissions are largely associated with the relative abundance of ciliates and fungi. The role of nitrate electron acceptors can be particularly important because this respiration mechanism directly competes with methanogenesis. Whole metagenome sequencing is necessary to jointly consider the relative abundance of Bacteria, Archaea, and Eukaryota in the statistical analyses. Nutritional and genetic strategies to reduce CH₄ emissions should focus on reducing the relative abundance of Alveolata and Fungi in the rumen. This experiment has generated the largest ONT ruminal metagenomic dataset currently available.

To date, there is a lack of research into the vaginal and sperm microbiome and its bearing on artificial insemination (AI) success in the ovine species. Using hypervariable regions V3–V4 of the 16S rRNA, we describe, for the first time, the combined effect of the ovine microbiome of both females (50 ewes belonging to five herds) and males (five AI rams from an AI center) on AI outcome. Differences in microbiota abundance between pregnant and non-pregnant ewes and between ewes carrying progesterone-releasing intravaginal devices (PRID) with or without antibiotic were tested at different taxonomic levels. The antibiotic treatment applied with the PRID only altered Streptobacillus genus abundance, which was significantly lower in ewes carrying PRID with antibiotic. Mageebacillus, Histophilus, Actinobacilllus and Sneathia genera were significantly less abundant in pregnant ewes. In addition, these genera were more abundant in two farms with higher AI failure. Species of these genera such as Actinobacillus seminis and Histophilus somni have been associated with reproductive disorders in the ovine species. These genera were not present in the sperm samples of AI rams, but were found in the foreskin samples of rams belonging to herd 2 (with high AI failure rate) indicating that their presence in ewes’ vagina could be due to prior transmission by natural mating with rams reared in the herd.

Molecular chaperones have long been understood to be preferentially transcribed in response to multiple perturbations of the cellular homeostasis. In this study, several polymorphisms in the gene encoding the inducible form of the cytoplasmic Hsp90 (HSP90AA1) were addressed in 24 sheep breeds reared in different climatic regions of Europe, Africa, and Asia. Significant differences in the genotype frequencies for a C/G single nucleotide polymorphism (SNP) located at position −660 in the HSP90AA1 5′flanking region were found between the different breeds. Regression analyses reflected significant correlations (from 0.41 to 0.62) between the alternative genotypes of this polymorphism and several climatic and geographic variables characteristic of the regions where these breeds are reared. Real-time analysis revealed that animals bearing the CC₋₆₆₀ genotype presented higher expression levels than those presenting the CG₋₆₆₀ or GG₋₆₆₀ in summer, but not in spring. Mutation at −660 site seems to affect HSP90AA1 transcription rates which could have important effects on the adaptation to different environmental conditions in sheep. Thus, the variability found in the genotype frequencies for the SNP at −660 in the ovine HSP90AA1 locus could be the result of the different environmental pressures occurring in the regions where these breed are maintained.

While physiological differences across skeletal muscles have been described, the differential gene expression underlying them and the discovery of how they interact to perform specific biological processes are largely to be elucidated. The purpose of the present study was, firstly, to profile by cDNA microarrays the differential gene expression between two skeletal muscle types, Psoas major (PM) and Flexor digitorum (FD), in beef cattle and then to interpret the results in the context of a bovine gene coexpression network, detecting possible changes in connectivity across the skeletal muscle system. Eighty four genes were differentially expressed (DE) between muscles. Approximately 54% encoded metabolic enzymes and structural-contractile proteins. DE genes were involved in similar processes and functions, but the proportion of genes in each category varied within each muscle. A correlation matrix was obtained for 61 out of the 84 DE genes from a gene coexpression network. Different groups of coexpression were observed, the largest one having 28 metabolic and contractile genes, up-regulated in PM, and mainly encoding fast-glycolytic fibre structural components and glycolytic enzymes. In FD, genes related to cell support seemed to constitute its identity feature and did not positively correlate to the rest of DE genes in FD. Moreover, changes in connectivity for some DE genes were observed in the different gene ontologies. Our results confirm the existence of a muscle dependent transcription and coexpression pattern and suggest the necessity of integrating different muscle types to perform comprehensive networks for the transcriptional landscape of bovine skeletal muscle.

The aim of this work was to investigate the association between polymorphisms located at the HSP90AA1 ovine gene promoter and gene expression rate under different environmental conditions, using a mixed model approach. Blood samples from 120 unrelated rams of the Manchega sheep breed were collected at three time points differing in environmental conditions. Rams were selected on the basis of their genotype for the transversion G/C located 660 base pairs upstream the gene transcription initiation site. Animals were also genotyped for another set of 6 SNPs located at the gene promoter. Two SNPs, G/C₋₆₆₀ and A/G₋₄₄₄, were associated with gene overexpression resulting from heat stress. The composed genotype CC₋₆₆₀-AG₋₄₄₄ was the genotype having the highest expression rates with fold changes ranging from 2.2 to 3.0. The genotype AG₋₅₂₂ showed the highest expression levels under control conditions with a fold change of 1.4. Under these conditions, the composed genotype CC₋₆₀₁-TT₋₅₂₄-AG₋₅₂₂-TT₋₄₆₈ is expected to be correlated with higher basal expression of the gene according to genotype frequencies and linkage disequilibrium values. Some putative transcription factors were predicted for binding sites where the SNPs considered are located. Since the expression rate of the gene under alternative environmental conditions seems to depend on the composed genotype of several SNPs located at its promoter, a cooperative regulation of the transcription of the HSP90AA1 gene could be hypothesized. Nevertheless epigenetic regulation mechanisms cannot be discarded.