Multispecies comparison of the casein gene loci and evolution of casein gene family

Variations in Casein Genes Are Associated with Milk Protein and Fat Contents in Sarda Goats (Capra hircus), with an Important Role of CSN1S2 for Milk Yield

This work aimed to assess the variability of casein genes in a population of 153 bucks and 825 lactating does of the Sarda breed, and to perform association analysis between polymorphic sites and milk yield and composition traits. To genotype the casein genes, we chose an SNP panel including 44 SNPs mapping to the four casein genes CSN1S1, CSN2, CSN1S2, and CSN3. Genotyping (made by KASP™ genotyping assay, based on competitive allele-specific PCR) revealed the high variability of the Sarda goat, and haplotype analysis revealed linkage disequilibrium (LD) between CSN1S1 and CSN2 genes, in addition to two LD blocks within the CSN1S2 and two LD blocks within the CSN3 gene, in bucks and does. Association analysis revealed that variability at all four casein genes was associated with milk protein content, total solids, and milk energy. The three Ca-sensitive casein genes were associated with lipid content, and CSN1S2 showed a unique pattern, with intron variants associated with milk yield, in addition to milk pH, NaCl, and SCS (Somatic Cell Score). This information might prove useful in selection schemes and in future investigations aiming to better understand the biology of lactation, and the direct link between genotype and phenotype.

Potential status of A1 and A2 variants of bovine beta-casein gene in milk samples of Indian cattle breeds

This study aimed to determine the polymorphism in 7th exon of beta-casein gene (CSN2) gene in seven domestic (Kosali, Tharparkar, Gangatiri, Sahiwal, Gir, Khariar, Motu) and two exotic cattle breeds (Jersey and Holstein-Friesian). Genomic DNA was extracted from 1000 milk samples, and the C > A polymorphism in CSN2 was determined using the tetra-primer amplification refractory mutation system-polymerase chain reaction method. In all Indigenous cattle breeds, the mean frequency of A1A2 and A2A2 genotypes was 0.19 and 0.80, respectively. The A1A1 genotype was absent in all seven domestic cattle breeds. The frequency of the A2A2 genotype was highest in the Gir breed (0.93). However, the Sahiwal, Tharparkar, and Motu breeds also had a higher frequency of A2A2 genotype compared to other breeds. In contrast, Gangatiri breed of India showed lowest frequency of A2A2 genotype. The mean A1 and A2 allele frequency was 0.09 and 0.91, respectively. In exotic breeds, the mean frequencies of the A1A1, A1A2, and A2A2 genotypes were 0.42, 0.55, and 0.03, respectively. Similarly, the mean A1 and A2 allele frequency was 0.69 and 0.31, respectively. This study suggests the high potential of Gir, Sahiwal, Tharparkar, and Motu cattle for A2 milk production since they carry a favorable A2 genotype.

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.

Identification and molecular evolution of the GLX genes in 21 plant species: a focus on the Gossypium hirsutum

BACKGROUND

The glyoxalase system includes glyoxalase I (GLXI), glyoxalase II (GLXII) and glyoxalase III (GLXIII), which are responsible for methylglyoxal (MG) detoxification and involved in abiotic stress responses such as drought, salinity and heavy metal.

RESULTS

In this study, a total of 620 GLX family genes were identified from 21 different plant species. The results of evolutionary analysis showed that GLX genes exist in all species from lower plants to higher plants, inferring that GLX genes might be important for plants, and GLXI and GLXII account for the majority. In addition, motif showed an expanding trend in the process of evolution. The analysis of cis-acting elements in 21 different plant species showed that the promoter region of the GLX genes were rich in phytohormones and biotic and abiotic stress-related elements, indicating that GLX genes can participate in a variety of life processes. In cotton, GLXs could be divided into two groups and most GLXIs distributed in group I, GLXIIs and GLXIIIs mainly belonged to group II, indicating that there are more similarities between GLXII and GLXIII in cotton evolution. The transcriptome data analysis and quantitative real-time PCR analysis (qRT-PCR) show that some members of GLX family would respond to high temperature treatment in G.hirsutum. The protein interaction network of GLXs in G.hirsutum implied that most members can participate in various life processes through protein interactions.

CONCLUSIONS

The results elucidated the evolutionary history of GLX family genes in plants and lay the foundation for their functions analysis in cotton.

Caseins: Versatility of Their Micellar Organization in Relation to the Functional and Nutritional Properties of Milk

The milk of mammals is a complex fluid mixture of various proteins, minerals, lipids, and other micronutrients that play a critical role in providing nutrition and immunity to newborns. Casein proteins together with calcium phosphate form large colloidal particles, called casein micelles. Caseins and their micelles have received great scientific interest, but their versatility and role in the functional and nutritional properties of milk from different animal species are not fully understood. Caseins belong to a class of proteins that exhibit open and flexible conformations. Here, we discuss the key features that maintain the structures of the protein sequences in four selected animal species: cow, camel, human, and African elephant. The primary sequences of these proteins and their posttranslational modifications (phosphorylation and glycosylation) that determine their secondary structures have distinctively evolved in these different animal species, leading to differences in their structural, functional, and nutritional properties. The variability in the structures of milk caseins influence the properties of their dairy products, such as cheese and yogurt, as well as their digestibility and allergic properties. Such differences are beneficial to the development of different functionally improved casein molecules with variable biological and industrial utilities.

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.

