Skip to main content

Mitochondrial DNA T7719G in tRNA-Lys gene affects litter size in Small-tailed Han sheep



In farm animals, mitochondrial DNA (mtDNA) effect on economic performance remains hot-topic for breeding and genetic selection. Here, 53 maternal lineages of Small-tailed Han sheep were used to investigate the association of mitochondrial DNA variations and the lambing litter size.


Sequence sweeping of the mitochondrial coding regions discovered 31 non-synonymous mutations, and the association study revealed that T7719G in mtDNA tRNA-Lys gene was associated with litter size (P < 0.05), manifesting 0.29 lambs per litter between the G and T carriers. Furthermore, using the mixed linear model, we assayed the potential association of the ovine litter size and haplogroups and multiple-level mtDNA haplotypes, including general haplotypes, assembled haplotypes of electron transport chain contained sequences (H-ETC), mitochondrial respiratory complex contained sequences (H-MRC) and mitochondrial genes (H-gene, including polypeptide-coding genes, rRNA genes and tRNA genes). The strategy for assembled mitochondrial haplotypes was proposed for the first time in mtDNA association analyses on economic traits, although none of the significant relations could be concluded (P > 0.05). In addition, the nuclear major gene BMPR1B was significantly correlated with litter size in the flock (P < 0.05), however, did not interact with mtDNA T7719G mutation (P > 0.05).


Our results highlight mutations of ovine mitochondrial coding genes, suggesting T7719G in tRNA-Lys gene be a potentially useful marker for selection of sheep litter size.


Mitochondria are responsible for ATP production in the electron transport chain (ETC) in cells. The ETC consists of five mitochondrial respiratory complexes (MRCs), of which complex I, complex III, complex IV and complex V are constructed by both mitochondrial and nuclear encoded proteins, and complex II is entirely encoded by nuclear genes. Mitochondrial genome codes 13 polypeptides, 2 rRNAs and 22 tRNAs [1]. Since the first report of mitochondrial DNA (mtDNA) effect on milk production traits in dairy cattle [2], mtDNA effects on various economic traits have been widely studied in livestock, including pigs [3, 4], dairy cattle [5, 6], beef cattle [7, 8], sheep [9, 10] and chickens [11, 12]. Litter size is one of the vital economic traits for animal breeding and production, which has been studied for decades. For sheep, the nuclear gene, BMPR1B was identified as one of causative genes for sheep prolificacy [13], and has been widely used in sheep breeding. For the mtDNA effect on litter size, researchers reported the association with ewe litter size among haplogroups in an Afec-Assaf flock, but did not found the interaction with BMPR1B effect [10]. However, previous studies reported poor mtDNA effects [1417], which made it necessary to uncover the genetic contribution of mtDNA for ewe litter size. In this study, Small-tailed Han sheep, a prolific breed of China, were used to investigate the association of litter size with mtDNA coding genes, in which non-synonymous mutations were considered as possible functional SNPs. Besides the non-synonymous mutation and haplogroup, assembled haplotypes of ETC-contained mtDNA sequences, MRC-contained mtDNA sequences, and mitochondrial genes were analyzed, respectively. The strategy for assembled mitochondrial haplotypes was proposed for the first time in association analyses. In addition, BMPR1B effect and interaction with mtDNA mutations were analyzed.



In total, 117 lambing Small-tailed Han sheep from 53 maternal lineages (families divided by female ancestors) of the same flock were performed blood sample collection, and recorded one or more times of litter size (Additional file 1: Table S1). All sheep were kept indoors year-round, and fed with mixed silage and hay to meet their nutritional requirements. Litter size (number of lambs born) was recorded, and full pedigree information was collected for all animals as reproductive management of the flock where ewes were mated with the selected rams after spontaneous estrus. All the ewes were genotyped for the BMPR1B at >2-month age.

Genotyping of BMPR1B and sequencing of mitochondrial coding genes

Genomic DNA was extracted using the standard phenol/chloroform method [18]. The BMPR1B was genotyped using PCR-RFLP assay with the Ava II restriction enzyme [19]. Mitochondrial complete coding sequences were amplified by 17 primer pairs [20], and PCR products were sequenced in the Sanger method.

Haplotype and haplogroup constitution

All non-synonymous mutations were used to constitute the haplotype and haplogroup. Considering the possible action for mitochondrial function, haplotypes were furthermore assembled by mtDNA sequence of ETC, MRCs and genes respectively. Haplotypes were determined by the online software FaBox [21], and haplogroups were constituted based on network analysis by Network [22].

