Assessment of the gene content of the chromosomal regions flanking bovine DGAT1

As a first step towards verifying the candidate status of DGAT1 as the causal gene for milk fat percentage in cattle, we constructed a bovine BAC contig spanning 576 kb of the chromosomal region containing DGAT1. High content of NotI sites (21 within the contig) indicated that the region is gene-rich. Twenty-three genes neighboring DGAT1 were mapped, including two bovine cDNA sequences that have no orthologous sequences within the NCBI sequence databases. On average, 2015 bp for each of the 23 neighboring genes were sequenced and entered into EMBL. Likewise, 10 new STS markers were established by BAC-end sequencing. Within the genes and STS markers, 55 polymorphisms were discovered. These will form the basis of future linkage disequilibrium studies to test whether any genes neighboring DGAT1 are associated with variation in milk fat percentage, thereby testing the candidate status of DGAT1.     

Effects of DGAT1 and GHR on milk yield and milk composition in the Chinese dairy population

Previous studies have demonstrated that the p.Lys232Ala substitution in the acylCoA: diacylglycerol acyltransferase ( DGAT1 ) gene and the p.Phe279Tyr mutation in the growth hormone receptor ( GHR ) gene are the causative quantitative trait loci underlying milk yield and composition on BTA14 and BTA20 respectively. To examine their applications in the genetic improvement of Chinese dairy cattle productivity, we herein investigated the effects of the DGAT1 p.Lys232Ala and GHR p.Phe279Tyr mutations on milk, fat and protein yield, as well as fat and protein percentage in the milk of 1222 Holstein cows. Genotyping was performed using PCR-RFLP for DGAT1 or primer-introduced restriction analysis (PCR-PIRA) for GHR . With a mixed animal model, the significant associations of the DGAT1 p.Lys232Ala substitution with 305-day milk, fat and protein yield were identified ( P  =   0.0001). The DGAT1 allele that encode lysine at position 232 was associated with increased 305-day milk fat yield, but with decreased 305-day milk and protein yield, whereas the GHR p.Phe279Tyr mutation was found to be significantly associated with protein percentage ( P  =   0.0014). The allele substitution effect of p.279Phe by p.279Tyr may lead to a significant increase in protein percentage. Our findings indicate that DGAT1 p.232Ala and GHR p.279Phe could be used to increase milk yield and protein yield of Chinese Holstein cows.     

The roles of segmental and tandem gene duplication in the evolution of large gene families in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Background

Most genes in Arabidopsis thaliana are members of gene families. How do the members of gene families arise, and how are gene family copy numbers maintained? Some gene families may evolve primarily through tandem duplication and high rates of birth and death in clusters, and others through infrequent polyploidy or large-scale segmental duplications and subsequent losses.

Results

Our approach to understanding the mechanisms of gene family evolution was to construct phylogenies for 50 large gene families in Arabidopsis thaliana, identify large internal segmental duplications in Arabidopsis, map gene duplications onto the segmental duplications, and use this information to identify which nodes in each phylogeny arose due to segmental or tandem duplication. Examples of six gene families exemplifying characteristic modes are described. Distributions of gene family sizes and patterns of duplication by genomic distance are also described in order to characterize patterns of local duplication and copy number for large gene families. Both gene family size and duplication by distance closely follow power-law distributions.

Conclusions

Combining information about genomic segmental duplications, gene family phylogenies, and gene positions provides a method to evaluate contributions of tandem duplication and segmental genome duplication in the generation and maintenance of gene families. These differences appear to correspond meaningfully to differences in functional roles of the members of the gene families.

     

Whole genome scan to detect quantitative trait loci for bovine milk protein composition