Assessment of A1 and A2 variants in the CNS2 gene of some cattle breeds by using ACRS-PCR method

The aim of this study is to reveal β-casein polymorphism of some cattle breeds and also the potential to produce A2 milk from existing animals and to develop strategies in this area. Therefore, a total of 400 cattle, 100 animals from each breed of Holstein, Brown Swiss, Jersey and Simmental raised commonly in Turkey, were obtained, and C > A polymorphism in 67th amino acid in the 7th exons of β-casein gene was determined by TaqI enzyme with PCR-ACRS method. Blood samples were collected from dairy cattle farms raising Holstein, Brown Swiss and Jersey breeds from Konya province and Simmental breed from Kütahya province in Turkey. A1 and A2 allele frequencies in Holstein, Brown Swiss, Jersey and Simmental cattle breeds were determined as 0.475/0.525, 0.370/0.630, 0.215/0.785 and 0.440/0.560, respectively. The highest A2 allele frequency (0.785) was found in Jersey breed. A1A1 genotypes in Holstein, Brown Swiss, Jersey and Simmental breeds were 0.240, 0.150, 0.030 and 0.160, respectively; A1A2 genotypes were 0.470, 0.440, 0.370 and 0.560, respectively; A2A2 genotypes were determined as 0.290, 0.410, 0.600 and 0.280, respectively. In these breeds, the highest A2A2 genotype frequency was found in Jersey (0.600), the lowest A1A1 genotype frequency (0.030) was found in Jersey and the highest A1A2 genotype frequency (0.560) was found in Simmental. Holstein, Brown Swiss, Simmental and Jersey populations were at the level of Hardy-Weinberg in terms of β-casein gene (p > 0.05). The average H_o, H_e and PIC values were calculated as 0.460, 0.469 and 0.605, respectively. In conclusion, the frequency of commonly reared cattles in Turkey especially Brown Swiss, and Jersey breeds in A2A2 genotype are satisfactory, but it can be said that the use of animals with A2 allele in selection is extremely important for increasing A2 milk producing potential in the future.

Genetic polymorphism of bovine beta-casein gene in Japanese dairy farm herds

The aim of this study was to investigate beta‐casein polymorphism among 320 Japanese cows sampled from eight dairy farms. We used a newly‐developed genotyping method that involved collecting DNA from hairs and a Cycleave polymerase chain reaction (PCR) assay to detect the A1, A2, and B variants. Results revealed the presence of five genotypes (A1A1, A2A2, A1A2, A1B, and A2B). We found that the most common genotype was A2A2 (0.42), followed by A1A2 (0.39) and A1A1 (0.11). The A1B and A2B genotypes were less frequent (<0.05). The frequencies of alleles A1, A2, and B were calculated to be 0.32, 0.64, and 0.04, respectively. Our study is the first to show the current status of beta‐casein polymorphisms in Japanese dairy farms. Given the adverse effects of A1 beta‐casein on human health, attempts have been made to develop herds consisting solely of A2A2 cows. Our study provides a reference for improving cow populations in Japanese dairy farms. The Cycleave PCR‐based assay we developed here can be used for rapid and reliable genotyping of bovine beta‐casein.

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.

The molecular basis of mammary gland development and epithelial differentiation

Our understanding of the molecular events underpinning the development of mammalian organ systems has been increasing rapidly in recent years. With the advent of new and improved next-generation sequencing methods, we are now able to dig deeper than ever before into the genomic and epigenomic events that play critical roles in determining the fates of stem and progenitor cells during the development of an embryo into an adult. In this review, we detail and discuss the genes and pathways that are involved in mammary gland development, from embryogenesis, through maturation into an adult gland, to the role of pregnancy signals in directing the terminal maturation of the mammary gland into a milk producing organ that can nurture the offspring. We also provide an overview of the latest research in the single-cell genomics of mammary gland development, which may help us to understand the lineage commitment of mammary stem cells (MaSCs) into luminal or basal epithelial cells that constitute the mammary gland. Finally, we summarize the use of 3D organoid cultures as a model system to study the molecular events during mammary gland development. Our increased investigation of the molecular requirements for normal mammary gland development will advance the discovery of targets to predict breast cancer risk and the development of new breast cancer therapies.

Comparative Genomics, Evolutionary and Gene Regulatory Regions Analysis of Casein Gene Family in Bubalus bubalis

