Molecular characterisation, genetic variability and detection of a functional polymorphism influencing the promoter activity of OXT gene in goat and sheep

A Comprehensive Analysis of CSN1S2 I and II Transcripts Reveals Significant Genetic Diversity and Allele-Specific Exon Skipping in Ragusana and Amiatina Donkeys

The αs2-casein is a phosphoprotein secreted in the milk of most mammals, and it is the most hydrophilic of all caseins. Contrary to genes found in ruminants, in donkeys two different encoding genes for donkey αs2-casein (CSN1S2 I and CSN1S2 II) have been identified. However, unlike in ruminants, the variability at these loci has not been characterized in detail in donkeys until now. In this study, we analyze the transcript profile of the donkey CSN1S2 I and CSN1S2 II genes, and we identify and describe the variability of these loci in the Ragusana and Amiatina breeds reared in Italy. The analysis of the CSN1S2 I Reverse Transcriptase-Polymerase Chain Reaction (RT-PCR) products and subsequent sequencing showed, in addition to correctly spliced mRNA, seven other minor mRNAs resulting from differential splicing events involving, in various combinations, entire exons (4, 5, 6, and 11), parts of exons (5' or 3' end of exon 17), or the recognition of intronic sequences as an exon (exon 12'). Similarly, the transcription analysis of the CSN1S2 II gene revealed a remarkable variability in splicing events, mainly concerning the alternative insertion of an extra exon 7 (named 7'); the first 33 bp of exon 13; or the alternative skipping of exons 9, 10, 11, 12, and 15, and their combinations. At the mRNA level for CSN1S2 I, seven SNPs were observed, five of which led to amino acid changes: p.T73>A, p.I109>V, p.I130>V, p.I146>T, and p.D217>Y. Similarly, nine SNPs were observed at the CSN1S2 II locus, seven of which are non-synonymous: p.L63>F, p.H70>Q, p.D90>N, p.129A>T, p.H131>Y, p.E144>G, and p.F157>S. In addition, the DNA sequencing of exon 17 and flanking introns of the CSN1S2 I gene revealed a G>A transition at the splice acceptor site of CSN1S2 I exon 17 (FM946022.1:c.375-1G>A), resulting in an allele-specific skipping of the first 15 nucleotides of this exon, which encode the peptide 176NKINQ180, and the recognition of an in-frame cryptic splicing acceptor site: arAACAAAATCAACCAG. A genotyping method based on restriction fragment length polymorphism (XbaI PCR-RFLP) was set up for this SNP. In the total population studied (105 Ragusana and 14 Amiatina donkeys), the A allele had a frequency of 0.2437 with no evidence of deviation from the Hardy-Weinberg equilibrium. This study adds new knowledge regarding the genetic variability of αs2-caseins in donkeys and may contribute significantly to the genetic improvement of milk production for this species.

Genetic variability among and within domestic Old and New World camels at the α-lactalbumin gene (LALBA) reveals new alleles and polymorphisms responsible for differential expression

α-Lactalbumin (α-LA), which is encoded by the LALBA gene, is a major whey protein that binds to Ca2+ and facilitates lactose synthesis as a regulatory subunit of the synthase enzyme complex. In addition, it has been shown to play central roles in immune modulation, cell-growth regulation, and antimicrobial activity. In this study, a multitechnical approach was used to fully characterize the LALBA gene and its variants in both coding and regulatory regions for domestic camelids (dromedary, Bactrian camel, alpaca, and llama). The gene analysis revealed a conserved structure among the camelids, but a slight difference in size (2,012 bp on average) due to intronic variations. Promoters were characterized for the transcription factor binding sites (11 found in total). Intraspecies sequence comparison showed 36 SNPs in total (2 in the dromedary, none in the Bactrian camel, 22 in the alpaca, and 12 in the llama), whereas interspecies comparison showed 86 additional polymorphic sites. Eight SNPs were identified as trans-specific polymorphisms, and 2 of them (g.112A>G and g.1229A>G) were particularly interesting in the New World camels. The first creates a new binding site for transcription factor SP1. An enhancing effect of the g.112G variant on the expression was demonstrated by 3 independent pGL3 gene reporter assays. The latter is responsible for the p.78Ile>Val AA replacement and represents novel allelic variants (named LALBA A and B). A link to protein variants has been established by isoelectric focusing (IEF), and bioinformatics analysis revealed that carriers of valine (g.1229G) have a higher glycosylation rate. Genotyping methods based on restriction fragment length polymorphism (PCR-RFLP) were set up for both SNPs. Overall, adenine was more frequent (0.54 and 0.76) at both loci. Four haplotypes were found, and the AA and GA were the most common with a frequency of 0.403 and 0.365, respectively. Conversely, a putative biological gain characterizes the haplotype GG. Therefore, opportunities for rapid directional selection can be realized if this haplotype is associated with favorable milk protein properties. This study adds knowledge at the gene and protein level for α-LA (LALBA) in camelids and importantly contributes to a relatively unexplored research area in these species.

