Molecular analysis of cytoplasmic genetic variation in Holstein cows

The role of mtDNA in oocyte quality and embryo development

The mitochondrial genome resides in the mitochondria present in nearly all cell types. The porcine (Sus scrofa) mitochondrial genome is circa 16.7 kb in size and exists in the multimeric format in cells. Individual cell types have different numbers of mitochondrial DNA (mtDNA) copy number based on their requirements for ATP produced by oxidative phosphorylation. The oocyte has the largest number of mtDNA of any cell type. During oogenesis, the oocyte sets mtDNA copy number in order that sufficient copies are available to support subsequent developmental events. It also initiates a program of epigenetic patterning that regulates, for example, DNA methylation levels of the nuclear genome. Once fertilized, the nuclear and mitochondrial genomes establish synchrony to ensure that the embryo and fetus can complete each developmental milestone. However, altering the oocyte's mtDNA copy number by mitochondrial supplementation can affect the programming and gene expression profiles of the developing embryo and, in oocytes deficient of mtDNA, it appears to have a positive impact on the embryo development rates and gene expression profiles. Furthermore, mtDNA haplotypes, which define common maternal origins, appear to affect developmental outcomes and certain reproductive traits. Nevertheless, the manipulation of the mitochondrial content of an oocyte might have a developmental advantage.

Genetic and phenotypic associations of mitochondrial DNA copy number, SNP, and haplogroups with growth and carcass traits in beef cattle

Mitochondrial DNA copy number (mtDNA CN) is heritable and easily obtained from low-pass sequencing (LPS). This study investigated the genetic correlation of mtDNA CN with growth and carcass traits in a multi-breed and crossbred beef cattle population. Blood, leucocyte, and semen samples were obtained from 2,371 animals and subjected to LPS that resulted in nuclear DNA (nuDNA) and mtDNA sequence reads. Mitochondrial DNA CN was estimated as the ratio of mtDNA to nuDNA coverages. Variant calling was performed from mtDNA, and 11 single nucleotide polymorphisms (SNP) were identified in the population. Samples were classified in taurine haplogroups. Haplogroup and mtDNA type were further classified based on the 11 segregating SNP. Growth and carcass traits were available for between 7,249 and 60,989 individuals. Associations of mtDNA CN, mtDNA haplogroups, mtDNA types, and mtDNA SNP with growth and carcass traits were estimated with univariate animal models, and genetic correlations were estimated with a bivariate animal model based on pedigree. Mitochondrial DNA CN tended (P-value ≤0.08) to be associated with birth weight and weaning weight. There was no association (P-value >0.10) between mtDNA SNP, haplogroups, or types with growth and carcass traits. Genetic correlation estimates of mtDNA CN were -0.30 ± 0.16 with birth weight, -0.31 ± 0.16 with weaning weight, -0.15 ± 0.14 with post-weaning gain, -0.11 ± 0.19 with average daily dry-matter intake, -0.04 ± 0.22 with average daily gain, -0.29 ± 0.13 with mature cow weight, -0.11 ± 0.13 with slaughter weight, -0.14 ± 0.13 with carcass weight, -0.07 ± 0.14 with carcass backfat, 0.14 ± 0.14 with carcass marbling, and -0.06 ± 0.14 with ribeye area. In conclusion, mtDNA CN was negatively correlated with most traits investigated, and the genetic correlation was stronger with growth traits than with carcass traits.

The ND1 and CYTB genes polymorphisms associated with in vitro early embryo development of Sanjabi sheep

This study aimed to investigate the association between polymorphisms of ND1 and CYTB genes and in vitro early embryo development of Sanjabi sheep. Blood and ovarian samples were collected from a local slaughterhouse. The cumulus-oocyte complexes with a diameter greater than 3 mm were aspirated from follicles, and in vitro maturation (IVM) and in vitro culture (IVC) rates of them were recorded. A respective 1200 bp and 980 bp fragments of ND1 and CYTB genes were genotyped using a modified single strand conformation polymorphism (SSCP) method. The results of this study revealed that four different patterns, named as A, B, C, and D were observed for both ND1 and CYTB genes. The ND1 gene polymorphisms had significant effects on the IVM and IVC rate (p < 0.05). The pattern C of the ND1 gene significantly increased the IVM rate compared to the patterns A, B and D. For the IVC, the highest and lowest means were related to the C and B patterns, respectively. The CYTB gene polymorphisms also had significant effects on IVC (p < 0.01), but the IVM did not affected (p = 0.07). Here, the pattern D had the highest and the pattern C had the lowest means for both IVM and IVC rates.

