Genome-wide association studies for reproductive seasonality traits in Rasa Aragonesa sheep breed

Complex Evolutionary History With Extensive Ancestral Gene Flow in an African Primate Radiation

Understanding the drivers of speciation is fundamental in evolutionary biology, and recent studies highlight hybridization as an important evolutionary force. Using whole-genome sequencing data from 22 species of guenons (tribe Cercopithecini), one of the world's largest primate radiations, we show that rampant gene flow characterizes their evolutionary history and identify ancient hybridization across deeply divergent lineages that differ in ecology, morphology, and karyotypes. Some hybridization events resulted in mitochondrial introgression between distant lineages, likely facilitated by cointrogression of coadapted nuclear variants. Although the genomic landscapes of introgression were largely lineage specific, we found that genes with immune functions were overrepresented in introgressing regions, in line with adaptive introgression, whereas genes involved in pigmentation and morphology may contribute to reproductive isolation. In line with reports from other systems that hybridization might facilitate diversification, we find that some of the most species-rich guenon clades are of admixed origin. This study provides important insights into the prevalence, role, and outcomes of ancestral hybridization in a large mammalian radiation.

Whole-genome resequencing of the native sheep provides insights into the microevolution and identifies genes associated with reproduction traits

BACKGROUND

Sheep genomes undergo numerous genes losses, gains and mutation that generates genome variability among breeds of the same species after long time natural and artificial selection. However, the microevolution of native sheep in northwest China remains elusive. Our aim was to compare the genomes and relevant reproductive traits of four sheep breeds from different climatic environments, to unveil the selection challenges that this species cope with, and the microevolutionary differences in sheep genomes. Here, we resequenced the genomes of 4 representative sheep breeds in northwest China, including Kazakh sheep and Duolang sheep of native breeds, and Hu sheep and Suffolk sheep of exotic breeds with different reproductive characteristics.

RESULTS

We found that these four breeds had a similar expansion experience from ~ 10,000 to 1,000,000 years ago. In the past 10,000 years, the selection intensity of the four breeds was inconsistent, resulting in differences in reproductive traits. We explored the sheep variome and selection signatures by F_ST and θπ. The genomic regions containing genes associated with different reproductive traits that may be potential targets for breeding and selection were detected. Furthermore, non-synonymous mutations in a set of plausible candidate genes and significant differences in their allele frequency distributions across breeds with different reproductive characteristics were found. We identified PAK1, CYP19A1 and PER1 as a likely causal gene for seasonal reproduction in native sheep through qPCR, Western blot and ELISA analyses. Also, the haplotype frequencies of 3 tested gene regions related to reproduction were significantly different among four sheep breeds.

CONCLUSIONS

Our results provide insights into the microevolution of native sheep and valuable genomic information for identifying genes associated with important reproductive traits in sheep.

Ram Semen Cryopreservation for Portuguese Native Breeds: Season and Breed Effects on Semen Quality Variation

The semen quality is one of the determinant factors of ram semen cryopreservation. The present retrospective study aimed to characterize the seasonal ram pattern during the year for ten Portuguese local sheep breeds, hypothesizing that the breed and season had low effects on the main spermatozoa traits. A total of 1471 ejaculates were used and evaluated (fresh semen) from 85 rams between 2004 and 2020 and re-evaluated after thawing (thawed semen). The effect of breed, season, and sperm cryopreservation on nine semen traits were evaluated. The volume per ejaculate, spermatozoa (SPZ) concentration, and total number of SPZ per ejaculate, were affected by breed (p < 0.001) but not by season (p > 0.05). As expected, the semen processing was the most significant (p < 0.001) factor of variation on seminal parameters. Moreover, breed and interactions between breed × semen processing, modulated the response of alive SPZ, abnormal morphology, head, and intermediate piece defects. In fresh semen, season only affected the intermediate piece defects due to the highest percentage observed between February and April period in some breeds. Overall, and despite the mentioned particularities, there were similarities among the ten local breeds. We also concluded that the seasonal effect on ejaculate and SPZ traits is not significant in our region. These local ram breeds have low seasonality and can be employed in natural mating as well as semen donors for cryopreservation and assisted reproductive biotechnologies during the whole year at our latitude.

Detection of selection signatures in the genome of a farmed population of anadromous rainbow trout (Oncorhynchus mykiss)

