The β-casein in camels: molecular characterization of the CSN2 gene, promoter analysis and genetic variability

Review on camel genetic diversity: ecological and economic perspectives

Camels, known as the "Ship of the Desert," play a vital role in the ecosystems and economies of arid and semi-arid regions. They provide meat, milk, transportation, and other essential services, and their resilience to harsh environments makes them invaluable. Despite their similarities, camel breeds exhibit notable differences in size, color, and structure, with over 40 million camels worldwide. This number is projected to increase, underscoring their growing significance. Economically, camels are crucial for food production, tourism, and trade, with camel racing being particularly significant in Arab countries. Their unique physiological traits, such as low disease susceptibility and efficient water conservation, further enhance their value. Camel products, especially meat and milk, offer substantial nutritional and therapeutic benefits, contributing to their high demand. Genetic diversity studies have advanced our understanding of camels' adaptation to extreme environments. Functional genomics and whole-genome sequencing have identified genes responsible for these adaptations, aiding breeding programs and conservation efforts. High-throughput sequencing has revealed genetic markers linked to traits like milk production and disease resistance. The development of SNP chips has revolutionized genetic studies by providing a cost-effective alternative to whole-genome sequencing. These tools facilitate large-scale genotyping, essential for conserving genetic diversity and improving breeding strategies. To prevent the depletion of camel genetic diversity, it is crucial to streamline in situ and ex situ conservation efforts to maintain their ecological and economic value. A comprehensive approach to camel conservation and genetic preservation, involving advanced genomic technologies, reproductive biotechniques, and sustainable management practices, will ensure their continued contribution to human societies.

Whole-Genome Resequencing Analysis of the Camelus bactrianus (Bactrian Camel) Genome Identifies Mutations and Genes Affecting Milk Production Traits

Milk production is an important trait that influences the economic value of camels. However, the genetic regulatory mechanisms underlying milk production in camels have not yet been elucidated. We aimed to identify candidate molecular markers that affect camel milk production. We classified Junggar Bactrian camels (9-10-year-old) as low-yield (<1.96 kg/d) or high-yield (>2.75 kg/d) based on milk production performance. Milk fat (5.16 ± 0.51 g/100 g) and milk protein (3.59 ± 0.22 g/100 g) concentrations were significantly lower in high-yielding camels than those in low-yielding camels (6.21 ± 0.59 g/100 g, and 3.93 ± 0.27 g/100 g, respectively) (p < 0.01). There were no apparent differences in gland tissue morphology between the low- and high-production groups. Whole-genome resequencing of 12 low- and 12 high-yield camels was performed. The results of selection mapping methods, performed using two methods (FST and θπ), showed that 264 single nucleotide polymorphism sites (SNPs) overlapped between the two methods, identifying 181 genes. These genes were mainly associated with the regulation of oxytocin, estrogen, ErbB, Wnt, mTOR, PI3K-Akt, growth hormone synthesis/secretion/action, and MAPK signaling pathways. A total of 123 SNPs were selected, based on significantly associated genomic regions and important pathways for SNP genotyping, for verification in 521 additional Bactrian camels. This analysis showed that 13 SNPs were significantly associated with camel milk production yield and 18 SNPs were significantly associated with camel milk composition percentages. Most of these SNPs were located in coding regions of the genome. However, five and two important mutation sites were found in the introns of CSN2 (β-casein) and CSN3 (κ-casein), respectively. Among the candidate genes, NR4A1, ADCY8, PPARG, CSN2, and CSN3 have previously been well studied in dairy livestock. These observations provide a basis for understanding the molecular mechanisms underlying milk production in camels as well as genetic markers for breeding programs aimed at improving milk production.

A discussion on A1-free milk: Nuances and comments beyond implications to the health

