Comparative analysis of gene sequence of goat CSN1S1 F and N alleles and characterization of CSN1S1 transcript variants in mammary gland

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.

Sequencing and Characterization of αs2-Casein Gene (CSN1S2) in the Old-World Camels Have Proven Genetic Variations Useful for the Understanding of Species Diversification

The CSN1S2 gene encodes αs2-casein, the third most abundant protein in camel milk. Despite its importance in foals, human nutrition, and dairy processing, the CSN1S2 gene in camels has received little attention. This study presents the first complete characterization of the CSN1S2 gene sequence in Old-World camels (Camelus bactrianus and Camelus dromedarius). Additionally, the gene promoter, consisting of 752 bp upstream of exon 1, was analyzed. The entire gene comprises 17 exons, ranging in length from 24 bp (exons 4, 8, 11, and 13) to 280 bp (exon 17). Interesting was the identification of the exon 12 in both species. The promoter analysis revealed 24 putative binding sites in the Bactrian camel and 22 in dromedary camel. Most of these sites were typical elements associated with milk protein, such as C/EBP-α, C/EBP-β, Oct-1, and AP1. The SNP discovery showed relatively high genetic diversity compared to other camel casein genes (CSN1S1, CSN2, and CSN3), with a total of 34 polymorphic sites across the two species. Particularly noteworthy is the transition g.311G>A in the CSN1S2 promoter, creating a new putative consensus binding site for a C/EBP-β in the Bactrian camel. At the exon level, two novel variants were found. One was detected in exon 6 of the Bactrian camel (g.3639C>G), resulting in an amino acid replacement, p.36Ile>Met. The second variant was found in noncoding exon 17 of dromedary CSN1S2 (g.1511G>T). Although this mutation occurs in the 3'-UnTranslated Region, it represents the first example of exonic polymorphism in the CSN1S2 for this species. This SNP also affects the binding sites of different microRNAs, including the seed sequence of the miRNA 4662a-3p, highlighting its role as a regulatory factor for CSN1S2 gene. A PCR-RFLP was set up for genotyping a dromedary Tunisian population (n = 157), and the minor allele frequency was found to be 0.27 for the G allele, indicating a potential yield improvement margin. The interspersed elements (INEs) analysis revealed 10 INEs covering 7.34% and 8.14% of the CSN1S2 sequence in the Bactrian and dromedary camels, respectively. Furthermore, six elements (A, B, F, H, I, and L) are shared among cattle and camels and are partially found in other ruminants, suggesting a common ancestral origin of these retrotransposons. Conversely, elements C, D, E, and G are specific to camels.

A Genotyping Method for Detecting Foreign Buffalo Material in Mozzarella di Bufala Campana Cheese Using Allele-Specific- and Single-Tube Heminested-Polymerase Chain Reaction

Mozzarella di Bufala Campana (MdBC) cheese is a Protected Designation of Origin (PDO) product that is important for the economy and cultural heritage of the Campania region. Food fraud can undermine consumers' trust in this dairy product and harm the livelihood of local producers. The current methods for detecting adulteration in MdBC cheese due to the use of buffalo material from foreign countries could exhibit limitations associated with the required use of expensive equipment, time-consuming procedures, and specialized personnel. To address these limits here, we propose a rapid, reliable, and cost-effective genotyping method that can detect foreign buffalo milk in a counterpart from the PDO area and in MdBC cheese, ensuring the quality and authenticity of the latter dairy product. This method is based on dedicated allele-specific and single-tube heminested polymerase chain reaction procedures. By using allele-specific primers that are designed to detect the nucleotide g.472G>C mutation of the CSN1S1B^bt allele, we distinguished an amplicon of 330 bp in the amplification product of DNA when extracted from milk and cheese, which is specific to the material originating from foreign countries. By spiking foreign milk samples with known amounts of the counterpart from the PDO area, the sensitivity of this assay was determined to be 0.01% v/v foreign to PDO milk. Based on a rough estimate of its simplicity, reliability, and cost, this method could be a valuable tool for identifying adulterated buffalo PDO dairy products.

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.

Molecular Characterization and Phylogenetic Analysis of Casein Gene Family in Camelus ferus

Camel milk is known for its exceptional medical uses. It has been used since ancient times to treat infant diarrhea, hepatitis, insulin-dependent diabetes (IDDM), lactose intolerance, alcohol-induced liver damage, allergies, and autism. It has the power to treat several diseases, with cancer being the most significant. This study investigated the evolutionary relationship, physiochemical characteristics, and comparative genomic analysis of the casein gene family (CSN1S1, CSN2, CSN1S2, and CSN3) in Camelus ferus. Molecular phylogenetics showing the camelid species clustered casein nucleotide sequences into four groups: CSN1S1, CSN2, CSN1S2, and CSN3. The casein proteins from camels were evaluated and found to be unstable, thermostable, and hydrophilic. CSN1S2, CSN2, and CSN3 were acidic, but CSN1S1 was basic. CSN1S1 showed positive selection for one amino acid (Q), CSN1S2 and CSN2 for three (T, K, Q), and CSN3 showed no positive selection. We also compared high-milk-output species such as cattle (Bos Tarus) and low-milk-yield species such as sheep (Ovies Aries) with camels (Camel ferus) and discovered that YY1 sites are more frequent in sheep than in camels and very low in cattle. We concluded that the ratio of YY1 sites in these species may affect milk production.

A New AS-PCR Method to Detect CSN201 Allele, Genotyping at Ca-Sensitive Caseins Loci and Milk Traits Association Studies in Autochthonous Lazio Goats

Calcium-sensitive caseins are the main protein component of milk. In the goat, they are encoded by three genes (CSN1S1, CSN2, and CSN1S2) located on chromosome 6. A high number of alleles has been discovered for these genes in the goat species, responsible for changes in the milk’s qualitative and quantitative characteristics. This study aimed to develop an Allele-Specific PCR (AS-PCR), which allowed us to unequivocally detect goat carriers of the CSN201 allele. Subsequently, the calcium-sensitive casein loci genotype was investigated in three native goat breeds of the Lazio Region (Bianca Monticellana, Capestrina, and Ciociara Grigia). No individuals were carriers of the CSN1S101, CSN1S1E, CSN201, CSN1S2D, and CSN1S20 alleles, while a high frequency of the alleles CSN1S1F and CSN1S1A*,B* was observed. Association analyses between the different genotypes at the CSN1S1 locus and some milk traits, namely the fat and protein yielded and the fat, protein, solids-not-fat, and casein percentages without an effect on the milk yield, were observed.

