Genomic and phenotypic analysis of invasive Streptococcus suis isolated in Spain reveals genetic diversification and associated virulence traits.

Streptococcus suis is a zoonotic pathogen that causes a major health problem in the pig production industry worldwide. Spain is one of the largest pig producers in the world. This work aimed to investigate the genetic and phenotypic features of invasive S. suis isolates recovered in Spain. A panel of 156 clinical isolates recovered from 13 Autonomous Communities, representing the major pig producers, were analysed. MLST and serotyping analysis revealed that most isolates (61.6%) were assigned to ST1 (26.3%), ST123 (18.6%), ST29 (9.6%), and ST3 (7.1%). Interestingly, 34 new STs were identified, indicating the emergence of novel genetic lineages. Serotypes 9 (27.6%) and 1 (21.8%) prevailed, followed by serotypes 7 (12.8%) and 2 (12.2%). Analysis of 13 virulence-associated genes showed significant associations between ST, serotype, virulence patterns, and clinical features, evidencing particular virulence traits associated with genetic clusters. The pangenome was generated, and the core genome was distributed in 7 Bayesian groups where each group included a variable set of over- and under-represented genes of different categories. The study provides comprehensive data and knowledge to improve the design of new vaccines, antimicrobial treatments, and bacterial typing approaches.

Blood transcriptome of Rasa Aragonesa rams with different sexual behavior phenotype reveals CRYL1 and SORCS2 as genes associated with this trait.

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.

Analyzing ram sexual behavior via blood transcriptome profiling can help to identify in vivo sexual behavior biomarkers as an innovative alternative to invasive and time-consuming methods in farms. Using RNA-sequencing technique, we compared 12 Rasa Aragonesa rams with different sexual behavior (6 sexually active and 6 nonactive) to identify differentially expressed genes (DEGs) in peripheral blood putatively responsible of libido differences between rams. Comparative analysis revealed four candidate genes and several signaling pathways related to sexual behavior such as lysosome, and positive regulation of the extracellular signal-regulated kinase 1/2 (ERK1 and ERK2) cascade. This data will be helpful for further investigations to understand the differences of sheep sexual behavior.

Exploring differentially expressed genes in hypothalamic, pars tuberalis and pineal gland transcriptomes in different sexual behavior phenotypes in rams using RNA-Seq.

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 like 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.

Male sexual behavior is a key factor in reproduction, especially in seasonal breeders such as sheep. The identification of differentially expressed genes (DEGs) in brain regions involved in male reproduction and sexual behavior between rams with different sexual activity by RNA high-throughput sequencing can provide useful information to the sheep meat industry. This work aimed to determine the possible molecular mechanisms underlying the sexual behavior of Rasa Aragonesa rams by investigating transcriptional changes in the hypothalamus (HT), pars tuberalis (PT) and pineal gland (PG) between active (A) and nonactive (NA) rams. Comparative analysis revealed 103 and 12 DEGs between the A vs. NA comparison in the HT and the PG, respectively, but no DEGs were found in the PT. Gene ontology (GO) enrichment analysis of DEGs in HT samples revealed significant pathways, associated mainly with neuroactive ligand-receptor interactions, and the GnRH and cAMP signaling pathways. Furthermore, gene set enrichment analysis (GSEA) detected many overrepresented pathways related to sexual behavior via an interaction network within the hypothalamic-pituitary-gonadal axis. These data will be helpful for further investigations to look for mutations or functional single nucleotide polymorphisms (SNPs) that may be used for genetic assisted selection to improve sexual behavior in sheep.

Characterization of the pars tuberalis and hypothalamus transcriptome in female sheep under different reproductive stages.

For understanding the molecular events underlying the follicular (F) and luteal (L) phases of estrous cycle, and anestrous (A) phase, the pars tuberalis (PT), and hypothalamus (HT) transcriptomes of 21 ewes were studied. In HT, 72 and 3 differential expression genes (DEGs) were found when comparing F vs. A and L vs. A, respectively. In PT, 6 and 4 DEGs were found in F vs. A and L vs. A comparisons, respectively. Enrichment analysis for DEGs between the F and A phases in the HT revealed significant clusters, mainly associated with actin-binding, and cytoskeleton, that are related to neural plasticity modulated by gonadal steroid hormones, as well as with oxytocin signaling. We found that DEGs in PT had higher differences in expression levels than those found in HT. In this sense, the ITLN was highly upregulated in the F and L vs. A phases, being MRPL57 and IRX4 highly downregulated in L vs. A comparison. The DDC gene in PT, related to LH regulation, was upregulated in the F phase. The gene set enrichment analysis (GSEA) revealed multiple pathways related to neurotransmission and neuronal plasticity. Our study reveals new candidate genes involved in the reproductive stages' transitions in seasonal sheep.