Cattle phenotypes can disguise their maternal ancestry

Background

Cattle are bred for, amongst other factors, specific traits, including parasite resistance and adaptation to climate. However, the influence and inheritance of mitochondrial DNA (mtDNA) are not usually considered in breeding programmes. In this study, we analysed the mtDNA profiles of cattle from Victoria (VIC), southern Australia, which is a temperate climate, and the Northern Territory (NT), the northern part of Australia, which has a tropical climate, to determine if the mtDNA profiles of these cattle are indicative of breed and phenotype, and whether these profiles are appropriate for their environments.

Results

A phylogenetic tree of the full mtDNA sequences of different breeds of cattle, which were obtained from the NCBI database, showed that the mtDNA profiles of cattle do not always reflect their phenotype as some cattle with Bos taurus phenotypes had Bos indicus mtDNA, whilst some cattle with Bos indicus phenotypes had Bos taurus mtDNA. Using D-loop sequencing, we were able to contrast the phenotypes and mtDNA profiles from different species of cattle from the 2 distinct cattle breeding regions of Australia. We found that 67 of the 121 cattle with Bos indicus phenotypes from NT (55.4%) had Bos taurus mtDNA. In VIC, 92 of the 225 cattle with Bos taurus phenotypes (40.9%) possessed Bos indicus mtDNA. When focusing on oocytes from cattle with the Bos taurus phenotype in VIC, their respective oocytes with Bos indicus mtDNA had significantly lower levels of mtDNA copy number compared with oocytes possessing Bos taurus mtDNA (P < 0.01). However, embryos derived from oocytes with Bos indicus mtDNA had the same ability to develop to the blastocyst stage and the levels of mtDNA copy number in their blastocysts were similar to blastocysts derived from oocytes harbouring Bos taurus mtDNA. Nevertheless, oocytes originating from the Bos indicus phenotype exhibited lower developmental potential due to low mtDNA copy number when compared with oocytes from cattle with a Bos taurus phenotype.

Conclusions

The phenotype of cattle is not always related to their mtDNA profiles. MtDNA profiles should be considered for breeding programmes as they also influence phenotypic traits and reproductive capacity in terms of oocyte quality.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0523-5) contains supplementary material, which is available to authorized users.

The relationship between mitochondrial DNA haplotype and the reproductive capacity of domestic pigs (Sus scrofa domesticus)

Background

The maternally inherited mitochondrial genome encodes key proteins of the electron transfer chain, which produces the vast majority of cellular ATP. Mitochondrial DNA (mtDNA) present in the mature oocyte acts as a template for all mtDNA that is replicated during development to meet the specific energy requirements of each tissue. Individuals that share a maternal lineage cluster into groupings known as mtDNA haplotypes. MtDNA haplotypes confer advantages and disadvantages to an organism and this affects its phenotype. In livestock, certain mtDNA haplotypes are associated with improved milk and meat quality, whilst, other species, mtDNA haplotypes have shown increased longevity, growth and susceptibility to diseases. In this work, we have set out to determine whether mtDNA haplotypes influence reproductive capacity. This has been undertaken using a pig model.

Results

To determine the genetic diversity of domestic pigs in Australia, we have sequenced the D-loop region of 368 pigs, and identified five mtDNA haplotypes (A to E). To assess reproductive capacity, we compared oocyte maturation, fertilization and development to blastocyst, and found that there were significant differences for maturation and fertilization amongst the haplotypes. We then determined that haplotypes C, D and E produced significantly larger litters. When we assessed the conversion of developmentally competent oocytes and their subsequent developmental stages to offspring, we found that haplotypes A and B had the lowest reproductive efficiencies. Amongst the mtDNA haplotypes, the number of mtDNA variants harbored at >25 % correlated with oocyte quality. MtDNA copy number for developmentally competent oocytes positively correlated with the level of the 16383delC variant. This variant is located in the conserved sequence box II, which is a regulatory region for mtDNA transcription and replication.