This chapter provides an overarching view of the multifaceted aspects of milk β-casein, focusing on its genetic variants A1 and A2. The work examines the current landscape of A1-free milk versus regular milk, delving into health considerations, protein detection methods, technological impacts on dairy production, non-bovine protein, and potential avenues for future research. Firstly, it discussed ongoing debates surrounding categorizing milk based on A1 and A2 β-casein variants, highlighting challenges in establishing clear regulatory standards and quality control methods. The chapter also addressed the molecular distinction between A1 and A2 variants at position 67 of the amino acid chain. This trait affects protein conformation, casein micelle properties, and enzymatic susceptibility. Variations in β-casein across animal species are acknowledged, casting doubt on non-bovine claims of "A2-like" milk due to terminology and genetic differences. Lastly, this work explores the burgeoning field of biotechnology in milk production.

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.

Genetic polymorphism of Β-casein gene and its association with milk production and composition in Azi-Kheli buffalo

The aim of this study was to find out the genetic polymorphism in β-casein gene CSN2 in Azi-Kheli buffaloes found in district Swat. Blood samples from 250 buffaloes were collected and processed in lab for sequencing to see the genetic polymorphism in CSN2 gene on 67 position of exon7. The β-casein is a milk second abundant protein having some variants, wherein A1 and A2 are the most common. After performing sequence analysis, it was found that Azi-Kheli buffaloes were homozygous for only A2 type variant. The amino acid change (proline to histadine) on 67 position of exon 7 was not found; however, three other novel SNPs at loci g.20545A > G, g.20570G > A, and g.20693C > A were identified in the study. Amino acid change due to SNPs were found as SNP1, valine > proline; SNP2, leucin > phenylalanine; and SNP3, threonine > valine. Allelic and genotypic frequencies' analysis exhibited that all three SNPs were following the Hardy-Weinberg equilibrium (HWE: P < 0.05). All the three SNPs showed medium PIC value and gene heterozygosity. The SNPs located on different position of exon 7 of CSN2 gene exhibited associations with some of the performance traits and milk composition. Higher daily milk yield of 9.86 ± 0.43 L and the peak milk yield of 13.80 ± 0.60 L were found in response to SNP3 followed by SNP2 and SNP1. The percentage of milk fat and protein was found significantly higher (P ≤ 0.05) in relation to SNP3 followed by SNP2 and SNP1 given as 7.88 ± 0.41, 7.48 ± 0.33, and 7.15 ± 0.48 for fat% and 4.00 ± 0.15, 3.73 ± 0.10 and 3.40 ± 0.10 for protein%. It was concluded that Azi-Kheli buffalo milk contains A2 genetic variant along with other useful novel variants indicating quality milk for human health. Genotypes of SNP3 should be given preference in selection both in indices and nucleotide polymorphism.

Do non-bovine domestic animals produce A2 milk?: an in silico analysis

A2 milk is an easily digestible product since it has only A2 beta-casein. In cattle, the A1 and A2 alleles are found in the population and the A2 milk is produced from A2A2 animals. Little is known about these alleles in other domestic dairy species. The present study aims to analyze sequence of genetic material available on public databases and quantify the animals genotyped. Eight domestic species were analyzed. There is strong evidence that domestic non-bovine species only carry A2 beta-casein. The data reported here for goats already confirm it due to the large number of animals genotyped as well as buffaloes. It means that they naturally produce A2 milk and no selection must be done. Thus, the fact that A2 milk is easier to digest can be used to add value to dairy product of these species. It helps to conquer new markets. It also improves people's health and breeder profitability.

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.

Tunisian camel casein gene characterization reveals similarities and differences with Sudanese and Nigerian populations

Milk is a primary protein source that has always played a role in mammalian health. Despite the intensification of research projects on dromedary and the knowledge of the genetic diversity at the casein loci, the genetic structure of the Tunisian camel population still needs exploration. This study sought to determine the genetic diversity of 3 casein gene variants in 5 Tunisian camel ecotypes: c.150G>T at CSN1S1 (α_S1-casein), g.2126A>G at CSN2 (β-casein), and g.1029T>C at CSN3 (κ-casein). The obtained results were compared with data published on Sudanese and Nigerian camels to establish the level of differentiation within and between populations. A total of 159 blood samples were collected from 5 Tunisian camel ecotypes and the extracted DNA was genotyped by PCR-RFLP. A streamlined genotyping protocol was also developed for CSN3. Results indicated that allele T was quite rare (0.06) at CSN1S1 for all ecotypes. Minor allele frequency was found for G (0.462) in CSN2 except for Ardhaoui Medenine ecotype who deviated from the average CSN2 allele frequency of the total population. Allele C showed minor allele frequency of 0.384 in CSN3. Among the Tunisian population, GAT (0.343) was the most represented haplotype in all ecotypes except for Ardhaoui Medenine, where GGC (0.322) was the most frequent one. Significant differences in heterozygosity and local inbreeding were observed across the Tunisian, Sudanese, and Nigerian populations, although the global fixation index indicated that only 2.2% of the genetic variance is related to ecotype differences. Instead, phylogenetic analysis revealed a closer link between the Tunisian and Sudanese populations through a clade subdivision with 3 main branches among the ecotypes. This study represents the first attempt to understand casein gene variability in Tunisian camels; with further study, milk traits and genetic differentiation among populations can be associated with the history of camel domestication.