The objective of this study was to perform a whole genome scan to detect quantitative trait loci (QTL) for milk protein composition in 849 Holstein–Friesian cows originating from seven sires. One morning milk sample was analysed for the major milk proteins using capillary zone electrophoresis. A genetic map was constructed with 1341 single nucleotide polymorphisms, covering 2829 centimorgans (cM) and 95% of the cattle genome. The chromosomal regions most significantly related to milk protein composition ( P _genome < 0.05) were found on Bos taurus autosomes (BTA) 6, 11 and 14. The QTL on BTA6 was found at about 80 cM, and affected α_S1-casein, α_S2-casein, β-casein and κ-casein. The QTL on BTA11 was found at 124 cM, and affected β-lactoglobulin, and the QTL on BTA14 was found at 0 cM, and affected protein percentage. The proportion of phenotypic variance explained by the QTL was 3.6% for β-casein and 7.9% for κ-casein on BTA6, 28.3% for β-lactoglobulin on BTA11, and 8.6% for protein percentage on BTA14. The QTL affecting α_S2-casein on BTA6 and 17 showed a significant interaction. We investigated the extent to which the detected QTL affecting milk protein composition could be explained by known polymorphisms in β-casein , κ -casein , β-lactoglobulin and DGAT1 genes. Correction for these polymorphisms decreased the proportion of phenotypic variance explained by the QTL previously found on BTA6, 11 and 14. Thus, several significant QTL affecting milk protein composition were found, of which some QTL could partially be explained by polymorphisms in milk protein genes.     

Exploring polymorphisms and associations of the bovine MOGAT3 gene with growth traits

Monoacylglycerol acyltransferase (MGAT3, also known as MOGAT3) catalyzes the synthesis of diacylglycerol (DAG) using 2-monoacylglycerol and fatty acyl coenzyme A. This enzymatic reaction is believed to be an essential and rate-limiting step for the absorption of dietary fat in the small intestine. However, similar research for the bovine MOGAT3 gene is lacking. Therefore, in this paper, polymorphisms of the bovine MOGAT3 gene were detected in 1145 individuals from five cattle breeds by DNA pooling, PCR-RFLP, and DNA sequencing methods. The results showed that 26 novel SNPs were identified, which included 16 mutations in the coding region and the others in the introns. Additionally, association analysis between two missense mutations, g.A229G and g.G1627A, and growth traits in Nanyang cattle up to 2 years of age and adult Qinchuan cattle was performed. The results indicated that polymorphisms were significantly associated with Nanyang cattle, but no convincing associations were observed for Qinchuan cattle for the studied traits.     

Genome Wide Association Study Identifies 20 Novel Promising Genes Associated with Milk Fatty Acid Traits in Chinese Holstein

Detecting genes associated with milk fat composition could provide valuable insights into the complex genetic networks of genes underling variation in fatty acids synthesis and point towards opportunities for changing milk fat composition via selective breeding. In this study, we conducted a genome-wide association study (GWAS) for 22 milk fatty acids in 784 Chinese Holstein cows with the PLINK software. Genotypes were obtained with the Illumina BovineSNP50 Bead chip and a total of 40,604 informative, high-quality single nucleotide polymorphisms (SNPs) were used. Totally, 83 genome-wide significant SNPs and 314 suggestive significant SNPs associated with 18 milk fatty acid traits were detected. Chromosome regions that affect milk fatty acid traits were mainly observed on BTA1, 2, 5, 6, 7, 9, 13, 14, 18, 19, 20, 21, 23, 26 and 27. Of these, 146 SNPs were associated with more than one milk fatty acid trait; most of studied fatty acid traits were significant associated with multiple SNPs, especially C18:0 (105 SNPs), C18 index (93 SNPs), and C14 index (84 SNPs); Several SNPs are close to or within the DGAT1, SCD1 and FASN genes which are well-known to affect milk composition traits of dairy cattle. Combined with the previously reported QTL regions and the biological functions of the genes, 20 novel promising candidates for C10:0, C12:0, C14:0, C14:1, C14 index, C18:0, C18:1n9c, C18 index, SFA, UFA and SFA/UFA were found, which composed of HTR1B, CPM, PRKG1, MINPP1, LIPJ, LIPK, EHHADH, MOGAT1, ECHS1, STAT1, SORBS1, NFKB2, AGPAT3, CHUK, OSBPL8, PRLR, IGF1R, ACSL3, GHR and OXCT1. Our findings provide a groundwork for unraveling the key genes and causal mutations affecting milk fatty acid traits in dairy cattle.     

Detection of QTL for milk production traits in cattle by application of a specifically developed marker map of BTA6