Buffalo is a luxurious genetic resource with multiple utilities (as a dairy, draft, and meat animal) and economic significance in the tropical and subtropical regions of the globe. The excellent potential to survive and perform on marginal resources makes buffalo an important source for nutritious products, particularly milk and meat. This study was aimed to investigate the evolutionary relationship, physiochemical properties, and comparative genomic analysis of the casein gene family (CSN1S1, CSN2, CSN1S2, and CSN3) in river and swamp buffalo. Phylogenetic, gene structure, motif, and conserved domain analysis revealed the evolutionarily conserved nature of the casein genes in buffalo and other closely related species. Results indicated that casein proteins were unstable, hydrophilic, and thermostable, although αs1-CN, β-CN, and κ-CN exhibited acidic properties except for αs2-CN, which behaved slightly basic. Comparative analysis of amino acid sequences revealed greater variation in the river buffalo breeds than the swamp buffalo indicating the possible role of these variations in the regulation of milk traits in buffalo. Furthermore, we identified lower transcription activators STATs and higher repressor site YY1 distribution in swamp buffalo, revealing its association with lower expression of casein genes that might subsequently affect milk production. The role of the main motifs in controlling the expression of casein genes necessitates the need for functional studies to evaluate the effect of these elements on the regulation of casein gene function in buffalo.

Beta-Casein Gene Polymorphism in Serbian Holstein-Friesian Cows and Its Relationship with Milk Production Traits

The most common types of beta-casein in cow’s milk are A1 and A2, which differ in one amino acid. During the gastrointestinal proteolysis of A1 beta-casein in humans this difference results in the release of beta-casomorphin-7, an opioid which may lead to severe effects on human health, causing various ailments (type-1 diabetes mellitus, ischemic heart disease, arteriosclerosis, sudden infant death syndrome, autism, schizophrenia, gastrointestinal digestive discomfort, as well as increased gastrointestinal transit time). By contrast, A2 beta-casein cannot exert these effects owing to its different composition and metabolism. Furthermore, studies have shown that it can influence milk productivity traits. Our research aimed to screen the frequency of A1 and A2 alleles of beta-casein gene in a population of Serbian Holstein-Friesian cows and to detect how the genotypes influence milk production, and milk protein and fat yields. Out of 106 animals, 13 (12.26%) were of A1A1 genotype, 58 (54.72%) of A1A2, and 35 (33.02%) of A2A2 genotype. Milk yield was significantly (P<0.01) higher in A2A2 compared to both A1A1 and A1A2 genotypes. Milk protein concentrations were significantly (P<0.01) higher in A2A2 compared to A1A2 genotype, while milk fat concentrations were significantly (P<0.01) higher in A2A2 compared to both A1A1 and A1A2 genotypes.

An improved LC-MS method to profile molecular diversity and quantify the six main bovine milk proteins, including genetic and splicing variants as well as post-translationally modified isoforms

•

Identification of the main milk proteins isoforms are inferred from a mass database.

•

Quantification of the six main milk proteins is achieved from corrected UV at 214 nm.

•

Multiple isoforms can be quantified from mass signal, within each protein family.

•

Glycosylation isoforms of κ-casein which impact micelle stability are quantified.

Here we describe a method based on Liquid Chromatography coupled with Mass Spectrometry (LC-MS) that provides an accurate determination of the six main bovine milk proteins, including allelic and splicing variants, as well as isoforms resulting from post-translational modifications, with an unprecedented level of resolution. Proteins are identified from observed molecular masses in comparison with theoretical masses of intact proteins indexed in an “in-house” database that includes nearly 3000 entries. Quantification was performed either from UV (214 nm) or mass signals. Thus, up to one hundred molecules, derived from the six major milk proteins, can be identified and quantified from an individual milk sample. This powerful and reliable method, initially developed as an anchoring method to estimate the composition of the six main bovine milk proteins from MIR spectra, is transferable to several mammalian species, including small ruminants, camels, equines, rabbits, etc., for which specific mass databases are available.

Frequencies Evaluation of β-Casein Gene Polymorphisms in Dairy Cows Reared in Central Italy

Simple Summary

Bovine milk contains several β-casein variants, with the A1 and A2 variants occurring most frequently. The presence of some variants, such as A1, B, and C, is considered a risk factor for disease in humans who consume milk. These variants are probably involved in intolerance to milk and some human diseases due to the production of a bioactive peptide with opioid activity during digestion, β-casomorphin 7 (BCM-7). In contrast, the A2 variant is not involved in pathogenetic mechanisms; thus, its presence in milk is a desirable feature. The difference between the A1 and A2 variants is a mutation at position 67 of the β-casein gene (CSN2), which causes an amino acid to change from histidine (in the A1, B, and C variants) to proline (in the A2 variant). To select dairy cows on the basis of the presence of the β-casein variant A2, allele frequencies of CSN2 variants were evaluated in Italian dairy cows reared in central Italy. The results of this study may help with the selection of animals with the β-casein gene variant A2 to produce a more digestible milk that only contains the β-casein variant A2.

Abstract