Conclusions

We have identified five mtDNA haplotypes in Australian domestic pigs indicating that genetic diversity is restricted. We have also shown that there are differences in reproductive capacity amongst the mtDNA haplotypes. We conclude that mtDNA haplotypes affect pig reproductive capacity and can be used as a marker to complement current selection methods to identify productive pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0375-4) contains supplementary material, which is available to authorized users.

The control of mtDNA replication during differentiation and development

BACKGROUND

Mitochondrial DNA (mtDNA) is important for energy production as it encodes some of the key genes of electron transfer chain, where the majority of cellular energy is generated through oxidative phosphorylation (OXPHOS). MtDNA replication is mediated by nuclear DNA-encoded proteins or enzymes, which translocate to the mitochondria, and is strictly regulated throughout development. It starts with approximately 200 copies in each primordial germ cell and these copies undergo expansion and restriction events at various stages of development.

SCOPE OF REVIEW

I describe the patterns of mtDNA replication at key stages of development. I explain that it is essential to regulate mtDNA copy number and to establish the mtDNA set point in order that the mature, specialised cell acquires the appropriate numbers of mtDNA copy to generate sufficient adenosine triphosphate (ATP) through OXPHOS to undertake its specialised function. I discuss how these processes are dependent on the controlled expression of the nuclear-encoded mtDNA-specific replication factors and that this can be modulated by mtDNA haplotypes. I discuss how these events are altered by certain assisted reproductive technologies, some of which have been proposed to prevent the transmission of mutant mtDNA and others to overcome infertility. Furthermore, some of these technologies are predisposed to transmitting two or more populations of mtDNA, which can be extremely harmful.

MAJOR CONCLUSIONS

The failure to regulate mtDNA replication and mtDNA transmission during development is disadvantageous.

GENERAL SIGNIFICANCE

Manipulation of oocytes and embryos can lead to significant implications for the maternal-only transmission of mtDNA. This article is part of a Special Issue entitled Frontiers of mitochondrial research.

The phylogenetic status of typical Chinese native pigs: analyzed by Asian and European pig mitochondrial genome sequences

China is one of the most diverse countries, which have developed 88 indigenous pig breeds. Several studies showed that pigs were independently domesticated in multiple regions of the world. The purpose of this study was to investigate the origin and evolution of Chinese pigs using complete mitochondrial genomic sequences (mtDNA) from Asian and European domestic pigs and wild boars. Thirty primer pairs were designed to determine the mtDNA sequences of Xiang pig, Large White, Lantang, Jinhua and Pietrain. The phylogenetic status of Chinese native pigs was investigated by comparing the mtDNA sequences of complete coding regions and D-loop regions respectively amongst Asian breeds, European breeds and wild boars. The analyzed results by two cluster methods contributed to the same conclusion that all pigs were classified into two major groups, European clade and Asian clade. It revealed that Chinese pigs were only recently diverged from each other and distinctly different from European pigs. Berkshire was clustered with Asian pigs and Chinese pigs were involved in the development of Berkshire breeding. The Malaysian wild boar had distant genetic relationship with European and Asian pigs. Jinhua and Lanyu pigs had more nucleotide diversity with Chinese pigs although they all belonged to the Asian major clade. Chinese domestic pigs were clustered with wild boars in Yangtze River region and South China.

Association between mitochondrial DNA haplotype compatibility and increased efficiency of bovine intersubspecies cloning

Reconstructed embryos derived from intersubspecies somatic cell nuclear transfer (SCNT) have poorer developmental potential than those from intrasubspecies SCNT. Based on our previous study that Holstein dairy bovine (HD) mitochondrial DNA (mtDNA) haplotype compatibility between donor karyoplast and recipient cytoplast is crucial for SCNT embryo development, we performed intersubspecies SCNT using HD as donor karyoplast and Luxi yellow heifer (LY) as recipient cytoplast according to mtDNA haplotypes determined by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis. The results demonstrated that intersubspecies mtDNA homotype SCNT embryos had higher pre- and post-implantation developmental competence than intrasubspecies mtDNA heterotype embryos as well as improved blastocyst reprogramming status, including normal H3K9 dimethylation pattern and promoter hypomethylation of pluripotent genes such as Oct4 and Sox2, suggesting that intersubspecies SCNT using LY oocytes maintains HD cloning efficiency and may reprogram HD nuclei to develop into a normal cloned animal ultimately. Our results indicated that karyoplast-cytoplast interactions and mtDNA haplotype compatibility may affect bovine intersubspecies SCNT efficiency. This study on bovine intersubspecies SCNT is valuable for understanding the mechanisms of mtDNA haplotype compatibility between karyoplast and cytoplast impacting the bovine SCNT efficiency, and provides an alternative and economic resource for HD cloning.