Review: Genetic and protein variants of milk caseins in goats

The milk casein genes in goats, are highly polymorphic genes with numerous synonymous and non-synonymous mutations. So far, 20 protein variants have been reported in goats for alpha-S1-casein, eight for beta-casein, 14 for alpha-S2-casein, and 24 for kappa-casein. This review provides a comprehensive overview on identified milk casein protein variants in goat and non-coding DNA sequence variants with some affecting the expression of the casein genes. The high frequency of some casein protein variants in different goat breeds and geographical regions might reflect specific breeding goals with respect to milk processing characteristics, properties for human nutrition and health, or adaptation to the environment. Because protein names, alongside the discovery of protein variants, go through a historical process, we linked old protein names with new ones that reveal more genetic variability. The haplotypes across the cluster of the four genetically linked casein genes are recommended as a valuable genetic tool for discrimination between breeds, managing genetic diversity within and between goat populations, and breeding strategies. The enormous variation in the casein proteins and genes is crucial for producing milk and dairy products with different properties for human health and nutrition, and for genetic improvement depending on local breeding goals.

Comparative Analysis of the Polymorphism of the Casein Genes in Camels Bred in Kazakhstan

Caseins play an important role in determining the technological properties and quantitative characteristics of camel milk. To date, only a few studies on the genetic polymorphism of casein genes have been reported in the camel populations of Kazakhstan. Therefore, this work aimed to identify the genetic polymorphism level of casein genes among camel populations of the Almaty region of Kazakhstan. The PCR-RFLP method was used for this purpose and the following genotypes were revealed as a result: CSN3 gene—CC, CT, TT, where the T allele predominated in all populations, with a frequency of 0.60; CSN2 gene—AA, AG, GG, with the predomination of A allele (0.64); and CSN1S1 gene—GG and GT, with the predomination of G allele (0.94). Statistical analysis was carried out using the POPGENE and GenAlEx software. The χ2 values were equal to 12.1 (CSN3), 8.6 (CSN2), and 14.5 (CSN1S1). As a result, three out of 53 animals were designated as the “core” of the population—animals with the desired genotypes: CC genotype for the CSN3 gene and AA genotype for the CSN2 gene. Such animals can be selected for further use with an increase in the number of livestock with high productivity rates.

Genetic Diversity in Casein Gene Cluster in a Dromedary Camel (C. dromedarius) Population from the United Arab Emirates

Genetic polymorphisms, causing variation in casein genes (CSN1S1, CSN1S2, CSN2, and CSN3), have been extensively studied in goats and cows, but there are only few studies reported in camels. Therefore, we aimed to identify alleles with functional roles in the United Arab Emirates dromedary camel (Camelus dromedarius) population to complement previous studies conducted on the same species. Using targeted next-generation sequencing, we sequenced all genes in the casein gene cluster in 93 female camels to identify and characterize novel gene variants. Most variants were found in noncoding introns and upstream sequences, but a few variants showed the possibility of functional impact. CSN2 was found to be most polymorphic, with total 91 different variants, followed by CSN1S1, CSN3 and CSN1S2. CSN1S1, CSN1S2 and CSN2 each had at least two variants while CSN3 had only one functional allele. In future research, the functional impact of these variants should be investigated further.

Capture Sequencing to Explore and Map Rare Casein Variants in Goats

Genetic variations in the four casein genes CSN1S1, CSN2, CSN1S2, and CSN3 have obtained substantial attention since they affect the milk protein yield, milk composition, cheese processing properties, and digestibility as well as tolerance in human nutrition. Furthermore, milk protein variants are used for breed characterization, biodiversity, and phylogenetic studies. The current study aimed at the identification of casein protein variants in five domestic goat breeds from Sudan (Nubian, Desert, Nilotic, Taggar, and Saanen) and three wild goat species [Capra aegagrus aegagrus (Bezoar ibex), Capra nubiana (Nubian ibex), and Capra ibex (Alpine ibex)]. High-density capture sequencing of 33 goats identified in total 22 non-synonymous and 13 synonymous single nucleotide polymorphisms (SNPs), of which nine non-synonymous and seven synonymous SNPs are new. In the CSN1S1 gene, the new non-synonymous SNP ss7213522403 segregated in Alpine ibex. In the CSN2 gene, the new non-synonymous SNPs ss7213522526, ss7213522558, and ss7213522487 were found exclusively in Nubian and Alpine ibex. In the CSN1S2 gene, the new non-synonymous SNPs ss7213522477, ss7213522549, and ss7213522575 were found in Nubian ibex only. In the CSN3 gene, the non-synonymous SNPs ss7213522604 and ss7213522610 were found in Alpine ibex. The identified DNA sequence variants led to the detection of nine new casein protein variants. New variants were detected for alpha S1 casein in Saanen goats (CSN1S1^∗C1), Bezoar ibex (CSN1S1^∗J), and Alpine ibex (CSN1S1^∗K), for beta and kappa caseins in Alpine ibex (CSN2^∗F and CSN3^∗X), and for alpha S2 casein in all domesticated and wild goats (CSN1S2^∗H), in Nubian and Desert goats (CSN1S2^∗I), or in Nubian ibex only (CSN1S2^∗J and CSN1S2^∗K). The results show that most novel SNPs and protein variants occur in the critically endangered Nubian ibex. This highlights the importance of the preservation of this endangered breed. Furthermore, we suggest validating and further characterizing the new casein protein variants.

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.

Casein Gene Cluster in Camelids: Comparative Genome Analysis and New Findings on Haplotype Variability and Physical Mapping

