Abstract

Traditional analyses of a QTL on Bota 19 implicate a surfeit of candidates, but each is of marginal significance in explaining the deposition of healthy, low melting temperature fat within marbled muscle of Wagyu cattle. As an alternative approach, we have used genomic, multigenerational segregation to identify 14 conserved, ancestral 20 Mb haplotypes. These determine the degree and rate of marbling in Wagyu and other breeds of cattle. The melting temperature of intramuscular fat is highly heritable and traceable by haplotyping. Fortunately, for the production of healthy beef, some of these haplotypes are sufficiently penetrant to be expressed in heterozygous crossbreds, thereby allowing selection of sires which will improve the healthiness of beef produced under even harsh climatic conditions. The region of Bota 19 is syntenic to a region of Hosa 17 known to be important in muscle metabolism and in determining susceptibility to a form of human muscular dystrophy.

1. Introduction

Conserved ancestral haplotypes have been highly informative in revealing the functional significance of polymorphic sequences especially when extensive enough to encode regulatory regions as well as structural elements [13].

In Homo sapiens (Hosa), there are now many examples of conserved haplotypes regulating the degree of expression of a trait or disease such as an enzyme deficiency, immune competence, drug reactions, and disease severity. The best known examples are the extended haplotypes of the major histocompatibility complex [1, 2, 47], but the principle applies throughout the genome, [8] with the important caveat that the haplotypes are identified and tagged correctly [9, 10]. Many of these haplotypes were discovered through the finding of genetic associations with quantitative traits, such as diseases. It became obvious that there are tens or even hundreds of alternative haplotypes encoding functionally important regions of the genome and that these interact with each other positively or negatively. A potent example is IgA deficiency [2, 6, 1113] which is associated with the 8.1 haplotype inter alia. The precise mechanisms underlying this regulatory defect are not clear from analysis of SNPs and coding sequences [14].

For the breeding of livestock, it is traditional [15, 16] to follow the precepts of the infinitesimal model which assumes there are independent genes separated by free recombination [17]. However, it is quite clear that this assumption is not correct [3, 17]. The genomic architecture of cattle [1821], sheep [22, 23], goats [24], dogs [25, 26], mice [27], and pigs [28, 29] is remarkably similar to that present in Hosa, where there are kilobase and even megabase blocks of polymorphism at which there are retroviral and other indels, segmental duplications, and differences in copy number in addition to nucleotide substitutions [1, 2, 30, 31].

In cattle, there are many traits of proven importance to the production of healthy food. There have been numerous attempts to find useful genetic markers if only for selection of superior breeding pairs for commercial gain. Results have been disappointing for several reasons generally reflecting an oversimplified understanding of the genomic structure. Thus, conserved ancestral haplotypes can be confounding if not identified. For example, desirable traits such as increased milk production can be associated with disadvantages such as infertility [3]. In other instances, the genetic marker may appear useful in one breed but not in another, as illustrated by polling [3].

Here, we address these disappointments by identifying megabase haplotypes of critical interest to healthy food and public health [32, 33] and potential relevance to the metabolic changes underlying limb girdle muscular dystrophy.

True or fine marbling describes the deposition of lipid between and within muscle fibres after supplementary feeding for a year or more occurs. The genesis is very complex, including differentiation of satellite or progenitor cells, and is obviously tightly regulated, [34], but the ultimate effect is to improve taste and juiciness thereby increasing the value of beef several folds. More importantly, marbling increases the content of the healthy, monounsaturated fatty acid (MUFA) and is therefore “cholesterol lowering” and statin-sparing [35, 36]. These features are characteristic of Japanese cattle, known as Wagyu, black and red, but may be found to lesser degrees in other breeds and crossbreeds. It is important to distinguish true or fine marbling from the deposition of subcutaneous and interfascicular fat as a characteristic of European breeds.

For more than 20 years, initially in Japan, there have been attempts to find predictors of (i)how much a given animal will marble;(ii)which matings are preferable;(iii)whether the very costly feeding can be reduced.

At face value, the results are contradictory but also tantalising since the inconsistencies may be due to the complexity of the marbling pathway and the sheer number of genes involved.

