A genome-wide analysis of copy number variation in Murciano-Granadina goats.

BACKGROUND: In this work, our aim was to generate a map of the copy number variations (CNV) segregating in a population of Murciano-Granadina goats, the most important dairy breed in Spain, and to ascertain the main biological functions of the genes that map to copy number variable regions. RESULTS: Using a dataset that comprised 1036 Murciano-Granadina goats genotyped with the Goat SNP50 BeadChip, we were able to detect 4617 and 7750 autosomal CNV with the PennCNV and QuantiSNP software, respectively. By applying the EnsembleCNV algorithm, these CNV were assembled into 1461 CNV regions (CNVR), of which 486 (33.3% of the total CNVR count) were consistently called by PennCNV and QuantiSNP and used in subsequent analyses. In this set of 486 CNVR, we identified 78 gain, 353 loss and 55 gain/loss events. The total length of all the CNVR (95.69 Mb) represented 3.9% of the goat autosomal genome (2466.19 Mb), whereas their size ranged from 2.0 kb to 11.1 Mb, with an average size of 196.89 kb. Functional annotation of the genes that overlapped with the CNVR revealed an enrichment of pathways related with olfactory transduction (fold-enrichment = 2.33, q-value = 1.61 × 10-10), ABC transporters (fold-enrichment = 5.27, q-value = 4.27 × 10-04) and bile secretion (fold-enrichment = 3.90, q-value = 5.70 × 10-03). CONCLUSIONS: A previous study reported that the average number of CNVR per goat breed was ~ 20 (978 CNVR/50 breeds), which is much smaller than the number we found here (486 CNVR). We attribute this difference to the fact that the previous study included multiple caprine breeds that were represented by small to moderate numbers of individuals. Given the low frequencies of CNV (in our study, the average frequency of CNV is 1.44%), such a design would probably underestimate the levels of the diversity of CNV at the within-breed level. We also observed that functions related with sensory perception, metabolism and embryo development are overrepresented in the set of genes that overlapped with CNV, and that these loci often belong to large multigene families with tens, hundreds or thousands of paralogous members, a feature that could favor the occurrence of duplications or deletions by non-allelic homologous recombination.

A genome-wide association analysis for body, udder, and leg conformation traits recorded in Murciano-Granadina goats.

Morphological traits are of great importance to dairy goat production given their effect on phenotypes of economic interest. However, their underlying genomic architecture has not yet been extensively characterized. Herein, we aimed to identify genomic regions associated with body, udder, and leg conformation traits recorded in 825 Murciano-Granadina goats. We genotyped this resource population using the GoatSNP50 BeadChip (Illumina Inc., San Diego, CA) and performed genome-wide association analyses using the GEMMA software. We found 2 genome-wide significant associations between markers rs268273468 [Capra hircus (CHI) 16:69617700] and rs268249346 (CHI 28:18321523) and medial suspensory ligament. In contrast, we did not detect any genome-wide significant associations for body and leg traits. Moreover, we found 12, 19, and 7 chromosome-wide significant associations for udder, body, and leg traits, respectively. Comparison of our data with previous studies revealed a low level of positional concordance between regions associated with morphological traits. In addition to technical factors, this lack of concordance could be due to a substantial level of genetic heterogeneity among breeds or to the strong polygenic background of morphological traits, which makes it difficult to detect genetic factors that have small phenotypic effects.

Assessing the levels of intraspecific admixture and interspecific hybridization in Iberian wild goats (Capra pyrenaica).

Iberian wild goats (Capra pyrenaica, also known as Iberian ibex, Spanish ibex, and Spanish wild goat) underwent strong genetic bottlenecks during the 19th and 20th centuries due to overhunting and habitat destruction. From the 1970s to 1990s, augmentation translocations were frequently carried out to restock Iberian wild goat populations (very often with hunting purposes), but they were not systematically planned or recorded. On the other hand, recent data suggest the occurrence of hybridization events between Iberian wild goats and domestic goats (Capra hircus). Augmentation translocations and interspecific hybridization might have contributed to increase the diversity of Iberian wild goats. With the aim of investigating this issue, we have genotyped 118 Iberian wild goats from Tortosa-Beceite, Sierra Nevada, Muela de Cortes, Gredos, Batuecas, and Ordesa and Monte Perdido by using the Goat SNP50 BeadChip (Illumina). The analysis of genotypic data indicated that Iberian wild goat populations are strongly differentiated and display low diversity. Only three Iberian wild goats out from 118 show genomic signatures of mixed ancestry, a result consistent with a scenario in which past augmentation translocations have had a limited impact on the diversity of Iberian wild goats. Besides, we have detected eight Iberian wild goats from Tortosa-Beceite with signs of domestic goat introgression. Although rare, hybridization with domestic goats could become a potential threat to the genetic integrity of Iberian wild goats; hence, measures should be taken to avoid the presence of uncontrolled herds of domestic or feral goats in mountainous areas inhabited by this iconic wild ungulate.

Identification of circular RNAs in porcine sperm and evaluation of their relation to sperm motility.

Circular RNAs (circRNAs) are emerging as a novel class of noncoding RNAs which potential role as gene regulators is quickly gaining interest. circRNAs have been studied in different tissues and cell types across several animal species. However, a thorough characterization of the circRNAome in ejaculated sperm remains unexplored. In this study, we profiled the sperm circRNA catalogue using 40 porcine ejaculates. A complex population of 1,598 circRNAs was shared in at least 30 of the 40 samples. Generally speaking, the predicted circRNAs presented low abundances and were tissue-specific. Around 80% of the circRNAs identified in the boar sperm were reported as novel. Results from abundance correlation between circRNAs and miRNAs together with the prediction of microRNA (miRNA) target sites in circRNAs suggested that circRNAs may act as miRNA sponges. Moreover, we found significant correlations between the abundance of 148 exonic circRNAs and sperm motility parameters. Two of these correlations, involving ssc_circ_1458 and ssc_circ_1321, were confirmed by RT-qPCR using 36 additional samples with extreme and opposite sperm motility values. Our study provides a thorough characterization of circRNAs in sperm and suggests that circRNAs hold potential as noninvasive biomarkers for sperm quality and male fertility.

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

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.