Mitochondrial DNA single nucleotide polymorphism associated with weight estimated breeding values in Nelore cattle (Bos indicus)

ANIMAL GENETICS

RESEARCH ARTICLE

Mitochondrial DNA single nucleotide polymorphism associated with weight estimated breeding values in Nelore cattle (Bos indicus)

Fernando Henrique Biase^I; Flávio Vieira Meirelles^II; Ricardo Gunski^III; Pedro Alejandro Vozzi^I; Luiz A.F. Bezerra^I; Reginaldo A. Vila^I; Artur J.M. Rosa^IV; Raysildo B. Lôbo^I; Lúcia Martelli^I

^IDepartamento de Genética, Faculdade de Medicina de Ribeirão Preto, Universidade de São Paulo, Ribeirão Preto, São Paulo, Brazil

^IIDepartamento de Ciências Básicas, Faculdade de Zootecnia e Engenharia de Alimentos, Universidade de São Paulo, Pirassununga, São Paulo, Brazil

^IIICentro de Ciências Rurais de São Gabriel, Universidade Federal do Pampa, São Gabriel, RS, Brazil

^IVDepartment of Animal and Range Sciences, South Dakota State University, Brookings, South Dakota, United States of America

^{Send correspondence to}

ABSTRACT

We sampled 119 Nelore cattle (Bos indicus), 69 harboring B. indicus mtDNA plus 50 carrying Bos taurus mtDNA, to estimate the frequencies of putative mtDNA single nucleotide polymorphisms (SNPs) and investigate their association with Nelore weight and scrotal circumference estimated breeding values (EBVs). The PCR restriction fragment length polymorphism (PCR-RFLP) method was used to detect polymorphisms in the mitochondrial asparagine, cysteine, glycine, leucine and proline transporter RNA (tRNA) genes (tRNA^asn, tRNA^cys, tRNA^gly, tRNA^leu and tRNA^pro). The 50 cattle carrying B. taurus mtDNA were monomorphic for all the tRNA gene SNPs analyzed, suggesting that they are specific to mtDNA from B. indicus cattle. No tRNA^cys or tRNA^gly polymorphisms were detected in any of the cattle but we did detect polymorphic SNPs in the tRNA^asn, tRNA^leu and tRNA^pro genes in the cattle harboring B. indicus mtDNA, with the same allele observed in the B. taurus sequence being present in the following percentage of cattle harboring B. indicus mtDNA: 72.46% for tRNA^asn, 95.23% for tRNA^leu and 90.62% for tRNA^pro. Analyses of variance using the tRNA^asn SNP as the independent variable and EBVs as the dependent variable showed that the G ® T SNP was significantly associated (p < 0.05) with maternal EBVs for weight at 120 and 210 days (p < 0.05) and animal's EBVs for weight at 210, 365 and 455 days. There was no association of the tRNA^asn SNP with the scrotal circumference EBVs. These results confirm that mtDNA can affect weight and that mtDNA polymorphisms can be a source of genetic variation for quantitative traits.

Key words: bovine, mitochondria, mtDNA, SNP, weight.

Introduction

Mitochondria are eukaryotic cell organelles involved in various cellular functions, including cell proliferation, apoptosis and, mostly important, energy production (Birch-Machin, 2006) by oxidative phosphorylation (Taanman, 1999). These organelles are responsible for approximately 90% of the energy produced by the mammal cell (Boettcher et al., 1996b). Mitochondria have their own genome, which is maternally inherited in mammals and is an important source of cytoplasmic genetic variation (Gibson et al., 1997; Birky, 2001). In cattle, direct maternal effects have been discussed since Wagner (1972) but the contribution of cytoplasmic effects estimated for growth and milk production are not consistent between populations (Gibson et al., 1997).

Boettcher et al. (1996b) demonstrated bias in the heritability, permanent environmental variance, and accuracy estimation in the animal model when the cytoplasmic effects were ignored. Moreover, the effect of cytoplasmic inheritance has been studied in respect to milk production (Bell et al., 1985; Tess et al., 1987; Schutz et al., 1992; Boettcher et al., 1996a) and beef cattle growth traits (Tess and Robison, 1990; Northcutt et al., 1991; Tess and MacNeil, 1994; Quintanilla et al., 1999).