Association analysis

Association analyses were carried out in the following mixed model by MIXED procedure in SAS software version 9.2 (SAS Institute Inc., Cary, North Carolina, USA).

$$ l s= y s+ parity+ ram+ BMPR 1 B+ mutations+ BMPR 1 B\times mutations+ ID+ E P+ e $$

In the model, the effects of lambing year-season (ys), parity number (parity), service ram (ram), BMPR1B genotype (BMPR1B), mtDNA mutations (mutations, including the effects of mtDNA non-synonymous mutations, haplotypes and haplogroups), the interaction between BMPR1B and mtDNA mutations (BMPR1B × mutations), the polygenic effect (ID), the permanent environmental effect (EP), and the random residual (e) were included. The response variable was the ewe litter size (ls). Each cell of these effects contained observations. The polygenic effect corrected the genetic background by the additive genetic relationship matrix, i.e. the pedigree information. The permanent environmental effect dealed with the repeated measurement data.

Inference about the interaction effect was made. If non-significant, that effect was dropped from the model and inference was made about the main effects. If significant, the cell means for the interaction became of interest.

When a set of statistical inferences were simultaneously considered, multiple comparisons were conducted by the false discovery rate (FDR) in the R project (R version 3.2.5) [23].


Mutations in mitochondrial coding genes

In total, 95 mutations in mitochondrial coding genes were discovered (Additional file 1: Table S2), including 64 synonymous and 31 non-synonymous mutations (19 missense mutations in protein coding genes, 8 mutations in rRNAs, and 4 mutations in tRNAs respectively), which were illustrated in Table 1.

Table 1 Non-synonymous mutations in mitochondrial coding genes and corresponding effects on litter size

BMPR1B genotypes

For BMPR1B gene, 117 ewes were genotyped, including ++ genotype of 18 ewes, B+ genotype of 87 ewes and BB genotype of 12 ewes.

Effects of mitochondrial haplotype and haplogroup on ovine litter size

The 31 mitochondrial non-synonymous mutations assigned the Small-tailed Han sheep flock to 44 haplotypes (Additional file 1: Table S3), which were clustered into 4 haplogroups (Fig. 1). Using the mixed linear model, no significant effect on any haplotype or haplogroup was associated with litter size of Small-tailed Han sheep (P > 0.05), and the interaction between BMPR1B and haplotype or haplogroup also did not significantly affect litter size (P > 0.05), while the BMPR1B was positively associated with litter size (P < 0.05).

Fig. 1

Neighbor-joining tree of different maternal pedigrees for haplogroup assignments by the 31 non-synonymous mutations

Effects of missense mutations and haplotypes in mitochondrial protein coding genes on ovine litter size

Totally 35 H-ETCs were constituted by 19 missense mutations in protein coding sequences, while 15 missense mutations in MRCI assembled 33 H-MRCIs, and both 2 mutations in MRCIV and MRCV constituted 3 H-MRCIVs and 3 H-MRCVs respectively. Intensively, H-genes were assembled by mutations of individual genes (Table 2 and Additional file 1: Table S4, S5 and S6). With mixed linear model analyses, no significant association was detected between litter size and the 19 non-synonymous mutations (P > 0.05) (Table 1), nor was any haplotype (P > 0.05) (Table 2). BMPR1B was strongly associated with litter size (P < 0.05) (Table 3), but the interaction of BMPR1B and mtDNA mutations was not remarkable (P > 0.05).

Table 2 Haplotypes constituted by mitochondrial non-synonymous mutations in multiple levels and significant effects on litter size
Table 3 Effects of BMPR1B on litter size

Effects of mutations and haplotypes in mitochondrial rRNA genes on ovine litter size

The 3 mutations in 12SrRNA sorted the sheep flock into 4 haplotypes (H-12SrRNA), and the 5 mutations in 16SrRNA constituted 15 haplotypes (H-16SrRNA) (Table 2 and Additional file 1: Table S7). Association analyses revealed that neither mutation nor haplotype in rRNA genes affected litter size (P > 0.05) (Tables 1 and 2). The BMPR1B genotype was remarkably associated with litter size (P < 0.05) (Table 3), however, was inconspicuously correlated to rRNA mutations (P > 0.05).