Interval mapping was carried out to identify quantitative trait loci (QTL) for milk production traits in five granddaughter design families of the German Holstein population. Fourteen randomly generated markers spanning the whole of BTA6 and six targeted microsatellite markers from BTA6q21-31 were included in the analysis. In one family a QTL with effects on milk fat yield and milk protein yield was mapped to the interval TGLA37-FBN13 (3 CM proximal to FBN13, lodscore 3.22) in the middle part of the chromosome. Although there are several reports about QTL with effects on milk production traits on BTA6 in the literature, a QTL with effects on milk fat and milk protein yield has not been previously described.     

Identification of novel single nucleotide polymorphisms in the DGAT1 gene of buffaloes by PCR-SSCP

Diacylglycerol O-acyltransferase 1 (DGAT1) is a microsomal enzyme that catalyzes the final step of triglyceride synthesis. The DGAT1 gene is a strong functional candidate for determining milk fat content in cattle. In this work, we used PCR-SSCP (polymerase chain reaction-single-strand conformation polymorphism) and DNA sequencing to examine polymorphism in the region spanning exon 7 to exon 9 of the DGAT1 gene in Murrah and Pandharpuri buffaloes. Three alleles (A, B and C) and four novel single-nucleotide polymorphisms were identified in the buffalo DGAT1 gene. The frequencies of the alleles differed between the two buffalo breeds, with allele C being present in Murrah but not in Pandharpuri buffalo. The allele variation detected in this work may influence DGAT1 expression and function. The results described here could be useful in examining the association between the DGAT1 gene and milk traits in buffalo.     

Association of the DGAT1 gene polymorphism in stud bulls with milk productivity in cows

The substitution of lysine for alanine (K232A) in the acyl-CoA-diacylglycerol acyltransferase, which is encoded by the DGAT1 gene, was tested for the significance for breeding evaluation of stud bulls of the holsteinized Black-and-White breed. The breeding value (deviation of milk productivity in daughters compared with cows of the same age) was estimated by the DYDC (daughter yield deviation to contemporary) method with modification. The frequency of allele 232K in the bulls examined was 0.28, lower than in Holstein bulls (0.4-0.6). The greatest effect of the A232K substitution was observed for the percent fat content in milk (1.4 sigma) and milk yield (0.76 sigma), and a lower effect was established for the milk protein yield (0.47 sigma) and percent protein content in milk (0.44 sigma). In the case of milk fat yield, the effect was nonsignificant. A method was proposed for converting the data on fat yield in order to obtain significant results in this case as well. The effect of the A232K substitution was estimated at 154 kg for milk yield, 2.8 kg for protein yield, 0.079% for percent fat content, and 0.015% for percent protein content. The results are discussed in the context of multiple pleiotropic effects of the K232A substitution in the DGAT1 gene. It is proposed that the K232A substitution of the DGAT1 gene be used as a golden standard in comparisons of the effect on milk productivity for the total gene set. This approach will allow a meta-analysis of the gene effects in spite of the differences in dairy cattle breeds and methods used to analyze their breeding value. In view of more than 30-year experience of using sperm of Holstein stud bulls, including those bred in North America, it was noted that the effect of the A232K substitution on milk productivity traits agreed well with the data reported for the North American commercial population of Holstein cattle.     

Patterns of gene duplication in the plant SKP1 gene family in angiosperms: evidence for multiple mechanisms of rapid gene birth