Recently, molecular biology has furthered our understanding of the function and inheritance of the mitochondrial genome and its importance on livestock development and production (Smith and Alcivar, 1993; Smith et al., 2000). Since Anderson et al. (1982) published the complete mitochondrial DNA (mtDNA) sequence, nucleotide variants in coding and non-coding regions have been studied to associate molecular markers with the production of Bos taurus breeds (Ron et al., 1993; Suzuki et al., 1993; Mannen et al., 1998; Mannen et al., 2003). Pegoraro et al. (1996) identified sequence alterations in mitochondrial transfer RNA (tRNA) genes and in the origin of light strand (Ori L) replication of Bos indicus Nelore cattle as compared to B. taurus and described a single nucleotide polymorphism in the Nelore Ori L region and asparagine tRNA (tRNA^asn) gene. The polymorphisms were later confirmed by the complete genomic sequence of the B. indicus mitochondrial DNA (mtDNA) and deposited in the GenBank under accession number AY126697.

Meirelles et al. (1999) reported the existence of cattle carrying the B. taurus or B. indicus mitochondrial genome in Brazilian Nelore cattle. These animals, with same breed composition but harboring different mitochondrial genomes, constitute an interesting model to investigate the influence of mtDNA polymorphisms on quantitative traits. During the work described in the present paper we tested the hypothesis that cattle with different specific nucleotide variations differ in estimated breeding values (EBVs) for growth and reproductive traits. The aims of this study were to estimate the frequency of single nucleotide polymorphisms of the mitochondrial genes for asparagine tRNA (tRNA^asn), cysteine tRNA (tRNA^cys) (Pegoraro et al., 1996), glycine tRNA (tRNA^gly), leucine tRNA (tRNA^leu) and proline tRNA (tRNA^pro) in B. indicus and B. taurus mtDNA from public GenBank sequences AY126697 and AY526085 in Brazilian Nelore cattle and to investigate the effects of these polymorphisms on weight and scrotal circumference estimated breeding values at different ages.

Material and Methods

Cattle sample

The sample (n = 119) consisted of 69 purebred adult Nelore cattle registered as pure origin imported (POI) in the Brazilian Association of Zebu Breeder's Herdbook and identified as carriers of B. indicus mtDNA plus a further 50 purebred B. indicus Nelore cattle which were carriers of B. taurus mtDNA and not registered as POI. A blood sample was collected from each animal and total DNA extracted according to Sambrook et al. (2001). The presence of B. indicus and B. taurus mtDNA was confirmed by the polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP) method by the absence (B. indicus) or the presence (B. taurus) of HindIII restriction site at nucleotide 12178 on the NADH dehydrogenase subunit 5 gene (Meirelles et al., 1999).

Detection of polymorphisms

The single nucleotide substitutions (SNP) in five transfer RNA (tRNA) mitochondrial genes were detected by PCR-RFLP using endonuclease enzymatic digestion of the amplicons. Primers were designed using the B. indicus mtDNA sequence as a reference (GenBank AY126697). The PCR was performed using 150 ng of total DNA, 1x PCR buffer (20 mM Tris-HCl, 50 mM KCl, Invitrogen, Brazil), 3 mM MgCl₂, 0.2 mM of each dNTP, 0.2 µM of each primer (Table 1) and 1.5 units of Taq DNA polymerase (Invitrogen, Brazil) in a final volume of 50 µL. Cycling was set for 5 min at 95 °C followed by 35 cycles of 40 s at 95 °C, 30 s at the specific temperature of each primer and 40 s at 72 °C, with a final 72 °C extension for 5 min. The PCR products (15 µL) were digested with 1 unit of Bsr I, Mnl I, Spe I, Dra I or Bsm AI restriction endonuclease specific for each tRNA gene sequence for one hour at the temperature recommended by the supplier (New England Biolabs, USA). The mitochondrial tRNA PCR amplicons, specific primer sequences, annealing temperatures, SNP positions and restriction enzymes used in this study are presented in Table 1. Electrophoresis was carried out on 2% (w/v) agarose gel and DNA was stained with ethidium bromide and the fragments visualized under ultraviolet light. Mitochondrial transfer RNA SNP frequencies were obtained by direct counting of the PCR-RFLP results.