A novel deleterious oxytocin variant is associated with the lower twinning ratio in Awassi ewes

This study aimed to assess the possible association of oxytocin (OXT) gene with reproductive traits in two groups of Awassi ewes that differ in their reproductive potentials. Sheep were genotyped using PCR-single-stranded conformation polymorphism approach. Three genotypes were detected in exon 2, CC, CA, and AA, and a novel SNP was identified with a missense effect on oxytocin (c.188C > A → p.Arg55Leu). A significant (p < 0.01) association of p.Arg55Leu with the twinning rate was found as ewes with AA and CA genotypes exhibited, respectively a lower twinning ratio than those with the wild-type CC genotype. The deleterious impact of p.Arg55Leu was demonstrated by all in silico tools that were utilized to assess the effect of this variant on the structure, function, and stability of oxytocin. Molecular docking showed that p.Arg55Leu caused a dramatic alteration in the binding of oxytocin with its receptor and reduced the number of interacted amino acids between them. Our study suggests that ewes with AA and CA genotypes showed a lower reproductive performance due to the presence of p.Arg55Leu, which caused damaging impacts on oxytocin and is binding with the OXT receptor. The utilization of the p.Arg55Leu could be useful for improving Awassi reproductive potential.

The relationship between OXT gene polymorphisms and reproductive hormones in pregnant and lactating Awassi Ewes

BACKGROUND

Numerous genetic loci interact intricately to control reproduction in mammals. The oxytocin gene (OXT) is a promising candidate for reproductive traits in mammals. Previously, sheep and goats have been studied for the presence of the OXT polymorphism. As of yet, no polymorphisms have been identified in the OXT gene of Awassi sheep. Thus, this study was conducted to determine the effects of OXT polymorphism and litter size on reproductive hormones in pregnant and lactating Awassi ewes.

METHODS AND RESULTS

This study evaluated 232 ewes aged 3 and 4 years (123 single-progeny ewes and 109 twin-producing ewes). Serum was collected to measure reproductive hormones using ELISA kits manufactured by ELK Biotechnology. DNA was extracted from sheep blood for genotyping and sequencing to identify variations in OXT gene (exon 2, 266 bp). Genotyping analysis revealed three genotypes within 266 bp: CC, CA, and AA. Sequence analysis revealed a novel mutation in exon 2: 188 C > A. Statistical analysis showed significant associations between the 188 C > A SNP and phenotypic traits. Twin-pregnant ewes carrying CC genotypes had higher estrogen, progesterone, and follicle-stimulating hormone/luteinizing hormone levels (65.86 ± 3.87) (pg/mL), (6.51 ± 0.39) (ng/mL), and (20.22 ± 1.27) (ng/mL)/( 23.37 ± 2.14) (ng/mL) respectively, compared to CA and AA genotypes in the fourth month of twin-pregnant ewes compared to single-pregnant ewes.

CONCLUSIONS

This study found that the 188 C > A SNP negatively affected reproductive hormone levels in Awassi sheep. These findings provide breeders with a new insight into the sheep OXT gene, useful for future breeding.