The majority of proteins in cow’s milk are caseins, which occur in four groups (α-s1, α-s2, β, and k) encoded by different genes (CSN1S1, CSN1S2, CSN2, and CSN3, respectively). In this study, we focused on the β-casein allele variants A1 and A2 due to their influence on milk’s technological characteristics and human health. Digestion of the β-casein variant A1 leads to the formation of β-casomorphin 7 (BCM-7), a bioactive peptide that has been suggested to be a possible cause of various human diseases and associated with low milk digestibility. The potential negative role of the β-casein variant A1 in human health has stimulated the planning of cattle breeding programs based on genetic selection to increase the frequency of the A2 variant, which is associated with increased milk digestibility. The aim of this work was to evaluate the frequencies of the different β-casein variants in Italian Holstein Friesian dairy cows from cattle farms located in central Italy to select a population of A2 homozygous animals. β-casein genotypes were identified by evaluating the presence of single nucleotide polymorphisms (SNPs) of the CSN2 gene using PCR and sequencing analysis. The frequency of the desirable β-casein variant A2 in the studied bovine population was 0.61. The frequency of the undesirable A1 variant in the studied bovine population was 0.30. The frequency of the A2 allele was higher than expected for the breed; therefore, genetic selection for the A2 variant in these animals could be achieved in a fairly short time using A2 homozygous bulls.

Genetic potential of Sindhi cattle for A2 milk production

ContextCow milk is consumed in many countries. However, its consumption can be associated with gastrointestinal disorders in humans. The genetic profile of β-casein is one of the factors responsible for these disorders and the A1 variant is a causal agent. The β-casein A1 leads to the formation of bioactive peptides, which trigger intestinal reactions. The same is not observed for the A2 allele. A2 milk is obtained from cows with the A2A2 genotype. The identification of breeds and selection of animals for A2 milk production is important for public health and animal production. The Sindhi cattle (Bos taurus indicus) possesses characteristics such as adaptation to hot climate, resistance to parasites and utilisation of low-quality forage that are important for rearing in the tropics. AimsThe objective of the present study was to genotype the single-nucleotide polymorphism (SNP) that forms A1 and A2 alleles of the β-casein in Sindhi cattle. MethodsIn total, 348 Sindhi animals were used. DNA was extracted from hair follicles and submitted to PCR-sequencing for genotyping. Key resultsAllele frequency for A2 was 0.94 and for A1 it was 0.06. Genotype frequencies were 0.90 (A2A2), 0.09 (A1A2) and 0.01 (A1A1). The frequency of the A2A2 genotype in Sindhi cattle was similar to or higher than those found in zebu breeds and higher than those reported for specialised taurine breeds. ConclusionsThis high frequency of the A2A2 genotype indicates a genetic potential of the breed for A2 milk production and easy, assisted selection (if that is the intention). ImplicationsThis information gives importance to the breed and permits to add value to dairy product derived from it. It may increase producer’s income through producing and marketing A2 milk.

Implications of kappa-casein evolutionary diversity for the self-assembly and aggregation of casein micelles

Milk alpha-, beta- and kappa-casein proteins assemble into casein micelles in breast epithelial cells. The glycomacropeptide (GMP) tails of kappa-casein that extend from the surface of the micelle are key to assembly and aggregation. Aggregation is triggered by stomach pepsin cleavage of GMP from para-kappa-casein (PKC). While one casein micelle model emphasizes the importance of hydrophobic interactions, another focuses on polar residues. We performed an evolutionary analysis of kappa-casein primary sequence and predicted features that potentially impact on protein interactions. We noted more rapid change in the earlier period (166 to 60 Ma). Pepsin and plasmin cleavage sites were avoided in the GMP, which may partly explain its amino acid composition. Short tandem repeats have led to modest expansions of PKC, and to large GMP expansions, suggesting the GMP is less length constrained. Amino acid compositional constraints were assessed across species. Polarity and hydrophobicity properties were insufficient to explain differences between PKC and GMP. Among polar residues, threonine dominates the GMP, compared to serine, probably reflecting its preference for O-glycosylation over phosphorylation. Glutamine, enriched in the bovine PQ-rich region, is not positionally conserved in other species. Among hydrophobic residues, isoleucine is clearly preferred over leucine in the GMP, and patches of hydrophobicity are not markedly positionally conserved. PKC tyrosine and charged residues showed stronger conservation of position, suggesting a role for pi-interactions, seen in other structurally dynamic protein membraneless assemblies. Independent acquisitions of cysteines are consistent with a trend of increasing stabilization of multimers by covalent disulphide bonds, over evolutionary time. In conclusion, kappa-casein compositional and positional constraints appear to be influenced by modification preferences, protease evasion and protein-protein interactions.

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.

Genomic analysis of the origins of extant casein variation in goats

The variation in the casein genes has a major impact on the milk composition of goats. Even though many casein polymorphisms have been identified so far, we do not know yet whether they are evolutionarily ancient (i.e., they existed before domestication) or young (i.e., they emerged after domestication). Herewith, we identified casein polymorphisms in a data set of 106 caprine whole-genome sequences corresponding to bezoars (Capra aegagrus, the ancestor of domestic goats) and 4 domestic goat (Capra hircus) populations from Europe, Africa, the Far East, and the Near East. Domestic and wild goat populations shared a substantial number of casein SNP, from 36.1% (CSN2) to 55.1% (CSN1S2). The comparison of casein variation among bezoars and the 4 domestic goat populations demonstrated that more than 50% of the casein SNP are shared by 2 or more populations, and 18 to 44% are shared by all populations. Moreover, the majority of casein alleles reported in domestic goats also segregate in the bezoar, including several alleles displaying significant associations with milk composition (e.g., the A/B alleles of the CSN1S1 and CSN3 genes, the A allele of the CSN2 gene). We conclude that much of the current diversity of the caprine casein genes comes from ancient standing variation segregating in the ancestor of modern domestic goats.

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.