Donor-host mitochondrial compatibility improves efficiency of bovine somatic cell nuclear transfer

Background

The interaction between the karyoplast and cytoplast plays an important role in the efficiency of somatic cell nuclear transfer (SCNT), but the underlying mechanism remains unclear. It is generally accepted that in nuclear transfer embryos, the reprogramming of gene expression is induced by epigenetic mechanisms and does not involve modifications of DNA sequences. In cattle, oocytes with various mitochondrial DNA (mtDNA) haplotypes usually have different ATP content and can further affect the efficiency of in vitro production of embryos. As mtDNA comes from the recipient oocyte during SCNT and is regulated by genes in the donor nucleus, it is a perfect model to investigate the interaction between donor nuclei and host oocytes in SCNT.

Results

We investigated whether the in vitro development of reconstructed bovine embryos produced by SCNT would be influenced by mtDNA haplotype compatibility between the oocytes and donor cells. Embryos from homotype A-A or B-B showed significantly higher developmental ability at blastocyst stages than the heterotype A-B or B-A combinations. Post-implantation development ability, pregnancy rate up to day 90 of gestation, as well as percent of term births were higher in the homotype SCNT groups than in the heterotype groups. In addition, homotype and heterotype SCNT embryos showed different methylation patterns of histone 3-lysine 9 (H3K9) genome-wide and at pluripotency-related genes (Oct-4, Sox-2, Nanog).

Conclusion

Both histone and DNA methylation show that homotype SCNT blastocysts have a more successful epigenetic asymmetry pattern than heterotype SCNT blastocysts, which indicates more complete nuclear reprogramming. This may result from variability in their epigenetic patterns and responses to nuclear reprogramming. This suggests that the compatibility of mtDNA haplotypes between donor cells and host oocytes can significantly affect the developmental competence of reconstructed embryos in SCNT, and may include an epigenetic mechanism.

Ovine mitochondrial DNA: restriction enzyme analysis, mapping and sequencing data

Restriction endonuclease fragment patterns of mitochondrial DNA (mtDNA) in sheep were analysed with 11 enzymes. Four breeds (Merinolandschaf, Rhoenschaf, Schwarzkoepfiges Fleischschaf and Skudde) of domestic sheep and European Mouflon were examined. A restriction map with 28 cleavage sites of seven enzymes was established. KpnI and PstI do not cut ovine mtDNA. Two EcoRI fragments of Merinolandschaf, Rhoenschaf and Mouflon each were cloned and partially sequenced. Intraspecific nucleotide sequence differences within 1.101 kb ranged from 0.09 to 0.27%. Hybridization analysis with a fragment of porcine mtDNA along with sequencing data from cloned fragments was used for orientation of the restriction map along the bovine sequence. Ovine mtDNA sequences encompassing parts of the Cyt.b-, ND5-, CoIII- and ATPase6 genes were compared with the corresponding sequences of the bovine mtDNA.

Effect of oocyte mitochondrial DNA haplotype on bovine somatic cell nuclear transfer efficiency