Gene duplication plays important roles in organismal evolution, because duplicate genes provide raw materials for the evolution of mechanisms controlling physiological and/or morphological novelties. Gene duplication can occur via several mechanisms, including segmental duplication, tandem duplication and retroposition. Although segmental and tandem duplications have been found to be important for the expansion of a number of multigene families, the contribution of retroposition is not clear. Here we show that plant SKP1 genes have evolved by multiple duplication events from a single ancestral copy in the most recent common ancestor (MRCA) of eudicots and monocots, resulting in 19 ASK (Arabidopsis SKP1-like) and 28 OSK (Oryza SKP1-like) genes. The estimated birth rates are more than ten times the average rate of gene duplication, and are even higher than that of other rapidly duplicating plant genes, such as type I MADS box genes, R genes, and genes encoding receptor-like kinases. Further analyses suggest that a relatively large proportion of the duplication events may be explained by tandem duplication, but few, if any, are likely to be due to segmental duplication. In addition, by mapping the gain/loss of a specific intron on gene phylogenies, and by searching for the features that characterize retrogenes/retrosequences, we show that retroposition is an important mechanism for expansion of the plant SKP1 gene family. Specifically, we propose that two and three ancient retroposition events occurred in lineages leading to Arabidopsis and rice, respectively, followed by repeated tandem duplications and chromosome rearrangements. Our study represents a thorough investigation showing that retroposition can play an important role in the evolution of a plant gene family whose members do not encode mobile elements.     

DGAT1 K232A quantitative trait nucleotide polymorphism in Polish Black-and-White cattle

The diacylglycerol o-acyltransferase 1 gene (DGAT1) was investigated in Polish Black-and-White cattle. The frequency of the K allele was 0.60, 0.68 and 0.48 for AI sires (n=150), young bulls (n=139) and cows (n=213), respectively. The method of selective genotyping for identification of the quantitative trait nucleotide was verified through identification of DGAT1 effect on milk production traits. Daughters of six heterozygous bulls were selectively genotyped based on their milk traits. The genotypic frequencies differed between high and low yield groups representing milk and fat contents. The Kruskal-Wallis test revealed a highly significant effect of DGAT1 K232A in cows with extremely low fat content and a significant effect in cows with extremely high protein content of milk. No significant effect of AI sires' genotypes on their breeding value was found.     

FISH-mapping of LEP and SLC26A2 genes in sheep, goat and cattle R-banded chromosomes: comparison between bovine, ovine and caprine chromosome 4 (BTA4/OAR4/CHI4) and human chromosome 7 (HSA7)

Sheep (OAR), goat (CHI) and cattle (BTA) R-banded chromosome preparations, obtained from synchronized cell cultures, were used to FISH-map leptin (LEP) and solute carrier family 26 member 2 (SLC26A2) genes on single chromosome bands. LEP maps on OAR4q32 and CHI4q32, being the first assignment of this gene to these two species. SLC26A2 maps on BTA7q24, OAR5q24 and CHI7q24. This gene, too, was assigned for the fist time to both sheep and goat chromosomes, while it was more precisely localized on a single chromosome band in cattle. Improved cytogenetic maps of BTA4/OAR4/CHI4 were constructed and compared with HSA7 revealing five main conserved segments and complex chromosome rearrangements, including a centromere repositioning, differentiating HSA7 and BTA4/OAR4/CHI4.     

Revisiting the quantitative trait loci for milk production traits on BTA6

A parallel association study was performed in two independent cattle populations based on 41 validated, targeted single nucleotide polymorphisms (SNPs) and four microsatellite markers to re-evaluate the multiple quantitative trait loci (QTL) architecture for milk performance on bovine chromosome 6 (BTA6). Two distinct QTL located in the vicinity of the middle region of BTA6, but differing unambiguously regarding their effects on milk composition and yield traits were validated in the German Holstein population. A highly significant association of the protein variant ABCG2 p.Tyr581Ser with milk composition traits reconfirmed the causative molecular relevance of the ABCG2 gene in QTL region 1, whereas in QTL region 2, significant and tentative associations between gene variants RW070 and RW023 (located in the promoter region and exon 9 of the PPARGC1A gene for milk yield traits) were detected. For the German Fleckvieh population, only RW023 showed a tentative association with milk yield traits, whereas those loci with significant effects in German Holsteins (ABCG2 p.Tyr581Ser, RW070) showed fixed alleles. Even though our new data highlight two variants in the PPARGC1A gene (RW023, RW070) in QTL region 2, based on the results of our study, currently no unequivocal conclusion about the causal background of this QTL affecting milk yield traits can be drawn. Notably, the German Holstein and Fleckvieh populations, known for their divergent degree of dairy type, differ substantially in their allele frequencies for the growth-associated NCAPG p.Ile442Met locus.     

Association of the <Emphasis Type="Italic">DGAT1</Emphasis> gene polymorphism in bulls with cow milk performance