Rapid autophagic regression of the milk gland during involution is critical for maximizing tsetse viviparous reproductive output

Tsetse flies are important vectors of human and animal trypanosomiasis. Ability to reduce tsetse populations is an effective means of disease control. Lactation is an essential component of tsetse’s viviparous reproductive physiology and requires a dramatic increase in the expression and synthesis of milk proteins by the milk gland organ in order to nurture larval growth. In between each gonotrophic cycle, tsetse ceases milk production and milk gland tubules undergo a nearly two-fold reduction in width (involution). In this study, we examined the role autophagy plays during tsetse fly milk gland involution and reproductive output. Autophagy genes show elevated expression in tissues associated with lactation, immediately before or within two hours post-parturition, and decline at 24-48h post-parturition. This expression pattern is inversely correlated with that of the milk gland proteins (lactation-specific protein coding genes) and the autophagy inhibitor fk506-bp1. Increased expression of Drosophila inhibitor of apoptosis 1, diap1, was also observed in the milk gland during involution, when it likely prevents apoptosis of milk gland cells. RNAi-mediated knockdown of autophagy related gene 8a (atg8a) prevented rapid milk gland autophagy during involution, prolonging gestation, and reducing fecundity in the subsequent gonotrophic cycle. The resultant inhibition of autophagy reduced the recovery of stored lipids during the dry (non-lactating) periods by 15–20%. Ecdysone application, similar to levels that occur immediately before birth, induced autophagy, and increased milk gland involution even before abortion. This suggests that the ecdysteroid peak immediately preceding parturition likely triggers milk gland autophagy. Population modeling reveals that a delay in involution would yield a negative population growth rate. This study indicates that milk gland autophagy during involution is critical to restore nutrient reserves and allow efficient transition between pregnancy cycles. Targeting post-birth phases of reproduction could be utilized as a novel mechanism to suppress tsetse populations and reduce trypanosomiasis.

Tsetse flies are vectors for trypanosomes that cause both African sleeping sickness in humans and Nagana in animals. The reduction of tsetse populations is the most efficient way to reduce the prevalence of this economically important disease with current control methods including pesticide application, traps, and sterile insect techniques. Tsetse pregnancy and milk production represent a species-specific target for population control and milk gland transition during each larval growth cycle could represent a novel target for tsetse control. Within one day after birth, the milk gland organ, essential for provisioning nutrients to the intrauterine larva, undergoes involution marked by an ecdysone driven increase in autophagy that allows breakdown of this gland. Inhibiting the process of autophagy prevents the timely transition from the lactation phase to the dry phase, triggering a delay in subsequent pregnancy cycle. This misregulation of milk gland involution leads to an overall decrease in the number of offspring that each female can produce per lifetime. This study has determined the molecular components of this process, and reveals new targets of interference for vector control.

Evaluation of bovine beta casein polymorphism in two dairy farms located in northern Italy

Bovine beta casein A1 is one of the most common variants in dairy cattle breeds; it is considered a risk factor in milk intolerance and in other important human diseases, because of the bioactive peptide beta casomorphin-7 (BCM7) produced by raw or processed A1-milk, but not by A2-milk, during digestion. The aim of this study was to perform a cheap and rapid method to investigate beta casein polymorphism in copious animals. The study included 2 dairy farms with a totally of 1230 cows. Beta casein genotypes were estimated evaluating Exon 7 region of bovine beta casein gene (CSN2) by sequences analysis. In the population included in the study 5 variants (A1, A2, B, F, I) and 13 genotypes (A1A1, A1A2, A1B, A1F, A1I, A2A2, A2B, A2F, A2I, BB, BF, BI, FI) were detected. The method showed high sensibility and specificity, resulted low-cost and few time consuming.

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

Distinct roles of prolactin, epidermal growth factor, and glucocorticoids in β-casein secretion pathway in lactating mammary epithelial cells

Beta-casein is a secretory protein contained in milk. Mammary epithelial cells (MECs) synthesize and secrete β-casein during lactation. However, it remains unclear how the β-casein secretion pathway is developed after parturition. In this study, we focused on prolactin (PRL), epidermal growth factor (EGF), and glucocorticoids, which increase in blood plasma and milk around parturition. MECs cultured with PRL, EGF and dexamethasone (DEX: glucocorticoid analog) developed the β-casein secretion pathway. In the absence of PRL, MECs hardly expressed β-casein. EGF enhanced the expression and secretion of β-casein in the presence of PRL and DEX. DEX treatment rapidly increased secreted β-casein concurrent with enhancing β-casein expression. DEX also up-regulated the expression of SNARE proteins, such as SNAP-23, VAMP-8 and Syntaxin-12. Furthermore, PRL and DEX regulated the expression ratio of α_s1-, β- and κ-casein. These results indicate that PRL, EGF and glucocorticoids have distinct roles in the establishment of β-casein secretion pathway.