The structure of casein genes has been fully understood in llamas, whereas in other camelids, this information is still incomplete. In fact, structure and polymorphisms have been identified in three (CSN1S1, αs1-CN; CSN2, β-CN; CSN3, κ-CN) out of four casein genes, whereas controversial information is available for the CSN1S2 (αs2-CN) in terms of structure and genetic diversity. Data from the genome analysis, whose assembly is available for feral camel, Bactrian, dromedary, and alpaca, can contribute to a better knowledge. However, a majority of the scaffolds available in GenBank are still unplaced, and the comparative annotation is often inaccurate or lacking.Therefore, the aims of this study are 1) to perform a comparative genome analysis and synthesize the literature data on camelids casein cluster; 2) to analyze the casein variability in two dromedary populations (Sudanese and Nigerian) using polymorphisms at CSN1S1 (c.150G > T), CSN2 (g.2126A > G), and CSN3 (g.1029T > C); and 3) to physically map the casein cluster in alpaca. Exon structures, gene and intergenic distances, large insertion/deletion events, SNPs, and microsatellites were annotated. In all camelids, the CSN1S2 consists of 17 exons, confirming the structure of llama CSN1S2 gene. The comparative analysis of the complete casein cluster (∼190kb) shows 12,818 polymorphisms. The most polymorphic gene is the CSN1S1 (99 SNPs in Bactrian vs. 248 in dromedary vs. 626 in alpaca). The less polymorphic is the CSN3 in the Bactrian (22 SNPs) and alpaca (301 SNPs), whereas it is the CSN1S2 in dromedary (79 SNPs). In the two investigated dromedary populations, the allele frequencies for the three markers are slightly different: the allele C at CSN1S1 is very rare in Nigerian (0.054) and Sudanese dromedaries (0.094), whereas the frequency of the allele G at CSN2 is almost inverted. Haplotype analysis evidenced GAC as the most frequent (0.288) and TGC as the rarest (0.005). The analysis of R-banding metaphases hybridized with specific probes mapped the casein genes on chromosome 2q21 in alpaca. These data deepen the information on the structure of the casein cluster in camelids and add knowledge on the cytogenetic map and haplotype variability.

Alternative splicing events expand molecular diversity of camel CSN1S2 increasing its ability to generate potentially bioactive peptides

In a previous study on camel milk from Kazakhstan, we reported the occurrence of two unknown proteins (UP1 and UP2) with different levels of phosphorylation. Here we show that UP1 and UP2 are isoforms of camel αs2-CN (αs2-CNsv1 and αs2-CNsv2, respectively) arising from alternative splicing events. First described as a 178 amino-acids long protein carrying eight phosphate groups, the major camel αs2-CN isoform (called here αs2-CN) has a molecular mass of 21,906 Da. αs2-CNsv1, a rather frequent (35%) isoform displaying a higher molecular mass (+1,033 Da), is present at four phosphorylation levels (8P to 11P). Using cDNA-sequencing, αs2-CNsv1 was shown to be a variant arising from the splicing-in of an in-frame 27-nucleotide sequence encoding the nonapeptide ENSKKTVDM, for which the presence at the genome level was confirmed. αs2-CNsv2, which appeared to be present at 8P to 12P, was shown to include an additional decapeptide (VKAYQIIPNL) revealed by LC-MS/MS, encoded by a 3'-extension of exon 16. Since milk proteins represent a reservoir of biologically active peptides, the molecular diversity generated by differential splicing might increase its content. To evaluate this possibility, we searched for bioactive peptides encrypted in the different camel αs2-CN isoforms, using an in silico approach. Several peptides, putatively released from the C-terminal part of camel αs2-CN isoforms after in silico digestion by proteases from the digestive tract, were predicted to display anti-bacterial and antihypertensive activities.

The single nucleotide polymorphism g.133A>C in the stearoyl CoA desaturase gene (SCD) promoter affects gene expression and quali-quantitative properties of river buffalo milk

The stearoyl-CoA desaturase (SCD) gene has been investigated in depth in ruminants because of its effect on milk fat composition. In river buffalo, the single nucleotide polymorphism (SNP) g.133A>C in the gene promoter has been associated with milk quality and yield. However, the biological reason for such effects remains unexplored. In this study, we combined mRNA profile analysis, an electromobility shift assay, and quantitative PCR to elucidate the role of this SNP on gene transcription and its effects on milk fat traits. A preliminary genotyping of g.133A>C was carried out on a group of 303 river buffaloes to choose individuals for the downstream applications. Analysis of allele frequencies showed an increase in the minor allele C (0.25) compared with previous findings (0.16). Six animals (2 for each genotype) were chosen for cloning and 216 positive cDNA recombinant clones for SCD (72 per genotype) were analyzed by PCR. All clones showed the same length on agarose gel; therefore, random clones were chosen for sequencing. No qualitative differences were found and all gene transcripts assembled correctly. An electrophoretic mobility shift assay was performed to evaluate the binding of the transcription factor Sp1 to DNA sequences including g.133A>C. Genotype CC showed a higher binding (mean ± standard error of the mean) than genotype AA in 2 different conditions [Enzo buffer (EB), Enzo Life Science Inc., Farmingdale, NY: 201.77 ± 4.06 vs. 141.65 ± 3.77 band intensity values and Poletto buffer (PB): 95.90 ± 1.15 vs. 67.30 ± 2.14 band intensity values]. The subsequent quantitative PCR confirmed the upregulation of the CC genotype compared with the AA and AC genotypes. The association study with milk fat traits revealed a favorable effect of allele C. The heterozygous genotype had the highest values for monounsaturated fatty acids, oleic acid (C18:1 cis-9), polyunsaturated fatty acids, and odd- and branched-chain fatty acids, and the lowest values for saturated fatty acids and atherogenic and thrombogenic indices; the heterozygous genotype differed significantly from the AA genotype. The AC genotype has previously been associated with higher milk yield. Therefore, the g.133A>C SNP is a marker with dual effects and is an interesting candidate for assisted selection programs in river buffalo. These data clarified the biological role of the SNP g.133A>C in the SCD promoter and how it affects gene function, providing important knowledge on the genetic background of lipid metabolism, including the future possibility of selecting alleles with quantitatively or qualitatively favorable effects.

A novel indel within goat casein alpha S1 gene is significantly associated with litter size