Genome-Wide Association Studies of Somatic Cell Count in the Assaf Breed.

A genome-wide association study (GWAS) was performed to identify new single nucleotide polymorphisms (SNPs) and genes associated with mastitis resistance in Assaf sheep by using the Illumina Ovine Infinium® HD SNP BeadChip (680K). In total, 6173 records from 1894 multiparous Assaf ewes with at least three test day records and aged between 2 and 7 years old were used to estimate a corrected phenotype for somatic cell score (SCS). Then, 192 ewes were selected from the top (n = 96) and bottom (n = 96) tails of the corrected SCS phenotype distribution to be used in a GWAS. Although no significant SNPs were found at the genome level, four SNPs (rs419096188, rs415580501, rs410336647, and rs424642424) were significant at the chromosome level (FDR 10%) in two different regions of OAR19. The SNP rs419096188 was located in intron 1 of the NUP210 and close to the HDAC11 genes (61 kb apart), while the other three SNPs were totally linked and located 171 kb apart from the ARPP21 gene. These three genes were related to the immune system response. These results were validated in two SNPs (rs419096188 and rs424642424) in the total population (n = 1894) by Kompetitive Allele-Specific PCR (KASP) genotyping. Furthermore, rs419096188 was also associated with lactose content.

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

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 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.

A new allele in the BMP15 gene (FecXRA) that affects prolificacy co-segregates with FecXR and FecXGR in Rasa aragonesa sheep.

A FecX-mutated allele called FecXR in the BMP15 gene has been described in Rasa aragonesa sheep. FecXR causes increased prolificacy when heterozygous and sterility when homozygous in ewes. However, highly prolific ewes without the FecXR allele have been found in this breed. Therefore, a genome-wide association study (GWAS) was performed in 158 ewes (tail H: N = 73, mean prolificacy ± standard deviation = 2.14 ± 0.26; tail L: N = 85, mean prolificacy = 1.06 ± 0.08) with the Ovine HD SNP BeadChip. In this analysis, the FecXGR allele was found to have genome-wide significance associated with prolificacy, first described in the Grivette sheep breed. We also identified a novel polymorphism in exon 2 of BMP15 in 9 high prolific ewes by Sanger sequencing. This new mutation, called FecXRA, is a SNP (Oar3.1_X: g. 50970948C > T; NM_001114767.1: c.1172C > T) that produces an amino acid substitution (ENSOART00000010201: p.T400I) that is predicted to be deleterious and to alter the predicted secondary structure of the mature protein. To confirm if this SNP had any the effect on prolificacy, we genotyped sires with known EBVs (Estimated Breeding Values), finding one hemizygous sire for the FecXRA allele with the highest EBV in the breeding program (effect on litter size at + 0.39 lamb per lambing). A very low frequency, ranging from 0.13 to 2%, was found for the FecXGR and FecXRA alleles in 3428 animals belonging to four different flocks. Finally, an association study was performed to validate and quantify the effects of the FecXGR and FecXRA alleles. Significant increased prolificacy of 0.52 ± 0.05, 0.42 ± 0.05 and 0.32 ± 0.01 were found when comparing FecXGR, FecXRA and FecXR heterozygous ewes to wild type homozygous ones. These effects are of the same order of magnitude as the effect of most of other known major genes for prolificacy. Only significant differences between FecXGR and FecXR were found among the three alleles associated with increased prolificacy. However, we cannot confirm the effect of the FecXRA allele at homozygous state because we did not find any homozygous ewes. These results confirm that these three alleles in the BMP15 gene that affect prolificacy co-segregate in Rasa aragonesa sheep.

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

Sheep breeds from Mediterranean area show reproductive seasonal patterns of oestrous behaviour and ovulatory activity, mainly regulated by variation in the photoperiod. Maximal reproductive activity is associated with short days from August to March. The aim of this study therefore was, to identify new SNPs and genes associated to reproductive seasonality in sheep by using the Illumina OvineSNP50 Beadchip. A total of 239 adult Rasa Aragonesa breed ewes from one flock were controlled from January to August. Three reproductive seasonality traits were considered: the total days of anoestrus (TDA), based on weekly individual plasma progesterone levels and defined as the sum of days in anoestrus, considering anoestrus those periods with three or more consecutive P4 concentrations lower than 0.5 ng/ml; the progesterone cycling months (P4CM), defined for each ewe as the rate of cycling months between January and August based on progesterone determinations and the oestrus cycling months (OCM), defined for each ewe as the rate of months cycling between January and August based on oestrus records. Genotyping of 123 ewes was performed with the OvineSNP50 Infinium Beadchip. After the quality control (QC) performed on the raw genotypes, a total of 47,206 SNPs distributed over the 27 ovine chromosomes and 110 ewes were included in subsequent analyses. Principal component analysis revealed a substructure within the total dataset and identified 4 principal clusters in the experimental flock. None of the SNPs overcame the genome-wide significance level (P = 1.06 × 10-6). However, the SNPs OAR4_66002395 (9.41E-6), and OAR8_25877010 (1.86E-5) reached the genome-wide suggestive significance level (set to 2.32 × 10-5) for TDA and P4CM traits, respectively, while OAR23_14608581 was significant for both TDA (2.02E-5) and P4CM (1.05E-5) traits. Five SNPs evidenced association at chromosome-wise level: SNPs OAR4_66002395, OAR23_14608581 and s20800 (DTA), and OAR8_25877010, OAR23_14608581 and s48474 (P4CM). Several genes related to circadian and circannual rhythms were found close to these SNPs: NPSR1 (SNP OAR4_66002395), HS3ST5 (SNP OAR8_25877010), RPTOR (SNP s48474), and NPTX1 (SNP s48474) and could be considered as candidate gene related to TDA and P4CM traits.