Domestication processes and artificial selection are likely to leave signatures that can be detected at a molecular level in farmed rainbow trout (Oncorhynchus mykiss). These signatures of selection are genomic regions that contain functional genetic variants conferring a higher fitness to their bearers. We genotyped 749 rainbow trout from a commercial population using a rainbow trout Axiom 57 K SNP array panel and identified putative genomic regions under selection using the pcadapt, Composite Likelihood Ratio (CLR) and Integrated Haplotype Score (iHS) methods. After applying quality-control pipelines and statistical analyses, we detected 12, 96 and 16 SNPs putatively under selection, associated with 96, 781 and 115 candidate genes, respectively. Several of these candidate genes were associated with growth, early development, reproduction, behavior and immune system traits. In addition, some of the SNPs were found in interesting regions located in autosomal inversions on Omy05 and Omy20. These findings could represent a genome-wide map of selection signatures in farmed rainbow trout and could be important in explaining domestication and selection for genetic traits of commercial interest.

The GWAS Analysis of Body Size and Population Verification of Related SNPs in Hu Sheep

Body size is an important indicator of growth and health in sheep. In the present study, we performed Genome-Wide Association Studies (GWAS) to detect significant single-nucleotide polymorphisms (SNPs) associated with Hu sheep's body size. After genotyping parental (G1) and offspring (G2) generation of the nucleus herd for meat production of Hu sheep and conducting GWAS on the body height, chest circumference, body length, tail length, and tail width of the two groups, 5 SNPs associated with body height and 4 SNPs correlated with chest circumference were identified at the chromosomal significance level. No SNPs were significantly correlated to body length, tail length, and width. Four out of the 9 SNPs were found to be located within the 4 genes. KITLG and CADM2 are considered as candidate functional genes related to body height; MCTP1 and COL4A6 are candidate functional genes related to chest circumference. The 9 SNPs found in GWAS were verified using the G3 generation of the nucleus herd for meat production. Nine products were amplified around the 9 sites, and 29 SNPs were found; 3 mutation sites, G > C mutation at 134 bp downstream of s554331, T > G mutation at 19 bp upstream of s26859.1, and A > G mutation at 81 bp downstream of s26859.1, were significantly correlated to the body height. Dual-luciferase reporter gene experiments showed that the 3 SNPs could significantly impact dual-luciferase and gene transcription activity.

Genome-Wide Association Study Demonstrates the Role Played by the CD226 Gene in Rasa Aragonesa Sheep Reproductive Seasonality

Simple Summary

To elucidate the genetic basis of reproductive seasonality in Rasa Aragonesa sheep breed, we performed a genome-wide association study (GWAS) in order to detect single nucleotide polymorphisms (SNPs) or regions associated with traits related to ovarian function and behavioural signs of estrous. The GWAS included 205 ewes with genotypes for 583882 SNPs. Only one SNP overcame the genome-wide significance level. Nine potential SNPs overcame the chromosome-wise significance level (FDR 10%). Gene annotation demonstrated that CD226molecule (CD226) and neuropeptide Y (NPY) genes that could be involved in reproductive seasonality were close to the significant SNPs. To validate the results, we sequenced the entire coding region of the NPY gene and four exons of the CD226 gene to search for polymorphisms that could be involved in the phenotypes studied. Two synonymous and two nonsynonymous SNPs in the NPY and CD226 genes, respectively, were genotyped in the whole population. We demonstrated that the AA genotype of the SNP rs404360094 located in exon 3 of the CD226 gene was associated with higher and lower total days of anoestrus and oestrous cycling months, respectively. Therefore, this SNP could be utilized as a genetic marker for assisted selection marker to reduce seasonality.

Abstract

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.

Genome-Wide Analyses Reveal Genetic Convergence of Prolificacy between Goats and Sheep

The litter size of domestic goats and sheep is an economically important trait that shows variation within breeds. Strenuous efforts have been made to understand the genetic mechanisms underlying prolificacy in goats and sheep. However, there has been a paucity of research on the genetic convergence of prolificacy between goats and sheep, which likely arose because of similar natural and artificial selection forces. Here, we performed comparative genomic and transcriptomic analyses to identify the genetic convergence of prolificacy between goats and sheep. By combining genomic and transcriptomic data for the first time, we identified this genetic convergence in (1) positively selected genes (CHST11 and SDCCAG8), (2) differentially expressed genes (SERPINA14, RSAD2, and PPIG at follicular phase, and IGF1, GPRIN3, LIPG, SLC7A11, and CHST15 at luteal phase), and (3) biological pathways (genomic level: osteoclast differentiation, ErbB signaling pathway, and relaxin signaling pathway; transcriptomic level: the regulation of viral genome replication at follicular phase, and protein kinase B signaling and antigen processing and presentation at luteal phase). These results indicated the potential physiological convergence and enhanced our understanding of the overlapping genetic makeup underlying litter size in goats and sheep.