The ruminant casein gene family (CSNs, link as CSN1S1-CSN2-CSN1S2-CSN3) is characterized by diverse variations and has been extensively studied for ruminant milk traits, however, studies on insertion/deletion (indel) mutations within this gene family and their effects on prolificacy are extremely limited. This study aimed to detect possible novel indels within CSNs in an indigenous Chinese goat breed-Shaanbei White Cashmere goat (SBWC, n = 3047) and four other Chinese goat breeds (n = 1136) with varied litter size rates (105%-283%) and different estrus types (seasonal vs. perennial), as well as exploring the association between these potential indels and litter size. Only one novel 11-bp indel within the CSN1S1 gene was found. The association analyses uncovered that this novel indel was related to the first-birth litter size of SBWC population (n = 2690) (P < 1.0 E-8). Individuals with the II genotype (n = 676) had the best litter size when compared with those ID genotype (n = 1098) and DD genotype (n = 916) individuals. Animals with the II genotype were found to have higher relative expression level of CSN1S1 gene in the ovary (P < 0.01). Besides, Chi-square tests for different litter size and estrous cycle breeds showed that perennial-estrus breeds and multi-kids breeds had higher "I" allelic frequencies and "II" genotypic frequencies. These findings suggest the 11-bp indel within the CSN1S1 gene is significantly associated with reproduction traits and can be an effective molecular marker for litter size of goat breeding.

Transcript analysis at DGAT1 reveals different mRNA profiles in river buffaloes with extreme phenotypes for milk fat

Buffalo DGAT1 (diacylglycerol O-acyltransferase 1) was mainly investigated for the characterization of the gene itself and for the identification of the K232A polymorphism, similar to what has been accomplished in cattle, although no information has been reported so far at the mRNA level. The importance of DGAT1 for lipid metabolism led us to investigate the transcript profiles of lactating buffaloes characterized as high (9.13 ± 0.23) and low (7.94 ± 0.29) for milk fat percentage, and to explore the genetic diversity at the RNA and DNA level. A total of 336 positive clones for the DGAT1 cDNA were analyzed by PCR and chosen for sequencing according to the differences in length. The clone assembling revealed a very complex mRNA pattern with a total of 21 transcripts differently represented in the 2 groups of animals. Apart from the correct transcript (17 exons long), the skipping of exon 12 is the most significant in terms of distribution of clones with 11.6% difference between the 2 groups, whereas a totally different mRNA profile was found in approximately 12% of clones. The sequencing of genomic DNA allowed the identification of 10 polymorphic sites at the intron level, which clarify, at least partially, the genetic events behind the production of complex mRNA. Genetic diversity was found also at the exon level. The single nucleotide polymorphism c.1053C>T represents the first example of polymorphism in a coding region for the DGAT1 in the Italian Mediterranean breed. To establish whether this polymorphism is present in other buffalo breeds, a quick method based on PCR-RFLP was set up for allelic discrimination in the Italian Mediterranean and the Romanian Murrah (200 animals in total). The alleles were equally represented in the overall population, whereas the analysis of the 2 breeds showed different frequencies, likely indicating diverse genetic structure of the 2 breeds. The T allele might be considered as the ancestral condition of the DGAT1 gene, being present in the great part of the sequenced species. These data add knowledge at the transcript and genetic levels for the buffalo DGAT1 and open the opportunity for further investigation of other genes involved in milk fat metabolism for the river buffalo, including the future possibility of selecting alleles with quantitative or qualitative favorable effects (or both).

Characterization of the liver X receptor-dependent regulatory mechanism of goat stearoyl-coenzyme A desaturase 1 gene by linoleic acid

Stearoyl-coenzyme A desaturase 1 (SCD1) is a key enzyme in the biosynthesis of palmitoleic and oleic acid. Although the transcriptional regulatory mechanism of SCD1 via polyunsaturated fatty acids (PUFA) has been extensively explored in nonruminants, the existence of such mechanism in ruminant mammary gland remains unknown. In this study, we used goat genomic DNA to clone and sequence a 1,713-bp fragment of the SCD1 5' flanking region. Deletion assays revealed a core region of the promoter located between -415 and -109 bp upstream of the transcription start site, and contained the highly conserved PUFA response region. An intact PUFA response region was required for the basal transcriptional activity of SCD1. Linoleic acid reduced endogenous expression of SCD1 and sterol regulatory element binding factor-1 (SREBF1) in goat mammary epithelial cells. Further analysis indicated that both the sterol response element (SRE) and the nuclear factor Y (NF-Y) binding site in the SCD1 promoter were responsible for the inhibition effect by linoleic acid, whereas the effect was abrogated once NF-Y was deleted. In addition, SRE and NF-Y were partly responsible for the transcriptional activation induced via the liver X receptor agonist T 4506585 (Sigma-Aldrich, St. Louis, MO). When goat mammary epithelial cells were cultured with linoleic acid, addition of T 4506585 markedly increased SCD1 transcription in controls, but had no effect on cells with a deleted SRE promoter. These results demonstrated that linoleic acid can regulate SCD1 expression at the transcriptional level through SRE and NF-Y in a liver X receptor-dependent fashion in the goat mammary gland.

Lama glama αS1-casein: Identification of new polymorphisms in the CSN1S1 gene

South American camelids have been poorly genetically investigated and little information is available in llamas (Lama glama) regarding the diversity of the caseins at the protein and gene level. Exon skipping and duplication events previously reported in the α_S1-casein gene (CSN1S1) led us to investigate the genetic variability at this locus. Seventy-two positive clones for the α_S1-casein transcripts were analyzed and randomly sequenced. The comparative analysis of the sequences revealed 2 transitions, c.366A>G and c.690T>C, at the 10th nucleotide of exon 12 and 94 bp of exon 19, respectively. These SNP are responsible for 2 amino acid changes, Ile→Val in position 86 and Tyr→His in position 194 of the mature protein. Both polymorphisms clarify the genetic events behind the protein variants A and B. This result was confirmed by isoelectric focusing analysis of llama milk samples. Quick methods based on PCR-RFLP and allele-specific PCR were set up for allelic discrimination in a population of 128 animals. Based on genotyping results, 4 haplotypes were observed and the estimated frequencies indicated B as the most common haplotype (0.629) in the investigated population. These data add knowledge to the genetic variability of a species little investigated, and open opportunity for new investigation in the field of milk protein for South American camelids, including the possibility, in the future, to select alleles with favorable characteristics.

Thyroid hormone responsive (THRSP) promotes the synthesis of medium-chain fatty acids in goat mammary epithelial cells