The substitution of lysine for alanine (K232A) in the acyl-CoA:diacylglycerol acyltransferase, which is encoded by the DGAT1 gene, was tested for the significance for breeding evaluation of bulls of the holsteinized Black-and-White breed. The breeding value was estimated by the DYDC (daughter yield deviation to contemporaries) method with modification. The frequency of allele 232K in the bulls examined was 0.28, lower than in Holstein bulls (0.4–0.6). The greatest effect of the A232K substitution was observed for the percent fat milk yield (1.4σ) and milk yield (0.76σ), and a lower effect was established for the milk protein yield (0.47σ) and percent protein milk yield (0.44σ). In the case of milk fat yield, the effect was nonsignificant. A method was proposed for converting the data on fat yield in order to obtain significant results in this case as well. The effect of the A232K substitution was estimated at 154 kg for milk yield, 2.8 kg for milk protein yield, 0.079% for percent milk fat, and 0.015% for percent milk protein yield. The results are discussed in the context of multiple pleiotropic effects of the K232A substitution in the DGAT1 gene. It is proposed that the K232A substitution of the DGAT1 gene may be used as a golden standard in comparisons of the effect on milk production traits for the total gene set. This approach will allow a meta-analysis of the gene effects in spite of the different dairy cattle breeds and methods used to analyze their breeding value. In view of more than 30-year experience of using sperm of Holstein bulls, including those breeds in North America, it was noted that the effect of the A232K substitution on milk production traits agreed well with the data reported for the North American commercial population of Holstein cattle.     

Comparative map alignment of BTA27 and HSA4 and 8 to identify conserved segments of genome containing fat deposition QTL

Quantitative trait loci (QTL) associated with fat deposition have been identified on bovine Chromosome 27 (BTA27) in two different cattle populations. To generate more informative markers for verification and refinement of these QTL-containing intervals, we initiated construction of a BTA27 comparative map. Fourteen genes were selected for mapping based on previously identified regions of conservation between the cattle and human genomes. Markers were developed from the bovine orthologs of genes found on human Chromosomes 1 (HSA1), 4, 8, and 14. Twelve genes were mapped on the bovine linkage map by using markers associated with single nucleotide polymorphisms or microsatellites. Seven of these genes were also anchored to the physical map by assignment of fluorescence in situ hybridization probes. The remaining two genes not associated with an identifiable polymorphism were assigned only to the physical map. In all, seven genes were mapped to BTA27. Map information generated from the other seven genes not syntenic with BTA27 refined the breakpoint locations of conserved segments between species and revealed three areas of disagreement with the previous comparative map. Consequently, portions of HSA1 and 14 are not conserved on BTA27, and a previously undefined conserved segment corresponding to HSA8p22 was identified near the pericentromeric region of BTA8. These results show that BTA27 contains two conserved segments corresponding to HSA8p, which are separated by a segment corresponding to HSA4q. Comparative map alignment strongly suggests the conserved segment orthologous to HSA8p21-q11 contains QTL for fat deposition in cattle. Received: 25 February 2000 / Accepted: 30 March 2000     

DGAT1, a new positional and functional candidate gene for intramuscular fat deposition in cattle

Intramuscular fat content, also assessed as marbling of meat, represents an important beef quality trait. Recent work has mapped a quantitative trait locus (QTL) with an effect on marbling to the centromeric region of bovine chromosome 14, with the gene encoding thyroglobulin (TG) being proposed as a positional and functional candidate gene for this QTL. Recently, the gene encoding diacylglycerol O-acyltransferase (DGAT1), which also has been mapped within the region of the marbling QTL, has been demonstrated to affect the fat content of milk. In the present study, the effects of a 5'-polymorphism of TG and of a lysine/alanine polymorphism of DGAT1 on the fat content of musculus (m.) semitendinosus and m. longissimus dorsi in 55 bovine animals (28 German Holstein and 27 Charolais) has been investigated. Significant effects were found for both candidate genes in both the breeds. These effects seem to be independent of one another because the alleles of the two polymorphisms showed no statistically significant disequilibrium. The DGAT1 effect is mainly on the m. semitendinosus. The TG polymorphism only affects m. longissimus dorsi. However, both intramuscular fat enhancing effects seem to be recessive. The possibility of two linked loci, acting recessively on intramuscular fat content, will require special strategies when selecting for higher marbling scores.     