Effects of mutations in mitochondrial tRNA genes on ovine litter size

There were 4 mutations in tRNA genes, and only one variation was observed in each of them (Table 1), therefore no haplotype was constituted. Notably, T7719G in tRNA-Lys affected litter size by 0.29 lambs per litter between the G (1.79) and T (1.50) carriers (P < 0.05). Even though the BMPR1B was in prominent correlation with litter size (P < 0.05) (Table 3), there was no interaction with T7719G (P > 0.05) (Additional file 1: Table S8).


In human diseases, mutations in mitochondrial coding genes led to changes of oxidative phosphorylation enzyme complexes [24]. It is rational to image the possible effects of mtDNA on livestock traits. In this study, the Small-tailed Han sheep were used to explore the correlation of mtDNA and ewe litter size.

In order to investigate the overall mtDNA effects, a novel strategy of assembling mitochondrial haplotypes based on biology functions is proposed. Extensively, there are four levels of mitochondrial haplotypes for mitochondrial biological functions. The general haplotype is assembled by all non-synonymous mutations to reflect the integrated characteristics of mtDNA coding regions. Subsequently, the ETC based haplotype (H-ETC) represents the general feature of all ETC-contained mtDNA sequences. The MRC based haplotype infers the integrated signal of particular MRC-contained mtDNA sequences. H-MRC includes four types, i.e. MRC-I, MRC-III, MRC-IV and MRC-V. The last is the gene level, which indicates the information of a particular gene-contained mtDNA sequence.

In the mixed linear model, unnecessary factors interfere with the estimation of interested factors, as the degree of freedom is wasted. Therefore, based on the thoughts of multiple models [10] and variable selection in a mixed linear model [25], a two-step method of environmental variable selection was put forward in the association analyses, in which non-significant environmental factors were ignored to improve accuracy of genetic estimation. In this study, the environmental factors included lambing month, parity of dam, and service ram, and the genetic factors, which we were interested in, included BMPR1B genotype, mtDNA mutations and their interaction. We started a full model (Model I) with all the environmental and genetic factors to test if environmental factors were noises. Results revealed that the effects of parity of dam and service ram were not significant on litter size. Subsequently, these two factors were excluded in the aims to construct an optimized model (Model II), which was used to test the significance of genetic factors, especially the effects of mtDNA mutations on litter size.

Many mitochondrial tRNA mutations in human were reported to be associated with wide range of pathological conditions [26], for example, deafness was correlated with mitochondrial tRNA-Asp A7551G mutation [27], hypertension was associated with mitochondrial tRNA-Ile A4263G mutation [28], and mitochondrial tRNA mutations might decrease in carcinoma hepatocyte [29]. In sheep, A7755G in tRNA-Lys was reported to affect litter size in an Afec-Assaf flock [10], while our study revealed that T7719G in tRNA-Lys was significantly associated with ewe litter size in Small-tailed Han sheep, with the means of 1.79 and 1.50 for the G and T carriers, manifesting a difference of 0.29 lambs per litter. Furthermore, the mutation was predicted to produce a transversion of U to G at DHU loop on 2D cloverleaf of the tRNA-Lys structure (Fig. 2).

Fig. 2

The predicted secondary structures of the tRNA-Lys for U7719G mutation


As a conclusion, we summarize the research procedures for mtDNA effects. Firstly, sweeping mitochondrial variations on maternal lineages; secondly, constituting the haplotype, haplogroup and assembled haplotypes; lastly, analyzing the association between mtDNA mutations (individual mutations, haplotypes and haplogroups) and interested traits. The present study discovered the mtDNA T7719G was linked with ewe litter size. For traits of low heritability, the marker assisted selection could increase the accuracy of breeding and selection. For the further study, the mtDNA T7719G should be put into post-association validation, and may become a genetic marker in the sheep breeding programs.



Mitochondrial DNA


Electron transport chain


Mitochondrial respiratory complex


The haplotype constituted by all non-synonymous mutations in mitochondrial coding sequences


The haplogroup classified by H


The assembled haplotypes of mtDNA electron transport chain contained sequences


The assembled haplotypes of mtDNA mitochondrial respiratory complex contained sequences, including H-MRCI, H-MRCIII, H-MRCIV, and H-MRCV, corresponding to complex I, III, IV and V, respectively

H-gene :