Estimated breeding values

The cattle selected for mtDNA genotyping were part of the Nelore cattle breeding program (NCBP) at the University of São Paulo (ANCP, http://www.ancp.org.br/ sumarios) which estimated breeding values for approximately 900,000 cattle. The estimated breeding values used in this study for association analysis were kindly provided by the NCBP.

The traits analyzed were maternal effects for weight (MW) at 120 and 210 days (MW₁₂₀ and MW₂₁₀), animal weight (AW) at 120, 210, 365, and 450 days (AW₁₂₀, AW₂₁₀, AW₃₆₅, and AW₄₅₀) and scrotal circumference (SC) at 365, 450, and 550 days (SC₃₆₅, SC₄₅₀, and SC₅₅₀).

Briefly, breeding values for each animal were estimated using an animal model applying the best linear unbiased predictor (Boldman et al., 1995). The model used for genetic parameter estimations was y = Xb + Z1a + Z2m + Z3p + e, where y is a N x 1 vector of records, b the vector of fixed effects (herd-year-season, sex, dam age class), X is the matrix associating b with y, a is the vector of estimated breeding values for direct genetic effects, Z1 is the matrix associating a with y, m is the vector of estimated breeding values for maternal genetic effects, Z2 is the matrix associating m with y, p is the vector of permanent environmental non-additive genetic effects contributed by dams to the records of their progeny, Z3 is the matrix associating p with y and e is the vector of residual effects (Gunski et al., 2001).

Association analysis

Maternal and animal estimated breeding values were independently compared between animals with and without the Bsr I restriction site (G or T nucleotide) on the tRNA^asn gene by analysis of variance (ANOVA). The linear model used was: y_ij = µ+N_i+e_ij, where y_ij is the breeding value for each trait of the ij^th animal, N is the fixed effect of the i^th tRNA^asn nucleotide (G or T) and e_ij is the random error effect associated with the ij^th observation. Means were compared using the Tukey Studentized Range Test in the SAS program (SAS Institute Inc., 2001). Differences with a < 0.05 were considered statistically significant.

Nelore Heard-book genealogy data incorporated in the breeding value estimations allowed us to obtain the estimated breeding values for all parents and relatives from each genotyped animal. However, only those ascending the animals harboring B. indicus mtDNA were considered for further analysis. Because mitochondrial DNA is inherited maternally, females and their offspring from the same matrilineal lineage were assumed to have the same mitochondrial tRNA^asn nucleotide at position 5501. This allowed us to group females from the same lineage with previously genotyped animals. The new data set was based on 345 observations for each characteristic analyzed and was used for ANOVA and the Tukey test.

Results and Discussion

All the cattle genotyped showed no sign of heteroplasmy for the five single nucleotide polymorphisms (SNPs) analyzed, suggesting maternal inheritance and homoplasmic distribution within tissues (Attardi, 1985; Smith et al., 2000;Birky, 2001). The 50 cattle genotyped as carrying B. taurus mtDNA were monomorphic for the five mitochondrial tRNA gene SNPs analyzed in this study, suggesting that the SNPs evaluated are specific to mtDNA from B. indicus cattle.

In this Nelore sample (n = 119) no mitochondrial tRNA^cys or tRNA^gly mutations were found using PCR-RFLP. However, we did detect polymorphic SNPs in the tRNA^asn, tRNA^leu and tRNA^pro mitochondrial genes of the cattle harboring B. indicus mtDNA (n = 69) where the B. taurus sequence (Anderson et al., 1982) was present in 72.46% of animals for tRNA^asn, 95.23% for tRNA^leu and 90.62% for tRNA^pro, the remaining animals having a variant tRNA pattern similar to the B. indicus mtDNA GenBank sequence AY126697. Tracing the matrilineal genealogy, we found that each of the 19 cattle genotyped as having B. indicus tRNA^asn were descendants of different females imported from India, suggesting that the nucleotide mutations in the tRNA^asn, tRNA^leu and tRNA^pro mitochondrial genes, characteristic of B. indicus cattle, occurred in India and were brought to Brazil with imported Indian cows.

Because the mitochondrial tRNA^cys and tRNA^gly genes were not polymorphic for the SNPs tested in our sample and the mitochondrial tRNA^leu, and tRNA^pro genes lacked sufficient variability, we used the SNP on the tRNA^asn gene to compare the estimated breeding values of the cattle harboring B. indicus mtDNA in our sample.

As stated above, maternal and animal estimated breeding values were independently compared between the cattle in our sample which had, or lacked, the Bsr I restriction site (G or T nucleotide). The guanidine to thymine substitution was significantly associated (p < 0.05) with changes in maternal breeding value for weight at 120 and 210 days and the animal breeding value for weight at 210, 365 and 455 days (Table 2). Cattle with the mitochondrial tRNA^asn guanidine SNP had higher estimated breeding values compared to those with the thymine SNP (Table 3), although there were no differences for the scrotal circumference estimated breeding values between cattle with different nucleotides at this position (Tables 4 and 5).

Gunski et al. (2001) studied mitochondrial genome effects to compare pre-weaning growth traits between were in Nelore cattle harboring B. taurus mtDNA and B. indicus mtDNA, but found no significant (p > 0.05) differences for maternal and animal estimated breeding values.

Our results suggest an association between polymorphism in the mitochondrial tRNA^asn gene with growth traits in B. indicus Nelore cattle. The differences for maternal and direct weight estimated breeding values were not caused by a B. taurus nuclear genome effect since this sample were known to harbor the B. indicus mtDNA and were progeny of registered POI sires and cows.

There are reports associating cytoplasmic lineage and growth and milk traits effects in various cattle breeds. Tess et al. (1987) reported that cytoplasmic effects influenced pre-weaning and milk traits in Hereford cattle but later results did not confirm this effect (Tess and Robison, 1990), while Schutz et al., (1993) found that maternal lineage was responsible for 4% to 15% of the phenotypic variation in milk traits of Holstein-Friesian cattle. Boettcher et al., (1996b) reported that the inclusion of the maternal lineage effect in a simulated animal model gave better estimations of heritability and permanent environment variance and also accurate estimated breeding values. However, due to small number of different lineages in our data set we could not test if cytoplasmic lineage was associated with breeding value variation.

Quantitative traits have also been associated with mtDNA polymorphisms. D-loop variation has been associated to carcass traits such as longissimus muscle area and beef marbling score (Mannen et al., 1998), milk production as measured by yield, fat content and estimated milk energy (Schutz et al., 1994), and calving rates (Sutarno et al., 2002). Furthermore, polymorphic sites on rRNA genes have also been associated with milk traits (Boettcher et al., 1996a).

The effect of variation in the mitochondrial genome on milk and beef traits is thought to be a consequence of changes in metabolic rate or the energy available for milk production and muscle development, the differences reported by us in this study tend to support this hypothesis.

We investigated scrotal circumference (SC) because this has been correlated with reproductive traits in Nelore cattle (Martins-Filho and Lôbo, 1991) and the heritability of scrotal circumference is considered moderate to high at 0.36 to 0.47 (Garnero et al., 2001; Pereira et al., 2002; Silveira et al., 2004). However, no studies have suggested that the mitochondrial genome may influence fertility traits and the polymorphism tested by us was not associated with scrotal circumference estimated breeding values.

In humans, there are approximately 108 disorders that are associated with SNPs of mitochondrial tRNA genes. Three of these disorders involve the tRNA^asn gene (Brandon et al., 2005), indicating that SNPs of mitochondrial tRNA genes have an important role in mitochondria. Mitochondrial tRNA^asn gene SNP (G5501 ® T) could affect maternal and direct weight estimated breeding values due to a direct effect on the functional efficiency of the tRNA^asn gene. Otherwise, the tRNA^asn gene SNP genotyped in our study may be associated to other functional polymorphisms or a group of polymorphisms (haplotype). While our results are not conclusive concerning the effect of the tRNA^asn gene SNP on Nelore estimated breeding values, they provide important evidence that nucleotide variation within the mitochondrial genome should be considered as genetic markers for the assisted selection of quantitative traits.

The B. taurus and B. indicus mitochondrial genomes have been completely sequenced but the existing polymorphisms between these two breeds need to be studied to further characterize the existing genetic variability in these breeds and elucidate the effects of such polymorphisms on production traits. The association of weight estimated breeding values with a tRNA^asn gene SNP confirms that mitochondrial DNA is a source of genetic variation that influences growth during the bovine pre-weaning and post-weaning periods.

Acknowledgments

The authors would like to thank the Brazilian agencies CAPES, CNPq-PRONEX and FAPESP (proc.: 98/16573-8) for financial support and also the Breeders belonging to the Nelore Cattle Breeding Program – USP for providing blood samples from the cattle used in this study.

Internet Resource

Nelore cattle breeding program (NCBP) at the University of São Paulo, under the auspices of the Brazilian National Association of Breeders and Researchers of Nelore Cattle Associação Nacional de Criadores e Pesquisadores – ANCP) http://www.ancp.org.br/sumarios.