A whole genome scan for QTL affecting milk protein percentage in Italian Holstein cattle, applying selective milk DNA pooling and multiple marker mapping in a daughter design

We report on a complete genome scan for quantitative trait loci (QTL) affecting milk protein percentage (PP) in the Italian Holstein-Friesian cattle population, applying a selective DNA pooling strategy in a daughter design. Ten Holstein-Friesian sires were chosen, and for each sire, about 200 daughters, each from the high and low tails of estimated breeding value for PP, were used to construct milk DNA pools. Sires and pools were genotyped for 181 dinucleotide microsatellites covering all cattle autosomes. Sire marker allele frequencies in the pools were obtained by shadow correction of peak height in the electropherograms. After quality control, pool data from eight sires were used for all subsequent analyses. The QTL heterozygosity estimate was lower than that of similar studies in other cattle populations. Multiple marker mapping identified 19 QTL located on 14 chromosomes (BTA1, 2, 3, 4, 5, 6, 8, 9, 12, 14, 17, 20, 23 and 27). The sires were also genotyped for seven polymorphic sites in six candidate genes (ABCG2, SPP1, casein kappa, DGAT1, GHR and PRLR) located within QTL regions of BTA6, 14 and 20 found in this study. The results confirmed or excluded the involvement of some of the analysed markers as the causative polymorphic sites of the identified QTL. The QTL identified, combined with genotype data of these candidate genes, will help to identify other quantitative trait genes and clarify the complex QTL patterns observed for a few chromosomes. Overall, the results are consistent with the Italian Holstein population having been under long-term selection for high PP.     

Comparative FISH-mapping of the survival of motor neuron gene (SMN) in domestic bovids

A comparative fluorescence in situ mapping of the SMN gene was performed on R-banded chromosome preparations of cattle (Bos taurus, BTA, 2n = 60), river buffalo (Bubalus bubalis, BBU, 2n = 50), sheep (Ovis aries, OAR, 2n = 54) and goat (Capra hircus, CHI, 2n = 60), as well as on those of a calf from Piedmont breed affected by arthrogryposis. SMN was located on BTA20q13.1, OAR16q13.1, CHI20q13.1 and BBU19q13. These chromosomes and chromosome bands are believed to be homeologous, confirming the high degree of chromosome homeologies among bovids. The position of SMN was refined in cattle, compared to the two previous localizations, while it is a new gene assignment in the other three bovids. A comparative fiber-FISH performed on extended chromatin of both normal cattle and calf affected by arthrogryposis revealed more extended FITC signals in the calf, compared to the normal cattle (control), suggesting a possible duplication of the SMN gene in the calf affected by arthrogryposis. .     

Molecular evolution of glycinin and β-conglycinin gene families in soybean (Glycine max L. Merr.)

There are two main classes of multi-subunit seed storage proteins, glycinin (11S) and β-conglycinin (7S), which account for approximately 70% of the total protein in a typical soybean seed. The subunits of these two protein classes are encoded by a number of genes. The genomic organization of these genes follows a complex evolutionary history. This research was designed to describe the origin and maintenance of genes in each of these gene families by analyzing the synteny, phylogenies, selection pressure and duplications of the genes in each gene family. The ancestral glycinin gene initially experienced a tandem duplication event; then, the genome underwent two subsequent rounds of whole-genome duplication, thereby resulting in duplication of the glycinin genes, and finally a tandem duplication likely gave rise to the Gy1 and Gy2 genes. The β-conglycinin genes primarily originated through the more recent whole-genome duplication and several tandem duplications. Purifying selection has had a key role in the maintenance of genes in both gene families. In addition, positive selection in the glycinin genes and a large deletion in a β-conglycinin exon contribute to the diversity of the duplicate genes. In summary, our results suggest that the duplicated genes in both gene families prefer to retain similar function throughout evolution and therefore may contribute to phenotypic robustness.