Demographic pattern of A1/A2 beta casein variants indicates conservation of A2 type haplotype across native cattle breeds (Bos indicus) of India

Genetic variations of the beta casein gene hold importance because of their probable association with human health. Comparative sequence analysis of β-casein gene across Indian native, crossbred and exotic breeds in India revealed 15 SNPs and 4 INDELs corresponding to 14 haplotypes. The frequency of A2 type haplotype was maximum (0.941) across all Indian native breeds. Among the 15 variants reported for taurine breeds, only three (A1, A2 and B) were observed in analysed populations. Allelic profiling of A1/A2 β-casein variants in ~ 4000 animals belonging to three cattle types and breeding bulls also revealed the predominance of A2 allele (0.95) in Indian cattle. The high proportion of A2 allele/haplotype indicates that Indian native cattle are the best suited to meet the demands for A2 milk globally. However, a higher percentage of heterozygous genotype (A1A2) in breeding bulls warrants the need to screen sire lines so as to drift the herd towards A2.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03232-0.

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.

Identification and in silico characterization of structural and functional impacts of genetic variants in milk protein genes in the Zebu breeds Guzerat and Gyr

Whole genome sequencing of bovine breeds has allowed identification of genetic variants in milk protein genes. However, functional repercussion of such variants at a molecular level has seldom been investigated. Here, the results of a multistep Bioinformatic analysis for functional characterization of recently identified genetic variants in Brazilian Gyr and Guzerat breeds is described, including predicted effects on the following: (i) evolutionary conserved nucleotide positions/regions; (ii) protein function, stability, and interactions; (iii) splicing, branching, and miRNA binding sites; (iv) promoters and transcription factor binding sites; and (v) collocation with QTL. Seventy-one genetic variants were identified in the caseins (CSN1S1, CSN2, CSN1S2, and CSN3), LALBA, LGB, and LTF genes. Eleven potentially regulatory variants and two missense mutations were identified. LALBA Ile60Val was predicted to affect protein stability and flexibility, by reducing the number the disulfide bonds established. LTF Thr546Asn is predicted to generate steric clashes, which could mildly affect iron coordination. In addition, LALBA Ile60Val and LTF Thr546Asn affect exonic splicing enhancers and silencers. Consequently, both mutations have the potential of affecting immune response at individual level, not only in the mammary gland. Although laborious, this multistep procedure for classifying variants allowed the identification of potentially functional variants for milk protein genes.

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.

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 β-casein gene polymorphism with milk composition traits of Egyptian Maghrebi camels (Camelus dromedarius)