A simple polymorphism of the SCD (delta-9 desaturase on Bota 26) appears to affect the amount of marbling presumably by encouraging desaturation and therefore the production of MUFA and particularly oleic acid. Polymorphisms of growth hormone (GH) and fatty acid synthase (FASN) amongst many others may also have an effect. The evidence for some of these individual genes is reviewed and presented in Table 1 [3754]. From these 14 studies, the only fair conclusion is that, although the genomic region controls fat deposition, multiple genes and therefore haplotypes must be involved.

Our interest increased when it became apparent that a regulator of SCD, SREBF1, maps to the Bota 19 chromosome together with GH and FASN. Taken with the pleiotropic regulatory role of SREBF1, in humans at least [55], and with the evidence for the influence of Bota 19 on fatty acid composition in Japanese black cattle [48, 49], we asked whether extensive, multimegabase haplotypes could be identified. Previous investigations [3754] have suggested that some intervening and adjacent genes could be relevant. Such apparent clustering of plausible candidates is often illusory, but in this case, we were encouraged by the fact that Japanese breeders are sufficiently convinced to be prepared to pay handsomely for SCD and GH typing even though it seems clear that these genes are unlikely to be involved directly.

Accordingly, we have been defining conserved ancestral haplotypes within the region. The many conserved megabase haplotypes have very different frequencies when breeds are compared [56]. The need for a simple summary of breed differences led to the development of the W plots which separate haplotypes common to all cattle from those which are found in, and potentially defining for, one breed but not in the other [57, 58]. Of course, since Wagyu cattle marble, it follows that all Wagyu-specific haplotypes must be associated with marbling irrespective of their function; all too often, such passive associations have been misinterpreted as causal. Common haplotypes may be implicated unless nonmarbling breeds are compared to determine which are specific for marbling breeds. The challenge is to assess the Wagyu haplotypes in different backgrounds such as can be obtained by crossing Wagyu with other breeds which are relatively resistant to marbling.

As a strategy, we asked whether Wagyu-specific haplotypes, once defined, are penetrant when heterozygous in F1 crossbreds. If so, it would be fair to conclude that these haplotypes have direct effects on the marbling process. We also asked whether non-Wagyu haplotypes affect the rate of marbling alone or in combination. Therefore, we started the investigation by determining whether fine marbling of progeny is a function of parentage.

In order to undertake such studies, it is critical to have a reliable measurement of marbling. In spite of many attempts, including the use of photographic and ultrasound images, there is still no agreed international benchmark, largely because the deposition of the preferred fat is very fine and barely visible. The most objective approach is the cumbersome slip test for the measurement of melting temperature (Tm) of the lipid fraction and effectively the predominance of oleic acid [34, 59]. In Wagyu, Tms below body temperature (37°C) correlate with fine marbling thereby explaining the cholesterol lowering and statin-sparing effects [34].

We have introduced a more convenient and highly reproducible variation using the same thermocycler as used for the polymerase chain reaction [60].

With this method, elsewhere we will show that low Tm is highly heritable as expected given its extreme expression in only the Wagyu breed. Furthermore, we show that the effect is determined by multiple, conserved, megabase haplotypes. These can be used to predict marbling in crossbreds. There is evidence for cis-interactions within and trans-interactions between haplotypes. We also describe a strategy for improving the quality and healthiness of beef by judicious crossing.

The step-wise approach used here can be applied to other challenges in livestock genetics especially where regulatory ancestral haplotypes are suspected. Systematically, we address (i)synteny as a guide to functions and interactions already defined in Hosa;(ii)breed effects as an initial screen for heritability;(iii)heritability as the only variable;(iv)genomic markers of inheritance;(v)haplotypes in crossbreds;(vi)mechanisms of haplotype effects.

The conclusions of all these complementary studies lead to a model of the genetic control of fatty acid composition.

2. Results

2.1. Synteny as a Guide

In Figure 1, we show the regions of Hosa chr17 which correspond to SREBF1 to FASN on Bota chr 19. Note also the regulatory effects of SREBF1 on remote sites many of which map closely with locations implicated in limb girdle muscular dystrophy.

In the dot plot of GH to FASN, it is apparent that there are two separate components in Hosa; the blue and green sections are distinct but note how they remain contiguous even though the green has been inverted. Similar phenomena are present elsewhere in the synteny plot. Together, these suggest that there are multiple interacting haplotypes from SREBF1 to FASN.