The development capability of reconstructed bovine embryos via ovum pick-up (OPU)-somatic cell nuclear transfer (SCNT) technique has been influenced by the maternal lineage of oocyte cytoplasm, but the underlying mechanism remains unclear. Since mitochondria are the richest maternal-inherited organelle, in this study, we intended to clarify the effect of mtDNA haplotypes on cloning efficiency. By PCR-RFLP method, we identified mtDNA haplotypes A and B, differing in six restriction sites. Reconstructed embryos with haplotype A cytoplast achieved better fusion and blastocyst formation rate (64.6% and 39.4%), as compared with haplotype B (53.6% and 26.3%; P < 0.05). To further evaluate the role of mitochondria, the quantity of mtDNA, ATP content, and mRNA level of mtDNA-encoded COXI, COXIII in both oocytes were measured. Our data indicated that mtDNA copy number in haplotype A oocyte was significantly higher than that in haplotype B oocyte, both at the GV (10(5.03 +/- 0.69) vs. 10(4.81 +/- 0.86) copies/oocyte) and MII stages (10(5.31 +/- 0.71) vs. 10(5.13 +/- 0.63) copies/oocyte; logarithmically transformed values; P < 0.05). ATP content in type A oocyte was also greater at the GV (1.67 +/- 0.09 vs. 1.27 +/- 0.1 pmol) and MII stages (5.18 +/- 0.07 vs. 2.68 +/- 0.03 pmol; P < 0.05). Similarly, the mRNA expression level of mtDNA-encoded COXI and COXIII in haplotype A oocyte was significantly higher comparing to haplotype B oocyte (3.3 +/- 2.0 x 10(3) vs. 0.68 +/- 0.45 x 10(3); 24.9 +/- 10.5 x 10(3) vs. 9.4 +/- 3.3 x 10(3), respectively; P < 0.05). The data suggest that mitochondrial structure, quantity, and function may significantly affect the developmental competence of reconstructed embryos.

Mitochondrial DNA Polymorphism and Determination of Effects on Reproductive Trait in Pigs

The purposes of this study were to investigate single strand conformation polymorphisms (SSCP) in the D-loop region of pig mitochondrial DNA (mtDNA) and to determinate their association with the reproductive traits of meishan pigs. A total of four types of band patterns, designed SSCP band pattern A, B, C and D, were identified. A type of SSCP band pattern was present in all European-American breeds, but not in East Asian breeds. This result showed the diversified sequence in the D-loop region between European-American and East Asian populations. Two types of band patterns, B and C, were found in Meishan pigs. The average body weight at day 21 of piglets from B type dams was significantly heavier than the body weight of C type (p < 0.05). We also tested whether the SSCP patterns would be suitable for paternity testing in a family group and found that bands of all the offspring were derived from their maternal parent. Therefore, we conclude that SSCP may be a marker for identification of maternal ancestors.

Transmission of mitochondrial DNA following assisted reproduction and nuclear transfer

Mitochondria are the organelles responsible for producing the majority of a cell's ATP and also play an essential role in gamete maturation and embryo development. ATP production within the mitochondria is dependent on proteins encoded by both the nuclear and the mitochondrial genomes, therefore co-ordination between the two genomes is vital for cell survival. To assist with this co-ordination, cells normally contain only one type of mitochondrial DNA (mtDNA) termed homoplasmy. Occasionally, however, two or more types of mtDNA are present termed heteroplasmy. This can result from a combination of mutant and wild-type mtDNA molecules or from a combination of wild-type mtDNA variants. As heteroplasmy can result in mitochondrial disease, various mechanisms exist in the natural fertilization process to ensure the maternal-only transmission of mtDNA and the maintenance of homoplasmy in future generations. However, there is now an increasing use of invasive oocyte reconstruction protocols, which tend to bypass mechanisms for the maintenance of homoplasmy, potentially resulting in the transmission of either form of mtDNA heteroplasmy. Indeed, heteroplasmy caused by combinations of wild-type variants has been reported following cytoplasmic transfer (CT) in the human and following nuclear transfer (NT) in various animal species. Other techniques, such as germinal vesicle transfer and pronuclei transfer, have been proposed as methods of preventing transmission of mitochondrial diseases to future generations. However, resulting embryos and offspring may contain mtDNA heteroplasmy, which itself could result in mitochondrial disease. It is therefore essential that uniparental transmission of mtDNA is ensured before these techniques are used therapeutically.

Identification of mitochondrial DNA substitutions related to meat quality in Japanese Black cattle