The objectives of this study were to detect the polymorphism of 2126A/G SNP in the β-casein (CSN2) gene among Egyptian Maghrebi camels and to investigate the association of 2126A/G SNP genotypes, parity, lactation stage, and temperature–humidity index (THI) with the milk composition traits of Maghrebi camels. Sixty-eight hair samples were collected from three different populations of Maghrebi camels for DNA extraction. Fat, protein, total solids, solids-not-fat, and lactose percentages were determined in Maghrebi camel milk using an automatic milk analyzer device. Three different genotypes – A/A, A/G, and G/G – were identified in the 5′ flanking region of β-casein gene by using PCR-RFLP method with the A/G genotype showing the highest frequency. Association among these three genotypes with milk composition traits suggests a positive effect of A/A genotype on acidity and protein percentage. Higher protein and acidity values were observed in the milk of individuals carrying the A/A genotype. The protein percentage of this study significantly increased from the first till the fourth parity and then decreased. Fat and total solid percentages were significantly higher in the late stage of lactation, while lactose showed a decreasing trend from the early till the late stages of lactation. Fat and protein percentages were highest in the low THI class. Our results encourage the utilization of Maghrebi camel milk for cheese and butter processing at the late lactation stages of the middle parities of their productive life. Moreover, the A/G SNP of the CSN2 gene may be used as a DNA marker in selection programs for the improvement of camel milk composition. Further studies are needed in order to fully explore the variation in the chemical composition of camel milk due to the effect of CSN2 gene, parity, lactation stage, and THI factors.

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.

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.

From Desert to Medicine: A Review of Camel Genomics and Therapeutic Products

Camels have an important role in the lives of human beings, especially in arid regions, due to their multipurpose role and unique ability to adapt to harsh conditions. In spite of its enormous economic, cultural, and biological importance, the camel genome has not been widely studied. The size of camel genome is roughly 2.38 GB, containing over 20,000 genes. The unusual genetic makeup of the camel is the main reason behind its ability to survive under extreme environmental conditions. The camel genome harbors several unique variations which are being investigated for the treatment of several disorders. Various natural products from camels have also been tested and prescribed as adjunct therapy to control the progression of ailments. Interestingly, the camel employs unique immunological and molecular mechanisms against pathogenic agents and pathological conditions. Here, we broadly review camel classification, distribution and breed as well as recent progress in the determination of the camel genome, its size, genetic distribution, response to various physiological conditions, immunogenetics and the medicinal potential of camel gene products.

5'-flanking variants of equine casein genes (CSN1S1, CSN1S2, CSN2, CSN3) and their relationship with gene expression and milk composition

Genes encoding casein proteins are important candidates for milk composition traits in mammals. In the case of the domestic horse, our knowledge of casein genes is limited mainly to coding sequence variants. This study involved screening for polymorphism in 5’-flanking regions of four genes encoding equine caseins (CSN1S1, CSN1S2, CSN2, and CSN3) and making a preliminary assessment of their effect on the gene expression (on the mRNA and protein levels) and milk composition traits in selected horse breeds. Altogether, 23 polymorphisms (21 described previously SNPs and two novel InDels) were found in the studied sequences, the majority of which are common in various horse breeds. Statistical analysis revealed that some are putatively associated with gene expression or milk composition — for example, the c.-2047_-2048insAT polymorphism (CSN1S1) turns out to be related to the total milk protein content in Polish Primitive Horse (p < 0.05), whereas c.-2105C>G SNP (CSN2) is related to beta-casein relative mRNA level and milk lactose concentration in the Polish Coldblood Horse breed (p < 0.05). We have also found significant effects of horse breed and lactation time-point on gene expression and mare’s milk composition. Our study indicates that the 5’-regulatory regions of genes encoding casein proteins are interesting targets for functional studies of their expression and the composition traits of mare’s milk.

Variability of Some Milk-Associated Genes and Proteins in Several Breeds of Saudi Arabian Camels

To gain knowledge on the molecular basis of diversity of several clans of Saudi camel (Camelus dromedarius) characterization of these animals was conducted at both genetic and protein levels. To this end, blood and milk samples were collected from several camel breeds at different Saudi Arabia locations (northern Jeddah, Riyadh, and Alwagh governorates). Genomic DNA was extracted from blood of four Saudi camel breeds (Majahem, Safra, Wadha, and Hamara), and DNA fragments of the casein and α-lactalbumin genes were amplified. The retrieved DNA sequences were analyzed for genetic variability. The inter-simple sequence repeat technique was used for confirming the relationships among the analyzed camel breeds, and the PCR-RFLP with two restriction enzymes was utilized for exploring their molecular variations. The number of haplotypes, gene diversity, nucleotide diversity, average number of nucleotide differences, and sequence conservation were calculated for all the analyzed DNA sequences. These analyses revealed the presence of several single nucleotide polymorphisms in the analyzed DNA sequences. A group of neighbor joining trees was built for inferring the evolutionary variations among the studied animals. Protein profiling of milk from different camel clans was also conducted, and differences between and within the Saudi camel clans were easily found based on the isoelectric focusing (IEF) profiles using ampholytes with different IEF range. This study revealed that analyzed camel breeds show low levels of genetic differences. This may be a reflection of the evolutionary history of C. dromedarius that was domesticated based on a highly homogeneous ancestor ecotype.