Received: August 15, 2006; Accepted: March 8, 2007.

Associate Editor: Luiz Lehmann Coutinho

• Anderson S, Bruijn MHL, Coulson AR, Eperon IC, Sanger F and Young IG (1982) Complete sequence of bovine mitochondrial DNA: Conserved features of the mammalian mitochondrial genome. J Mol Biol 156:683-717.
• Attardi G (1985) Animal mitochondrial DNA: An extreme example of genetic economy. Int Rev Cytol 93:93-145.
• Bell BR, McDaniel BT and Robison OW (1985) Effects of cytoplasmic inheritance on production traits of dairy cattle. J Dairy Sci 68:2038-2051.
• Birch-Machin MA (2006) The role of mitochondria in ageing and carcinogenesis. Clin Exp Dermatol 31:548-552.
• Birky CW (2001) The inheritance of genes in mitochondria and chloroplasts: Laws, mechanisms, and models. Annu Rev Genet 35:125-148.
• Boettcher PJ, Freeman AE, Johnston SD, Smith RK, Beitz DC and McDaniel BT (1996a) Relationships between polymorphism for mitochondrial deoxyribonucleic acid and yield traits of Holstein cows. J Dairy Sci 79:647-654.
• Boettcher PJ, Kuhn MT and Freeman AE (1996b) Impacts of cytoplasmic inheritance on genetic evaluations. J Dairy Sci 79:663-675.
• Boldman K, Van Vleck L, Kriese L and Kachman S (1995) MTDFREML: User's Guide. Departament of Agriculture/Agricultural Research Service, Lincon.
• Brandon MC, Lot MT, Nguyen KC, Spolim S, Navathe SB, Baldi P and Wallace DC (2005) MITOMAP: A human mitochondrial genome database-2004 update. Nucleic Acids Res 33:D611-613.
• Garnero AV, Lôbo RB, Bezzerra LAF and Oliveira HN (2001) Comparação de alguns critérios de seleção para crescimento na raça nelore. Rev Bras Zootec 30:714-718. Abstract in English.
• Gibson JP, Freeman AE and Boettcher PJ (1997) Cytoplasmic and mitochondrial inheritance of economic traits in cattle. Livest Prod Sci 47:115-124.
• Gunski RJ, Bezzerra LAF, Meirelles FV and Lôbo RB (2001) Efeito do DNA mitocondrial sobre características de crescimento pré-desmama na raça Nelore. In: Anais da 38Ş Reunião Anual da Sociedade Brasileira de Zootecnia, Piracicaba, pp 23-26. Abstract in English.
• Mannen H, Kojima T, Oyama K, Mukai F, Ishida T and Tsuji, S (1998) Effect of mitochondrial DNA variation on carcass traits of Japanese Black cattle. J Anim Sci 76:36-41.
• Mannen H, Morimoto ML, Oyamat K, Mukai F and Tsuji S (2003) Identification of mitochondrial DNA substitutions related to meat quality in Japanese Black cattle. J Anim Sci 81:68-73.
• Martins-Filho R and Lôbo RB (1991) Circunferência escrotal como critério de seleção em bovinos tendo em vista o melhoramento de características reprodutivas em ambos os sexos. Ciênc Animal 1:83-91. Abstract in English.
• Meirelles FV, Rosa AJM, Lôbo RB, Garcia JM, Smith LC and Duarte FM (1999) Is the American Zebu really Bos indicus? Genet Mol Biol 22:543-546.
• Northcutt SL, Willham RL and Wilson DE (1991) Genetic parameters for nuclear and non-nuclear inheritance in three synthetic lines of beef cattle differing in mature size. J Anim Sci 69:4745-4753.
• Pegoraro L, Yang Z, Samake S, Meirelles FV, Bordingnon V, Moquin LC and Smith LC (1996) Sequence comparison of mitochondrial tRNA genes and origin of light strand replication in Bos taurus and Nelore (Bos indicus) breeds. Anim Genet 27:91-94.
• Pereira E, Eler JP and Ferraz JP (2002) Análise genética de características reprodutivas na raça Nelore. Pesq Agropec Bras 37:703-708. Abstract in English.
• Quintanilla R, Varona L, Pujol MR and Piedrafita J (1999) Maternal animal model with correlation between maternal environmental effects of related dams. J Anim Sci 77:2904-2917.
• Ron M, Yoffe O and Weller JI (1993) Sequence variation in D-loop mtDNA of cow lineages selected for high and low maternal effects on milk production. Anim Genet 24:183-186.
• Sambrook J, Fritch EF and Maniatis T (2001) Molecular Cloning: A Laboratory Manual. 3 edition. Cold Spring Harbor Laboratory Press, New York.
• SAS Institute Inc. (2001) Statistical Analysis Software. 8.2 edition. SAS Institute Inc., Cary.
• Schutz MM, Freeman AE, Beitz DC and Mayfield JE (1992) The importance of maternal lineage on milk yield traits of dairy cattle. J Dairy Sci 75:1331-1341.
• Schutz MM, Freeman AE, Lindberg GL and Beitz DC (1993) Effects of maternal lineages grouped by mitochondrial genotypes on milk yield and composition. J Dairy Sci 76:621-629.
• Schutz MM, Freeman AE, Lindberg GL, Koehler CM and Beitz DC (1994) The effect of mitochondrial DNA on milk production and health of dairy cattle. Livestock Production Science 37:283-295.
• Silveira JC, McManus C, Mascioli AS, Silva LOC, Silveira AC, Garcia JAS and Louvandini H (2004) Study of genetic and environmental factors on production and reproduction traits in a Nelore herd in Mato Grosso do Sul state. Rev Bras Zootec 33:1432-1444.
• Smith LC and Alcivar AA (1993) Cytoplasmic inheritance and its effects on development and performance. J Reprod Fertil Suppl 48:31-43.
• Smith LC, Bordignon V, Garcia JM and Meirelles FV (2000) Mitochondrial genotype segregation and effects during mammalian development: Applications to biotechnology. Theriogenology 53:35-46.
• Sutarno CG, Cummins JM, Greeff J and Lymbery AJ (2002) Mitochondrial DNA polymorphisms and fertility in beef cattle. Theriogenology 57:1603-1610.
• Suzuki R, Kemp SJ and Teale AJ (1993) Polymerase chain reaction analysis of mitochondrial DNA polymorphism in N'Dama and Zebu cattle. Anim Genet 24:339-343.
• Taanman JW (1999) The mitochondrial genome: Structure, transcription, translation and replication. Biochim Biophys Acta 1410:103-123.
• Tess MW and MacNeil MD (1994) Evaluation of cytoplasmic genetic effects in Miles City Line 1 Hereford cattle. J Anim Sci 72:851-856.
• Tess MW, Reodecha C and Robison OW (1987) Cytoplasmic genetic effects on preweaning growth and milk yield in Hereford cattle. J Anim Sci 65:675-684.
• Tess MW and Robison OW (1990) Evaluation of cytoplasmic genetic effects in beef cattle using an animal model. J Anim Sci 68:1899-1909.
• Wagner RP (1972) The role of maternal effects in animal breeding II: Mitochondria and animal inheritance. J Anim Sci 35:1280-1287.

Publication Dates

• Publication in this collection
  13 Dec 2007
• Date of issue
  2007

History

• Received
  15 Aug 2006
• Accepted
  08 Mar 2007