A genome-wide association study (GWAS) was used to identify genomic regions influencing seasonality reproduction traits in Rasa Aragonesa sheep. Three traits associated with either ovarian function based on blood progesterone levels (total days of anoestrus and progesterone cycling months) or behavioral signs of oestrous (oestrous cycling months) were studied. The GWAS included 205 ewes genotyped using the 50k and 680k Illumina Ovine Beadchips. Only one SNP associated with the progesterone cycling months overcame the genome-wide significance level (rs404991855). Nine SNPs exhibited significant associations at the chromosome level, being the SNPs rs404991855 and rs418191944, that are located in the CD226 molecule (CD226) gene, associated with the three traits. This gene is related to reproductive diseases. Two other SNPs were located close to the neuropeptide Y (NPY) gene, which is involved in circadian rhythms. To validate the GWAS, partial characterization of both genes by Sanger sequencing, and genotyping of two synonymous and two nonsynonymous SNPs in the NPY and CD226 genes, respectively, were performed. SNP association analysis showed that only SNP rs404360094 in the exon 3 of the CD226 gene, which produces an amino acid substitution from asparagine (uncharged polar) to aspartic acid (acidic), was associated with the three seasonality traits. Our results suggest that the CD226 gene may be involved in the reproductive seasonality in Rasa Aragonesa.
The aim of this study was to characterize and identify causative polymorphisms in the leptin receptor (LEPR) gene responsible for the seasonal variation of reproductive traits in sheep. Three reproductive seasonality traits were studied: the total days of anoestrous (TDA), the progesterone cycling months (P4CM) and the oestrous cycling months (OCM). In total, 18 SNPs were detected in 33 ewes with extreme values for TDA and OCM. Six SNPs were non-synonymous substitutions and two of them were predicted in silico as deleterious: rs596133197 and rs403578195. These polymorphisms were then validated in 239 ewes. The SNP rs403578195, located in exon 8 and leading to a change of alanine to glycine (Ala284Gly) in the extracellular domain of the protein, was associated with the OCM trait, being the G allele associated with a decrease of 12 percent of the OCM trait. Haplotype analyses also suggested the involvement of other non-synonymous SNP located in exon 20 (rs405459906). This SNP also produces an amino acid change (Lys1069Glu) in the intracellular domain of the protein and segregates independently of rs403578195. These results confirm for the first time the role of the LEPR gene in sheep reproductive seasonality.
Reproductive seasonality is a limiting factor in sheep production. Sexual behavior is a key element in reproductive efficiency, and this function is regulated by the hypothalamus-pituitary-gonadal (HPG) axis. To understand the mechanisms of sexual behavior, transcriptomic sequencing technology was used to identify differentially expressed genes (DEGs) in the hypothalamus (HT), pars tuberalis (PT) and pineal gland (PG) in Rasa Aragonesa rams with different sexual behavior. Bioinformatics analysis of the 16,401 identified genes by RNA-Seq revealed 103 and 12 DEGs in the HT and the PG, respectively, at a false discovery rate (FDR) of 5% with an absolute value of expression ≥ 1 (log2FC). However, no DEGs were found in the PT. Functional annotation and pathway enrichment analysis showed that DEGs of HT were enriched mainly in neuroactive ligand-receptor interactions and signaling pathways, including notable candidate genes such as MTNR1A, CHRNA2, FSHB, LHB, GNRHR, AVP, PRL, PDYN, CGA, GABRD and TSHB, which play a crucial role in sexual behavior. The GnRH and cAMP signaling pathways were also highlighted. In addition, gene set enrichment analysis (GSEA) identified potential pathways, dominated mainly by biological process category, that could be responsible for the differences in sexual behavior observed in rams. The intracellular protein transport and pattern specification process were enriched within the PT and the transcription factor binding and protein ubiquitination pathways for the PG. Thus, these pathways together, may play an important role in the regulation of the sexual behavior in Rasa Aragonesa rams through the hypothalamic-pituitary-gonadal axis. The validation of 5 DEGs using reverse transcription quantitative polymerase chain reaction (RT- qPCR) showed expression patterns similar to the found with RNA-Seq. Overall, these results contribute to understanding the genomic basis of sexual behavior in rams. Our study demonstrates that multiple networks and pathways orchestrate sexual behavior in sheep.
Reproductive fitness of rams is seasonal, showing the highest libido during short days coinciding with the ovarian cyclicity resumption in the ewe. However, the remarkable variation in sexual behavior between rams impair farm efficiency and profitability. Intending to identify in vivo sexual behavior biomarkers that may aid farmers to select active rams, transcriptome profiling of blood was carried out by analyzing samples from 6 sexually active (A) and 6 nonactive (NA) Rasa Aragonesa rams using RNA-Seq technique. A total of 14,078 genes were expressed in blood but only four genes were differentially expressed (FDR < 0.10) in the A vs. NA rams comparison. The genes, acrosin inhibitor 1 (ENSOARG00020023278) and SORCS2 were upregulated (log2FC > 1) in active rams, whereas the CRYL1 and immunoglobulin lambda-1 light chain isoform X47 (ENSOARG00020025518) genes were downregulated (log2FC < -1) in this same group. Gene set Enrichment Analysis (GSEA) identified 428 signaling pathways, predominantly related to biological processes. The lysosome pathway (GO:0005764) was the most enriched, and may affect fertility and sexual behavior, given the crucial role played by lysosomes in steroidogenesis, being the SORCS2 gene related to this signaling pathway. Furthermore, the enriched positive regulation of ERK1 and ERK2 cascade (GO:0070374) pathway is associated with reproductive phenotypes such as fertility via modulation of hypothalamic regulation and GnRH-mediated production of pituitary gonadotropins. Furthermore, external side of plasma membrane (GO:0009897), fibrillar center (GO:0001650), focal adhesion (GO:0005925) and lamellipodium (GO:0030027) pathways were also enriched, suggesting that some molecules of these pathways might also be involved in rams’ sexual behavior. These results provide new clues for understanding the molecular regulation of sexual behavior in rams. Further investigations will be needed to confirm the functions of SORCS2 and CRYL1 in relation to sexual behavior.
Haplotypes from sequencing data may improve the prediction accuracy in genomic evaluations as haplotypes are in stronger linkage disequilibrium with quantitative trait loci than markers from SNP chips. This study focuses first, on the creation of haplotypes in a population sample of 450 Holstein animals, with full-sequence data from the 1000 bull genomes project; and second, on incorporating them into the whole genome prediction model. In total, 38,319,258 SNPs (and indels) from Next Generation Sequencing were included in the analysis. After filtering variants with minor allele frequency (MAF< 0.025) 13,912,326 SNPs were available for the haplotypes extraction with findhap.f90. The number of SNPs in the haploblocks was on average 924 SNP (166,552 bp). Unique haplotypes were around 97% in all chromosomes and were ignored leaving 153,428 haplotypes. Estimated haplotypes had a large contribution to the total variance of genomic estimated breeding values for kilogram of protein, Global Type Index, Somatic Cell Score and Days Open (between 32 and 99.9%). Haploblocks containing haplotypes with large effects were selected by filtering for each trait, haplotypes whose effect was larger/lower than the mean plus/minus 3 times the standard deviation (SD) and 1 SD above the mean of the haplotypes effect distribution. Results showed that filtering by 3 SD would not be enough to capture a large proportion of genetic variance, whereas filtering by 1 SD could be useful but model convergence should be considered. Additionally, sequence haplotypes were able to capture additional genetic variance to the polygenic effect for traits undergoing lower selection intensity like fertility and health traits.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.