The relationship between muscle α-tocopherol concentration and meat oxidation in light lambs fed vitamin E supplements prior to slaughter.

BACKGROUND: The use of concentrates supplemented with α-tocopherol in animals is an effective method to reduce the oxidative processes that occur in meat products. The high cost of α-tocopherol requires accurate feeding, so it is necessary to define the minimum period of α-tocopherol concentrate supplementation that will ensure an acceptable meat quality. Indoor concentrate-fed light lambs (n = 35) were supplemented with 500 mg dl-α-tocopheryl acetate (VE) kg(-1) concentrate for a period of between 4 and 28 days before being slaughtered at 22-24 kg body weight. Control lambs (n = 12) were not supplemented with α-tocopherol. RESULTS: The α-tocopherol content in both plasma and muscle tissues increased significantly with the length of supplementation (P < 0.001). The thiobarbituric acid-reactive substance (TBARS) concentration in meat decreased exponentially when the muscle α-tocopherol concentration was increased to 0.61-0.90 mg kg(-1) fresh meat (P < 0.05). After 7 days of display, the formation of metmyoglobin (MMb) decreased significantly as the α-tocopherol content increased to 0.31-0.60 mg kg(-1) meat (P < 0.05). CONCLUSION: A range of 0.61-0.90 mg α-tocopherol kg(-1) fresh meat protected fresh lamb meat from lipid oxidation and MMb formation. This level can be achieved by supplementation with 500 mg VE kg(-1) concentrate for a period of 7-14 days before slaughter.

Investigation of the genetic diversity among native Turkish sheep breeds using mtDNA polymorphisms.

A total of 135 unrelated sheep from nine Turkish native sheep breeds (Daglic, Kivircik, Imroz, Chios, Morkaraman, Ivesi, Hemsin, Karayaka and Akkaraman) were investigated to determinate the maternal genetic diversity using a sequence of a 531-bp segment of the mtDNA control region. Analysis of the mtDNA control region sequence revealed 63 haplotypes and 53 polymorphic sites. Haplotype diversity, nucleotide diversity and the average number of nucleotide differences were estimated to be 0.9496 ± 0.011, 0.01407 ± 0.00060 and 7.456, respectively. The sequence analysis also revealed high level of genetic diversity among the native Turkish breeds. These breeds were grouped into three major maternal haplogroups: A, B and C, with one animal belonging from the Akkaraman breed to the rare haplogroup E. Irregular shape of mismatch distribution of haplogroup C could be an indicator that haplogroup C may represent different haplogroups. Contrarily to previous studies carried out on Turkish native breeds, majority of animals grouped in haplogroup A in the present study. This result and the irregular shape of mismatch curve of haplogroup C indicate that genetic structure of Turkish native sheep breeds could be more complicated than it is thought.

Influence of lairage time on some welfare and meat quality parameters in pigs.

A total of 150 Large White cross Landrace pigs (110-120 kg) of both sexes were used to investigate the effects of three different lairage times (0 h, 3 h, 9 h). Blood samples were collected at exanguination and cortisol, glucose, lactate, muscle enzymes and haematological parameters were determined. Post-mortem measurements of muscle pH were taken at 20 min, 2 h and 24 h from Longissimus thoracis and Semimembranosus. Lairage time showed a significant effect on pH24, internal muscle reflectance using the fibre optic probe (FOP24), red blood cells, neutrophils and lymphocytes, glucose and enzymatic activities. Changes in blood profile and meat quality parameters indicated that three hours of lairage in Spanish commercial conditions may reduce the amount of stress exhibited by pigs at slaughter and better meat quality can be obtained. No lairage or an excessively long lairage period without food may compromise animal welfare and meat quality.