In nonruminants, thyroid hormone responsive (THRSP) is a crucial protein for cellular de novo lipogenesis. However, the role of THRSP in regulating the synthesis of milk fatty acid composition in goat mammary gland remains unknown. In the present study, we compared gene expression of THRSP among different goat tissues. Results revealed that THRSP had the highest expression in subcutaneous fat, and expression was higher during lactation compared with the dry period. Overexpression of THRSP upregulated the expression of fatty acid synthase (FASN), stearoyl-coenzyme A desaturase 1 (SCD1), diacylglycerol acyltransferase 2 (DGAT2), and glycerol-3-phosphate acyltransferase (GPAM) in goat mammary epithelial cells. In contrast, overexpression of THRSP led to downregulation of thrombospondin receptor (CD36) and had no effect on the expression of acetyl-coenzyme A carboxylase α (ACACA) and sterol regulatory element binding transcription factor1 (SREBF1). In addition, overexpressing THRSP in vitro resulted in a significant increase in triacylglycerol (TAG) concentration and the concentrations of C12:0 and C14:0. Taken together, these results highlight an important role of THRSP in regulating lipogenesis in goat mammary epithelial cells.

Characterization of equine CSN1S2 variants considering genetics, transcriptomics, and proteomics

Currently, research interest is increasing in horse milk composition and its effect on human health. Despite previously published studies describing the presence of intra- and interbreed variability of equine milk components, no investigations have focused on the genetic background of this variation. Among horse caseins and the genes encoding them, least is known about the structure and expression of the α-S2 casein gene, CSN1S2. Herein, based on direct sequencing of the equine CSN1S2 coding sequence, we describe the presence of 51-bp insertion-deletion (in/del) polymorphism, which significantly changes the protein sequence (lack or presence of 17-amino acid serine-rich peptide). Bioinformatic analysis revealed that the observed in/del polymorphism spanned exactly 2 exons; therefore, we hypothesized that we were observing different CSN1S2 splicing isoforms. However, further investigation indicated that the detected sequence variation was caused by a large (1.3-kb) deletion in the genomic DNA. We found that the polymorphic forms (A, longer; B, shorter; KP658381 and KP658382 GenBank records, respectively) were unevenly distributed among different horse breeds (the highest frequency of variant B was observed in coldblood horses and Haflingers). We propose that the analyzed polymorphism is associated with CSN1S2 expression level (the highest expression was recorded for individuals carrying the BB genotype), which was much more pronounced for milk CSN1S2 protein content than for relative transcript abundance (measured in milk somatic cells). Our results provide insight into the equine CSN1S2 structure and lay a foundation for further functional analyses regarding, for example, allergenicity or physiochemical properties of the observed CSN1S2 variants.

Molecular Characterization of the Llamas (Lama glama) Casein Cluster Genes Transcripts (CSN1S1, CSN2, CSN1S2, CSN3) and Regulatory Regions

In the present paper, we report for the first time the characterization of llama (Lama glama) caseins at transcriptomic and genetic level. A total of 288 casein clones transcripts were analysed from two lactating llamas. The most represented mRNA populations were those correctly assembled (85.07%) and they encoded for mature proteins of 215, 217, 187 and 162 amino acids respectively for the CSN1S1, CSN2, CSN1S2 and CSN3 genes. The exonic subdivision evidenced a structure made of 21, 9, 17 and 6 exons for the αs1-, β-, αs2- and κ-casein genes respectively. Exon skipping and duplication events were evidenced. Two variants A and B were identified in the αs1-casein gene as result of the alternative out-splicing of the exon 18. An additional exon coding for a novel esapeptide was found to be cryptic in the κ-casein gene, whereas one extra exon was found in the αs2-casein gene by the comparison with the Camelus dromedaries sequence. A total of 28 putative phosphorylated motifs highlighted a complex heterogeneity and a potential variable degree of post-translational modifications. Ninety-six polymorphic sites were found through the comparison of the lama casein cDNAs with the homologous camel sequences, whereas the first description and characterization of the 5'- and 3'-regulatory regions allowed to identify the main putative consensus sequences involved in the casein genes expression, thus opening the way to new investigations -so far- never achieved in this species.

Mediterranean river buffalo CSN1S1 gene: search for polymorphisms and association studies

The aim of the present study was to investigate the variability at CSN1S1 locus of the Italian Mediterranean river buffalo and to study possible allele effects on milk yield and its composition. Effects of parity, calving season and month of production were also evaluated. Three single-nucleotide polymorphisms were detected. The first mutation, located at position 89 of the 17th exon (c.628C>T), is responsible for the amino acid change p.Ser178 (B allele)/Leu178 (A allele). The other two polymorphisms, detected at the positions 144 (c.882G>A) and 239 (c.977A>G) of 19th exon, respectively, are silent (3ʹ UTR, untranslated region). Associations between the CSN1S1 genotypes and milk production traits were investigated using 4122 test day records of 503 lactations from 175 buffalo cows. Milk yield, fat and protein percentages were analysed using a mixed linear model. A significant association between the c.628C>T SNP and the protein percentage was found. In particular, the CC genotype showed an average value ~0.04% higher than the CT and TT genotypes. The allele substitution effect of cytosine into thymine was –0.014, with a quite low (0.3%) protein percentage contribution to total phenotypic variance. A large dominance effect was detected. Characterisation of the CSN1S1 transcripts and a method based on MboI amplification created restriction site PCR for a rapid genotyping of c.628C>T are provided.

Lipoprotein lipase, tissue expression and effects on genes related to fatty acid synthesis in goat mammary epithelial cells

Lipoprotein lipase (LPL) serves as a central factor in hydrolysis of triacylglycerol and uptake of free fatty acids from the plasma. However, there are limited data concerning the action of LPL on the regulation of milk fat synthesis in goat mammary gland. In this investigation, we describe the cloning and sequencing of the LPL gene from Xinong Saanen dairy goat mammary gland, along with a study of its phylogenetic relationships. Sequence analysis showed that goat LPL shares similarities with other species including sheep, bovine, human and mouse. LPL mRNA expression in various tissues determined by RT-qPCR revealed the highest expression in white adipose tissue, with lower expression in heart, lung, spleen, rumen, small intestine, mammary gland, and kidney. Expression was almost undetectable in liver and muscle. The expression profiles of LPL gene in mammary gland at early, peak, mid, late lactation, and the dry period were also measured. Compared with the dry period, LPL mRNA expression was markedly greater at early lactation. However, compared with early lactation, the expression was lower at peak lactation and mid lactation. Despite those differences, LPL mRNA expression was still greater at peak, mid, and late lactation compared with the dry period. Using goat mammary epithelial cells (GMEC), the in vitro knockdown of LPL via shRNA or with Orlistat resulted in a similar degree of down-regulation of LPL (respectively). Furthermore, knockdown of LPL was associated with reduced mRNA expression of SREBF1, FASN, LIPE and PPARG but greater expression of FFAR3. There was no effect on ACACA expression. Orlistat decreased expression of LIPE, FASN, ACACA, and PPARG, and increased FFAR3 and SREBF1 expression. The pattern of LPL expression was similar to the changes in milk fat percentage in lactating goats. Taken together, results suggest that LPL may play a crucial role in fatty acid synthesis.