Perinatal exposure to bisphenol A modifies the transcriptional regulation of the β-Casein gene during secretory activation of the rat mammary gland

With the aim to analyze whether bisphenol A (BPA) modifies β-Casein (β-Cas) synthesis and transcriptional regulation in perinatally exposed animals, here, pregnant F0 rats were orally exposed to 0, 0.6 or 52 μg BPA/kg/day from gestation day 9 until weaning. Then, F1 females were bred and mammary glands were obtained on lactation day 2. Perinatal BPA exposure decreased β-Cas expression without modifying the activation of prolactin receptor. It also decreased the expression of glucocorticoid receptor in BPA52-exposed dams and β1 and α6 integrins as well as dystroglycan in both BPA groups. In addition, BPA exposure altered the expression of histone-modifying enzymes and induced histone modifications and DNA methylation in the promoter, enhancer and exon VII of the β-Cas gene. An impaired crosstalk between the extracellular matrix and lactogenic hormone signaling pathways and epigenetic modifications of the β-Cas gene could be the molecular mechanisms by which BPA decreased β-Cas expression.

Translation attenuation via 3' terminal codon usage in bovine csn1s2 is responsible for the difference in αs2- and β-casein profile in milk

The rate of secretion of αs2-casein into bovine milk is approximately 25% of that of β-casein, yet mammary expression of their respective mRNA transcripts (csn1s2 and csn2) is not different. Our objective was to identify molecular mechanisms that explain the difference in translation efficiency between csn1s2 and csn2. Cell-free translational efficiency of csn2 was 5 times that of csn1s2. Transcripts of csn1s2 distributed into heavier polysomes than csn2 transcripts, indicating an attenuation of elongation and/or termination. Stimulatory and inhibitory effects of the 5' and 3' UTRs on translational efficiency were different with luciferase and casein sequences in the coding regions. Substituting the 5' and 3' UTRs from csn2 into csn1s2 did not improve csn1s2 translation, implicating the coding region itself in the translation difference. Deletion of a 28-codon fragment from the 3' terminus of the csn1s2 coding region, which displays codons with low correlations to cell fitness, increased translation to a par with csn2. We conclude that the usage of the last 28 codons of csn1s2 is the main regulatory element that attenuates its expression and is responsible for the differential translational expression of csn1s2 and csn2.

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.

Cross-species comparison of mammalian saliva using an LC-MALDI based proteomic approach

Despite the importance of saliva in the regulation of oral cavity homeostasis, few studies have been conducted to quantitatively compare the saliva of different mammal species. Aiming to define a proteome signature of mammals' saliva, an in-depth SDS-PAGE-LC coupled to MS/MS (GeLC-MS/MS) approach was used to characterize the saliva from primates (human), carnivores (dog), glires (rat and rabbit), and ungulates (sheep, cattle, horse). Despite the high variability in the number of distinct proteins identified per species, most protein families were shared by the mammals studied with the exception of cattle and horse. Alpha-amylase is an example that seems to reflect the natural selection related to digestion efficacy and food recognition. Casein protein family was identified in all species but human, suggesting an alternative to statherin in the protection of hard tissues. Overall, data suggest that different proteins might assure a similar role in the regulation of oral cavity homeostasis, potentially explaining the specific mammals' salivary proteome signature. Moreover, some protein families were identified for the first time in the saliva of some species, the presence of proline-rich proteins in rabbit's saliva being a good example.

An evolutionarily conserved non-coding element in casein locus acts as transcriptional repressor

In mammals, the casein locus consists of stretches of non-coding DNA, the functions of most of which are unknown. These regions are believed to harbour elements responsible for spatio-temporally regulated expression of genes in this locus and so far, only a few such elements have been identified. In this study, we report a novel regulatory element in the casein locus. Comparative analysis of genomic DNA sequences of casein loci from different mammals identified a 147bp long evolutionarily conserved region (ECR) upstream of Odam, a gene in this locus. The ECR was found in close proximity of Odam gene in all the mammals examined. In-silico analysis predicted the ECR as a potential regulatory element. Functional analysis in different cell lines identified it as a unidirectional repressor element. From our findings we speculate that the ECR may be involved in the repression of the Odam expression in the mammary gland during lactation.

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.

Investigating the genetic polymorphism of sheep milk proteins: a useful tool for dairy production

Sheep is the second most important dairy species after cow worldwide, and especially in the Mediterranean and Middle East regions. In some countries, the difficult environmental conditions require a peculiar adaptation and, in these contexts, sheep are able to provide higher quality protein than cattle. In the least-developed countries, the amount of dairy sheep and ovine milk production is progressively increasing. In order to improve dairy productions, in particular those with local connotations, it is necessary to obtain in-depth information regarding milk quality and rheological properties. The genetic polymorphisms of milk proteins are often associated with quantitative and qualitative parameters in milk and are potential candidate markers that should be included in breeding strategies similar to those already available for cattle. Due to the current and growing interest in this topic and considering the large amount of new information, the aim of this study was to review the literature on sheep milk protein polymorphisms with a particular emphasis on recent findings in order to give scientists useful support. Moreover, the effects of different protein variants on milk yield and composition are discussed.