The assembled haplotypes of mitochondrial genes, including polypeptide-coding genes, rRNA genes and tRNA genes


  1. 1.

    Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson AR, Drouin J, et al. Sequence and organization of the human mitochondrial genome. Nature. 1981;290:457–65.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Bell BR, Mcdaniel BT, Robison OW. Effects of cytoplasmic inheritance on production traits of dairy-cattle. J Dairy Sci. 1985;68:2038–51.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Yen NT, Lin CS, Ju CC, Wang SC, Huang MC. Mitochondrial DNA polymorphism and determination of effects on reproductive trait in pigs. Reprod Domest Anim. 2007;42:387–92.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Fernandez AI, Alves E, Fernandez A, de Pedro E, Lopez-Garcia MA, Ovilo C, et al. Mitochondrial genome polymorphisms associated with longissimus muscle composition in Iberian pigs. J Anim Sci. 2008;86:1283–90.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Qin YH, Chen SY, Lai SJ. Polymorphisms of mitochondrial Atpase 8/6 genes and association with milk production traits in Holstein cows. Anim Biotechnol. 2012;23:204–12.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Paneto JCC, Ferraz JBS, Balieiro JCC, Bittar JFF, Ferreira MBD, Leite MB, et al. Bos indicus or Bos taurus mitochondrial DNA: comparison of productive and reproductive breeding values in a Guzerat dairy herd. Genet Mol Res. 2008;7:592–602.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Zhang B, Chen H, Hua LS, Zhang CL, Kang XT, Wang XZ, et al. Novel SNPs of the mtDNA ND5 gene and their associations with several growth traits in the Nanyang cattle breed. Biochem Genet. 2008;46:362–8.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Auricélio AM, Bittar JFF, Bassi PB, Ronda JB, Bittar ER, Panetto JCC, et al. Influence of endogamy and mitochondrial DNA on immunological parameters in cattle. BMC Vet Res. 2014;10:79–88.

    Article  Google Scholar 

  9. 9.

    Henry BA, Loughnan R, Hickford J, Young IR, St. John JC, Clarke I. Differences in mitochondrial DNA inheritance and function align with body conformation in genetically lean and fat sheep. J Anim Sci. 2015;93:2083–93.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Reicher S, Seroussi E, Weller JI, Rosov A, Gootwine E. Ovine mitochondrial DNA sequence variation and its association with production and reproduction traits within an Afec-Assaf flock. J Anim Sci. 2012;90:2084–91.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Li S, Zadworny D, Aggrey SE, Kuhnlein U. Mitochondrial PEPCK: a highly polymorphic gene with alleles co-selected with Marek’s disease resistance in chickens. Anim Genet. 1998;29:395–7.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Lu WW, Hou LL, Zhang WW, Zhang PF, Chen W, Kang X, et al. Study on heteroplasmic variation and the effect of chicken mitochondrial ND2. Mitochondrial DNA Part A, DNA mapping, sequencing, and analysis. 2016;27:2303-9.

  13. 13.

    Mulsant P, Lecerf F, Fabre S, Schibler L, Monget P, Lanneluc I, et al. Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Merino ewes. Proc Natl Acad Sci U S A. 2001;98:5104–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Hanford KJ, Snowder GD, Vleck LDV. Models with nuclear, cytoplasmic, and environmental effects for production traits of Columbia sheep. J Anim Sci. 2003;81:1926–32.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Van Vleck LD, Snowder GD, Hanford KJ. Models with cytoplasmic effects for birth, weaning, and fleece weights, and litter size, at birth for a population of Targhee sheep. J Anim Sci. 2003;81:61–7.

    Article  PubMed  Google Scholar 

  16. 16.

    Snowder GD, Hanford KJ, Van Vleck LD. Comparison of models including cytoplasmic effects for traits of Rambouillet sheep. Livest Prod Sci. 2004;90:159–66.

    Article  Google Scholar 

  17. 17.

    Van Vleck LD, Hanford KJ, Snowder GD. Lack of evidence for cytoplasmic effects for four traits of Polypay sheep. J Anim Sci. 2005;83:552–6.

    Article  PubMed  Google Scholar 

  18. 18.

    Sambrook J, Russell D. Molecular Cloning: A Laboratory Manual. 3rd ed. Cold Spring Harbor, New York: Cold Spring Harbor Laboratory Press; 2001.

    Google Scholar 

  19. 19.

    Liu FL, Liu YB, Wang F, Wang R, Tian CY, Liu MX. Study on the polymorphism of Bone Morphogenetic Protein Receptor IB in China partial sheep. Acta Agr Boreali-Sinica. 2007;22:151–4.

    Google Scholar 

  20. 20.

    Meadows JRS, Hiendleder S, Kijas JW. Haplogroup relationships between domestic and wild sheep resolved using a mitogenome panel. Heredity. 2011;106:700–6.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Villesen P. FaBox: an online toolbox for fasta sequences. Mol Ecol Notes. 2007;7:965–8.

    CAS  Article  Google Scholar 

  22. 22.

    Network Accessed 5 Sept 2016.

  23. 23.

    R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013. URL

    Google Scholar 

  24. 24.

    Taylor RW, Turnbull DM. Mitochondrial DNA mutations in human disease. Nat Rev Genet. 2005;6:389–402.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Haapalainen E, Laurinen P, Roning J, Kinnunen H. Estimation of exercise energy expenditure using a Wrist-Worn accelerometer: a linear mixed model approach with fixed-effect variable selection. Seventh International Conference on Machine Learning and Applications, Proceedings. 2008. p. 796–801.

    Google Scholar 

  26. 26.

    A Human Mitochondrial Genome Database. Accessed 4 May 2012.

  27. 27.

    Wang M, Peng Y, Zheng J, Zheng B, Jin X, Liu H, et al. A deafness-associated tRNAAsp mutation alters the m1G37 modification, aminoacylation and stability of tRNAAsp and mitochondrial function. Nucleic Acids Res. 2016;44:10974-10985.

  28. 28.

    Chen X, Zhang Y, Xu B, Cai Z, Wang L, Tian J, et al. The mitochondrial calcium uniporter is involved in mitochondrial calcium cycle dysfunction: Underlying mechanism of hypertension associated with mitochondrial tRNA(Ile) A4263G mutation. Int J Biochem Cell Biol. 2016;78:307–14.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Li G, Duan YX, Zhang XB, Wu F. Mitochondrial tRNA mutations may be infrequent in hepatocellular carcinoma patients. Genet Mol Res. 2016;15.