The LEPR Gene Is Associated with Reproductive Seasonality Traits in Rasa Aragonesa Sheep

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.

Genome-Wide Association Study Reveals Candidate Genes for Litter Size Traits in Pelibuey Sheep

The Pelibuey sheep has adaptability to climatic variations, resistance to parasites, and good maternal ability, whereas some ewes present multiple births, which increases the litter size in farm sheep. The litter size in some wool sheep breeds is associated with the presence of mutations, mainly in the family of the transforming growth factor β (TGF-β) genes. To explore genetic mechanisms underlying the variation in litter size, we conducted a genome-wide association study in two groups of Pelibuey sheep (multiparous sheep with two lambs per birth vs. uniparous sheep with a single lamb at birth) using the OvineSNP50 BeadChip. We identified a total of 57 putative SNPs markers (p < 3.0 × 10^-3, Bonferroni correction). The candidate genes that may be associated with litter size in Pelibuey sheep are CLSTN2, MTMR2, DLG1, CGA, ABCG5, TRPM6, and HTR1E. Genomic regions were also identified that contain three quantitative trait loci (QTLs) for aseasonal reproduction (ASREP), milk yield (MY), and body weight (BW). These results allowed us to identify SNPs associated with genes that could be involved in the reproductive process related to prolificacy.

Genomic Approaches Identify Novel Gene Associations with Out of Season Lambing in Sheep

Sheep are seasonally polyestrous, traditionally breeding when the day length shortens in the autumn. The changing photoperiod stimulates reproductive hormones through a series of chemical pathways, ultimately leading to cyclicity. Some breeds of sheep, such as the Polypay and Dorset, have been selected for reduced seasonality and can lamb year-round. Despite this selection, there is still variation within these breeds in the ability to lamb out of season. The identification of out of season lambing quantitative trait loci has the potential to improve genetic progress using genomic selection schemes. Association studies, fixation index (FST), and runs of homozygosity (ROH) were evaluated to identify regions of the genome that influence the ability of ewes to lamb out of season. All analyses used genotypic data from the Illumina Ovine HD beadchip. Genome-wide associations were tested both across breeds in 257 ewes and within the Dorset and Polypay breeds. FST was measured across breeds and between UK and US Dorsets to assess population differences. ROH were estimated in ewes to identify homozygous regions contributing to out of season lambing. Significant associations after multiple testing correction were found through these approaches, leading to the identification of several candidate genes for further study. Genes involved with eye development, reproductive hormones, and neuronal changes were identified as the most promising for influencing the ewe’s ability to lamb year-round. These candidate genes could be advantageous for selection for improved year-round lamb production and provide better insight into the complex regulation of seasonal reproduction.

Assessing genetic architecture and signatures of selection of dual purpose Gir cattle populations using genomic information

Gir is one of the main cattle breeds raised in tropical South American countries. Strong artificial selection through its domestication resulted in increased genetic differentiation among the countries in recent years. Over the years, genomic studies in Gir have become more common. However, studies of population structure and signatures of selection in divergent Gir populations are scarce and need more attention to better understand genetic differentiation, gene flow, and genetic distance. Genotypes of 173 animals selected for growth traits and 273 animals selected for milk production were used in this study. Clear genetic differentiation between beef and dairy populations was observed. Different criteria led to genetic divergence and genetic differences in allele frequencies between the two populations. Gene segregation in each population was forced by artificial selection, promoting isolation, and increasing genetic variation between them. Results showed evidence of selective forces in different regions of the genome. A total of 282 genes were detected under selection in the test population based on the fixation index (Fst), integrated haplotype score (iHS), and cross-population extend haplotype homozygosity (XP-EHH) approaches. The QTL mapping identified 35 genes associated with reproduction, milk composition, growth, meat and carcass, health, or body conformation traits. The investigation of genes and pathways showed that quantitative traits associated to fertility, milk production, beef quality, and growth were involved in the process of differentiation of these populations. These results would support further investigations of population structure and differentiation in the Gir breed.