Effects of CSN1S2 Genotypes on Economic Traits in Chinese Dairy Goats

The aim of this study was to investigate allele frequencies at the CSN1S2 locus in two Chinese dairy goat breeds and the effects of its variation on dairy goat economic traits. Seven hundred and eight goats from Xinong Saanen (XS, n = 268) and Guanzhong (GZ, N = 440) breeds were selected. The milk samples of 268 XS goats were collected during the middle of lactation, body size parameters (708 goats) and daily milk yield (202 goats) were registered. The RFLP (restriction fragment length polymorphism) and SSCP (single strand conformation polymorphism) were used to detect the polymorphisms in CSN1S2. The Hardy-Weinberg (HW) equilibrium and the associations between body size, milk yield and composition and the genotypes were calculated. The results revealed that only A and F CSN1S2 alleles were found in the two Chinese dairy goat breeds. Allelic frequencies of A and F were 0.795, 0.205 and 0.739, 0.261 in Xinong Saanen and Guanzhong population respectively. Xinong Saanen breed was in Hardy-Weinberg equilibrium, while Guanzhong breed deviated from Hardy-Weinberg equilibrium (p<0.05). The association of polymorphism with economic traits indicated that the goats with FF genotype have higher milk fat and total solid concentration than those with AA and AF genotypes (p<0.05).

Cloning and Molecular Characterization of Porcine β-casein Gene (CNS2)

The production of therapeutic proteins from transgenic animals is one of the most important successes of animal biotechnology. Milk is presently the most mature system for production of therapeutic proteins from a transgenic animal. Specifically, β-casein is a major component of cow, goat and sheep milk, and its promoter has been used to regulate the expression of transgenic genes in the mammary gland of transgenic animals. Here, we cloned the porcine β-casein gene and analyzed the transcriptional activity of the promoter and intron 1 region of the porcine β-casein gene. Sequence inspection of the 5′-flanking region revealed potential DNA elements including SRY, CdxA, AML-a, GATA-3, GATA-1 and C/EBP β. In addition, the first intron of the porcine β-casein gene contained the transcriptional enhancers Oct-1, SRY, YY1, C/EBP β, and AP-1, as well as the retroviral TATA box. We estimated the transcriptional activity for the 5′-proximal region with or without intron 1 of the porcine β-casein gene in HC11 cells stimulated with lactogenic hormones. High transcriptional activity was obtained for the 5′-proximal region with intron 1 of the porcine β-casein gene. The β-casein gene containing the mutant TATA box (CATAAAA) was also cloned from another individual pig. Promoter activity of the luciferase vector containing the mutant TATA box was weaker than the same vector containing the normal TATA box. Taken together, these findings suggest that the transcription of porcine β-casein gene is regulated by lactogenic hormone via intron 1 and promoter containing a mutant TATA box (CATAAAA) has poor porcine β-casein gene activity.

Goat casein genotypes are associated with milk production traits in the Sarda breed

The aim of the current work was to analyze, in the Sarda breed goat, genetic polymorphism within the casein genes and to assess their influence on milk traits. Genetic variants at the CSN1S1, CSN2, CSN1S2 and CSN3 gene loci were investigated using PCR-based methods, cloning and sequencing. Strong alleles prevailed at the CSN1S1 gene locus and defective alleles also were revealed. Null alleles were evidenced at each calcium-sensitive gene locus. At the CSN3 gene locus, we observed a prevalence of the CSN3 A and B alleles; the occurrence of rare alleles such as CSN3 B'', C, C', D, E and M; and the CSN3 S allele (GenBank KF644565) described here for the first time in Capra hircus. Statistical analysis showed that all genes, except CSN3, significantly influenced milk traits. The CSN1S1 BB and AB genotypes were associated with the highest percentages of protein (4.41 and 4.40 respectively) and fat (5.26 and 5.34 respectively) (P < 0.001). A relevant finding was that CSN2 and CSN1S2 genotypes affected milk protein content and yield. The polymorphism of the CSN2 gene affected milk protein percentage with the highest values recorded in the CSN2 AA goats (4.35, at P < 0.001). The CSN1S2 AC goats provided the highest fat (51.02 g/day) and protein (41.42 g/day) (P < 0.01) production. This information can be incorporated into selection schemes for the Sarda breed goat.

Effects of Arginine concentration on the in vitro expression of Casein and mTOR pathway related genes in mammary epithelial cells from dairy cattle