Combining different proteomic approaches to resolve complexity of the milk protein fraction of dromedary, Bactrian camels and hybrids, from different regions of Kazakhstan

Nutritional suitability of milk is not only related to gross composition, but is also strongly affected by the microheterogeniety of the protein fraction. Hence, to go further into the evaluation of the potential suitability of non-bovine milks in human/infant nutrition it is necessary to have a detailed characterization of their protein components. Combining proven proteomic approaches (SDS-PAGE, LC-MS/MS and LC-ESI-MS) and cDNA sequencing, we provide here in depth characterization of the milk protein fraction of dromedary and Bactrian camels, and their hybrids, from different regions of Kazakhstan. A total 391 functional groups of proteins were identified from 8 camel milk samples. A detailed characterization of 50 protein molecules, relating to genetic variants and isoforms arising from post-translational modifications and alternative splicing events, belonging to nine protein families (κ-, α_s1-, α_s2-, β-; and γ-CN, WAP, α-LAC, PGRP, CSA/LPO) was achieved by LC-ESI-MS. The presence of two unknown proteins UP1 (22,939 Da) and UP2 (23,046 Da) was also reported as well as the existence of a β-CN short isoform (946 Da lighter than the full-length β-CN), arising very likely in both genetic variants (A and B) from proteolysis by plasmin. In addition, we report, for the first time to our knowledge, the occurrence of a α_s2-CN phosphorylation isoform with 12P groups within two recognition motifs, suggesting thereby the existence of two kinase systems involved in the phosphorylation of caseins in the mammary gland. Finally, we demonstrate that genetic variants, which hitherto seemed to be species- specific (e.g. β-CN A for Bactrian and β-CN B for dromedary), are in fact present both in Camel dromedarius and C. bactrianus.

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.

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.

Promoter analysis and tissue expression of the chicken ASB15 gene

1. This study was conducted to explore the promoter region of the chicken ASB15 gene by detecting the activities of the dual luciferase reporter gene and to assess expression profiles of the ASB15 gene in 10 different tissues from Gushi chickens. 2. Five dual luciferase reporter gene vectors were constructed and transfected into DF1 cells. The activities of recombined plasmids were measured and the core promoter was confirmed by bioinformatic analysis. Total RNA was extracted and the relative expression of the ASB15 gene was examined. 3. Data analysis indicated that the promoter was located from -955 to -212 bp. Results showed that the chicken ASB15 gene was expressed in heart, breast muscle and leg muscle. 4. This study has confirmed the promoter region and the expression profile of the chicken ASB15 gene, which provides a foundation for further exploring its transcriptional regulation and function.

The history of Old World camelids in the light of molecular genetics

Old World camels have come into the focus as sustainable livestock species, unique in their morphological and physiological characteristics and capable of providing vital products even under extreme environmental conditions. The evolutionary history of dromedary and Bactrian camels traces back to the middle Eocene (around 40 million years ago, mya), when the ancestors of Camelus emerged on the North American continent. While the genetic status of the two domestic species has long been established, the wild two-humped camel has only recently been recognized as a separate species, Camelus ferus, based on molecular genetic data. The demographic history established from genome drafts of Old World camels shows the independent development of the three species over the last 100,000 years with severe bottlenecks occurring during the last glacial period and in the recent past. Ongoing studies involve the immune system, relevant production traits, and the global population structure and domestication of Old World camels. Based on the now available whole genome drafts, specific metabolic pathways have been described shedding new light on the camels' ability to adapt to desert environments. These new data will also be at the origin for genome-wide association studies to link economically relevant phenotypes to genotypes and to conserve the diverse genetic resources in Old World camelids.

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.