Complete sequences of mitochondrial (mt) genomes of eight Japanese Black cattle were determined to investigate the relationships between mt deoxyribonucleic acid (DNA) displacement loop (D-loop) types and other mtDNA regions and to identify the variation in the coding region that may influence the economic traits. The survey of mitochondrial sequences in the encoding region revealed 14 substitutions including six antonymous substitutions and one in 16S ribosomal ribonucleic acid (rRNA). Three methods of polymorphic DNA analyses (polymerase chain reaction [PCR]-restriction fragment length polymorphism [RFLP], mismatch PCR-RFLP, PCR-single-strand conformation polymorphism [SSCP]) were performed on these seven candidate substitutions (base pair [bp] 2,232, 12,158, 12,908, 13,310, 14,122, 14,140, and 14,565) for 202 Japanese Black cattle. The substitution of bp 13,310 was observed in all samples, but not in the reference sequence, indicating that this is a minor substitution or a sequencing mistake in the reference sequence. The substitutions at bp 14,122, 14,140, and 14,565 were observed in only a few samples, suggesting that these were also minor substitutions. The substitutions at bp 2,232 (16S rRNA), 12,158, and 12,908 (reduced nicotinamide adenine dinucleotide-ubiquinone oxidoreductase chain-5) were closely related to mitochondrial D-loop types that have previously been related to differences in the carcass traits of Japanese Black cattle. Evaluation of the effects on six carcass traits with mixed model procedures suggests that the bp 2,232 substitution affects longissimus muscle area and beef marbling score. The substitution at bp 2,232 is a strong candidate for the mitochondrial effect on meat quality.

Mitochondrial DNA polymorphisms and fertility in beef cattle

Four hundred and twenty-two beef cattle of two different breeds (purebred Hereford and composite multibreed) were characterized by polymerase chain reaction-restriction fragment length polymorphism, using the restriction enzymes ApaI, AvaII, HindIII, PstI, SpeI, SspI and TaqI in two regions (the D-loop and the ND-5 gene) of mitochondrial DNA. The association between molecular haplotypes and records on calving rate, defined as the mean number of live calves born per year over 4 years, were examined by analysis of variance. A significant association was found between calving rate and mitochondrial polymorphisms in both breeds. This may have implications for genetically improving cow fertility.

Genetic variation in functionally important domains of the bovine mtDNA control region

DNA sequences of the mitochondrial control region (CR) of 32 unrelated Austrian cattle were analysed in order to determine the extent of variability in functionally important domains. Using sequencing of PCR products, allele-specific PCR (AS-PCR) and primer introduced restriction analysis (PIRA), 43 differences were observed. They included 33 transitions, five transversions, one deletion and four differences in the number of consecutive cytosines. Twenty-three of these polymorphisms have not been reported before. In addition, we analysed all available European cattle sequences for this region. The transcriptional start sites, the conserved sequence block CSB 1 and both binding sites for the mitochondrial transcription factor mtTFA were highly conserved. We found a transition in each of the inter-specifically conserved Mt4 and Mt5 elements, three nucleotide substitutions in the termination-associated sequence TAS-A and six polymorphisms in the conserved sequence block CSB 2+3, a region which has been implicated in mitochondrial RNA processing.

Variation in the control region sequence of the sheep mitochondrial genome

The DNA sequences of the control region of the mitochondrial genome of fifty unrelated sheep were determined in order to ascertain the extent and distribution of its variability. A consensus sequence was derived, and 1081 differences from it were observed amongst the fifty animals. Some constant groups of differences were observed that were held in common by a number of animals, which thus fell into two main groups, although neither group was typical of any of the breeds sampled. The consensus sequence also allowed comparison between the control region sequences of sheep and other mammals. The sequence contains four tandem repeats of a 75 base-pair motif that accounts for the difference in its size from the cattle control region, to which it is otherwise very similar. Comparison with the cattle sequence allowed the determination of the homologues of various functionally important sites. The homologues of the transcription promoters, the origin of replication and the central conserved sequence block were all identified by this method.

Mapping the sheep genome: practice, progress and promise

Over 100 loci have now been mapped on the sheep genome, double that of 3 years ago. Achieved through converging developments in DNA technologies and genetics, this rate of gene mapping is increasing. Its goal is to identify genes for important traits. This will elucidate the biological basis of these traits, paving the way for their manipulation, and accelerate the genetic improvement of sheep and quality of their products by marker-assisted selection. The aim of this review is to provide an introduction to the strategies and methods of gene mapping, highlighting the value of comparative genome analysis. These themes are then elaborated in the description of the ongoing project to map the sheep genome. Opportunities are discussed for the application of the fruits of this research to the understanding and treatment of diseases, and to the development of tools for the genetic improvement of sheep by marker-assisted selection.