2.2. Breed Effects and Heritability

Most Wagyu have Tm of intramuscular fat below human body temperature of 37°C and distinctly lower than fat from Simmentals which is mostly subcutaneous and generally above 40°C [60]. The length of feeding was much greater for Wagyu due to the fact that the Simmentals grow rapidly and fatten quickly through deposition of largely subcutaneous fat. Wagyu are much slower to both grow and fatten as they deposit intramuscular rather than subcutaneous fat.

We compared Wagyu cattle with various European breeds and crossbreeds which were fed the same ration under the same environmental conditions, but with days on feed varied to achieve the desired finish. The results showed that the Tm of the subcutaneous fat was lower in animals with higher proportions of Wagyu ancestry (data not shown here). Differences were noticeable in animals with only 25% Wagyu ancestry indicating that the breed effect can be highly penetrant.

Wagyu steers were fed together for 300 days and assessed identically. The only significant variable was parentage which was confirmed by DNA haplotyping. All sires are considered to be elite in terms of their breeding, reputation, and estimated breeding value (EBV) when available. The progeny groups differ significantly, with mean Tm for the progeny of two sires differing by nearly 2 degrees. Importantly, this effect was not seen with the visual marbling score. There is an unequivocal genetic effect detectable within one generation. Thus, in spite of the undoubted complexity of the inheritance, sires can be ranked according to heritability of low Tm.

Thus, breed differences, individual sires, and presumably identifiable genetic factors, determine the type, rate, and location of lipid deposition. However, it is also clear that there is a major effect of environment in that Tm falls with increasing DOF irrespective of breed, and this fact compounded by imprecision may have contributed to past confusion.

In the analysis of genetic differences presented here, we define alleles and haplotypes common in Wagyu but rare in other breeds that may explain some of the differences in the type, rate, and location of lipid deposition.

2.3. Genomic Markers of Inheritance

The genomic structure of the relevant region of Bota 19 contains the previously implicated SREBF1 and GH, as shown in Figures 1 and 2. Note the extensive duplication as found in other genomic regions of polymorphism including differences in copy number. The degree of polymorphism (at the insensitive gel level) ranges from at least 7 alleles at the most polymorphic locus (MPRIP) to only 2 at the least (TCAP). In the SREBF1 to GH (S-G) segment, there are 3 × 5 × 7 × 2 × 3 = 630 possible haplotypes on each paternal or maternal chromosome. Therefore, there should be some thousands of genotypes if the alleles at each locus are segregating randomly due to free recombination as specified in the infinitesimal model of population genetics.

Here, we report the actual frequencies in Wagyu and in some of the breeds with which they have been crossed. Only 14 haplotypes occur commonly, say greater than 5%, in the breeds studied here. Importantly, as shown in Figures 3 and 4, these are the conserved ancestral haplotypes. Four of the fourteen are present in the three-generation family which also shows the segregation of the (a) haplotype in three genotypes. None of the haplotypes could be identified through linkage disequilibrium; they must be observed empirically. They are not recent recombinants or mutants. They are not those expected from the allele frequencies and random assortment. Random segregation and free recombination can be excluded at least in this 20-megabase region.

As shown in Figure 2, haplotypes are designated by the alleles at MPRIP, TCAP, SREBF1, and NT5M in that order so as to emphasize the conservation over the 5 megabases between MPRIP and TCAP (M-T segment). Note in Figure 1 that MPRIP and TCAP are adjacent in Bota but not in Hosa.

The 60.10.S.10 is the most common Wagyu-specific haplotype. Other haplotypes are found in all breeds (e.g., 30.20.L.20 and 40.20.L.20) and yet others are characteristic of a breed such as Simmental (e.g., 30.10.L.10). These breed differences are shown as a W plot in Figure 4.

It can be seen that the conserved ancestral haplotypes extend from SREBF1 to FASN. The breed differences are especially striking at GH and FASN. Note that Wagyu-specific haplotypes have L at FASN whereas Simmental-specific haplotypes have S (Figure 4(a)). This remarkable difference might suggest that FASN is critical, but similar breed specificity is apparent at GH. Further, an association with 60.10.S.10.B.L or 30.20.S.20.C.L in Wagyu could be misinterpreted as implicating the shared alleles SREBF1 S and FASN L which are megabases apart. Note also that the remarkable but unpredictable associations between alleles at SREBF1, GH, and FASN have the potential to explain why there has been 20 years of confusion and argument as to which particular gene is important. As also shown in Figure 4(a), the five haplotypes in the family are also present in the W plot. Three are essentially Wagyu specific whereas two are common to cattle generally implying conservation over at least thousands of generations.