Download references


The authors thank J. C. Jiang (China Agricultural University) for data analysis. The authors also thank Dr. J. M. Zhong and Ms. D. Miller (Auburn University) for language editing.


The research was supported by the Project of Science and Technology of Hebei Province (15226308D) and the National Key Basic Research Program of China (2014CB138500). The former collected the sample and fed the experimental sheep, and the latter designed the study, performed the experiments, analyzed the data, interpreted the data, and drafted the manuscript.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its supplementary information files.

Authors’ contributions

XZ designed the study. WD and XC collected the sample and fed the experimental sheep. the XC, HX, DB and TY performed the experiments. DW, HX and XC analyzed the data. DW, HX, XC and XZ interpreted the data and drafted the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval

All experimental sheep were maintained under the standard farm management practices according to the guidelines of the experimental animal management of China Agricultural University (CAU), and the study was approved by the Experimental Animal Care and Use Committee of CAU. Blood samples were collected followed the guideline of Animal Ethical Committee of China.

Author information



Corresponding author

Correspondence to Xingbo Zhao.

Additional file

Additional file 1: Table S1.

Detail information of used sheep. Table S2. Information of variations in mtDNA coding region of Small-tailed Han sheep. Table S3. The haplotype constituted by the 31 mitochondrial non-synonymous mutations in mitochondrial coding regions (H). Table S4. Assembled haplotypes of mitochondrial electron transport chain contained sequences(H-ETC). Table S5. Assembled haplotypes of mitochondrial respiratory complex contained sequences (H-MRCI, H-MRCIV, and H-MRCV). Table S6. Assembled haplotypes of protein coding genes (H-ND4L, H-ND4, H-ND5, H-ND6, H-ATP6). Table S7. Assembled haplotypes of rRNA coding genes (H-12SrRNA, H-16SrRNA). Table S8. Statistics for the number of observations in each cell of the interaction between BMPR1B and T7719G in tRNA-Lys and corresponding effects on litter size. (XLSX 33 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chen, X., Wang, D., Xiang, H. et al. Mitochondrial DNA T7719G in tRNA-Lys gene affects litter size in Small-tailed Han sheep. J Animal Sci Biotechnol 8, 31 (2017).

Download citation


  • Association
  • Haplotype
  • Mitochondria
  • Non-synonymous mutation
  • Reproduction
  • Sheep