Cloning and comparative analysis of gene structure in promoter site of alpha-s1 casein gene in Naeinian goat and sheep

The 5' end or alpha-S1 casein promoter has a significant role in milk protein gene expression. The understanding of the translation process of alpha-S1 casein mutants will provide us an opportunity to make the best selection in livestock providing more proteins in milk. Blood samples were taken from three hundred of Naeinian goats and sheep, and DNA extraction was done using modified salting out method. Polymerase chain reactions (PCR) were carried out using a specific primer pairs for amplification a fragment of 1133 bp from part of 5'-UTR and exon 1 of alpha s1 casein gene. The AluI and HinfI restriction enzyme treatment of all samples provided the same homozygous AA genotype in both species. Subsequently, one sample of each species was selected and cloned, and the final sequences were analyzed by BioEdit, CLC genomic, Mega4 and DNASIS MAX software. Several polymorphisms are recognized between Naeinian goat and sheep that are presented on motif sites. In this research, the interested location, including exon I and a part of 5', was analyzed, and genetic element comparisons were done between Naeinian goat and sheep. The number and location of probable binding sites can have a crucial role as a result of antagonistic and synergistic effects on gene regulation activities.

The Application of Genomic Technologies to Investigate the Inheritance of Economically Important Traits in Goats

Goat genomics has evolved at a low pace because of a lack of molecular tools and sufficient investment. Whilst thousands and hundreds of quantitative trait loci (QTL) have been identified in cattle and sheep, respectively, about nine genome scans have been performed in goats dealing with traits as conformation, growth, fiber quality, resistance to nematodes, and milk yield and composition. In contrast, a great effort has been devoted to the characterization of candidate genes and their association with milk, meat, and reproduction phenotypes. In this regard, causal mutations have been identified in the αS1-casein gene that has a strong effect on milk composition and the PIS locus that is linked to intersexuality and polledness. In recent times, the development of massive parallel sequencing technologies has allowed to build a reference genome for goats as well as to monitor the expression of mRNAs and microRNAs in a broad array of tissues and experimental conditions. Besides, the recent design of a 52K SNP chip is expected to have a broad impact in the analysis of the genetic architecture of traits of economic interest as well as in the study of the population structure of goats at a worldwide scale.

Newly identified mutations at the CSN1S1 gene in Ethiopian goats affect casein content and coagulation properties of their milk

Very high casein content and good coagulation properties previously observed in some Ethiopian goat breeds led to investigating the αs1-casein (CSN1S1) gene in these breeds. Selected regions of the CSN1S1 gene were sequenced in 115 goats from 5 breeds (2 indigenous: Arsi-Bale and Somali, 1 exotic: Boer, and 2 crossbreeds: Boer × Arsi-Bale and Boer × Somali). The DNA analysis resulted in 35 new mutations: 3 in exons, 3 in the 5' untranslated region (UTR), and 29 in the introns. The mutations in exons that resulted in an amino acid shift were then picked to evaluate their influence on individual casein content (αs1-, αs2-, β-, and κ-CN), micellar size, and coagulation properties in the milk from the 5 goat breeds. A mutation at nucleotide 10657 (exon 10) involved a transversion: CAG→CCG, resulting in an amino acid exchange Gln77→Pro77. This mutation was associated with the indigenous breeds only. Two new mutations, at nucleotide 6072 (exon 4) and 12165 (exon 12), revealed synonymous transitions: GTC→GTT in Val15 and AGA→AGG in Arg100 of the mature protein. Transitions G→A and C→T at nucleotides 1374 and 1866, respectively, occurred in the 5' UTR, whereas the third mutation involved a transversion T→G at nucleotide location 1592. The goats were grouped into homozygote new (CC), homozygote reference (AA), and heterozygote (CA) based on the nucleotide that involved the transversion. The content of αs1-CN (15.32g/kg) in milk samples of goats homozygous (CC) for this newly identified mutation, Gln77→Pro77 was significantly higher than in milks of heterozygous (CA; 9.05g/kg) and reference (AA; 7.61g/kg) genotype animals. The αs2-, β-, and κ-CN contents showed a similar pattern. Milk from goats with a homozygous new mutation had significantly lower micellar size. Milk from both homozygote and heterozygote new-mutation goats had significantly shorter coagulation rate and stronger gel than the reference genotype. Except the transversion, the sequence corresponded to allele A and presumably derived from it. Therefore, this allele is denoted by A3. All goats from the reference genotype (AA) were homozygous for the allele at nucleotide position 1374 and 1866, whereas all mutations in the 5' UTR existed in a heterozygous form in both heterozygous (CA) and the new mutation (CC) genotype. The newly identified mutation (CC) detected in some of the goat breeds is, therefore, important in selection for genetic improvement and high-quality milk for the emerging goat cheese-producing industries. The finding will also benefit farmers raising these goat breeds due to the increased selling price of goats. Further studies should investigate the effect of this amino acid exchange on the secondary and tertiary structure of the αs1-CN molecule and on the susceptibility of peptide hydrolysis by digestive enzymes.

Identification of an intronic regulatory mutation at the buffalo αS1-casein gene that triggers the skipping of exon 6