By comparing the haplotypes in the red and black Wagyu, it can be seen that the 60.10 segment at M-T is shared but the SREBF1 to NTM5 (S-N) segments differ (see Figure 2). Both types of Wagyu are known for their fine marbling suggesting that M-T is crucial.

2.4. Haplotypes in Crossbreds

A key question is which Wagyu haplotypes can have an effect in crossbreds when the genetic background is foreign. Samples were taken from exhibits submitted to a steak competition intended to compare production systems rather than genetics. The results of Tm measurements are shown in Figure 5. The full bloods are below 36°C in 7/9 as expected. Although the crossbreds, taken as a group, have higher Tm, remarkably, 4/16 are below 36°C. Of the 11/16 crossbreds with Tm below 38°C, 9/11 have one or the other of two Wagyu-specific haplotypes. Thus, these two haplotypes (60.10.S.10 and 30.20.S.20) are associated with, if not directly responsible for, the dramatic fall in Tm obtained by crossing a high Tm breed with a black Wagyu which transmits one of these haplotypes. Note that these two haplotypes differ in the M-T segment but share S at SREBF1.

So as to examine the effect of haplotypes in another setting, we compared marbling scores in crossbreds differing by whether they inherited haplotypes classified as specific for either Wagyu or Bos indicus. As shown in Table 2, the former group had greater marbling.

2.5. Mapping Mechanisms to Haplotypes

It is unrealistic to expect to map single functions to ancestral haplotypes as if these were just a string of biallelic coding genes. The sequence conservation is so extensive, so complex, and yet, so poorly understood! However, it is reasonable to expect that the differences will be quantitative and that penetrance will be dependent on multiple variables. At minimum, there will be multiple cis- and trans-interactions between the various polymorphisms. It follows that we expect to find associations at multiple markers. The results given above appear to implicate sequences around S and within the M-T segment of Wagyu haplotypes irrespective of effects marked by GH and FASN (see Figure 2).

In order to evaluate such conclusions and so as to assess other markers and haplotypes, we have taken advantage of the availability of two very different data sets:

Mayurafixed environment including 300 DOF, 100% black Wagyu, permitting examination of genetic differences within black Wagyu.

Melaleukafixed environment but variable DOF and breed, permitting examination of all haplotypes, Wagyu haplotypes within crossbreds and effects on the rate of marbling

Table 3 gives some of the results of the Mayura set. As already shown, see Figure 4, this is the set which revealed the heritability of Tm confirming its value in dissecting complex genetics. For convenience, we show a cut-off of 37°C, giving approximately 50% in the above and below groups. By comparing alleles at four of the loci shown in Figure 2, it is clear that the TCAP 20 allele is associated with lower Tm. The MPRIP 60 is associated with higher Tm. Although SREBF1 S is a marker for low Tm in some earlier studies and patents, here we find the opposite. At the NT5M locus, the 22 allele is associated with low Tm.

The inconsistencies between past and present results and between locus and haplotype analyses can be illustrated as follows. The prototypic Wagyu ancestral sequence, 60.10.S10, can be recognized by the 60.10.S.10 combination, but these alleles, taken individually, are not haplospecific. As shown in Figure 2, each is found on more than one haplotype.

When the alleles are examined without regard to whether they are carried by the ancestral haplotype, Tm would be higher at MPRIP, TCAP, SREB-P, and NT5M (see Table 3). However, as shown in Figure 5, the 60.10.S.10 ancestral haplotype actually lowers the Tm of crossbreds. Thus, there may be multiple cis- and trans-interactions contributing to the effect of a particular haplotype.

As a further approach to address the complexity, we asked whether there is any influence of trans-interaction or homozygosity. Since it would be necessary to study many hundreds of thousands of individuals to obtain sufficient homozygotes of even the most common haplotypes, we chose to examine the least polymorphic loci such as TCAP with only two alleles: 10 and 20. In the Mayura long-fed Wagyu data set, all 8 with Tm below 36°C are 20,20 homozygotes (see Figure 6). Using 37°C as a cut-off, 60% of 20,20 are lower compared with only 33% of the 10,10 homozygotes.