A novel allelic donkey β-lactogobulin I protein isoform generated by a non-AUG translation initiation codon is associated with a nonsynonymous SNP

β-Lactoglobulin I (β-LG I) is one of the most important whey proteins in donkey milk. However, to our knowledge, there has been no study focusing on the full nucleotide sequences of this gene (BLG I). Current investigation of donkey BLG I gene is very limited with only 2 variants (A and B) characterized so far at the protein level. Recently, a new β-LG I variant, with a significantly higher mass (+1,915 Da) than known variants has been detected. In this study, we report the whole nucleotide sequence of the BLG I gene from 2 donkeys, whose milk samples are characterized by the β-LG I SDS-PAGE band with a normal electrophoretic mobility (18,514.25 Da, β-LG I B1 form) the first, and by the presence of a unique β-LG I band with a higher electrophoretic mobility (20,428.5 Da, β-LG I D form) the latter. A high genetic variability was found all over the 2 sequenced BLG I alleles. In particular, 16 polymorphic sites were found in introns, one in the 5' flanking region, 3 SNPs in the 5' untranslated region and one SNP in the coding region (g.458G > A) located at the 40th nucleotide of exon 2 and responsible for the AA substitutions p.Asp28 > Asn in the mature protein. Two SNPs (g.920-922CAC > TGT and g.1871G/A) were genotyped in 93 donkeys of 2 Italian breeds (60 Ragusana and 33 Amiatina, respectively) and the overall frequencies of g.920-922CAC and g.1871A were 0.3065 and 0.043, respectively. Only the rare allele g.1871A was observed to be associated with the slower migrating β-LG I. Considering this genetic diversity and those found in the database, it was possible to deduce at least 5 different alleles (BLG I A, B, B1, C, D) responsible for 4 potential β-LG I translations. Among these alleles, B1 and D are those characterized in the present research, with the D allele of real novel identification. Haplotype data analysis suggests an evolutionary pathway of donkey BLG I gene and a possible phylogenetic map is proposed. Analyses of mRNA secondary structure showed relevant changes in the structures, as consequence of the g.1871G > A polymorphism, that might be responsible for the recognition of an alternative initiation site providing an additional signal peptide. The extension of 19 AA sequence to the mature protein, corresponding to the canonical signal peptide with an additional alanine residue, is sufficient to provide the observed molecular weight of the slower migrating β-LG I encoded by the BLG I D allele.

Complementary Role of Oxytocin and Vasopressin in Cardiovascular Regulation

The neurons secreting oxytocin (OXY) and vasopressin (AVP) are located mainly in the supraoptic, paraventricular, and suprachiasmatic nucleus of the brain. Oxytocinergic and vasopressinergic projections reach several regions of the brain and the spinal cord. Both peptides are released from axons, soma, and dendrites and modulate the excitability of other neuroregulatory pathways. The synthesis and action of OXY and AVP in the peripheral organs (eye, heart, gastrointestinal system) is being investigated. The secretion of OXY and AVP is influenced by changes in body fluid osmolality, blood volume, blood pressure, hypoxia, and stress. Vasopressin interacts with three subtypes of receptors: V1aR, V1bR, and V2R whereas oxytocin activates its own OXTR and V1aR receptors. AVP and OXY receptors are present in several regions of the brain (cortex, hypothalamus, pons, medulla, and cerebellum) and in the peripheral organs (heart, lungs, carotid bodies, kidneys, adrenal glands, pancreas, gastrointestinal tract, ovaries, uterus, thymus). Hypertension, myocardial infarction, and coexisting factors, such as pain and stress, have a significant impact on the secretion of oxytocin and vasopressin and on the expression of their receptors. The inappropriate regulation of oxytocin and vasopressin secretion during ischemia, hypoxia/hypercapnia, inflammation, pain, and stress may play a significant role in the pathogenesis of cardiovascular diseases.

Complete CSN1S2 Characterization, Novel Allele Identification and Association With Milk Fatty Acid Composition in River Buffalo