The characterization of casein polymorphism is an essential step in order to understand the genetic basis of milk quality in dairy ruminants. In this work, we report the identification of a regulatory mutation at the buffalo αs1-casein (CSN1S1) gene that alters the normal processing of the primary transcript. Sequencing of CSN1S1 cDNA from individuals harbouring this new variant revealed that its most distinctive feature is the loss of exon 6 that encodes eight amino acids between positions 35-42 of mature protein. In an effort to map the causal mutation, we sequenced a genomic region spanning exons 5-7 of the buffalo CSN1S1 gene. This experiment allowed us to establish that exon 6-skipping is produced by a G to C substitution at the first position of intron 6 that inactivates the donor splice site. This mutation can be typed by PCR-RFLP by using either TaaI or Bpu10I diagnostic restriction enzymes, and it has a frequency of 0.18 in Romanian buffaloes. This exon skipping phenomenon is the first one described in buffalo CSN1S1 locus.

PPAR γ Regulates Genes Involved in Triacylglycerol Synthesis and Secretion in Mammary Gland Epithelial Cells of Dairy Goats

To explore the function of PPARγ in the goat mammary gland, we cloned the whole cDNA of the PPARγ gene. Homology alignments revealed that the goat PPARγ gene is conserved among goat, bovine, mouse, and human. Luciferase assays revealed that rosiglitazone enhanced the activity of the PPARγ response element (PPRE) in goat mammary epithelial cells (GMECs). After rosiglitazone (ROSI) treatment of GMECs, there was a significant (P < 0.05) increase in the expression of genes related to triacylglycerol synthesis and secretion: LPL, FASN, ACACA, PLIN3, FABP3, PLIN2, PNPLA2, NR1H3, SREBF1, and SCD. The decreases in expression observed after knockdown of PPARγ relative to the control group (Ad-NC) averaged 65%, 52%, 67%, 55%, 65%, 58%, 85%, 43%, 50%, and 24% for SCD, DGAT1, AGPAT6, SREBF1, ACACA, FASN, FABP3, SCAP, ATGL, and PLIN3, respectively. These results provide direct evidence that PPARγ plays a crucial role in regulating the triacylglycerol synthesis and secretion in goat mammary cells and underscore the functional importance of PPARγ in mammary gland tissue during lactation.

Effect of CSN1S1 gene polymorphism and stage of lactation on milk yield and composition of extensively reared goats

The effect of CSN1S1 genotype and lactation stage on milk yield and composition were investigated in 80 extensively reared goats. Milk yield was recorded in early, mid and late lactation and individual milk samples were collected to determine: fat, protein, lactose and casein content, pH, freezing point, somatic cell count (SCC) and total microbic mesophilic count (TMC). Relative casein composition and amino acid profile were quantified by HPLC. Fatty acid profile was measured by gas-chromatography. Genotype did not affect milk yield, while this trait was significantly affected by lactation stage (P < 0·01). CSN1S1 BB goats produced significantly higher protein and casein percentages (P < 0·05). αs1-casein (CN) was significantly higher in BB and AB goats than AF and BF, showing intermediate values in AA goats (P < 0·01). The protein percentage and the αs1 and αs2-CN fractions were not affected by lactation stage, while the casein content and the β and κ-CN significantly increased throughout lactation (P < 0·01). C4 : 0 and C6 : 0 were not affected by genotype, while C8 : 0 and C10 : 0 were higher in the AA goats than BB; most of the long chain FA were higher in BB than AA goats. MUFA and PUFA increased in late lactation. In addition, BB goats showed higher essential amino acids, resulting in an optimal composition from the nutritional point of view, when compared with AA goats. The increase of MUFA, PUFA, essential and cis-FA in late lactation indicate that the lipid composition of goat's milk, with the progress of lactation, tends to improve its nutritional value.

Genetic polymorphism at the CSN1S1 gene in Girgentana dairy goat breed

The aim of this work was to evaluate the variability of the αs1-casein locus in the endangered Girgentana dairy goat breed in order to define genetic improvement and a conservation program for this breed. The study was performed on 200 dairy goats by means of different PCR protocols. The most frequent alleles were A (0.590) and F (0.290) followed by B (0.065) and N (0.047). CSN1S1 E allele was identified with a very low frequency (0.008). The most common genotype was AF (0.365) followed by AA (0.340). The high frequency of the strong genotypes is associated with the production of milk with high fat and protein content and with optimal technological properties. In Girgentana goat breed, the CSN1S1 genotype information could be utilised in selection strategies for milk protein content and milk yield, in order to select genetic lines for the production of ‘drinking milk’ using weak and null genotypes, and for niche products using strong genotypes.

Additive and dominance effects of the α(s1)-casein locus on milk yield and composition traits in dairy goats

The objective of this study was to evaluate the effects of the CSN1S1 locus polymorphism on 305-d records of milk, fat, protein, lactose and total solids yields, fat, protein, lactose and total solids contents in Mexican dairy goats. A total of 514 lactation records belonging to Alpine (n=60), Saanen (n=105) and Toggenburg (n=74) goats, born from 2003 to 2006 in three herds were used. Discrimination between alleles E, F, N, A* (CSN1S1 A, G, H, I, O1 and O2) and B* (CSN1S1 B1, B2, B3, B4, C and L) were made by amplification of fragments of the gene CSN1S1 and digestion with the restriction endonuclease XmnI. In order to estimate additive and dominance effects, data sets including (1) all genotypes, and (2) only homozygote genotypes, were analysed using linear mixed models. The allele A*, had significant additive effects for protein content (0·21±0·07%; P=0·002) and total solids content (0·66±0·23%; P=0·005) when compared with allele F. An unfavourable additive effect of allele A* on milk yield was found in the Alpine breed (-81·4±40·2; P=0·046) when compared with allele F. Favourable dominance effects were found for some genotypes (P<0·05) for milk yield (A*N and B*N), fat yield (A*N and B*E), protein yield (A*N and B*E), lactose yield (A*N) and total solids yield (A*N). Also, unfavourable dominance effects were found (P<0·05) for protein content (A*B* and A*N) and total solids content (A*B*, A*N, and A*F). Allele A* was the only one with a positive effect for protein content. Significant allele-year interaction effects were also observed. The presence of significant dominance effects, estimated between specific pairs of alleles, challenged the purely additive nature of the genetic effect at the CSN1S1 locus. Implications from use of CSN1S1 effects in goat breeding programmes are presented.

Molecular characterization and genetic variability at κ-casein gene (CSN3) in camels