The effect of TCAP homozygosity is clear in a genetically very diverse group fed for 450 days. Only subcutaneous fat from over the rump was available. Although not comparable to other samples reported above, TCAP 10,10 homozygotes have higher Tm than 20,20 homozygotes (Table 4).

So as to address the issue in a very different data set, we analysed the Melaleuka Euro and crossbred cattle and found similar results. As shown in Figure 7, TCAP 10,10 and 20,20 homozygotes diverge after 50–100 DOF. After 100 DOF, 9 of 12 with TCAP 10,10 remain above 39°C compared with 3 of 12 with TCAP 20,20 (). The 20,20 homozygotes assume a faster trajectory to healthier fat indicating a quantitative regulatory effect (Figure 8).

A similar pattern is obtained when a model is constructed based on the slow induction of SCD (Figure 9).

3. Discussion

The history of Japanese Wagyu or “Kobe” beef is intriguing. According to legend [3], Wagyu cattle were used for heavy burden although fed only with rice husks, suggesting that they may have been selected for efficient food conversion combined with an ability to utilise fatty acids as a direct energy source. Much later, after the eating of meat was permitted, selection (according to the Shogun’s preferences) must have led to enhanced flavour as found with fine marbling and low melting temperature intramuscular fat. Another consequence has been a tendency to muscle weakness attributed to the replacement of myocytes with adipocytes. Whatever the selection criteria, Wagyu have conserved the capacity to marble from their remote North Asian ancestors and seemingly still do so after crossing with diverse breeds now found worldwide. This history suggested that there must be highly conserved and penetrant haplotypes rather than recent random mutations.

The present study began with the seemingly simple aim of improving the healthiness of beef by crossing elite Japanese Wagyu with composite breeds able to survive in harsh conditions. The literature suggests that transfer of a Wagyu gene might improve the fat content thereby providing cholesterol-lowering, statin-sparing, tasty beef at reasonable cost. Experience has shown that improvement can occur but inconsistently and only after expensive grain-feeding thereby restricting worldwide access to potential benefits. Our interpretation is that the outcome depends on the particular mix of conserved haplotypes which exert quantitative control over all the multiple functions which shift the balance from muscle to lipid and especially MUFA.

As shown elsewhere [60], it was necessary to develop a simple method of measuring the amount of beneficial fat. It transpires that heritability of the preferred fat can be demonstrated with Tm but not the standard visual method of quantifying marbling. Hitherto, therefore, selection of sires must have suffered. This is important because, using Tm, the sire effect is evident within one generation, meaning that success or failure can occur rapidly, leading to the remarkable scatter as demonstrated in Figure 2. Fortunately, Tm measurements should lead to progressive improvement.

Given accurate and precise measurement of the endpoints, we have used four data sets.

Samples from a tasting competition reveal that some crossbreds can be as desirable as the most elite Wagyu. Here, we offer a strategy based on the observation that transmission of a certain Wagyu haplotype, such as 60.10.S.10, may dramatically improve quality. If confirmed, a homozygous bull will benefit all progeny and remove some of the existing inconsistency and expense. In an ongoing multicentre study, we hope to be able to evaluate the benefits of sires which are homozygous for 30.20.S.20 and 60.10.S.10 in their anticipated order of benefit. However, because of relatively low haplotype frequencies of say 0.05, the natural frequencies of homozygotes will be less than 0.0025, meaning that the value of such homozygotes will only be clear if all bulls are haplotyped. The solution to the lack of homozygotes for use today is being addressed through artificial insemination and embryo transfer. In the meanwhile, we are developing some hypotheses which recognize that evaluating all combinations of multimegabase haplotypes (as diplotypes) is impractical. The n values required for the comparisons are greater than the number of cattle alive today.

One of these hypotheses arose from the ability to compare two totally different data sets. The Mayura set standardizes essentially all variables, other than the genetic differences between Wagyu. The Melaleuka set uses Tm/DOF as a measure to compare crossbreds.

Both show a remarkable effect of TCAP 10,10 versus 20,20. It should be emphasized that differences at TCAP itself cannot explain all findings such as the superiority of Wagyu over other breeds and the benefit of 60.10.S.10 in crossbreds.