The αs2-casein is one of the phosphoproteins secreted in all ruminants' milk, and it is the most hydrophilic of all caseins. However, this important gene (CSN1S2) has not been characterized in detail in buffaloes with only two alleles detected (reported as alleles A and B), and no association studies with milk traits have been carried out unlike what has been achieved for other species of ruminants. In this study, we sequenced the whole gene of two Mediterranean river buffalo homozygotes for the presence/absence of the nucleotide C (g.7539G>C) realized at the donor splice site of exon 7 and, therefore, responsible for the skipping of the same exon at mRNA level (allele B). A high genetic variability was found all over the two sequenced CSN1S2 alleles. In particular, 74 polymorphic sites were found in introns, six in the promoter, and three SNPs in the coding region (g.11072C>T, g.12803A>T, and g.14067A>G) with two of them responsible for amino acid replacements. Considering this genetic diversity, those found in the database and the SNP at the donor splice site of exon 7, it is possible to deduce at least eight different alleles (CSN1S2 A, B, B1, B2, C, D, E, and F) responsible for seven different possible translations of the buffalo αs2-casein. Haplotype data analysis suggests an evolutionary pathway of buffalo CSN1S2 gene consistent with our proposal that the published allele CSN1S2 A is the ancestral αs2-CN form, and the B2 probably arises from interallelic recombination (single crossing) between the alleles D and B (or B1). The allele CSN1S2 C is of new identification, while CSN1S2 B, B1, and B2 are deleted alleles because all are characterized by the mutation g.7539G>C. Two SNPs (g.7539G>C and g.14067A>G) were genotyped in 747 Italian buffaloes, and major alleles had a relative frequency of 0.83 and 0.51, respectively. An association study between these SNPs and milk traits including fatty acid composition was carried out. The SNP g.14067A>G showed a significant association (P < 0.05) on the content of palmitic acid in buffalo milk, thus suggesting its use in marker-assisted selection programs aiming for the improvement of buffalo milk fatty acid composition.

Association of GH polymorphisms with growth traits in buffaloes

Members of the somatotrophic axis, especially GH and IGF1, are essential for growth. The association between GH polymorphisms and growth traits was numerously studied in cattle; however, no data are available for such association studies in buffalo. Therefore, this study was conducted to screen for polymorphisms in the GH gene and to study their putative association with growth traits in 200 Egyptian buffaloes. Polymerase chain reaction single-strand conformation polymorphism and sequencing were applied to look for polymorphisms in 3 loci spanning all exons and introns of buffalo GH. The C (MspI+) >T (MspI-) SNP in intron3, which is well known in cattle, was not detected in the examined buffaloes. However, 2 missense mutations were detected in exon5: one previously detected p.Leu153Val SNP, with very low frequencies for the mutant (Val) allele and one novel p.Asn174His SNP. At weaning age, the p.Leu153Val SNP was significantly associated with weaning body weight and gain with the positive effect of the wild allele (Leu) and higher GH serum, mRNA, and protein levels in animals with Leu/Asn and Leu/His haplotypes. At yearling age, the 2 SNPs associated with yearling weight and gain with positive effect for the mutant (Val and His) alleles with increased GH, and IGF1 serum, mRNA, and protein and GHR mRNA and protein levels in animals with Val/Asn haplotype. Therefore, the selection of Egyptian buffaloes with the Val/Asn haplotype could improve the growth traits of Egyptian buffaloes at yearling age which is the target age for perfect growing.

Sequencing of lipoprotein lipase gene in the Mediterranean river buffalo identified novel variants affecting gene expression