κ-casein is a glycosilated protein belonging to a family of phosphoproteins (αs1, β, αs2, κ) that represent the major protein component in mammalian milk. κ-casein plays an essential role in the casein micelle stabilization, determining the size and the specific function. In the present paper, we report for the first time the characterization of the nucleotide sequence of the whole κ-casein-encoding gene (CSN3) plus 1045 nucleotides at the 5' flanking region in Camelus dromedarius. The promoter region and the complete cDNA were also provided for the first time in Camelus bactrianus. The gene is spread over 9.3kb and consists of 5 exons varying in length from 33bp (exon 3) to 494bp (exon 4), and 4 introns from 1200bp (intron 3) to 2928bp (intron 2). Highly conserved sequences, located in the 5' flanking region, have been found. The regulatory regions of camels seems to be more related to equids than to other compared species. 17 polymorphic sites have been detected, one of these (g.1029T>C) is responsible for the creation of a new putative consensus sequence for the transcription factor HNF-1. In general, these SNPs are the first reported in camels for casein loci. Finally, seven interspersed repeated elements were also identified at intronic level.

Genetic testing for phenotype-causing variants in sheep and goats

This review gives an overview on ovine and caprine defects/disorders, disease predispositions, production traits and coat colours for which causal gene variants are known. Most phenotypes are inherited autosomal-recessive or dominant and in the majority are caused by single nucleotide substitutions or deletions. Causative sequence variants mainly were identified by sequencing candidate genes in the past, and recently also by whole genome analysis using the ovine 50k SNP chip. While PCR-fragment length polymorphism analyses were developed for the majority of causative sequence variants, other low- to medium-throughput PCR-based methods as PCR-single strand conformation analysis and allele-specific PCR were also established frequently. For processing large sample numbers, high-throughput methods as MALDI-ToF MS or real-time PCR are available for some gene variants. Further progress in development of ovine and caprine genome sequences and SNP chips will be beneficial for the discovery of additional causative variants in these two species.

Casein SNP in Norwegian goats: additive and dominance effects on milk composition and quality

Background

The four casein proteins in goat milk are encoded by four closely linked casein loci (CSN1S1, CSN2, CSN1S2 and CSN3) within 250 kb on caprine chromosome 6. A deletion in exon 12 of CSN1S1, so far reported only in Norwegian goats, has been found at high frequency (0.73). Such a high frequency is difficult to explain because the national breeding goal selects against the variant's effect.

Methods

In this study, 575 goats were genotyped for 38 Single Nucleotide Polymorphisms (SNP) located within the four casein genes. Milk production records of these goats were obtained from the Norwegian Dairy Goat Control. Test-day mixed models with additive and dominance fixed effects of single SNP were fitted in a model including polygenic effects.

Results

Significant additive effects of single SNP within CSN1S1 and CSN3 were found for fat % and protein %, milk yield and milk taste. The allele with the deletion showed additive and dominance effects on protein % and fat %, and overdominance effects on milk quantity (kg) and lactose %. At its current frequency, the observed dominance (overdominance) effects of the deletion allele reduced its substitution effect (and additive genetic variance available for selection) in the population substantially.

Conclusions

The selection pressure of conventional breeding on the allele with the deletion is limited due to the observed dominance (overdominance) effects. Inclusion of molecular information in the national breeding scheme will reduce the frequency of this deletion in the population.

The CSN1S1 N and F alleles identified by PCR-SSCP and their associations with milk yield and composition in Chinese dairy goats

A method was depicted to identify null allele CSN1S1 N and low allele CSN1S1 F of the CSN1S1 gene of goat using PCR-SSCP (polymerase chain reaction-single strand conformation polymorphism). First, primer A was designed to amplify the exon 9 of CSN1S1 gene which produced three genotypes AA, AB, and BB. Among these three genotypes, only AA and AB individuals had a cytosine deletion at exon 9 after DNA sequencing, which cannot be used to identify the N and F alleles. Therefore, primer B was used to amplify intron 14 of CSN1S1 of described AA and AB individuals. Genotypes FF, FN and NN were detected within AA individuals and genotypes FO and NO were detected in the above AB individuals. The frequencies of F and N alleles in 708 samples from Xinong Saanen (XS) and Guanzhong (GZ) dairy goat breeds were 0.1139, 0.0927, and 0.2376, 0.1193, respectively. In 268 XS samples, the individuals with NN genotype contained a significant lower protein content than that of other genotypes (P<0.01). Individuals of FF genotype had significant higher milk yield than that of NO genotype in the first milk lactation of 202 XS individuals (P<0.05). Therefore, the variability at CSN1S1 locus contains enough genetic diversity to be potentially useful in improving the quality and production of milk in Chinese dairy goat breeds.

Characterization of casein gene complex and genetic diversity analysis in Indian goats

Milk protein polymorphism plays an important role in genetic diversity analysis, phylogenetic studies, establishing geographical diversity, conservation decision, and improving breeding goals. Milk protein polymorphism in Indian goat breeds has not been well studied; therefore, an investigation was carried out to analyze the genetic structure of the casein gene and milk protein diversity at six milk protein loci in nine Indian goat breeds/genetic groups from varied agro-climatic zones. Milk protein genotyping was carried out in 1098 individual milk samples by SDS-PAGE at alphaS1-CN (CSN1S1), beta-CN (CSN2), alphaS2-CN (CSN1S2), kappa-CN (CSN3), beta-LG, and alpha-LA loci. Indian goats exhibited alphaS1-casein A allele in higher frequency in the majority of breeds except Ganjam and local goats. The alphaS1-casein A allele frequencies varied from 0.45 to 0.77. A total of 16 casein haplotypes were observed in seven breeds and breed specific haplotypes were observed with respect to geographic region. The average number of alleles was lowest in Ganjam (1.66 +/- 0.81) and highest in Sirohi goats (2.50 +/- 1.05). Expected heterozygosity at six different loci demonstrated genetic diversity and breed fragmentation. Neighbor-Joining tree was built basing on Nei's distance. There was about 16.95% variability due to differences between breeds, indicating a strong subdivision. Principal component analysis was carried out to highlight the relationship among breeds. The variability among goat breeds was contributed by alphaS2-CN, beta-LG and alphaS1-CN. The Indian goats exhibited alphaS1-CN (CSN1S1) A allele in higher frequency in all the breeds indicating the higher casein yield in their milk.