Arginine (Arg) is a conditionally-essential amino acid that is taken up by bovine mammary gland in excess of its output in milk. In this study we evaluated the effects of Arg concentration on the expression of casein and signaling pathway-related genes in mammary epithelial cells. The treatments (applied for 24 h) were designed to be devoid of Arg 0X (control; 0.00 mg/L), resemble the profile of Arg in casein (Arg 1X; 278.00 mg/L), be deficient [Arg 0.25X (69.50 mg/L) and Arg 0.5X (139.00 mg/L)], or be in excess of the amount in casein [Arg 2X (556.00 mg/L), Arg 4X (1,112 mg/L), and Arg 8X (2,224 mg/L)]. The expression of CSN1S, CSN3 and mTOR in the experimental groups was higher than those of the control group (P<0.05). Except for Arg 0.25X and Arg 8X (P>0.05), the expression of CSN1S2, CSN2 and JAK2 in other experimental groups was higher (P<0.05) than those in the control group. Except for Arg 8X (P>0.05), the expression of STAT5 in the other experimental groups was higher than those of the control (P<0.05). It also was observed that except for Arg 0.5X, the S6K expression was higher in other experimental groups than the control (P<0.05). In contrast, except for Arg 0.25X the other experimental groups resulted in lower 4EBP1 expression than the control (P<0.05). Among groups, the expression of CSN1S1, CSN1S2, CSN2, CSN3, JAK2, STAT5, mTOR and S6K gene was highest with Arg 2X (P<0.05); the reverse was true for 4EBP1 gene, with the lowest expression in this group (P<0.05). Taken together, Arg appears to play an important role in the transcriptional regulation of casein genes and mTOR-related genes in bovine mammary epithelial cells.

A novel highly divergent protein family identified from a viviparous insect by RNA-seq analysis: a potential target for tsetse fly-specific abortifacients

In tsetse flies, nutrients for intrauterine larval development are synthesized by the modified accessory gland (milk gland) and provided in mother's milk during lactation. Interference with at least two milk proteins has been shown to extend larval development and reduce fecundity. The goal of this study was to perform a comprehensive characterization of tsetse milk proteins using lactation-specific transcriptome/milk proteome analyses and to define functional role(s) for the milk proteins during lactation. Differential analysis of RNA-seq data from lactating and dry (non-lactating) females revealed enrichment of transcripts coding for protein synthesis machinery, lipid metabolism and secretory proteins during lactation. Among the genes induced during lactation were those encoding the previously identified milk proteins (milk gland proteins 1-3, transferrin and acid sphingomyelinase 1) and seven new genes (mgp4-10). The genes encoding mgp2-10 are organized on a 40 kb syntenic block in the tsetse genome, have similar exon-intron arrangements, and share regions of amino acid sequence similarity. Expression of mgp2-10 is female-specific and high during milk secretion. While knockdown of a single mgp failed to reduce fecundity, simultaneous knockdown of multiple variants reduced milk protein levels and lowered fecundity. The genomic localization, gene structure similarities, and functional redundancy of MGP2-10 suggest that they constitute a novel highly divergent protein family. Our data indicates that MGP2-10 function both as the primary amino acid resource for the developing larva and in the maintenance of milk homeostasis, similar to the function of the mammalian casein family of milk proteins. This study underscores the dynamic nature of the lactation cycle and identifies a novel family of lactation-specific proteins, unique to Glossina sp., that are essential to larval development. The specificity of MGP2-10 to tsetse and their critical role during lactation suggests that these proteins may be an excellent target for tsetse-specific population control approaches.

Amelioration of reproduction-associated oxidative stress in a viviparous insect is critical to prevent reproductive senescence

Impact of reproductive processes upon female health has yielded conflicting results; particularly in relation to the role of reproduction-associated stress. We used the viviparous tsetse fly to determine if lactation, birth and involution lead to damage from oxidative stress (OS) that impairs subsequent reproductive cycles. Tsetse females carry an intrauterine larva to full term at each pregnancy cycle, and lactate to nourish them with milk secretions produced by the accessory gland ( = milk gland) organ. Unlike most K-strategists, tsetse females lack an apparent period of reproductive senescence allowing the production of 8-10 progeny over their entire life span. In a lactating female, over 47% of the maternal transcriptome is associated with the generation of milk proteins. The resulting single larval offspring weighs as much as the mother at birth. In studying this process we noted an increase in specific antioxidant enzyme (AOE) transcripts and enzymatic activity at critical times during lactation, birth and involution in the milk gland/fat body organ and the uterus. Suppression of superoxide dismutase (sod) decreased fecundity in subsequent reproductive cycles in young mothers and nearly abolished fecundity in geriatric females. Loss of fecundity was in part due to the inability of the mother to produce adequate milk to support larval growth. Longevity was also impaired after sod knockdown. Generation of OS in virgin females through exogenous treatment with hydrogen peroxide at times corresponding to pregnancy intervals reduced survival, which was exacerbated by sod knockdown. AOE expression may prevent oxidative damage associated with the generation of nutrients by the milk gland, parturition and milk gland breakdown. Our results indicate that prevention of OS is essential for females to meet the growing nutritional demands of juveniles during pregnancy and to repair the damage that occurs at birth. This process is particularly important for females to remain fecund during the latter portion of their lifetime.