Lipoprotein lipase (LPL) is a key enzyme for lipid metabolism, playing a fundamental role in the composition of fat in adipose tissue and milk. The LPL gene has been seldom investigated in dairy ruminants and barely studied in river buffalo (Bubalus bubalis). The aim of this work was to explore the genetic diversity of LPL and its promoter and to identify functional mutations, using a combined approach based on sequencing, dual-color electrophoretic mobility shift assay, and quantitative PCR. Thirteen consensus sequences for transcription factors were found in the promoter. Eleven SNP were detected, and the attention was focused on the SNP with potential functional effects: g.-446A>G, because the presence of G created a consensus motif for the transcription factor Sp1, and g.107A>G, which was the only exonic SNP. We developed PCR-RFLP methods for genotyping the 2 SNP and calculated the allele frequencies. A strong linkage disequilibrium (D' = 1; r² = 0.903) was found between the 2 SNP. The dual-color electrophoretic mobility shift assay demonstrated that only genotype g.-446GG allowed the binding of the Sp1 transcription factor, resulting in overexpression of the gene (~2.5 fold), as confirmed by the quantitative PCR results. Haploinsufficiency is proposed as a regulation mechanism. This study adds further knowledge on the structure of the LPL gene and its expression in river buffalo, with potential effects on milk qualitative and quantitative production.

Genetic characterization of the oxytocin-neurophysin I gene (OXT) and its regulatory regions analysis in domestic Old and New World camelids

Oxytocin is a neurohypophysial peptide linked to a wide range of biological functions, including milk ejection, temperament and reproduction. Aims of the present study were a) the characterization of the OXT (Oxytocin-neurophysin I) gene and its regulatory regions in Old and New world camelids; b) the investigation of the genetic diversity and the discovery of markers potentially affecting the gene regulation. On average, the gene extends over 814 bp, ranging between 825 bp in dromedary, 811 bp in Bactrian and 810 bp in llama and alpaca. Such difference in size is due to a duplication event of 21 bp in dromedary. The main regulatory elements, including the composite hormone response elements (CHREs), were identified in the promoter, whereas the presence of mature microRNAs binding sequences in the 3'UTR improves the knowledge on the factors putatively involved in the OXT gene regulation, although their specific biological effect needs to be still elucidated. The sequencing of genomic DNA allowed the identification of 17 intraspecific polymorphisms and 69 nucleotide differences among the four species. One of these (MF464535:g.622C>G) is responsible, in alpaca, for the loss of a consensus sequence for the transcription factor SP1. Furthermore, the same SNP falls within a CpG island and it creates a new methylation site, thus opening future possibilities of investigation to verify the influence of the novel allelic variant in the OXT gene regulation. A PCR-RFLP method was setup for the genotyping and the frequency of the allele C was 0.93 in a population of 71 alpacas. The obtained data clarify the structure of OXT gene in domestic camelids and add knowledge to the genetic variability of a genomic region, which has received little investigation so far. These findings open the opportunity for new investigations, including association studies with productive and reproductive traits.

Remarkable genetic diversity detected at river buffalo prolactin receptor (PRLR) gene and association studies with milk fatty acid composition

Prolactin is an anterior pituitary peptide hormone involved in many different endocrine activities and is essential for reproductive performance. This action is mediated by its receptor, the prolactin receptor, encoded by the PRLR gene. In this study, we sequenced and characterized the Mediterranean river buffalo PRLR gene (from exon 3 to 10), and we found remarkable genetic diversity. In particular, we found 24 intronic polymorphisms and 13 exonic SNPs, seven of which were non-synonymous. Furthermore, the polymorphisms identified in the 3'-UTR were investigated to establish their possible influence on microRNA binding sites. Considering all the amino acid changes and the observed allelic combinations, it is possible to deduce at least six different translations of the buffalo prolactin receptor and, consequently, the presence at the PRLR gene of at least six alleles. Furthermore, we identified a deletion of a CACTACC heptamer between nucleotides 1102 and 1103 of exon 10 (3'-UTR), and we developed an allele-specific PCR to identify the carriers of this genetic marker. Finally, the SNP g.11188A>G, detected in exon 10 and responsible for the amino acid replacement p.His328Arg, was genotyped in 308 Italian Mediterranean river buffaloes, and an association study with milk fat traits was carried out. The statistical analysis showed a tendency that approached significance for the AA genotype with higher contents of odd branched-chain fatty acids. Thus, our results suggest that the PRLR gene is a good candidate for gene association studies with qualitative traits related to buffalo milk production.