Rather, taken together, the results argue for the importance of ancestral haplotypes rather than individual loci. This is illustrated in Figure 4. The alleles at GH and FASN are different depending upon the megabase haplotype, meaning that effects attributed to these individual alleles could reflect sequences more than 10 megabases away. Therefore, it is important to define the megabase haplotypes before sequencing or SNP analyses. No doubt, this explains why progress has been slow. Genomic analysis of polygenic traits like QTLs must consider the implications of synteny and the potential importance of paralogy, duplication, copy number, and retroviral regulators to list a few mechanisms [2, 3, 30, 31, 55]. As proposed here, duplication and copy number are more relevant than coding region differences; so too is regulation over extensive distances. Thus, SREBF1 is a master regulator which controls gene expression within the SREBF1 to FASN region via targeting transcriptional repressors of MPRIP, TCAP, and other muscle–related genes (BECN1, CACNB1, MYH2, MYH3, and MYH13); one of these repressors BHLHB 2 maps near FASN [55]. Note also that muscular dystrophy, limb girdle type 2G is associated in humans with TCAP polymorphism (OMIM#604488) and the histopathology resembles marbling in some respects. Lipid accumulation has been observed in limb girdle muscular dystrophy [61]. Another unexpected recent finding is the location of UTS2R which is associated with development of type II diabetes and fat deposition [44].

The balance between myogenesis and adipogenesis may be regulated in part by the TCAP-myostatin pathway in Wagyu as first suggested by Shibata et al. [62].

We conclude that genomic differences regulate the amount of MUFA such as oleic acid in beef and determine the rate at which Tm falls with DOF. We regard this finding as indicating, possibly for the first time in the setting of genomics, that the genetic dissection confirms the expectation that the important differences are conserved and quantitative, rather than qualitative, in keeping with our earlier studies on regulation of immunoglobulin and antibodies by MHC ancestral haplotypes [2, 1113].

Genomic markers for favorable rates of marbling will be of practical as well as theoretical interest. Reducing the time and cost on feed brings the potential to provide healthy beef to a wider sector of the world’s population with the secondary benefits of reducing cardiovascular disease and the billions of dollars spent on cholesterol lowering statins.

One of the most important lessons from the present study is that a haplotype sequence of multiple genes and regulators can have complex associations and confounding effects. The 60.10.S.10 haplotypes is Wagyu specific, beneficial in upgrading non-Wagyu in the tasting competition, directly or indirectly associated with low Tm but only the 3rd or 4th best in the ranking of the Wagyu haplotypes. We believe that further data will show that, as in the case of the human MHC, diplotypes and their quantitative cis- and trans-interactions will be critical in explaining apparent contradictions.

In conclusion, we recommend combining an accurate measurement of the phenotype with a genomic description of the genotype allowing identification of megabase ancestral haplotypes.

4. Materials and Methods

4.1. Full-Blood Wagyu with Identified Sires

Two cohorts of Wagyu steers () and heifers () were fed for approximately 300 ± 20 days on a proprietary ration within a commercial feedlot. One-gram samples of meat and intramuscular fat were taken from each carcass between the 10th and 11th rib. AUS-MEAT marbling score (MS) was scored between the 10th and 11th rib, with an average of 7.6 and a range from 2 to 11. Steers chosen for the sire comparison had their paternity confirmed via DNA testing [55].

4.2. Branded Beef Competition

Samples were obtained from 9 full blood, 16 crossbred, and 4 grass-fed entries to the Australian Wagyu Association’s “Branded Beef Competition” from 2014 to 2015. Crossbred entries were F1, F2, F3, and F4 Wagyu. Feed composition and days on feed varied, but all entries had more than 25% intramuscular fat, as assessed by camera image analysis [64].

4.3. Short-Fed European and Crossbreds

The herd at Melaleuka Stud has a variety of European breeds including Simmental, Gelbvieh, and Angus. This herd was selected to produce high-quality beef on pasture, finished within 2 to 4 months of supplemental feeding. Black Wagyu have been introduced into the herd as full blood or pure bred and for crossing with some of the European breeds. Those designated 50% are the F1 progeny of a mating between Wagyu and European, whereas 25% and 75% are the second cross progeny of F1s mated with European and Wagyu, respectively.

Akaushi sires have recently been used over some European breed cows with European breed bulls continuing to be used over the remaining cows.

Calves stay on milk until 4 months of age when they are weaned, and male calves are castrated. After weaning, they continue grazing Kikuyu and Ryegrass pasture until they reach 300 kg. Their feed is then supplemented with 9 mm EasyBeef pellets (Milne Feeds, Perth, Australia) ad libitum. The main ingredients of the EasyBeef pellets are lupins, barley, oats, wheat, and triticale. The nutritional composition, based on dry matter, is crude protein (min) 14.5%, metabolizable energy (est.) 11.0 MJ/kg, crude fiber (max) 20.0%, urea (max) 1.5%, and monensin 26.6 ppm.

The feeders are considered ready for slaughter when they reach a weight of 400 kg and are slaughtered to match demand. Some animals were kept on feed longer to test the effect of increased feeding on Tm and meat quality. The average live weight at slaughter for animals in this study was 461 kg, average age at slaughter was 15.4 months (range 8 to 23), and the average days on feed was 104 days (range 17 to 288). Body numbers from abattoirs were matched to farm records and pedigrees via their RFID tags, where possible identity was confirmed by in-house proprietary DNA testing [56].

Subcutaneous fat samples from the sirloin of these cattle were collected after boning and wet aging.

4.4. Long-Fed Wagyu Crossbreeds with Bos indicus Ancestry

Fifty samples of subcutaneous fat were taken from the rump of crossbred Wagyu steers fed for 450 days in a commercial feedlot. Samples were collected from the abattoir and frozen immediately. DNA was extracted and tested for MPRIP markers, and triglycerides which separated from the adipose tissue during DNA extraction were tested for Tm. Wagyu ancestry ranged from 48% to 96% with the majority being F1 and average Wagyu ancestry at 57%. 19 of the 50 animals had some Bos indicus content. Other dam breeds included Shorthorn and Angus.

4.5. DNA Extraction and C19 Haplotype Testing

Genomic DNA was extracted from all meat samples using the standard salting out method. Alleles of SREB, NT5M, MPRIP, TCAP, and GH markers were determined by PCR and capillary electrophoresis using the primers and method described in Lecomte et al. [55]. Alleles of FASN were determined using the primer and method described in Williamson et al. [56]. The FASN marker is more correctly known as SCT-FSN since the present marker is adjacent to the coding region for FASN, within the segmental duplication containing SECTM1. Haplotypes in heterozygous individuals were determined from segregation or alleles in related animals. The W plots were based on haplotype frequencies measured in Australian herds of registered Wagyu.

Conflicts of Interest

Collectively, the authors associated with the CY O’Connor ERADE Village Foundation as publication number MS1505 have an interest in the work described in this manuscript as it forms part of the Foundation’s intellectual property.

Authors’ Contributions

Roger L. Dawkins and Sally S. Lloyd conceived the paper’s outline, data collection, and writing plan. Roger L. Dawkins wrote the paper. Sally S. Lloyd, Edward J. Steele, and Roger L. Dawkins contributed to the analysis and correction of the various manuscript versions. Edward J. Steele was also responsible for the review of the literature summarized in Table 1. Jose L. Valenzuela contributed the samples from Melaleuka Stud.

Acknowledgments

This research was supported by the CY O’Connor ERADE Village Foundation as publication number MS1505. The authors would like to thank J.R. Dawkins and the Australian Wagyu Association for the samples and associated data. For beef and fat samples, the authors would also like to thank Scott de Bruin (Mayura Station, SA), Peter Gilmour (Irongate Wagyu, WA), Darren Hamblin (Strathdale Wagyu), and Dom Bayard (Global Reproductive Solutions, Goorambatt, VIC). The authors would like to thank Professor Frank Mastaglia and Sue Lester for their helpful comments.

Supplementary Materials

Supplementary Figure 1: The three dotplots show the similarity between regions of the bovine chromsome 19 on the horizontal axis and the human chromosome 17 on the vertical axis. Diagonal lines from top left to bottom right show homologous regions, while diagonals in the opposite regions show homologous regions that have been through a process of inversion. Colour shading matches to the regions shown in Figure 1.

  1. Supplementary Material
  2. Supplementary Material
  3. Supplementary Material